# Soft Sensor for Oxide Scales on the Steam Side of Superheater Tubes under Uneven Circumferential Load

Authors: Qing Wei Li; Gui Huan Yao
Journal: Advances in Materials Science and Engineering
 (2015)

---

## Abstract

A soft sensor for oxide scales on the steam side of superheater tubes of utility boiler under uneven circumferential loading is proposed for the first time. First finite volume method is employed to simulate oxide scales growth temperature on the steam side of superheater tube. Then appropriate time and spatial intervals are selected to calculate oxide scales thickness along the circumferential direction. On the basis of the oxide scale thickness, the stress of oxide scales is calculated by the finite element method. At last, the oxide scale thickness and stress sensors are established on support vector machine (SMV) optimized by particle swarm optimization (PSO) with time and circumferential angles as inputs and oxide scale thickness and stress as outputs. Temperature and stress calculation methods are validated by the operation data and experimental data, respectively. The soft sensor is applied to the superheater tubes of some power plant. Results show that the soft sensor can give enough accurate results for oxide scale thickness and stress in reasonable time. The forecasting model provides a convenient way for the research of the oxide scale failure.


---

## Body

## 1. Introduction

Oxide scale formed on the steam side of superheaters and reheaters of utility boilers in power plant brings many problems like tube clogging, overheating, and erosion of turbine blade [1–3]. It is estimated that 10% of all power plant breakdowns are caused by creep fractures of boiler tubes due to the scales formation [4]. Oxide failure often occurs at boiler’s start-up and stop stages. During these stages, large stress in oxide scales would be generated [5]. Many researchers have been devoted to studying the stress generated in oxide scales. For theoretical analysis, Sabau and Wright [6] studied the effects of creep, growth stress, and thermal stress on the internal stress state of the oxide film. Feng et al. [7] corrected the stoney formula for the nonuniform temperature and multilayered oxides to calculate the average stress combined with the deformation curvature of the layers. As for the experiment aspect, Galerie et al. [8] measured the average thermal stress of oxides on ferritic alloy with flat structure in three ways, namely, Raman spectroscopy, X-ray diffraction, and single-sided bending test. They pointed out that the results measured by Raman spectroscopy were accurate enough without the modifications by the X-ray diffraction and single-sided bending test. Luzin et al. [9] measured the thermal stress in Cu/Al coating systems of the flat structure using neutron diffraction method, which can give results reflecting the distribution along the thickness direction. For numerical aspect, Bian et al. [10] simulated the effects of temperature perturbations of gas and steam on the stress states of tubes and oxides using ANSYS finite element software.As the load along the circumferential direction of the superheater tube is uneven, the oxide growth temperature is different circumferentially and so is the thickness of the oxide scales. Sabau et al. [11] pointed out that the steam temperature gradually increases along the flow direction of the tube and calculated the growth temperature distribution in the axial direction.Support vector machine (SVM) introduced by Vapnik is a useful tool for data mining, especially in the fields of pattern recognition and regression. During the past few years, its solid theoretical foundation and good behaviors have attracted a number of researchers [12, 13]. However, the selection of parameters of SVM is still a problem. Grid search method [14] and optimization methods like particle swarm optimization [15], genetic algorithm [16], and ant colony optimization [17] have been introduced to solve the parameter selection problem.PSO is an evolutionary computation technique developed by Dr. Eberhart and Dr. Kennedy in 1995, inspired by social behavior of bird flocking. It is getting more and more attention as it is simple, fast, and easy to implement [18]. Like other intelligent algorithms, PSO is troubled with prematurity and convergence problem. So some improved PSOs have been proposed [19–21].It is not convenient to get oxide scales thickness and stress values through theoretical method, experimental method, and simulation method for online application. In this paper, a soft sensor for oxide scales based on SVM optimized by PSO is proposed. The rest of the paper is organized as follow. Section2 introduces the calculation of the thickness and stress of oxide scales. Then Section 3 provides a soft sensor for oxide scales, after which the temperature and stress calculation method are validated by the operation data and experiment data, respectively. An application is given to show the performance of the soft sensor in Section 5. Finally, Section 6 summarizes the contribution of this paper.


## 2. Calculation of the Thickness and Stress of Oxide Scales

### 2.1. Physical Model

The model area is divided into four regions, that is, steam region, scale region, steel region, and gas region from inside to outside, as shown in Figure1. Define the angle of the windward side as 0° and take counterclockwise as positive.Figure 1
Geometric model.
(a)
Cross section
(b)
Longitudinal sectionThe steam is fully developed turbulent flow and convection coefficient between steam and oxide scale is expressed as(1)
h
s
=
0.023
k
s
D
R
e
s
0.8
P
r
s
0.4
,where R
e
s is Reynolds number defined in (2); P
r
s is Prandtl number defined in (3); k
s is thermal conductivity of the steam, W/(m·K); D is the inner diameter of the tube, m:(2)
R
e
s
=
4
m
s
π
D
μ
s
,where m
s is mass flow rate of the steam, kg/s; μ
s is steam viscosity, kg/(m·s):(3)
P
r
s
=
μ
s
Cp
s
k
s
,where Cp
s is specific heat of the steam, J/(kg·K).Heat transfer between gas and steel is forced convection and convection coefficient is given as(4)
h
g
=
0.33
12
k
s
d
0
R
e
g
0.6
P
r
g
0.33
,where R
e
g is Reynolds number that is defined as in (5); P
r
g is Prandtl number that is defined as in (6); k
g is thermal conductivity of flue gas, W/(m·K); d
0 is outer diameter of the tube, m:(5)
R
e
g
=
G
d
0
12
μ
g
,where μ
g is steam viscosity, kg/(m·s); G is gas mass velocity, kg/s:(6)
P
r
g
=
μ
g
Cp
g
k
g
,where C
p
g is specific heat of the flue gas, J
/
(
kg
·
K
).The temperature can be solved by the following equation [25]:(7)
ρ
Cp
∂
T
∂
τ
+
u
∂
T
∂
x
+
v
∂
T
∂
y
=
λ
∂
2
T
∂
x
2
+
∂
2
T
∂
y
2
,where ρ is the density, kg/m3; Cp is the specific heat, J/(kg·K); τ is time, s; λ is thermal conductivity, W/(m·K); x and y are the space variables, m.


### 2.2. Mathematical Model of Oxide Thickness Calculation

Oxide scale thickness is the basis of the calculation of stress. For the case where protective oxide scale is formed, the rate of scale growth is controlled by diffusion of ions through the oxide, and the growth of the oxide scales can be described by [5] (8)
d
ox
2
=
2
×
k
p
×
t
k
p
=
A
×
e
-
Q
/
R
T
,where d
ox is the thickness of the oxide scales, μm; k
p is the oxidation rate, μm2/hr; t is time, hr; A is Arrhenius constant, μm2/hr, Q is activation energy, kJ/mole; R is universal gas constant, kJ/(mole·K); T is absolute growth temperature, K.Take the temperature between the tube and oxide scales as the growth temperature. The growth temperature changes with the thickness of the oxide scales, so the thickness of the oxide scales should be calculated through [26](9)
d
ox
,
k
2
=
d
ox
,
k
-
1
2
+
d
d
ox
2
d
d
ox
2
=
2
Δ
t
k
p
T
gr
t
k
,where d
ox
,
k is the thickness of the oxides at the end of the kth time intervals, μm; Δ
t is the time interval, hr; T
gr
t
k is the growth temperature in the kth interval, K; k
p
T
gr
t
k is the oxidation rate in the kth interval, μm2/hr.


### 2.3. Mathematical Model of the Stress in Oxide Scales

The stress in oxide scales mainly consists of thermal stress, which can be calculated by the following equations [27]:(10)
ε
met
=
∫
T
a
T
b
β
met
(
T
)
d
T
ε
oxide
=
∫
T
c
T
d
β
oxide
(
T
)
d
T
,where ε
met and ε
oxide are the free strain of tube and oxide scale, respectively.The displacement and stress of tube and oxide scales are calculated according to mechanics of elasticity as follows:(11)
u
M
=
1
+
ν
M
r
M
∫
r
a
r
b
1
+
ν
M
ε
M
r
d
r
+
C
M
r
M
+
D
M
r
M
σ
r
,
M
=
E
M
-
1
r
M
2
∫
r
a
r
b
(
1
+
ν
M
)
ε
M
r
d
r
vvvvvvv
+
C
M
1
-
ν
M
-
D
M
(
1
+
ν
M
)
r
M
2
σ
φ
,
M
=
E
M
1
r
M
2
∫
r
a
r
b
(
1
+
ν
M
)
ε
M
r
d
r
vvvvvvv
-
(
1
+
ν
M
)
ε
M
+
C
M
1
-
ν
M
+
D
M
(
1
+
ν
M
)
r
M
2
σ
z
,
M
=
ν
M
σ
r
,
M
+
σ
φ
,
M
-
E
M
ε
M
,where ε
met is the metal free thermal strain; ε
oxide is the oxide free thermal strain; β
met is the linear expansion coefficient of the metal, m/K; β
oxide is the linear expansion coefficient of the oxide scales, m/K; T
a and T
b are the inside and outside temperature of metal, respectively, K; T
c and T
d are the inside and outside temperature of oxide scales, respectively, K; M represents metal or oxide scales; u is the radial displacement, m; ν is the Poisson ration; r
M is the average radial coordinate, m; r is the radial argument, m; r
a and r
b are the inside radius and outside radius, m; ε
M is the free thermal strain; C
M and D
M are undetermined parameters and are obtained with boundary conditions; σ
r
,
M is the radial stress, MPa; E
M is modulus of elasticity, Mpa; σ
φ
,
M is circumferential stress, MPa; σ
z
,
M is axial stress, MPa.The inlet gas temperature and steam temperature are given. The radial boundary stress of the gas side and steam side equals the gas pressure and steam pressure, respectively. The radial stress and temperature are continuous at the interface of tube and oxide scales and the interfaces of different oxide layers.


## 2.1. Physical Model

The model area is divided into four regions, that is, steam region, scale region, steel region, and gas region from inside to outside, as shown in Figure1. Define the angle of the windward side as 0° and take counterclockwise as positive.Figure 1
Geometric model.
(a)
Cross section
(b)
Longitudinal sectionThe steam is fully developed turbulent flow and convection coefficient between steam and oxide scale is expressed as(1)
h
s
=
0.023
k
s
D
R
e
s
0.8
P
r
s
0.4
,where R
e
s is Reynolds number defined in (2); P
r
s is Prandtl number defined in (3); k
s is thermal conductivity of the steam, W/(m·K); D is the inner diameter of the tube, m:(2)
R
e
s
=
4
m
s
π
D
μ
s
,where m
s is mass flow rate of the steam, kg/s; μ
s is steam viscosity, kg/(m·s):(3)
P
r
s
=
μ
s
Cp
s
k
s
,where Cp
s is specific heat of the steam, J/(kg·K).Heat transfer between gas and steel is forced convection and convection coefficient is given as(4)
h
g
=
0.33
12
k
s
d
0
R
e
g
0.6
P
r
g
0.33
,where R
e
g is Reynolds number that is defined as in (5); P
r
g is Prandtl number that is defined as in (6); k
g is thermal conductivity of flue gas, W/(m·K); d
0 is outer diameter of the tube, m:(5)
R
e
g
=
G
d
0
12
μ
g
,where μ
g is steam viscosity, kg/(m·s); G is gas mass velocity, kg/s:(6)
P
r
g
=
μ
g
Cp
g
k
g
,where C
p
g is specific heat of the flue gas, J
/
(
kg
·
K
).The temperature can be solved by the following equation [25]:(7)
ρ
Cp
∂
T
∂
τ
+
u
∂
T
∂
x
+
v
∂
T
∂
y
=
λ
∂
2
T
∂
x
2
+
∂
2
T
∂
y
2
,where ρ is the density, kg/m3; Cp is the specific heat, J/(kg·K); τ is time, s; λ is thermal conductivity, W/(m·K); x and y are the space variables, m.


## 2.2. Mathematical Model of Oxide Thickness Calculation

Oxide scale thickness is the basis of the calculation of stress. For the case where protective oxide scale is formed, the rate of scale growth is controlled by diffusion of ions through the oxide, and the growth of the oxide scales can be described by [5] (8)
d
ox
2
=
2
×
k
p
×
t
k
p
=
A
×
e
-
Q
/
R
T
,where d
ox is the thickness of the oxide scales, μm; k
p is the oxidation rate, μm2/hr; t is time, hr; A is Arrhenius constant, μm2/hr, Q is activation energy, kJ/mole; R is universal gas constant, kJ/(mole·K); T is absolute growth temperature, K.Take the temperature between the tube and oxide scales as the growth temperature. The growth temperature changes with the thickness of the oxide scales, so the thickness of the oxide scales should be calculated through [26](9)
d
ox
,
k
2
=
d
ox
,
k
-
1
2
+
d
d
ox
2
d
d
ox
2
=
2
Δ
t
k
p
T
gr
t
k
,where d
ox
,
k is the thickness of the oxides at the end of the kth time intervals, μm; Δ
t is the time interval, hr; T
gr
t
k is the growth temperature in the kth interval, K; k
p
T
gr
t
k is the oxidation rate in the kth interval, μm2/hr.


## 2.3. Mathematical Model of the Stress in Oxide Scales

The stress in oxide scales mainly consists of thermal stress, which can be calculated by the following equations [27]:(10)
ε
met
=
∫
T
a
T
b
β
met
(
T
)
d
T
ε
oxide
=
∫
T
c
T
d
β
oxide
(
T
)
d
T
,where ε
met and ε
oxide are the free strain of tube and oxide scale, respectively.The displacement and stress of tube and oxide scales are calculated according to mechanics of elasticity as follows:(11)
u
M
=
1
+
ν
M
r
M
∫
r
a
r
b
1
+
ν
M
ε
M
r
d
r
+
C
M
r
M
+
D
M
r
M
σ
r
,
M
=
E
M
-
1
r
M
2
∫
r
a
r
b
(
1
+
ν
M
)
ε
M
r
d
r
vvvvvvv
+
C
M
1
-
ν
M
-
D
M
(
1
+
ν
M
)
r
M
2
σ
φ
,
M
=
E
M
1
r
M
2
∫
r
a
r
b
(
1
+
ν
M
)
ε
M
r
d
r
vvvvvvv
-
(
1
+
ν
M
)
ε
M
+
C
M
1
-
ν
M
+
D
M
(
1
+
ν
M
)
r
M
2
σ
z
,
M
=
ν
M
σ
r
,
M
+
σ
φ
,
M
-
E
M
ε
M
,where ε
met is the metal free thermal strain; ε
oxide is the oxide free thermal strain; β
met is the linear expansion coefficient of the metal, m/K; β
oxide is the linear expansion coefficient of the oxide scales, m/K; T
a and T
b are the inside and outside temperature of metal, respectively, K; T
c and T
d are the inside and outside temperature of oxide scales, respectively, K; M represents metal or oxide scales; u is the radial displacement, m; ν is the Poisson ration; r
M is the average radial coordinate, m; r is the radial argument, m; r
a and r
b are the inside radius and outside radius, m; ε
M is the free thermal strain; C
M and D
M are undetermined parameters and are obtained with boundary conditions; σ
r
,
M is the radial stress, MPa; E
M is modulus of elasticity, Mpa; σ
φ
,
M is circumferential stress, MPa; σ
z
,
M is axial stress, MPa.The inlet gas temperature and steam temperature are given. The radial boundary stress of the gas side and steam side equals the gas pressure and steam pressure, respectively. The radial stress and temperature are continuous at the interface of tube and oxide scales and the interfaces of different oxide layers.


## 3. Soft Sensor for the Thickness and Stress of Oxide Scales

Though the stresses can be calculated by simulation method, it is time consuming. Operators are eager to monitor the oxide thickness and stress online. So a soft sensor is proposed.


### 3.1. Introduction of SVM

Train samples of SVM are{
(
x
1
,
y
1
)
,
…
,
(
x
n
,
y
n
)
}
⊂
R
n
×
R, which are mapped into high-dimensional space through nonlinear function Φ and regressed linearly as follows:(12)
f
x
=
w
Φ
x
+
b
,where w
∈
R
n; b
∈
R.The learning process is converted into an optimization problem according to structural risk minimization principle:(13)
min
⁡
1
2
w
2
+
C
∑
i
=
1
n
(
ξ
i
+
ξ
*
)
,
s
.
t
.
y
i
-
f
(
x
i
)
≤
ε
+
ξ
i
,
vvvvv
f
(
x
i
)
-
y
i
≤
ε
+
ξ
i
*
vvvvv
ξ
i
,
ξ
i
*
≥
0
,
i
=
1,2
,
…
,
n
,where C is the penalty factor; ξ
i and ξ
i
* are slack variables.Linear insensitive loss function is defined withε:(14)
L
ε
x
,
y
,
f
=
max
⁡
o
,
y
-
f
x
-
ε
.Equation (13) can be transformed into the following form through the dual form of Lagrange polynomial:(15)
max
⁡
∑
i
=
1
n
y
i
(
α
i
*
-
α
i
)
-
ε
∑
i
=
1
n
(
α
i
*
+
α
i
)
vvvvvv
-
1
2
∑
i
,
j
=
1
n
(
α
i
*
-
α
i
)
(
α
j
*
-
α
j
)
K
(
x
i
,
x
j
)
s
.
t
.
∑
i
n
α
i
*
=
∑
i
n
α
i
,
vvvvv
C
≥
α
i
*
,
α
i
≥
0
,
i
=
1,2
,
…
,
n
,where α
i
* and α
i are Lagrange multipliers; K
(
x
i
,
x
j
) is kernel function meeting the conditions of mercer. Radial basis kernel function is selected in this study:(16)
K
x
i
,
x
j
=
exp
⁡
-
x
i
-
x
j
2
σ
2
.Due to the sparsity, only some samples’ coefficients are not 0 in quadratic programming (15), which are the support vector machines. Assuming the number of support vector machines is g, the regression function can be given as follows: (17)
f
x
=
∑
i
=
1
g
α
i
*
-
α
i
K
x
,
x
i
+
b
.


### 3.2. Introduction of PSO

In PSO, each feasible solution of the optimization problem is seen as a “particle” of the solution space. Each particle searches the solution space following the optimal particle of the whole group. Suppose the group searches in aD dimensional solution space. The particle i is expressed as x
(
x
i
1
,
x
i
2
,
…
,
x
i
d
). Particles search for new feasible solutions by constantly adjusting their positions. P
i
(
P
i
1
,
P
i
2
,
…
,
P
i
d
) is the optimal position of particle i. P
g
(
P
g
1
,
P
g
2
,
…
,
P
g
d
) is the optimal position of the whole group. V
i
(
V
i
1
,
V
i
2
,
…
,
V
i
d
) is position change rate of particle i. Each particle uses the following equations to adjust its velocity and position: (18)
V
i
d
t
+
1
=
w
V
i
d
t
+
c
1
×
rand
1
×
p
i
d
-
x
i
d
t
+
c
2
×
rand
2
×
p
g
d
-
x
i
d
t
x
i
d
t
+
1
=
x
i
d
t
+
V
i
d
t
,
1
≤
i
≤
N
,
1
≤
d
≤
D
,where w is inertia weight; D is the dimension of the space; N is the particle number in the group; c
1 and c
2 are acceleration coefficients, where c
1 is known as the individual cognitive factor; c
2 is known as social learning factor; rand1 and rand2 are two separately generated uniformly distributed random numbers in the range (
0,1
).


### 3.3. Optimization of SVM’s Parameters

C and σ have great influence on the learning ability and generalization ability. C compromises structural risks and sampling error. σ is decided by the width or the scope of the input space.In this paper, support vector machine parameters are optimized by PSO. Each particle represents a potential set ofC and σ. The fitness function is defined as(19)
MSE
=
1
N
∑
i
=
1
N
f
x
i
-
y
i
2
,where f
(
x
i
) is the forecasting value and y
i is the actual value.Initialize PSO and iterate to get proper set ofC and σ to minimize the predicting error, namely, the fitness value evaluated by (19), to promote the predicting accuracy.


### 3.4. Soft Sensor for Thickness and Stress of Oxide Scales

Traditional methods like theoretical method, experimental method, and simulation method cannot provide accurate enough oxide scale thickness and stress in reasonable time. With the development of intelligent algorithms, the soft sensor can be qualified to do this work. In our study, SVM is employed to predict oxide scale thickness and stress.(1)Soft Sensor for Oxide Scale Thickness. Take time and angles as inputs and oxide scale thickness as the output, as shown in Figure 2(a). Train samples and test samples are obtained from simulations. Take 6 out of the samples as train samples and the one left is test sample. Leave-one-out method is employed to train SVM parameters [16].Figure 2
Soft sensor for oxide scale thickness and stress.
(a)
Soft sensor for oxide scale thickness
(b)
Soft sensor for oxide scale stress(2)Soft Sensor for Oxide Scale Stress. Take time, oxide scale thickness, and angles as inputs and oxide scale stress as the output, as shown in Figure 2(b). Train samples and test samples are obtained from simulations. The first 18 samples are train samples, while the last 3 samples are test samples. The training method is the same with thickness soft sensor. Leave-one-out method is employed to train SVM parameters.The steps of constructing the soft sensors are detailed as follows.Step 1.
Collect calculation parameters of oxide scales and superheater tubes.Step 2.
Program thermal calculation procedure according to 73 boiler thermal calculation standards to get the convection coefficient between steam and oxide scales and convection coefficient between gas and substrate and steam temperature and gas temperature.Step 3.
Divide operating time and spatial space into several intervals reasonably considering computational cost and accuracy of the results.Step 4.
Calculate the temperature distribution by finite volume method with the boundaries got in Step2.Step 5.
Calculate the oxide scale thickness according to (8)-(9) with intervals in Step 3 and temperature distribution in Step 4. Repeat Steps 4 and 5 until all the intervals have been done.Step 6.
Calculate the stress of oxide scales according to (10)-(11) on the platform of ANSYS with given temperature drop and oxide scale thickness in corresponding direction.Step 7.
Optimize SVM parameters for the oxide scale thickness soft sensor and stress soft sensor with leave-one-out method.Step 8.
Construct thickness soft sensor and stress soft sensor with the data collecting from Steps5 and 6 and the parameters optimized in Step 7.


## 3.1. Introduction of SVM

Train samples of SVM are{
(
x
1
,
y
1
)
,
…
,
(
x
n
,
y
n
)
}
⊂
R
n
×
R, which are mapped into high-dimensional space through nonlinear function Φ and regressed linearly as follows:(12)
f
x
=
w
Φ
x
+
b
,where w
∈
R
n; b
∈
R.The learning process is converted into an optimization problem according to structural risk minimization principle:(13)
min
⁡
1
2
w
2
+
C
∑
i
=
1
n
(
ξ
i
+
ξ
*
)
,
s
.
t
.
y
i
-
f
(
x
i
)
≤
ε
+
ξ
i
,
vvvvv
f
(
x
i
)
-
y
i
≤
ε
+
ξ
i
*
vvvvv
ξ
i
,
ξ
i
*
≥
0
,
i
=
1,2
,
…
,
n
,where C is the penalty factor; ξ
i and ξ
i
* are slack variables.Linear insensitive loss function is defined withε:(14)
L
ε
x
,
y
,
f
=
max
⁡
o
,
y
-
f
x
-
ε
.Equation (13) can be transformed into the following form through the dual form of Lagrange polynomial:(15)
max
⁡
∑
i
=
1
n
y
i
(
α
i
*
-
α
i
)
-
ε
∑
i
=
1
n
(
α
i
*
+
α
i
)
vvvvvv
-
1
2
∑
i
,
j
=
1
n
(
α
i
*
-
α
i
)
(
α
j
*
-
α
j
)
K
(
x
i
,
x
j
)
s
.
t
.
∑
i
n
α
i
*
=
∑
i
n
α
i
,
vvvvv
C
≥
α
i
*
,
α
i
≥
0
,
i
=
1,2
,
…
,
n
,where α
i
* and α
i are Lagrange multipliers; K
(
x
i
,
x
j
) is kernel function meeting the conditions of mercer. Radial basis kernel function is selected in this study:(16)
K
x
i
,
x
j
=
exp
⁡
-
x
i
-
x
j
2
σ
2
.Due to the sparsity, only some samples’ coefficients are not 0 in quadratic programming (15), which are the support vector machines. Assuming the number of support vector machines is g, the regression function can be given as follows: (17)
f
x
=
∑
i
=
1
g
α
i
*
-
α
i
K
x
,
x
i
+
b
.


## 3.2. Introduction of PSO

In PSO, each feasible solution of the optimization problem is seen as a “particle” of the solution space. Each particle searches the solution space following the optimal particle of the whole group. Suppose the group searches in aD dimensional solution space. The particle i is expressed as x
(
x
i
1
,
x
i
2
,
…
,
x
i
d
). Particles search for new feasible solutions by constantly adjusting their positions. P
i
(
P
i
1
,
P
i
2
,
…
,
P
i
d
) is the optimal position of particle i. P
g
(
P
g
1
,
P
g
2
,
…
,
P
g
d
) is the optimal position of the whole group. V
i
(
V
i
1
,
V
i
2
,
…
,
V
i
d
) is position change rate of particle i. Each particle uses the following equations to adjust its velocity and position: (18)
V
i
d
t
+
1
=
w
V
i
d
t
+
c
1
×
rand
1
×
p
i
d
-
x
i
d
t
+
c
2
×
rand
2
×
p
g
d
-
x
i
d
t
x
i
d
t
+
1
=
x
i
d
t
+
V
i
d
t
,
1
≤
i
≤
N
,
1
≤
d
≤
D
,where w is inertia weight; D is the dimension of the space; N is the particle number in the group; c
1 and c
2 are acceleration coefficients, where c
1 is known as the individual cognitive factor; c
2 is known as social learning factor; rand1 and rand2 are two separately generated uniformly distributed random numbers in the range (
0,1
).


## 3.3. Optimization of SVM’s Parameters

C and σ have great influence on the learning ability and generalization ability. C compromises structural risks and sampling error. σ is decided by the width or the scope of the input space.In this paper, support vector machine parameters are optimized by PSO. Each particle represents a potential set ofC and σ. The fitness function is defined as(19)
MSE
=
1
N
∑
i
=
1
N
f
x
i
-
y
i
2
,where f
(
x
i
) is the forecasting value and y
i is the actual value.Initialize PSO and iterate to get proper set ofC and σ to minimize the predicting error, namely, the fitness value evaluated by (19), to promote the predicting accuracy.


## 3.4. Soft Sensor for Thickness and Stress of Oxide Scales

Traditional methods like theoretical method, experimental method, and simulation method cannot provide accurate enough oxide scale thickness and stress in reasonable time. With the development of intelligent algorithms, the soft sensor can be qualified to do this work. In our study, SVM is employed to predict oxide scale thickness and stress.(1)Soft Sensor for Oxide Scale Thickness. Take time and angles as inputs and oxide scale thickness as the output, as shown in Figure 2(a). Train samples and test samples are obtained from simulations. Take 6 out of the samples as train samples and the one left is test sample. Leave-one-out method is employed to train SVM parameters [16].Figure 2
Soft sensor for oxide scale thickness and stress.
(a)
Soft sensor for oxide scale thickness
(b)
Soft sensor for oxide scale stress(2)Soft Sensor for Oxide Scale Stress. Take time, oxide scale thickness, and angles as inputs and oxide scale stress as the output, as shown in Figure 2(b). Train samples and test samples are obtained from simulations. The first 18 samples are train samples, while the last 3 samples are test samples. The training method is the same with thickness soft sensor. Leave-one-out method is employed to train SVM parameters.The steps of constructing the soft sensors are detailed as follows.Step 1.
Collect calculation parameters of oxide scales and superheater tubes.Step 2.
Program thermal calculation procedure according to 73 boiler thermal calculation standards to get the convection coefficient between steam and oxide scales and convection coefficient between gas and substrate and steam temperature and gas temperature.Step 3.
Divide operating time and spatial space into several intervals reasonably considering computational cost and accuracy of the results.Step 4.
Calculate the temperature distribution by finite volume method with the boundaries got in Step2.Step 5.
Calculate the oxide scale thickness according to (8)-(9) with intervals in Step 3 and temperature distribution in Step 4. Repeat Steps 4 and 5 until all the intervals have been done.Step 6.
Calculate the stress of oxide scales according to (10)-(11) on the platform of ANSYS with given temperature drop and oxide scale thickness in corresponding direction.Step 7.
Optimize SVM parameters for the oxide scale thickness soft sensor and stress soft sensor with leave-one-out method.Step 8.
Construct thickness soft sensor and stress soft sensor with the data collecting from Steps5 and 6 and the parameters optimized in Step 7.


## 4. Verification

### 4.1. Validation of Oxide Scale Thickness Calculation Method

Two reheater tube samples of power plant in site are taken for validation [28, 29]. Two different cases with different tube diameters from two different locations are used. Operating steam temperature of both tubes is 576°C. The flue gas temperatures were reported ranging from 800 to 900°C. The estimated scale thickness and the actual data are plotted in Figure 3. It can be seen that the estimated scale thickness is closed to the actual data.Figure 3
The estimated scale thickness and the actual data.


### 4.2. Validation of Oxide Scale Stress Calculation Method

The average internal stress of oxides on AISI 441 is measured by Raman spectroscopy in [30]. The dimensions of the specimen are 2 mm × 1.5 mm. The thickness of the alloy and oxide scales is 1.2 mm and 0.001 mm, respectively. This literature assumed that the composition of the oxides was chromia only. The specimen was cooled to room temperature after being oxidized at 800°C. The pressure is 1 atm. The average stress of the oxide scales is simulated with finite element method. The simulated result is compressive and the value is 1.126 Gpa. The error is 6.17% compared with 1.2 Gpa provided in [30]. So the results based on (10)-(11) by finite element method are credible.


## 4.1. Validation of Oxide Scale Thickness Calculation Method

Two reheater tube samples of power plant in site are taken for validation [28, 29]. Two different cases with different tube diameters from two different locations are used. Operating steam temperature of both tubes is 576°C. The flue gas temperatures were reported ranging from 800 to 900°C. The estimated scale thickness and the actual data are plotted in Figure 3. It can be seen that the estimated scale thickness is closed to the actual data.Figure 3
The estimated scale thickness and the actual data.


## 4.2. Validation of Oxide Scale Stress Calculation Method

The average internal stress of oxides on AISI 441 is measured by Raman spectroscopy in [30]. The dimensions of the specimen are 2 mm × 1.5 mm. The thickness of the alloy and oxide scales is 1.2 mm and 0.001 mm, respectively. This literature assumed that the composition of the oxides was chromia only. The specimen was cooled to room temperature after being oxidized at 800°C. The pressure is 1 atm. The average stress of the oxide scales is simulated with finite element method. The simulated result is compressive and the value is 1.126 Gpa. The error is 6.17% compared with 1.2 Gpa provided in [30]. So the results based on (10)-(11) by finite element method are credible.


## 5. An Application Case

### 5.1. Description on the Calculation Related Data

The physical parameters of oxide scales and superheater material are listed in Tables1-2. Geometrical parameters listed in Table 3 are taken from the superheater of some power plant. Take operational parameters as boundary conditions; see Table 4.Table 1
The main parameters for oxide scales [22, 23].
Name
Elasticity modulus/×1011 [Pa]
The Poisson ration
Fe-Cr-spinel
2.3
0.31
Magnetite
2.1
0.29
Hematite
2.2
0.19Table 2
The parameters of TP347 [24].
Temperature [°C]
Elasticity modulus ×1011 [Pa]
The linear expansion coefficient ×10−6 [K−1]
The Poisson
Conduction coefficient [W/(m·K)]
Arrhenius constantA [μm2/h]
Activation energy [KJ/mole]
25
2
15.27
0.3
23.2
5 × 109
−171
100
1.96
17.3
0.3
23.2
5 × 109
−171
200
1.88
17.5
0.3
23.2
5 × 109
−171
300
1.81
17.7
0.3
23.2
5 × 109
−171
400
1.72
18.2
0.3
23.2
5 × 109
−171
500
1.63
18.6
0.3
23.2
5 × 109
−171
600
1.56
18.9
0.3
23.2
5 × 109
−171Table 3
Geometrical parameters.
Inner diameter [mm]
Outer diameter [mm]
Axial length [mm]
Length of gas area [mm]
Width of gas area [mm]
Operating time [h]
36.7
44.5
15
1000
500
10000Table 4
Boundary conditions.
Inlet temperature [°C]
Inlet velocity [m/s]
Outlet pressure [MPa]
Gas
950
9.6
0.09
Steam
560
6.9
25


### 5.2. Growth Temperature and Thickness of the Oxide Scales

Growth temperature decreases along the circumferential direction and changes rapidly at approximately 90°. Reduced gas flow area enhances flow speed, which enlarges the convective heat transfer coefficient. At first, the growth temperature of the oxide scales increases fast as the difference of thermal conductivity between the metal and the oxide scales is large; see Figure4. The thickness of the oxide scales is given in Figure 5.Figure 4
The growth temperature along the tube.Figure 5
The oxide thickness along the tube.


### 5.3. Forecasting of the Oxide Scales Thickness

The predicting results for oxide scale thickness are shown in Figure6. We can see the forecasting values are closed to the actual values for train samples and test samples in all directions, showing that the results given by oxide scale thickness soft sensor are accurate enough.Figure 6
The forecasting results.
(a)
0°, 30°, and 60°
(b)
90°, 120°, 150°, and 180°


### 5.4. Forecasting of the Oxide Scale Stress

Assume that the temperature drop of steam is 100°C while all other parameters remain the same. The stresses are calculated with (10)-(11). The optimization progress of SVM parameters for the radial stress in 0° is shown in Figure 7. The best c and g are found after 38 iterations. The best c and best g are 100 and 0.87, respectively.Figure 7
The optimization progress of SVM parameters (Bestc = 100 g = 0.86866).Radial stresses, axial stresses, and circumferential stresses in 0°, 90°, and 180° are shown in Figures8, 9, and 10. The first 18 samples are training samples and the last 3 samples are testing samples. From the figures, it can be seen that the forecasting results provided by soft senor are accurate enough. The first 3 training samples’ forecasting results are not ideal, which is caused by two reasons: the samples are not enough and the stress rises too dramatically.Figure 8
Forecasting results of radial stresses.
(a)
0° direction
(b)
90° direction
(c)
180° directionFigure 9
Forecasting results of axial stress.
(a)
0° direction
(b)
90° direction
(c)
180° directionFigure 10
Forecasting results of circumferential stress.
(a)
0° direction
(b)
90° direction
(c)
180° directionThe corresponding average errors and maximum errors are shown in Tables5-6. The maximum forecasting error comes from the first training sample. The average forecasting errors are small enough to show the stress level in the oxide scales, which can be further reduced by adding samples in the early stage. To get the stress, the simulation method takes more than 10 minutes; however, the soft sensor only needs few seconds, cutting down the calculation time greatly.Table 5
The average forecasting errors for the stress.
0°
90°
180°
Axial stress average error/[MPa]
1.35
1.41
2.26
Circumferential stress average error/[MPa]
0.40
0.50
0.66
Radial stress average error/[MPa]
0.02
0.02
0.15Table 6
The maximum forecasting errors for the stress.
0°
90°
180°
Axial stress maximum error/[MPa]
76.85
40.99
24.28
Circumferential stress maximum error/[MPa]
58.38
43.45
41.32
Radial stress maximum error/[MPa]
4.39
4.46
1.99


## 5.1. Description on the Calculation Related Data

The physical parameters of oxide scales and superheater material are listed in Tables1-2. Geometrical parameters listed in Table 3 are taken from the superheater of some power plant. Take operational parameters as boundary conditions; see Table 4.Table 1
The main parameters for oxide scales [22, 23].
Name
Elasticity modulus/×1011 [Pa]
The Poisson ration
Fe-Cr-spinel
2.3
0.31
Magnetite
2.1
0.29
Hematite
2.2
0.19Table 2
The parameters of TP347 [24].
Temperature [°C]
Elasticity modulus ×1011 [Pa]
The linear expansion coefficient ×10−6 [K−1]
The Poisson
Conduction coefficient [W/(m·K)]
Arrhenius constantA [μm2/h]
Activation energy [KJ/mole]
25
2
15.27
0.3
23.2
5 × 109
−171
100
1.96
17.3
0.3
23.2
5 × 109
−171
200
1.88
17.5
0.3
23.2
5 × 109
−171
300
1.81
17.7
0.3
23.2
5 × 109
−171
400
1.72
18.2
0.3
23.2
5 × 109
−171
500
1.63
18.6
0.3
23.2
5 × 109
−171
600
1.56
18.9
0.3
23.2
5 × 109
−171Table 3
Geometrical parameters.
Inner diameter [mm]
Outer diameter [mm]
Axial length [mm]
Length of gas area [mm]
Width of gas area [mm]
Operating time [h]
36.7
44.5
15
1000
500
10000Table 4
Boundary conditions.
Inlet temperature [°C]
Inlet velocity [m/s]
Outlet pressure [MPa]
Gas
950
9.6
0.09
Steam
560
6.9
25


## 5.2. Growth Temperature and Thickness of the Oxide Scales

Growth temperature decreases along the circumferential direction and changes rapidly at approximately 90°. Reduced gas flow area enhances flow speed, which enlarges the convective heat transfer coefficient. At first, the growth temperature of the oxide scales increases fast as the difference of thermal conductivity between the metal and the oxide scales is large; see Figure4. The thickness of the oxide scales is given in Figure 5.Figure 4
The growth temperature along the tube.Figure 5
The oxide thickness along the tube.


## 5.3. Forecasting of the Oxide Scales Thickness

The predicting results for oxide scale thickness are shown in Figure6. We can see the forecasting values are closed to the actual values for train samples and test samples in all directions, showing that the results given by oxide scale thickness soft sensor are accurate enough.Figure 6
The forecasting results.
(a)
0°, 30°, and 60°
(b)
90°, 120°, 150°, and 180°


## 5.4. Forecasting of the Oxide Scale Stress

Assume that the temperature drop of steam is 100°C while all other parameters remain the same. The stresses are calculated with (10)-(11). The optimization progress of SVM parameters for the radial stress in 0° is shown in Figure 7. The best c and g are found after 38 iterations. The best c and best g are 100 and 0.87, respectively.Figure 7
The optimization progress of SVM parameters (Bestc = 100 g = 0.86866).Radial stresses, axial stresses, and circumferential stresses in 0°, 90°, and 180° are shown in Figures8, 9, and 10. The first 18 samples are training samples and the last 3 samples are testing samples. From the figures, it can be seen that the forecasting results provided by soft senor are accurate enough. The first 3 training samples’ forecasting results are not ideal, which is caused by two reasons: the samples are not enough and the stress rises too dramatically.Figure 8
Forecasting results of radial stresses.
(a)
0° direction
(b)
90° direction
(c)
180° directionFigure 9
Forecasting results of axial stress.
(a)
0° direction
(b)
90° direction
(c)
180° directionFigure 10
Forecasting results of circumferential stress.
(a)
0° direction
(b)
90° direction
(c)
180° directionThe corresponding average errors and maximum errors are shown in Tables5-6. The maximum forecasting error comes from the first training sample. The average forecasting errors are small enough to show the stress level in the oxide scales, which can be further reduced by adding samples in the early stage. To get the stress, the simulation method takes more than 10 minutes; however, the soft sensor only needs few seconds, cutting down the calculation time greatly.Table 5
The average forecasting errors for the stress.
0°
90°
180°
Axial stress average error/[MPa]
1.35
1.41
2.26
Circumferential stress average error/[MPa]
0.40
0.50
0.66
Radial stress average error/[MPa]
0.02
0.02
0.15Table 6
The maximum forecasting errors for the stress.
0°
90°
180°
Axial stress maximum error/[MPa]
76.85
40.99
24.28
Circumferential stress maximum error/[MPa]
58.38
43.45
41.32
Radial stress maximum error/[MPa]
4.39
4.46
1.99


## 6. Conclusions

The thickness and stress of oxide scales under uneven circumferential loading are calculated by finite volume method and finite element method, which are supplied as the samples. A soft sensor for oxide thickness and stress is proposed for the first time. An application shows that the soft sensor can give enough accurate results for oxide scale thickness and stress while the calculation time is greatly cut down. The soft sensor provides a convenient way to study oxide scale failure. In future, more operating parameters, tube diameters, and tube materials will be studied and incorporated into the predicting model.



---

*Source: 275340-2015-03-03.xml*