Estimate Simulation Time

It is challenging to estimate the exact real time required for a DEM simulation, as each simulation and each computer is different.

Time Step

One of the key numbers in DEM simulation is the Rayleigh Time Step which is the time taken for a shear wave to propagate through a solid particle. It is, therefore, a theoretical maximum Time Step for a DEM simulation of a quasi-static particulate collection in which the coordination number (total number of contacts per particle) for each particle remains above 1 and is defined as:

TR=πRρG0.1631υ+0.8766

where R is the particle’s radius, ρ its density, G the Shear modulus, and v the Poisson’s ratio. This formula assumes that the relative velocity between contacting particles is very small. Other than for quasi-static systems, in practice some fraction of this maximum value is used, and for high coordination numbers (4 and above) a typical Time Step of 0.2TR (20%) is recommended. For lower coordination numbers, 0.4TR (40%) is more appropriate.

Hertzian Contact

While the Rayleigh Time Step is a suitable starting point for quasi-static simulations, a shorter Time Step is required for systems undergoing flow. Consider two elements approaching each other at a speed v. In one Time Step, t, the maximum possible overlap is defined as:

dmax=υt

In DEM, particles undergoing elastic (Hertzian) contact are treated as overlapping and this overlap is equated to a surface compression. t in the above equation must be such that the maximum overlap is lower than the theoretical maximum overlap for Hertzian contact. In practice, to get a good numerical integral to the contact graph, at least six time points should occur - three during approach and three during separation (though ten is more desirable) .

Elastic Impact

From the Hertz theory of elastic collision, the total time of contact is defined as:

Th=2.87(m*2R*E*2Vz)15
where
1m*=1mi+1mj
where mi and mj are the masses of the two elements.
1R*=1Ri+1Rj
where Ri and Rj are the radii of the two elements.
1E*=1υ2iEi+1υ2jEj
where νi and νj are the Poisson’s ratios of the two materials.
Vz=ViVj
where Vz is the relative velocity and Vi and Vj are the velocities of the two elements.