Time Stepping Schemes

edited December 2016 in DualSPHysics_v4.0
Hello,

I am having some trouble understanding the Symplectic time stepping scheme described in the User Guide (p. 19f). I compared the equations 28-29 to the ones in Monaghan (2005), which is cited. I found that they match equations 5.39-5.41 in Monaghans paper, except for one point in equation 29 in the user guide, where the v(t=n+1) = v(t=n+1/2)+... . In Monaghan this equation is given as v(t=n+1) = v(t=n)+... , withouth having a predicted velocity at t=n+1/2.

So is it possible that the "+1/2" in the UserGuide is wrong? Or is this scheme a different one? If so, how to calculate v(t=n+1/2)?

Regards

Comments

  • edited December 2016
    And I have another question, concerning the Verlet scheme. In other publications I found similar equations, just that they always divided the force by m (mass). Why is that not necessary here?
Sign In or Register to comment.