Measurements of Local Heat Flux and Water-Side Heat Transfer Coefficient in Water Wall Tubes
Measurements of Local Heat Flux and Water-Side Heat Transfer Coefficient in Water Wall Tubes
Measurements of Local Heat Flux and Water-Side Heat Transfer Coefficient in Water Wall Tubes
Wall Tubes
Jan Taler1 and Dawid Taler2
[1]
[2]
Institute of Heat Transfer Engineering and Air Protection, Cracow University of Technology, Cracow, Poland
1. Introduction
Measurements of heat flux and heat transfer coefficient are subject of many current studies. A proper understanding of combustion and heat transfer in furnaces and
heat exchange on the water-steam side in water walls requires accurate measurement of heat flux which is absorbed by membrane furnace walls. There are three
broad categories of heat flux measurements of the boiler water-walls: (1) portable heat flux meters inserted in inspection ports [1], (2) Gardon type heat flux meters
welded to the sections of the boiler tubes [1-4], (3) tubular type instruments placed between two adjacent boiler tubes [5-14]. Tubular type and Gardon meters
strategically placed on the furnace tube wall can be a valuable boiler diagnostic device for monitoring of slag deposition. If a heat flux instrument is to measure the
absorbed heat flux correctly, it must resemble the boiler tube as closely as possible so far as radiant heat exchange with the flame and surrounding surfaces is
concerned. Two main factors in this respect are the emissivity and the temperature of the absorbing surface, but since the instrument will almost always be coated
with ash, it is generally the properties of the ash and not the instrument that dominate the situation. Unfortunately, the thermal conductivity can vary widely.
Therefore, accurate measurements will only be performed if the deposit on the meter is representative of that on the surrounding tubes. The tubular type instruments
known also as flux-tubes meet this requirement. In these devices the measured boiler tube wall temperatures are used for the evaluation of the heat fluxqm. The
measuring tube is fitted with two thermocouples in holes of known radial spacing r1 and r2. The thermocouples are led away to the junction box where they are
connected differentially to give a flux related electromotive force.
The use of the one dimensional heat conduction equation for determining temperature distribution in the tube wall leads to the simple formula
qm=k(f1f2)roln(r1/r2)qm=k(f1f2)roln(r1/r2)
Options
(1)
The accuracy of this equation is very low because of the circumferential heat conduction in the tube wall.
However, the measurement of the heat flux absorbed by water-walls with satisfactory accuracy is a challenging task. Considerable work has been done in recent years
in this field.Previous attempts to accurately measure the local heat flux to membrane water walls in steam boilers failed due to calculation of inside heat transfer
coefficients. The heat flux can only be determined accurately if the inside heat transfer coefficient is measured experimentally.
New numerical methods for determining the heat flux in boiler furnaces, based on experimentally acquired interior flux-tube temperatures, will be presented. The
tubular type instruments have been designed to provide a very accurate measurement of absorbed heat flux qm, inside heat transfer coefficient hin, and water steam
temperature Tf.
Two different tubular type instruments (flux tubes) were developed to identify boundary conditions in water wall tubes of steam boilers.
The first meter is constructed from a short length of eccentric bare tube containing four thermocouples on the fire side below the inner and outer surfaces of the tube.
The fifth thermocouple is located at the rear of the tube on the casing side of the water wall tube. First, formulas for the view factor defining the heat flux distribution
at the outer surface of the flux tube were derived. The exact analytical expressions for the view factor compare very well with approximate methods for determining
view factor which are used by the ANSYS software. The meter is constructed from a short length of eccentric tube containing four thermocouples on the fireside
below the inner and outer surfaces of the tube. The fifth thermocouple is located at the rear of the tube (on the casing side of the water-wall tube). The boundary
conditions on the outer and inner surfaces of the water flux-tube must then be determined from temperature measurements at the interior locations. Four K-type
sheathed thermocouples, 1 mm in diameter, are inserted into holes, which are parallel to the tube axis. The thermal conduction effect at the hot junction is minimized
because the thermocouples pass through isothermal holes. The thermocouples are brought to the rear of the tube in the slot machined in the tube wall. An austenitic
cover plate with the thickness of 3 mm welded to the tube is used to protect the thermocouples from the incident flame radiation. A K-type sheathed thermocouple
with a pad is used to measure the temperature at the rear of the flux-tube. This temperature is almost the same as the water-steam temperature.
The non-linear least squares problem was solved numerically using the LevenbergMarquardt method. The temperature distribution at the cross section of the flux
tube was determined at every iteration step using the method of separation of variables.The heat transfer conditions in adjacent boiler tubes have no impact on the
temperature distribution in the flux tubes.
The second flux tube has two longitudinal fins which are welded to the eccentric bare tube. In contrast to existing devices, in the developed flux-tube fins are not
welded to adjacent water-wall tubes. Temperature distribution in the flux-tube is symmetric and not disturbed by different temperature fields in neighboring tubes.
The temperature dependent thermal conductivity of the flux-tube material was assumed. An inverse problem of heat conduction was solved using the least squares
method. Three unknown parameters were estimated using the Levenberg-Marquardt method. At every iteration step, the temperature distribution over the crosssection of the heat flux meter was computed using the ANSYS CFX software. Test calculations were carried out to assess accuracy of the presented method. The
uncertainty in determined parameters was calculated using the variance propagation rule by Gauss. The presented method is appropriate for membrane water-walls.
The developed meters have one particular advantage over the existing flux tubes to date.The temperature distribution in the flux tube is not affected by the water wall
tubes, since the flux tube is not connected to adjacent waterwall tubes with metal bars, referred to as membrane or webs. To determine the unknown parameters only
the temperature distribution at the cross section of the flux tube must be analysed.
FIGURE 1.
The heat flux tube placed between two water wall tubes, a flux tube, b water wall tube, c thermal insulation
The boundary conditions at the outer and inner surfaces of the water flux-tube must then be determined from temperature measurements at the interior locations. Four
K-type sheathed thermocouples, 1 mm in diameter, are inserted into holes, which are parallel to the tube axis. The thermal conduction effect at the hot junction is
minimized because the thermocouples pass through isothermal holes. The thermocouples are brought to the rear of the tube in the slot machined in the protecting pad.
An austenitic cover plate with the thickness of 3 mm welded to the tube is used to protect the thermocouples from the incident flame radiation. A K-type sheathed
thermocouple witha pad is used to measure the temperature at the rear of the flux-tube. This temperature is almost the same as the water-steam temperature. A method
for determining fireside heat flux, heat transfer coefficient on the inner surface and temperature of water-steam mixture in water-wall tubes is developed. The
unknown parameters are estimated based on the temperature measurements at a few internal locations from the solution of the inverse heat conduction problem. The
non-linear least squares problem is solved numerically using the LevenbergMarquardt method. The diameter of the measuring tube can be larger than the water-wall
tube diameter. The view factor defining the distribution of the heat flux on the measuring tube circumference is determined using exact analytical formulas and
compared with the results obtained numerically using ANSYS software. The method developed can also be used for an assessment of scale deposition on the inner
surfaces of the water wall tubes or slagging on the fire side. The presented method is suitable for water walls made of bare tubes as well as for membrane water walls.
The heat transfer conditions in adjacent boiler tubes have no impact on the temperature distribution in the flux tubes.
2.1. VIEW FACTOR FOR RADIATION HEAT TRANSFER BETWEEN HEAT FLUX TUBE AND FLAME
The heat flux distribution in the flux tube depends heavily on the heat flux distribution on its outer surface. To determine the heat flux distribution q as a function of
angular coordinate , the analytical formulas for the view factor , defining radiation interchange between an infinitesimal surface on the outer flux tube
circumference and the infinite flame or boiler surface, will be derived. The heat flux absorbed by the outer surface of the heat flux tube q() is given by
q()=qm().q()=qm().
(2)
Options
The specific thermal load of the water wall qmis defined as the ratio of the heat transfer rate absorbed by the waterwall to the projected surface area of the water wall.
The view factor is the fraction of the radiation leaving the surface element located on the flux tube surface that arrives at the flame surface. The view factor can be
computed from
=12(sin1+sin2).=12(sin1+sin2).
Options
(3)
The angles 1 and 2 are formed by the normal to the flux tube at and the tangents to the flux tube and adjacent water-wall tube (Figures 2,4,6). Positive values of
1 are measured clockwise with respect to the normal while positive values of 2 are measured counterclockwise with respect to the normal. The radial
coordinate ro of the flux tube outer surface measured from the center 0 (Figure 2) is
ro=ecos+b2e2(sin)2.ro=ecos+b2e2(sin)2.
(4)
Options
(5)
Options
1=arcsin[(ecos+b2(esin)2)sinb],21.1=arcsin[(ecos+b2(esin)2)sinb],21.
(6)
Options
First, the view factor for the angle interval 011,l 1011,l
1 was determined
=1+cos12,01 1,l 1.=1+cos12,01
1,
l 1.
(7)
Options
(8)
Options
(Figure3).
The view factor is computed from Eq.(2), taking into account that (Figure 4)
1=2,2=2(+1),=+2,xi=bsin1,xi=bcos1,=arcsinc(txi)2+(yi+e)2,=arcsintxi(txi)2+(yi+e)2,l 1 1 l 2,1=2,2=2(+
(9
1),=+2,xi=bsin1,xi=bcos1,=arcsinc(txi)2+(yi+e)2,=arcsintxi(txi)2+(yi+e)2,l
1
1
l 2,
Options
FIGURE 2.
Determination of view factor in the angle interval
0 1 1,l10 1
1,l1
FIGURE 3.
Limit angles
1,l11,l1
and
1,l21,l2
FIGURE 4.
Determination of view factor in the angle interval
1,l 1 1 1,l21,l
1
1 1,l2
FIGURE 5.
Limit angle
1,l31,l3
The limit angle1,l31,l3
(Figure 5) can be expressed as
1,l3=2++,1,l3=2++,
(10)
Options
(11)
Options
=arccosb+ct2+e2.=arccosb+ct2+e2.
(12)
Options
(13)
Options
where
1=2,1=2,
Options
(14)
2=+12,2=+12,
(15)
Options
=+2=+2
(16)
Options
=arcsinc(txi)2+(yi+e)2,=arcsinc(txi)2+(yi+e)2,
(17)
Options
=arcsintxi(txi)2+(yi+e)2,=arcsintxi(txi)2+(yi+e)2,
(18)
Options
xi=bsin1,xi=bsin1,
(19)
Options
yi=bcos1.yi=bcos1.
Options
(20)
FIGURE 6.
Determination of view factor in the angle interval
1,l2 11,l31,l2
11,l3
FIGURE 7.
Determination of mean view factor bs for boiler setting over tube pitch t using the crossed string method
Radiation leaving the flame reaches also the boiler setting.The view factor for the radiation heat exchange between boiler setting and rear side of the measuring tube
can be calculated in similar way as for the forward part. The mean heat flux qbsresulting from the radiation heat transfer between the flame and the boiler setting can
be determined using the crossed-string method [20-21].
The mean value of the view factor bs over the pitch length t is calculated from (Figure7)
bs=12t[(FC+BG)(FG+BC)]bs=12t[(FC+BG)(FG+BC)]
(21)
Options
After substituting the lengths of straight FC and BG and circular segments FG and BC into Eq. (21), the mean value of the view factor over the boiler setting can be
expressed as:
bs=b+ct(tan).bs=b+ct(tan).
(22)
Options
qbs=qmbs.qbs=qmbs.
(23)
Options
If the diameters of the heat flux and waterwall tubes are equal, then Eq.(24) simplifies to
(24)
tan=(t2c)21.tan=(t2c)21.
(25)
Options
The view factor for the radiation heat exchange between boiler setting and rear side of the measuring tube can be calculated in similar way as for the forward part.
The view factor in the angle interval 1,l411,l51,l411,l5
(26)
Options
(Figure 8)
1,l4=2+.1,l4=2+.
Options
(27)
FIGURE 8.
Limit angles
1,l41,l4
and
1,l5=1,l2=(/2)+arcsin(e/b)1,l5=1,l2=(/2)+arcsin(e/b)
FIGURE 9.
Determination of view factor in the angle interval
1,l4 11,l51,l4
11,l5
The angles 1and 2 are (Figure 9)
1=2+1,1=2+1,
Options
(28)
2=2,2=2,
(29)
Options
where
=+2,=+2,
(30)
Options
=arcsinc(txi)2+(yi+e)2,=arcsinc(txi)2+(yi+e)2,
(31)
Options
=arcsintxi(txi)2+(yi+e)2,=arcsintxi(txi)2+(yi+e)2,
(32)
Options
xi=bsin1,xi=bsin1,
(33)
Options
yi=bcos1.yi=bcos1.
(34)
Options
where 1,l5=1,l21,l5=1,l2,
is given by
=bs12(sin1+sin2),1,l5 ,=bs12(sin1+sin2),1,
l5 ,
Options
(35)
where
1=12,1=12,
(36)
Options
2=2,2=2,
(37)
Options
=2(),=2(),
(38)
Options
=arcsinc(txi)2+(yi+e)2,=arcsinc(txi)2+(yi+e)2,
(39)
Options
=arcsintxi(txi)2+(yi+e)2,=arcsintxi(txi)2+(yi+e)2,
(40)
Options
xi=bcos(12),xi=bcos(12),
(41)
Options
yi=bsin(12).yi=bsin(12).
Options
(42)
FIGURE 10.
Determination of view factor in the angle interval
1,l5 11,l5
1
The total view factor accounts for the radiation heat exchange between the heat flux tube and flame and between the heat flux tube and the boiler setting.
heat transfer coefficient at the inner surface of the measuring tube does not vary on the tube circumference,
rear side of the water wall, including the measuring tube, is thermally insulated,
diameter of the eccentric flux tube is larger than the diameter of the water wall tubes,
the outside surface of the measuring flux tube is irradiated by the flame, so the heat absorption on the tube fire side is non-uniform.
FIGURE 11.
Approximation of the boundary condition on the outer tube surface
The temperature distribution in the eccentric heat flux tube is governed by heat conduction
1rr(krr)+1r(kr)=01rr(krr)+1r(kr)=0
(43)
Options
kn|r=ro=qm()kn|r=ro=qm()
(44)
krr=a=h|r=akr|r=a=h|r=a
(45)
Options
Options
(46)
Options
The second term in Eq. (46) can be neglected since it is very small and the boundary condition (44) simplifies to
krr=ro=qm()cos(1)kr|r=ro=qm()cos(1)
(47)
Options
The heat flux over the tube circumference can be approximated by the Fourier polynomial
qm()cos(1)=q0+n=1qncos(n)qm()cos(1)=q0+n=1qncos(n)
(48)
Options
where
q0=10qm()cos(1)d,qn=20qm()cos(1)cos(n)d,n=1,...q0=10qm()cos(1)d,qn=20qm()cos(1)cos(n)d,n=1,...
(49)
Options
The boundary value problem (43, 45, 47) was solved using the separation of variables to give
(r,)=A0+B0lnr+n=1(Cnrn+Dnrn)cosn.(r,)=A0+B0lnr+n=1(Cnrn+Dnrn)cosn.
(50)
Options
where
A0=q0ro()k(1Bilna),A0=q0ro()k(1Bilna),
(51)
Options
B0=q0ro()k,B0=q0ro()k,
(52)
Cn=qnro()k1nun(Bi+n)1anBi(u2n+1)+n(u2n1),Cn=qnro()k1nun(Bi+n)1anBi(u2n+1)+n(u2n1),
(53)
Dn=qnro()k1nun(Bin)anBi(u2n+1)+n(u2n1).Dn=qnro()k1nun(Bin)anBi(u2n+1)+n(u2n1).
(54)
Options
Options
Options
The ratio of the outer to inner radius of the eccentric flux tube: u = u()= ro()/a depends on the angle, since the outer radius of the tube flux
ro=ecos+b2(esin)2ro=ecos+b2(e
sin)2
(55)
Options
absorbed heat flux referred to the projected furnace wall surface: x1= qm,
heat transfer coefficient on the inner surface of the boiler tube: x2= h,
These parameters appear in boundary conditions (44) and (45) and will be determined based on the wall temperature measurements at m internal points (ri,i)
T(ri,i)=fi,i=1,...,m,m3.T(ri,i)=fi,i=1,...,m,m3.
(56)
Options
In a general case, the unknown parameters: x1, , xn are determined by minimizing sum of squares
S=(fTm)T(fTm),S=(fTm)T(fTm),
(57)
Options
where f = (f1, , fm)T is the vector of measured temperatures, and Tm = (T1, , Tm)T the vector of computed temperatures Ti = T(ri,i), i = 1, , m.
The parameters x1... xn, for which the sum (34) is minimum are determined using the Levenberg-Marquardt method [23,25]. The parameters, x, are calculated by the
following iteration
x(k+1)=x(k)+(k),k=0,1,....x(k+1)=x(k)+(k),k=0,1,....
(58)
Options
where
(k)=[(J(k)m)TJ(k)m+(k)In]1(J(k)m)T[fTm(x(k))].(k)=[(Jm(k))TJm(k)+(k)In]1(Jm(k))T[fTm(x(k))].
(59)
Options
where (k)(k)is the multiplier and In is the identity matrix. The LevenbergMarquardt method is a combination of the GaussNewton method ((k)0) and the
steepest-descent method ((k)). The mx nJacobian matrix of T(x(k), ri) is given by
J(k)=T(x)xTx=x(k)=T1x1Tmx1T1xnTmxnx=x(k),m=5,
(60)
n=3,J(k)=T(x)xT|x=x(k)=[T1x1T1xnTmx1Tmxn]|x=x(k),m=5,n=3,
Options
The symbol In denotes the identity matrix of n n dimension, and (k) the weight coefficient, which changes in accordance with the algorithm suggested by
Levenberg and Marquardt. The upper index Tdenotes the transposed matrix. Temperature distribution T(r,, x(k)) is computed at each iteration step using Eq. (50).
After a few iterations we obtain a convergent solution.
(61
(xikk)2+j=1m(xirjrj)2+j=1m(xijj)2]1/2,i= 1,2,3
Options
The 95% uncertainty in the estimated parameters can be expressed in the form
xi=xi2xi,xi=xi*2xi,
Options
where
xi,i=1,2,3xi*,i=1,2,3
(62)
represent the value of the parameters obtained using the least squares method. The sensitivity coefficients
xi/fj,xi/k,xi/fj,xi/k,
xi/rjxi/rj
, and
xi/jxi/j
in the expression (61) were calculated by means of the numerical approximation using central difference quotients:
xifj=xi(f1,f2,...,fj+,...,fm)xi(f1,f2,...,fj,...,fm)2,xi
fj=xi(f1,f2,...,fj+,...,fm)xi(f1,f2,...,fj,...,fm)2,
Options
Firstly, a computational example will be presented. Experimental data are generated artificially using the analytical solution (50).
(63)
FIGURE 12.
View factor associated with radiation heat exchange between elemental surface on the boiler setting or flux tube and flame: (a) view factor for radiation heat
transfer between flame and boiler setting, (b) 1 - total view factor accounting radiation from furnace and boiler setting, 2 - approximation by the Fourier polynomial
of the seventh degree, 3 - exact view factor for furnace radiation, 4- view factor from boiler setting
The view factor distributions on the outer surface of the flux-tube and boiler setting were calculated analytically and numerically by means of the finite element
method (FEM) [22]. The changes of the view factor over the pitch length and tube circumference are illustrated in Figures 12 and 13.
FIGURE 13.
Comparison of total view factor calculated by exact and FEM method
The agreement between the temperatures of the outer and inner tube surfaces which were calculated analytically and numerically is also very good (Figures 14
and 15). The small differences between the analytical and FEM solutions are caused by the approximate boundary condition (47). The temperature distribution in the
flux tube cross section is shown in Figure 14.
FIGURE 14.
Computed temperature distribution in oC in the cross section of the heat flux tube; qm = 200000 W/m2, h = 30000 W/(m2K), Tf=318 oC
FIGURE 15.
Temperature distribution at the inner and outer surfaces of the flux tube calculated by the analytical and finite element method
The following input data is generated using Eq. (50):
2fj= 0.2K ,j= 1, ,5,2k= 0.5 W/(mK),2rj=0.05mm,2j=0.5o,j=1,,5.2fj= 0.2K ,j= 1, ,5,2k= 0.5 W/
(mK),2rj=0.05mm,2j=0.5o,j=1,,5.
The uncertainties (95% confidence interval) of the coefficients xi were determined using the error propagation rule formulated by Gauss.
The calculation using Eq. (61) yielded the following results: x1 =
30 001.5630 001.56
2698.81 W/(m2 K), x3 = 318.0 0.11 oC. The accuracy of the obtained results is very satisfactory.There is only a small difference between the estimated parameters and the input values. The
highest temperature occurs at the crown of the flux-tube ( Figures 14 and 15). The temperature of the inner surface of the flux tube is only a few degrees above the saturation temperature of the
water-steam mixture. Since the heat flux at the rear side of the tube is small, the circumferential heat flow rate is significant. However, the rear surface thermocouple indicates temperatures of 24 oC above the saturation temperature. Therefore, the fifth thermocouple can be attached to the unheated side of the tube so as to measure the temperature of the water-steam mixture flowing
through the flux tube.
In the second example, experimental results will be presented. Measurements were conducted at a 50MW pulverized coal fired boiler. The temperatures indicated by
the flux tube at the elevation of 19.2 m are shown in Figure 16. The heat flux tube is of 20G low carbon steel with temperature dependent thermal conductivity
k(T)=53.260.02376224T,k(T)=53.260.02376224T,
(64)
Options
FIGURE 16.
Measured flux tube temperatures; marks denote measured temperatures taken for the inverse analysis
FIGURE 17.
Estimated parameters: absorbed heat flux qm, heat transfer coefficient h, and temperature of water-steam mixture Tf
Flux tubes can also be used as a local slag monitor to detect a build up of slag. The presence of the scale on the inner surface of the tube wall can also be detected.
FIGURE 18.
The cross-section of the membrane wall in the combustion chamber of the steam boiler
Test calculations were carried out to assess accuracy of the presented method. The uncertainty in determined parameters was calculated using the Gauss variance
propagation rule. The presented method is appropriate for membrane water walls (Figure 18). The new method has advantages in terms of simplicity and flexibility.
3.1. THEORY
The furnace wall tubes in most modern units are welded together with steel bars (fins) to provide membrane wall panels which are insulated on one side and exposed
to a furnace on the other, as shown schematically in Figure 18.
In a heat conduction model of the flux-tube the following assumptions are made:
the heat transfer coefficient hin and the scale thickness ds is uniform over the inner tube surface.
[k(T)T]=0,[k(T)T]=0,
Options
where is the vector operator, which is called nabla (gradient operator), and in Cartesian coordinates is defined by = i/x +j/y + k/z +. The unknown
boundary conditions may be expressed as
(65)
[k(T)Tn]s=q(s),[k(T)Tn]|s=q(s),
(66)
Options
where q(s) is the radiation heat flux absorbed by the exposed fluxtube and membrane wall surface. The local heat flux q(s) is a function of the view factor (s)
(Figure 19)
q(s)=qm(s),q(s)=qm(s),
(67)
Options
where qm is measured heat flux (thermal loading of heating surface). The view factor (s) from the infinite flame plane to the differential element on the membrane
wall surface can be determined graphically[7], or numerically [22].
In this chapter, (s) was evaluated numerically using the finite element program ANSYS [22], and is displayed in Figure 19 as a function of the extended
coordinate s.Because of the symmetry, only the representative water-wall section illustrated in Figure 20 needs to be analyzed. The convective heat transfer from the
inside tube surfaces to the water-steam mixture is described by Newtons law of cooling
[k(T)Tn]sin=hin(T|sinTf),[k(T)Tn]|sin=hin(T|sinTf),
(68)
Options
where T/n is the derivative in the normal direction, hin is the heat transfer coefficient and Tf denotes the temperature of the watersteam mixture.
The reverse side of the membrane water-wall is thermally insulated.In addition to the unknown boundary conditions, the internal temperature measurements fiare
included in the analysis
Te(ri)=fi,i=1,,m,Te(ri)=fi,i=1,,m,
(69)
Options
where m = 5 denotes the number of thermocouples (Figure 18). The unknown parameters: x1 = qm, x2= hin, and x3 = Tf were determined using the least-squares
method. The symbol rin denotes the inside tube radius, and k(T) is the temperature dependent thermal conductivity. The object is to choose x = (x1, , xn)T for n = 3
such that computed temperatures T(x, ri) agree within certain limits with the experimentally measured temperatures fi.
This may be expressed as
T(x,ri)fi0,i=1,,m,m=5.T(x,ri)fi0,i=1,,m,m=5.
Options
(70)
FIGURE 19.
View factor distribution on the outer surface of water-wall tube
FIGURE 20.
Temperature distribution in the flux tube cross-section for: qm = 150000 W/m2, hin = 27000 W/(m2K) and Tf = 317C
The least-squares method is used to determine parameters x. The sum of squares
S=i=1m[fiT(x,ri)]2,m=5,S=i=1m[fiT(x,ri)]2,m=5,
(71)
Options
The material of the heat flux-tube is 20G steel. The composition of the 20G mild steel is as follows: 0.170.24% C, 0.71.0% Mn, 0.150.40% Si, Max 0.04% P, Max
0.04% S, and the remainder is iron Fe. The heat flux-tube thermal conductivity is assumed to be temperature dependent (Table 1).
FIGURE 21.
Solution of the inverse problem for the exact data:f1 = 419.66C, f2 = 417.31C, f3 = 374.90C, f4= 373.19C, f5 = 318.01C ; (a) - temperature distribution in the
flux-tube, (b) - iteration number for the temperature T1
FIGURE 22.
Solution of the inverse problem for the perturbed data: f1 = 420.16C, f2 = 416.81C, f3 = 375.40C, f4 = 372.69C, f5 = 318.01C;(a) - temperature distribution in the
flux-tube, (b) - iteration number for the temperature T1
TemperatureT, C
Thermal conductivityk, W/(mK)
TABLE 1.
Thermal conductivity k(T) of steel 20G as a function of temperature
To demonstrate that the maximum temperature of the fin tip is lower than the allowable temperature for the 20G steel, the flux tube temperature was computed using
ANSYS CFX package [22]. Changes of the view factor on the flux tube, weld and fin surface were calculated with ANSYS CFX. The temperature distribution shown
in Figure 20 was obtained for the following data: absorbed heat flux,qm = 150000 W/m2, temperature of the water-steam mixture, Tf = 317C, and heat transfer
coefficient at the tube inner surface, hin = 27000 W/(m2 K). An inspection of the results shown in Figure 20indicates that the maximum temperature of the fin does
not exceed 375C.
Next, to illustrate the effectiveness of the presented method, test calculations were carried out. The measured temperatures fi, i = 1, 2, , 5 were generated
artificially by means of ANSYS CFX for:qm = 250000 W/m2, hin = 30000 W/(m2 K) and Tf = 318C. The following values of measured temperatures were
obtained f1 = 419.66C, f2 = 417.31C, f3 = 374.90C, f4 = 373.19C, f5 = 318.01C. The temperature distribution in the flux tube cross-section, reconstructed on the
basis of five measured temperatures is depicted in Figure 21a.
The proposed inverse method is very accurate since the estimated parameters: qm = 250000.063 W/m2,hin = 30000.054 W/(m2 K)and Tf = 318.0C differ
insignificantly from the input values. In order to show the influence of the measurement errors on the determined parameters, the 95% confidence intervals were
estimated. The following uncertainties of the measured values were assumed (at 95% confidence interval):
2fj=0.5 K,2fj=0.5 K,
j = 1, 2, , 5,
2rj2rj
=0.05mm,
2j2j
=0.5o,j=1,,5. The uncertainties (95% confidence interval) of the coefficients xi were determined using the error propagation rule formulated by Gauss [ 23-26]. The calculated uncertainties
are: 6% for qm, 33% for hin and 0.3% for Tf. The accuracy of the results obtained is acceptable.
Then, the inverse analysis was carried out for perturbed data: f1 = 420.16C, f2 = 416.81C, f3 = 375.40C, f4 = 372.69C, f5 = 318.01C. The reconstructed temperature
distribution illustrates Figure 22a.
The obtained results are: qm = 250118.613 W/m2, hin = 30050.041 W/(m2 K) and Tf = 317.99C. The errors in the measured temperatures have little effect on the
estimated parameters. The number of iterations in the Levenberg-Marquardt procedure is small in both cases (Figures 21b and 22b).
4. Conclusions
Two different tubular type instruments (flux tubes) were developed to identify boundary conditions in water wall tubes of steam boilers.The first measuring device is
an eccentric tube. The ends of the four thermocouples are located at the fireside part of the tube and the fifth thermocouple is attached to the unheated rear surface of
the tube. The meter presented in the paper has one particular advantage over the existing flux tubes to date. The temperature distribution in the flux tube is not
affected by the water wall tubes, since the flux tube is not connected to adjacent waterwall tubes with metal bars, referred to as membrane or webs. To determine the
unknown parameters only the temperature distribution at the cross section of the flux tube must be analyzed.
The second flux tube has two longitudinal fins. Fins attached to the flux tube are not welded to the adjacent water-wall tubes, so the temperature distribution in the
measuring device is not affected by neighboring water-wall tubes. The installation of the flux tube is easier because welding of fins to adjacent water-wall tubes is
avoided. Based on the measured flux tube temperatures the non-linear inverse heat conduction problem was solved. A CFD based method for determining heat flux
absorbed water wall tubes, heat transfer coefficient at the inner flux tube surface and temperature of the water-steam mixture has been presented. The proposed flux
tube and the inverse procedure for determining absorbed heat flux can be used both when the inner surface of the heat flux tube is clean and when scale or corrosion
deposits are present on the inner surface what can occur after a long time service of the heat flux tube.
The flux tubes can work for a long time in the destructive high temperature atmosphere of a coal-fired boiler.
NOMENCLATURE
a inner radius of boiler tube and flux-tube (m)
b outer radius of flux-tube (m)
Bi Biot number, Bi =ha/k
c outer radius of boiler tube (m)
e eccentric (m)
fi measured wall temperature at the i-th location (oC or K)
f vector of measured wall temperatures
h heat transfer coefficient (W/(m2K))
In identity matrix
J Jacobian matrix of T
k thermal conductivity (W/(mK))
l arbitrary length of boiler tube (m)
m number of temperature measurement points
n number of unknown parameters
qm heat flux to be determined (absorbed heat flux referred to the projected furnace water wall surface) (W/m 2)
r coordinate in cylindrical coordinate system or radius (m)
ri radial coordinate of the i-th thermocouple (m)
rin inner radius of the flux-tube ( m)
ro outer radius of the flux-tube ( m)
r position vector
s extended coordinate along the fireside water-wall surface (m)
S sum of the temperature difference squares (K2)
t pitch of the water wall tubes (m)
T temperature (oCor K)
Tf fluid temperature(C or K)
Ti calculated temperature at the location(ri,i) (C or K)
Tm m - dimensional column vector of calculated temperatures
u() ratio of the outer to the inner radius of the tube, u()= ro/a
xi unknown parameter
GREEK SYMBOLS
, , , 1, 2, angles (rad)
temperature excess over the fluid temperature,= T - Tf (K)
angular coordinate (rad)
i angular coordinate of the i-th thermocouple (rad)
multiplier in the Levenberg-Marquardt algorithm
view factor
SUBSCRIPTS
in inner
o outer
i i-thtemperature measurement point
f fluid