Abstract:We summarize and numerically compare two approaches for modeling and simulating the dynamics of dry granular matter. The first one, the discrete-element method via penalty (DEM-P), is commonly used in the soft matter physics and geomechanics communities; it can be traced back to the work of Cundall and Strack [P. Cundall, Proc. Symp. ISRM, Nancy, France 1, 129 (1971); P. Cundall and O. Strack, Geotechnique 29, 47 (1979)GTNQA80016-850510.1680/geot.1979.29.1.47]. The second approach, the discrete-element method via complementarity (DEM-C), considers the grains perfectly rigid and enforces nonpenetration via complementarity conditions; it is commonly used in robotics and computer graphics applications and had two strong promoters in Moreau and Jean [J. J. Moreau, in Nonsmooth Mechanics and Applications, edited by J. J. Moreau and P. D. Panagiotopoulos (Springer, Berlin, 1988), pp. 1-82; J. J. Moreau and M. Jean, Proceedings of the Third Biennial Joint Conference on Engineering Systems and Analysis, Montpellier, France, 1996, pp. 201-208]. The DEM-P and DEM-C are manifestly unlike each other: They use different (i) approaches to model the frictional contact problem, (ii) sets of model parameters to capture the physics of interest, and (iii) classes of numerical methods to solve the differential equations that govern the dynamics of the granular material. Herein, we report numerical results for five experiments: shock wave propagation, cone penetration, direct shear, triaxial loading, and hopper flow, which we use to compare the DEM-P and DEM-C solutions. This exercise helps us reach two conclusions. First, both the DEM-P and DEM-C are predictive, i.e., they predict well the macroscale emergent behavior by capturing the dynamics at the microscale. Second, there are classes of problems for which one of the methods has an advantage. Unlike the DEM-P, the DEM-C cannot capture shock-wave propagation through granular media. However, the DEM-C is proficient at handling arbitrary grain geometries and solves, at large integration step sizes, smaller problems, i.e., containing thousands of elements, very effectively. The DEM-P vs DEM-C comparison is carried out using a public-domain, open-source software package; the models used are available online.

DEM Contact Model for Spherical and Polyhedral Particles Based on Energy Conservation

Discrete Element Modeling of Polyhedral Representation of Granular Materials

Probability-based Contact Algorithm for Non-Spherical Particles in DEM

Investigation of Modeling Non-Spherical Particles by Using Spherical Discrete Element Model with Rolling Friction

Multi-super-ellipsoid model for non-spherical particles in DEM simulation

Comparative Study of Discrete Element Modeling of Tablets Using Multi-Spheres, Multi-Super-ellipsoids, and Polyhedrons

Dem Investigation Of Sag Mill With Spherical Grinding Media And Non-Spherical Ore Based On Polyhedron-Sphere Contact Model

Discrete element modeling and simulation of non-spherical particles using polyhedrons and super-ellipsoids

The Influence of Particle-Size Distribution on Critical State Behavior of Spherical and Non-Spherical Particle Assemblies

A GPU-based DEM framework for simulation of polyhedral particulate system

Study of discrete element model for non-sphere particles base on super-quadrics

A deformable spheropolygon-based discrete element method

Development and Verification of Coarse‐grain CFD‐DEM for Nonspherical Particles in a Gas–solid Fluidized Bed

Geometrically Exact Discrete-Element-method (DEM) Simulation on the Flow and Mixing of Sphero-Cylinders in Horizontal Drums

GPU based discrete element modeling for convex polyhedral shape particles: Development and validation

A new discrete element method for small adhesive non-spherical particles

An Experimentally Validated Coarse-Grain DEM Study of Monodisperse Granular Mixing

Compliant contact versus rigid contact: A comparison in the context of granular dynamics

Packing of Cylindrical Particles: DEM Simulations and Experimental Measurements

CFD-DEM Modeling of Rod-Like Particles in a Fluidized Bed with Complex Geometry

A poly‐superellipsoid‐based approach on particle morphology for DEM modeling of granular media