Effect of Milling Parameters on DEM Modeling of a Planetary Ball Mill

Mohsen Mhadhbi

Laboratory of Useful Materials, National Institute of Research and Physicochemical Analysis, Ariana, Tunisia.

DOI: 10.4236/ampc.2023.134004                

Abstract

The effects of the milling parameters involving shape of powder particles, rotation speed, and ball-to-powder diameter (BPDR) on DEM modeling in the planetary ball mill were investigated. BPDR was varied from 1 to 10. The results revealed that the size and shape of the powder particles do not give a significant change in simulation results when BPDR attains maximum value of 10. The increasing of BPDR leads to the increase of simulation time and size. Hence, the effect of change of the powder particle shape on the calculated data size is not significant. The results also revealed that the increasing rotation speed increases impact energy between powder particles.

Keywords

DEM Modeling, Milling Parameters, Planetary Ball Mill, Particles Shape

1. Introduction

The discrete element method (DEM) has been commonly used by engineers and researchers to simulate granular materials in various fields. A key application of DEM modeling was the construction of the model and the calibration of the input parameters [1] [2] [3] [4] .

However, a number of research have been developed to investigate the balls motion inside ball mill based on the DEM simulations [5] [6] [7] [8] . For example, Li et al. [9] developed a scale-up model based on DEM modeling to study the performance of tumbling ball mills. They showed that the particle-ball contacts were the principal breakage mechanism of particles during milling. They also revealed that the estimated milling rates of the particles under various conditions were consistent with the experimental results. Jayasundara and Zhu [10] used DEM simulations and data driven approach to simulate the impact energy of particles in ball mills under different operating conditions. The obtained mode was validated by comparing the simulated results with the experimental data. The results showed that the milling performance decreases with increasing mill size and the impact energy was influenced by the operating conditions of the mill. Kim et al. [11] investigated the effects of the shape of the powder particles and the ball-to-powder diameter ratio on EDEM simulations and time in a planetary ball mill. They concluded that the simulation time and date size increase with increases in ball-to-powder diameter ratio. Kalpanadevi et al. [12] investigated the influence of four rice varieties on rice milling. They revealed that the larger grain length shows the elevated degree of milling. Xie et al. [13] used DEM simulation to study the effect of spherical and irregular polyhedral ore on the performance of SAG mill. They concluded that the liner wear of the polyhedron-sphere milling system is harder than pure spherical. Hence, Bibak and Banisi [14] investigated the influence of particle shape on load movement in tumbling mills by using physical and DEM modeling approach. They concluded that the cubical particles contributed 5% more in the high-energy impact action in comparison to the spherical particles. In addition, the simulation time increased by 35 times with changes in particles shape from spherical to cubical. Mio et al. [15] investigated the effect of rotation-to-revolution speed ratio and rotational direction on the milling performance using three kinds of planetary ball mill. They showed that the counter rotational direction is very powerful for enhancing the milling performance rather than the normal rotational direction. Hence, the optimum rotation-to-revolution speed ratio can be used to determine an effective milling condition. In this regard, Mori et al. [16] carried out the ball mill simulation in wet grinding using a tumbling mill and discussed the correlation between the experimental and simulation results. Moreover, it was found that the specific impact energy determined from the simulation correlates with the milling rate of the sample, showing that the energy is the most influential factor used in the control of wet milling.

Additionally, Cleary [17] analyzed the effect of rock shape representation in DEM on flow and energy utilization by using a pilot SAG mill. The results showed that the rock shape has a greater influence on the big size classes in the charge. In the research of Liu et al. [18] , the electrical system and multi-body dynamics model (EMBD) was combined with DEM method to explore the performance of a semi-autogenous (SAG) mill under various conditions. It was shown that the fluctuation of motor load torque, current, and mill speed in the initial stage is very greater than that in the stable running stage. Li et al. [19] studied the performance of vertical roller mill (VRM) by changing the material properties and operating conditions based on DEM simulation and experiment. They revealed that the crushing characteristics of VRM are very powerful for the operation parameters and optimization of the equipment. Oliveira et al. [20] conducted DEM simulation on a ball mill to study the effect of mill geometry, milling medium size, milling speed, and percent solids on apparent breakage rates. They concluded that the breakage rates increase considerably with stirrer speed and the increase in the milling media lead to an increase of breakage rates. Based on DEM simulation combining particle replacement model, Zeng et al. [21] investigated the modeling of grain breakage of in a vertical rice mill. It was found that the results obtained from the simulation are in agreement with those from experiments. Orozco et al. [22] used DEM simulations to study the breakage of particles in a 2D ball mill. They showed that the change in balls size has little effect on the evolution of milling. They also concluded that the milling rate increases with the number of balls of the same size. Lv et al. [23] studied the effect of various initial modes on the performance of a ball mill based on the electromechanical coupling model. It was found that the simulation results are consistent with the experimental data. Chimwani and Bwalya [24] used DEM simulation to investigate the effect of mill fill rate, end liner configuration, shape, and ball size on the distribution of forces acting on the liner. They revealed that the forces (radially and tangentially) increase with the increase of mill filling and with bigger balls. Based on DEM simulation, Cleary [25] investigated the effect of changes in mill operating conditions and particle properties on both the charge properties and mill power. It was found a sensitivity of mill behavior to the charge properties and the crucial significance of different hypothesis used DEM. Ji et al. [26] studied the influence of blockiness, particle aspect ratio (AR), and rotation speed on the mixing rate based on super-quadratic DEM simulations. It was found that the rotating speed has a great influence on the mixing rate. Jiang et al. [27] investigated the effects of number of lifers, filling degree, drum diameter, and rotation speed on the mixing uniformity of irregular gravel and sand materials in a rotating drum. They concluded that the rotation speed was the very important factor that affecting mixing, whereas the diameter of the drum was the lowest important. Bor et al. [28] used DEM simulation to study the effects of rotation speed, ball material, and milling time on the coating characteristics using two kinds of ball mills. The results showed that the impact forces have less affect in the coating than the contact number. In addition, the ball material properties influence no more difference in the coating.

The objective of this work is to investigate the effect of milling parameters, including shape of powder particles, rotation speed, and ball-to-powder diameter (BPDR) on DEM simulation results in a planetary ball mill and to develop a method to minimize the calculation cost during simulation.

2. DEM Model and Simulation Conditions

2.1. Apparatus Description

A high-energy planetary ball mill (Fritsch Pulverisette 7) was used in the simulation. This ball mill consists of supporting disc, two vials with volume of 15 ml, and balls with diameter of 15 mm. The operating principle of the planetary ball mill is illustrated in Figure 1.

Figure 1. The operating principle of the planetary ball mill.

2.2. Simulation Conditions

In this work, the simulation was started by mixing of powders with balls. The mill then rotated at a fixed speed to lift the ball-particles mixture. After that, the flow reached a steady state. Details on the discrete element method and simulation parameters can be found elsewhere [29] [30] .

The simulations were carried out using DEM simulation software EDEM Altair 2021.2 [31] . This software gives various methods for theory and simulation. Thus, Hertez-Mindlin contact model, based on Hertzian contact theory and Mindlin Deresiewicz model, was used in DEM simulations because of its useful and precise force calculation. Therefore, the shapes of the powder particles are illustrated in Figure 2. The volume of particle and balls was assumed to equal the volume of a sphere. However, the simulation method has been validated in a previous work [1] .

Figure 2. The different shapes of the powder particles used during simulation: (a) Single sphere; (b) Dual spheres; (c) Triple spheres.

3. Results and Discussion

In this study, the effect of the milling parameters on DEM modeling of a planetary ball mill was investigated at a laboratory scale. The results are detailed below.

3.1. Effect of the Shape of the Powder Particles

Figure 3 shows the snapshots of particles and milling balls at various particle shapes. As shown in the figure, there is no great difference in the motion patterns that all have cascading regime. This suggests that the shape of powder particles does not having a great effect on the motion patterns inside planetary ball mill. Feng et al. [32] concluded that the size and shape of the particles do not give a big effect in the motion pattern and simulation results when BPDR is over 20/3.

However, some works [33] [34] [35] assumed that the process of the planetary ball mill is very complex and the optimal parameters strongly depend on milling materials and operation conditions.

Figure 3. Snapshots of particles and milling balls at various particle shapes.

3.2. Effect of the Ball-to-Powder Diameter Ratio

Figure 4 presents the snapshots of particles and milling balls at various ball-to-powder diameter ratio (BPDR = 1, 3, 5, and 10). In this part, we assumed that the particle shape was a single-sphere. We can observe a slight difference in the motion patterns with increasing of BPDR from 1 to 10. For the case of BPDR equal to 1 and 3, we noted a cataracting regime. For the case of BPDR equal to 5 and 10, we noted a cascading regime, which the ball-to-ball contacts could be limited because of the contacts between the powder and the balls are retained for a prolonged milling time. On the other hand, the increase in BPDR leads to an increase in the number of ball-powder contacts. This indicates that the size of particles does not affect the motion pattern in the vial when BPDR reach a max value of 10. Thus, this phenomenon was also reported by several researchers [7] [36] [37] [38] .

Figure 4. Snapshots of particles and milling balls at various ball-to-powder diameter ratio: (a) BPDR = 1; (b) BPDR = 3; (c) BPDR = 5; (d) BPDR = 10.

3.3. Effect of the Rotation Speed

Figure 5 illustrates the snapshots of particles and milling balls at various rotation speed (100, 300, and 500 rpm). As shown in the figure, there is a significant difference in the motion patterns with increasing rotation speed. However, the rotation speed has a significant effect on DEM simulation in the planetary ball mill. Indeed, the figure demonstrates an acceptable agreement between the experimental and simulation results. On the other hand, the particles are represented by different colors based on their speeds. Thus, these results were similar to those in the ref. [7] . As explained by Li et al. [39] , the rotation speed can affect the average packing density: the density is greater at small rotation speed. In their work, Jayasundara et al. [40] investigated the effects of rotation speed and media loading on particle flow and milling performance in a horizontal stirred

Figure 5. Snapshots of particles and milling balls at various milling speeds: (a) 100; (b) 300; (c) 500 rpm.

mill. They revealed that the increase in mill speed resulted as an increase in the velocities of particles, which induces greater collisions between particles. However, Xie and Zhao [41] revealed that the effect of mill speeds is related to the rotation-revolution radius. They also concluded that the rotation speed is strongly related to the rotation-revolution speed ratio. Additionally, Zeng et al. [42] used DEM simulation to investigate the particle flow in a vertical rice mill. They showed that the influence of aspect ratio on the milling strongly depend on rotation speed.

4. Conclusion

Numerical modeling based on DEM approach has been performed on planetary ball mill. The effect of the shape of the powder particles, the ball-to-powder diameter ration, and the rotation speed of a planetary ball mill was investigated by simulating the behavior of the particles and milling balls through DEM modeling. The shape of powder particles does not affect significantly the motion patterns and simulation results. The simulation time increase with increasing the ball-to-powder diameter ratio. The increase of rotation speed greatly affected the motion patterns and simulation time. The results suggest that the ball-to-powder diameter ratio and rotation speed should be adjusted appropriately before milling.

Acknowledgements

The author would like to thank ALTAIR DEM for the EDEM software through the Academic License Program. The author also thanks Ms. Y. F. Zhang (Yavonne) for her technical support.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1]	Mhadhbi, M. (2022) Calibration of DEM Parameters to Simulate a Planetary Ball Mill. Advances in Materials Physics and Chemistry, 12, 155-162. https://doi.org/10.4236/ampc.2022.127011
[2]	Grima, A.P. and Wypych, P.W. (2011) Development and Validation of Calibration Methods for Discrete Element Modelling. Granular Matter, 13, 127-132. https://doi.org/10.1007/s10035-010-0197-4
[3]	Rackl, M. and Hanley, K.J. (2017) A Methodical Calibration Procedure for Discrete Element Models. Powder Technology, 307, 73-83. https://doi.org/10.1016/j.powtec.2016.11.048
[4]	Yan, Z., Wilkinson, S.K., Stitt, E.H. and Marigo, M. (2015) Discrete Element Modelling (DEM) Input Parameters: Understanding Their Impact on Model Predictions Using Statistical Analysis. Computational Particle Mechanics, 2, 283-299. https://doi.org/10.1007/s40571-015-0056-5
[5]	Mishra, B.K., and Rajamani, R.K. (1992) The Discrete Element Method for the Simulation of Ball Mills. Applied Mathematical Modelling, 16, 598-604. https://doi.org/10.1016/0307-904X(92)90035-2
[6]	Tavares, L.M. (2017) A Review of Advanced Ball Mill Modelling. KONA Powder and Particle Journal, 34, 106-124. https://doi.org/10.14356/kona.2017015
[7]	Rodriguez, V.A., de Carvalho, R.M. and Tavares, L.M. (2018) Insights into Advanced Ball Mill Modelling through Discrete Element Simulations. Minerals Engineering, 127, 48-60. https://doi.org/10.1016/j.mineng.2018.07.018
[8]	He, K., Jia, M., Sun, F., Yang, Q., Wu, B., Li, X. and Meng, C. (2022) Simulation of Grinding Media Motion and Collisions in Wet-Operated Accelerator-Mill Using DEM Coupled with CFD. Minerals, 12, 341. https://doi.org/10.3390/min12030341
[9]	Li, Y., You, Y., Gou, D., Yu, A. and Yang, R. (2022) A DEM Based Scale-Up Model for Tumbling Ball Mills. Powder Technology, 409, Article ID: 117854. https://doi.org/10.1016/j.powtec.2022.117854
[10]	Jayasundara, C.T. and Zhu, H.P. (2022) Impact Energy of Particles in Ball Mills Based on DEM Simulations and Data-Driven Approach. Powder Technology, 395, 226-234. https://doi.org/10.1016/j.powtec.2021.09.063
[11]	Kim, K., Jiang, T., Kim, N. and Kwon, C. (2022) Effects of Ball-to-Powder Diameter Ratio and Powder Particle Shape on EDEM Simulation in a Planetary Ball Mill. Journal of the Indian Chemical Society, 99, Article ID: 100300. https://doi.org/10.1016/j.jics.2021.100300
[12]	Kalpanadevi, C., Singh, V. and Subramanian, R. (2019) Impact of Physicochemical Properties on Duration and Head Rice Yield during Abrasive and Friction Milling of Rice. Journal of Food Science and Technology, 56, 3900-3909. https://doi.org/10.1007/s13197-019-03861-z
[13]	Xie, C.H., Ma, H.Q., Song, T. and Zhao, Y.Z. (2021) DEM Investigation of SAG Mill with Spherical Grinding Media and Non-Spherical Ore Based on Polyhedron-Sphere Contact Model. Powder Technology, 386, 154-165. https://doi.org/10.1016/j.powtec.2021.03.042
[14]	Bibak, Z. and Banisi, S. (2021) A Combined Physical and DEM Modelling Approach to Investigate Particle Shape Effects on Load Movement in Tumbling Mills. Advanced Powder Technology, 32, 916-930. https://doi.org/10.1016/j.apt.2021.01.034
[15]	Mio, H., Kano, J., Saito, F. and Kaneko, K. (2004) Optimum Revolution and Rotational Directions and Their Speeds in Planetary Ball Milling. International Journal of Mineral Processing, 74, S85-S92. https://doi.org/10.1016/j.minpro.2004.07.002
[16]	Mori, H., Mio, H., Kano, J. and Saito, F. (2004) Ball Mill Simulation in Wet Grinding Using a Tumbling Mill and Its Correlation to Grinding Rate. Powder Technology, 143, 230-239. https://doi.org/10.1016/j.powtec.2004.04.029
[17]	Cleary, P.W. (2019) Effect of Rock Shape Representation in DEM on Flow and Energy Utilisation in a Pilot SAG Mill. Computational Particle Mechanics, 6, 461-477. https://doi.org/10.1007/s40571-019-00226-3
[18]	Liu, Z., Wang, G., Guan, W., Guo, J., Sun, G. and Chen, Z. (2022) Research on Performance of a Laboratory-Scale SAG Mill Based on DEM-EMBD. Powder Technology, 406, Article ID: 117581. https://doi.org/10.1016/j.powtec.2022.117581
[19]	Liu, C., Chen, Z., Zhang, W., Mao, Y., Xu, P. and Xie, Q. (2022) Analysis of Vertical Roller Mill Performance with Changes in Material Properties and Operating Conditions Using DEM. Minerals Engineering, 182, Article ID: 107573. https://doi.org/10.1016/j.mineng.2022.107573
[20]	de Oliveira, A.L.R., de Carvalho, R.M. and Tavares, L.M. (2021) Predicting the Effect of Operating and Design Variables in Grinding in a Vertical Stirred Mill Using a Mechanistic Mill model. Powder Technology, 387, 560-574. https://doi.org/10.1016/j.powtec.2021.04.057
[21]	Zeng, Y., Mao, B., Jia, F., Han, Y. and Li, G. (2022) Modelling of Grain Breakage of in a Vertical Rice Mill Based on DEM Simulation Combining Particle Replacement Model. Biosystems Engineering, 215, 32-48. https://doi.org/10.1016/j.biosystemseng.2021.12.022
[22]	Orozco, L.F., Nguyen, D.H., Delenne, J.Y., Sornay, P. and Radjai, F. (2020) Discrete-Element Simulations of Comminution in Rotating Drums: Effects of Grinding Media. Powder Technology, 362, 157-167. https://doi.org/10.1016/j.powtec.2019.12.014
[23]	Lv, W., Wang, G. and Tian, H. (2018) Electromechanical Dynamic Behaviour and Start-Up Evaluation of Tumbling Ball Mills. Mathematical Problems in Engineering, 2018, Article ID: 3515308. https://doi.org/10.1155/2018/3515308
[24]	Chimwani, N. and Bwalya, M.M. (2020) Exploring the End-Liner Forces Using DEM Software. Minerals, 10, 1047. https://doi.org/10.3390/min10121047
[25]	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, 79-114. https://doi.org/10.1016/S0301-7516(01)00037-0
[26]	Ji, S., Wang, S. and Zhou, Z. (2020) Influence of Particle Shape on Mixing Rate in Rotating Drums Based on Super-Quadric DEM Simulations. Advanced Powder Technology, 31, 3540-3550. https://doi.org/10.1016/j.apt.2020.06.040
[27]	Jiang, S., Ye, Y., He, M., Duan, C., Liu, S., Liu, J., Xiao, X., Zhang, H. and Tan, Y. (2019) Mixing Uniformity of Irregular Sand and Gravel Materials in a Rotating Drum with Determination of Contact Model Parameters. Powder Technology, 354, 377-391. https://doi.org/10.1016/j.powtec.2019.06.005
[28]	Bor, A., Jargalsaikhan, B., Lee, J. and Choi, H. (2020) Effect of Different Milling Media for Surface Coating on the Copper Powder Using Two Kinds of Ball Mills with Discrete Element Method Simulation. Coatings, 10, 898. https://doi.org/10.3390/coatings10090898
[29]	Mhadhbi, M. (2021) Modelling of the High-Energy Ball Milling Process. Advances in Materials Physics and Chemistry, 11, 31-44. https://doi.org/10.4236/ampc.2021.111004
[30]	Mhadhbi, M. (2021) Simulation of a Laboratory Scale Ball Mill via Discrete Element Method Modelling. Advances in Materials Physics and Chemistry, 11, 167-175. https://doi.org/10.4236/ampc.2021.1110016
[31]	Altair’s 2021.2 Simulation Software. https://www.altair.com/simulation-2021-2
[32]	Feng, Y., Han Kook and Owen, D.R.J. (2004) Discrete Element Simulation of the Dynamics of High Energy Planetary Ball Milling Processes. Materials Science and Engineering: A, 375-377, 815-819. https://doi.org/10.1016/j.msea.2003.10.162
[33]	Christine, F.B. and Arno, K. (2013) Process Engineering with Planetary Ball Mills. Chemical Society Reviews, 42, 7660-7667. https://doi.org/10.1039/c3cs35455e
[34]	Marica, B., Luca, G., Mirco, D., Cristy, A.R., Nicola, P. and Scardi, P. (2015) Modeling of the Planetary Ball-Milling Process: The Case Study of Ceramic Powders. Journal of the European Ceramic Society, 36, 2205-2212.
[35]	Giraud, M., Gatumel, C., Vaudez, S., Nos, J., Gervais, T., Bernard-Granger, G. and Berthiaux, H. (2021) Investigating Grinding Mechanisms and Scaling Criteria in a Ball Mill by Dimensional Analysis. Advanced Powder Technology, 32, 2988-3001. https://doi.org/10.1016/j.apt.2021.06.016
[36]	Khasenova, R., Komarov, S., Ishihara, S., Kano, J. and Zadorozhnyy, V.Y. (2018) Discrete Element Method Simulations of Mechanical Plating of Composite Coatings on Aluminum Substrates. Surface and Coatings Technology, 349, 949-958. https://doi.org/10.1016/j.surfcoat.2018.04.056
[37]	Ashrafizadeh, H. and Ashrafizaadeh, M. (2012) Influence of Processing Parameters on Grinding Mechanism in Planetary Mill by Employing Discrete Element Method. Advanced Powder Technology, 23, 708-716. https://doi.org/10.1016/j.apt.2011.09.002
[38]	Kim, K.C., Jiang, T., Xu, Y., Kim, N.I., Jin, H.S. and Kim, J.C. (2021) Application of Discrete Element Simulation in Mechanical Activation of Boron Concentrate. Powder Technology, 382, 441-453. https://doi.org/10.1016/j.powtec.2020.12.031
[39]	Li, M., et al. (2021) DEM Simulation on Mixing Characteristics and Macroscopic/Microscopic Flow Behaviors of Different-Shaped Sphero-Cylinders in a Rotating Drum. Industrial & Engineering Chemistry Research, 60, 8874-8887. https://doi.org/10.1021/acs.iecr.1c00962
[40]	Jayasundara, C.T., Yang, R.Y., Yu, A.B. and Rubenstein, J. (2010) Effects of Disc Rotation Speed and Media Loading on Particle Flow and Grinding Performance in a Horizontal Stirred Mill. International Journal of Mineral Processing, 96, 27-35. https://doi.org/10.1016/j.minpro.2010.07.006
[41]	Xie, C. and Zhao, Y. (2022) Investigation of the Ball Wear in a Planetary Mill by DEM Simulation. Powder Technology, 398, Article ID: 117057. https://doi.org/10.1016/j.powtec.2021.117057
[42]	Zeng, Y., Mao, B., Li, A., Han, Y. and Jia, F. (2022) DEM Investigation of Particle Flow in a Vertical Rice Mill: Influence of Particle Shape and Rotation Speed. Powder Technology, 399, Article ID: 117105. https://doi.org/10.1016/j.powtec.2021.117105