Imaginary Coating Algorithm Approaching Dense Accumulation of Granular Material in Simulations with Discrete Element Method

Abstract

The difficulty of obtaining a densely packed granular material as an initial condition is a very common problem in numerical simulations of granular materials. In this article, an Imaginary Coating Algorithm (ICA) is introduced. To avoid unreasonable particle deformation when using a longer time step and a lower Young’s modulus, the radii used in calculating the action forces in a binary collision are slightly larger than the real values. In other words, an imaginary coat is added to each particle or element. To validate this algorithm, simulations were carried out by using ÅDEM, and a A Discrete Element Method (DEM) software program was developed. Compared with traditional Simulated Annealing Algorithms (SAA), this technique can approach the densely packed state with less CPU/GPU time and is easy to operate.