A project by DCS Computing and CFDEMresearch
Velocity Verlet integration assumes that acceleration only depends on position , and does not depend on velocity. Then is it logical to perform Velocity Verlet integration on granular where there is velocity dependent terms in formula of force and consequently acceleration depends on the velocity?
Thanks
This website uses cookies to ensure you get the best experience on our website. By continuing to browse the site you are agreeing to our use of cookies. This website uses cookies. By continuing to browse the site you are agreeing to our use of cookies. More info
ckloss | Tue, 03/06/2012 - 11:07
Many people in DEM use
Many people in DEM use Velocity Verlet. If you simulate glass beads, i.e. coefficient of restitution = 0.98, this will not harm you. If the cor gets really low, this will "impair" the order of the Velocity Verlet integration scheme, i.e. it will be 1 < order < 2 for granular systems. If you think another integrator would fit better for your application, there are many other 2nd order integrators available, such as Adams-Bashforth and others
See eg "Selection of an appropriate time integration scheme for the discrete element method (DEM)" Computers & Chemical Engineering
Volume 32, Issue 10, 17 October 2008, Pages 2263–2279 for an overview. Any of these should be pretty easy to implement in LIGGGHTS.
Cheers, Christoph