Optimum Cell Size for Contact Detection in the Algorithm of the Discrete Element Method

The objective of this paper is to find out the optimum cell size for contact detection in the algorithm of DEM (Discrete Element Method) to improve the calculation speed. The computing time of contact detection was measured by using a cubic box, which contains particles, to investigate the effect of cell size. The contact detection was inefficient when too fine grid cell or large one was introduced, the optimum cell size was given around 1.5–2.0 times of particle radius. The theoretical formula for the optimum cell size was proposed by considering the effects of the number of searched cells and contact checks. A significant correlation between the computing time of DEM work and the theoretical equation was seen in all conditions of volume fraction. Furthermore, it was found that the optimum cell size was the size, which each cell has about 0.7–0.8 particles in it, and this result was applied to a modelling of random-sized particles. Therefore, the calculation speed of DEM would be improved by using the optimum cell size for the contact detection.

[1]  B. Alder,et al.  Studies in Molecular Dynamics. I. General Method , 1959 .

[2]  Thomas M. Liebling,et al.  Parallel implementation of a distinct element method for granular media simulation on the Cray T3D , 1996 .

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

[4]  Paul Langston,et al.  Validation tests on a distinct element model of vibrating cohesive particle systems , 2002 .

[5]  Jpk Seville,et al.  The influence of DEM simulation parameters on the particle behaviour in a V-mixer , 2002 .

[6]  Paul W. Cleary,et al.  DEM modelling of industrial granular flows: 3D case studies and the effect of particle shape on hopper discharge , 2002 .

[7]  A. Yu,et al.  Microdynamic analysis of particle flow in a horizontal rotating drum , 2003 .

[8]  J. Kano,et al.  Estimation of mechanochemical dechlorination rate of poly(vinyl chlorde). , 2002, Environmental Science and Technology.

[9]  Mojtaba Ghadiri,et al.  Effect of the impact angle on the breakage of agglomerates: a numerical study using DEM , 2003 .

[10]  Junya Kano,et al.  Correlation of grinding rate of gibbsite with impact energy in tumbling mill with mono-size balls , 2001 .

[11]  Paul W. Cleary,et al.  A parallel discrete element method for industrial granular flow simulations , 1999 .

[12]  A. J. Matchett,et al.  Vibrating powder beds: a comparison of experimental and Distinct Element Method simulated data , 2000 .

[13]  Garry H. Rodrigue,et al.  A High Performance Chemical Kinetics Algorithm for 3-D Atmospheric Models , 1999, Int. J. High Perform. Comput. Appl..

[14]  S. Masson,et al.  Effect of particle mechanical properties on silo flow and stresses from distinct element simulations , 2000 .

[15]  Mridul Gautam,et al.  Computer simulation of bubbles in large-particle fluidized beds , 1998 .

[16]  S. Nakatsuka,et al.  The Terzaghi-Voigt combined model for constant-pressure consolidation of filter cakes and homogeneous semi-solid materials , 1986 .