A numerical study of viscoelastic phase separation in polymer solutions