DEM Simulation to Determine the Influence on the Experimental Results of Tests of Iron Pellets When the Dimensions of the Test Device Are Varied

The current study is based on the DEM computer simulation of three experimental test devices with different dimensions to determine the difference in the results of the formation of shear and repose angles that the particles experience when grouped under the action of the gravitational force. In this respect, the experimental test devices with different height, width, and depth were geometrically modeled with iron pellet particles using morphology and a granulometric variation from 6 mm to 9 mm of equivalent diameter in its spherical shape. Depending on the results obtained, a reliable size of the experimental test device will be available to obtain the necessary data for a cor-rect adjustment of the calibration parameters for the DEM simulation of min-ing-metallurgical processes that use granulated material of iron pellet.


Introduction
The growing use of computer simulation through the discrete elements method (DEM) in mining, food, processing and farming industry to characterize and improve their processes has allowed an efficient development of them, allowing in turn to know the impact of their processes on the environment.
The method developed by Cundall, P. A., & Strack [1] at the beginning of the 70's and its subsequent current evolution [2] [3] have allowed together with the development calculation processes using graphics processors (GPU) [4], a sus-How to cite this paper: Aguilera-Carvajal, Y., Robledo, Y.T. and Cortes, S.P. Given the relevance of benefits associated with the development of DEM Simulations, it is where the calibration of its behavior parameters plays a fundamental role to generate simulations that adjust to reality. The importance of calibrations lies in the fact that the results of the simulations must be developed under the physical phenomena associated with the displacement dynamics of the particles and the interaction between them [12], due to this, it is necessary to be able to count on the characterization of the test devices to obtain reliable experimental parameters that will be used for the configuration of the adhesive model in computer simulations using the discrete elements method.
The fundamental concept of the discrete element method is based on the interaction of granular materials through a finite number of elements called particles which interact with each other when they come into contact. The interaction between them is defined by the contact models. Each contact model has a variety of different parameters, which must be assigned numerical values. These parameters, as well as other properties of the particles, influence the macroscopic behavior of the simulated material. Since there is no computer approach to determine numerical parameter values directly from a system response, the models must be calibrated. Calibration is an iterative process, which tries to identify a set of parameters whose resulting response from the system corresponds as closely as possible to the desired one. Due to the long calculation times of the DEM simulation, the number of calculations required should be kept as low as possible.
The iterative process of searching for adequate parameter values is highly inefficient and in some cases very difficult to find, so the experimental obtaining of the calibration values is a more efficient method and with high probability than the values entered into the models of adhesion give results adjusted to the phenomena of actual particle flow. In this sense, progress in the improvement of calibration processes has been developed by different researches [3] [13] but as indicated, the main disadvantage of iteration methods is the extensive computational time required, and in some cases can reach thousands of iterations before finding calibration values satisfactory as stated [14] in the research of a methodical calibration process for discrete element models The experiences carried out in accordance with establishing the variations of the results of the calibration tests in experimental test devices have shown that the dimensions of the test devices vary according to the different length scales used as indicated by the experience of [15].
Although advances in calibration methods have evolved notably, such as the use of artificial intelligence exposed by [16] based on genetic algorithms as well as other more advanced techniques in the use of these advanced calculation tools Y. Aguilera-Carvajal et al. Int. J. Modern Nonlinear Theory and Application exposed by [17] on models of surface response and neural networks such as the work carried out by [18], they require necessarily a satisfactory tuning and iterative search for calibration values in addition to the coupling of other external numerical simulation software.
From what has been said about the iterative techniques and based on the scarce development in the search for experimental validation methodologies, the current research focuses on the differences that occur when using parameter determination devices for DEM simulation in a simple way and with a high percentage of accuracy, which allows substantial savings in the calculation processes in search of satisfactory results for the calibration of simulation models. To obtain a comparative model and obtain their differences, three devices of different dimensions were geometrically modeled with spherical particles whose granulometry varied from 6 mm to 9 mm under identical material input conditions and simulation time.
This research work methodically establishes the variables that must be considered for the construction of a device model that allows results of experimentation for calibration of reliable DEM parameters of granular material in iron pellets, avoiding costly time to perform calibrations with the iterative method in addition to indicating which are the variables that generate significant differences in the results of the tests.

Formulation of the Discrete Element Method
The basic foundations of the discrete elements method are based on the fact that the particles basically describe two types of movements, the first is translational and the second is rotational, these happen when the interaction between them occurs and when a particle interacts with the walls of the system that contains it and/or with a circulating fluid in the means of analysis. When the movement of the particles occurs, an exchange of moment and energy is generated. Through the interactions that occur between particles it can determine the resulting force on it [19].
This approach allows using Newton's second law of motion to obtain the motion of individual particles, assuming that the system contains n discrete elements, the particles i have a mass (m i ) and a moment of inertia (l i ), the equations for the translational and rotational motion can be written as: where v i correspond to the speed of translation, while w i is the component of the angular velocity of the particle i, the contact force and torque exerted by a particle j or walls on the particle i are represented as F c , M ij , respectively, F c , M ij is the force of F c , M ij non-contact caused by particle k on particle i, while f is a component of the force produced by the particle-fluid interaction on particle i and F g in the gravitational force. In Figure 1 the forces and pairs involved in a DEM simulation [20] are shown schematically.

Contact Forces between Particles
The contact between two particles generally does not occur in a single point, due to the geometric deformation that the contact zone possesses is defined by a finite area, which is equivalent to the contact that occurs between two rigid bodies that slightly overlap in the DEM. In Figure 2 the contact force distribution over an area can be decomposed into two forces, one normal (F n ) and the other tangential (F t ) [21].
The normal unit vector (n) that allows the union of the center of the particles with a direction outside the particle i. Using a constitute model it is possible to obtain the contact between spherical particles shown in Figure 3.   At the contact interface k n and k t are the normal and tangential stiffness respectively, k n the friction coefficient and C n the contact damping coefficient.
The normal component of the norm force is the result of the sum of the elastic normal force (f ne ) and the damping force (f nd ). n ne nd

Elastic Normal Contac Force
The elastic normal contact force (f ne ) is proportional to the value of the normal stiffness (k n ) and to that of the penetration of the particle surfaces (u rn ).
where d corresponds to the distance between centers that is generated when contact between particles occurs, r i and r j is the radius of the particle i, j respectively. For the case studies that the particles do not have cohesion with each other, there will be no normal elastic contact force (or tension). What causes that when 0 rn u ≤ the value of the normal tension force 0 ne f = is set as null.

Normal Damping Contact Force
The normal damping contact force ( nd f ) is required to cause a decrease in the oscillations of the contact forces, in addition to dissipating the kinetic energy of the particles. For this type of force it must be assumed that it is viscous [22].
nd n rn f c v = (6) v rn is the relative velocity of the normal component between the centers of the particles.
( ) Considering that the damping coefficient (c n ) is contained within the critical damping (c cr ) [23], we have:

DEM Calibration
The calibration in the discrete elements method is based on obtaining optimal parameters at the time that an interaction between solid bulk particles occurs, which are entered as input values in the simulation of any process for obtaining Int. J. Modern Nonlinear Theory and Application solid material at in bulk, the data sought to be obtained by means of the calibration are the static and dynamic friction coefficients, coefficient of restitution and rolling coefficient [24].

Calibration by Angle of Repose Test
Obtaining the angle of repose consists in determining the angle of inclination with respect to the horizontal when solid bulk material of different dimensions is stacked, obtaining this value depends on the method used in experimental tests or computer simulation, and there are four ways to obtain this result [25].
The techniques of calibration by testing the angle of repose are pouring, draining and dynamic repose.

Calibration by Inclined Plane Test
The inclined plane test is used to determine the coefficients of sliding and rolling friction between the particle and a flat surface of different material. For the determination of the coefficient of sliding friction, a grain of material having a flat side is selected. The granular particle is placed on a horizontal flat wall, to keep the granulated material in a stationary state and thus avoid a possible rolling.
This approach ensures that the particle always slides before it starts rolling so that only the slip component can be measured. Then gradually the plane is inclined until the particle begins to slide [26].

Calibration by Coefficient of Restitution Test
This test is used to measure the coefficient of restitution of the particle (e), by means of the ratio between the relative velocities of the particle after and before the impact, when spherical particle of the granulated material is dropped on a horizontal surface. To facilitate obtaining the coefficient, the ratio of the square root of the rebound and the height of fall is usually used, and not the velocity ratio. The physically measured value of the coefficient of restitution is directly used in the DEM model so that the particle in the method has the same dissipation energy property as the real particle [27].

Results and Discussion
The design of experimentation (DOE) of this research work was based on the series of tests to determine the variations of the results against different dimensions of the same test device design to obtain DEM calibration parameters, for the simulation of the calibration process the Rocky DEM software was used to generate the virtualized tests of the tests of different devices. The material used is based on the properties and morphology of mineral iron pellets, where its main values are shown in Table 1. shown in Figure 4. The test process consists of filling the upper part of the device with particles, when the filling is finished, the opening of the gate that has the central division is configured, which causes a large part of the granular material to pass from the upper part of the device to the lower part, due to the effect of gravity, which begins to stack since the lower horizontal cover prohibits its displacement, with this an angle is generated between the lower cover and the material deposited on it, called the angle of repose, while the remaining material remains stored in the upper part since it fails to break its stationary state, it generates an angle with central division, this take the name shear angle.   In the discrete elements method, spherical shaped particles are used [30], as can be seen in Figure 6, for its approach in this study there are particles whose dimensions are within the range of 3 to 6.6 mm, at in turn, there will be 3 test devices similar to Figure 4, each with variations in their measurements.
The shape and main dimensions of the iron pellet particles are shown in Figure 7, where from the dimensions in accordance with the granulometry obtained, the equivalent volumes of the particles in their spherical shape were specified.    Table 3. of bulk material in the process it is 7.5 (ton/h) and for the test an adhesive distance of 0.0001 (m) will be used, with these parameters we obtain the results shown in Table 4.

Results
When there is a variation in the depth dimension of a test device as observed in test mechanism 2, in which the physical dimension in the depth of device 2 is decreased with respect to device 1, it is observed in the shear angle both on the right and on the left as observed in Table 5, which generates a standard deviation around 2% in each case, so the drawdown angle is not very sensitive to the change in the depth dimension.
The height dimension is modified as it happens in device 3 in which with respect to mechanism 2 the width and depth measurement is kept constant, generating an increase in its height, this produces a variation in the drawdown angles as observed in Table 5, the value of the angle generated on the right has a decreased of 16% in the case of the left side, the value decreases by 13% with respect to the values obtained from device 2.    In the case of the angles of repose, the value of the left side presents a decrease of 1˚ as seen in Table 6, for the case of the angle generated on the right side it is observed that it decreases by 2˚ with respect to the value obtained from device 1, with which a deviation of 6% is generated on the right side and 3% on the left sector.
In the values presented in Table 6, it is observed that the resting angles obtained in the test of device 3 suffer an increase with respect to the second, for the case of the angle that is formed to the right and to the left it is observed a growth of 18%.
There is a decrease in the calculation time, with device 1 being the one that takes the longest to process to obtain the final results as presented in Table 7, this is due to the fact that this mechanism has a greater quantity of particles and longer fill time when testing starts.
When decreasing one of the dimensions of the test device, the calculation time decreases, on the contrary when the dimensions of the test mechanism are increased, the calculation time suffers an increase, this is observed since device 3 suffers an increase in the height with respect to device 2, which generates an increase in the processing time as can be seen in Table 7.
When reducing one of the dimensions of the test device, the calculation time decreases, on the contrary, when the dimensions of the test mechanism are increased, the calculation time suffers an increase, this occurs when device 3 suffers and increase in height with respect to device 2, which generates an increase in processing time as can be seen in Table 7.

Conclusions
From the results obtained, it can be established that calibration is an important factor to be able to obtain reliable and realistic results of simulations using the discrete elements method of solid granular material composed of iron pellets for the simulation of industrial processes that process this type of material since it allows to define precisely the parameters that must be entered as input variables in them, that is why defining the appropriate dimensions will allow the calibra-  It should be noted that the granulometry used in the study must be respected so that the results of the DEM simulations adhere to the phenomenology that the particles experience in their real performance as well as the real geometric shape of the particle. In this study, the equivalent volumes of the sizes of the real particles were used to approximate their behavior using a spherical shape. The changes that occur when using morphology similar to the real one will vary according to the number of sides and different volumetric dimensions of them. It is to be expected that using granulometries in the ranges mentioned and with geometric shapes identical to the real particles, results are obtained that are adjusted to the phenomena that actually occur.
Future research should be focused on conducting tests with physical devices for experimental validation and with geometric morphology of the real particle.
Another research objective is to be able to determine the differences that occur in the experimental results due to the roughness of the surfaces and their effect on the static and dynamic friction coefficients of the particles.
From the above, it is a valuable contribution to the research to be able to determine with small modifications in the design of the values of the device of coefficients of friction and restitution with a single test procedure, in this case, it is necessary to investigate the governing and binding equations that describe what the particles experience during the granular material calibration testing process.
Obtaining reliable results of calibration values for DEM simulations will allow elucidating many of the phenomena that occur during the processes associated with the transport and processing of iron pellet materials with realistic and reliable simulations. The use of the methodology of discrete elements and its ability to couple other methodologies such as the elements and finite volumes will allow an important development in the mining industry, not only to verify and establish levels of processes but also in the effects that impact the useful life of assets by improving their maintenance strategies with a binding approach of condition monitoring techniques