The discrete element method (DEM) is a family of numerical methods for computing the motion of a large number of particles like molecules or grains of sand. The method is sometimes called molecular dynamics (MD), even when the particles are not molecules. Two prominent members of the DEM family are the distinct element method proposed by Cundall in 1971, and the discontinuous deformation analysis (DDA) proposed by Shi in 1988.
Discrete element methods are processor intensive and this limits either the length of a simulation or the number of particles. Advances in the software are beginning to take advantage of parallel processing capabilities (shared or distributed systems) to scale up the number of particles or length of the simulation. An alternative to treating all particles separately is to treat the material as a continuum and use computational fluid dynamics.
Applications
The fundamental assumption of the method is that the material consists of separate, discrete particles. These particles may have different shapes and properties. Some examples are:
- liquids and solutions, for instance of sugar or proteins;
- bulk materials in storage silos, like cereal;
- granular matter, like sand;
- powders, like toner.
Typical industries using DEM are:
- Mining
- Pharmaceutical
- Oil and gas
- Agriculture and food handling
- Chemical
Outline of the method
A DEM-simulation is started by putting all particles in a certain position and giving them an initial velocity. Then the forces which act on each particle are computed from the initial data and the relevant physical laws.
The following forces may have to be considered in macroscopic simulations:
- friction, when two particles touch each other;
- recoil, when two particles collide;
- gravity (the force of attraction between particles due to their mass), which is only relevant in astronomical simulations.
On a molecular level, we may consider
- the Coulomb force, the electrostatic attraction or repulsion of particles carrying electric charge;
- Pauli repulsion, when two atoms approach each other closely;
- van der Waals force.
All these forces are added up to find the total force acting on each particle. An integration method is employed to compute the change in the position and the velocity of each particle during a certain time step from Newton's laws of motion. Then, the new positions are used to compute the forces during the next step, and this loop is repeated until the simulation ends.
Typical integration methods used in a discrete element method are:
- the Verlet algorithm,
- velocity Verlet,
- the leapfrog method.
Long-range forces
When long-range forces (typically gravity or the Coulomb force) are taken into account, then the interaction between each pair of particles needs to be computed. The number of interactions, and with it the cost of the computation, increases quadratically with the number of particles. This is not acceptable for simulations with large number of particles. A possible way to avoid this problem is to combine some particles, which are far away from the particle under consideration, into one pseudoparticle. Consider as an example the interaction between a star and a distant galaxy: The error arising from combining all the stars in the distant galaxy into one point mass is negligible. So-called tree algorithms are used to decide which particles can be combined into one pseudoparticle. These algorithms arrange all particles in a tree, a quadtree in the two-dimensional case and an octree in the three-dimensional case.
However, simulations in molecular dynamics divide the space in which the simulation take place into cells. Particles leaving through one side of a cell are simply inserted at the other side (periodic boundary conditions); the same goes for the forces. The force is no longer taken into account after the so-called cut-off distance (usually half the length of a cell), so that a particle is not influenced by the mirror image of the same particle in the other side of the cell. One can now increase the number of particles by simply copying the cells.
Algorithms to deal with long-range force include:
Bibliography
- P.A. Cundall, O.D.L. Strack, A distinct element model for granular assemblies. Geotechnique, 29:47–65, 1979.
- Griebel, Knapek, Zumbusch, Caglar: Numerische Simulation in der Molekulardynamik. Springer, 2004. ISBN 3-540-41856-3.
- Bicanic, Ninad, Discrete Element Methods in Stein, de Borst, Hughes Encyclopedia of Computational Mechanics, Vol. 1. Wiley, 2004. ISBN 0-470-84699-2.
- Shi, G, Discontinuous deformation analysis - A new numerical model for the statics and dynamics of deformable block structures, 16pp. In 1st U.S. Conf. on Discrete Element Methods, Golden. CSM Press: Golden, CO, 1989.
- Kawaguchi, T., Tanaka, T. and Tsuji, Y., "Numerical simulation of two-dimensional fluidized beds using the discrete element method (comparison between the two- and three-dimensional models), Powder Technology, 96(2):129–138, 1998.
Software
Open source and non-commercial software:
- BALL & TRUBAL (1979-1980) distinct element method (FORTRAN code), originally written by P.Cundall and currently maintained by C.Thornton.
- SDEC Spherical Discrete Element Code.
- YADE Yet Another Dynamic Engine, second incarnation of SDEC written from ground-up, GPL license.
Commercially available DEM software packages include PFC3D and EDEM:
- PFC2D and PFC3D (Particle Flow Code in 2 Dimensions; Particle Flow Code in 3 Dimensions), PFC2D uses BALL codebase, PFC3D uses TRUBAL codebase.
- EDEM (DEM Solutions Ltd.)
- GROMOS 96
- ELFEN