Discrete element method

Last updated

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.

Contents

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.

The DEM family

The various branches of the DEM family are the distinct element method proposed by Peter A. Cundall and Otto D. L. Strack in 1979, [5] the generalized discrete element method ( Williams, Hocking & Mustoe 1985 ) [6] , the discontinuous deformation analysis (DDA) ( Shi 1992 ) and the finite-discrete element method concurrently developed by several groups (e.g., Munjiza and Owen). The general method was originally developed by Cundall in 1971 to problems in rock mechanics. Williams [7] showed that DEM could be viewed as a generalized finite element method, allowing deformation and fracturing of particles. Its application to geomechanics problems is described in the book Numerical Methods in Rock Mechanics( Williams, Pande & Beer 1990 ). The 1st, 2nd and 3rd International Conferences on Discrete Element Methods have been a common point for researchers to publish advances in the method and its applications. Journal articles reviewing the state of the art have been published by Williams and O'Connnor, [8] Bicanic, and Bobet et al. (see below). A comprehensive treatment of the combined Finite Element-Discrete Element Method is contained in the book The Combined Finite-Discrete Element Method. [9]

Discrete-element simulation with particles arranged after a photo of Peter A. Cundall. As proposed in Cundall and Strack (1979), grains interact with linear-elastic forces and Coulomb friction. Grain kinematics evolve through time by temporal integration of their force and torque balance. The collective behavior is self-organizing with discrete shear zones and angles of repose, as characteristic to cohesionless granular materials. Cundall DEM.gif
Discrete-element simulation with particles arranged after a photo of Peter A. Cundall. As proposed in Cundall and Strack (1979), grains interact with linear-elastic forces and Coulomb friction. Grain kinematics evolve through time by temporal integration of their force and torque balance. The collective behavior is self-organizing with discrete shear zones and angles of repose, as characteristic to cohesionless granular materials.

Applications

The fundamental assumption of the method is that the material consists of separate, discrete particles. These particles may have different shapes and properties that influence inter-particle contact. Some examples are:

Typical industries using DEM are:

Outline of the method

A DEM-simulation is started by first generating a model, which results in spatially orienting all particles and assigning an initial velocity. The forces which act on each particle are computed from the initial data and the relevant physical laws and contact models. Generally, a simulation consists of three parts: the initialization, explicit time-stepping, and post-processing. The time-stepping usually requires a nearest neighbor sorting step to reduce the number of possible contact pairs and decrease the computational requirements; this is often only performed periodically.

The following forces may have to be considered in macroscopic simulations:

On a molecular level, we may consider:

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:

Thermal DEM

The discrete element method is widely applied for the consideration of mechanical interactions in many-body problems, particularly granular materials. Among the various extensions to DEM, the consideration of heat flow is particularly useful. Generally speaking in Thermal DEM methods, the thermo-mechanical coupling is considered, whereby the thermal properties of an individual element are considered in order to model heat flow through a macroscopic granular or multi-element medium subject to a mechanical loading. [12] Interparticle forces, computed as a part of classical DEM, are used to determined areas of true interparticle contact and thus model the conductive transfer of heat from one solid element to another. A further aspect that is considered in DEM is the gas phase conduction, radiation and convection of heat in the interparticle spaces. To facilitate this, properties of the inter-element gaseous phase need to be considered in terms of pressure, gas conductivity and the mean-free path of gas molecules. [13]

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. Both the number of interactions and cost of computation increase 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:

Combined finite-discrete element method

Following the work by Munjiza and Owen, the combined finite-discrete element method has been further developed to various irregular and deformable particles in many applications including pharmaceutical tableting, [14] packaging and flow simulations, [15] and impact analysis. [16]

Advantages and limitations

Advantages

Disadvantages

See also

Related Research Articles

<span class="mw-page-title-main">Computational fluid dynamics</span> Analysis and solving of problems that involve fluid flows

Computational fluid dynamics (CFD) is a branch of fluid mechanics that uses numerical analysis and data structures to analyze and solve problems that involve fluid flows. Computers are used to perform the calculations required to simulate the free-stream flow of the fluid, and the interaction of the fluid with surfaces defined by boundary conditions. With high-speed supercomputers, better solutions can be achieved, and are often required to solve the largest and most complex problems. Ongoing research yields software that improves the accuracy and speed of complex simulation scenarios such as transonic or turbulent flows. Initial validation of such software is typically performed using experimental apparatus such as wind tunnels. In addition, previously performed analytical or empirical analysis of a particular problem can be used for comparison. A final validation is often performed using full-scale testing, such as flight tests.

Computational science, also known as scientific computing, technical computing or scientific computation (SC), is a division of science that uses advanced computing capabilities to understand and solve complex physical problems. This includes

Numerical methods for partial differential equations is the branch of numerical analysis that studies the numerical solution of partial differential equations (PDEs).

<span class="mw-page-title-main">Mesh generation</span> Subdivision of space into cells

Mesh generation is the practice of creating a mesh, a subdivision of a continuous geometric space into discrete geometric and topological cells. Often these cells form a simplicial complex. Usually the cells partition the geometric input domain. Mesh cells are used as discrete local approximations of the larger domain. Meshes are created by computer algorithms, often with human guidance through a GUI, depending on the complexity of the domain and the type of mesh desired. A typical goal is to create a mesh that accurately captures the input domain geometry, with high-quality (well-shaped) cells, and without so many cells as to make subsequent calculations intractable. The mesh should also be fine in areas that are important for the subsequent calculations.

Computational mechanics is the discipline concerned with the use of computational methods to study phenomena governed by the principles of mechanics. Before the emergence of computational science as a "third way" besides theoretical and experimental sciences, computational mechanics was widely considered to be a sub-discipline of applied mechanics. It is now considered to be a sub-discipline within computational science.

<span class="mw-page-title-main">Fluid–structure interaction</span>

Fluid–structure interaction (FSI) is the interaction of some movable or deformable structure with an internal or surrounding fluid flow. Fluid–structure interactions can be stable or oscillatory. In oscillatory interactions, the strain induced in the solid structure causes it to move such that the source of strain is reduced, and the structure returns to its former state only for the process to repeat.

<span class="mw-page-title-main">Meshfree methods</span> Methods in numerical analysis not requiring knowledge of neighboring points

In the field of numerical analysis, meshfree methods are those that do not require connection between nodes of the simulation domain, i.e. a mesh, but are rather based on interaction of each node with all its neighbors. As a consequence, original extensive properties such as mass or kinetic energy are no longer assigned to mesh elements but rather to the single nodes. Meshfree methods enable the simulation of some otherwise difficult types of problems, at the cost of extra computing time and programming effort. The absence of a mesh allows Lagrangian simulations, in which the nodes can move according to the velocity field.

Discontinuous deformation analysis (DDA) is a type of discrete element method (DEM) originally proposed by Shi in 1988. DDA is somewhat similar to the finite element method for solving stress-displacement problems, but accounts for the interaction of independent particles (blocks) along discontinuities in fractured and jointed rock masses. DDA is typically formulated as a work-energy method, and can be derived using the principle of minimum potential energy or by using Hamilton's principle. Once the equations of motion are discretized, a step-wise linear time marching scheme in the Newmark family is used for the solution of the equations of motion. The relation between adjacent blocks is governed by equations of contact interpenetration and accounts for friction. DDA adopts a stepwise approach to solve for the large displacements that accompany discontinuous movements between blocks. The blocks are said to be "simply deformable". Since the method accounts for the inertial forces of the blocks' mass, it can be used to solve the full dynamic problem of block motion.

<span class="mw-page-title-main">Computational engineering</span>

Computational Engineering is an emerging discipline that deals with the development and application of computational models for engineering, known as Computational Engineering Models or CEM. Computational engineering uses computers to solve engineering design problems important to a variety of industries. At this time, various different approaches are summarized under the term Computational Engineering, including using computational geometry and virtual design for engineering tasks, often coupled with a simulation-driven approach In Computational Engineering, algorithms solve mathematical and logical models that describe engineering challenges, sometimes coupled with some aspect of AI, specifically Reinforcement Learning.

The material point method (MPM) is a numerical technique used to simulate the behavior of solids, liquids, gases, and any other continuum material. Especially, it is a robust spatial discretization method for simulating multi-phase (solid-fluid-gas) interactions. In the MPM, a continuum body is described by a number of small Lagrangian elements referred to as 'material points'. These material points are surrounded by a background mesh/grid that is used to calculate terms such as the deformation gradient. Unlike other mesh-based methods like the finite element method, finite volume method or finite difference method, the MPM is not a mesh based method and is instead categorized as a meshless/meshfree or continuum-based particle method, examples of which are smoothed particle hydrodynamics and peridynamics. Despite the presence of a background mesh, the MPM does not encounter the drawbacks of mesh-based methods which makes it a promising and powerful tool in computational mechanics.

In contact mechanics, the term unilateral contact, also called unilateral constraint, denotes a mechanical constraint which prevents penetration between two rigid/flexible bodies. Constraints of this kind are omnipresent in non-smooth multibody dynamics applications, such as granular flows, legged robot, vehicle dynamics, particle damping, imperfect joints, or rocket landings. In these applications, the unilateral constraints result in impacts happening, therefore requiring suitable methods to deal with such constraints.

The CFD-DEM model, or Computational Fluid Dynamics / Discrete Element Method model, is a process used to model or simulate systems combining fluids with solids or particles. In CFD-DEM, the motion of discrete solids or particles phase is obtained by the Discrete Element Method (DEM) which applies Newton's laws of motion to every particle, while the flow of continuum fluid is described by the local averaged Navier–Stokes equations that can be solved using the traditional Computational Fluid Dynamics (CFD) approach. The interactions between the fluid phase and solids phase is modeled by use of Newton's third law.

A CFD-DEM model is suitable for the modeling or simulation of fluid-solids or fluid-particles systems. In a typical CFD-DEM model, the phase motion of discrete solids or particles is obtained by the Discrete Element Method (DEM) which applies Newton's laws of motion to every particle and the flow of continuum fluid is described by the local averaged Navier–Stokes equations that can be solved by the traditional Computational Fluid Dynamics (CFD). The model is first proposed by Tsuji et al. The interactions between the fluid phase and solids phase is better modeled according to Newton's third law.

Particle damping is the use of particles moving freely in a cavity to produce a damping effect.

<span class="mw-page-title-main">Extended discrete element method</span> Granular material interaction simulation technique

The extended discrete element method (XDEM) is a numerical technique that extends the dynamics of granular material or particles as described through the classical discrete element method (DEM) by additional properties such as the thermodynamic state, stress/strain or electro-magnetic field for each particle. Contrary to a continuum mechanics concept, the XDEM aims at resolving the particulate phase with its various processes attached to the particles. While the discrete element method predicts position and orientation in space and time for each particle, the extended discrete element method additionally estimates properties such as internal temperature and/or species distribution or mechanical impact with structures.

Model order reduction (MOR) is a technique for reducing the computational complexity of mathematical models in numerical simulations. As such it is closely related to the concept of metamodeling, with applications in all areas of mathematical modelling.

Goma is an open-source, parallel, and scalable multiphysics software package for modeling and simulation of real-life physical processes, with a basis in computational fluid dynamics for problems with evolving geometry. It solves problems in all branches of mechanics, including fluids, solids, and thermal analysis. Goma uses advanced numerical methods, focusing on the low-speed flow regime with coupled phenomena for manufacturing and performance applications. It also provides a flexible software development environment for specialty physics.

In granular mechanics, the μ(I) rheology is one model of the rheology of a granular flow.

<span class="mw-page-title-main">Numerical modeling (geology)</span> Technique to solve geological problems by computational simulation

In geology, numerical modeling is a widely applied technique to tackle complex geological problems by computational simulation of geological scenarios.

The Cohesion number (Coh) is a useful dimensionless number in particle technology by which the cohesivity of different powders can be compared. This is especially useful in DEM simulations of granular materials where scaling of the size and stiffness of the particles are inevitable due to the computationally demanding nature of the DEM modelling.

References

  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.
  5. Cundall, Peter. A.; Strack, Otto D. L. (1979). "Discrete numerical model for granular assemblies" (PDF). Géotechnique. 29 (1): 47–65. doi:10.1680/geot.1979.29.1.47.
  6. {{cite journal |last1=Williams |first1=J. R. |last2=Hocking |first2=G. |last3=Mustoe |first3=G. G. W. |title=The Theoretical Basis of the Discrete Element Method |journal=NUMETA 1985, Numerical Methods of Engineering, Theory and Applications |publisher=A.A. Balkema |location=Rotterdam |date=January 1985|url=https://docs.google.com/document/d/1ljujwjib2h2NwYksdh9wONZhEpNljGQdAmehXANFJw4%7C
  7. ( Williams, Hocking & Mustoe 1985 )
  8. Williams, J. R.; O'Connor, R. (December 1999). "Discrete element simulation and the contact problem". Archives of Computational Methods in Engineering. 6 (4): 279–304. CiteSeerX   10.1.1.49.9391 . doi:10.1007/BF02818917. S2CID   16642399.
  9. Munjiza, Ante (2004). The Combined Finite-Discrete Element Method. Chichester: Wiley. ISBN   978-0-470-84199-0.
  10. Alizadeh, Mohammadreza; Hassanpour, Ali; Pasha, Mehrdad; Ghadiri, Mojtaba; Bayly, Andrew (2017-09-01). "The effect of particle shape on predicted segregation in binary powder mixtures" (PDF). Powder Technology. 319: 313–322. doi:10.1016/j.powtec.2017.06.059. ISSN   0032-5910.
  11. Behjani, Mohammadreza Alizadeh; Motlagh, Yousef Ghaffari; Bayly, Andrew; Hassanpour, Ali (2019-11-07). "Assessment of blending performance of pharmaceutical powder mixtures in a continuous mixer using Discrete Element Method (DEM)". Powder Technology. 366: 73–81. doi:10.1016/j.powtec.2019.10.102. ISSN   0032-5910. S2CID   209718900. Archived from the original on 21 Feb 2020.
  12. Gan, Yixiang; Hernandez, Francisco; Hanaor, Dorian; Annabattula, Ratna; Kamlah, Marc; Pereslavtsev, Pavel (2014). "Thermal Discrete Element Analysis of EU Solid Breeder Blanket Subjected to Neutron Irradiation". Fusion Science and Technology. 66 (1): 83–90. arXiv: 1406.4199 . Bibcode:2014FuST...66...83G. doi:10.13182/FST13-727. S2CID   51903434.
  13. Tsory, Tal; Ben-Jacob, Nir; Brosh, Tamir; Levy, Avi (2013). "Thermal DEM–CFD modeling and simulation of heat transfer through packed bed". Powder Technology. 244: 52–60. doi:10.1016/j.powtec.2013.04.013.
  14. Lewis, R. W.; Gethin, D. T.; Yang, X. S.; Rowe, R. C. (2005). "A combined finite-discrete element method for simulating pharmaceutical powder tableting". International Journal for Numerical Methods in Engineering. 62 (7): 853. arXiv: 0706.4406 . Bibcode:2005IJNME..62..853L. doi:10.1002/nme.1287. S2CID   122962022.
  15. Gethin, D. T.; Yang, X. S.; Lewis, R. W. (2006). "A two dimensional combined discrete and finite element scheme for simulating the flow and compaction of systems comprising irregular particulates". Computer Methods in Applied Mechanics and Engineering. 195 (41–43): 5552. Bibcode:2006CMAME.195.5552G. doi:10.1016/j.cma.2005.10.025.
  16. Chen, Y.; May, I. M. (2009). "Reinforced concrete members under drop-weight impacts". Proceedings of the ICE - Structures and Buildings. 162: 45–56. doi:10.1680/stbu.2009.162.1.45.
  17. S.J. Antony (2007). "Link between single-particle properties and macroscopic properties in particulate assemblies: role of structures within structures". Philosophical Transactions of the Royal Society of London, Series: A. 356: 2879-2891.
  18. S. J. Antony, D. Chapman, J. Sujatha and T. Barakat (2015). "Interplay between the inclusions of different sizes and their proximity to the wall boundaries on the nature of their stress distribution within the inclusions inside particulate packing". Powder Technology. 286: 286, 98-106.{{cite journal}}: CS1 maint: multiple names: authors list (link)
  19. He, Yi; Bayly, Andrew E.; Hassanpour, Ali; Muller, Frans; Wu, Ke; Yang, Dongmin (2018-10-01). "A GPU-based coupled SPH-DEM method for particle-fluid flow with free surfaces". Powder Technology. 338: 548–562. doi: 10.1016/j.powtec.2018.07.043 . ISSN   0032-5910.

Bibliography

Book

Periodical

Proceedings