Measurements of Local Heat Flux and Water-Side Heat Transfer Coefficient in Water Wall Tubes

Download as docx, pdf, or txt
Download as docx, pdf, or txt
You are on page 1of 61

Measurements of Local Heat Flux and Water-Side Heat Transfer Coefficient in Water

Wall Tubes
Jan Taler1 and Dawid Taler2
[1]

Department of Thermal Power Engineering, Cracow University of Technology, Cracow, Poland

[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.

2. Tubular type heat flux meter made of a bare tube


Heat flux meters are used for monitoring local waterwall slagging in coal and biomass fired steam boilers [5-19].
The tubular type instruments (flux tubes) [10-14,19] and other measuring devices [15-18] were developed to identify boundary conditions in water wall tubes of
steam boilers. The meter is constructed from a short length of eccentric 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.

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

where: e eccentric (Figure 2), b outer radius of flux-tube.


The angle11can be expressed in terms of the angle , flux tube outer radius b, and eccentric e (Figure 2)
1=arcsin[(ecos+b2(esin)2)sinb],12,1=arcsin[(ecos+b2(esin)2)sinb],12,

(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

The limit angle 1,l 11,


l 1(Figures 2 and 3) is given by
1,l 1=arccosceb,1,l
1=arccosceb,

(8)

Options

where cis the outer radius of the boiler tube.


Next the view factor in the angle interval 1,l 111,l 21,
l 111,

l 2will be determined. The limit angle1,l 21,


l 2is:1,l2=
1(=/2)=(/2)+arcsin(e/b)1,l2=1(=/2)=(/2)+arcsin(e/b)

(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

where t is the pitch of the water wall tubes.


Next the view factor () is determined in the angle interval 1,l211,l31,l211,l3
(Figures 3 and 5).

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

where the angles i are given by


=arctanbct,=arctanbct,

(11)

Options

=arccosb+ct2+e2.=arccosb+ct2+e2.

(12)

Options

The view factor in the interval 1,l211,l 31,


l 211,l

3 is calculated from the following expression (Figure 6)
=12(sin2sin1),1,l211,l 3,=12(sin2sin1),1,l
211,l

3,

(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

The mean heat flux over the setting surface is

qbs=qmbs.qbs=qmbs.

(23)

Options

The angle is determined from


tan=e2+t2(b+c)2b+c,tan=e2+t2(b+c)2b+c,
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

(Figure 8), accounting for the setting radiation, is given by


=bs12(sin2sin1),1,l411,l5=bs12(sin2sin1),1,l411,l5

(26)

Options

where the limit angle 1,l41,l4is

(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

The view factor in the interval 1,l51,


l5,

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.

2.2. THEORY OF THE INVERSE PROBLEM


At first, the temperature distribution at the cross section of the measuring tube will be determined, i.e. the direct problem will be solved. Linear direct heat conduction
problem can be solved using an analytical method. The temperature distribution will also be calculated numerically using the finite element method (FEM). In order
to show accuracy of a numerical approach, the results obtained from numerical and analytical methods will be compared. The following assumptions have been
made:

thermal conductivity of the flux tube material is constant,

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.

The cylindrical coordinate system is shown in Figure11.

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

subject to the following boundary conditions

kn|r=ro=qm()kn|r=ro=qm()

(44)

krr=a=h|r=akr|r=a=h|r=a

(45)

Options

Options

The left side of Eq. (44) can be transformed as follows (Figure11)


kn|r=ro=(qr+q)nr=ro==[kTrcos(1)+krTsin(1)]r=rokn|r=ro=(qr+q)n|r=ro==[kTrcos(1)+krTsin(1)]|r=ro

(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

is the function of the angle .


Eq. (50) can be used for the temperature calculation when all the boundary conditions are known. In the inverse heat conduction problem three parameters are to be
determined:

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,

fluid bulk temperature: x3=Tf.

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.

2.3. THE UNCERTAINTY OF THE RESULTS


The uncertainties of the determined parameters x* will be estimated using the error propagation rule of Gauss [23-26]. The propagation of uncertainty in the
independent variables: measured wall temperatures fj, j=1, m, thermal conductivity k, radial and angular positions of temperature sensorsrj, j, j=1, m is
estimated from the following equation
2xi=[j=1m(xifjfj)2+(xikk)2+j=1m(xirjrj)2+j=1m(xijj)2]1/2,i= 1,2,32xi=[j=1m(xifjfj)2+

(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

where is a small number.

2.4.COMPUTATIONAL. AND BOILER TESTS

Firstly, a computational example will be presented. Experimental data are generated artificially using the analytical solution (50).

Consider a water-wall tube with the following parameters (Figure1.):

outer radius b = 35 mm,

inner radius a = 25 mm,

pitch of the water-wall tubes t = 80 mm,

thermal conductivity k = 28.5 W/(m K),

(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

absorbed heat flux qm = 200000 W/m2,

heat transfer coefficient h = 30000 W/(m2 K),

fluid temperature Tf= 318 oC.

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):

f1=437.98 oC, f2=434.47 oC,f1=437.98 oC, f2=434.47 oC,


f3=383.35 oC, f4=380.70 oC, f5=321.58oC.f3=383.35 oC, f4=380.70 oC, f5=321.58oC.
The following values were obtained using the proposed method:

qm=200000.35 W/m2, h=30001.56 W/(m2K),Tf=318.00 oC.qm*=200000.35

W/m2, h*=30001.56 W/(m2K),Tf*=318.00 oC.


In order to show the influence of the measurement errors on the determined thermal boundary parameters, the 95% confidence intervals were calculated. The
following uncertainties of the measured values were assumed (at a 95% confidence interval):

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 =

200 000.35200 000.35


3827.72 W/m2, x2 =

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

where the temperature T is expressed in oC and thermal conductivity in W/(m K).


The unknown parameters were determined for eight time points which are marked in Figure16.
The inverse analysis was performed assuming the constant thermal conductivity k(T)k(T) which was obtained from Eq.(64) for the average
temperature: T=(T1+T2+T3+T4)/4T=(T1+T2+T3+T4)/4.
The estimated parameters: heat flux qm, heat transfer coefficient h, and the water-steam mixture Tf are depicted in Figure17.The developed flux tube can work for a
long time in the destructive high temperature atmosphere of a coal-fired boiler.

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.

3. Tubular type heat flux meter made of a finned tube


In this section, a numerical method for determining the heat flux in boiler furnaces, based on experimentally acquired interior flux-tube temperatures, is presented.
The tubular type instrument has been designed (Figure 18) to provide a very accurate measurement of absorbed heat flux qm, inside heat transfer coefficient hin, and
water steam temperature Tf. The number of thermocouples is greater than three because the additional information can help enhance the accuracy of parameter
determining. In contrast to the 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. 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 in 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. An inverse problem of heat conduction was solved using
the least squares method. Three unknown parameters were estimated using the Levenberg-Marquardt method [23, 25]. At every iteration step, the temperature
distribution over the cross-section of the heat flux meter was computed using the ANSYS CFX software

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:

temperature distribution is two-dimensional and steady-state,

the thermal conductivity of the flux-tube and membrane wall,

may be dependent of temperature,

the heat transfer coefficient hin and the scale thickness ds is uniform over the inner tube surface.

The temperature distribution is governed by the non-linear partial differential equation

[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

is minimized usingthe LevenbergMarquardt method [23, 25].


The uncertainties of the determined parameters x* will be estimated using the error propagation rule of Gauss [23-26].

3.2. TEST COMPUTATIONS


The flux-tubes were manufactured in the laboratory and then securely welded to the water-wall tubes at different elevations in the furnace of the steam boiler. The
coal fired boiler produces 58.3 kg/s superheated steam at 11 MPa and 540C.

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,

2k=1 W/(mK)2k=1 W/(mK)


,

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

x n-dimensional column vector of unknown parameters

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

You might also like