1. Introduction
The discrete element method (DEM) was first developed by Cundal et al. (1979) [
1] and is a numerical technique for simulating the mechanical behavior of granular assemblies. It has the advantage of data trackings, such as the trajectories, velocities, and transient forces of all particles at any stage of the test. Therefore, the DEM has been successfully applied to simulate the bulk behavior of agricultural seeds, granular particles, and pharmaceutical tablets and to describe chemical reactions and heat transfer at the individual particle scale [
2]. The accuracy of bulk behavior simulation in the DEM depends on the DEM model and input properties, including the material properties and interaction properties [
3]. Farmers often use granular fertilizers to apply to the soil using a variety of machines. However, there is a lack of information on the granular fertilizer interaction properties to simulate in the DEM. Therefore, the first objective of this research was to calibrate particle-particle and particle-material interaction properties to investigate the granular fertilizer behavior with various devices.
A common DEM model for predicting the bulk behavior of non-cohesive materials is the Hertz–Mindlin model. However, this model is computationally time-consuming when many particles are generated [
4]. The Hertz-Mindlin model is used in conjunction with the Johnson-Kendall-Roberts (JKR) model for cohesive bulk materials if the tensional forces between particles have to be considered [
5]. Since the moisture content of granular fertilizer does not exceed 2%, it is regarded as a non-cohesive material. In this case, the Hertz-Mindlin (no-slip) model, a simplified form of the Hertz-Mindlin model, describes the cohesionless granular fertilizers [
6]. Therefore, in this study, the Hertz-Mindlin (no-slip) model was selected to simulate the particle-particle and particle-geometry interactions.
The granular fertilizer properties, such as the particle shape, size, particle size distribution, Poisson’s ratio, shear modulus, and density, are the primary input properties of the simulations. The objective of this paper is to validate the calibration strategy to determine the interaction properties of the granular fertilizers. Six interaction properties are required for the Hertz-Mindlin (no-slip) model: the coefficient of restitution, static friction coefficient, and rolling friction coefficient for the particle-particle and particle-material interactions. These interaction properties can be determined by measurements [
7], virtual calibration using the DEM [
8], or a combination of both [
9]. When the interaction properties of a single particle are measured with laboratory tests, the DEM simulation results do not match the experimental results because the physical meaning of the input properties is unknown due to the particle shape and size [
7]. If the bulk particle shape and size are accurately modeled by the DEM, the bulk behavior of the particles and their contacts can be accurately modeled [
10]. However, accurate modeling of the particle shape and size is challenging unless the particles are spherical [
11]. Nevertheless, the mismatch between the actual and simulated particle shape and size can be compensated for by performing a calibration of the interaction properties, which have a much stronger influence than the particle shape [
12]. Therefore, robust calibration procedures to determine the interaction properties (which are also called micro-properties) are required to reproduce the granular fertilizer behavior accurately so that the macro response simulated by the DEM matches the laboratory test results.
Researchers devoted to studying calibration methods of the particles have been developing laboratory tests for many years. A swing-arm slump test was developed by Grima and Wypych (2011) [
13] to calibrate the particle-particle static friction and rolling friction coefficients of polyethylene pellets to minimize the particle-material interactions. The swing-arm slump test creates a rapid flow of bulk material when the left and right arms attached to the symmetrically divided cylinder filled with polyethylene granules are instantly extended. Tekeste et al. (2018) [
14] used an anchor-Lift test to calibrate the interparticle static and rolling friction coefficients by measuring the repose angle of a conical corn pile; the anchor was lifted from a corn-filled cylinder to form the conical corn pile in the anchor. Derakhshani et al. (2015) [
15] and Do et al. (2018) [
16] performed a sandglass test to calibrate the interparticle static and rolling friction coefficients. Quartz sand in the upper transplant chamber flowed from the top part when the plug was removed from the sandglass neck. The discharge time (DT) and the angle of repose (AOR) were measured to determine the stabilization time of the sand particles, which formed a conical sand pile in the upper chamber. Similar to the sandglass test, the draw-down test was used to calibrate interparticle static and rolling friction coefficients of the coarse gravel by Richter et al. (2020) [
17], the corn grains by Coetzee (2020) [
9], and the particles with different diameters varying from 2 to 10mm by Ye et al. [
8], measuring the conical pile angle of repose of the upper box and the bottom box. The particles in the upper box of the draw-down test flowed to the bottom box after the outlet was opened, creating conical coarse gravel piles in the upper and bottom boxes. Ye et al. [
8] used a rotating drum test to calibrate the particle-particle and particle-material friction coefficients. The dynamic angles of repose of the upper and lower flows of the particles in a rotating 300mm diameter drum were measured. An Anton Paar powder cell was used by Salehi et al. [
18] to calibrate the interparticle interaction properties of powder by measuring the rotation torque of the impeller immersed in powder. Various impellers were used to compare the results. Cabiscol et al. [
19] used three laboratory tests: a pouring test to calibrate the interparticle friction coefficients; a rotating drum with a diameter of 130 mm to calibrate the particle-particle and particle-material sliding and rolling friction coefficients; and a drop test to calibrate the particle-particle and particle-material coefficient of restitution of cylindrical tablets. A cylinder with a diameter of 40 mm and a height of 255 mm was used in the pouring test. The cylinder was filled with the tablets at a constant speed of 3 mm/s. When the tablets formed a conical pile, the AOR between the horizontal surface and the pile’s lateral profile was measured. A single pharmaceutical tablet was released from the suction handle in the drop test from a height of 100 mm, and the rebound height was extracted with a high-speed camera (Promon 501, AOS Technologies, Switzerland) at an acquisition frequency of 1459 fps. The coefficient of restitution was determined by dividing the fall height by the bounce height. The drop test was used to determine the coefficient of restitution (COR) of a single particle. However, an average COR input value was used to simulate the bulk behavior using the DEM, even though the COR value depends on the particle’s weight, size and shape. The accuracy of calibrating the COR coefficient depends on the laboratory test device. In a pouring test with coke provided by Zu et al. (2018) [
20], the influence of the particle-particle COR was not significant because of the negligible influence of the particle-particle collision on the AOR of the conical coke pile. Nevertheless, the results of the draw-down test provided by Wei et al. (2017) [
21] and Geer et al. (2018) [
22] showed that the particle-particle COR affected the slope angle of the pile because of the high probability of the particle to slide and jump along the free surface. Xia et al. (2019) [
23] designed a slide plate test to calibrate the particle-material COR. The particles in a hopper fell onto the inclined collision plate located under the hopper and collided with other particles. The collision distance of the coal particle from the inclined collision plate and the AOR of the formed conical coal pile were measured.
However, there is a lack of information comparing the responses that best describe the interaction properties. Therefore, the second objective of this research was to design laboratory tests to calibrate the interaction properties and compare the responses best describing interaction properties. A potential problem of calibration is that more than one micro-property may influence the bulk response of the particles in the numerical experiment. It would be ideal if every input property could be determined separately in laboratory tests. However, it is challenging to measure the rolling and static friction of granular particles. A desirable strategy would be to calibrate the particle-particle interaction properties separately from the particle-material interaction properties using various laboratory tests to compare the macro responses. Performing separate calibrations of the particle-particle and particle-material interaction properties with two laboratory tests allows us to limit the number of factors influencing the interaction properties, thereby increasing the calibration accuracy. The calibration accuracy of the interaction properties depends on the design of the laboratory test and the experimental method.
The designed laboratory test is adequate only for granular non-cohesive materials with a particle diameter range of 1-5 mm. In this study, we determine the interaction properties of granular fertilizer particles for simulating applications with fertilizer metering devices.
2. The Experimental Framework
In this study, the granular fertilizer properties were determined experimentally, including the particle shape, size, and density. The material properties of the granular fertilizer, such as the Poisson’s ratio and shear modulus, were obtained from the literature because they have a negligible effect on the bulk behavior [
24]. The measured particle shape, size, and density were compared with the simulated values obtained from the DEM. If the bulk mass and volume match the experiment results, the particle-particle interaction properties are calibrated using the DEM. The sandglass test was modified to eliminate the particle-material interaction properties. After the calibration of the particle-particle micro-properties, a discharge-drop laboratory test was performed to calibrate the interaction properties of particles with different geometries (
Figure 1). The calibration accuracy of the particle-particle interaction properties has to be high because it affects the results of the particle-material interaction properties. Calibration is typically performed using a “trial and error” approach or optimization algorithms to reduce the number of simulations. Here, we adopt the central composite design (CCD) as an optimization algorithm.
2.1. Determination of the Shape, Size, and Density of the Granular Fertilizer
Choosing the appropriate particle size and shape is a critical part of the calibration process due to the substantial influence on the macro response of the bulk material. The particle shape influences the rolling coefficient parameters and the computational complexity. A dial caliper was used to measure three orthogonal dimensions of 80 randomly selected granular fertilizer particles [
25]. The largest of the three dimensions was designated as the length, the second largest was the width, and the smallest was the thickness. The sphericity of the particles was calculated using the following equation [
26]:
where
L,
W, and
T are the length, width, and thickness of the fertilizer granule (mm);
is the sphericity of the granule (dimensionless).
In addition to the particle shape, the particle size also significantly influences the macro response [
27]. Therefore, a 1000 g sample was used to determine the particle size and distribution using a sieve. The sieve opening diameter was changed from 6 to 1 mm with an interval of 0.5 mm. The granular fertilizer particle size ranged from 1.5 to 5 mm. The experiment was conducted in triplicate, and the results were averaged.
The liquid displacement method was used to determine the density of the granular fertilizer [
28]. A cylinder with an inner radius of 47.5 mm and a height of 45 mm was loaded with compound fertilizers, and the weight was determined with an electronic scale (accuracy of 0.01 g). Since granular fertilizers absorb water, a wetting agent (polyethylene glycol tert-octylphenyl ether) was added to the water at a concentration of 1.5 g/L. The wetting agent prevents the granular fertilizer from absorbing water by reducing the surface tension. The water was poured into the cylinder filled with the granular fertilizer and completely covered the granular fertilizer. The volume of the poured water was considered as free space between the granular fertilizer. The density of the granular fertilizer was determined as:
where
is the mass of the granular fertilizer in the cylinder, g;
is the volume of the cylinder, mm
3;
is the free space between the granular fertilizer particles, mm
3. The granular fertilizer density was 1575 kg m
−3. The granular fertilizer density depends on its moisture, which was 2.5%, with a standard deviation of 0.3%. The moisture of the granular fertilizer was measured before each test with a moisture analyzer.
2.2. EDEM Input Parameters
The simulation was carried out in the EDEM software. The values of the material properties used in the Hertz-Mindlin no-slip numerical model for the DEM simulations were obtained from the literature (except for the granular fertilizer density) (
Table 1). The values of the shear modulus and the Poisson’s ratio of the granular fertilizer obtained from the literature as the determined granular fertilizer density were identical to our experimental results. The particle size distribution has to be known to determine the particle-particle and particle-material interaction properties because it significantly affects the experimental results. The boundary conditions were not specified, since we were not interested in how the particles behave when leaving the domain. In the Hertz-Mindlin no-slip model the Euler was selected as a time integration method. The Raleigh time-step was 30%. The estimated cell radius of the simulator grid was 3 mm.
The low, mid and high levels of the interaction coefficients were determined to develop the CCD using the Design-Expert software (
Table 2). Choosing the low and high levels is essential. If the range between the minimum and maximum values is large, there may not be a sufficient number of points. If the range is small, some points may be outside of the range. In this study, the range was not large to ensure accurate calibration results. In this research, low and high levels were chosen according to pre-provided simulations in the DEM.
In the CCD, twenty configurations of the interaction properties were generated for the EDEM simulation using the sandglass test with three factors: the interparticle static friction coefficient (A), the interparticle rolling friction coefficient (B), and the interparticle restitution coefficient (C). The particle-material interaction parameters, i.e., the particle-material static friction coefficient (D), the particle-material rolling friction coefficient (E), and the particle-material restitution coefficient (F), were kept at the mid-levels. The interparticle interaction properties calibrated by the sandglass test were used to calibrate the particle-material interaction properties of the discharge-drop test. In the CCD, twenty configurations of the interaction parameters were generated for the EDEM simulation using the discharge-drop test with three factors: the particle-material static friction coefficient (D), the particle-material rolling friction coefficient (E), and the particle-material restitution coefficient (F). The optimization of the interaction properties according to the simulation results was conducted using the Design-Expert software.
2.3. Sandglass Test
2.3.1. Sandglass Test Process
The sandglass test was used with some modifications that did not take into account the particle-material interaction properties due to design (
Figure 2). First, a hole in the center of the bottom of the upper vessel (1) was surrounded by a ring with the same internal radius, with a height of 5 mm. The ring was added to prevent the bottom layer of granular fertilizers from moving when in contact with the material. The top of the ring was tilted outward at 45°. Second, the upper part of the inner vessel (4), which was attached to the cylinder (3) by a micro load cell (5), was tilted inward at 45° to prevent particles from interacting with the material. Therefore, the interaction of particles with material in the upper vessel and in the inner vessel was neglected, since the bottom was filled with immobile particles. Left and right plugs (2) were used to close the central hole of the upper vessel (1).
In the EDEM software, forty thousand particles were generated in the upper vessel in 4 s (
Figure 3a). The particles did not move until the left and right plugs were removed at a speed of 0.5 m/s, opening the hole in the center (
Figure 3b). The opening speed was 0.5 m/s because this was the speed of removing the plugs by hand in experiments, as determined by a high-speed camera. The particles were discharged from the upper vessel through the hole in the center and filled the inner vessel, and excess particles fell downward (
Figure 3c,d).
The sandglass test was designed to investigate the effect of the particle-particle interaction parameters on the responses, such as the discharge mass (DM) and the discharge time (DT) of the granular fertilizers from the upper vessel, and the AOR and the height (H) of the conical granular fertilizer pile. The sandglass test was repeated five times, and the simulations with the different configurations of the interaction properties were conducted once. The average of the five tests was the target value for the calibration of the EDEM interaction parameters.
2.3.2. Measurement of the Discharge Mass of the Granular Fertilizer
The
DM of the granular fertilizer was calculated by subtracting the final weight from the original weight. The original mass of the forty thousand particles simulated in EDEM in the upper vessel was 1837 g. Therefore, 1837 g of the granular fertilizer was loaded into the upper vessel in the real experiments. After discharging the granular fertilizer by removing the left and right plugs, the plugs were closed, and the granular fertilizer mass left in the upper vessel was weighed (
Figure 3d). The average of five replicates was the target
DM of the granular fertilizer for calibration of the interparticle properties. The fertilizer weight was measured with an electronic scale with an accuracy of 0.01 g.
2.3.3. Measurement of the Discharge Time of the Granular Fertilizer
The DT was measured using two methods, i.e., using a micro load cell (4) and a high-speed video camera (5). The results of the two methods were compared (
Figure 4). In the first method, the inner vessel (7) was connected to a micro load cell (4), which sent an analog signal to the Arduino Uno board (2) through the HX711 load cell amplifier (3). The Arduino Uno board converted the analog signal to a digital signal and sent the data to the serial monitor of the laptop (1), where the data and received data time were stored and later were processed to determine the
DT. The data in the monitor was zero when the right and left plugs blocked the central orifice of the upper vessel. The measurement time started when the right (9) and left plugs (10) were removed, and the granular fertilizers started to fall into the inner acrylic vessel (7). The measurement time ended when the granular fertilizer stopped discharging from the upper vessel. The data in the monitor were fluctuating when the particles were falling dawn and were constant when the particles stopped falling. The time when the data was fluctuating was considered the
DT. In the second method, as all pieces of the equipment (6–11) were acrylic, the discharge process could be monitored using a high-speed video camera (5), which captured images with an exposure time of less than 1000 s
−1. After recording, the images could be played back in slow motion, frame by frame. The measurement time started when the granular fertilizer began to discharge from the upper vessel and ended when the last granular fertilizer particle had left the upper vessel. The time between the start time and finish time was considered the
DT. The
DT results obtained from the load cell and the high-speed video camera were compared. The average of five replicates was considered the target value of the granular fertilizer’s
DT to calibrate the interparticle properties.
2.3.4. Measurement of the Angle of Repose of the Conical Granular Fertilizer Pile
Images of the four sides of the vessel were captured to determine the
AOR of the pile in Cartesian coordinates (
Figure 5a). It was determined that the area from 20 to 50 mm from the edge best described the
AOR of the conical granular fertilizer pile in the inner vessel. Therefore, sample points in this area were randomly chosen in the images and linearly fitted to determine the
AOR of the pile (
Figure 5b). The results showed that the area of interest for obtaining the
AOR was located at a distance of up to 30 mm from the edge of the vessel since it provided the best linear regression results. It should be noted that linear regression depends on the number of points sampled along the border. In this study, the number of sampling points was considered sufficient because the distance between the sampling points was 0.05 mm; thus, the slope of the line represents the
AOR of the particle pile (
Figure 6). It should be noted that all the conditions to determine the
AOR were identical.
Since the experiment was repeated five times and the pictures were taken from four sides, the average of twenty values of the AOR in the experiment was considered the target value to calibrate the interparticle interaction properties.
2.3.5. Measurement of the Height of the Conical Granular Fertilizer Pile
The five highest points from the top of the vessel were measured in Cartesian coordinates, and the average was used as the height of the pile (
Figure 7). The average was considered the target
H of the conical granular fertilizer pile for calibration of the interparticle properties.
2.4. Discharge-Drop Test
2.4.1. Discharge-Drop Test Process
The discharge-drop test was used to calibrate the particle-material interactions. The discharge-drop test equipment consisted of a container (1), drum (2), plug (3), collision plate (4), and collision plate attachment (5) located between the two acrylic plates (
Figure 8). The distance between the two acrylic plates was 50 mm. All the parts were acrylic and were attached to each other, except for the plug, drum, and collision plate attachment. The collision plate attachment was tilted 30 degrees from the horizontal. The drum and the collision plate can be made from various materials. In this research, polylactic acid (PLA) and acrylic were compared. The PLA drum and PLA collision plates were printed with a 3D printer, and the acrylic drum and acrylic collision plate were cut from an acrylic tube and acrylic board, respectively. The drum rotated clockwise at 10 rpm. The plug moved upward and downward, opening and closing the rectangular hole between the container and drum.
2.4.2. Measuring the Discharge Time of the Granular Fertilizer
When measuring the
DT, the high-speed camera was filming the upper left side of the discharge-drop test. The plug was closed when the container was empty (
Figure 9a) and was filled with 1640 g of granular fertilizer particles. In the EDEM, 35,800 particles were generated in the container to obtain a mass of 1640 g (
Figure 9b). Before the drum started to rotate, the plug was opened (
Figure 9c). In a real experiment, the plug was removed manually, and the opening speed did not affect the results of the experiment. In the EDEM, the plug was moved upward until the open hole was free. The drum started to rotate at 10 rpm, discharging the particles (
Figure 9d). The
DT started when the drum began to turn and ended when the level of the discharged granular fertilizer reached the open hole (
Figure 9e). In a real experiment, the moment when the discharged granular fertilizer reached the open hole could be monitored using a high-speed camera. The discharged granular fertilizer exiting the rotating drum fell by gravity onto the collision plate. The use of a fertilizer granular discharge drum to calibrate particle-material interaction properties is a novelty in this study.
2.4.3. Measuring the Collision Height
The discharged particles fell and collided with the inclined plate. The highest collision heights (
CH) were measured. The average of the measured 36 highest collision heights (
CH) was considered the target value to calibrate the particle-material interaction properties (
Figure 10). The advantages of using an inclined plate allow us to measure the height of the collision, since particles that collide with the collision plate bounce in one direction.