Numerical Investigation of the Characteristics of Erosion in a Centrifugal Pump for Transporting Dilute Particle-Laden Flows

: Erosion in centrifugal pumps for transporting ﬂows with dilute particles is a main pump failure problem in many engineering processes. A numerical model combining the computational ﬂuid dynamics (CFD) and Discrete Element Method (DEM) is applied to simulate erosion in a centrifugal pump. Different models of the liquid-solid inter-phase forces are implemented, and the particle-turbulence interaction is also deﬁned. The inertial particles considered in this work are monodisperse and have ﬁnite size. The numerical results are validated by comparing the results with a series of experimental data. Then, the effects of particle volume fraction, size, and shape on the pump erosion are estimated in the simulations. The results demonstrate that severe erosive areas are located near the inlet and outlet of the pressure side of the impeller blade, the middle region of the blade, the corners of the shroud and hub of the impeller adjoining to the pressure side of the blade, and the volute near the pump tongue. Among these locations, the maximum erosion occurs near the inlet of the pressure side of the blade. Erosion mitigation occurs under the situation where more particles accumulate in the near-wall region of the eroded surface, forming a buffering layer. The relationship between the particle size and the erosion is nonlinear, and the 1 mm particle causes the maximum pump erosion. The sharp particles cause more severe erosion in the pump because both the frequency of particle-wall collisions and the impact angle increase with the increasing sharpness of the particle.


Introduction
Transport of flows with dilute coarse particles by pumps is widely encountered in many industries, such as coal and mining, as well as metallurgical and chemical processes [1,2]. The collisions between hard solid particles and rotating hydraulic parts of the pump can cause surface damage or erosive wear, resulting in fracture, pitting, and even breakdown of the hydraulic parts. Recently, the activities of exploiting natural resources (e.g., oil and gas, geothermal energy, and minerals) from the deep earth and ocean require the development of multiphase flow pumps with high performance and abrasion resistance [3]. Large amounts of particles with different shapes, sizes, and properties will be conveyed by the liquids in the pumps. Compared with pipes and fittings, pumps are more vulnerable equipment in the conveying system because of the fierce collisions of particles and high rotating components. The pump erosion will become the primary concern for designing such systems in deep wells or undersea because repairing or changing pumps will be expensive. Poor design in the capability of resisting erosion may reduce the pump service life. Therefore, the mechanisms and characteristics of the erosion process in the pump caused by coarse particles should be well revealed in order to design new pumps for harsh environments.
Early studies of pump erosion were mainly based on erosion experiments that observed and measured material weight losses in the tested pumps [4][5][6]. However, a limitation of this method is that the results occur in real conditions with a vast number of influencing factors acting simultaneously, and the erosive mechanism of each factor is intractably studied from the experiment. Thanks to computational fluid dynamics (CFD), two methods based on the CFD technique are developed to model multiphase flows. In the Eulerian-Eulerian (E-E) framework, both the liquid and solid phases are solved as a continuum in the Navier-stokes equations. Several models were developed to model the inter-phase interactions with different assumptions and applicable scopes. Zhao et al. [7] simulated the liquid-solid two-phase flow in a double-channel centrifugal pump by using the mixture model (a special type of the E-E method) and found that the particle size and concentration have greater impacts on the erosion of the hydraulic components. Zhang et al. [8] investigated the characteristics of liquid-solid mixing flow in a centrifugal pump with different particle concentrations and diameters based on the mixture model. However, this model is not accurate when the particles are large and distributed unevenly in the fluid domain. In such conditions, the particles should be treated as individual elements. The Eulerian-Lagrangian (E-L) formulation is introduced in the CFD simulation. The particles are solved as discrete elements, and the particle motion is governed by Newton's Second Law. Lin et al. [9] used the discrete phase model (DPM) to simulate the slurry flow in the pipe and found good agreement with the experiment in terms of the distributions of particle concentration and velocity. Pagalthivarthi et al. [10] used the DPM method to simulate the wear on the volute of a pump transporting dilute slurry. Shen et al. [11] studied the influence of particle size on the wear rate and position in a centrifugal pump by using the DPM model. Huang et al. compared the effects of three different particle diameters on the wear position of the pump blades by using an unsteady calculation with the DPM method [12].
The existing drawbacks of the DPM method are the one-way coupling from fluid to solid and the absence of particle-particle collisions. The discrete element method (DEM) is then introduced and coupled with the CFD simulation, which is also viewed as the E-L method. Although the DEM simulation has been used in many engineering processes (e.g., fluidized bed, particle sediment) for a long time [13][14][15], the application of coupling CFD and DEM simulations in erosion prediction has been recently realized in both the pneumatic and hydraulic conveying systems [16,17]. The erosions of the bends with different diameters, curvatures, and orientations in the pneumatic system were predicted using the CFD-DEM simulations [18], and several experiments of the erosion in a 90 • bend were conducted to calibrate the numerical method [18,19]. Erosion in the bends, conveying liquid-solid flow, was also simulated using the CFD-DEM method [20]; the results were normally validated with experimental data measured by Zeng et al. [21]. There are few studies on erosion predictions in more complicated machines, such as pumps.
In fact, the accuracy of erosion predictions in the fluid-particle flow systems is determined by the predictions of various characteristics of the two-phase flow, such as particle distribution, particle velocity, and fluid-particle interactions in both momentum and turbulence. It is argued that the solid volume fraction should be considered in the fluid governing equations when the solid volume fraction is larger than 10 −6 [22]. Moreover, the fluid-particle interaction force should be considered in the fluid, which is known as the two-way coupling method [16]. Compared with that used in gas-solid flows, the lift force model is paid a special interest in the liquid-solid flow, and a recently developed model, i.e., the Loth model, is used in the simulation instead of the conventional Saffman and Magnus models [17]. The turbulence-particle interaction was ignored in the previous investigation but was considered using the discrete random walk (DRW) model in others [16]. However, the effect that the turbulent modification due to particles was not considered in the simulation or described in detail in the model. From the literature review, turbulence plays an important role in determining particle distribution and helping particle cluster formation [23][24][25][26]. Since turbulence in the centrifugal pump is usually very high, this effect on the particle distribution and associated erosion should be revealed.
In the present work, a delicate numerical model based on the CFD-DEM coupling is developed for modeling the dilute particle-laden flow and resulting erosion in a centrifugal pump. Since the turbulent intensity is normally high in pumps, the two-way coupling between turbulence and particles is considered in the simulation. The descriptions of the physical models of the pump and particle and the mathematic formulas of the numerical model are presented in Section 2. The authenticity of the numerical model is validated by comparing the simulated results with a series of previous experiments in Section 3. The effects of different factors, such as particle size, concentration, and shape, as well as the working flow rate on the erosion of different pump components are studied numerically. The main conclusions are summarized in the final section.

Description of Pump Model
The studied pump in this work was designed as a typical centrifugal pump, and the prototype of the pump can be referred to in [27]. The reason for reproducing their pump was to compare the predicted characteristics of the liquid-solid two-phase flows with the experimental results. The 3D geometry of the pump is presented in Figure 1, and the main parameters of the pump model are listed in Table 1. In order to simplify the geometry and save computational resources, the very small gaps between the impeller's shroud and pump casing were not simulated. Although the erosion in these gaps (due to recirculating flow mixed with small particles) plays an important role in the pump damage, most coarse particles can barely flow into the gaps. The particle motion and erosion in the flow passages of the hydraulic parts are of primary interest to us. ered in the simulation or described in detail in the model. From the literature review, turbulence plays an important role in determining particle distribution and helping particle cluster formation [23][24][25][26]. Since turbulence in the centrifugal pump is usually very high, this effect on the particle distribution and associated erosion should be revealed.
In the present work, a delicate numerical model based on the CFD-DEM coupling is developed for modeling the dilute particle-laden flow and resulting erosion in a centrifugal pump. Since the turbulent intensity is normally high in pumps, the two-way coupling between turbulence and particles is considered in the simulation. The descriptions of the physical models of the pump and particle and the mathematic formulas of the numerical model are presented in Section 2. The authenticity of the numerical model is validated by comparing the simulated results with a series of previous experiments in Section 3. The effects of different factors, such as particle size, concentration, and shape, as well as the working flow rate on the erosion of different pump components are studied numerically. The main conclusions are summarized in the final section.

Description of Pump Model
The studied pump in this work was designed as a typical centrifugal pump, and the prototype of the pump can be referred to in [27]. The reason for reproducing their pump was to compare the predicted characteristics of the liquid-solid two-phase flows with the experimental results. The 3D geometry of the pump is presented in Figure 1, and the main parameters of the pump model are listed in Table 1. In order to simplify the geometry and save computational resources, the very small gaps between the impeller's shroud and pump casing were not simulated. Although the erosion in these gaps (due to recirculating flow mixed with small particles) plays an important role in the pump damage, most coarse particles can barely flow into the gaps. The particle motion and erosion in the flow passages of the hydraulic parts are of primary interest to us.

Particle Modeling
Modeling of particle motion was carried out in the software EDEM, where the particle topology can be defined by assembling a number of elementary balls to approximately represent the particle configuration, as shown in Figure 2. The equivalent diameter of a non-spherical particle is determined to be equal to the diameter of a spherical particle with the same volume. Material properties of the particle can be defined in the software. According to [28], a parameter Ψ is defined to indicate the sphericity of the modeled particle. Sphericity is a measure of how close an object's shape is to a perfect sphere. It is defined as the ratio of the surface area of a spherical particle (having the same volume as a given particle) to the surface area of the particle and is calculated as: where V p is the volume of the particle and A p , the surface area of the particle. Three particle shapes are studied, and the corresponding Ψ values are 0.67, 0.84, and 1, respectively. The particle of Ψ = 0.67 is defined with relatively sharp vertices and edges, which are expected to cause the most severe surface erosion. The particle of Ψ = 0.84 has moderate vertices and edges and is more similar to the particle in the real working conditions. The particle of Ψ = 1 is a purely spherical shape and can represent an ideal situation for comparison.

Particle Modeling
Modeling of particle motion was carried out in the software EDEM, where the particle topology can be defined by assembling a number of elementary balls to approximately represent the particle configuration, as shown in Figure 2. The equivalent diameter of a non-spherical particle is determined to be equal to the diameter of a spherical particle with the same volume. Material properties of the particle can be defined in the software. According to [28], a parameter Ψ is defined to indicate the sphericity of the modeled particle. Sphericity is a measure of how close an object's shape is to a perfect sphere. It is defined as the ratio of the surface area of a spherical particle (having the same volume as a given particle) to the surface area of the particle and is calculated as: where Vp is the volume of the particle and Ap, the surface area of the particle. Three particle shapes are studied, and the corresponding Ψ values are 0.67, 0.84, and 1, respectively. The particle of Ψ = 0.67 is defined with relatively sharp vertices and edges, which are expected to cause the most severe surface erosion. The particle of Ψ = 0.84 has moderate vertices and edges and is more similar to the particle in the real working conditions. The particle of Ψ = 1 is a purely spherical shape and can represent an ideal situation for comparison.

Governing Equations
The CFD-DEM simulation consists of two coupled steps. The first step involves the modeling of particle motion using the DEM approach, and in the second step, the information of particle position and velocity is introduced into the CFD modeling and is coupled with the fluid simulation in the pump's fluid zone. After that, the forces due to the fluid flow imposed on the particle are returned to the DEM modeling of particle motion. A home-made interface code was written to couple the solvers of CFD and DEM and implement the associated inter-phase interactions, such as lift force, minor forces, and turbulence-particle interactions.

Governing Equations
The CFD-DEM simulation consists of two coupled steps. The first step involves the modeling of particle motion using the DEM approach, and in the second step, the information of particle position and velocity is introduced into the CFD modeling and is coupled with the fluid simulation in the pump's fluid zone. After that, the forces due to the fluid flow imposed on the particle are returned to the DEM modeling of particle motion. A home-made interface code was written to couple the solvers of CFD and DEM and implement the associated inter-phase interactions, such as lift force, minor forces, and turbulence-particle interactions.

Governing Equations of Liquid Phase
The RANS equations were used in modeling the liquid phase since they are still the most efficient model for engineering problems. The two-way coupling method is considered for the liquid-particle interaction. Therefore, the liquid volume fraction is solved in the governing equations, and the interaction between the liquid and particles were modeled as a body force in the conservation equations of momentum, which are expressed as: where ρ f is the fluid density, u i is the Reynolds-averaged velocity in the i direction, p is the pressure, µ is the fluid dynamic viscosity, and µ t , the turbulent viscosity. F p-f denotes the volumetric interaction force from the particles to the fluid, which equals the opposite of the total force exerted by all the particles in the considered CFD cell divided by the cell volume V cell . n is the number of particles in the considered CFD cell. The volume fraction α f is calculated as: where V p,i is the volume of the particle i. The turbulence was simulated by using the RNG k-ε model. A model different from the default model in the used commercial software was defined to simulate the turbulence modulation due to the particle presence since this effect may play an important role in determining the particle motion and causing the erosion in the bend. S k and S ε are the source terms of the turbulent kinetic energy and turbulent dissipation rate to represent the modifications in the turbulence by the particles and read [29]: where v f and v f are the average and fluctuating fluid velocity vectors, and v p denotes the instant particle velocity vector. f k is the drag factor, or the ratio of the drag coefficient to Stokes drag, which is defined as f k = C D Re p 24 . CD is obtained from Equations (12) and (13) and Re P is calculated from Equation (26). dp is the particle diameter. C ε3 = C ε2 C εp Re p 1.416 in which C ε2 ≈ C ε2 = 1.92 and C εp = 0.058 [29]. The source terms are defined in the interface code between Fluent and EDEM and linked to the turbulence equations by the user-defined function in Fluent.

Governing Equations of Solid Phase
The equations of DEM solved in EDEM are expressed as: where F drag is the drag force, F g is the gravitational force, F Lift is the lift force, F p is the pressure gradient force, F vm is the virtual mass force, F c is the contact force of particleparticle/wall, and F b is the buoyant force. Ω p is the angular velocity of a particle and I p = 0.1m p d 2 p is the moment of inertia for a sphere. T c is the contact torque, and T f is the torque caused by the fluid.

Drag Force Model
The particle drag force was calculated using the Di Felice [30] model, where the effect of the fluid volume fraction on the drag force is considered. It is thus expressed as: is the particle Reynolds number considering the fluid volume fraction. The drag coefficient C D of the spherical particle was calculated as follows: The drag coefficient of the non-spherical particle was calculated as follows [31]: where Re p,α is the particle Reynolds number, and Ψ is the Sphericity of the particle as defined in Equation (1).

Lift Force Model
The most widely used model of the lift force is the Saffman and Magnus lift force, which is developed from the gas-solid two-phase flow [32,33]. However, its validation in the simulation of liquid-solid flow is argued, although it was used in some works [34]. In this work, a more recent model developed specifically for liquid-solid flow is written in the home-made interface code to model the lift force in such flow. The Loth lift force is expressed as [35]: The lift coefficient reads: where the contained parameters are defined as follows [35]: where ω f is the fluid vorticity and Ω p , the particle angular velocity.

Models of Pressure Gradient and Virtual Mass Forces
The pressure gradient force and the virtual mass force are considered in the liquidsolid two-phase flow in that the densities of both phases are close. The virtual mass force can be calculated as: where C vm is the virtual mass force coefficient, which is equal to 0.5 as usual.
The pressure gradient force was also included in the model and is obtained using the following expression: The two forces were defined in the interface code in order to facilitate the on/off switch in the simulations.

Particle-Particle Collision Model
The particle-particle collision was modeled directly using the default model in the software EDEM. The contact force F c includes both the normal and tangential components, which were obtained through the linear spring-dashpot model proposed by Cundall and Strack [36].
where n and t stand for the normal and tangential unit vectors, respectively; k n,ij and k t,ij are the normal and tangential elastic stiffness; δ n,ij , δ t,ij are the displacements of colliding particles, and γ n,ij , γ t,ij are the damping coefficients. v r = v i − v j is the relative velocity vector of particle i and particle j; r is the vector from the mass center to the contacting point. Notice that this model can equally be applied to describe the interaction between a particle and a wall. The equations for calculating the k, δ, and γ can be found in the literature [17]. It is worth mentioning that the tangential contact force can also be calculated by the Coulomb friction model when sliding friction occurs: where µ s is the sliding friction coefficient.

Turbulence-Particle Interaction Model
The influence of turbulence on the particle motion can be simulated by using the discrete random walk (DRW) model [37]. The turbulent velocity fluctuation obeys the Gaussian distribution and is obtained via the following relations: where ξ i , ξ j , ξ k are random numbers obeying the Gaussian distribution. Due to the simplicity of the RANS model, the turbulence was assumed isotropic, and the root mean square (RMS) of the fluctuating velocity was calculated as: with k being the turbulence kinetic energy, obtained using the RNG k-ε model. The randomly fluctuating velocity was then used as an input into the expression for the fluid velocity, and the turbulence effect was also considered when calculating the particle motion. The particle-turbulence interaction is two-way coupling using Equations (6), (7), and (33) and yields more accurate results.

Erosion Model
Zhang's model [19] obtained from controlled experimental tests for solid particles carried in a liquid flow was used here. The correlation proposed by Zhang et al. is: where ER is the erosion ratio, defined as the amount of mass lost by the wall material due to particle impacts divided by the mass of the particles impacting; BH is the Brinell hardness of the wall material; F s is the particle shape coefficient, F s = 1.0 for sharp (Ψ = 0.67), 0.53 for semi-rounded (Ψ = 0.84), or 0.2 for fully rounded (Ψ = 1) sand particles; V p is the particle impact speed in m/s; θ is the impact angle in radians; n = 2.41 and C = 2.17 × 10 −7 are empirical constants. Values of A i for i = 1-5 are 5.40, −10.11, 10.93, −6.33, 1.42, respectively.

Boundary Conditions and Model Setups
The commercial CFD software ANSYS Fluent was used for liquid modeling, and the software EDEM was coupled with ANSYS Fluent. The RNG k-ε turbulence model was selected for modeling the turbulence in the bend and pump because of its correction for swirling flow in the model [38]. The mesh in the impeller region was set to rotate with the angular velocity of the pump revolution speed. The flow information was exchanged on the interface defined on the meshes of the static and rotating zones. A mass flow rate of 4.17 kg/s (Q d = 15 m 3 /h) was defined at the inlet of the pump model to simulate the nominal flow rate, and the outlet is set as free outflow because of the long extended pipe at the pump outlet. The turbulence parameters at the inlet can be defined as 5% for the turbulent intensity, and the hydraulic diameter is 0.065 m. A straight pipe with 8D length was created in front of the pump inlet to allow a fully developed inflow. The SIMPLEC algorithm was chosen to couple the velocity and pressure terms, and the second-order upwind scheme was employed in the discretizations of pressure, momentum, turbulent kinetic energy, and specific dissipation rate.
In EDEM, all surfaces were defined as stainless steel, of which the Poisson's ratio is 0.3, Young's modulus is 2 × 10 11 Pa, and the density is 8000 kg/m 3 . For the particleparticle collision, the coefficient of restitution is 0.95, the coefficient of static friction is 0.004, and the coefficient of rolling friction is 0.5 [34]. For the particle-wall collision, the coefficient of restitution is 0.74, the coefficient of static friction is 0.2, and the coefficient of rolling friction is 0.3 [34]. Since the impeller rotates at high speed, the contact model chosen in the simulation was the Hertz-Mindlin (no slip). The particles are uniformly introduced in the extended pipe from the pump inlet and removed when they leave out the CFD zone. Depending on the simulated particle diameter or particle volume (for non-spherical particles), the equivalent numbers of particles per second are defined in the zone of particle creation in order to obtain the simulated particle volume fractions. The model of erosion prediction was plugged in the EDEM through the user-defined program for specific purposes. A home-made interface code was developed to couple the calculations between the two tools. Except for the drag force, the various minor forces, i.e., Loth lift force, pressure gradient force, and virtual mass force, were all defined in the code to calculate the particle motion. The opposite force of the total force acted on each particle was summed through all of the particles present in the same CFD cell, and the summing force was added into the momentum equation of the fluid phase as the source term. In addition, the influence of the turbulence on the particle motion was also considered by decomposing the fluid velocity into the RANS-based mean velocity from the results of the CFD solver and the fluctuating velocity sampled by the DRW model. At the same time, the influence of particles on fluid turbulence was also considered by adding appropriate source terms as defined in Equations (6) and (7) in the transport equations of turbulence.
Transient simulation is performed in both Fluent and EDEM. The time step for the liquid modeling was 1 × 10 −4 s. The time step in EDEM should be defined as 20-40% of the Rayleigh time calculated in EDEM. Therefore, the time step in EDEM was equal to 1 × 10 −5 s. The total simulation time was 1.33 s, and the impeller accomplished 20 revolutions. The default criterion was set for the simulation convergence in each time step, i.e., 10 −4 for mass, momentum, and turbulent parameters, and the maximum iteration number of each time step was 40. Each simulation takes approximately 100 h with 6 cores of 2.1 GHz.

Mesh Verification
Three meshes with different element numbers were calculated, and several parameters at the monitored points (MPs) were compared. The results are all presented in Table 2. It is noticed that the discrepancy of all compared parameters is relatively large between mesh #1 and mesh #2, and this discrepancy became marginal between mesh #2 and mesh #3. Therefore, mesh #2 with 1.3 million elements was used in the following study. The chosen mesh is shown in Figure 3.

Validation of Liquid-Solid Two-Phase Flow in the Pump
The accuracy of the CFD-DEM model for predicting particle motion in the pump, especially in the rotating region of the pump, should be verified with experimentation. Shi et al. [27] proposed a measuring method based on PIV and successfully measured the two-phase velocity fields of a water-glass bead mixture and a water-rapeseed mixture in a centrifugal pump. Their results were used as a reference for the quantitative verification of the particle motion in the pump predicted by the CFD-DEM model in this work. The centrifugal pump model is reproduced in a CFD model, as shown in Figure 1. A straight pipe with 8D length was added at the inlet of the impeller for fully developing the mixture of fluid and particles. To be consistent with the experimental parameters, the mass flow rate at the inlet was 4.17 kg/s, and the input particle number was 3,184,713 per second, which corresponds to the 500 µm particles with 5% volume fraction. The impeller speed was 900 r/min, and free outflow was defined at the model outlet. The particle density was 2500 kg/m 3 , and the Poisson's ratio and shear modulus of the particles were 0.4 and 1.3 × 10 8 Pa, respectively. A Hertz-Mindlin (no-slip) model for the collisions of particleparticle and particle-wall was used, and the time step in EDEM was set as 1 × 10 −5 s. The simulation was performed for 1.33 s, corresponding to 20 revolutions. Figure 4a,b shows the contours of absolute water velocity in the zone of the PIV camera between the experimental and simulation results. The visualized plane is located in the middle of the blade span. The fluid was clean water. Figure 4b shows the velocity field on both sides of the blade when the impeller was turned at the same position as the PIV shoot. It demonstrates that the CFD result is consistent with the PIV result. Less than 0.1 m/s difference was found at the compared positions. Figure 5 shows the fluid and particle velocities in the vicinity of the pressure side of a blade. The glass beads were 500 µm in diameter and 5% in volume fraction. The scattered symbols show the PIV measured velocities of the two phases, while the simulation results are represented in solid lines. It clearly illustrates that the simulation results are in good agreement with the experimental results, and the tendency of the velocity variations of the two phases along the blade length was well captured by the CFD-DEM model. diameter and 5% in volume fraction. The scattered symbols show the PIV measured velocities of the two phases, while the simulation results are represented in solid lines. It clearly illustrates that the simulation results are in good agreement with the experimental results, and the tendency of the velocity variations of the two phases along the blade length was well captured by the CFD-DEM model.

Validation of the Erosion Prediction in a 90° Bend
Erosion in a 90° bend was experimentally studied by Zeng et al. [21]. The diameter of the bend was D = 50 mm, and the ratio of bend radius to its diameter The conveying particle concentration was 1.2 wt% (and approximately 0.5% in fraction), and the particle size was within 400-500 μm. The simulated bend and m shown in Figure 6a,b. A horizontal pipe of 20 D length was added in front of t The downstream vertical pipe was 10 D in length, and constant pressure was de the boundary condition at the outlet. The developed numerical model was applie ulate the erosion process in the bend by using the spherical particles with a 450 μ eter and 0.5% in concentration. The fluid velocity was set as 4 m/s at the inlet. In t the RNG k-ε turbulence model is employed since a strong secondary flow occu bend. The Hertz-Mindlin (no-slip) contact model was chosen in the particle sol the erosion model was defined in the module of the particle-wall contact model u user-defined programs for specific purposes. The parameters employed in the sim are summarized in Table 3.

Validation of the Erosion Prediction in a 90 • Bend
Erosion in a 90 • bend was experimentally studied by Zeng et al. [21]. The internal diameter of the bend was D = 50 mm, and the ratio of bend radius to its diameter was 1.5. The conveying particle concentration was 1.2 wt% (and approximately 0.5% in volume fraction), and the particle size was within 400-500 µm. The simulated bend and mesh are shown in Figure 6a,b. A horizontal pipe of 20 D length was added in front of the bend. The downstream vertical pipe was 10 D in length, and constant pressure was defined as the boundary condition at the outlet. The developed numerical model was applied to simulate the erosion process in the bend by using the spherical particles with a 450 µm diameter and 0.5% in concentration. The fluid velocity was set as 4 m/s at the inlet. In this case, the RNG k-ε turbulence model is employed since a strong secondary flow occurs in the bend. The Hertz-Mindlin (no-slip) contact model was chosen in the particle solver, and the erosion model was defined in the module of the particle-wall contact model using the user-defined programs for specific purposes. The parameters employed in the simulation are summarized in Table 3.
eter and 0.5% in concentration. The fluid velocity was set as 4 m/s at the inlet. In this case, the RNG k-ε turbulence model is employed since a strong secondary flow occurs in the bend. The Hertz-Mindlin (no-slip) contact model was chosen in the particle solver, and the erosion model was defined in the module of the particle-wall contact model using the user-defined programs for specific purposes. The parameters employed in the simulation are summarized in Table 3.
(a) (b) Figure 6. Schematics of the simulated 90° bend (a) and the mesh (b) employed in the simulation. Figure 7 shows the predicted erosion rates along the curves at different circumferential positions, as depicted in Figure 6a. The curve with the azimuthal angle of 180° corresponds to the extrados of the bend. The simulated erosion rates are compared with Zeng's experiment. It is clearly demonstrated that the E/CRC's model predicted the erosion rates accurately. The E/CRC erosion model is therefore employed in the following simulations.

Validation of the Erosion Prediction on the Hydraulic Components of Pump
Zhang's erosion model was proved to predict a reasonable erosion rate in the bend. However, there is a lack of experimental data for the erosion on hydraulic parts of pumps in complex flow conditions. Fortunately, the studied pump was operated under circumstances very similar to that in Zeng's bend experiment: the particle diameter was 500 µm, the particle volume fraction was 0.5%, and the average fluid velocity in the impeller was approximately 4 m/s. The predicted erosion rate on the hydraulic parts of the pump was compared with that in the bend. The computational domain and boundary condition settings of the pump model were the same as those described in Section 2.5. The simulated particles are of 500 µm and 0.5% volume fraction. The flow rate was 15 m 3 /h, at which the average fluid velocity in each channel of the pump impeller is approximately 4 m/s. The predicted erosion rate is presented in the unit of mm/y. Figure 8 shows the annual erosion rates along the pressure side of one blade of the pump impeller and the extrados of the bend in Zeng's experiment. The X-axis indicates the non-dimensional position on both entities. The full length of the blade and bend extrados is represented by 1.0. It was seen that the predicted erosion rate on the impeller blade is in the same order of magnitude as the experimental value under similar particle and hydrodynamic conditions. The severely erosive areas were located near the inlet and outlet of the blade, as well as the middle region. The maximum erosion rate was very close to that in the bend erosion. The results validate the numerical model in predicting the erosion on the hydraulic parts of the pump.

Influence of Particle Concentration on the Pump Erosion
The influence of the concentration of loading particles on pump erosion is e in this section. The shape of the simulated particles is spherical, and the particle d is 1 mm. At the inlet of the model, the number of simulated particles is appro defined in order to obtain the simulated volume fractions of 0.1%, 0.5%, and 2.5% tively. The former two values are often found in the transporting medium of mining engineering and chemical engineering, and the latter one occurs in the h transporting system of liquid-solid mixed flows. Since the modeled pump posse typical configuration of a centrifugal pump, the simulated results of pump eros guide the design and optimization of anti-erosion centrifugal pumps.
The particles are tracked in a specific flow passage of the pump under the si particle concentrations. At the same time, the influence of the working flow rat estimated at 0.6Qd, 1.0Qd, and 1.4Qd. The results are shown in the figure below. T tures of particle transport can be obtained from the perspective of particle dist Part of the particles at the impeller inlet are decelerated, which is attributed to with the different components of the impeller as well as particle-particle collisi blocking effect becomes significant when the particle concentration increases Once the particles come into the passage, they move along the pressure side of th At the small particle concentration, the particles cannot always be constrained on sure side; they flux into the bulk flow at the outlet of the passage. As the particle

Influence of Particle Concentration on the Pump Erosion
The influence of the concentration of loading particles on pump erosion is estimated in this section. The shape of the simulated particles is spherical, and the particle diameter is 1 mm. At the inlet of the model, the number of simulated particles is appropriately defined in order to obtain the simulated volume fractions of 0.1%, 0.5%, and 2.5%, respectively. The former two values are often found in the transporting medium of coal and mining engineering and chemical engineering, and the latter one occurs in the hydraulic transporting system of liquid-solid mixed flows. Since the modeled pump possesses the typical configuration of a centrifugal pump, the simulated results of pump erosion will guide the design and optimization of anti-erosion centrifugal pumps.
The particles are tracked in a specific flow passage of the pump under the simulated particle concentrations. At the same time, the influence of the working flow rate is also estimated at 0.6Q d , 1.0Q d , and 1.4Q d . The results are shown in the figure below. Two features of particle transport can be obtained from the perspective of particle distribution. Part of the particles at the impeller inlet are decelerated, which is attributed to impacts with the different components of the impeller as well as particle-particle collisions. The blocking effect becomes significant when the particle concentration increases to 2.5%. Once the particles come into the passage, they move along the pressure side of the blade. At the small particle concentration, the particles cannot always be constrained on the pressure side; they flux into the bulk flow at the outlet of the passage. As the particle concentration increases, more particles are constrained on the pressure side and move out from the outlet of the blade. It can be expected that erosion on the surfaces of different components of the impeller near the outlet is increased with an increasing particle concentration. The particles in the flow passage are also tracked at different flow rates, and the particle distribution is very similar, except that the number of particles in the passage increases with increasing flow rate. The corresponding figures are ignored for the sake of brevity.
The erosion rates on the surfaces of different components of the pump are shown in the figures below for three particle concentrations. On the surface of the impeller shroud (as shown in Figure 9a), the erosion mainly occurs near the corner adjoining to the pressure side of the blades. The erosion rate is increased as the particle concentration increases. The maximum annual erosion rates are 0.89, 6.96, and 153.50 mm/year under the particle volume fractions of 0.1%, 0.5%, and 2.5%. It is shown the erosion is not distributed evenly in different flow passages at the low particle volume fraction. However, the distribution of the erosion becomes more uniform among the different passages with the increasing volume fraction. The erosion on the impeller hub is shown in Figure 9b. The maximum annual erosion rates are comparable to those on the impeller shroud. The regions of severe erosion mainly occur near the inlet of the impeller. This region is elongated along the streamwise direction and near the corner adjoining to the pressure side of the blade with the increasing particle volume fraction. The erosion on the impeller blades is shown in Figure 9c. It is seen that severe erosion mainly occurs on the pressure side of the blades. The region near the impeller inlet is always eroded by the incoming particles, and the region is extended over one-third of the blade length. The erosion in the region near the blade outlet is significantly increased with the increasing particle volume fraction because more particles are accumulated in the near-wall region of the pressure side, and more collisions occur between the particles and wall. It is found that the erosion rate is maximum on the blades than on the surfaces of the other components. The erosion on the volute is shown in Figure 9d. The results show that severe erosions mainly occur on the internal surface of the volute near the tongue. The tongue plays the role of blocking the conveying particles out of the volute and increasing the collisions between the rebounded particles and the internal surface of the volute, resulting in an augmented erosion.
The influence of working flow rate on the pump erosion is also estimated for the different particle concentrations, as shown in Figure 10. The contours of erosion rate on different components are replaced by the maximum annual erosion rate on each component. The simulated flow rates are 0.6Q d , 1.0Q d , and 1.4Q d , respectively. It is found from Figure 10a that the maximum erosion occurs on the pressure side of the blades at the low (0.6Q d ) and nominal (1.0Q d ) flow rates under 0.1% particle volume fraction. The maximum annual erosion rate is 1.12 mm/year and 1.97 mm/year, respectively. At the large (1.4Q d ) flow rate, the maximum erosion occurs on the internal surface of the volute with an annual erosion rate of 3.45 mm/year. The increasing rate of erosion is largest on the volute surface with the increasing flow rate, as seen in Figure 10a. When the particle volume fraction increases to 0.5%, the maximum erosion occurs on the blades at all simulated flow rates. At the low and nominal flow rates, the maximum erosion rates on the shroud and hub of the impeller and the volute are comparable, while at the large flow rate, the erosion on the volute is significantly increased to the level comparable to that on the blades. When the particle volume fraction increases to 2.5%, the maximum erosion occurs on the blades at all simulated flow rates, and the annual erosion rate exceeds over 100 mm/year, which indicates that the hydraulic part of the pump cannot work normally in such harsh conditions. The general features of the erosion on the pump components can be concluded as follows. The maximum erosion always occurs on the pressure side of the blades at all simulated scenarios except for the case of 0.1% and 1.4Q d , in which the erosion on the volute is larger. The rate of increase in the erosion rate is largest on the internal surface of the volute with respect to the flow rate. The erosion rates on the shroud and hub of the impeller are comparable at all simulated conditions. The erosions that occurred in the extended pipes of the inlet and outlet are much smaller than the other components.
collisions occur between the particles and wall. It is found that the erosion rate is maximum on the blades than on the surfaces of the other components. The erosion on the volute is shown in Figure 9d. The results show that severe erosions mainly occur on the internal surface of the volute near the tongue. The tongue plays the role of blocking the conveying particles out of the volute and increasing the collisions between the rebounded particles and the internal surface of the volute, resulting in an augmented erosion. (d) Volute The influence of working flow rate on the pump erosion is also estimated for the different particle concentrations, as shown in Figure 10. The contours of erosion rate on different components are replaced by the maximum annual erosion rate on each component. The simulated flow rates are 0.6Qd, 1.0Qd, and 1.4Qd, respectively. It is found from tions. The general features of the erosion on the pump components can be concluded as follows. The maximum erosion always occurs on the pressure side of the blades at all simulated scenarios except for the case of 0.1% and 1.4Qd, in which the erosion on the volute is larger. The rate of increase in the erosion rate is largest on the internal surface of the volute with respect to the flow rate. The erosion rates on the shroud and hub of the impeller are comparable at all simulated conditions. The erosions that occurred in the extended pipes of the inlet and outlet are much smaller than the other components. The relationships between maximum annual erosion rates on different components and the conveying particle concentration are revealed in the Figure 11 at different working flow rates. It is found that the erosion rate is increased significantly with the increasing particle concentration while the rate firstly increases and then decreases. The sharp increase in erosion rate occurs in the range of 0.5~1.5% particle volume fraction. The reason may be attributed to the formation of a buffering layer of particle aggregation on the impinging surfaces [39]. This phenomenon was also observed in the studies of erosion in bend [40]. If the particles accumulate near the eroded surface, the particles in the inner layer are rebounded by the particles in the bulk flow, and their impact velocity and angle are reduced. The surface erosion tends to be saturated, although the particle concentration in the bulk flow still increases. From Figure 12, it is seen that a large number of particles are aggregated on the pressure side of the blade at 2.5% particle volume fraction, and the buffering layer has been forming with the increasing particle concentration. Figure 11. Curves of the erosion rates on the impeller blade and volute with the varying particle concentration under different working flow rates. The relationships between maximum annual erosion rates on different components and the conveying particle concentration are revealed in the Figure 11 at different working flow rates. It is found that the erosion rate is increased significantly with the increasing particle concentration while the rate firstly increases and then decreases. The sharp increase in erosion rate occurs in the range of 0.5~1.5% particle volume fraction. The reason may be attributed to the formation of a buffering layer of particle aggregation on the impinging surfaces [39]. This phenomenon was also observed in the studies of erosion in bend [40]. If the particles accumulate near the eroded surface, the particles in the inner layer are rebounded by the particles in the bulk flow, and their impact velocity and angle are reduced. The surface erosion tends to be saturated, although the particle concentration in the bulk flow still increases. From Figure 12, it is seen that a large number of particles are aggregated on the pressure side of the blade at 2.5% particle volume fraction, and the buffering layer has been forming with the increasing particle concentration.
layer are rebounded by the particles in the bulk flow, and their impact velocity an are reduced. The surface erosion tends to be saturated, although the particle conce in the bulk flow still increases. From Figure 12, it is seen that a large number of are aggregated on the pressure side of the blade at 2.5% particle volume fraction, buffering layer has been forming with the increasing particle concentration.

Influence of Particle Size on the Pump Erosion
The particle size is another important factor in the liquid-solid two-phase flow systems. Since coarse particles are of particular interest in this work, the particle diameters of 0.5 mm, 1 mm, and 2 mm were simulated, and the particle volume fraction was 0.5%. Two important parameters relating to the particle size are the particle Reynolds number and the particle Stokes number. The particle Reynolds numbers, in this case, were approximately 250, 300, and 600, respectively. The Stokes number is defined as the ratio of particle relaxation time and fluid characteristic time is a dimensionless number for representing the curvilinear movement of the solid particle, which is defined as Here, D is the characteristic length, which can be set as the equivalent hydraulic diameter of the flow channel in the impeller. The corresponding Stokes numbers of different particle sizes are calculated as approximately 2.3, 11.3, and 45.3. It is shown that the Stokes number is greater than 1, especially for 2 mm particles and the inertial effect plays a significant role in the particle motion. The larger the particle size is, the weaker the fluid constraint is. Figure 13 shows the particle distribution in a flow passage of the pump impeller for three particle diameters. The larger the particle diameter is, the fewer the particles in the pump. It was seen that the particles are distributed more uniformly in the passage for 2 mm particles, which indicates that the highly inertial particles do not follow the fluid trajectories very closely and do not cluster. As the particle diameter reduces, the Stokes number decreases, particle motion is more affected by the underlying fluid flow, and clustering is observed. More particles are presented in the passage, and the particles intensively move

Influence of Particle Size on the Pump Erosion
The particle size is another important factor in the liquid-solid two-phase flow systems. Since coarse particles are of particular interest in this work, the particle diameters of 0.5 mm, 1 mm, and 2 mm were simulated, and the particle volume fraction was 0.5%. Two important parameters relating to the particle size are the particle Reynolds number and the particle Stokes number. The particle Reynolds numbers, in this case, were approximately 250, 300, and 600, respectively. The Stokes number is defined as the ratio of particle relaxation time and fluid characteristic time is a dimensionless number for representing the curvilinear movement of the solid particle, which is defined as St = ρ p d 2 p u 18µD . Here, D is the characteristic length, which can be set as the equivalent hydraulic diameter of the flow channel in the impeller. The corresponding Stokes numbers of different particle sizes are calculated as approximately 2.3, 11.3, and 45.3. It is shown that the Stokes number is greater than 1, especially for 2 mm particles and the inertial effect plays a significant role in the particle motion. The larger the particle size is, the weaker the fluid constraint is. Figure 13 shows the particle distribution in a flow passage of the pump impeller for three particle diameters. The larger the particle diameter is, the fewer the particles in the pump. It was seen that the particles are distributed more uniformly in the passage for 2 mm particles, which indicates that the highly inertial particles do not follow the fluid trajectories very closely and do not cluster. As the particle diameter reduces, the Stokes number decreases, particle motion is more affected by the underlying fluid flow, and clustering is observed. More particles are presented in the passage, and the particles intensively move along the pressure side of the blade. The results demonstrate that the particle motion in the pump impeller is firstly dominated by the collision between the particles and the pressure side of the blade. If the number of particles increases, the collisions among particles affect the particle distribution. The major part of the particles is located near the pressure side of the blade.
It was seen that the particles are distributed more uniformly in the passage for 2 mm par ticles, which indicates that the highly inertial particles do not follow the fluid trajectories very closely and do not cluster. As the particle diameter reduces, the Stokes number de creases, particle motion is more affected by the underlying fluid flow, and clustering is observed. More particles are presented in the passage, and the particles intensively move along the pressure side of the blade. The results demonstrate that the particle motion in the pump impeller is firstly dominated by the collision between the particles and the pres sure side of the blade. If the number of particles increases, the collisions among particles affect the particle distribution. The major part of the particles is located near the pressure side of the blade. dp = 0.5 mm dp = 1 mm dp = 2 mm The contours of erosion on different components of the pump are figured, and the features on the erosion pattern are very similar to those analyzed in the above section. The severe erosions on the shroud and hub of the impeller occur in the jointed regions with the pressure side of the blade. The severe erosion on the blade occurs in the regions next to the inlet and outlet of the blade. If the particle size is small and more particles are presented, the erosive area is extended toward the middle of the blade. The severe erosion on the volute, as observed in Figure 10, occurs near the pump tongue. Figure 14 shows the maximum annual erosion rates on different components of the pump for different particle sizes. It is seen that the erosions on the shroud and hub of the impeller, the volute, and the extensions of the inlet and outlet are not very sensitive to the particle size, while the erosion on the pressure side of the blade is strongly affected by the particle size. Figure 15 shows the relationships between the maximum annual erosion rate and the particle diameter. On both the blade and volute, the erosion rate is a nonlinear relationship with particle size, and the maximum erosion occurs in the case of the 1 mm particle. For the 2 mm particles, each particle has a larger mass and kinetic energy. However, the presented particles are significantly fewer at the fixed volume fraction, and the collisions with the pump surfaces are scarce, causing the lowest annual erosion rate. For the 1mm particles, although the particle mass and kinetic energy are reduced, the number of particles is significantly increased, as well as the collision frequency with the pump surface. The overall effect is to cause the largest erosion rate. For the 0.5 mm particles, the particle mass and kinetic energy are further reduced while the number of particles is greatest. However, it is found that a buffering layer is formed on the surface when the presented particles are intensively accumulated in the vicinity of the surface. The inner particles are rebounded by the outer particles and move along the surface. The impact velocity and angle are significantly reduced, resulting in the mitigation of surface erosion. The results demonstrate that erosion is determined by both the particle size and concentration near the eroded surface. Moreover, if the particles near the surface accumulate, the mitigation of surface erosion may occur due to the buffering effect mentioned above.
presented particles are intensively accumulated in the vicinity of the surface. T particles are rebounded by the outer particles and move along the surface. The velocity and angle are significantly reduced, resulting in the mitigation of surface The results demonstrate that erosion is determined by both the particle size and tration near the eroded surface. Moreover, if the particles near the surface accumu mitigation of surface erosion may occur due to the buffering effect mentioned ab

Effect of Particle Shape on the Pump Erosion
The effect of particle shape on pump erosion was studied since particles a non-spherical in real working conditions. Three shapes, as shown in Figure 2, we eled, and the resulting erosion on the hydraulic parts of the pump was predicted section, particles of different shapes share the same particle volume so that the energy of the particles is constant at the same velocity. The simulated particle fraction is 0.5%. The specific model of the drag force expressed by Equation (13) for the non-spherical particles in the simulation.
The particle distributions in one flow passage of the pump impeller are pres Figure 16. It is seen that the particles of Ψ = 0.67 tend to move along the pressure the blade, and fewer particles are distributed in the bulk flow of the passage. How the particle sphericity increases, the particles become more diffused in the pass the spherical particles, more particles are distributed in the bulk flow. It was fou the coefficient of the drag force is increased with the reducing particle spheric more irregular particles bear a larger drag force and stronger interaction with th Therefore, the particles follow better with the flow streamlines. Figure 15. Curves of erosion rates on the impeller blade and volute with varying particle sizes.

Effect of Particle Shape on the Pump Erosion
The effect of particle shape on pump erosion was studied since particles are often nonspherical in real working conditions. Three shapes, as shown in Figure 2, were modeled, and the resulting erosion on the hydraulic parts of the pump was predicted. In this section, particles of different shapes share the same particle volume so that the kinetic energy of the particles is constant at the same velocity. The simulated particle volume fraction is 0.5%. The specific model of the drag force expressed by Equation (13) is used for the non-spherical particles in the simulation.
The particle distributions in one flow passage of the pump impeller are presented in Figure 16. It is seen that the particles of Ψ = 0.67 tend to move along the pressure side of the blade, and fewer particles are distributed in the bulk flow of the passage. However, as the particle sphericity increases, the particles become more diffused in the passage. For the spherical particles, more particles are distributed in the bulk flow. It was found that the coefficient of the drag force is increased with the reducing particle sphericity. The more irregular particles bear a larger drag force and stronger interaction with the fluid. Therefore, the particles follow better with the flow streamlines. the blade, and fewer particles are distributed in the bulk flow of the passage. However, as the particle sphericity increases, the particles become more diffused in the passage. For the spherical particles, more particles are distributed in the bulk flow. It was found that the coefficient of the drag force is increased with the reducing particle sphericity. The more irregular particles bear a larger drag force and stronger interaction with the fluid. Therefore, the particles follow better with the flow streamlines.  Figure 17 shows the maximum annual erosion rates on different components of the pump suffering from particles of different shapes. It is seen that the erosion caused by non-spherical particles is larger than that by spherical ones, especially for the pressure side of the impeller blade. The erosion rate on the blade is 3~4 times larger than those on the other components. It is difficult to say the particles with the sharpest shape cause the most severe erosion since the erosion rates caused by the particles of Ψ = 0.84 on the impeller shroud and blade are the largest.  Figure 17 shows the maximum annual erosion rates on different components of the pump suffering from particles of different shapes. It is seen that the erosion caused by non-spherical particles is larger than that by spherical ones, especially for the pressure side of the impeller blade. The erosion rate on the blade is 3~4 times larger than those on the other components. It is difficult to say the particles with the sharpest shape cause the most severe erosion since the erosion rates caused by the particles of Ψ = 0.84 on the impeller shroud and blade are the largest. In order to estimate the effect of particle shape on the erosion of the pressure side o the blade, the erosion rate, collision frequency of the particle-wall collision, and the impac angle are statistically obtained by averaging these parameters along three streamwis curves on the pressure side of the blade. As shown in Figure 18, it is seen that the collisio frequency and impact angle of the non-spherical particles are significantly larger tha those of the spherical ones near the inlet and outlet of the blade. Since the erosion is in creased with the increasing impact angle when the angle is below 50°, the erosion rate i larger, caused by the irregular particles. However, the differences between the two non spherical particles are not distinct, and the difference in the resulting erosion rate is les significant. In order to estimate the effect of particle shape on the erosion of the pressure side of the blade, the erosion rate, collision frequency of the particle-wall collision, and the impact angle are statistically obtained by averaging these parameters along three streamwise curves on the pressure side of the blade. As shown in Figure 18, it is seen that the collision frequency and impact angle of the non-spherical particles are significantly larger than those of the spherical ones near the inlet and outlet of the blade. Since the erosion is increased with the increasing impact angle when the angle is below 50 • , the erosion rate is larger, caused by the irregular particles. However, the differences between the two non-spherical particles are not distinct, and the difference in the resulting erosion rate is less significant. frequency and impact angle of the non-spherical particles are significantly larger th those of the spherical ones near the inlet and outlet of the blade. Since the erosion is creased with the increasing impact angle when the angle is below 50°, the erosion rate larger, caused by the irregular particles. However, the differences between the two no spherical particles are not distinct, and the difference in the resulting erosion rate is l significant.

Conclusions
In this work, a CFD-DEM numerical model was applied in a centrifugal pump to simulate erosion on different pump components. Different models of fluid-solid interphase forces were considered in the numerical model. The particle-turbulence interaction was also implemented in a home-made interface code and combined in the simulation. Different particle shapes were modeled in the particle modeling. The authenticity of the numerical model was systematically validated by comparing the liquid-solid two-phase field with the PIV measurement in the same pump. The predicted erosion rates in a 90 • bend and the hydraulic part of the pump were validated with the measured erosion rate in a 90 • bend under similar particle and hydrodynamic conditions. The validated numerical model was then employed to simulate erosion in the studied pump. The effects of particle volume fraction, particle size, and particle shape on pump erosion were estimated by simulations. It was found that severe erosion occurs near the inlet and outlet of the pressure side of the impeller blade, the middle region of the blade, the joint corners of the shroud and hub of the impeller with the blade, and the volute near the pump tongue. Among these locations, the maximum erosion rate always occurs near the inlet of the pressure side of the impeller blade. The increase in the particle volume fraction brings in more particles in the impeller and pushes the particles along the pressure side of the blade. Mitigation of erosion rate is found in the case of a large particle volume fraction since more particles are accumulated in the near-wall region of the eroded surface and form a buffering layer to reduce the erosion. The particle size barely influences the particle distribution in the impeller. However, it is found that the relationship between the particle size and erosion is nonlinear, and the 1 mm particle causes the largest erosion rate compared to the 0.5 and 2 mm particles at the same volume fraction. Finally, the particle shape can influence the erosion rate. The sharper particles receive a larger drag force and tend to move along the pressure side of the blade. The resulting erosion rates on different components are larger than those caused by spherical particles. The collision frequency and impact angle of the non-spherical particles are also larger than those of the spherical ones. The difference between the two irregular particles is small, and the difference of the resulting erosion is less significant.