Influence of Inter-Particle Friction and Damping on the Dynamics of Spherical Projectile Impacting Onto a Soil Bed

This study investigates the dynamics of a spherical projectile impact onto a granular bed via numerical simulations by discrete element method (DEM). The granular bed is modeled as an assembly of polydisperse spherical particles and the projectile is represented by a rigid sphere. The DEM model is used to investigate the cratering process, including the dynamics of the projectile and energy transformation and dissipation. The cratering process is illustrated by tracking the motion of the projectile and granular particles in the bed. The numerical results show that the dynamics of the projectile follows the generalized Poncelet law that the final penetration depth is a power-law function of the falling height. The numerical results can match well the experimental data reported in the literature, demonstrating the reliability of the DEM model in analyzing the impact of a spherical projectile on a granular bed. Further analyses illustrate that the impact process consists of three main stages, namely the impact, penetration and collapse, as characterized by the evolution of projective velocity, strong force chains and crater shape. The initial kinetic and potential energy of the projectile is dissipated mainly by inter-particle friction which governs the projectile dynamics. The stopping time of projectile decreases as the initial impact velocity increases. The final penetration depth scales as one-third the power of total falling height and is inversely proportional to the macroscopic granular friction coefficient.


INTRODUCTION
The impact of projectiles on granular media are widespread phenomena in nature, such as asteroids colliding onto planetary surfaces (Senft and Stewart, 2009), raindrops falling into soil (Marston et al., 2010), people walking on sand beaches (Uehara et al., 2003) and rockfalls impacting onto soil buffering layers (Wang and Cavers, 2008;Calvetti and di Prisco, 2012;Su et al., 2018;Shen et al., 2019). The related research contributes to a better understanding of the formation of impact craters and the design of efficient shock absorbers. Though it has been studied via experiments, numerical simulations and analytical theories (Newhall and Durian, 2003;Wada et al., 2006;Crassous et al., 2007;Katsuragi and Durian, 2007;Clark et al., 2015;Li et al., 2016;Ye et al., 2016;Horabik et al., 2018;Ye et al., 2018;Zhang et al., 2021) in the past several decades, the understanding of impact process is still limited owing to the complexities of such events (Wada et al., 2006;Tiwari et al., 2014;Horabik et al., 2018).
The dynamic interaction between projectiles and granular media is complex because it depends on the kinematics of the projectile (e.g. size, density, velocity and shape) and the mechanical properties of the granular media (e.g. bulk density, inter-particle friction, inter-particle damping and particle size). A part of these factors have been systematically analyzed by a series of experimental tests (Newhall and Durian, 2003;Uehara et al., 2003;Katsuragi and Durian, 2013). The focus of these tests is on the ejection process, the crater morphology and the penetration dynamics, aiming to find a scaling law for crater size and penetration depth. For a spherical projectile impacting onto granular media, it is well established that the crater diameter scales with the power of 1/4 the kinetic energy of the projectile and the final penetration depth scales with the power of 1/3 the falling height (Uehara et al., 2003). The observed deceleration of the impacting sphere can be interpreted by the generalized Poncelet force law (Uehara et al., 2003). The stopping time of a spherical projectile has also been studied in different experiments (Ciamarra et al., 2004;Katsuragi and Durian, 2007). It was found to be a decreasing function of the impact velocity with an asymptotic plateau at high impact velocities. Newhall and Durian (2003) investigated the effect of projectile shape on the penetration of a projectile into a granular layer. Their results indicate that for low velocity impacts, the projectile shape plays a crucial role such that the sharper/ elongated objects can penetrate more deeply. Recently, with the help of photo-elastic particles and high speed camera, the microscopic inter-particle force networks in the granular media have been clearly identified by Clark et al. (2012). Based on the obtained results of force networks, the propagation and topology of the force networks were investigated (Clark et al., 2015;Takahashi et al., 2018). However, researchers are still far from a comprehensive understanding of the cratering mechanisms and the mechanical responses of the granular media (Wada et al., 2006;Tiwari et al., 2014), because it is difficult to control some physical parameters without changing other crucial parameters during the experiments. In particular, the influence of inter-particle friction and damping on the dynamics of projectile has not been investigated (Clark et al., 2015).
The aforementioned micro-and macro-mechanical responses of granular media can also be addressed by the discrete element method (DEM) (Cundall and Strack, 1979), especially for their loose and discontinuous natures (Wada et al., 2006;Li et al., 2016;Shen et al., 2018). With careful calibrations, the DEM modeling allows researchers to quantitatively analyze some physical processes that are nearly impossible to obtain in experiments (Utili et al., 2015;Zhang et al., 2015;Gao G. and Meguid M., 2018;Gao G. and Meguid M. A., 2018;Shen et al., 2021a), e.g. force wave propagation, energy evolution. In addition, the DEM allows a parametric analysis on one factor without altering other factors (Seguin et al., 2009). In the literature, both two-and three-dimensional (i.e., 2D and 3D) DEM models have been employed to study the projectile impacts on granular media (Tanaka et al., 2002;Ciamarra et al., 2004;Wada et al., 2006;Seguin et al., 2009;Kondic et al., 2012;Tiwari et al., 2014;Li et al., 2016;Horabik et al., 2018). Bourrier et al. (2010) and Zhang et al. (2017) investigated the impactinduced force chain evolution and its relation to the global mechanical response of granular media by 2D DEM. Wada et al. (2006) and Li et al. (2016) investigated the impact cratering processes of granular materials by 3D DEM models, reproducing the total mass and the velocity distribution of the impact-induced ejecta. The above studies show that the DEM is an effective and efficient method to investigate the impact of projectiles onto granular media. It is worth to note that the inter-particle friction and damping can be changed easily but not change the geometrical and mechanical properties of granular particles by using DEM.
The present study employs a 3D DEM model to investigate the dynamics of projectile impact onto granular media, with detailed parametric analyses of inter-particle friction and material damping. The paper is organized as follows: Section 2 presents a brief introduction of the DEM theory and the numerical model configurations. Section 3 summarizes the numerical results and discusses the dynamics of projectile and the energy transformation and dissipation during the impact process. Finally, conclusions on the projectile impact process are provided in Section 4.

Discrete Element Method Theory
The open-source DEM code ESyS-Particle (Wang and Mora, 2009;Weatherley et al., 2014) is employed to run all simulations presented herein. This code has been widely employed to study the mechanical behavior of soil and rock (Utili et al., 2015; Shen et al., 2018;Zhao et al., 2018;Shen et al., 2021b;Gao and Meguid, 2021). In DEM, the granular material is mimicked as an assembly of closely packed rigid spherical particles. The translational and rotational motions of each particle are governed by the Newton's second law of motion as expressed in Eqs 1, 2.
where F i is the resultant force acting on particle i; r i is the position of its centroid; m i is the particle mass; T i is the resultant moment acting on the particle; ω i is the angular velocity and I i is the moment of inertia.
In DEM, the inter-particle interactions are computed by the cohesionless frictional model. This employs a linear elastic-spring model to calculate the normal and tangential contact forces ( Figure 1). In order to replicate the energy dissipation by particle asperities being sheared off and the plastic deformations at contacts, a linear viscous damping model (dashpot model) is employed in the normal direction. Thus, the normal contact force (F n ) is calculated as, where u n is the overlapping length between two particles in the normal contact direction; k n is the normal contact stiffness and F d n is the normal damping force. The stiffness of particle normal contact is defined as k n πE(R A + R B )/4 with E being the Young's modulus of particles, R A and R B being the radii of the two contacting particles, respectively.
The damping force (F d n ) is calculated as where β is the damping ratio; m A and m B are the mass of the two contacting particles, and v n is the relative velocity between the particles in the normal direction. The shear force at the current time step (F n s ) is calculated incrementally as, where F n−1 s is the shear force at the previous iteration time step. ΔF s1 is calculated as △u s k s with k s being the shear contact stiffness and △u s being the incremental shear displacement. The shear stiffness is calculated as k s πE(R A + R B )/(8(1 + υ)) with υ being the particle Poisson's ratio. ΔF s2 is the shear force related to the rotation of particle contact plane. A detailed description of these two shear force terms can be found in Wang and Mora (2009).
The magnitude of the shear force is limited by the Coulomb's friction law as, where μ is the friction coefficient of particle. The shear induced moment is computed as:

Numerical Model Configurations
The numerical model configuration of a projectile colliding onto a granular bed is shown in Figure 2. The projectile is modeled as a rigid sphere of diameter (D) 2.54 cm. The granular bed is represented by an assembly of rigid spherical particles with radius uniformly ranging from 1.6 to 3.0 mm. It is prepared by randomly generating a loose packing of spherical particles in a cubic container space confined by rigid boundary walls. These particles are then released to settle downward under gravity until the kinetic energy of the system becomes nil. The acceleration due to gravity (g) is equal to 9.81 m/s 2 . After the gravitational deposition, the granular bed has dimensions of 22.1 cm in thickness (T), 20.0 cm in length (L) and width (W). The solid volume fraction is 0.58 and the bulk density (ρ g ) is 1.54 g/cm 3 . Consequently, it can be calculated T/D = 8.7, L/D = 7.9 and W/ D = 7.9. According to Seguin et al. (2008), the influence of model size on granular penetration can be neglected if T/D > 2.6, L/D > 5.0 and W/D > 5.0. Hence, the boundary conditions employed in the present study are acceptable. The input parameters of the DEM model are listed in Table 1. The material properties of the projectile and particles in the granular bed are chosen according to the well-documented experimental data reported in Katsuragi and Durian (Katsuragi and Durian, 2007;Katsuragi and Durian, 2013). For the granular particles, the interparticle friction coefficient (μ) is obtained by measuring the material angle of repose (θ r ) in the DEM to match approximately that of glass beads commonly used in experiments (≈22°). For the projectile, the Young's modulus and friction angle are set the same as those of the bed particles. The density of the projectile is set the same as that of a steel ball. During simulations, the projectile is positioned in the middle and just above the granular bed with an initial downward velocity (v 0 ). Each simulation lasts around 20 h on a standard desktop computer (Intel ® Core ™ i7 CPU, 4.00 GHz × 8, and 16 GB RAM).

RESULTS
The DEM model has been employed to investigate the dynamic response of the projectile and granular bed during the impact process. A series of simulations are computed under conditions of various initial impact velocities (v 0 ). The obtained numerical results are compared with the well-documented experimental and numerical data reported in the literature (Uehara et al., 2003;Katsuragi and Durian, 2007;Tiwari et al., 2014), in terms of the formation of crater, the dynamics of projectile, the penetration depth and the stopping time of projectile. In addition, the energy transformation and dissipation during the impact are analyzed quantitatively. Figure 3 illustrates the dynamic cratering process of the projectile colliding on and penetrating into the granular bed for test with v 0 = 5.0 m/s. For visualization purpose, the granular bed has been divided into seven equal-sized sub-layers with distinct colors at t = 0 ( Figure 3). To plot the velocity field, the bed has been divided into 20 × 21 equal-sized grid cells along a vertical plane and then the average velocity of all the particles in each grid cell can be obtained. From Figure 3, it can be observed that after t = 0, the projectile gradually penetrates the granular bed. The granular particles around the projectile are significantly disturbed. In the meantime, the deformation of granular bed is still small, and no obvious crater can be observed. As the simulation continues, particles beneath the projectile move downward and the side particles move laterally, forming a bowl-shaped crater. After a few milliseconds, the projectile enters the granular bed completely ( Figure 3). Meanwhile, the particles above the projectile continually move laterally and the diameter of the crater increases gradually (Figure 3). After the projectile penetrates into the granular bed, the top and lateral particles are still in dynamic motion with relatively small velocities (Figure 3). At  50 ms, the projectile stops moving, and an instantaneous deep crater is formed (Figure 3). Then, this transient crater begins to collapse. The collapsed particles gradually bury the projectile, forming a final stable crater (Figure 3). The numerical observations of the cratering process reveal three distinct stages: impact, penetration and side collapse. The impact stage is featured by the initial collision of a projectile onto the granular medium with very small bed deformation (e.g., Figure 3). The penetration stage is characterized by the rapid downward movement of the projectile in the granular media and the gradual expansion of the crater size (e.g., Figures 3a3-a5).

Formation of Crater
The collapse stage involves the toppling of the transient crater and the formation of a final stable crater (e.g. Figures 3a6, b6). These observations agree with some well-documented experimental and numerical results reported in the literature (Ciamarra et al., 2004;Wada et al., 2006;Li et al., 2016). The corresponding compressive force wave propagation within the impacted granular bed can be represented by the evolution of force chains, as shown in Figure 4. Here, the force chain is defined as a network of straight lines connecting the centers of contacting particles. The thickness of these lines is proportional to the magnitude of normal contact force. According to Bourrier et al. (2008), the propagation of force chains within an assembly of granular materials can be used to analyze the propagation of compressive stress waves. From Figure 4, it can be observed that before impact, the force chains are uniformly distributed within the granular bed, as determined by the gravity force ( Figure 4A). At impact, large contact forces occur immediately beneath the projectile and small contact forces distribute in the propagation front ( Figure 4B), indicating the propagation of compressive wave. Over time, a clear and intact force chain network is formed in the granular bed ( Figure 4C). After a certain time, the force chains stop propagating and begin to destruct ( Figures 4D, E). As the projectile penetrates into the granular bed, the force chain beneath the projectile is almost totally broken, but a small number of force chains still exist at the bottom edge of the projectile (Figures 4F-H). During the collapse stage, the large contact force chains (red lines) totally disappear ( Figures 4I, J). Figures 5A, B show the evolution of velocity (v) and penetration depth (y) of the projectile during the cratering process, respectively, for simulations of different initial velocities. The positive directions of the position and velocity are defined vertically down. Once the projectile impacts on the granular bed, the velocity firstly exhibits a rapid decrease and then slows down gradually to zero. According to Tiwari et al. (2014), the rapid decrease of the projectile velocity is due to the formation of strong force chain network in the granular bed. During the impact, the projectile must overcome the resistant force exerted by the bed materials (i.e. the force chain network), resulting in the rapid decrease of its velocity. In general, the time at which the projectile stops moving is defined as the stopping time (t s ). The stopping time decreases with the initial impact velocity. An impact at a high velocity exhibits a short stopping time. This phenomenon has also been observed by some experimental and numerical studies (Katsuragi and Durian, 2007;Seguin et al., 2009). For the evolution of projectile penetration depth ( Figure 5B), it increases gradually to the final penetration depth (h f ). As expected, it increases with the impact velocity. The impact at low velocity can lead to a relatively shallow crater, and specifically, the final penetration depth may be less than the projectile diameter (i.e. 0.0254 m). By contrast, the projectile with a relatively high speed can be entirely submerged in the granular bed, as implied in Figure 5B. According to Uehara et al. (2003), the final penetration depth of a projectile can be estimated by an empirical formula h f (0.14/μ′)(ρ p /ρ g ) 1/2 D 2/3 H 1/3 with μ′ being the macroscopic friction coefficient of bed particles and H being the falling height. μ′ is calculated as the tangent of the repose angle (θ r ). H is representative falling height calculated as the sum of the final penetration depth and the free fall height (H h f + v 2 0 /2g). Figure 6 shows the relationship between the final penetration depth and the total falling height from the current DEM simulations. The black curve shows the analytical results predicted by the empirical formula by Uehara et al. (2003). As shown in this figure, the numerical results can match well the empirical formula, demonstrating the accuracy of the empirical formula.

Projectile Dynamics
In the literature, several studies reported that the dynamics of a projectile can be assumed to follow the generalized Poncelet law (Katsuragi and Durian, 2007;Seguin et al., 2009) where h is the penetration depth of the projectile, m is the projectile mass; g is the gravity acceleration of value 9.8 m/s 2 ; d 1 is the material parameter of unit length; k is the elastic constant of contact. The solution of this equation shows the relationship between the velocity and penetration depth, as:    Figure 7 illustrates the relationships between the velocity and penetration depth for numerical simulations with different impact velocities. These curves exhibit the trend of a progressive change from a concave-down to a convex-up shape. Generally, there is a rapid decrease of velocity to zero at the end of penetration stage. These behaviors have also been observed in experiments by Katsuragi and Durian (2013) and DEM simulations by Tiwari et al. (2014). In addition, as shown in Figure 7, the numerical results can be well fitted by Eq. 9. This indicates that the dynamics of projectile in the present DEM study obeys well the generalized Poncelet law.
According to Katsuragi and Durian (2013), the fitting parameters d 1 and k in Eq. 8 can be approximately calculated as d 1 D 0.25 μ′ ρ p ρ g (10) k 12mgD −1 μ′ ρ g ρ p 1/2 (11) Thus, d 1 and k are estimated as 0.0831 m and 56.049 N/m, respectively. By substituting these values into Eq. 8, the final penetration depth and stopping time of the projectiles with different impact velocities can be predicted and compared with the numerical results ( Figure 8). It can be seen that the final penetration depths predicted by Eq. 8 using the parameters d 1 and k calculated from Eqs 10, 11 are larger than the numerical results. However, the increasing trend of the data obtained in numerical modeling is in accordance with the prediction of Eq. 8. The stopping time firstly shows a rapid decrease as the impact velocity increases from 0 to 1.0 m/s, and then it gradually approaches a stable value. For low impact velocities (v 0 < 1.0 m/ s), the stopping time can be well predicted by Eq. 8, while for higher impact velocities, Eq. 8 would underestimate the stopping time. However, the general decreasing trend of the stopping time matches well the prediction of Eq. 8. In addition, according to Katsuragi and Durian (2007), the ultimate stopping time can be estimated using D/g . This theoretical estimation is plotted as a gray dashed line in Figure 8B.

Energy Transformation and Dissipation
Analyses of the energy evolution, transformation and dissipation are important for a comprehensive understanding of the interaction process between the projectile and granular bed. The analyses also allow a quantitative evaluation of the role of inter-particle friction and damping on the projectile dynamics. The total energy (E T ) of the granular system consists of potential energy (E P ), kinetic energy (E K ), elastic strain energy (E S ) and the energy loss due to interparticle friction (E F ) and local contact viscous damping (E D ). The method to calculate these energy components in the context of DEM has been detailed in (Shen et al., 2018). In this study, all these energy components were recorded during the simulations and subsequently analyzed regarding the energy evolution, transformation and dissipation, as shown in Figure 9. The kinetic energy of the projectile (E p K ) and the granular bed (E g K ), the potential energy of the projectile (E p P ), the energy dissipated by friction between the projectile and granular bed (E p F ), the energy dissipated by damping (E g D ) and friction (E g F ) are recorded separately. The potential energy (E p P ) is defined with respect to the bottom of the granular bed. The strain energy (E S ) is not plotted because its value is negligibly small. All these energy components are normalized by the initial mechanical energy of the projectile (E 0 : the initial kinetic energy plus potential energy). As illustrated in Figure 9, the whole system involves a cascade of energy evolution which begins when the projectile impacts onto the granular bed. After the collision occurs (t = 0), the kinetic energy of the projectile decreases rapidly, while that of the granular bed increases rapidly to the peak value at t = 1.5 ms. During this period, only a small amount of energy is dissipated by inter-particle friction and damping. After t = 1.5 ms, the kinetic energy of both the projectile and granular bed decreases slowly. However, the decreasing rate of kinetic energy of the granular bed is smaller than that of the projectile. This is because during the penetration stage of the projectile, a certain portion of projectile kinetic energy can still be transferred to the granular bed.
To limit the computation time, the calculation is terminated when the projectile stops. Hence, at t = 51.8 ms, the kinetic energy of the projectile vanishes, while the granular bed still has a little of kinetic energy. From Figure 9, it also can be seen that nearly 70% of the mechanical energy of the projectile is dissipated by the inter-particle friction in the granular bed. The energy consumed by the friction of projectile and the damping of granular bed is less than 10% of E 0. Thus, from the perspective of energy consumption, it can be speculated that compared with the particle damping, the interparticle friction plays a dominant role in the mechanical response of granular bed against the projectile impact. However, in addition to the energy consumption, it is also necessary to study the influences of inter-particle friction and damping on the projectile dynamics.

Influence of Inter-Particle Friction on Projectile Dynamics
To investigate the influence of inter-particle friction (μ b ), a series of impact tests are performed with the bed particle friction coefficient (μ b ) ranging from 0.1 to 0.6. In these numerical tests, the granular bed is also generated by gravitational deposition as described in Section 2.2. Since the porosity of the granular bed depends on the inter-particle friction, the value of μ b is set as 0.1 during the gravitational deposition, such that the inter-particle friction can be varied to investigate its influence on the projectile dynamics without changing the initial stable state and especially the porosity of the granular bed. After the gravitational deposition, the volume fraction and bulk density (ρ g ) of the granular bed are 0.61 and 1.61 g/ cm 3 , respectively. Figure 10 shows the evolutions of velocities and penetration depths of the projectiles for the numerical tests with various interparticle friction coefficients. Generally, the velocity first shows a rapid decrease at the initial impact and then decreases gradually to zero. As the inter-particle friction increases, the stopping time of projectile decreases. For μ b ≥ 0.4, the evolutions of velocities are almost identical and the projectiles stop at the same time. As for the penetration depth, it increases gradually to a peak value after the collision is initiated. The final penetration depth decreases as the inter-particle friction increases. For the tests with μ b ≥ 0.4, the evolution of the projectile penetration depth is nearly identical and the projectiles almost stop at the same position. These results indicate that the interparticle friction indeed has a significant influence on the projectile dynamics. However, it has a relatively negligible influence on the velocity and penetrating depth of the projectile after the inter-particle friction is beyond 0.4. The dependence of the final penetration depth and stopping time on the friction coefficient of particles are illustrated in Figures 11A, B, respectively. It can be seen that both the final penetration depth and stopping time decrease as the inter-particle friction increases. After the inter-particle friction increases to 0.4, both the final penetration depth and stopping time become relatively stable. This is in accord with the general trend of the prediction by Eq. 8. Actually, such a phenomenon is related to the dependence of the macroscopic friction coefficient of the granular bed (μ′) on the microscopic inter-particle friction. For spherical particles, if particle rotation is allowed, the macroscopic friction (μ′) increases with the inter-particle friction (μ b ) and reaches the peak value after μ b increases to 0.45 (Suiker and Fleck, 2004). According to Eqs 10, 11, d 1 and k tend to be constant values when the macroscopic friction coefficient reaches its limit. Thus, the dynamics of projectile is not sensitive to μ b when it increases beyond 0.45.

Influence of Particle Damping on Projectile Dynamics
In this study, the linear dashpot model Eq. (4) is used in DEM to simulate the energy dissipation due to elastic wave propagation through a solid particle. This kind of model has been widely applied to numerical simulations of granular materials (Mao et al., 2004;Zhou et al., 2016;Gao G. and Meguid M. A., 2018;Li and Zhao, 2018;Shen et al., 2018). The damping ratio (β) in this model quantifies the plastic properties of particles. In the laboratory, the coefficient of restitution (COR) is commonly used as an index of this plastic property of solid materials. It can be measured in the laboratory from a series of drop tests (Imre et al., 2008). The relationship between the COR and the damping ratio (β) is depicted in Figure 12. The derivation of this equation can be found in (Gao G. and Meguid M., 2018). Thus, by varying the damping ratio, the influence of material damping of the granular bed on the projectile dynamics can be investigated. Figure 13 gives the evolution of velocity and penetration depth of the projectiles for the numerical tests at various damping ratios (i.e. 0.01, 0.1, 0.2 and 0.3). The corresponding COR values are 0.97, 0.73, 0.53 and 0.37, respectively. From Figure 13A, it can be observed that the damping ratio has little influence on the projectile velocity. For the tests with different damping ratios, the evolutions of velocity of the projectiles are nearly the same and the projectiles nearly stop at the same time. From Figure 13B, it can be observed that in the cases of β < 0.2, the final penetration depth decreases with the increase of the damping ratio. However, an increase of β from 0.01 to 0.2 leads to only 7.3% decrease of the final penetration depth. In addition, for the cases of β ≥ 0.2, the evolution of penetration depth of is nearly the same. Thus, it can be concluded that the main factor determining the dynamics of projectile is the inter-particle friction rather than the particle damping of a granular material.

CONCLUSIONS
The impact of a spherical projectile into a granular bed is analyzed via three-dimensional numerical modeling by the discrete element method. This model is validated by comparing the numerical and experimental results reported in the literature. The validated model is then used to investigate the influence of inter-particle friction and material damping on the dynamics of projectile. The main conclusions of this study are as follows: 1. Based on the numerical modeling, three key interaction stages, namely the impact, penetration and collapse, are identified. The impact stage is characterized by the initial collision of projectile onto the granular bed, the rapid decrease of projectile velocity and the formation of strong force chain networks. The penetration stage is featured by the movement of the projectile inside the granular bed and the crater expansion. The collapse stage involves the toppling of the deepest transient crater and the formation of a final static crater. 2. During the whole impact process, the initial kinetic and potential energy of the projectile is dissipated mainly by inter-particle friction, while only a small part of energy loss is induced by projectile-particle friction and inter-particle damping. 3. The parametric study on the inter-particle friction and damping shows that the inter-particle friction is the main factor determining the projectile dynamics.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.

AUTHOR CONTRIBUTIONS
TZ, GC, and FD contributed to conception and design of the study. WS performed the numerical simulations and organized the database. WS wrote the first draft of the manuscript. TZ, WS, GC, FD, and GD contributed to manuscript revision, read, and approved the submitted version.