Abstract

To investigate the impact of lifting velocity and container shape on angle of repose, the fixed-base cylinder method was performed using three types of container shape. The container shape was lifted a series of lifting velocities. Six size fractions of iron ore particles ranging from coarser to fine particles were used as the test materials. And the sand-pile calibration method was then used to calibrate the contact parameters of iron ore particles. Results show angle of repose decreased exponentially with the lifting velocity, while it appeared approximately to be invariant to particle shape, for all size fractions. The sand pile highly depends on the container shape at a low lifting velocity but appears to be invariant to particle size for a high lifting velocity. And then a predictive equation is established and a very close agreement between the predicted and measured angle of repose is attained. Finally, a series of DEM simulations considering the irregular particle shape are conducted by means of sphere clump method to calibrate the contact parameters and are in good visual agreement with the experimental results, indicating the “tuned” contact parameters as well as the applicability of the predicted equation.

1. Introduction

Size reduction of minerals by grinding hitherto is a significant step in whole stage of mineral processing in which feeds are reduced in size to liberate impurities and to change the particle size distribution. The operating costs associated with size reduction are a significant discipline in mineral industry [1, 2]. Grinding is known to be an energy-intensive operation process, since it consumes approximately two-thirds of the total energy of the mineral processing plants [3]. The literature reported that the power estimation of tumbling mills plays an important role in determining the grinding efficiency and mill performance [4, 5]. To this end, establishing power equations of tumbling mill has attracted much attention in recent years. At present, researchers presented that mill power strongly affected by the flow properties of materials in mills [68]. Besides, with the rapid development of computing power and advanced contact algorithm, the Discrete Element Method (DEM) has been extensively used as a leading tool in tumbling mills. DEM simulations hitherto have been demonstrated extremely desirable to predict what happens in reality as well as the quantitatively accurate information representation inside mills, while the accuracy of outcomes for DEM simulations highly depends on the input parameters [914]. So far, the sand-pile calibration has been employed exclusively in determining the input parameters. The aim of this approach was to make an approximate comparison between angle of repose of experimental results and the one obtained by DEM numerical simulations. Thus, the investigation of angle of repose is of importance in studying the performance of tumbling mills.

The angle of repose is one of the most significant macroscopic parameters in characterizing the flow behavior of powders and granular materials. It is closely related to many phenomena, including packing [15, 16], avalanching [17, 18], and segregation and stratification [1921], and it is therefore a research focus in recent years. The angle of repose, as a basic characteristic of granular materials in nature, has been indicated varying from 25° for smooth particles to 45° for irregular particles [9, 22, 23]. Currently, The literature reported that angle of repose was closely related to physical properties of granular materials such as particle shape and size [2430] and materials properties such as shear modulus and particle moisture content and friction coefficients [12, 3133]. In the works of scholars who have studied the impact of particle shape and size on angle of repose in the literature [2430], the angle of repose was proportional to friction coefficient and was inversely proportional to particle size. Besides, the particle shapes used were smooth, and the container shapes used were usually a cylinder. Moreover, a small lifting velocity was used to determine the impact of different factors on angle of repose [3436]. Since it was difficult to form the sand pile in their works for high and low enough lifting velocities, however, the particle shapes used in this study were highly irregular and the particle size distribution was evidently heterogeneous. It is essential to recognize that angle of repose not only highly depends on the experimental method involved but also is strongly sensitive to container shape, due to different contact behavior [37]. To this end, the impacts of lifting velocity and container shape on angle of repose were needed.

The purpose of this study was to determine the impact of container shape and lifting velocity on the angle of repose. The test materials were iron ore samples. A series of experiments were performed on different container shapes and lifting velocities. Then the contact parameters were determined by means of the comparison of experiments and DEM simulations. Our achievements will be of great benefit for determining the flow property of iron ore particles in industrial applications.

2. Materials and Methods

2.1. Experimental Setup

In this study, the sand-pile test, namely, fixed-base cylinder method, was used to determine the angle of repose. This method is in accordance with the most common forming process of sand pile in nature [34]. Shown in Figure 1, the stainless steel plate used was 0.425 m in length and 0.45 m in width and was fixed horizontally on the bottom of the aluminum frame. Two fixed pulleys were installed at the top of the aluminum frame. The container was placed perpendicularly on the top of the stainless steel plate. The grooved wheel, which was coupled to a motor shaft, was used to lift the container at different lifting velocities. Finally, the iron ore particles were discharged out of the container under the natural force of gravitation forming a pile of materials and accumulated from the center of the steel plate and gradually expanded to the edge.

To study the impact of container shape on angle of repose, three kinds of containers with typical cross-sections were employed, as shown in Figure 2. The stainless steel container, opened at both ends, was 0.27 m in height and 0.02 m in thickness. The cross-sectional areas were approximately equal as well as the volume of container (Table 1).

2.2. Materials

Six different size fractions ranging from coarse particles to fine ones were used as the test materials. The iron ore samples with 67.46% grade were obtained from the Xuzhou iron mine (located in Jiangsu Province, China). The larger lumps were broken by an industrial jaw crusher followed by a laboratory-scale jaw crusher to approximately −50 mm. The crushed particles were then sieved carefully (20 min) by a vibrating screen to obtain the required size fractions (Table 2). The weight of iron ore particles used in this study was 2.28 kg.

2.3. Angle of Repose

As the container was lifted at a constant velocity, the iron ore particles spontaneously formed a sand pile. Angle of repose, as a basic parameter of particle property, defined as the angle between the inclined surface of the pile and the base surface, is known to qualitatively analyze the flow characteristic of granular particles, shown in Figure 3. For high angle, the iron ore particle starts to flow. For low angle, the iron ore particles stop flowing and become stable gradually. After sand-pile formation, the pile angles were recorded carefully with a high-precision digital protractor, and the average value was obtained and expressed as follows:where are the measured angle in the orthogonal axes with a unit of rad.

2.4. Procedure

Firstly, the prepared iron ore particles were poured into the container by a hopper. And then the container was lifted at different velocities, for example, 0.003 m/s, 0.007 m/s, 0.01 m/s, 0.03 m/s, 0.07 m/s, 0.1 m/s, 0.3 m/s, and 0.5 m/s. Finally, the sand piles with different slope angles were formed, and therefore the corresponding angles of reposes were determined. Each test was repeated three times and the average values of repose angle were obtained.

3. Experimental Results and Discussions

3.1. Effect of Lifting Velocity on Angle of Repose

Figure 4 shows the angle of repose plotted against the lifting velocity for three different containers. As seen in Figure 4, it changes approximately exponentially with the lifting velocity for all size factions. For lifting velocity below 0.1 m/s, the angle of repose decreases rapidly with the increase of lifting velocity, while it decreases steadily above 0.1 m/s. The evident expanded effect of edge appears when a lifting velocity is higher enough. A high lifting velocity will cause the fierce collision in particle systems because of the high impact velocity, and there is therefore large amount of energy consumed. Thus, it is hard to form a high sand pile. For low lifting velocity, the consumed collision energy decreases. As a result, it will have high potential to form a high sand pile. Most of particles used in the literature commonly have regular shape and smooth surface such as glass beads and chocolate candies. It is difficult to form the sand piles at a low enough lifting velocity, due to the reduction on friction coefficient between particle-particle and particle-wall. Moreover, the iron ore particle shapes are highly irregular and the tendency of rotational motion of iron ore particles is harder than that of regular particles. Hence, the sand-pile formation of iron ore particles is markedly different from others. Additionally, it is clear that the lifting velocity has a stronger effect on angle of repose than particle size.

3.2. Effect of Particle Size on Angle of Repose

Figure 5 shows the variation of angle of repose as a function of particle size for different size fractions. The experimental errors are inevitable due to the experimental replication, and a maximum error of 1.47 and a minimum error of 0.12 are obtained. For most studies, researchers presented that increasing particle size will decrease angle of repose. As shown in Figure 5, the angle of repose slightly increases with the increase of particle size. As the lifting velocity increases, the variation on angle of repose decreases up to a minimum value and then increases slightly. For lifting velocity above 0.3 m/s, the effect of lifting velocity on angle of repose can be neglected. The effect of particle size on angle of repose is mainly due to the balance between interparticle forces and interaction forces. When particles are small, the interparticle force will play the dominant role on angle of repose, and therefore the particles will be considered as a cohesive group [22]. When the particle size is large enough, the interparticle force no longer plays a dominant role in angle of repose so that the particles behave as individual ones [22]. In addition, the shapes of coarser particles are more irregular than that of fine particles so that the friction coefficient of coarse particle is higher, resulting in high potential to form the sand pile. Moreover, the collisions of particle-particle and particle-wall are relatively gentle at a low lifting velocity, but that will be subjected to fierce collision at a high lifting velocity, resulting in more energy consumed. To sum up, it is clear that angle of repose is only a weak variation with the particle size so that the particle size effect can be neglected. Therefore it is reasonable in the practice of DEM simulation to assume a constant value of contact parameters for various particle size distributions.

3.3. Effect of Container Shape on Angle of Repose

Figure 6 shows the impact of container shape on angle of repose for three different lifting velocities, that is, a low velocity of 0.01 m/s, an intermediate velocity of 0.1 m/s, and a high velocity of 0.5 m/s. For a low lifting velocity, angle of repose formed by rectangle container (RE) is higher than that formed by round container (R). For an intermediate velocity, angle of repose formed by rectangle container (RE) is lower than that formed by round container (R). However, angles of reposes formed by rectangle container (RE) and round container (R) are obviously greater than that formed by triangle container (T). For a high velocity, angle of repose appears approximately to be invariant to particle size. It is owing to the fact that the container shape plays significant role in the contact behavior. Broadly speaking, ball mill is fitted with various types of lifters such as lifter bar, wave lifter, and triangle lifter. Hence, the contact behaviors between the particles and lifters are markedly different, resulting in forming diverse load behavior of charge.

In this research, the coarse particles of 4 × 7 mesh and the fine particles of 14 × 20 mesh, lifted at a low velocity of 0.01 m/s and a high velocity of 0.1 m/s, are used to further study the boundary effect, as shown in Figure 7. From the top view, the bottom shapes of sand piles are profoundly affected by the cross-section of each container, for the lifting velocity of 0.01 m/s. However, the boundary shapes become the same for the lifting velocity of 0.1 m/s, and therefore the boundary effect disappears.

3.4. Formulation of a Predictive Equation

The above results show that the angle of repose is highly sensitive to the lifting velocity, while it appears to be invariant to particle size. The measured angle of repose is fitted using a nonlinear equation, and the equation is expressed as follows:where are model constants, is lifting velocity, and is average particle size

Figure 8 shows the measured angle of repose plotted against the predicted values using a nonlinear model. The data listed in Figure 8 demonstrated that good agreements between measured and predicted angle of repose indicate the applicability of this equation in this work. This equation can enormously shorten the measurement time of angle of repose and improve the efficiency for studying the input parameters of DEM simulations.

4. DEM Simulations

To determine the DEM input parameters, the sand-pile calibration method is used to make an approximate comparison between angle of repose of experimental results and the one obtained by DEM numerical simulations. A series of DEM simulations are conducted on the samples of 4 × 7 mesh due to the weak sensitivity to particle size. Considering the rationality of lifting velocity, the round container is lifted at a constant velocity of 0.01 m/s. In this study, twenty-eight particles were used to generate the 3D digitized and meshed shapes by a high-accuracy 3D scanner (TEXU-BLU, 2M pixels). Sphericity can define the shape of the particle to a very good extent with a simple method and has been used extensively in powder technology [38].where is the real surface area of the iron ore particle, is the surface area of sphere determined by the same volume of iron ore particle, and is the real volume of iron ore particle. The real physical properties of twenty-eight particles are calculated as shown in Table 3.

The geometrical dimensions of iron ore particles are determined on the basis of the 3D digitized and meshed shapes. Considering the equivalent sphere with the same volume of the real particle, the equivalent surface areas are calculated in Table 3, as shown in Figure 9. The equivalent surface area versus the real surface area is fitted by the method of least squares. The fitting value has a sphericity of 0.758 and a correlation coefficient of 0.95 which are slightly smaller than is the case for the average value (0.782), but both are less than the sphericity of sphere. It is demonstrated that the iron ore particle is highly irregular comparing with the sphere so that the simplifications and assumptions of irregular particle as simple spherical shape are unrealistic.

Hence, to allow the irregular shape to be accurately determined, the sphere clump method is employed in modelling the particle shape. An example of a modelled iron ore particle is shown in Figure 10. The no-slip Hertz–Mindlin model is employed to figure out the contact force between particle systems within the EDEM software package (DEM—Solutions Ltd., Edinburgh, UK). To achieve computational accuracy without increasing the amount of calculation, the geometrical model used in the work is a cluster of thirty overlapping spheres. In EDEM software any kind of irregular geometrical shape of particles can be generated as a clump composed of several touching or overlapping spheres. Therefore, the contact detection between the sphere clumps is sphere-based and, therefore, the discrete element algorithm of the sphere clumps is fully available for calculating the contact forces.

The mechanical properties of iron ores and walls are conducted using hydraulic MTS machine. The contact parameters between particle and wall are determined by means of drop test and pin-on-disk tribometer tester. The interparticle contact parameters are determined previously by authors. It is useful to determine which contact parameter has the greatest impact on angle of repose. A series of DEM simulations present that friction coefficient has the primary impact on angle of repose rather than restitution coefficient. Then adjustments could be made to the friction coefficient until DEM simulations match the experimental results. The simulated test shows a good visual agreement with the experimental results, shown in Figure 11.

Using the appropriate parameters used in the slump tests, a series of DEM simulations on angle of repose are conducted. Figure 11 compares the experimental and simulated slump test on iron ore particles. It indicates that the DEM simulations using the sphere clump methods are in good visual agreement with the experimental results, implying the “tuned” contact parameters. In order to test the obtained contact parameters, the comparisons between DEM simulations and experimental results are studied at the lifting velocity of 0.01 m/s and 0.1 m/s, using three containers, as shown in Figure 12.

Figure 12 shows the relationship among measured and calculated and DEM simulation results. As seen in Figure 12, the DEM simulations considering the irregular shape are in good visual agreement with the experimental results as well as the calculated results. The standard deviation is approximately 1.6 degrees, which almost within the measurement error.

5. Conclusions

The results presented above contribute to a better understanding of the impact of lifting velocity and particle size on angle of repose of iron ore particles. Three types of container shape were used to conduct the experiments on six size fractions ranging from coarser particles to fine ones. And the sand-pile calibration method was then used to calibrate the contact parameters. Results have shown that angle of repose decreased exponentially with the lifting velocity, while it appeared approximately to be invariant to particle shape. Hence, the impact of lifting velocity was more profound than that of particle size. For a low lifting velocity, the sand piles highly depended on the container shape. For a high lifting velocity, it appeared to be invariant to the particle size. Subsequently, a predictive equation is established to quantitatively study the impact of lifting velocity and particle size on angle of repose of iron ore particles. Finally, a sphere clump method was used to generate the geometrical model of real particles using a cluster of thirty overlapping spheres. Based on the sand-pile and DEM simulations, a very close agreement between the predicted and measured angle of repose is attained. The contact parameters determined by DEM simulations show a good approximation to the experimental and calculated results, indicating the “tuned” contact parameters as well as the applicability of the predicted equation.

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

Yuxing Peng proposed the research topic; Zhangfa Yu contributed materials to the work; Zhencai Zhu guided and supervised the entire project; Tongqing Li wrote the manuscript and analyzed the data; Yuxing Peng and Zixin Yin made the revisions of the manuscript.

Acknowledgments

The research reported here was supported by the National Nature Science Foundation of China (Grant no. 51475458), Key Project of National Natural Science Foundation of China (Grant no. U1510205), and the Program for Changjiang Scholars and Innovative Research Team in University (IRT_16R68). The authors also wish to thank the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD) and the Top-Notch Academic Programs Project of Jiangsu Higher Education Institutions (TATP).