Abstract

A ductile erosion model embodying the mechanisms of erosion involving cutting wear and repeated plastic deformation has been developed to predict erosion rates of boiler grade steels. The issue of erosion-oxidation interaction has also been addressed to further predict the mass loss resulted from this composite mechanism. A deterministic formalism for the kinetics of oxide-scale growth and a probabilistic approach to characterize the material loss are employed to describe simultaneous actions of high-temperature oxidation and mechanical erosion. The model predictions are in good agreement with the published data.

1. Introduction

Erosion is a process by which material is removed from the layers of a surface impacted by a stream of abrasive particles. The magnitude of the wear is quantified by the volume or mass of the material that is removed by the action of the impacting particles. In coal-fired boilers, ash particles collide with the surface of the boiler steel components (air preheater, economizer, and super-heater tubes) resulting in considerable erosion of metallic materials. Such erosion, together with the processes of blocking, fouling, and corrosion, shortens the service life of the boiler components. The resulting penalty is not only the cost of replacing the components but also the cost of stoppage of power production. It is desirable, therefore, to be able to predict the rate of erosion of the boiler components in order to systematically plan the maintenance and/or replacement schedules of these components to avoid forced outages [1–5].

High-temperature erosion, which involves the conjoint action of high-temperature oxidation and damage due to solid particle impact, is a complex process. Material damage occurs through the interaction of both oxidation and erosion, with metal recession rates often enhanced as a result of the synergetic effect of these two processes. In this paper, mathematical investigations have been carried out on the basis of a ductile erosion model to characterize and quantify the wear behavior and mechanism of various boiler grade steels. The erosion-oxidation synergetic impact is also modeled to characterize the material damage in a high-temperature-oxidizing environment based on a composite deterministic-probabilistic modeling framework.

2. Mathematical Simulation

There are three important phenomena by which metal is removed by the impacting particles at elevated temperature:

(i)removal of material due to cutting wear,(ii)removal of material due to repeated plastic deformation,(iii)effect of temperature on tensile properties of the material. The first two phenomena are applicable to erosion at room temperature.

2.1. Mechanism of Cutting Wear

An ash particle, that strikes the surface at an acute angle and at a velocity greater than the critical velocity needed for the penetration of the material’s surface, removes some material in a process similar to the cutting action of a machine tool (micromachining action). At the impact location, the particle loses a fraction of its kinetic energy to the target material in the form of heat and energy for deformation of the surface. When the shear strain exceeds the elastic strain limit of the target material, the particle penetrates the surface of the material and ploughs along the surface, removing material.

In the case of cutting wear, the material's resistance to deformation needs to be overcome by the impacting particles to penetrate the surface of the material. The differential equation of motion for the depth of penetration is given as [1–3] π‘šπ‘π‘‘2𝐿𝑑𝑑2𝑑=βˆ’πœ‹π‘2πΏπ‘πœŽπ‘¦,(1) where L is the depth of penetration, π‘šπ‘ is the mass of the particle, 𝑑𝑝 is the diameter of the particle, t is the time, πœŽπ‘¦ is the yield stress of the target material, and c is a particle shape factor. For sphere particles, shape factor is equal to 3.

The erosion rate due to cutting wear πœ€π‘, defined as the ratio of the mass of the material eroded (π‘š), from the target surface to mass of the impacting particle (π‘šπ‘), is given by the following equation [1–3]: πœ€π‘=π‘šπ‘šπ‘=π‘˜π‘πœŒπ‘šπœŒπ‘3/2𝑑3𝑝𝑣3sin3𝛽33πœŽπ‘¦3/2ξ‚€πœ‹πœŒπ‘ξ‚€π‘‘3𝑝=π‘˜/61πœŒπ‘šπœŒπ‘1/2𝑣3sin3π›½πœŽπ‘¦3/2,(2) where 𝑉 and 𝛽 are the impacting velocity and angle of the particle, respectively. πœŒπ‘š and πœŒπ‘ are the densities of the target material, and the particle 𝐾𝑐 and 𝐾1 are constants [1–3].

2.2. Mechanism of Plastic Deformation Wear

During particle impact, the loss of material from an eroding surface may occur as a consequence of a composite sequential mechanism similar to a combined process of extrusion-forging mechanism. The extruded platelets from shallow craters made by the impacting particle are further forged into a strained condition. Because of the high-strain rates, adiabatic shear heating occurs in the surface region immediate to the impact site. The kinetic energy of the impacting particles is sufficient to generate considerably greater force being imparted to the metal than is required to create platelets on the surface and as a consequence, a work hardened zone formed, and steady state erosion begins.

Particle kinetic energy, responsible for extrusion-forging process during plastic deformation, is given as [1–3] 𝐸1=12πœ‹π‘‘3𝑝6πœŒπ‘π‘‰2sin2πœ‹π›½=𝜌12𝑝𝑑3𝑝𝑉2sin2𝛽,(3) where 𝑑𝑝 and πœŒπ‘ are the particle diameter and density, respectively, and 𝑉 and 𝛽 are the particle incident velocity and angle, respectively. The final expression for the erosion rate, due to plastic deformation, πœ€π‘, is [1–3] πœ€π‘=π‘šπ‘šπ‘=πΎπ‘πœŒπ‘šπœŒπ‘3/2𝑑3𝑝𝑉3sin3𝛽𝐻𝑉3/2ξ‚€πœŒπ‘πœ‹π‘‘3𝑝=𝐾/62πœŒπ‘šπœŒπ‘1/2𝑉3sin3𝛽𝐻𝑉3/2,(4) where 𝐾2 is a constant [1–3].

2.3. Overall Erosion Rate

It is difficult to predict accurately the proportions contributed by each of the two mechanisms to the overall material loss. However, the overall erosion rate, combining the cutting and plastic deformation wear mechanisms, may be then given by the following equation [1–3]: πΎπœ€=3πœŒπ‘šπœŒπ‘1/2𝑉3sin3π›½πœŽπ‘¦3/2,(5) where 𝐾3 is a constant pertaining to erosion modelling of ductile metal surfaces impacted by fly ash.

2.4. Model for Erosion-Oxidation Interaction

Material removal from a brittle oxide surface on a ductile substrate at high temperature is due to a composite mechanism, involving fracture, chipping, oxide spallation at the metal-oxide interface, cutting wear, and plastic deformation of metal substrate. Following the general oxidation law, the thickness of the oxide layer 𝛿” formed during time ξ…’πœƒβ€ is given as [6, 7]𝛿=πΎπœƒπ‘›,(6) where K and n are mechanism dependent kinetic parameters. K generally being related to temperature, through Arrhenius relationship: 𝐾=𝐴𝑒(βˆ’π‘„/𝑅𝑇),(7) where A is an experimentally determined Arrhenius constant (preexponential factor), Q is the activation energy, R is the universal gas constant, and T is the absolute temperature.

Erosion can be described as a stochastic phenomenon. In the spatial domain, individual erosion events take place in a random manner. It is assumed that the instantaneous number of events on the entire surface can increase deterministically as a linear function of time. The average oxide-scale thickness, βŸ¨π›ΏβŸ©, under the influence of erosion, is given by the following dimensionless equation [6, 7]: βŸ¨π›ΏβŸ©=π‘›πΎξ‚€πΉπœ‹πœŽ2𝑛𝛾𝑛,πΉπœ‹πœŽ2πœƒξ‚,(8) where πΉπœ‹πœŽ2πœƒ represents dimensionless time, F is the average dimensional flux of impinging erodent, and 𝜎 is the radius of erosion footprint. 𝛾” is the incomplete gamma function and is given by [6, 7] ξ€œπ›Ύ(πœ‰,π‘₯)=π‘₯0π‘’πœ‰βˆ’1exp(βˆ’π‘’)𝑑𝑒.(9) The following mass loss equation is derived from the above equation to quantify the erosion under the influence of oxidation: Δ𝑀(πœƒ)=1βˆ’π‘“π‘šξ‚πœŒ0βŸ¨π›ΏβŸ©βˆ’π‘“π‘šπΉπœ‹πœŽ2𝜌0ξ€œπœƒ0βŸ¨π›ΏβŸ©π‘‘πœƒ,(10) where Δ𝑀(πœƒ), π‘“π‘š, and 𝜌0 are total erosion loss, mass fraction of oxide associated with the steel/alloy, and densities of individual oxides (NiO, Fe2O3, and Cr2O3), respectively.

3. Results

The erosion rates have been computed as function of ash particle impingement velocity, impingement angle, percentage silica in the ash sample, average density of ash particles, density of the steel component, yield stress of the steel component, and temperature of the steel component. The computations have been carried out for three boiler grade steels, namely, carbon steel, 1.25 Cr-1Mo-V steel, and alloy steel 800.

Some typical results of model predictions are discussed here. Figure 1 depicts the erosion rate (mg/kg of erodent) as a function of impact velocity for a given impingement angle of 30Β° and at room temperature for 1.25Cr-1Mo-V steel. The simulated result is validated with the published literature [1, 2] and is found to be in good agreement. Figure 2 shows the erosion rate as a function of impingement angle for carbon steel at different temperatures of the substrate, that is, 303, 573, and 873 K, respectively. The prediction of temperature dependent erosion is also validated with the published data [1, 2]. It is observed that for low values of impingement angle, the erosion rate increases with an increase in the impingement angle, and the maximum erosion rate occurs with impingement angle between 20°–40Β°. Thereafter, the erosion rate decreases. This has also been validated with the published data. The erosion rate with variation of impact velocity for the three steels is shown in Figure 3, while Figure 4 shows the change in erosion behavior with change in the temperature of the substrate. It is observed that 1.25 Cr-1Mo-V steel shows least erosion rate (or maximum erosion resistance), whereas the alloy steel 800 has the maximum erosion rate (or least erosion resistance). Figure 5 shows the erosion behavior with respect to silica content of 55% and 70% in the ash. It is seen that an increase of 15% silica in the ash results in 2–3 times enhancement in the erosion rate. A mass loss graph is shown in Figure 6 depicting dimensionless mass change as a function of dimensionless time. Erosion losses of three different oxides, namely, chromium, iron, and nickel oxides for alloy steel 800 are shown under solid particles impact. The quantitative mass loss prediction is verified with the published data [6, 7].

4. Conclusion

A mathematical model has been developed to predict the erosion rate for fly ash particles impingement on typical boiler grade steels. Stochastic approach is employed to model erosion-oxidation interaction phenomena. The following conclusions are drawn:

(i)erosion rate increases monotonically with an increase in impact velocity for a given impingement angle and temperature of the substrate;(ii)erosion rate is maximum for impact angles in the range of 20°–40Β° for a given velocity and temperature;(iii)the erosion rate is significantly enhanced with an increase in silica content of the ash;(vi)erosion rate of nickel oxide is faster with regard to iron and chromium oxides, which is attributed to the relatively faster growth of nickel oxide with respect to other oxides.