# Modeling and Numerical Analysis of the Solid Particle Erosion in Curved Ducts

Authors: Ke Sun; Lin Lu; Hanhui Jin
Journal: Abstract and Applied Analysis
 (2013)

---

## Abstract

This paper presents a modeling and computational study on particle erosion in curved ducts. It is found that the average erosion rates per impact range from4.2×10−7 to 9.5×10−3 mm3/g under current conditions. For each doubled inlet velocity, the increases of erosion rates per impact are 2–14 times. The erosion rate per impact varies with particle diameter with “√” shape through bends, which is similar to the particle deposition behavior in duct flows. The erosion rate curves per injected particle show the shapes of a 90-degree anticlockwise rotated “S” and a wide open “V,” respectively, for three larger and smaller inlet velocities. The average erosion rates per injected particle are 1.4–18.9 times those rates per impact due to huge amounts of impacting, especially for those depositing particles. It is obvious that the erosion rate distribution per impact is similar to a “fingerprint” with five clear stripes and a lower “cloud” along the bend deflection angle for the three largest particles; yet, for other smaller particles, the erosion rate distributions are much like an entire “cloud.”


---

## Body

## 1. Introduction

Particulate flow is a significant phenomenon in environmental, industrial, medical, and lifetime applications. For example, the conveying and ventilation system is severely affected by particle flow and its erosion [1–3]. The aerosol deposition, accumulation, and soiling on solar panels and glazing glass will erode their surfaces obviously [4]. These applications also include a large amount of straight and curved duct flow. However, the particle flow and erosion in bends have not been fully studied [5, 6].El-Behery et al. [7] studied the penetration rate distribution for estimating solid particle erosion in curved 90° and 180° ducts. Sun et al. [8] studied the particle penetration and deposition in and behind bends. Chen et al. [9] proposed a comprehensive procedure to estimate the erosion in elbows mainly for bubbly or droplet flow. Zhang et al. [10] adopted Computational Fluid Dynamics (CFD) with near-wall and volume improvement to predict particle impact in a sharp bend. Macchini et al. [11] investigated the influence of particle size, density, and concentration on bend erosive wear under high particle concentrations. Although some studies have been conducted to investigate the particle erosion in bends, prediction research on the erosion distribution is limited.Therefore, present study focuses on the erosion status modeling and analysis of curved walls. The airflow and particle flow are modeled and predicted to obtain the erosion information in typical 90-degree bends. The behaviors of erosion rates varying with particle diameter or Stokes number, inlet velocity, and deflection angles are analyzed in detail.


## 2. Method

### 2.1. Fluid and Particle Flow Model

Before solving the particle flow in a fluid, the fluid flow conservation equations were deducted and predicted by CFD tools [12]. Together with near-wall two-layer model, the Reynolds stress model (RSM) was utilized to solve the fluid flow system in the curved duct due to its ability to predict the swirling, turbulent, and near-surface flow [8]. Details of the RSM model are given as follows:
(1)∂u-i∂xi=0,∂u-i∂t+u-j∂u-i∂xj=-1ρ∂p-∂xi+1ρ∂∂xj(μ∂u-i∂xj-ρui′uj′¯),
where u-i is the mean velocity, r and m are the density and molecular viscosity of the fluid, respectively, p- is the mean pressure, and -ρui′uj′¯ is the Reynolds stress. The general transport equations for the Reynolds stresses could be estimated by
(2)∂∂t(ρui′uj′¯)+∂∂xk(ρuk¯ui′uj′¯)=-∂∂xk[ρui′uj′uk′¯+p(δkjui′+δikuj′)¯]+∂∂xk[μ∂∂xk(ui′uj′¯)]-ρ(ui′uk′¯∂uj¯∂xk+uj′uk′¯∂ui¯∂xk)-ρβ(giuj′θ¯+gjui′θ¯)+p(∂ui′∂xj+∂uj′∂xi)¯-2μ∂ui′∂xk∂uj′∂xk¯-2ρΩk(uj′um′¯Eikm+ui′um′¯Ejkm)+Suser,
where the 1st and 4th–6th terms on the right of this equation have detailed models to close this equation as that in [12]. Neumann outlet flow conditions were employed and the walls were set as smooth. Staggered numerical scheme was adopted for pressure equation and second-order discretization methods were applied for other variables. These variable convergences were achieved when their residuals were 10-5 or less.Based on the modeling and computation of fluid flow field, the particle flow was determined by tracking the path of each particle with the Lagrangian method. When using this method, particles were assumed to be spherical solid ones with diluted flow. The Lagrangian equation can be expressed as(3)dupidt=FD(ui-upi)+gi(1-ρaρp)+Fai,
where upi is particle velocity in the ith direction (m/s) and ui represents the air velocity. The 1st term on the formula right represents the particle drag force in the ith direction (m/s2), and the 2nd one is the gravitational force. The 3rd term stands for other possible forces [13]. To simplify the modelling of particle forces, this article only adopts lift force to predict the particle flow through curved duct [14]. The random fluctuation effects of the fluid flow on particles were modeled by “Eddy lifetime” method [15].


### 2.2. Particle Wall Interaction and Erosion Models

When particles flow along constrained domains, they will interact with different surfaces. For present solid particles, they would deposit on, rebound from, or reimpact with wall surfaces. To model this process, an algebraic particle-wall impact model [16] was utilized to describe the interaction process between particle and wall as follows:
(4)e=-v2,npv1,np=Rpw(1-ρpw),v2,np=-ev1,np,v2,tp=v1,tp-μ1(1+e)v1,np,
where e is the restitution coefficient, v1,np and v2,np are particle normal incident and reflected velocities, respectively, Rpw stands for the coefficient of restitution without adhesion, ρpw is the coefficient of adhesion, and μ1 is the ratio of tangential to normal impulse. This model was integrated into the main particle computation by in-house codes.When particle-wall impact happens, the wall will be eroded by impacted particles gradually. The wall surface erosion rate or wastage rate is commonly defined as the ratio of wall eroded mass or volume to that of impacted particles. The erosion rate is usually determined by particle impacting velocity and angle, wall and particle materials, and their temperatures. In this paper, a convenient model developed by Menguturk et al. [17, 18] was introduced into the curved duct flow to analyse the wall erosion condition on the curved duct walls as follows:
(5)Ev=1.63×10-6(v1pcosα)2.5sin(πα2αc)+4.68×10-7(v1psinα)2.5(α≤αc),Ev=1.63×10-6(v1pcosα)2.5+4.68×10-7(v1psinα)2.5(α>αc),
where the unit of erosion rate Ev is mm3/g, v1p is the incident impact velocity, α is the incident impact angle near the wall, and αc is the critical incident angle with a constant value of 22.7°. In these formulas, the erosion rate Ev has a power relationship with index 2.5 with impact velocity v1p. This model was incorporated into the particle-wall interaction model and main particle flow model to obtain the impact and erosion information in the curved duct flow.


## 2.1. Fluid and Particle Flow Model

Before solving the particle flow in a fluid, the fluid flow conservation equations were deducted and predicted by CFD tools [12]. Together with near-wall two-layer model, the Reynolds stress model (RSM) was utilized to solve the fluid flow system in the curved duct due to its ability to predict the swirling, turbulent, and near-surface flow [8]. Details of the RSM model are given as follows:
(1)∂u-i∂xi=0,∂u-i∂t+u-j∂u-i∂xj=-1ρ∂p-∂xi+1ρ∂∂xj(μ∂u-i∂xj-ρui′uj′¯),
where u-i is the mean velocity, r and m are the density and molecular viscosity of the fluid, respectively, p- is the mean pressure, and -ρui′uj′¯ is the Reynolds stress. The general transport equations for the Reynolds stresses could be estimated by
(2)∂∂t(ρui′uj′¯)+∂∂xk(ρuk¯ui′uj′¯)=-∂∂xk[ρui′uj′uk′¯+p(δkjui′+δikuj′)¯]+∂∂xk[μ∂∂xk(ui′uj′¯)]-ρ(ui′uk′¯∂uj¯∂xk+uj′uk′¯∂ui¯∂xk)-ρβ(giuj′θ¯+gjui′θ¯)+p(∂ui′∂xj+∂uj′∂xi)¯-2μ∂ui′∂xk∂uj′∂xk¯-2ρΩk(uj′um′¯Eikm+ui′um′¯Ejkm)+Suser,
where the 1st and 4th–6th terms on the right of this equation have detailed models to close this equation as that in [12]. Neumann outlet flow conditions were employed and the walls were set as smooth. Staggered numerical scheme was adopted for pressure equation and second-order discretization methods were applied for other variables. These variable convergences were achieved when their residuals were 10-5 or less.Based on the modeling and computation of fluid flow field, the particle flow was determined by tracking the path of each particle with the Lagrangian method. When using this method, particles were assumed to be spherical solid ones with diluted flow. The Lagrangian equation can be expressed as(3)dupidt=FD(ui-upi)+gi(1-ρaρp)+Fai,
where upi is particle velocity in the ith direction (m/s) and ui represents the air velocity. The 1st term on the formula right represents the particle drag force in the ith direction (m/s2), and the 2nd one is the gravitational force. The 3rd term stands for other possible forces [13]. To simplify the modelling of particle forces, this article only adopts lift force to predict the particle flow through curved duct [14]. The random fluctuation effects of the fluid flow on particles were modeled by “Eddy lifetime” method [15].


## 2.2. Particle Wall Interaction and Erosion Models

When particles flow along constrained domains, they will interact with different surfaces. For present solid particles, they would deposit on, rebound from, or reimpact with wall surfaces. To model this process, an algebraic particle-wall impact model [16] was utilized to describe the interaction process between particle and wall as follows:
(4)e=-v2,npv1,np=Rpw(1-ρpw),v2,np=-ev1,np,v2,tp=v1,tp-μ1(1+e)v1,np,
where e is the restitution coefficient, v1,np and v2,np are particle normal incident and reflected velocities, respectively, Rpw stands for the coefficient of restitution without adhesion, ρpw is the coefficient of adhesion, and μ1 is the ratio of tangential to normal impulse. This model was integrated into the main particle computation by in-house codes.When particle-wall impact happens, the wall will be eroded by impacted particles gradually. The wall surface erosion rate or wastage rate is commonly defined as the ratio of wall eroded mass or volume to that of impacted particles. The erosion rate is usually determined by particle impacting velocity and angle, wall and particle materials, and their temperatures. In this paper, a convenient model developed by Menguturk et al. [17, 18] was introduced into the curved duct flow to analyse the wall erosion condition on the curved duct walls as follows:
(5)Ev=1.63×10-6(v1pcosα)2.5sin(πα2αc)+4.68×10-7(v1psinα)2.5(α≤αc),Ev=1.63×10-6(v1pcosα)2.5+4.68×10-7(v1psinα)2.5(α>αc),
where the unit of erosion rate Ev is mm3/g, v1p is the incident impact velocity, α is the incident impact angle near the wall, and αc is the critical incident angle with a constant value of 22.7°. In these formulas, the erosion rate Ev has a power relationship with index 2.5 with impact velocity v1p. This model was incorporated into the particle-wall interaction model and main particle flow model to obtain the impact and erosion information in the curved duct flow.


## 3. Results and Discussion

In this paper, the background fluid was selected as air, and the flow domain was formed in the typical curved two-dimensional duct with a 90-degree bend as shown in Figure1. Since particles mainly deposit and impact on the outer bend wall, the erosions are also observed to happen at these places. The deflection angle θ is also demonstrated from 0° to 90° in this figure. Eight particle groups were injected into the duct flow from the duct inlet with even space distribution, which have the diameters dp=1,3,7,16,30,60,100, and 200 μm. At the duct inlet, 60000 particles were injected. The volume fraction is less than 1.02×10-6. Six inlet velocity uin conditions were adopted in this analysis, that is, 2.2,5.3,9,15,33, and 52 m/s. Totally, more than 48 different conditions were analyzed. The airflow and particle flow models were validated by experimental data through previous works [8, 19] in terms of the air and particle velocity profiles, and particle depositions and penetrations, which support the present investigations.Figure 1
Schematic chart of the model geometry.


### 3.1. Erosion Rate per Impact on Curved Walls

Based on the airflow and particle flow modeling and prediction, statistic results of particle-wall impact and erosion status are given in Figure1. This figure shows the average erosion rate per impacted particle on the curved wall for each inlet velocity and particle diameter conditions. Using the concept of erosion rate per impacted particle can have more physical meanings and higher application abilities.Generally, the average erosion rates shown in Figure2 range from 9.5×10-3 to 4.2×10-7 mm3/g for the inlet velocity from 52 m/s to 2.2 m/s. Based on these erosion rate results, the wall wastage can be calculated out considering the injected particle mass. The erosion rate changes with inlet velocity sharply as demonstrated in the figure. The inlet velocity conditions are roughly set around 2-time increase from one smaller inlet velocity to another bigger one. However, the erosion rate generally increases nearly 2 to 14 times. This phenomenon can be explained from the erosion rate equation in (5). Although the erosion rate has a power relationship with impact velocity with index 2.5, the erosion rate has much more complicated nonlinear correlation with inlet velocity due to the complex curved duct flow conditions.Comparison of the erosion rateEv per impact particle against the particle diameter and inlet flow velocity: (a) larger inlet velocities of 15,33, and 52 m/s and (b) smaller inlet velocities of 2.2, 5.3, and 9 m/s.
(a)
(b)In terms of the changing with particle diameter, the general trend of erosion rate is similar for particle diameter fromdp=1 to 200 μm for different inlet velocities. Firstly, the erosion rate decreases with particle diameter increase. When particle diameter dp is 7 or 16 μm, the erosion rate gets to the lowest value. After that, it increases with the particle diameter increase. For very large flow inlet velocity and particle diameter, like dp=200μm at uin=33 and 52 m/s, the erosion rate decreases a little due to the large inertia and gravity of these large particles with large Stokes number St. These phenomena are similar to that observed in the particle deposition behavior with “√” shape in duct flows [20, 21]. This phenomenon may be attributed to the similar impact mechanism during the deposition, rebounding, and reimpaction process. Smaller particles of dp=1 to 3 μm are affected more by diffusion, turbulence, and vortex in the bend and they may have higher impact angle, and thus higher erosion rates are observed. Coarse ones of dp=30 to 200 μm are controlled by inertia and gravity, and thus they would erode more wall masses.


### 3.2. Erosion Rate per Injected Particle on Curved Walls

Since particles may rebound and reimpact onto walls due to the near-wall forces or mechanism like turbulence and vortex, average erosion rate per impacted particle cannot give macroinformation about the injected particle erosion ability. Therefore, Figure3 shows the average erosion rate per injected particles statistically against the particle Stokes number. The Stokes number St is a key parameter to describe the particle response to fluid flow in bends. This dimensionless number is determined by the following formula:
(6)St=τpτb=τpuin(Dh/2)=Ccρpdp218μ·uin(Dh/2),
where τp is the particle relaxation time, τb is the time scale of the bend, uin is the air velocity at the duct inlet, Dh is the hydraulic diameter of the duct, Cc is the Cunningham slip correction factor for microparticles, ρp and dp are particle density and diameter, respectively, and μ is the air dynamic viscosity. The friction velocity, u*, shown in Figure 3 could be calculated by
(7)u*=uinf2=uin/2-3.6log[6.9/Re+(k′/3.7Dh)1.11]=-uin5.1log[6.9/Re+(k′/3.7Dh)1.11],
where f is the Fanning friction factor determined in straight duct by White [22], k′ is the average microscale roughness height of the rough wall, which is zero in this work for smooth walls, and Re is the Reynolds number.Comparison of the erosion rateEv per injected particle against the particle diameter and inlet flow velocity: (a) larger inlet velocities of 15, 33, and 52 m/s and (b) smaller inlet velocities of 2.2, 5.3, and 9 m/s.
(a)
(b)The erosion rate curves shown in Figure3 are generally smoother than those in Figure 2. The curve shapes in Figure 3(a) are much like a 90-degree anticlockwise rotated “S.” The erosion rates range from 9.2×10-5 to 2.0×10-2 mm3/g. The magnitudes of the values are generally much larger than those in Figure 2(a) due to huge amounts of impacting, especially for those depositing particles. For a similar condition, the erosion rates in Figure 3(a) are 1.4–18.9 times those in Figure 2(a). If the solid particles have high velocities and will deposit onto the bend surfaces, they may impact, rebound, and reimpact on the bend walls to gradually lower down its velocity for deposition. Furthermore, since particles will deposit more on ducts for larger inlet velocities or friction velocities as demonstrated in previous research [19, 20], the erosion rates in Figure 3(a) are generally high. In addition, smaller particles of Stokes number St=3.3×10-3–7.2×10-1 have higher impaction frequencies or amounts than the other larger ones as shown in the figure. This behavior may be due to their frequent interaction with airflows.For lower inlet velocity cases shown in Figure3(b), the erosion rate curves are much like wide open “V” shapes. If larger particles of dp>200μm are simulated, the curves may become the shapes demonstrated in Figure 3(a). The erosion rates range from 1.9×10-7 to 3.7×10-4 mm3/g. Most of these erosion rates are still larger than those shown in Figure 2(b) due to multi-impaction of a single depositing particle. The increasing degrees in this figure are smaller than those in Figure 3(a). For friction velocity u*=0.13 m/s, the increase degree is the smallest. Furthermore, three erosion rates are smaller than those in Figure 2(b) for the particle groups of St=0.02,0.11,0.37. This phenomenon may be attributed to the low particle impaction amount and low particle deposition for these particles with the smallest friction velocity.


### 3.3. Erosion Rate Distribution for Each Impact

Distribution of the erosion rate per impacted particle can further demonstrate the particle erosion picture on bend walls. Detailed erosion magnitude and locations can be interpreted. Particle depositions in bends mainly happen on the outer bend walls as shown in Figure1 and previous research results [23]. Therefore, the following analysis primarily states the erosions on different locations of the outer bend wall.Figures4–7 show the distribution of erosion rate Ev per impacted particle against the deflection angle as demonstrated in Figure 1. The angles from 0° to 90° represent bend inlet to outlet. The case of the inlet velocity uin=5.3 m/s is selected as an example to explain the erosion rate distribution. Figure 4 presents the erosion rate for the largest particles with diameters dp=100 and 200 μm. It is very obvious that the erosion rate is like a “fingerprint” with clear stripes. These stripes are probably caused by the clear impactions of particles. The top stripe is the first impact of particles on the outer bend wall with the largest erosion rate. Following the airflow, the deflection angles with first impaction are from 0° to 60°, where the peak erosion rate appears at around 20°. After this impaction, four much clear impactions are observed but with lower erosion rate due to moment or velocity reduction after each impaction. The lowest unclear cloud-like areas are the lightest impactions close to particle depositions. They are more randomly distributed due to turbulence, impaction, and diffusion. Deeper color represents that the erosion rates have a large amount at these areas. It can be seen that the particles in this figure have very high erosion rates from Ev=1.0×10-5 to 1.0×10-4 mm3/g. Furthermore, the “fingerprint” phenomenon in this figure indicates that further impactions would happen after the bend because the stripes do not have end tails at 90° deflection angle or bend outlet.Erosion rateEv distribution per impact particle against the deflection angle of curved duct at the inlet velocity 5.3 m/s: (a) particle diameter 200 μm and (b) particle diameter 100 μm.
(a)
(b)Erosion rateEv distribution per impact particle against the deflection angle of curved duct at the inlet velocity 5.3 m/s: (a) particle diameter 60 μm and (b) particle diameter 30 μm.
(a)
(b)Erosion rateEv distribution per impact particle against the deflection angle of curved duct at the inlet velocity 5.3 m/s: (a) particle diameter 16 μm and (b) particle diameter 7 μm.
(a)
(b)Erosion rateEv distribution per impact particle against the deflection angle of curved duct at the inlet velocity 5.3 m/s: (a) particle diameter 3 μm and (b) particle diameter 1 μm.
(a)
(b)Figure5 demonstrates the erosion rates for particles with diameters dp=60 and 30 μm. From Figure 4 to Figure 5(a), the stripe boundaries are gradually obscure, and they begin to widen their widths. From Figure 5(a) to Figure 5(b), stripes have disappeared, and the general erosion rates decrease. In Figure 5(b), most of the erosions are from Ev=1.0×10-8 to 1.0×10-6 mm3/g. From Figures 5 to 6, erosions decrease with impaction and deposition decreases because particle diameter decreases. From Figures 6 to 7, erosion amounts increase again as seen in Figures 2 and 3 due to diffusion-induced impactions of smaller particles. In Figure 7(b), most of the erosions are from Ev=5.0×10-7 to 1.0×10-5 mm3/g. These larger erosion values and amounts may be caused by larger incident angle impactions due to diffusion and turbulent eddy mixing.


## 3.1. Erosion Rate per Impact on Curved Walls

Based on the airflow and particle flow modeling and prediction, statistic results of particle-wall impact and erosion status are given in Figure1. This figure shows the average erosion rate per impacted particle on the curved wall for each inlet velocity and particle diameter conditions. Using the concept of erosion rate per impacted particle can have more physical meanings and higher application abilities.Generally, the average erosion rates shown in Figure2 range from 9.5×10-3 to 4.2×10-7 mm3/g for the inlet velocity from 52 m/s to 2.2 m/s. Based on these erosion rate results, the wall wastage can be calculated out considering the injected particle mass. The erosion rate changes with inlet velocity sharply as demonstrated in the figure. The inlet velocity conditions are roughly set around 2-time increase from one smaller inlet velocity to another bigger one. However, the erosion rate generally increases nearly 2 to 14 times. This phenomenon can be explained from the erosion rate equation in (5). Although the erosion rate has a power relationship with impact velocity with index 2.5, the erosion rate has much more complicated nonlinear correlation with inlet velocity due to the complex curved duct flow conditions.Comparison of the erosion rateEv per impact particle against the particle diameter and inlet flow velocity: (a) larger inlet velocities of 15,33, and 52 m/s and (b) smaller inlet velocities of 2.2, 5.3, and 9 m/s.
(a)
(b)In terms of the changing with particle diameter, the general trend of erosion rate is similar for particle diameter fromdp=1 to 200 μm for different inlet velocities. Firstly, the erosion rate decreases with particle diameter increase. When particle diameter dp is 7 or 16 μm, the erosion rate gets to the lowest value. After that, it increases with the particle diameter increase. For very large flow inlet velocity and particle diameter, like dp=200μm at uin=33 and 52 m/s, the erosion rate decreases a little due to the large inertia and gravity of these large particles with large Stokes number St. These phenomena are similar to that observed in the particle deposition behavior with “√” shape in duct flows [20, 21]. This phenomenon may be attributed to the similar impact mechanism during the deposition, rebounding, and reimpaction process. Smaller particles of dp=1 to 3 μm are affected more by diffusion, turbulence, and vortex in the bend and they may have higher impact angle, and thus higher erosion rates are observed. Coarse ones of dp=30 to 200 μm are controlled by inertia and gravity, and thus they would erode more wall masses.


## 3.2. Erosion Rate per Injected Particle on Curved Walls

Since particles may rebound and reimpact onto walls due to the near-wall forces or mechanism like turbulence and vortex, average erosion rate per impacted particle cannot give macroinformation about the injected particle erosion ability. Therefore, Figure3 shows the average erosion rate per injected particles statistically against the particle Stokes number. The Stokes number St is a key parameter to describe the particle response to fluid flow in bends. This dimensionless number is determined by the following formula:
(6)St=τpτb=τpuin(Dh/2)=Ccρpdp218μ·uin(Dh/2),
where τp is the particle relaxation time, τb is the time scale of the bend, uin is the air velocity at the duct inlet, Dh is the hydraulic diameter of the duct, Cc is the Cunningham slip correction factor for microparticles, ρp and dp are particle density and diameter, respectively, and μ is the air dynamic viscosity. The friction velocity, u*, shown in Figure 3 could be calculated by
(7)u*=uinf2=uin/2-3.6log[6.9/Re+(k′/3.7Dh)1.11]=-uin5.1log[6.9/Re+(k′/3.7Dh)1.11],
where f is the Fanning friction factor determined in straight duct by White [22], k′ is the average microscale roughness height of the rough wall, which is zero in this work for smooth walls, and Re is the Reynolds number.Comparison of the erosion rateEv per injected particle against the particle diameter and inlet flow velocity: (a) larger inlet velocities of 15, 33, and 52 m/s and (b) smaller inlet velocities of 2.2, 5.3, and 9 m/s.
(a)
(b)The erosion rate curves shown in Figure3 are generally smoother than those in Figure 2. The curve shapes in Figure 3(a) are much like a 90-degree anticlockwise rotated “S.” The erosion rates range from 9.2×10-5 to 2.0×10-2 mm3/g. The magnitudes of the values are generally much larger than those in Figure 2(a) due to huge amounts of impacting, especially for those depositing particles. For a similar condition, the erosion rates in Figure 3(a) are 1.4–18.9 times those in Figure 2(a). If the solid particles have high velocities and will deposit onto the bend surfaces, they may impact, rebound, and reimpact on the bend walls to gradually lower down its velocity for deposition. Furthermore, since particles will deposit more on ducts for larger inlet velocities or friction velocities as demonstrated in previous research [19, 20], the erosion rates in Figure 3(a) are generally high. In addition, smaller particles of Stokes number St=3.3×10-3–7.2×10-1 have higher impaction frequencies or amounts than the other larger ones as shown in the figure. This behavior may be due to their frequent interaction with airflows.For lower inlet velocity cases shown in Figure3(b), the erosion rate curves are much like wide open “V” shapes. If larger particles of dp>200μm are simulated, the curves may become the shapes demonstrated in Figure 3(a). The erosion rates range from 1.9×10-7 to 3.7×10-4 mm3/g. Most of these erosion rates are still larger than those shown in Figure 2(b) due to multi-impaction of a single depositing particle. The increasing degrees in this figure are smaller than those in Figure 3(a). For friction velocity u*=0.13 m/s, the increase degree is the smallest. Furthermore, three erosion rates are smaller than those in Figure 2(b) for the particle groups of St=0.02,0.11,0.37. This phenomenon may be attributed to the low particle impaction amount and low particle deposition for these particles with the smallest friction velocity.


## 3.3. Erosion Rate Distribution for Each Impact

Distribution of the erosion rate per impacted particle can further demonstrate the particle erosion picture on bend walls. Detailed erosion magnitude and locations can be interpreted. Particle depositions in bends mainly happen on the outer bend walls as shown in Figure1 and previous research results [23]. Therefore, the following analysis primarily states the erosions on different locations of the outer bend wall.Figures4–7 show the distribution of erosion rate Ev per impacted particle against the deflection angle as demonstrated in Figure 1. The angles from 0° to 90° represent bend inlet to outlet. The case of the inlet velocity uin=5.3 m/s is selected as an example to explain the erosion rate distribution. Figure 4 presents the erosion rate for the largest particles with diameters dp=100 and 200 μm. It is very obvious that the erosion rate is like a “fingerprint” with clear stripes. These stripes are probably caused by the clear impactions of particles. The top stripe is the first impact of particles on the outer bend wall with the largest erosion rate. Following the airflow, the deflection angles with first impaction are from 0° to 60°, where the peak erosion rate appears at around 20°. After this impaction, four much clear impactions are observed but with lower erosion rate due to moment or velocity reduction after each impaction. The lowest unclear cloud-like areas are the lightest impactions close to particle depositions. They are more randomly distributed due to turbulence, impaction, and diffusion. Deeper color represents that the erosion rates have a large amount at these areas. It can be seen that the particles in this figure have very high erosion rates from Ev=1.0×10-5 to 1.0×10-4 mm3/g. Furthermore, the “fingerprint” phenomenon in this figure indicates that further impactions would happen after the bend because the stripes do not have end tails at 90° deflection angle or bend outlet.Erosion rateEv distribution per impact particle against the deflection angle of curved duct at the inlet velocity 5.3 m/s: (a) particle diameter 200 μm and (b) particle diameter 100 μm.
(a)
(b)Erosion rateEv distribution per impact particle against the deflection angle of curved duct at the inlet velocity 5.3 m/s: (a) particle diameter 60 μm and (b) particle diameter 30 μm.
(a)
(b)Erosion rateEv distribution per impact particle against the deflection angle of curved duct at the inlet velocity 5.3 m/s: (a) particle diameter 16 μm and (b) particle diameter 7 μm.
(a)
(b)Erosion rateEv distribution per impact particle against the deflection angle of curved duct at the inlet velocity 5.3 m/s: (a) particle diameter 3 μm and (b) particle diameter 1 μm.
(a)
(b)Figure5 demonstrates the erosion rates for particles with diameters dp=60 and 30 μm. From Figure 4 to Figure 5(a), the stripe boundaries are gradually obscure, and they begin to widen their widths. From Figure 5(a) to Figure 5(b), stripes have disappeared, and the general erosion rates decrease. In Figure 5(b), most of the erosions are from Ev=1.0×10-8 to 1.0×10-6 mm3/g. From Figures 5 to 6, erosions decrease with impaction and deposition decreases because particle diameter decreases. From Figures 6 to 7, erosion amounts increase again as seen in Figures 2 and 3 due to diffusion-induced impactions of smaller particles. In Figure 7(b), most of the erosions are from Ev=5.0×10-7 to 1.0×10-5 mm3/g. These larger erosion values and amounts may be caused by larger incident angle impactions due to diffusion and turbulent eddy mixing.


## 4. Conclusions

This paper presented a modeling and computational study on particle erosion in curved ducts. An algebraic particle-wall impact model and a convenient erosion model were adopted in a typical 90-degree bend. The erosion rates per impacted particle and per injected particle were investigated against different particle diameter, inlet velocity, and Stokes number. The distribution of erosion rate was visualized and analyzed along bend deflection angles. The major findings are summarized as follows.(1)
Generally, the average erosion rates per impacted particles range from4.2×10-7 to 9.5×10-3 mm3/g for the inlet velocity from 2.2 m/s to 52 m/s. For each doubled inlet velocity, the increases of erosion rates per impact are 2–14 times. This phenomenon indicates a sharp enhancement of erosion rate due to inlet velocity increase.(2)
The erosion rate per impact changes with particle diameter with “√” shape in bend flows, which is similar to the particle deposition behavior in duct flows. When particle diameter dp is 7 or 16 μm, the erosion rate gets to the lowest value. Other smaller ones have higher erosion rates due to diffusion, turbulence, and vortex. Higher erosion rates for larger particles are caused by inertia and gravity.(3)
The erosion rate curves per injected particle are much like a 90-degree anticlockwise rotated “S” for three larger inlet velocities of 15, 33, and 52 m/s. The erosion rates range from9.2×10-5 to 2.0×10-2 mm3/g. The magnitudes of the values are 0.4-to-17.9 times larger than those per impact due to huge amounts of impacting, especially for those depositing particles. For the lower three inlet velocities, the erosion rate curves are much like wide open “V” shapes. Their erosion rates range from 1.9×10-7 to 3.7×10-4 mm3/g.(4)
It is obvious that the erosion rate is like a “fingerprint” with clear stripes for the three largest particles with diametersdp=60, 100, and 200 μm. The deflection angles with first impaction stripe are from 0° to 60°, where the peak erosion rate appears at around 20°. After this impaction, other four clear impaction stripes are observed but with lower value due to the velocity reduction after each impaction. The lowest cloud-like areas are the lightest impactions close to particle depositions.(5)
For other smaller particles, the erosion rate distribution is much like an entire “cloud” along the bend deflection angle. These “clouds” become smaller and then bigger when the particle diameter decreases. The deflection angles of “clouds” are mainly from 20° to 90°.



---

*Source: 245074-2013-11-19.xml*