Selection of an appropriate time integration scheme for the discrete element method (DEM)

2008 ◽  
Vol 32 (10) ◽  
pp. 2263-2279 ◽  
Author(s):  
H. Kruggel-Emden ◽  
M. Sturm ◽  
S. Wirtz ◽  
V. Scherer
1981 ◽  
Vol 103 (4) ◽  
pp. 657-664 ◽  
Author(s):  
H. U. Akay ◽  
A. Ecer

Analysis of transonic flow through a cascade of airfoils is investigated using the finite element method. Development of a computational grid suitable for complex flow structures and different types of boundary conditions is presented. An efficient pseudo-time integration scheme is developed for the solution of equations. Modeling of the shock and the convergence characteristics of the developed scheme are discussed. Numerical results include a 45 deg staggered cascade of NACA 0012 airfoils with inlet flow Mach number of 0.8 and angles of attack 1, 0, and −1 deg.


2019 ◽  
Vol 6 (4) ◽  
pp. 545-559 ◽  
Author(s):  
Joaquín Irazábal ◽  
Fernando Salazar ◽  
Miquel Santasusana ◽  
Eugenio Oñate

Author(s):  
G. F. Mathews ◽  
R. L. Mullen ◽  
D. C. Rizos

This paper presents the development of a semi-implicit time integration scheme, originally developed for structural dynamics in the 1970’s, and its implementation for use in Discrete Element Methods (DEM) for rigid particle interaction, and interaction of elastic bodies that are modeled as a cluster of rigid interconnected particles. The method is developed in view of ballast modeling that accounts for the flexibility of aggregates and the arbitrary shape and size of granules. The proposed scheme does not require any matrix inversions and is expressed in an incremental form making it appropriate for non-linear problems. The proposed method focuses on improving the efficiency, stability and accuracy of the solutions, as compared to current practice. A critical discussion of the findings of the studies is presented. Extended verification and assessment studies demonstrate that the proposed algorithm is unconditionally stable and accurate even for large time step sizes. It is demonstrated that the proposed method is at least as computationally efficient as the Central Difference Method. Guidelines for the implementation of the method to ballast modeling are discussed.


Minerals ◽  
2021 ◽  
Vol 11 (8) ◽  
pp. 878
Author(s):  
Manuel Moncada ◽  
Patricio Toledo ◽  
Fernando Betancourt ◽  
Cristian G. Rodríguez

Comminution by gyratory crusher is the first stage in the size reduction operation in mineral processing. In the copper industry, these machines are widely utilized, and their reliability has become a relevant aspect. In order to optimize the design and to improve the availability of gyratory crushers, it is necessary to calculate their power and torque accurately. The discrete element method (DEM) has been commonly used in several mining applications and is a powerful tool to predict the necessary power required in the operation of mining machines. In this paper, a DEM model was applied to a copper mining gyratory crusher to perform a comprehensive analysis of the loads in the mantle, the crushing torque, and crushing power. A novel polar representation of the radial forces is proposed that may help designers, engineers, and operators to recognize the distribution of force loads on the mantle in an easier and intuitive way. Simulations with different operational conditions are presented and validated through a comparison with nominal data. A calculation procedure for the crushing power of crushers is presented, and recommendations for the selection of the minimum resolved particle size are given.


2004 ◽  
Vol 10 (1) ◽  
pp. 3-14 ◽  
Author(s):  
Robertas Balevičius ◽  
Algis Džiugys

Application of discrete element method (DEM) to keel penetration in granular media is investigated. The basic relations for visco‐elastic granular media composed of spherical particles are presented, together with 5th order Gear predictor‐corrector scheme for time‐integration. The background version of DEM and numerical time integration algorithm are developed and implemented into DEMMAT code. The implementation of time‐integration algorithm is verified by simple tests concerning particle‐particle, particle‐wall interactions, for which analytical expressions exist. By limiting the size of the media domain, the three‐dimensional problem is reduced to particular case presented as two‐dimensional domain of spherical particles. The variation of keel reaction and distribution of the particle forces due to different material properties are investigated


Sign in / Sign up

Export Citation Format

Share Document