A numerical study of the SPH method for simulating transient viscoelastic free surface flows