Thermal conductivity *k*(*T*) of steel 20G as a function of temperature

## 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 flux *q*_{m}. The measuring tube is fitted with two thermocouples in holes of known radial spacing *r*_{1} and *r*_{2}. 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

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 *q*_{m}, inside heat transfer coefficient *h*_{in}, and water steam temperature *T*_{f}.

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 Levenberg–Marquardt 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 cross-section 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.

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 Levenberg–Marquardt 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

The specific thermal load of the water wall *q*_{m}is 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

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 *r*_{o} of the flux tube outer surface measured from the center 0 (Figure 2) is

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

The angle*φ*, flux tube outer radius *b*, and eccentric *e* (Figure 2)

First, the view factor for the angle interval

The limit angle

where *c*is the outer radius of the boiler tube.

Next the view factor in the angle interval

where *t* is the pitch of the water wall tubes.

Next the view factor *ψ* (*φ*) is determined in the angle interval

The limit angle

where the angles *κ* i *ω* are given by

The view factor *ψ* in the interval

where

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 *q*_{bs}resulting 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)

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:

The mean heat flux over the setting surface is

The angle *ω*is determined from

If the diameters of the heat flux and waterwall tubes are equal, then Eq.(24) simplifies to

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

where the limit angle

The angles *δ*_{1}and *δ*_{2} are (Figure 9)

where

The view factor *ψ* in the interval

where

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.

The temperature distribution in the eccentric heat flux tube is governed by heat conduction

subject to the following boundary conditions

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

The second term in Eq. (46) can be neglected since it is very small and the boundary condition (44) simplifies to

The heat flux over the tube circumference can be approximated by the Fourier polynomial

where

The boundary value problem (43, 45, 47) was solved using the separation of variables to give

where

The ratio of the outer to inner radius of the eccentric flux tube: *u* = *u*(*φ*)= *r*_{o}(*φ*)/*a* depends on the angle *φ*_{,} since the outer radius of the tube flux

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:

*x*_{1}=*q*_{m},heat transfer coefficient on the inner surface of the boiler tube:

*x*_{2}*= h*,fluid bulk temperature:

*x*_{3}*=T*_{f}.

These parameters appear in boundary conditions (44) and (45) and will be determined based on the wall temperature measurements at *m* internal points (*r*_{i},*φ*_{i})

In a general case, the unknown parameters: *x*_{1}, …, *x*_{n} are determined by minimizing sum of squares

where **f** = (*f*_{1}, …, *f*_{m})^{T} is the vector of measured temperatures, and **T**_{m} = (*T*_{1}, …, *T*_{m})^{T} the vector of computed temperatures *T*_{i} = *T*(*r*_{i},*φ*_{i}), *i* = 1, …, *m*.

The parameters *x*_{1}*... x*_{n}, 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

where

where **I**_{n} is the identity matrix. The Levenberg–Marquardt method is a combination of the Gauss–Newton method (*μ*^{(k)}→0) and the steepest-descent method (*μ*^{(k)}→∞). The *m* x *n*Jacobian matrix of *T*(**x**^{(k)}, **r**_{i}) is given by

The symbol **I**_{n} 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 *T* denotes 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 *f*_{j}, *j*=1, …*m*, thermal conductivity *k*, radial and angular positions of temperature sensors *r*_{j}, *φ*_{j}, *j*=1, …*m* is estimated from the following equation

The 95% uncertainty in the estimated parameters can be expressed in the form

where

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

absorbed heat flux

*q*_{m}= 200000 W/m^{2},heat transfer coefficient

*h*= 30000 W/(m^{2}K),fluid temperature

*T*_{f}= 318^{o}C.

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.

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.

The following input data is generated using Eq. (50):

The following values were obtained using the proposed method:

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

The uncertainties (95% confidence interval) of the coefficients *x*_{i} were determined using the error propagation rule formulated by Gauss.

The calculation using Eq. (61) yielded the following results: *x*_{1} = ^{2}, *x*_{2} = ^{2} K), *x*_{3} = 318.0 ± 0.11 ^{o}C. 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 2-4 ^{o}C 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

where the temperature *T* is expressed in ^{o}C 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

The estimated parameters: heat flux *q*_{m}, heat transfer coefficient *h*, and the water-steam mixture *T*_{f} 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.

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 *q*_{m}, inside heat transfer coefficient *h*_{in}, and water steam temperature *T*_{f}. 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

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

*h*_{in}and the scale thickness ds is uniform over the inner tube surface.

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

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

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)

where *q*_{m} 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 Newton’s law of cooling

where ∂*T*/∂*n* is the derivative in the normal direction, *h*_{in} is the heat transfer coefficient and *T*_{f} denotes the temperature of the water–steam mixture.

The reverse side of the membrane water-wall is thermally insulated.In addition to the unknown boundary conditions, the internal temperature measurements *f*_{i}are included in the analysis

where *m* = 5 denotes the number of thermocouples (Figure 18). The unknown parameters: *x*_{1} = *q*_{m}, *x*_{2} = *h*_{in}, and *x*_{3} = *T*_{f} were determined using the least-squares method. The symbol *r*_{in} denotes the inside tube radius, and *k*(*T*) is the temperature dependent thermal conductivity. The object is to choose **x** = (*x*_{1}, …, *x*_{n})^{T} for *n* = 3 such that computed temperatures *T*(**x, r**_{i}) agree within certain limits with the experimentally measured temperatures *f*_{i}.

This may be expressed as

The least-squares method is used to determine parameters **x**. The sum of squares

is minimized usingthe Levenberg–Marquardt 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 540ºC.

The material of the heat flux-tube is 20G steel. The composition of the 20G mild steel is as follows: 0.17–0.24% C, 0.7–1.0% Mn, 0.15–0.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).

TemperatureT, ºC | 100 | 200 | 300 | 400 |

Thermal conductivityk, W/(m·K) | 50.69 | 48.60 | 46.09 | 42.30 |

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, *q*_{m} = 150000 W/m^{2}, temperature of the water-steam mixture, *T*_{f} = 317ºC, and heat transfer coefficient at the tube inner surface, *h*_{in} = 27000 W/(m^{2} K). An inspection of the results shown in Figure 20 indicates that the maximum temperature of the fin does not exceed 375ºC.

Next, to illustrate the effectiveness of the presented method, test calculations were carried out. The “measured” temperatures *f*_{i}, *i* = 1, 2, …, 5 were generated artificially by means of ANSYS CFX for: *q*_{m} = 250000 W/m^{2}, *h*_{in} = 30000 W/(m^{2} K) and *T*_{f} = 318ºC. The following values of “measured” temperatures were obtained *f*_{1} = 419.66ºC, *f*_{2} = 417.31ºC, *f*_{3} = 374.90ºC, *f*_{4} = 373.19ºC, *f*_{5} = 318.01ºC. 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: *q*_{m} = 250000.063 W/m^{2}, *h*_{in} = 30000.054 W/(m^{2} K)and *T*_{f} = 318.0°C 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): *j* = 1, 2, …, 5, ^{o},*j*=1,…,5. The uncertainties (95% confidence interval) of the coefficients *x*_{i} were determined using the error propagation rule formulated by Gauss [23-26]. The calculated uncertainties are: ±6% for *q*_{m}, ±33% for *h*_{in} and ±0.3% for *T*_{f}. The accuracy of the results obtained is acceptable.

Then, the inverse analysis was carried out for perturbed data: *f*_{1} = 420.16ºC, *f*_{2} = 416.81ºC, *f*_{3} = 375.40ºC, *f*_{4} = 372.69ºC, *f*_{5} = 318.01ºC. The reconstructed temperature distribution illustrates Figure 22a.

The obtained results are: *q*_{m} = 250118.613 W/m^{2}, *h*_{in} = 30050.041 W/(m^{2} K) and *T*_{f} = 317.99°C. 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/kc outer radius of boiler tube (m)e eccentric (m)f_{i }measured wall temperature at the i-th location (^{o}C or K)f vector of measured wall temperaturesh heat transfer coefficient (W/(m^{2}·K))I _{n }identity matrixJ Jacobian matrix of T k thermal conductivity (W/(m·K))l arbitrary length of boiler tube (m)m number of temperature measurement pointsn number of unknown parametersqm 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)r_{i }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 (K^{2})t pitch of the water wall tubes (m)T temperature (^{o}Cor K)Tf fluid temperature(ºC or K)Ti calculated temperature at the location(ri,φi) (ºC or K)T m m - dimensional column vector of calculated temperaturesu(φ) ratio of the outer to the inner radius of the tube, u(φ)= ro/axi unknown parameterx n-dimensional column vector of unknown parameters |

### Greek symbols

α, β, γ, δ_{1}, δ_{2}, ε angles (rad)θ temperature excess over the fluid temperature,θ= T - T_{f} (K)φ angular coordinate (rad)φi angular coordinate of the i-th thermocouple (rad)μ multiplier in the Levenberg-Marquardt algorithm ψ view factor |

### Subscripts

in innero outeri i-thtemperature measurement pointf fluid |