Document Type: Original Research Paper


Faculty of Mining, Petroleum & Geophysics Engineering, Shahrood University of Technology, Shahrood, Iran


The number of lifters of mill shell liners, mill rotation speed, and filling percentage of grinding media are three of the most important parameters influencing the charge behavior and the trajectory of ball motion inside the SAG mills, and consequently, their performance. In this paper, the milling operation of pilot-scale SAG mills using the discrete element method (DEM) is investigated. First, a pilot-scale SAG mill with dimensions of 3.0 m × 1.5 m with no lifter is simulated. Then by adding, respectively, one, two, four, eight, sixteen, and thirty-two rectangle lifter(s), six other independent simulations are performed. The effects of the number of lifters on the two new parameters introduced by the authors, i.e. ‘head height’ and ‘impact zone length’ as well as on creation of cascading, cataracting, and centrifuging motions for balls at two different mill speeds, i.e. 70% and 80% of its critical speed (NC), are evaluated. Also in order to validate the simulation results, a laboratory-scale SAG mill is simulated. The results obtained indicate that the optimum number of lifters for pilot-scale SAG mills is between 16 and 32 lifters with medium thickness. Liners with the number of lifters in this range require less mill speed to create cataract motions. However, liners with the number of lifters less than this range require a higher mill speed. Also liners with the number of lifters beyond this range require less mill speed, and can cause centrifugal motions in the balls. Comparison of the simulations related to the laboratory-scale SAG mill with experimental results demonstrates a good agreement, which validates the DEM simulations and the software used.


Main Subjects

[1]. Rosales-Marín, G., Andrade, J., Alvarado, G., Delgadillo, J.A. and Tuzcu, E.T. (2019). Study of lifter wear and breakage rates for different lifter geometries in tumbling mill: Experimental and simulation analysis using population balance model. Minerals Engineering. 141: 105857.‏

[2]. Mishra, B.K. and Rajamani, R.K. (1992). The discrete element method for the simulation of ball mills. Applied Mathematical Modelling. 16 (11): 598-604.‏

[3]. Powell, M.S. and McBride, A.T. (2004). A three-dimensional analysis of media motion and grinding regions in mills. Minerals Engineering. 17 (11-12): 1099-1109.‏

[4]. Cleary, P.W. (1998). Predicting charge motion, power draw, segregation and wear in ball mills using discrete element methods. Minerals Engineering. 11 (11): 1061-1080.‏

[5]. Cleary, P.W. (2001). Charge behaviour and power consumption in ball mills: sensitivity to mill operating conditions, liner geometry and charge composition. International journal of mineral processing. 63 (2): 79-114.‏

[6]. Kalala, J.T., Bwalya, M. and Moys, M.H. (2005). Discrete element method (DEM) modelling of evolving mill liner profiles due to wear. Part II. Industrial case study. Minerals Engineering. 18 (15): 1392-1397.‏

[7]. Kalala, J.T., Bwalya, M. and Moys, M.H. (2005). Discrete element method (DEM) modelling of evolving mill liner profiles due to wear. Part II. Industrial case study. Minerals Engineering. 18 (15): 1392-1397.‏

[8]. Banisi, S. and Hadizadeh, M. (2007). 3-D liner wear profile measurement and analysis in industrial SAG mills. Minerals Engineering. 20 (2): 132-139.‏

[9]. Yahyaei, M., Banisi, S. and Hadizadeh, M. (2009). Modification of SAG mill liner shape based on 3-D liner wear profile measurements. International Journal of Mineral Processing. 91 (3-4): 111-115.‏

[10]. Cleary, P.W. (2001). Recent advances in DEM modelling of tumbling mills. Minerals Engineering. 14 (10): 1295-1319.‏

[11]. El-Rahman, M.A., Mishra, B.K. and Rajamani, R.K. (2001). Industrial tumbling mill power prediction using the discrete element method. Minerals engineering. 14 (10): 1321-1328.‏

[12]. Cleary, P.W., Morrisson, R. and Morrell, S. (2003). Comparison of DEM and experiment for a scale model SAG mill. International Journal of Mineral Processing. 68 (1-4): 129-165.‏

[13]. Mishra, B.K. (2003). A review of computer simulation of tumbling mills by the discrete element method: part I—contact mechanics. International journal of mineral processing. 71 (1-4): 73-93.‏

[14]. Morrison, R.D. and Cleary, P.W. (2004). Using DEM to model ore breakage within a pilot scale SAG mill. Minerals Engineering. 17 (11-12): 1117-1124.‏

[15]. Djordjevic, N., Shi, F.N. and Morrison, R. (2004). Determination of lifter design, speed and filling effects in AG mills by 3D DEM. Minerals Engineering. 17 (11-12): 1135-1142.‏

[16]. Djordjevic, N., Morrison, R., Loveday, B. and Cleary, P. (2006). Modelling comminution patterns within a pilot scale AG/SAG mill. Minerals Engineering. 19 (15): 1505-1516.‏

[17]. Maleki-Moghaddam, M., Ghasemi, A.R., Yahyaei, M. and Banisi, S. (2015). The impact of end-wall effect on the charge trajectory in tumbling model mills. International Journal of Mineral Processing. 144: 75-80.‏

[18]. Owen, P. and Cleary, P.W. (2015). The relationship between charge shape characteristics and fill level and lifter height for a SAG mill. Minerals Engineering. 83: 19-32.‏

[19]. Weerasekara, N.S., Liu, L.X. and Powell, M.S. (2016). Estimating energy in grinding using DEM modelling. Minerals Engineering. 85: 23-33.‏

[20]. Cleary, P.W. and Owen, P. (2018). Development of models relating charge shape and power draw to SAG mill operating parameters and their use in devising mill operating strategies to account for liner wear. Minerals Engineering. 117: 42-62.‏

[21]. Xu, L., Luo, K. and Zhao, Y. (2018). Numerical prediction of wear in SAG mills based on DEM simulations. Powder Technology, 329, 353-363.‏

[22]. Cleary, P.W. and Owen, P. (2019). Effect of operating condition changes on the collisional environment in a SAG mill. Minerals Engineering. 132: 297-315.‏

[23]. Hasankhoei, A.R., Maleki-Moghaddam, M., Haji-Zadeh, A., Barzgar, M.E. and Banisi, S. (2019). On dry SAG mills end liners: Physical modeling, DEM-based characterization and industrial outcomes of a new design. Minerals Engineering. 141: 105835.‏

[24]. Bian, X., Wang, G., Wang, H., Wang, S. and Lv, W. (2017). Effect of lifters and mill speed on particle behaviour, torque, and power consumption of a tumbling ball mill: Experimental study and DEM simulation. Minerals Engineering. 105: 22-35.‏

[25]. Pedrayes, F., Norniella, J.G., Melero, M.G., Menéndez-Aguado, J.M. and del Coz-Diaz, J.J. (2018). Frequency domain characterization of torque in tumbling ball mills using DEM modelling: Application to filling level monitoring. Powder Technology, 323, 433-444.‏

[26]. Chegeni, M.J. and Kolahi, S. An investigation on effect of shell liner type on performance of industrial SAG mills using DEM.‏

[27]. Cleary, P.W. (2004). Large scale industrial DEM modelling. Engineering Computations.‏

[28]. Cleary, P.W., Sinnott, M.D. and Morrison, R.D. (2008). DEM prediction of particle flows in grinding processes. International Journal for numerical methods in fluids. 58 (3): 319-353.‏

[29]. Cleary, P.W. (2009). Ball motion, axial segregation and power consumption in a full scale two chamber cement mill. Minerals Engineering. 22 (9-10): 809-820.‏

[30]. Kozicki, J. and Donze, F.V. (2009). YADE‐OPEN DEM: An open‐source software using a discrete element method to simulate granular material. Engineering Computations.‏

[31]. Chen, J., Huang, B., Chen, F. and Shu, X. (2012). Application of discrete element method to Superpave gyratory compaction. Road materials and pavement design. 13 (3): 480-500.‏

[32]. Zhang, L., Quigley, S.F. and Chan, A.H. (2013). A fast scalable implementation of the two-dimensional triangular Discrete Element Method on a GPU platform. Advances in Engineering Software. 60: 70-80.‏

[33]. Kruggel-Emden, H., Sturm, M., Wirtz, S. and Scherer, V. (2008). Selection of an appropriate time integration scheme for the discrete element method (DEM). Computers & Chemical Engineering. 32 (10): 2263-2279.‏

[34]. Nassauer, B., Liedke, T. and Kuna, M. (2013). Polyhedral particles for the discrete element method. Granular matter. 15 (1): 85-93.‏

[35]. Raji, A.O. and Favier, J.F. (2004). Model for the deformation in agricultural and food particulate materials under bulk compressive loading using discrete element method. I: Theory, model development and validation. Journal of food engineering. 64 (3): 359-371.‏

[36]. Nassauer, B. and Kuna, M. (2013). Contact forces of polyhedral particles in discrete element method. Granular Matter. 15 (3): 349-355.‏

[37]. Balevičius, R., Džiugys, A., Kačianauskas, R., Maknickas, A. and Vislavičius, K. (2006). Investigation of performance of programming approaches and languages used for numerical simulation of granular material by the discrete element method. Computer Physics Communications. 175 (6): 404-415.‏

[38]. Delaney, G.W., Cleary, P.W., Morrison, R.D., Cummins, S. and Loveday, B. (2013). Predicting breakage and the evolution of rock size and shape distributions in Ag and SAG mills using DEM. Minerals Engineering. 50: 132-139.‏

[39]. Ting, J.M., Khwaja, M., Meachum, L.R. and Rowell, J.D. (1993). An ellipse‐based discrete element model for granular materials. International Journal for Numerical and Analytical Methods in Geomechanics. 17 (9): 603-623.‏

[40]. Shmulevich, I. (2010). State of the art modeling of soil–tillage interaction using discrete element method. Soil and Tillage Research. 111 (1): 41-53.‏

[41]. Cleary, P.W. and Sawley, M.L. (2002). DEM modelling of industrial granular flows: 3D case studies and the effect of particle shape on hopper discharge. Applied Mathematical Modelling. 26 (2): 89-111.‏

[42]. Munjiza, A. and Cleary, P.W. (2009). Industrial particle flow modelling using discrete element method. Engineering Computations.‏

[43]. Cleary, P.W. and Sinnott, M.D. (2008). Assessing mixing characteristics of particle-mixing and granulation devices. Particuology. 6 (6): 419-444.‏

[44]. Cleary, P.W. (2010). DEM prediction of industrial and geophysical particle flows. Particuology. 8 (2): 106-118.‏

[45]. Cleary, P.W. and Morrison, R.D. (2009). Particle methods for modelling in mineral processing. International Journal of Computational Fluid Dynamics. 23 (2): 137-146.‏

[46]. Just, S., Toschkoff, G., Funke, A., Djuric, D., Scharrer, G., Khinast, J. and Kleinebudde, P. (2013). Experimental analysis of tablet properties for discrete element modeling of an active coating process. AAPS PharmSciTech. 14 (1): 402-411.‏

[47]. McBride, W. and Cleary, P.W. (2009). An investigation and optimization of the ‘OLDS’elevator using Discrete Element Modeling. Powder Technology. 193 (3): 216-234.‏

[48]. Goniva, C., Kloss, C., Deen, N.G., Kuipers, J.A. and Pirker, S. (2012). Influence of rolling friction on single spout fluidized bed simulation. Particuology. 10 (5): 582-591.‏

[49]. Goniva, C., Kloss, C., Hager, A. and Pirker, S. (2010, June). An open source CFD-DEM perspective. In Proceedings of OpenFOAM Workshop, Göteborg (pp. 22-24).‏

[50]. Kloss, C., Goniva, C., Aichinger, G. and Pirker, S. (2009). Comprehensive DEM-DPM-CFD simulations-model synthesis, experimental validation and scalability. In Proceedings of the seventh international conference on CFD in the minerals and process industries, CSIRO, Melbourne, Australia.‏

[51]. Weerasekara, N.S., Powell, M.S., Cleary, P.W., Tavares, L.M., Evertsson, M., Morrison, R.D. and Carvalho, R.M. (2013). The contribution of DEM to the science of comminution. Powder Technology, 248, 3-24.‏