Stable algorithm for event detection in event-driven particle dynamics

Abstract

Event-Driven Particle Dynamics is a fast and precise method to simulate particulate systems of all scales. In this work it is demonstrated that, despite the high accuracy of the method, the finite machine precision leads to simulations entering invalid states where the dynamics are undefined. A general event-detection algorithm is proposed which handles these situations in a stable and efficient manner. This requires a definition of the dynamics of invalid states and leads to improved algorithms for event-detection in hard-sphere systems.

References in corpus (6)

• Coefficient of Restitution for Viscoelastic Spheres: The Effect of Delayed Recovery
• Coefficient of Restitution as a Fluctuating Quantity
• Neighbor List Collision-Driven Molecular Dynamics Simulation for Nonspherical Particles. I. Algorithmic Details II. Applications to Ellipses and Ellipsoids
• Event-Driven Dynamics of Rigid Bodies Interacting via Discretized Potentials
• Event-driven Molecular Dynamics of Soft Particles
• Mapping continuous potentials to discrete forms

Cited by in corpus (7)

• Scalable Metropolis Monte Carlo for simulation of hard shapes
• Complex Crystals from Size-disperse Spheres
• Collective behavior of Vicsek particles without and with obstacles
• Hard sphere crystal nucleation rates: Reconciliation of simulation and experiment
• Evaluation of memory effects at phase transitions and during relaxation processes
• retQSS: A Novel Methodology for Efficient Modeling and Simulation of Particle Systems in Reticulated Geometries
• Stable algorithm for event detection in event-driven particle dynamics: logical states