Numerical simulations of granular dynamics. I. Hard-sphere discrete element method and tests

Abstract

We present a new particle-based (discrete element) numerical method for the simulation of granular dynamics, with application to motions of particles on small solar system body and planetary surfaces. The method employs the parallel N-body tree code pkdgrav to search for collisions and compute particle trajectories. Collisions are treated as instantaneous point-contact events between rigid spheres. Particle confinement is achieved by combining arbitrary combinations of four provided wall primitives, namely infinite plane, finite disk, infinite cylinder, and finite cylinder, and degenerate cases of these. Various wall movements, including translation, oscillation, and rotation, are supported. We provide full derivations of collision prediction and resolution equations for all geometries and motions. Several tests of the method are described, including a model granular "atmosphere" that achieves correct energy equipartition, and a series of tumbler simulations that show the expected transition from tumbling to centrifuging as a function of rotation rate.

54 manuscript pages, 8 figures including 4 in colour (online version only)

References in corpus (2)

• Microgravity experiments on the collisional behavior of Saturnian ring particles
• The S shape of a granular pile in a rotating drum

Cited by in corpus (19)

• Landslides and Mass Shedding on Spinning Spheroidal Asteroids
• Numerical predictions of surface effects during the 2029 close approach of asteroid 99942 Apophis
• The Brazil-nut effect and its application to asteroids
• Coupled Spin and Shape Evolution of Small Rubble-Pile Asteroids: Self-Limitation of the YORP Effect
• Tidal distortion and disruption of rubble-pile bodies revisited -- Soft-sphere discrete element analyses
• The Sands of Phobos: The Martian moon's eccentric orbit refreshes its surface
• Low-speed Impact Simulations into Regolith in Support of Asteroid Sampling Mechanism Design I.: Comparison with 1-g Experiments
• The spherical Brazil Nut Effect and its significance to asteroids
• Numerical Simulations of Collisional Disruption of Rotating Gravitational Aggregates: Dependence on Material Properties
• The role of fragment shapes in the simulations of asteroids as gravitational aggregates
• Validating N-body code Chrono for granular DEM simulations in reduced-gravity environments
• Numerical Simulations of Saturn's B Ring: Granular Friction as a Mediator between Self-Gravity and Viscous Overstability
• Crustal Failure on Icy Moons from a Strong Tidal Encounter
• Small-body deflection techniques using spacecraft: techniques in simulating the fate of ejecta
• The formation mechanism of 4179 Toutatis' elongated bi-lobed structure in a close Earth encounter scenario
• Numerical simulations of granular dynamics II. Particle dynamics in a shaken granular material
• Synchronization of Conservative Parallel Discrete Event Simulations on a Small-World Network
• Considering contact forces during the formation of planetesimals by gravitational collapse: mutual orbits, spin states, and shapes
• Space weathering of asteroids