Discrete element method

Not to be confused with finite element method.

A discrete element method (DEM), also called a distinct element method, is any of a family of numerical methods for computing the motion and effect of a large number of small particles. Though DEM is very closely related to molecular dynamics, the method is generally distinguished by its inclusion of rotational degrees-of-freedom as well as stateful contact, particle deformation and often complicated geometries (including polyhedra). With advances in computing power and numerical algorithms for nearest neighbor sorting, it has become possible to numerically simulate millions of particles on a single processor. Today DEM is becoming widely accepted as an effective method of addressing engineering problems in granular and discontinuous materials, especially in granular flows, powder mechanics, ice and rock mechanics. DEM has been extended into the Extended Discrete Element Method taking heat transfer,[1] chemical reaction[2] and coupling to CFD[3] and FEM[4] into account.

Discrete element methods are relatively computationally intensive, which limits either the length of a simulation or the number of particles. Several DEM codes, as do molecular dynamics codes, 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 average the physics across many particles and thereby treat the material as a continuum. In the case of solid-like granular behavior as in soil mechanics, the continuum approach usually treats the material as elastic or elasto-plastic and models it with the finite element method or a mesh free method. In the case of liquid-like or gas-like granular flow, the continuum approach may treat the material as a fluid and use computational fluid dynamics. Drawbacks to homogenization of the granular scale physics, however, are well-documented and should be considered carefully before attempting to use a continuum approach.

  1. ^ Peng, Z.; Doroodchi, E.; Moghtaderi, B. (2020). "Heat transfer modelling in Discrete Element Method (DEM)-based simulations of thermal processes: Theory and model development". Progress in Energy and Combustion Science. 79, 100847: 100847. doi:10.1016/j.pecs.2020.100847. S2CID 218967044.
  2. ^ Papadikis, K.; Gu, S.; Bridgwater, A.V. (2009). "CFD modelling of the fast pyrolysis of biomass in fluidised bed reactors: Modelling the impact of biomass shrinkage" (PDF). Chemical Engineering Journal. 149 (1–3): 417–427. doi:10.1016/j.cej.2009.01.036.
  3. ^ Kafui, K.D.; Thornton, C.; Adams, M.J. (2002). "Discrete particle-continuum fluid modelling of gas–solid fuidised beds". Chemical Engineering Science. 57 (13): 2395–2410. Bibcode:2002ChEnS..57.2395K. doi:10.1016/S0009-2509(02)00140-9.
  4. ^ Trivino, L.F.; Mohanty, B. (2015). "Assessment of crack initiation and propagation in rock from explosion-induced stress waves and gas expansion by cross-hole seismometry and FEM–DEM method". International Journal of Rock Mechanics & Mining Sciences. 77: 287–299. Bibcode:2015IJRMM..77..287T. doi:10.1016/j.ijrmms.2015.03.036.