Published in Communications in Numerical Methods in Engineering Vol. 25 (7), pp. 837-849, 2009
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 algorithm 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.