Further Study of Gas Temperature
Further Study of Gas Temperature
www.fuelfirst.com
Abstract
Gas temperature deviation in upper furnace is an important but a less reported issue in large-scale tangentially fired boilers, since they
endanger largely boilers operation. Simulations are conducted in this paper to study the deviation. Perfect agreement between the simulation
results and key boiler design values and available site operation records indicates that the calculations are reliable. Based on the simulations,
effect of some factors, including residual airflow swirling at furnace exit, super-heaters panels, coal particle trajectories and their combustion
histories, on temperature deviations are studied in details. The most important cause and how it affects the temperature deviation are located.
Two new methods, a nose on front-wall and re-arranged super-heater panels, are put forward unprecedentedly to alleviate the deviations.
q 2003 Elsevier Science Ltd. All rights reserved.
Keywords: Tangentially fired; Pulverized coal; Boiler; Gas deviation; Simulation
1. Introduction believed that they are related to steam flow deviation among
panels on the steam side and uneven gas temperature
Tangentially fired boiler is the most commonly used distribution in upper furnace and in crossover pass.
industrial coal combustion mode and, thus gains much As a special topic, the gas temperature deviation in large-
attention. With the development of comprehensive com- scale tangentially fired boilers was studied [12]. The
bustion models and their successful applications in modeled gas temperature deviations were compared with
industrial coal-fired boilers, as surveyed in Refs. [1,2], the measured values in a 600 MW boiler, and an
the complex phenomena in tangentially fired boilers, approximate agreement between calculations and measure-
including gas – solid flow, combustion, heat transfer and ments was found. It was concluded that residual airflow
NOx reduction, are widely studied with simulations and swirling and platen super-heaters were the two important
better understandings are achieved [3 –8]. In some cases, it causes for the gas temperature deviation.
replaces physical experiments with equivalent ‘numerical’ Following in the wake of earlier studies, this paper is
experiments. More recently, Fan et al. modeled system- devoted to investigating in details the effects of some
atically the combustion processes in a large-scale tangen- factors, including residual airflow swirling at furnace exit,
tial boiler (600 MW), including the flow, combustion and existence of super-heaters panels, coal particle trajectories
NOx emissions, and used their achievements to improve and their combustion histories inside and beyond combus-
boilers’ performance [9 –11]. tion section, on temperature deviations in tangentially fired
Compared with the above achievements, studies on gas boilers, based on which possible measures are put forward
temperature deviations in upper furnace in tangentially fired to alleviate the deviations.
boilers are less reported. In fact, it is also an important issue.
Tube overheating and explosion of super-heaters and/or re-
heaters rank in the first place among all the operational 2. The 609 MW utility boiler
accidents of large-scale tangentially fired boilers. It is
* Corresponding author. Tel.: þ 45-9635-9248; fax: þ 45-9815-1411. The 609 MW boiler is shown schematically in Fig. 1.
E-mail address: [email protected] (C. Yin). The lines, P1, P2, P3, P4, V1, V2 and V3, denote some
q
Published first on the web via Fuelfirst.com—http://www.fuelfirst.com cross-sections in upper part of the boiler. The dots L0 – L4
0016-2361/03/$ - see front matter q 2003 Elsevier Science Ltd. All rights reserved.
doi:10.1016/S0016-2361(02)00418-0
1128 C. Yin et al. / Fuel 82 (2003) 1127–1137
Table 1
Analysis data of the coal fired
Proximate analysis (%, dry) Ultimate analysis (%, DAF) Lower heating value (MJ/kg, DAF)
3. Simulation methods
Fig. 2. Concentric firing system in the furnace. The calculation domain is shown in Fig. 1. Since the flow
is still quite unmixed at the furnace exit, and the gas
temperature deviation in upper part of the boiler is the
emphasis of this study, so the upper furnace, the crossover
pass and part of the rear pass are included in the simulations.
Special attention has been paid to grids and grid design to
reduce numerical error as much as possible [15]. A
structured mesh is created with boundary-fitted coordinates,
in which the grid lines are approximately along the swirl
flow direction in furnace, as shown in Fig. 2. Fig. 4 also
shows the calculated flow field at a plane in the furnace,
from which one can see that the mesh resolves the major
flow property in the furnace. The flow in the crossover pass
and the rear pass is relatively simple and therefore, the
commonly used quadrilateral grid is kept. There are totally
454,776 hexahedral cells, with 42 £ 37 £ 234 ¼ 363,636
cells in the furnace. The grid is fine enough to give grid-
independent solutions [11]. First-order spatial discretization
scheme is used for all equations. Validating the simulations
by some key global design parameters of the boiler and
available site operational records assures reliable and
acceptable simulation results, as well as the residuals
control.
Table 2
Size distribution of pulverized coal
combustion, creating a source of fuel for reaction in gas the upper furnace, in the crossover pass and in the rear pass
phase. The two-competing-rates model is used to yield an are taken into account with the method described in Section
expression for the devolatilization. The heterogenerous 2. To investigate the effect of the panel super-heaters on the
combustion of char is calculated using kinetics/diffusion- gas temperature deviation, a second case, i.e. case 2, is also
limited model. The interaction between the gas and particle conducted, in which all platen heaters in the upper furnace,
phases is modeled by altering iterations between the flow in the crossover pass and in the rear pass are ignored. In this
calculation and the particle trajectory calculation, with 20 paper, case 1 is the default case for results presentation.
flow iterations per discrete phase iteration.
Species and chemical reactions are modeled using the
mixture-fraction/probability density function (PDF) 4. Results and discussions
approach and the full equilibrium chemistry, where the
turbulence-chemistry interaction is modeled using a double- 4.1. Calculation validation with boiler design values
delta PDF.
The discrete ordinates (DO) model is used to simulate The combustion calculation begins by solving the gas
radiation heat transfer, thanks to its ability to take into flow field equations assuming that the particles are absent.
account the particulate effects and localized heat sources, This simulation is well converged after more than 3000
the two important aspects of coal combustion simulation. iterations, as shown in Ref. [14], from which, full
All the models used in this study are general models, combustion simulations start.
widely used and well documented, and can be found in Refs. In the combustion calculations, the furnace exit par-
[2,16] for details. They are proven efficient to capture the ameters, such as gas temperature and O2 concentration on
complex phenomena in such a large-scale tangentially fired plane P1, are monitored, since they are very important for
boiler. The calculations in this paper just follow in the wake preserving the boiler operation and performance including
of the earlier studies [13,14]. emission, reliability and safety. For example, if furnace exit
gas temperature deviates from the design value, some
3.3. Calculation cases undesired condition might occur, such as, potential of
crossover pass tubes overheating, and increased slagging/
The basic case, i.e. case 1, is conducted on the basis of corrosion of water walls and super-heater tubes. The furnace
full-load operation. In this case, all the platen heaters in exit parameters are also used to validate the combustion
calculations. Table 3 shows the comparison of the
calculation (case 1) with some key design values, including
the furnace exit gas temperature and O2 mass fraction. The
designed furnace exit gas temperature and O2 mass fraction
are about 1600 K and 3.5%, respectively; and the total heat
transfer in the furnace is around 609 MW. The simulation
results in case 1 agree very well with the design values.
Here, it should be stated that the emissivity of the furnace
wall is set to 0.7 in this study for the radiation heat transfer
simulation, since the average temperature on plane P1 is
quite dependent on this parameter.
Table 3
Validation of combustion simulation with key global design values
Fig. 7. Gas velocity vectors on planes in the upper furnace: P1, P2, P3 and P4.
because of the residual airflow swirling. With the height the upper furnace and in the crossover pass: low in the left
increases, the gas in the left gradually reverses its and high in the right.
direction under effect of furnace walls and super-heaters From Fig. 7, large vx components from the left to the
panels, and flows forward directly into the crossover right are found in the gaps between the front-wall and the
pass; while the velocity component, vy, in the right has a first group of platen super-heaters (SH1) in all planes of
slight decrease because less gas goes from the left into different heights in the upper furnace. This portion of gas is
the right and then into the crossover pass. from two sources. Part of the gas comes directly from the
Fig. 7 also exhibits a typical uneven gas velocity in the near front-wall zone in the combustion section, which goes
upper furnace: low in the left and high in the right. upward from there at a certain z-velocity into the
Combined with the vy component distribution in plane P3, corresponding zone in the upper furnace; another part is
as shown in Fig. 8, the formation mechanism of the uneven from the left in the upper furnace, which goes firstly
gas velocity in the upper furnace may be explained in Fig. 9, backward toward the front-wall and then into this zone.
where M1 denotes the momentum of the residual airflow Both of them are caused by the residual airflow swirling at
swirling in upper furnace; M3 is the momentum introduced the furnace exit, which is probably one of the most
by induction fan, which can be supposed as constant here important causes for the gas temperature deviation in the
along the width of crossover pass; and M2 represents the upper furnace: high temperature in the left and low
composite momentum of gas in crossover pass. That is temperature in the right, and this temperature deviation
the mainstream of the uneven gas velocity distribution in increases with the furnace height.
C. Yin et al. / Fuel 82 (2003) 1127–1137 1133
Fig. 10. Two new methods to alleviate the gas temperature deviation in upper furnace.
swirling affects the temperature deviation in the upper the residual airflow swirling. This effect can be observed
furnace: the existence of high-speed gas flow from one side- from Fig. 11, in which the calculated deviations on lines of
wall to the opposite side-wall in the gap between the front- different heights in the two cases are shown.
wall and the platen super-heaters. The first is, a new nose The platen super-heaters attenuate the residual airflow
might be added on the front-wall near the furnace exit (P1 in swirling at furnace exit to some extent and reduce the
Fig. 1). In this way, the gas flow along the relatively cool velocity deviation between the two side-walls, as seen in
front-wall in the gap will be reduced under the directing of Fig. 12. In addition, super-heater panels can also attenuate
the nose. The second is, the super-heater panels can be re- the gas velocity imbalance in upper furnace and improve the
arranged to distribute uniformly the gas between two side- flow distribution.
walls. Taking this boiler as an example, they can be
described in Fig. 10. 4.5. Effect of particle trajectories and their combustion
histories on gas temperature deviation
4.4. Effect of platen super-heaters on gas
temperature deviation In tangentially fired boilers, the particles sweep around
the furnace volume, and it provides sufficient residence time
As stated in Ref. [12], the platen super-heaters in the for the very fine particles. During particle combustion, all of
upper furnace are another important cause for the gas the fuel species except char are consumed within 1 or 2 m of
temperature deviation. But, how the super-heater panels the burner throats, while char burns at a slower rate, and
affect the temperature deviation needs to be further studied. normally is not consumed until it reaches the furnace center.
From the above cause analysis of the temperature At the furnace exit cross-section (i.e. plane P1), the gas
deviation, it may be concluded that super-heater panels in temperature is fairly uniform between the two side-walls
the upper furnace and in the crossover pass can reduce the and the temperature deviation between the two side-walls is
temperature deviation to some extent since they attenuate very small, which was confirmed in Refs. [13,14] and can
C. Yin et al. / Fuel 82 (2003) 1127–1137 1135
Fig. 11. Gas temperature deviation on line L0, L2 and L3 in different cases.
also be seen from the temperature plot on lines L0 – L3 in section (i.e. above the furnace exit), the char concentration
Fig. 13, thanks to the reasonable design of the concentric is very low, with the calculated mass fraction less than
firing system in the furnace, and the distance between the 2 £ 10210, and thus the effect of char combustion in the
upmost burner and the furnace exit, and thus the pretty long upper furnace is very small, which can also be confirmed by
period of intensive mixing in the furnace. Therefore, it may the very low mass-exchange between two phases. So char
be concluded that one can ignore the influence of the combustion in the upper furnace is very limited and may be
particle trajectories and their combustion histories inside ignored, which can also be reflected from the monotonous
the combustion section at all when studying the formation of decrease of gas temperature with increase of section height
the gas temperature deviation in the upper furnace. (shown in Fig. 13).
Although char combustion might still exist somewhere Due to the quite complex flow in the upper furnace, the
in the upper furnace, the zones with char combustion are particle concentration and enthalpy change are here exam-
very small and the combustion rate is also very low: about ined instead of particles trajectories to study how the particles
two orders of magnitudes lower than those in the affect the gas temperature deviation in the upper furnace. The
combustion section. Moreover, beyond the combustion particle concentration is found a little higher in the right,
Fig. 14. Total enthalpy exchange from particles to gas in the upper furnace.
because gas velocity is a little higher in the right and particles the total mass flow rate in the boiler is about 660 kg/s. The
in upper furnace are very tiny. But this deviation is not accumulated enthalpy-exchange deviation itself in the upper
significant. Fig. 14 shows the exchange of total enthalpy from furnace is estimated to cause a low gas temperature deviation
the particles to the gas in the upper furnace. In fact, it is between the two side-walls, less than 15 K.
almost same as the sensible enthalpy exchange since the char So it may be concluded that the particle trajectories and
combustion is negligible. From the distribution of the combustion histories both inside and beyond the combustion
enthalpy exchange, one might be able to conclude that section do not have a significant effect on the temperature
the particles trajectories in the upper furnace will not cause a deviation in the upper furnace.
big gas temperature deviation, since both the transferred
enthalpy values and their deviations between the two side-
walls are fairly low, especially with increase of the furnace 5. Conclusions
height. If ignore all other factors, a 100 W deviation in the
enthalpy transferred from the particles to the gas will have Simulations are conducted to study the gas temperature
little effect on the gas temperature deviation since the specific deviations in a large-scale tangentially fired boiler,
heat of the gas is about 1300 J/(kg K) in the upper furnace and especially in the upper furnace. Good agreement between
C. Yin et al. / Fuel 82 (2003) 1127–1137 1137
the simulation results and key design values of the boiler proposed to reduce gas temperature deviation in
and available site operation records indicates that the upper furnace in tangentially fired boilers, with special
calculations are reliable for such an investigation. Generally aims at reducing hot gas flow going from one side-wall
speaking, the gas temperature deviation in this 609 MW to another along the cool front-wall and distributing
boiler is not serious, thanks to reasonable boiler design. this gas flow uniformly between the two side-walls.
Some special conclusions can be drawn as follows.