1. Introduction
Steam superheaters in subcritical and supercritical boilers are used to increase the efficiency of electricity generation in steam power plants. The calculation and testing of superheaters in steam boilers receive significant attention. The interest in superheaters is because of the high temperature of superheater tubes and their high failure rate. More than 40% of emergency boiler shutdowns result from superheater tube failures. The ash fouling of superheater tubes causes severe difficulties in operating steam boilers. In addition to significantly reducing the heat flow rates exchanged between the flue gas and steam, ash deposits cause superheater tubes’ corrosion. For this reason, many articles in the current literature address hydraulic, aerodynamic, thermal, and strength calculations. Examples of steam superheater failures and their causes are discussed extensively by French [
1]. Pronobis [
2] detailed superheater tubes’ erosion and ash fouling processes. The book [
2] contains rich experimental data on ash erosion of superheater tubes from the tests carried out on existing facilities and laboratory stands.
The literature on mathematical modelling and automatic control of supercritical power plants is extensive. The high interest in supercritical blocks is attributed to the large capacities of the power plants currently under construction and the greater efficiency of supercritical units. A very good review of the literature on the modelling and control of supercritical and ultra-supercritical steam power plants is presented by Omar [
3]. Publications presenting mathematical models of supercritical fossil-fuel-fired power plants based on thermal-flow process modelling are discussed. Models based on the laws of physics are more suitable for modelling significant power plant load changes and emergency conditions than the empirical-data-based models developed without any physical considerations. The authors highlight that the start-up optimisation process of a supercritical power plant has not been widely studied and appears to be an upcoming research proposal [
3].
A mathematical model of a supercritical boiler for simulating its transient operation was presented by K. Hedrick et al. in [
4]. Boiler load change from 100% to 60% was simulated, considering the change of boiler operation from subcritical to supercritical pressure. The boiler model was validated using measurements taken at an existing plant facility.
Haddad and Mohamed [
5] presented three modelling procedures for a once-through supercritical power unit. The first mathematical model of a 600 MW power unit is based on the conservation equations of mass, momentum, and energy, the second uses a linearised state–space model, and the third uses artificial neural networks. All three approaches were compared against each other.
Young et al. [
6] studied heat transfer at supercritical pressure with water flowing in a pipe. The pipe was heated at the outer surface. The heat flux varied from 100 to 1200 kW/m
2. The water flow regime varied from subcritical to supercritical. Heat transfer correlations were proposed to calculate the heat transfer coefficient at the pipe’s inner surface.
The method of determining the tube temperature and the thickness of the inner oxide scales at the superheater tubes was developed by Sun and Yan [
7]. The thickness of the oxide scale layer can be determined along the tube’s length as a function of time using the method proposed in [
7]. Local overheating of the tube walls in superheaters is a cause of tube creep and tube failure. Li et al. [
8] developed an online procedure for superheater tube monitoring based on Computational Fluid Dynamics (CFD) modelling and flue gas side measurements. The purpose of the system is to prevent overheating of the tube material. Oxide deposits on the inner surfaces of superheater tubes pose a significant threat to the durability of superheaters. Due to the increase in the operating temperature of the tube walls, there is a danger of the tube material overheating. In addition, stripping of the oxide layer from the pipe surfaces poses a severe threat to the life of the steam turbine. The impact of temperature on the oxide spallation intensity in ferritic–martensitic superheater tubes was studied by Qi et al. in [
9].
Analyses made by Balint et al. [
10] provide new insights into the melting and ageing behaviour of ash deposits, which play an essential role in the corrosion of superheater tubes. Cross sections of ash deposit layers on superheater tubes were analysed using scanning electron microscopy and energy-dispersive X-ray analysis. Failure analysis of the superheater tubes was carried out by Liu [
11] in a boiler fired with black liquor. Creep and thermal fatigue were identified as the main cause of superheater failure. The fractographic analysis identified two cracks at the heat-affected zone of the weld. The cracks initiated on the fireside surface of the tube and moved toward the inner surface of the tube. One of the analysed cracks propagated through the tube wall and caused damaging the superheater tube and steam flow out of it. An important phenomenon that significantly affects the life of superheater tubes is the uneven temperature distribution across the width of superheater platens. The non-uniformity of the mass flow rates in the individual tubes causes temperature differences in parallel tubes of the platens. In addition, the differences in the flue gas temperature in the duct cross-section can cause temperature differences across the width of the superheater. The steam temperature non-uniformity in an L-shaped platen superheater in a supercritical 600 MW circulating fluidised bed boiler (CFB) was studied by Ma et al. [
12]. Thermal flow calculations of the entire steam boiler with natural circulation are presented by Hossain et al. in [
13]. The paper [
13] is particularly interesting to boiler designers because it provides thermal calculations of the entire boiler. The effectiveness–number of transfer units (ε-NTU) method was used for superheater performance calculations, and the logarithmic mean temperature difference (LMTD) method for design calculations. Determining the superheater tube temperatures was not at the aim of this work. Steam superheaters in pulverised and fluidised bed boilers often fail due to the excessive temperature of the tube material. The computational determination of the maximum wall temperature of superheater tubes is very difficult due to the complexity of thermal and flow processes occurring in the boiler. Often, the only reliable method of determining the cause of superheater tube failure is to measure continuously the tube wall temperature during boiler operation. A new technique for measuring the tube wall temperature in CFB boilers was developed by Yao et al. [
14].
CFD methods are widely used in the analysis of flow and thermal phenomena occurring in pulverised coal-fired steam boilers or fluidised bed boilers. A three-dimensional (3D) mathematical model of combustion and heat transfer on the flue gas side and a one-dimensional mathematical model of heat transfer in the tubes of a hanging radiant superheater were developed by Laubscher and Rousseau [
15]. The steady-state and transient flow and thermal phenomena in subcritical boiler steam superheaters were analysed by Granda et al. [
16] using CFD modelling. Thermal and flow calculation results of the second-stage steam superheater located at the combustion chamber of the large-scale CFB boiler are presented by Madejski et al. [
17].
A significant group of calculation methods for steam superheaters are analytical and semi-analytical methods. Usually, a method using the LMTD between the flue gas and steam is used. The LMTD procedure is applied, for example, by Lokshin et al. [
18] in the standard method of hydraulic design for power boilers. In addition, Kagan et al. [
19] utilised the LMTD method for the thermal calculating of superheaters. The LMTD method makes it possible to calculate the heat exchange area of the superheater without considering the temperature of the tubes. Unlike low-temperature heat exchangers, a significant difficulty in calculating steam superheaters is the determination of the radiation heat transfer coefficient (HTC) on the flue gas side. An alternative method to [
19] for calculating the radiation and convection HTCs on the flue gas side was proposed by Olenev [
20]. A simple formula for calculating the radiation HTC at the superheater tube surface based on the P1 method of radiation heat transfer analysis in gases was presented by Taler and Taler [
21].
In large steam boilers with supercritical steam pressure, the wall thickness of superheater tubes is greater than in subcritical boilers. In addition, the steam temperature in modern steam superheaters is approximately 600 °C. In order not to increase excessively the thickness of tube walls and to ensure their resistance to high-temperature corrosion, high-alloy steels are used for superheater tubes. For steam temperatures not exceeding 600 °C, ferritic–martensitic steels are used, and at higher temperatures, austenitic steels are applied. Each superheater pass is made of a different steel grade adjusted to the operating temperature to reduce investment costs. It is necessary to know the temperature variation of steam along the flow path to select the suitable grade of steel for a given superheater pass. A hybrid numerical–analytical method for calculating cross-flow tube heat exchangers was presented by Taler et al. [
22]. The method also makes it possible to calculate quickly the temperature distribution of flue gas, steam, and tube walls, taking into account the temperature-dependent thermal properties of both fluids and the tube material. Flue gas and steam temperature are expressed in each finite volume by exact analytical formulas. A non-iterative calculation procedure for tube cross-flow heat exchangers was also developed by Węglarz et al. [
23]. The temperature of steam and flue gas at the outlet of the finite volume is calculated using analytical expressions. The calculation formulas are straightforward. Computer calculation times are short even with the more significant number of finite volumes required in this method compared with the method presented in the work [
22]. Both semi-analytical methods [
22,
23] are suitable for calculating cross-co-current and cross-counter-current superheaters with complex flow arrangements.
Most of the analytical methods used to calculate tubular heat exchangers, i.e., the LMTD, ε-NTU, and P-NTU methods, assume a constant specific heat throughout the heat exchanger. It should be added here that the ε-NTU method is also used in commercial CFD programs to calculate the average temperature of both fluids at the outlet of tubular heat exchangers, including steam superheaters.
A review of the work published to date shows that there is a lack of uncomplicated yet accurate superheater models to determine local steam, flue gas, and tube wall temperatures. The aim of this article is to develop a new numerical model of a four-pass cross-counter flow live steam superheater. The calculation procedure accounts for temperature-dependent steam, flue gas, and tube material properties. Superheater calculations were carried out for a third-stage live steam superheater installed in a supercritical boiler. The mass capacity of the boiler was 645 kg/s. The results of calculations using the proposed model were validated against the measured temperatures of steam and flue gas at the superheater outlet. The developed method is simple and can be used to calculate superheaters and other cross-flow tubular heat exchangers with complicated flow systems. The advantage of the method is that it requires little effort to build a mathematical model of the heat exchanger and has a very short computer calculation time. Because of the short calculation time, the superheater model developed in this paper can be used for online monitoring of superheater operation. By using the model, steam and wall temperature distributions can be determined along the length of the superheater tubes, which would be difficult to measure directly due to the high flue gas temperatures.
2. A Mathematical Modelling Method for Tube Cross-Flow Heat Exchanger
The calculation of the temperature distribution of the two fluids in the heat exchanger will first be presented, and then the procedure for determining the pressure drop inside the tubes will be described.
2.1. Determination of the Temperature Distribution of Fluids
The method presented in [
23] was used to simulate the third-stage steam superheater in a supercritical pressure boiler. The tubular heat exchanger was divided into finite volumes (
Figure 1). Formulas were obtained for calculating the temperature of the fluid
a and
w at the outlet of the finite volume. The energy conservation equations for steam flowing through the tubes and flue gas flowing perpendicular to the tube axis were used to derive formulas for outlet temperatures of both fluids.
The symbols
and
in
Figure 1 designate inlet temperatures of fluid
w and
a, and
and
are outlet temperatures of fluid
w and
a, respectively. The symbol
designates the mass flow rate of fluid
a per a single finite volume of length
. The symbol
represents the mass flow rate of fluid
w per a single tube. The
p2 symbol indicates the longitudinal tube spacing (pitch) in the heat exchanger.
If the inlet temperatures
and
are known, the temperature
is calculated using the following expression:
where
and
denote the number of heat transfer units (NTU) for fluid
a and
w, respectively.
The NTU values for the
i-th finite volume are calculated as follows:
The overall heat transfer coefficient refers to the external surface of the tube with an area of , where the symbol denotes the tube’s outer diameter.
The mean specific heat capacities
and
are defined as follows:
If the temperature
and the temperature
are known, the temperature
is calculated using the following formula:
Formula (1) applies when the temperature of the fluid w is calculated in the direction of flow. Formula (6) applies when the temperature of the fluid w is calculated opposite to its flow direction.
If the temperatures
and
are known, then the temperature
of the fluid
a at the outlet of the finite volume is calculated from the formula:
where the average fluid temperature
w on the finite volume length
is defined as
In the case of a uniform tube division in one pass of length L into n finite volumes, the finite volume length is . In order to avoid iteration, it is usually assumed , for Equation (1) and Equation (6).
Formulas (1), (6), and (7) apply to determine the outlet fluid temperatures of each finite volume. It should be noted that solving a system of algebraic equations is not necessary when determining the node temperatures in all finite volumes. It is an essential advantage of the method used.
The temperature
of the tube at the outer surface in the
i-th finite volume in the
j-th row is given by
where the heat flux at the tube outer surface is given by the expression:
The overall HTC
referred to the tube outer surface is calculated as follows:
The following designations are used in Formulas (9)–(11): i and j—the number of finite volumes and the number of tube rows,
—mean temperature of the fluid w across the length of the finite volume, °C,
—heat flux at the tube outer surface for the i-th finite volume and j-th tube row, W/m2,
, —inner and outer tube diameters, m,
—HTC at the inner surface of the tube (on the fluid side w), W/(m2·K),
—HTC at the outer surface of the tube (on the fluid side a), W/(m2·K),
—tube material thermal conductivity for i-th finite volume in j-th tube row, W/(m·K),
—number of tube rows in the heat exchanger.
The average flue gas temperature over the thickness of the
j-th row of tubes in the
i-th finite volume was calculated using the following formula:
The HTC at the tube’s inner surface
was calculated from the Dittus–Boelter formula [
18,
19]:
The HTC
can be calculated more precisely using the simple correlations given by Taler [
24] and Taler and Taler [
25], the form of which depends on the value of the Prandtl number.
2.2. Determination of Pressure Distribution
The momentum conservation equation for steady one-dimensional fluid flow in a channel shown in
Figure 2 is given by Taler [
26]:
where
w—fluid velocity, m/s,
pw—static pressure, Pa,
s—the coordinate passing through the centre of gravity of the duct cross-section, with the same direction as that of the velocity vector, m,
—density of the fluid, kg/m3,
g—acceleration of gravity, m/s2,
φ—inclination angle of the duct axis to the horizontal plane, rad,
—Darcy–Weisbach friction factor, -,
dh—hydraulic diameter, m.
In the component representing the friction pressure drop, the absolute value of the velocity is present to take into account the return of the velocity vector. It should be noted here that the pressure drop occurs in the direction of the flowing fluid. The momentum conservation Equation (14) can also be written as a function of the mass flow rate of the flowing fluid
:
where
A is the cross-sectional area of the channel, m
2.
Equation (14) can be written in the form:
The first term in square brackets in Equation (16) represents the pressure drop due to the variation in the velocity of the flowing fluid, the second the gravitational pressure loss, and the third the pressure drop due to friction. Equation (16) was solved using the finite difference method. The differential Equation (16) was approximated using the finite difference method:
The subscript
i refers to the inlet to the finite volume and index (
i + 1) to the outlet. The mean values of fluid density and flow velocity were assumed to be equal to the corresponding values at the inlet to the finite volume. The pressure at the outlet of the finite volume determined from the solution of Equation (17) has the following form:
The fluid flow velocity
wi at the inlet to the
i-th cell (finite volume) was calculated as follows:
It should be noted that in the steady state, .
Local pressure losses, such as pressure drops at bends, valves, and constrictions or expansions of the tube cross-section, are not included in the Equation (18).
Pressure drops at local resistances are usually taken into account at the boundary between two finite volumes. For example, for a modelled four-pass superheater, there are two 90° bends where steam flows from one pass to the other.
If the steam pressure
pw2,n+1 in the first tube row (in the fourth superheater pass) and the pressure
pw1,1 in the second tube row (the third superheater pass) are known at the finite volume nodes, the pressure drop across the two bends between the second and first tube row is given by the formula (
Figure 3). After taking into account that the pressure
pw2,n+1 is known from the calculation of the second tube row (third pass), the pressure
pw1,1 at the inlet to the first tube row (fourth pass) is calculated:
where the symbol
denotes the pressure lost coefficient for a 90° bend.
Similarly, we proceed with the determination of the pressure drop at the bends occurring between each row.
The values of the coefficients
can be found for a number of elements in the work of Lokshin et al. [
18] and in the handbook of Idelchik [
27]. The friction factor
can be read from Moody’s [
28] diagram. In computer calculations, it is more convenient to use explicit relations to determine the friction factor in tubes with rough surfaces. The explicit Churchill formula modified by Rennels and Hudson [
29] is very accurate for determining the friction factor:
where the coefficients
A and
B are given by
The symbol ε in Equation (22) designates relative roughness equal to , where the symbol Ra denotes the absolute roughness of the inner tube surface.
The method described above was applied to the determination of the steam and flue gas temperature distribution and pressure drop inside the tubes in a four-pass counter-flow superheater with supercritical steam pressure.
3. Mathematical Model of Cross-Counter-Current Superheater with Four Passes
The modelled third stage of the live steam superheater consisted of 88 four-pass panels. A side view of a single panel is shown in
Figure 3a. In a single panel (platen) with 4 passes, steam flowed parallel through 12 tubes. The height of a single panel was 5318 mm. There were
tubes (48 horizontal tube rows) per panel height. Considering that there were 88 panels in the superheater, the steam flowed parallel through
tubes. The transverse spacing of the panels was
mm, and the longitudinal tube spacing in the panel was
mm. The outer diameter of the superheater tubes was
, and the inner diameter was
. The average tube length in each superheater pass was
. The length L was calculated by dividing the total length of all tubes in the panel by
.
In the thermal–fluid calculations of the superheater, it was assumed that the length of the tubes in each superheater pass was the same. The equivalent length of a single pass was determined by dividing the total length of all tubes in the superheater by the number of parallel tubes in the superheater equal to .
The flow arrangement of the superheater under consideration with division into finite volumes is depicted in
Figure 4. Each superheater tube was divided into
n finite volumes. The superheater calculations were carried out for the given mass flow rates of steam
and flue gas
. At the superheater inlet, the pressure
and the steam temperature
were known. The flue gas temperature at the superheater inlet was uniform and was
, i.e.,
. The calculation of the steam temperature started with the outlet from the first row of tubes as the heat exchanger was cross-flow-counter.
The following designations are used in
Figure 4:
—flue gas temperature at the
i-th node upstream of the first tube row;
,
,
, and
—flue gas temperature downstream of the first, second, third, and fourth row of tubes, respectively;
,
,
, and
—the steam temperature at the
i-th node of the first, second, third, and fourth row of tubes, respectively.
The calculation of the superheater shown in
Figure 4 can be carried out in two different ways depending on the data available. If the outlet temperature of the steam
is known, it is possible to step in the opposite direction of the steam flow by calculating the temperatures
for
first. Then, the steam temperature at the outlet of the second tube row is known from the first tube row calculation, i.e., the following condition is assumed for the calculation
. The steam temperatures
at the nodes
are calculated sequentially by stepping from the right side of the superheater to the left. In the third row of pipes, it is assumed that
, where the temperature
was obtained from the calculations of the second tube row. For the fourth pass, the following boundary condition is assumed:
. By calculating the temperatures at nodes successively,
, one finds the steam temperature
at the superheater inlet. The HTC on the flue gas side is determined from the condition of equality of the calculated
and measured
steam temperature at the superheater outlet, i.e., from the condition
. The
temperature can be determined at each iterative step from the expression for heat exchanger efficiency
for a given HTC
on the flue gas side. The P-NTU method was used to determine the efficiency of the
superheater. The
efficiency was determined using the formula given in [
30].
In the second case, when the steam temperature
at the superheater outlet is not known, the determination of the steam temperature
at the superheater inlet proceeds as follows. First, the HTC
on the flue gas side is calculated. The formulae given in the standard method [
19] can be used to calculate the HTC. Formulae given in [
19] take into account the radiation [
19,
21] and convection [
18] components of the HTC. The determination of the
steam temperature is carried out iteratively. The temperature
is chosen so that the temperature
calculated by the proposed method and the set temperature
are equal to each other. The symbol
denotes the design or measured steam temperature at the superheater inlet.
Calculations of steam temperature, flue gas temperature, tube wall temperature, and steam pressure were carried out with a single superheater pass divided into
finite volumes. The physical properties of the steam were calculated using the equations given in the ASME steam tables [
30] and by Wagner and Kretzschmar [
31]. The equations given by Brandt [
32] were used to calculate the flue gas properties as a function of temperature. An in-house calculation programme in FORTRAN 90 was developed for the flow–thermal calculations of the superheater under analysis.
The boiler burned hard coal with the following mass proportions of the individual components: C = 52.6%, H2O = 20.8%, H2 = 3.64%, S = 0.9%, and N2 = 0.79%. The mass flow rate of hard coal was , and the excess air number at the combustion chamber outlet was equal to . The calculated and measured flue gas flow rates were equal: . The steam and flue gas temperatures determined during measurements carried out at full load of the boiler were as follows: the steam temperature at the superheater inlet , the steam temperature at the superheater outlet , average flue gas temperature before the superheater , and average flue gas temperature after the superheater . The steam pressures at the inlet and outlet of the superheater were equal to and , respectively. The absolute roughness was assumed to be . The thermal conductivity of the tube material was equal to 25.76 W/(m·K). The HTC at the inner surface of the tubes calculated from Equation (13) was . The HTC on the flue gas side was determined from the experimental data presented above so that the superheater effectiveness was equal to the effectiveness calculated using the measured steam and flue gas temperatures.
The pre-selected value of the HTC using the P-NTU method was improved using the proposed numerical superheater model. The value of was adjusted so that the steam temperature calculated and measured at the outlet of the superheater were equal to each other.
The effectiveness of the four-pass steam superheater shown in
Figure 3 can be calculated, assuming constant physical properties of the steam and flue gas, using the following formula given by Thulukkanam [
33] and by Nicole [
34]:
where the steam inlet temperature is
(
Figure 4). The symbols
and
designate the mean flue gas temperature before and after the superheater. The thermal efficiency of the superheater on the flue gas side is given by the formula [
25,
26]:
where
is defined as follows:
The symbol
in Formula (25) denotes the following expression:
where
K is calculated from the following relationship:
The HTC on the flue gas side determined from the equality condition of the superheater efficiency was calculated using the Formula (24), and the efficiency determined from experimental data was .
It can be seen (
Figure 5) that the highest steam temperature occurred in the first tube row (fourth pass) and the lowest in the fourth tube row (first pass). The average temperature differences between the outer surface of the tube and the steam along the length of the individual tube rows were
in the first tube row,
in the second tube row,
in the third tube row, and
in the fourth tube row. The temperature difference between the outer surface of the tube and the steam decreases in successive tube rows as the highest temperature difference between the flue gas and steam occurred in the first tube row and the lowest in the fourth tube row. The average temperature differences between the flue gas and the outer surface of the tube along the length of the individual tube rows were
in the first tube row,
in the second tube row,
in the third tube row, and
in the fourth tube row. The differences between the flue gas temperature and the tube wall temperature in the individual tube rows were much greater than the respective differences between the wall and steam temperatures. This is due to the significantly higher heat transfer coefficient on the steam side’s inner tube surfaces than on the outer tube surface on the flue gas side. This also reduced the heat flux at the outer surface of the tubes, which was transferred from the flue gas to the steam. The heat fluxes on the individual tube rows were
in the first tube row,
in the second tube row,
in the third tube row, and
in the fourth tube row.
The heat flux values were calculated at the centre of the tube length. The calculated and measured steam temperature at the outlet of the first tube row (fourth pass) was . The temperatures were the same because the steam temperature at the outlet of the first tube row was taken as the input for the calculation. The calculated steam temperature at the superheater inlet was equal to the measured value.
Figure 6 shows the flue gas temperature distributions along the length of the superheater downstream of the individual tube rows. From an analysis of the results shown in
Figure 6, it can be seen that the flue gas temperatures behind the different tube rows varied a little over the width of the superheater. The differences in flue gas temperature at points
and
behind the individual tube rows were minor: 3.46 K, (−0.03) K, 2.30 K, and (−0.04) K behind the first, second, third and fourth rows of tubes, respectively. In the calculations, the measured flue gas temperature of 731.65 °C was used as the flue gas temperature at the superheater inlet. The average flue gas temperature after the superheater calculated using the proposed method was
, while the measured temperature was 596.57 °C.
The heat flow rate
transferred in the superheater from the flue gas to the steam and the heat flow rate
taken up in the superheater by the steam were then compared to assess the accuracy of the developed model. To further assess the validity of the developed mathematical model of the superheater, the average outlet temperatures of the steam and flue gas were calculated using the P-NTU method. The relative difference
between the calculated and measured temperature difference of the flue gas behind the superheater to the measured flue gas temperature difference on the superheater was
. The heat flow rate calculated for the flue gas side was
while the heat flow rate calculated from the steam side was
In Formulas (29) and (30), the following designations were used:
, —the mass flow rate of flue gas and steam through the entire superheater,
, —average flue gas temperature before and after the superheater.
The relative difference between
and
was
The low value of the relative difference demonstrates the perfect accuracy of the numerical superheater model.
The superheater outlet steam temperature and the mean flue gas temperature after the superheater were also calculated using the P-NTU method [
33,
34]. When calculating the
Pa efficiency from Equation (25), it was considered that the average values of the number of heat transfer units were equal to NTU
w = 0.3826 and NTU
a = 0.8617 for steam and flue gas, respectively. The superheater efficiency calculated from Equation (25) was 0.5206. The average flue gas temperature after the superheater
was determined using the definition of the effectiveness given by Equation (24):
Using Equation (32), the following value
was obtained:
. The temperature of the steam at the superheater outlet was determined by the condition of equality of the heat flow rate transferred from the flue gas to the steam. The outlet steam temperature
is given by
After substituting the data into Equation (33), the following value of outlet steam temperature was obtained: . The values obtained for flue gas and steam temperatures downstream of the superheater agreed with the corresponding values obtained from the proposed superheater model despite assuming constant specific heat of steam and flue gas in the P-NTU method. However, the P-NTU method does not allow the determination of the tube wall temperature in each pass, which is needed when selecting the steel grade for a given pass. Due to the high cost of alloy steels, the individual superheater runs were made from different steel grades matched to the tube wall temperature in a given pass.
Figure 7 shows the variation of steam pressure in the superheater along the superheater steam flow path.
The steam pressure at the superheater inlet was
, and at the outlet
. The calculated steam pressure drop in the superheater was equal to the measured steam pressure drop
because the absolute surface roughness of the superheater tubes was chosen so that the calculated and measured steam outlet pressures were the same. From the analysis of the results shown in
Figure 5,
Figure 6 and
Figure 7, it can be seen that the developed numerical model of the four-pass counter-flow steam superheater has excellent accuracy.