An implementation of granular dynamics for simulating frictional elastic particles based on the DL_POLY code

We have modified Daresbury Laboratory's replicated data strategy (RDS) parallel molecular dynamics (MD) package DL_POLY (version 2.13) to study the granular dynamics of frictional elastic particles. DL_POLY [Smith and Forester, The DL_POLY_2 User Manual v2.13, 2001; Forester and Smith, The DL_POLY_2 Reference Manual v2.13, 2001] is a MD package originally developed to study liquid state and macromolecular systems by accounting for various molecular interaction forces. The particles of interest in this study are macroscopic grains in pharmaceutical powders, with sizes ranging from tens to hundreds of microns. We have therefore substituted the molecular interaction forces with contact forces (including linear-dashpot, HKK interaction forces and Coulombic friction) while taking advantage of the RDS scheme. In effect, we have created a parallel Discrete Element Simulation (DES) code. In this paper, we describe the modifications made to the original DL_POLY code and the results from the validation tests of the granular dynamics simulations for systems of monodisperse spherical particles settling under gravity. The code can also be utilized to study particle packings generated via uniaxial compaction and, in some cases, simultaneous application of shear, at constant strain.

[1]  William Smith,et al.  Molecular dynamics on hypercube parallel computers , 1991 .

[2]  Sam F. Edwards,et al.  The equations of stress in a granular material , 1998 .

[3]  M. Aulton.  Pharmaceutics : the science of dosage form design , 2002 .

[4]  D. Wolf,et al.  Force Schemes in Simulations of Granular Materials , 1996 .

[5]  P. Cundall,et al.  A discrete numerical model for granular assemblies , 1979 .

[6]  G. Kuwabara,et al.  Restitution Coefficient in a Collision between Two Spheres , 1987 .

[7]  S. Timoshenko,et al.  Theory of elasticity , 1975 .

[8]  William Smith,et al.  Molecular dynamics on distributed memory (MIMD) parallel computers , 1993 .

[9]  Kurt Kremer,et al.  Vectorized link cell Fortran code for molecular dynamics simulations for a large number of particles , 1989 .

[10]  D. Lin,et al.  Structure, stability and evolution of Saturn's rings , 1984, Nature.

[11]  Herrmann,et al.  Granular packings and fault zones , 2000, Physical review letters.

[12]  D. H. Zeuch,et al.  Mechanical properties and shear failure surfaces for two alumina powders in triaxial compression , 2000 .

[13]  Stefan Brode,et al.  An optimized MD program for the vector computer cyber 205 , 1986 .

[14]  H. Jaeger,et al.  Granular solids, liquids, and gases , 1996 .

[15]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[16]  I. Goldhirsch,et al.  A molecular dynamical study of granular fluids I: The unforced granular gas in two dimensions , 1993 .

[17]  C. Brooks.  Computer simulation of liquids , 1989 .

[18]  J. Barbera,et al.  Contact mechanics , 1999 .

[19]  Brian J. Briscoe,et al.  The internal form of compacted ceramic components: a comparison of a finite element modelling with experiment , 1996 .

[20]  Y-h. Taguchi,et al.  Powder turbulence: direct onset of turbulent flow , 1992 .

[21]  H. Herrmann,et al.  Physics of dry granular media , 1998 .

[22]  C. Brennen,et al.  Measurements of Solid Spheres Bouncing Off Flat Plates , 1990 .

[23]  Spahn,et al.  Model for collisions in granular gases. , 1996, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[24]  A. Cocks,et al.  Experimental investigation of yield behaviour of metal powder compacts , 2002 .

[25]  S. Alexander,et al.  Amorphous solids: their structure, lattice dynamics and elasticity , 1998 .

[26]  Gary S Grest,et al.  Geometry of frictionless and frictional sphere packings. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[27]  David Fincham,et al.  The Ewald Sum in Truncated Octahedral and Rhombic Dodecahedral Boundary Conditions , 1993 .