Determination of the normal spring stiffness coefficient in the linear spring–dashpot contact model of discrete element method

Navarro, Hélio Aparecido; Braun, Meire Pereira de Souza

doi:10.1016/j.powtec.2013.05.049

Cited by 60 publications

(20 citation statements)

References 67 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…To correctly capture this behaviour, Cummins et al [34] describe how to modify the contact duration calculation by using the dashpot coefficient as described by equation (A.8). For further discussion on the contact forces, additional contact laws, and specific DEM-related topics such as contact duration, etc., see [21,34,35,50,51].…”

Section: Abbreviationsmentioning

confidence: 99%

Partitioned Strong Coupling of Discrete Elements with Large Deformation Structural Finite Elements to Model Impact on Highly Flexible Tension Structures

Sautter

Teschemacher

Celigueta

et al. 2020

Advances in Civil Engineering

View full text Add to dashboard Cite

This article presents a staggered approach to couple the interaction of very flexible tension structures with large deformations, described with the finite element method (FEM), and objects undergoing large, complex, and arbitrary motions discretized with particle methods, in this case the discrete element method (DEM). The quantitative solution quality and convergence rate of this partitioned approach is highly time step dependent. Thus, the strong coupling approach is presented here, where the convergence is achieved in an iterative manner within each time step. This approach helps increase the time step size significantly, decreases the overall computational costs, and improves the numerical stability. Moreover, the proposed algorithm enables the application of two independent, standalone codes for simulating DEM and structural FEM as blackbox solvers. Systematic evaluations of the newly proposed iterative coupling scheme with respect to accuracy, robustness, and efficiency as well as cross comparisons between strong and weak FEM-DEM coupling approaches are performed. Additionally, the approach is validated against the rest position of an impacting object, and further examples with objects impacting highly flexible protection structures are presented. Here, the protection nets are described with nonlinear structural finite elements and the impacting objects as DEM elements. To allow the interested reader to independently reproduce the results, detailed code and algorithm descriptions are included in the appendix.

show abstract

Section: Abbreviationsmentioning

confidence: 99%

Partitioned Strong Coupling of Discrete Elements with Large Deformation Structural Finite Elements to Model Impact on Highly Flexible Tension Structures

Sautter

Teschemacher

Celigueta

et al. 2020

Advances in Civil Engineering

View full text Add to dashboard Cite

show abstract

“…Particle-particle and particle-wall collisional forces are handled through a soft-sphere spring-damper collision model, similar to the one originally proposed by Cundall and Strack (1979). In this work, the normal component of the model used by Navarro and De Souza Braun (2013) has been used. The collision force from the j th particle on the i th particle is composed of a conservative spring force (F S c;ij ) and dissipative damping force (F D c;ij ) by…”

Section: Lagrangian Descriptionmentioning

confidence: 99%

“…Þn ij is the normal relative contact velocity and h is the normal damping coefficient. Navarro and De Souza Braun (2013) rigorously derived the damping coefficient and the time in collision (t c ) as…”

Section: Lagrangian Descriptionmentioning

confidence: 99%

A scalable Euler–Lagrange approach for multiphase flow simulation on spectral elements

Zwick¹,

Balachandar

2019

The International Journal of High Performance Computing Applica

View full text Add to dashboard Cite

Multiphase flow can be difficult to simulate with high accuracy due to the wide range of scales associated with various multiphase phenomena. These scales may range from the size of individual particles to the entire domain of interest. Traditionally, large scale systems can only be simulated using averaging approaches that filter out the locations of individual particles. In this work, the Euler–Lagrange method is used to simulate large-scale dense particle systems in which each individual particle is tracked. In order to accomplish this, the highly scalable spectral element code nek5000 has been extended to handle the multiple levels of multiphase coupling in these systems. These levels include what has been called one-, two-, and four-way coupling. Here, each level has been separated to detail the computational impact of each stage. A binned ghost particle algorithm has also been developed to efficiently handle the challenges of two- and four-way coupling in a parallel processing context. The algorithms and their implementations are then shown to scale to 65,536 Message Passing Interface (MPI) ranks in both the strong and weak limits. After this, validation is performed through simulation of a small-scale fluidized bed. Lastly, a large-scale fluidized bed is simulated with 65,536 MPI ranks and is able to capture the unique physics of the onset of fluidization.

show abstract

“…Various methods for determining the value of this parameter have been proposed, most of which equate some aspect of the contact to a Hertzian contact. For example, the anticipated maximum contact overlap, the anticipated maximum elastic strain energy or the binary collision duration have all been used to determine the linear spring stiffness. In the current work,

k_{linear}

is chosen such that the overlap caused by the weight of N particles resting on top of the bottom particle is equal to the overlap that would occur if the linear spring was replaced by a Hertzian spring, assuming a circular contact area (Eq.…”

Section: Simulation Descriptionmentioning

confidence: 99%

Influence of normal contact force model on simulations of spherocylindrical particles

et al. 2018

View full text Add to dashboard Cite

How the choice of elastic normal contact force model affects predictions from discrete element method simulations of spherocylindrical particles is investigated in this article. Three force models were investigated: (1) a Hertzian force model (HFM) which assumes a circular contact area; (2) a linear force model (LFM) with a constant stiffness; and (3) a modified HFM (MFM) that accounts for various contact areas and contact transitions. With the MFM, transitions between contact area types must be accounted for otherwise discontinuities in the contact force can occur. It is found that simple force models (HFM, LFM) can be substituted for more accurate force models if only force data and bulk properties are of interest. However, if more detailed contact information, such as contact area, contact overlap, contact duration, or collision frequency, are needed, for example, in population balance models and transient liquid bridge modeling, then a more accurate force model should be used. The red band is the range of angles reported in the literature by B€ orzs€ onyi et al., 40 Guo et al., 41 and Hua et al. 4 [Color figure can be viewed at wileyonlinelibrary.com]

show abstract

Determination of the normal spring stiffness coefficient in the linear spring–dashpot contact model of discrete element method

Cited by 60 publications

References 67 publications

Partitioned Strong Coupling of Discrete Elements with Large Deformation Structural Finite Elements to Model Impact on Highly Flexible Tension Structures

Partitioned Strong Coupling of Discrete Elements with Large Deformation Structural Finite Elements to Model Impact on Highly Flexible Tension Structures

A scalable Euler–Lagrange approach for multiphase flow simulation on spectral elements

Influence of normal contact force model on simulations of spherocylindrical particles

Contact Info

Product

Resources

About