High-density sphere packing for discrete element method simulations

The first step in a discrete element simulation is the discretization of the domain into a set of particles. The cost of generating a good cylindrical or spherical packing has resulted in a great number of approaches during the last years. A new algorithin is proposed for high-density packing using a scheme that minimizes the distance between each particle. Using the Support of a finite element mesh, less time is needed in order to achieve a low porosity configuration. In addition, a boundary constraint is introduced. The application of the same optimization scheme is used as a condition to force a good surface definition. The results obtained show a high efficiency for the generation of low porosity packing, achieving values smaller than 10% in 2D cases and 30% in 3D cases. Copyright (C) 2008 John Wiley & Sons, Ltd.