Angle of repose and angle of marginal stability: molecular dynamics of granular particles

Abstract

The authors present an implementation of realistic static friction in molecular dynamics (MD) simulations of granular particles. In their model, to break contact between two particles, one has to apply a finite amount of force, determined by the Coulomb criterion. Using a two-dimensional model, they show that piles generated by avalanches have a finite angle of repose theta _R (finite slopes). Furthermore, these piles are stable under tilting by an angle smaller than a non-zero tilting angle theta _T, showing that theta _R is different from the angle of marginal stability theta _MS, which is the maximum angle of stable piles. These measured angles are compared to a theoretical approximation. They also measure theta _MS by continuously adding particles on the top of a stable pile.