Search Documents

Discrete Element Method - Computational Fluid Dynamics Modelling of Flow Behaviour in a Stirred Mill - Effect of Operation Conditions

Sep 13, 2010

By J. T. Kalala

Energy consumption and wear of grinding media and liners are the major cost factors in milling plants. Since the beginning of the last century, much effort has been directed at the prediction of the power draw as affected by the mill operating conditions. This paper explores the prediction of the liner wear using the discrete element method (DEM). A mathematical model taking into account adhesion, abrasion and impact wear which occur in dry mills is derived in order to predict the loss of material on liners. The mathematical model is implemented in the DEM code using an objective function in order to produce realistic simulated profiles. Our DEM simulated predictions are compared to published experimental and industrial data. Keywords: modelling; discrete element method; liner wear

Jan 1, 2004

By Barada K. Mishra

The potential merit of the discrete element method (DEM) as a numerical tool for analyzing process engineering problems involving particulate material is discussed. In this method, the dynamic equilibrium equations for each distinct and deformable body are solved. Therefore, it readily lends itself to modeling large numbers of discrete, interacting bodies as in mineral processing operations where a charge of particulate material is processed. Here some applications of DEM are reported: motion of particles in the ball mill, mixer, and jig. Although DEM has made very little headway in mineral processing research, it is shown here that a host of other areas of industrial interest in the context of flow of particulate materials can also be made amenable to analysis by this method.

Jan 1, 1993

By Qingwen Shi, Brijes Mishra

Laminated rock can delaminate and fail under certain stress states. Delamination-caused roof fall and cutter roof failure is a common failure mechanism in roadways in Pittsburgh coal seam due to the existence of laminated shale roof. Two-dimensional numerical models were developed at laboratory scale based on the bonded particle method in order to simulate the delamination with particle flow code. In these models, an assembly of bonded particles and parallel weak planes represents the laminated rock. This model calibrated the micro-parameters of the bonded particlematerial and weak planes with laboratory data. The mechanism of fracturing and delamination in a laminated specimen was investigated using unconfined and confined compression test and unsupported roof compression test. The numerical results show that, under a constant horizontal stress, delamination occurswhen the stress is much lower than its compressive strength. The confined compressive test indicates that delamination is restrained by applying confining stress since delamination-caused failure does not occur during any stage of the confined compressive test. Unsupported roof compression tests showed propagation of cutter roof failure with initiation of cracks and delamination at the intersections of the roof and ribs. An outer failed zone and inner failed zone form separately. The separated zones connect and form an opening in the roof. The presented numerical results provide additional insight into process of delamination.

Sep 14, 2020

By Dale S. Preece

A Discrete element computer program named DMC @istinct Motion Code) has been developed for modeling rock blasting. This program employs explicit time integration and uses spherical or cylindrical elements which are represented as circles in 2-D. DMC calculations have been compared with measurements on bench blasts in the field with relatively good comparison. Structural rock mass characteristics have a significant impact on any blast and DMC has not, until now, included these effects. This paper discusses a recently added DMC capability for treating joints and bedding planes in bench blast simulations.

Jan 1, 1995

By Dale S. Preece

A Discrete element computer program named DMC (Distinct Motion Code) has been developed for modeling rock blasting. This program employs explicit time integration and uses spherical or cylindrical elements which are represented as circles in 2-D. DMC calculations have been compared with measurements on bench blasts in the field with relatively good comparison. Structural geology and joint sets have a significant impact on any blast and DMC has not, until now, included these effects. This paper discusses a recently added DMC capability for treating joints and bedding planes in bench blast simulations.

Jan 1, 1995

By Graham G. W. Mustoe

Energy losses are important factors in a conveyor's performance and design. A significant energy loss due to belt sag within a conveyor system occurs when shearing motion between the bulk material particles occurs as the material moves along the belt. An initial analysis of these phenomena was performed using the discrete element method, which is a general-purpose numerical modeling technique for analyzing the mechanical behavior of particulate materials. The software used in this work was Overland Conveyor Company's, Chute Analyst Professional, a discrete element computer code designed specifically for analyzing bulk materials transportation and handling problems.

Jan 1, 2008

By A. Fakhimi

The effect of porosity on mechanical behavior of rock has been studied in the past using both experimental and numerical techniques. One important aspect of this study that has not received much attention is the role of void size and distribution on the mechanical behavior of rock. In this study, the CA2 computer program that is a hybrid discrete-finite element program for two-dimensional simulation of geomaterials is used. Numerical specimens with different macro-void sizes are loaded uniaxially. The Young?s modulus, Poisson?s ratio, uniaxial strength and crack initiation stress are studied. The numerical results are compared with the database of experimental tests on Plaster of Paris specimens containing Styrofoam spheres with different void sizes and porosities. It is shown that the macro void spatial distribution is partly responsible for the scatter of the experimental data. It is also shown that not-only the porosity but also the pore size can affect the deformability, crack initiation stress and rock strength.

May 1, 2009

By Ki-Bok Min, Bona Park

"In spite of the applications of DEM models to emulate the elastic modulus, Poisson’s ratio and strengths of isotropic rock in laboratory scale, extensions to anisotropic rock and larger engineering scale still have many challenges to overcome such as robust verification and validation against observed field case. In this study, a discrete element modeling was performed in order to represent a transversely isotropic rock such as shale, and the upscaled model was applied to foundation and borehole problems using PFC2D. The applicability of the upscaled rock mass model was verified by an analytical solution. The DEM modeling of stress redistribution induced by an applied surface line load compares well with Boussinesq’s isotropic solution and anisotropic solution given by Bray. The borehole model with a circular opening under biaxial stress was compared with both isotropic and anisotropic solutions resulting in good agreement. The abilities of DEM model to investigate the micromechanisms accounting for complex macroscopic behaviors and to predict the macroscopic behavior are discussed.The results demonstrate that DEM modeling has the potential to be applied to anisotropic rock mechanics problems on larger engineering scales such as wellbore stability problems especially in shale formations. Furthermore, the effect of anisotropy modeled through this study may be able to explain failure behavior of anisotropic rock encountered frequently in real cases.INTRODUCTIONAnisotropy is prevalent in various types of rocks, especially in shale that contain a set of parallel bedding planes. These weak planes in various forms are sources of large deformation and low shear strength and contribute significantly to the anisotropic behavior of rocks (Amadei, 1996; Cho et al., 2012; Kim et al., 2012). As elastic and strength behavior of anisotropic rock is substantially different from that of isotropic rock, proper consideration of anisotropy is of great importance. There is extensive of experimental evidence showing the anisotropy of shale. Niandou et al. (Niandou et al., 1997) analyzed the plastic deformation and failure mechanisms of Tournemire shale were highly dependent on the orientation of bedding planes. Cho et al. (Cho et al., 2012) presented various deformations and strengths of Boryeong shale obtained from laboratory experiments. Fjaer and Nes (Fjaer and Nes, 2013) performed uniaxial and triaxial compressive strength tests using Mancos shale, and they observed that the shear sliding along the planes of weakness before the occurrence of a major failure can reduce the strength of anisotropic rock. The recent research has demonstrated that the anisotropy of rock is an important consideration in rock engineering."

Jan 1, 2015

By D Johansson, C Yi

Blasting techniques are widely used to fragment rock masses into smaller pieces. Numerical modelling is an efficient method employed by many researchers in the blasting field. It is difficult for a conventional continuum-based approach such as the finite element method (FEM) to model the rock fragmentation by blasting and the expansion work on the rock by explosive and its detonation products. In this paper, the particle blast method (PBM) was employed to model the behaviour of the detonation and a bonded particle model (BPM) was used to model the brittle material to be blasted. A mortar cylinder with a centrally placed hole for the explosive was modelled and the results were compared to the experimental data. The blast process from crack initiation to fragment formation was analysed. The influence of coupling ratio on fragmentation was also investigated.CITATION:Yi, C and Johansson, D, 2015. Discrete element modelling of blast fragmentation of a mortar cylinder, in Proceedings 11th International Symposium on Rock Fragmentation by Blasting, pp 793–798 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Aug 24, 2015

By B. K. Mishra, R. K. Rajamani

Discrete element is a numerical technique for describing the collision dynamics and motion of an assembly of particles. In particular, calculations are done on full three-dimensional space with a force-contact law for collisions and Newton’s law for subsequent motion. It may be said that there is no empiricism in this approach. The paper describes the full scale simulation of large diameter semiautogenous grinding mill of size 40X24 ft. It is shown that the simulation predicts the power draft of plant size mill extremely accurately. Also these simulations have been adapted uniquely to assess the motion of pebbles and slurry in pulp lifters. The passage of slurry from pulp lifters into trunnion is extremely important since flow-back and carryover must be avoided at all cost. The dynamics of slurry flow through pulp lifters is also described in this paper.

Jan 1, 2003

By Jamie King

This paper describes the use of dynamic scheduling and discrete event control to provide coordinated control of multiple automated vehicles. In the past decade there have been significant advances in the field of robotics for underground mining. As this trend progresses toward a realization of a completely automated underground mine the need for a central control system, capable of providing coordination for multiple vehicles in a semistructured environment, becomes more evident. In the past year general research into multi-robot coordination has been applied to the underground mining scenario and demonstrated in a scaled model underground mine in a laboratory setting. This system shows significant potential to be adopted into a full-scale implementation for mine wide automation.

May 1, 2001

By T. Norvell, R. Gosine, S. O’Young

"This paper describes the use of dynamic scheduling and discrete event control to provide coordinated control of multiple automated vehicles. This work has been carried out within ongoing research into multi-robot coordination developed around the underground mining scenario. Demonstrations using scaled laboratory models of underground mines illustrate the potential for this technology for full-scale implementation for mine-wide automation. IntroductionIn the past decade, the mining industry has been improving the efficiency of underground mining operations. Robotics and automation have demonstrated the potential to further reduce costs. Tele-operation of underground vehicles from above ground and autonomous navigation have reduced travel delays and allowed for single operator control of multiple vehicles (Bloomquist and Hinton, 2001). Other recent projects include the automation of the digging process as well as the drilling and blasting operations. Infrastructureless guidance systems are also being developed (Bloomquist and Hinton, 2001).Mine-wide automation requires significant developments at two levels: (1) the development of mining robots capable of performing tasks either autonomously or semiautonomously; and (2) the development of methods for automating operations at a systems level. This paper focusses on the latter, and describes a general framework for the coordination of multiple mobile robots in an unstructured or semi-structured environment."

Jan 1, 2003

By Alessandro Navarra

Discrete event simulation provides a computational framework in which the use of smelting equipment can be analyzed and optimized, leading to the development of alternative modes of operation. Much of the complexity of copper and nickel smelters is centered on the converter aisle. These converters operate in parallel, produce batches, and set the tempo for the downstream operations. In practice, smelter staff schedule daily operations to coordinate the smelting furnaces and converters with downstream operations. Thus, a realistic simulation of smelter operations, over hundreds of days, should incorporate converter scheduling algorithms as presented in this paper.

Jan 1, 2014

By Sahin Furkan Sahiner, Onur Golbasi

Effective maintenance planning is vital for coal mining operations, where haul truck fleets directly determine production continuity. This study presents a discrete event simulation (DES) model for evaluating maintenance strategies under multiple failure modes and their interdependencies. The model integrates corrective, periodic preventive, and opportunistic maintenance, dynamically updating equipment lifetimes through probabilistic time-betweenfailure and time-to-repair distributions. Applied to a fivetruck coal haulage fleet over a one-year horizon, results show that periodic preventive maintenance delivers the highest availability (76%) under unconstrained conditions, slightly outperforming corrective-only maintenance (75%). Opportunistic and integrated approaches performed worse (68% and 69%), reflecting redundancy from overlapping interventions. These findings suggest that structured preventive programs provide the most reliable outcomes when costs and resources are not limiting. The DES framework serves as a baseline for future models incorporating cost, spare parts logistics, and workforce availability, supporting robust maintenance planning in coal mining operations.

Feb 22, 2026

By Kristopher John Shelswell

The work in this paper describes the theory, development and analysis of a computer model simulating a mine operation moving material using concurrent decline ramp truck haulage and skipping practices. This simulation was designed specifically to assess the performance of the truck haulage component and to expand on the scope and capabilities of traditional TKM fleet calculations. Haulage fleet simulation results were similar to TKM calculations for the lower to moderate production rate cases. However, comparisons of ramp-up and high production rate scenarios indicated that the simulation model was significantly more sensitive to operating practices and factors not accounted for by calculations derived when using a single TKM figure. Given the proper calibration of inputs and logic, the simulation model was additionally capable of providing information associated with general mine operations efficiency. Further benefits of the simulation model over a TKM approach included logistical assessments such as production schedule feasibility, productivity trade-off quantifications for proposed operating practices, and the validation of mine resource utilization scheduling. The simulation model was able to expand on traditional methods to approach planning and logistics by considering factors such as truck queuing, real-time traffic interactions, random events, and competition for resources on mine productivity.

May 1, 2013

By Greg Yuriy

This paper discusses research carried out by the Laurentian University Mining Automation Laboratory (LUMAL) in the field of discrete-event simulation of mine equipment systems and reliability assessment of these machinery using genetic algorithms. The AutoMod simulation package, from Brooks-PRI Automation Inc., is used in this research. It is a 3-D discrete-event simulation/animation software tool for creating realistic depictions of manufacturing and materials handling systems. Discrete-event simulation involves modeling of real-life systems based on variables and events which occur or change at distinct points in time. In our case, a two level, multi-drift sublevel stoping model was created in AutoMod for the simulation of the equipment fleet employed in the development cycle (drilling, explosives loading, blasting, materials hauling, ground support) of a typical underground mine in the nickel mining region of Sudbury in Ontario, Canada. The equipment primarily considered in the simulation is Load-Haul-Dump (LHD) vehicles (scooptrams), drills, explosives-loaders and roof bolters. The purpose of the simulation study is to assess the impact of equipment failures to mine?s production throughput. In more detail, using Mean Time Between Failures (MTBF) data of scooptrams derived from a reliability assessment model based on genetic algorithms, LHD failures are incorporated into the simulation model. These failures cause production interruptions to the development cycle and consequently loss of hauling capacity. This paper will indicate the significance of equipment failures to a mine?s production and suggest a methodology based on the combination of discrete-event simulation and a predictive reliability assessment model using genetic algorithms for the study and analysis of mine equipment systems.

May 1, 2003

By G. Goodall, M. Muinonen, R. Pandher, R. Jabbar, S. J. Muinonen

"The implementation of a new and highly efficient converting practice at the Glencore Sudbury Integrated Nickel Operations Smelter has opened the possibility to further capacity increases by using larger ladles. To quantify the impact of increasing ladle sizes from 200 to 300 ft3, XPS constructed a detailed discrete event simulation using the commercial software package Rockwell Arena®. The challenges associated with implementing a model on an intricate and tightly scheduled process, and the solutions needed to derive meaningful results will be presented. The aisle capacity is predicted to increase by approximately 17 %, with a reduction in furnace tapping demand of 8 %. Arena modeling capabilities and the key factors contributing to the overall capacity increase at the Sudbury INO Smelter will be reviewed and presented.INTRODUCTION The Glencore Sudbury Integrated Nickel Operations (INO) comprises two Ni-Cu-PGM mines, Strathcona Mill, and the Sudbury Smelter. The Smelter has operated continuously since 1929 using different processes to produce nickel-copper Bessemer matte. Feeds to the Sudbury INO Smelter are a blend of concentrate from Strathcona Mill, Raglan concentrate produced at the mine-mill complex in Raglan, Quebec, as well as third party concentrate feeds shipped to Sudbury from around the world. In 2016, the Sudbury INO Smelter produced 75,900 tonnes of nickel and 20,000 tonnes of copper in Bessemer matte. The product matte is shipped to Glencore’s Nikkelverk refinery in Kristiansand, Norway where it is refined to produce nickel and copper cathode. A schematic of the Sudbury INO Smelter’s flowsheet is shown in Figure 1."

Jan 1, 2017

By Bill Dershowitz

Consider a deep gold orebody with grade of 4 g/t (0.11 oz/st). At the time of this writing, with gold prices at historic highs, it could be a good, but not great, prospect. But if it is a vein deposit, and those veins are 200 g/t (5.8 oz/st), it would be significantly more attractive if there was a way to produce directly from those veins, leaving the rest of the rock in place. [Figure 1 The red around each of the four injector wells shows areas of high pressure, where the lixiviant fluid is pumped into the rock mass, with orange and yellow indicating dissipated pressure. The black lines indicate some of the fracture pathways through which the lixiviant, with absorbed metal, flows toward the central recovery well where low pressure is indicated in blue. Lixiviant is pumped out through the recovery well. Image courtesy of Golder Associates.] There is a way to do this, in situ leach (ISL) or solution mining. It has been a viable commercial technology for more than five decades.

Jan 1, 2011

By A L. de Ambrosis, F M. Weir, M Fowler

Traditional kinematic and statistical analyses of geological structure for geotechnical engineering applications, while considered industry standard practice, have some limitations and unrealistic assumptions. An improved approach involves a discrete fracture network (DFN), a numerical model that explicitly represents how defects may be spatially distributed through the rock mass. The DFN technique utilises the structural data in stochastic analyses to give a range of possible models used for visualisation and stability analyses.The early development of the DFN approach was largely associated with modelling fractured hydrocarbon reservoirs and nuclear waste isolation programs. With recent improvements to computational power it has become recognised internationally as a valuable tool for modelling some geomechanical problems encountered when designing excavations in fractured rock masses. For tunnelling or slope stability projects DFN modelling provides visualisation of a rock mass and identification of likely failure mechanisms.This paper presents a suite of DFN modelling undertaken for the design of a tunnel within the Ashfield Shale unit of the Sydney Basin, Australia. The necessary input parameters and derivation of these, along with the importance of both an accurate engineering geology model and validation of the synthetic fracture network are demonstrated. Multiple realisations of the DFN were analysed to provide an indication of the number and volume of unstable blocks for the proposed tunnel geometry, enhancing support design. One fracture network realisation was sliced along multiple sections for use in finite element modelling.CITATION:Weir, F M, de Ambrosis, A L and Fowler, M, 2014. Discrete fracture network modelling for tunnel design in the Sydney Basin, in Proceedings 15th Australasian Tunnelling Conference 2014 , pp 219–226 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Sep 17, 2014