Universal Approach for DEM Parameters Calibration of Bulk Materials

DEM parameters calibration is the most important step in preparing a DEM model. At the same time, the lack of a universal approach to DEM parameters calibration complicates this process. The paper presents the author’s approach to creating a universal calibration approach based on the physical meaning of the friction coefficients and conducting symmetrical experiments at full scale and in a simulation, as well as the implementation of the approach in the form of a physical test rig. Several experiments were carried out to determine the DEM parameters of six material–boundary pairs. The resulting parameters were adjusted using a refinement experiment. The results confirmed the adequacy of the developed approach, as well as its applicability in various conditions. The limitations of both the approach itself and its specific implementation in the form of a test rig were identified.


Introduction
The discrete element method is the most popular approach for computer modeling of bulk materials' behavior. The corresponding software that implements DEM in the user graphic interface is a highly effective tool for optimizing mining equipment. Lately, DEM has often been used in conjunction with CFD and other methods, which opens up the possibility of calculating complex multiphase processes [1][2][3][4][5].
A number of input parameters in DEM software directly represent material properties (shape and size of particles, density, etc.). Friction coefficients (Table 1) have a direct physical representation; however, in DEM software these parameters are integrated into DEM codes of contact models and may affect the behavior of bulk materials in different ways, that is, they are code dependent [6,7]. Since the values of the friction coefficients (DEM parameters) significantly affect the behavior of bulk materials, in order to build an adequate model, they have to be calibrated [8][9][10]. Many researchers offer their approaches and solutions for the DEM parameters calibration. The whole set of existing calibration methods can be divided into two groups: the bulk calibration approach (BCA) and the direct measuring approach (DMA) [11]. A collaborative approach is also often used. In BCA, a laboratory experiment is performed first (for 1. Code dependence (depending on DEM software and contact model).

2.
A possibility for several sets of DEM parameters to provide similar bulk responses in the simulation. 3.
The need to significantly simplify the model (including the shape and size of particles) to speed up the calculation. Thus, the use of calibrated values of DEM parameters is limited, as well as their dependence on a specific application (technological process).

4.
Imperfection and inaccuracy of modern measuring tools (including visual estimation of the bulk material responses using machine vision).
Nevertheless, a number of researchers have proposed approaches that can be called universal with a number of limitations [21][22][23]. Many studies are aimed at reducing the number of simulations required to achieve the desired result using optimization algorithms or at obtaining a unique set of calibrated DEM parameters either by introducing special criteria or by estimating all the possible factors that affect the bulk responses [24,25]. In this case, BCA is applied, which indicates a possible loss of the physical meaning of the obtained DEM parameters. The question also arises whether the obtained set of parameters is unique, that is, the bulk responses values are achieved only with this set of DEM parameters values.
The article describes the developed author's approach and a test rig for the bulk materials' DEM parameters calibration. The approach is based on the transfer of the friction coefficients' physical meaning into the measurement of bulk materials macro parameters using a high-speed camera and the calibration of the obtained DEM parameters set based on a refinement experiment.

Materials and Methods
The developed approach is a result of the research conducted by the authors since 2017. Approaches to using the BCA method directly using the design of the experiment were considered [26]. One way or another, this approach did not solve the main problems given in the introduction. In the course of the research, the relationship between the physical meaning of the friction coefficients and DEM parameters was confirmed.
The physical meaning of these coefficients is as follows: • For static friction (SF), this is the slope (tangent of the repose angle) at which the particle begins to slide over the surface (Figure 1a): where ∝-angle of shelf incline.
• For dynamic friction, the value is determined by the sliding time on a surface with a certain angle of repose ( Figure 1a): where ∝-angle of shelf incline, S-distance of particle slide movement, t-time of sliding, g-gravity acceleration.
• For the coefficient of restitution, the value is determined by the angles and velocities of rupture and reflection ( Figure 1b): Symmetry 2021, 13, 1088 3 of 14 lems given in the introduction. In the course of the research, the relationship between the physical meaning of the friction coefficients and DEM parameters was confirmed. The physical meaning of these coefficients is as follows: • For static friction (SF), this is the slope (tangent of the repose angle) at which the particle begins to slide over the surface (Figure 1a): where ∝-angle of shelf incline. • For dynamic friction, the value is determined by the sliding time on a surface with a certain angle of repose ( Figure 1a): where ∝-angle of shelf incline, -distance of particle slide movement, -time of sliding, -gravity acceleration.
• For the coefficient of restitution, the value is determined by the angles and velocities of rupture and reflection ( Figure 1b): The dynamic friction coefficient can be determined by the sliding time of the material along the inclined shelf. The coefficient of restitution can be determined by conducting experiments on the collision of particles with a vertical shelf [27][28][29]. Based on the physical meaning of the static friction coefficient, it can be determined by the angle of inclination of the shelf. A series of numerical experiments were carried out. In the first version, similar to [30], a box (without bottom) with particles was placed on a shelf, after which the shelf was slowly raised and the angle at which the box began to slide on the shelf was fixed. The second version involved the use of a counterweight for a box with particles standing on a horizontal plane. The counterweight force was slowly increased, and the force at which the box began to slide along the plane was recorded. However, the experiments carried out in the simulations did not show a correlation between SFPP, SFPB, and the angle of inclination of the shelf, as well as the force of the counterweight.
The approach was designed in such a way as to neutralize the main problems described in the introduction. For this, when developing, the authors started from the physical meaning of the coefficients, but transferring the meaning from a single particle to a portion of the investigated bulk material. It is possible to preserve the physical meaning of the calibrated DEM parameters, as well as reduce the number of possible parameter sets or even get a unique set. However, the approach has to include the determination of bulk material responses similar to the BCA method because of the problems with determining SF according to the physical meaning. This makes it possible to The dynamic friction coefficient can be determined by the sliding time of the material along the inclined shelf. The coefficient of restitution can be determined by conducting experiments on the collision of particles with a vertical shelf [27][28][29]. Based on the physical meaning of the static friction coefficient, it can be determined by the angle of inclination of the shelf. A series of numerical experiments were carried out. In the first version, similar to [30], a box (without bottom) with particles was placed on a shelf, after which the shelf was slowly raised and the angle at which the box began to slide on the shelf was fixed. The second version involved the use of a counterweight for a box with particles standing on a horizontal plane. The counterweight force was slowly increased, and the force at which the box began to slide along the plane was recorded. However, the experiments carried out in the simulations did not show a correlation between SF PP , SF PB , and the angle of inclination of the shelf, as well as the force of the counterweight.
The approach was designed in such a way as to neutralize the main problems described in the introduction. For this, when developing, the authors started from the physical meaning of the coefficients, but transferring the meaning from a single particle to a portion of the investigated bulk material. It is possible to preserve the physical meaning of the calibrated DEM parameters, as well as reduce the number of possible parameter sets or even get a unique set. However, the approach has to include the determination of bulk material responses similar to the BCA method because of the problems with determining SF according to the physical meaning. This makes it possible to keep the methodology flexible and applicable to various DEM software and contact models.
The developed approach scheme is shown in Figure 2. For the investigated bulk material Young's modulus, Poisson's ratio, and particle size distribution are determined using well-known techniques [31][32][33][34][35]. A certain amount of material (portion) is poured into the developed test rig, and then its rheological properties are studied. The flow of the portion is recorded with a pre-installed high-speed camera. After that, the recording from the camera is sent to the computing device. A video processing algorithm developed using machine vision is used; values and parameters characterizing the rheology of the material are recorded. For DF coefficients, these parameters are flow times along the inclined shelves (Equation (1)); for CoR, these are angles and velocities of rupture and reflection (Equation (2)) [36]. Then, the experiments are repeated in the DEM software within the same conditions, but a series of numerical experiments are carried out, where the DEM Symmetry 2021, 13, 1088 4 of 13 parameters vary in a given range. For each experiment, an animation is recorded and then processed using similar machine vision algorithms. The experiments in the simulation are symmetrical to the full-scale experiments. A functional relationship is built between the parameters obtained from the simulations and DEM parameters. The obtained function depends on the particle size distribution and particle shape of the material, as well as on the applied contact model and its specific implementation in the DEM software. This means that the input data must be determined in advance. As a result, the obtained dependence is substituted with the values from a full-scale experiment, which makes it possible to obtain specific values of DEM parameters. In this case, the obtained values require clarification; therefore, the repose angle and the flow time are additionally measured. Then, in DEM software, DEM parameters are iteratively varied in a narrow range, which results in fairly accurate results of DEM parameters values.
into the developed test rig, and then its rheological properties are studied. The flow of the portion is recorded with a pre-installed high-speed camera. After that, the recording from the camera is sent to the computing device. A video processing algorithm developed using machine vision is used; values and parameters characterizing the rheology of the material are recorded. For DF coefficients, these parameters are flow times along the inclined shelves (Equation (1)); for CoR, these are angles and velocities of rupture and reflection (Equation (2)) [36]. Then, the experiments are repeated in the DEM software within the same conditions, but a series of numerical experiments are carried out, where the DEM parameters vary in a given range. For each experiment, an animation is recorded and then processed using similar machine vision algorithms. The experiments in the simulation are symmetrical to the full-scale experiments. A functional relationship is built between the parameters obtained from the simulations and DEM parameters. The obtained function depends on the particle size distribution and particle shape of the material, as well as on the applied contact model and its specific implementation in the DEM software. This means that the input data must be determined in advance. As a result, the obtained dependence is substituted with the values from a full-scale experiment, which makes it possible to obtain specific values of DEM parameters. In this case, the obtained values require clarification; therefore, the repose angle and the flow time are additionally measured. Then, in DEM software, DEM parameters are iteratively varied in a narrow range, which results in fairly accurate results of DEM parameters values. The test rig at the current stage of research is shown in Figure 3. The test rig has the shape of a rectangular parallelepiped measuring 1 × 1.5 × 0.13 m with inclined shelves located inside, vertical partitions 1-2, bins for loading bulk material, and a funnel-shaped device for testing the angle of rupture and repose, as well as a system of dampers for controlling the flow of materials. The test rig at the current stage of research is shown in Figure 3. The test rig has the shape of a rectangular parallelepiped measuring 1 × 1.5 × 0.13 m with inclined shelves located inside, vertical partitions 1-2, bins for loading bulk material, and a funnel-shaped device for testing the angle of rupture and repose, as well as a system of dampers for controlling the flow of materials. The flow of the material portion in the test rig occurs under the gravity force and is regulated by mechanical extraction of the partitions 1-2. Determination of CoRPP and DFPP requires the preparation of shelf 3 with a uniformly poured bulk material under study. The bulk material is poured into an adhesive (for example, epoxy resin, liquid nails, ceramic glue). Particles of bulk material are mechanically embedded in the adhe-  The flow of the material portion in the test rig occurs under the gravity force and is regulated by mechanical extraction of the partitions 1-2. Determination of CoR PP and DF PP requires the preparation of shelf 3 with a uniformly poured bulk material under study. The bulk material is poured into an adhesive (for example, epoxy resin, liquid nails, ceramic glue). Particles of bulk material are mechanically embedded in the adhesive. After a while, the substance solidifies and, as a result, a shelf is formed of particles of the bulk material under study. The control of the uniformity of the formation of the shelf is carried out visually. Figure 4 shows an analogue of the shelf used in the simulations. The flow of the material portion in the test rig occurs under the gravity force and is regulated by mechanical extraction of the partitions 1-2. Determination of CoRPP and DFPP requires the preparation of shelf 3 with a uniformly poured bulk material under study. The bulk material is poured into an adhesive (for example, epoxy resin, liquid nails, ceramic glue). Particles of bulk material are mechanically embedded in the adhesive. After a while, the substance solidifies and, as a result, a shelf is formed of particles of the bulk material under study. The control of the uniformity of the formation of the shelf is carried out visually. Figure 4 shows an analogue of the shelf used in the simulations. After the completion of the preliminary preparation of the test rig, the material portions are divided into weighed portions for each type of testing. At the beginning of the experiment, portions of the investigated bulk material of the same mass are poured into the upper hoppers. The weight of the samples is determined using a laboratory balance and is about 700 g with an error of no more than 0.5 g. After mechanical removal (of the dampers 2), the bulk material begins to slide along shelves 3 and 4, which is fixed by the algorithm as the beginning of the experiment. Shelves are angled at 40 degrees with a sliding path of 800 mm. This makes it possible to visually distinguish bulk materials with different dynamic coefficients of friction. At the end of the flow of material from each shelf, the algorithm separately fixes the moment in time and calculates the flow time. The experiment is repeated several times. The results are converted to DFPP and DFPB, respectively. After the completion of the preliminary preparation of the test rig, the material portions are divided into weighed portions for each type of testing. At the beginning of the experiment, portions of the investigated bulk material of the same mass are poured into the upper hoppers. The weight of the samples is determined using a laboratory balance and is about 700 g with an error of no more than 0.5 g. After mechanical removal (of the dampers 2), the bulk material begins to slide along shelves 3 and 4, which is fixed by the algorithm as the beginning of the experiment. Shelves are angled at 40 degrees with a sliding path of 800 mm. This makes it possible to visually distinguish bulk materials with different dynamic coefficients of friction. At the end of the flow of material from each shelf, the algorithm separately fixes the moment in time and calculates the flow time. The experiment is repeated several times. The results are converted to DF PP and DF PB , respectively.
The full-scale experiment continues, but the shelf with particles 3 and 4 moves to position 5, as a result of which the flow from the shelf ends with the collision of particles with walls 5 and 6. CoR PP and CoR PB parameters are determined using equation: where β re f l , β rupt angles of reflection and rupture of bulk material flow, V re f l , V rupt the speed of reflection and rupture at the point of impact, k-coefficient that depends on the specifics of the experiment and is taken into account at the stage of refinement of the initially obtained values (by default k = 1) [34]. As a refinement experiment for the investigated bulk material, the angles of repose and rupture are determined using device 7. Bulk material is poured into the upper part. After that, partition 8 is pulled out and the angle of rupture is fixed. Then, partition 9 is pulled out, and the material is poured onto the lower shelf of the test rig. As a result, the repose angle and time of flow from the funnel are recorded. The computer vision system is implemented using LabVIEW software tools. The image from the camera is taken in perspective. Using the reference points on the test rig, the image is projected onto a vertical plane. The original image (Figure 5a) then goes through several stages. It is first filtered and binarized using image thresholding (Figure 5b), then reconstructed using morphological image processing (Figure 5c). Further, Symmetry 2021, 13, 1088 6 of 13 depending on the task, there is a search for geometric primitives, or the hit of particles in the region of interest (ROI) is recorded. For example, when determining the angle of the repose (Figure 5d), the boundaries of the object are determined, after which a geometric primitive is built-a triangle. The resulting left and right angles of the triangle are averaged.
To determine the flow time, the first entry of particles into the ROI is recorded next to shelves 3 and 4. It is recommended to use a camera with a speed of at least 100 frames per second.
After that, partition 8 is pulled out and the angle of rupture is fixed. Then, partition 9 is pulled out, and the material is poured onto the lower shelf of the test rig. As a result, the repose angle and time of flow from the funnel are recorded. The computer vision system is implemented using LabVIEW software tools. The image from the camera is taken in perspective. Using the reference points on the test rig, the image is projected onto a vertical plane. The original image (Figure 5a) then goes through several stages. It is first filtered and binarized using image thresholding ( Figure  5b), then reconstructed using morphological image processing (Figure 5c). Further, depending on the task, there is a search for geometric primitives, or the hit of particles in the region of interest (ROI) is recorded. For example, when determining the angle of the repose (Figure 5d), the boundaries of the object are determined, after which a geometric primitive is built-a triangle. The resulting left and right angles of the triangle are averaged. To determine the flow time, the first entry of particles into the ROI is recorded next to shelves 3 and 4. It is recommended to use a camera with a speed of at least 100 frames per second.

Results
As an example of the work of the developed approach, a series of experiments were carried out. For this, three different types of metal steels with different degrees of surface roughness were used as the boundary material. Iron ore and waste rock were used as bulk materials in the experiments. The task was to obtain DEM parameters for two bulk materials and three steels (six boundary-material pairs). Physical implementation of the developed test rig is shown in Figure 6.

Results
As an example of the work of the developed approach, a series of experiments were carried out. For this, three different types of metal steels with different degrees of surface roughness were used as the boundary material. Iron ore and waste rock were used as bulk materials in the experiments. The task was to obtain DEM parameters for two bulk materials and three steels (six boundary-material pairs). Physical implementation of the developed test rig is shown in Figure 6.

Experiment Setup
Before the start of the experiments, according to the d shelves with particles fixed with glue were prepared for ore and

Experiment Setup
Before the start of the experiments, according to the developed methodology, shelves with particles fixed with glue were prepared for ore and waste rock (Figure 7).

Experiment Setup
Before the start of the experiments, according to the developed methodology, shelves with particles fixed with glue were prepared for ore and waste rock (Figure 7).  In addition, a DEM model was prepared for simulations in Rocky DEM software. The parameters of bulk materials and model are presented in Table 2. For simplicity, the particles are represented in the model as a 10-sided polyhedron ( Figure 8).
Symmetry 2021, 13, 1088 8 of 14 For simplicity, the particles are represented in the model as a 10-sided polyhedron ( Figure 8).

Determining DFPB and DFPP
The first step is to determine the dynamic coefficient of friction. For this, for each boundary-material pair, experiments were carried out to determine the time of sliding on an inclined shelf. Each experiment was performed three times. The measurement results are presented in the Table 3.

Determining DF PB and DF PP
The first step is to determine the dynamic coefficient of friction. For this, for each boundary-material pair, experiments were carried out to determine the time of sliding on an inclined shelf. Each experiment was performed three times. The measurement results are presented in the Table 3. Next, a series of simulations were carried out in Rocky DEM. The DF PB and DF PP coefficients varied from 0.1 to 0.9 with a step of 0.1. The flow time was recorded for each experiment. According to the data, functional dependence was built (Figure 9).  Next, a series of simulations were carried out in Rocky DEM. The DFPB and DFPP coefficients varied from 0.1 to 0.9 with a step of 0.1. The flow time was recorded for each experiment. According to the data, functional dependence was built (Figure 9). From the obtained dependences, the DF coefficients were evaluated, presented in Table 4. Figure 10 shows an example of a visual comparison of a full-scale experiment with an ore-steel 1 pair and simulation with a DFPB value of 0.3. Visually, it is noticeable that the flow time is almost identical. Moreover, according to the obtained dependence, the value of the coefficient for this pair is 0.29. From the obtained dependences, the DF coefficients were evaluated, presented in Table 4. Figure 10 shows an example of a visual comparison of a full-scale experiment with an ore-steel 1 pair and simulation with a DF PB value of 0.3. Visually, it is noticeable that the flow time is almost identical. Moreover, according to the obtained dependence, the value of the coefficient for this pair is 0.29.

Determining CoRPB and CoRPP
After the DF coefficients have been obtained for all pairs, with the CoR coefficient determination. For this, the shelves (s cally. The angles of rupture and reflection, as well as the veloc collision of particles with the shelves, were determined by sof

Determining CoR PB and CoR PP
After the DF coefficients have been obtained for all pairs, it is possible to proceed with the CoR coefficient determination. For this, the shelves (steel) are mounted vertically. The angles of rupture and reflection, as well as the velocities before and after the collision of particles with the shelves, were determined by software algorithms of the computer vision system and recalculated into the values of the CoR PP and CoR PB coefficients according to Equation (4). Simultaneously, numerical experiments were carried out with varying CoR coefficients in the range from 0.1 to 0.9. Similarly, the values obtained by Equation (4) made it possible to refine the values of the coefficients obtained in full-scale experiments (coefficient k in formula 4). Figure 11 shows an example of a collision in a full-scale experiment.
Based on the results, four coefficients were obtained for all pairs. The results are presented in Table 4.

Determining SF PB and SF PP . Refinement of Results
The SF PB and SF PP coefficients are selected iteratively during the simulations of material flow from the funnel and the formation of rupture and repose angles using the bisection method. This takes into account the fact that SF < DF in most of the use cases. In addition, an error is included in the values of the coefficients in Table 4. Full-scale experiments were carried out to determine bulk responses for ore and waste rock. After that, a series of numerical experiments was launched, in which DF PB and DF PP were varied in the range of obtained value ± 0.05, SF PB in [0.1; DF PB ], and SF PP in [0.1; DF PP ].
The specified maximum tolerance was no more than four degrees or 10% of the obtained value in a full-scale experiment. The obtained values of the repose angles in full-scale experiments and in the model after calibration are presented in Table 5.
collision of particles with the shelves, were determined by sof computer vision system and recalculated into the values of the cients according to Equation (4). Simultaneously, numerical ex out with varying CoR coefficients in the range from 0.1 to 0.9. S tained by Equation (4) made it possible to refine the values of the full-scale experiments (coefficient k in formula 4). Figure 11 show sion in a full-scale experiment.   Based on the results of the refinement experiment, the values of the DEM parameters were recalculated. The calibrated values are presented in Table 6. An example of an original image obtained from a video camera in a full-scale experiment, and the result obtained in a calibrated model, are presented in Figure 12.

Discussion
The developed approach was combined. Four coefficients (DFPB, DFPP, CoRPB, and CoRPP) were determined by methods based on their physical meanings. However, then, the SF coefficients were determined using a refinement test with the determination of the macro parameters of the bulk material, where the values of other coefficients vary within a narrow range. This made it possible to obtain a unique combination of DEM parameters that do not lose their physical meaning and adequately reproduce the behavior of bulk material in the DEM model.
The obtained results shown in Table 4 make it possible to judge that with different surface roughness, DFPB parameters obtained using the developed approach have regular differences. Thus, steel 1 and steel 2 have almost the same roughness (steel 2 is slightly larger), and steel 3 has a much higher roughness. At the same time, DFPB for steel 3-ore and steel 3-waste rock pairs is higher than for other pairs. The DF coefficients obtained after the refinement experiment (Table 5) did not change compared to the initial ones (Table 4). This was due to the fact that the regression algorithm in the refinement experiment achieved the desired result immediately after the selection of the SF parameters. With a higher required accuracy, the DF coefficients could change in the range of obtained value ± 0.05. In general, possible differences in the values of the parameters before and after the refinement experiment are associated with the inaccuracy of the motion measurement using a video camera, as well as the image processing algorithms.
Although the approach was developed for the study of absolutely any bulk materials, the specific implementation of the approach in the form of a physical test rig has a number of limitations. These limitations arise due to the specific design features and the capabilities and accuracy of the measuring devices. First, the test rig is designed for particles no larger than 15 mm. In the opposite case, material sticking can form in narrow places of the test rig, as well as possible errors in the representation of macro parameters (for example, the angle of rupture in the refinement tests). At the same time, materials with particles more than 15 mm can be crushed to the required size, and experiments can be carried out on a test rig. For this, it is necessary that the shape of the particles in the initial and final forms is the same. That is, for particles with specific shapes (ball, cube, etc.), this option is not suitable. It is also recommended to examine materials with particles of at least 1 mm in size. On the one hand, particles can seep through the slits in the structure; on the other hand, a low-resolution video camera may not detect their movement (in particular, this concerns image processing algorithms). Secondly, the image processing algorithms require refinement in order to improve the accuracy of the unambiguousness of the results obtained. It is planned to consider other algorithms for determining the repose angle of bulk material and improve the accuracy of determining the flow time.

Discussion
The developed approach was combined. Four coefficients (DF PB , DF PP , CoR PB , and CoR PP ) were determined by methods based on their physical meanings. However, then, the SF coefficients were determined using a refinement test with the determination of the macro parameters of the bulk material, where the values of other coefficients vary within a narrow range. This made it possible to obtain a unique combination of DEM parameters that do not lose their physical meaning and adequately reproduce the behavior of bulk material in the DEM model.
The obtained results shown in Table 4 make it possible to judge that with different surface roughness, DF PB parameters obtained using the developed approach have regular differences. Thus, steel 1 and steel 2 have almost the same roughness (steel 2 is slightly larger), and steel 3 has a much higher roughness. At the same time, DF PB for steel 3-ore and steel 3-waste rock pairs is higher than for other pairs. The DF coefficients obtained after the refinement experiment (Table 5) did not change compared to the initial ones (Table 4). This was due to the fact that the regression algorithm in the refinement experiment achieved the desired result immediately after the selection of the SF parameters. With a higher required accuracy, the DF coefficients could change in the range of obtained value ± 0.05. In general, possible differences in the values of the parameters before and after the refinement experiment are associated with the inaccuracy of the motion measurement using a video camera, as well as the image processing algorithms.
Although the approach was developed for the study of absolutely any bulk materials, the specific implementation of the approach in the form of a physical test rig has a number of limitations. These limitations arise due to the specific design features and the capabilities and accuracy of the measuring devices. First, the test rig is designed for particles no larger than 15 mm. In the opposite case, material sticking can form in narrow places of the test rig, as well as possible errors in the representation of macro parameters (for example, the angle of rupture in the refinement tests). At the same time, materials with particles more than 15 mm can be crushed to the required size, and experiments can be carried out on a test rig. For this, it is necessary that the shape of the particles in the initial and final forms is the same. That is, for particles with specific shapes (ball, cube, etc.), this option is not suitable. It is also recommended to examine materials with particles of at least 1 mm in size. On the one hand, particles can seep through the slits in the structure; on the other hand, a low-resolution video camera may not detect their movement (in particular, this concerns image processing algorithms). Secondly, the image processing algorithms require refinement in order to improve the accuracy of the unambiguousness of the results obtained. It is planned to consider other algorithms for determining the repose angle of bulk material and improve the accuracy of determining the flow time.
In general, the calculation time for all simulations was about 4 h on 10 cores of the average processor in the Rocky DEM. With several GPUs, calculations can be significantly sped up. This means that the calibration of DEM parameters could be a relatively quick and easy process.

Conclusions
This paper describes the author's approach to determining and calibrating the DEM parameters of a wide range of bulk materials. Features of the approach are identification of each DEM parameter accordingly to its physical meaning and conducting symmetrical experiments in full-scale and simulation. The approach itself can be applied to any bulk material. The test rig developed as an implementation of the methodology includes several sections, where each section is responsible for determining a specific DEM parameter. The use of a refinement test for the rupture and repose angles ensures that the material with the obtained unique set of DEM parameters adequately reproduces the material rheology in the DEM model. The presence of restrictions on the size and shape of particles impose specific features of the test rig design, as well as cameras for recording material movement. However, the test rig has wide applicability in various industries, from mining and metallurgy to pharmaceutical.
It is planned to further develop the project and conduct additional tests of other particle-material pairs in order to identify weaknesses in the test rig design and the approach as a whole. Other possible implementations of the approach will also be considered. The introduction and use of the approach and the test rig in the modeling of specific technological processes will increase the adequacy of the behavior of the material in the DEM models and simplify the solution of engineering problems using the discrete elements method software.