This paper establishes a compact and practical model for a water-wall system comprising supercritical once-through boilers, which can be used for automatic control or simple analysis of the entire boiler-turbine system. Input and output variables of the water-wall system are defined, and balance equations are applied using a lumped parameter method. For practical purposes, the dynamic equations are developed with respect to pressure and temperature instead of density and internal energy. A comparison with results obtained using APESS, a practical thermal power plant simulator developed by Doosan Heavy Industries and Construction, is presented with respect to steady state and transient responses.
In spite of environmental issues, thermal power plants generate approximately 65% of the world’s power supply. In recent years, the construction of large-capacity thermal power plants with environmental facilities has been common
With respect to structure, the boilers of thermal power plants are classified into two types: the drum boiler and the once-through boiler (OTB)
shows a schematic of a supercritical once-through boiler-turbine system. A supercritical once-through boiler comprises several heat exchangers, such as economizers, a water wall, superheaters, and reheaters. As shown in
, the feedwater enters into the water wall through economizers. Then, the water is transformed into steam in the waterwall tube. The steam is superheated to generate electric power and circulated to the economizer again. Alternatively, a drum-type boiler system includes a drum, wherein saturated steam is separated from saturated fluid and provided to a superheater. The remaining saturated water re-enters the water-wall tubes through downcomers
Schematic of a supercritical once-through boilerturbine system
Currently, once-through boilers are constructed more commonly than drum boilers. Compared with drum boilers, once-through boilers can operate at higher pressure and temperature and thus allow greater energy efficiency. Because once-through boilers do not have drums or large-diameter downcomers, they exhibit less metal weight and smaller fluid storage capacity than drum boilers
. Therefore, although once-through boilers can respond rapidly to load changes, controlling them is more difficult than controlling drum boilers
In the system represented by
, the entire surface of the lower part of the furnace wall is surrounded by waterwall tubes. When the operation conditions of the boiler exceed the critical point (22.09 MPa, 374.14℃
), the unit is called a “supercritical unit.” In the water-wall tube of a supercritical once-through boiler, the phase of the water changes directly from liquid to vapor without undergoing saturation. This is a significant difference between the supercritical once-through boiler and other subcritical boilers. The water-wall system is one of the most important components affecting the dynamics of supercritical once-through boilers.
Although there are well-established models for drum boilers, such as those proposed by Bell and Åström
, standard models for once-through boilers are far less common. Because of the major difference in the water-wall system between the two types of boilers, a compact and effective model of the water-wall system is currently a relevant research topic.
There are many mathematical models of a water wall for subcritical once-through boilers
; however, there are comparatively few mathematical models of a water wall for super-critical once-through boilers.
Dumont and Heyen developed an abridged mathematical model for the entire once-through boiler system
. They modified internal heat transfer coefficients and pressure drop formulations and considered the changes in the flow pattern. Li and Ren describe a water-wall system using a moving boundary
. They used enthalpy to track the moving boundary location at supercritical pressure and used mass, energy, and momentum balances to obtain the length of each section. Pan and colleagues presented a detailed water-wall model for predicting the mass flux distribution and metal temperature in the water wall of an ultra-supercritical boiler
. They treated the water-wall system as a network comprising 178 circuits, 15 pressure grids, and 7 connecting tubes; the system can be described using 195 non-linear equations.
Recently, intelligent systems have been applied for modelling a once-through boiler. Chaibakhsh and colleagues developed a model for a subcritical once-through boiler whose parameters are adjusted on the basis of genetic algorithms
. Lee and colleagues established a model for a large-scale power plant based on the neural network method
, and Liu and colleagues described a supercritical once-through boiler using the fuzzy-neural network method
In the present study, we attempt to develop a compact and practical model of water-wall systems for supercritical boilers that can be used for automatic control, analysis, and modeling of entire boiler-turbine systems. The objective is to develop a relatively simple water-wall model with sufficient accuracy for analysis and control rather than to describe the detailed dynamics occurring inside the water-wall tube. We use pressure and temperature as state variables; both of these are practical variables in industrial applications.
First, we establish input and output variables of waterwall systems and apply fundamental laws of physics, i.e., mass, energy, and momentum balance equations, using a lumped parameter method. Then, complicated equations and variables are approximated by adopting reasonable and applicable assumptions. To change the state variables with pressure and temperature, enthalpy and density are approximated as functions of pressure and temperature using a steam table. To verify the proposed model, a model of the water-wall system obtained using APESS, a practical thermal power plant simulator
developed by Doosan Heavy Industries and Construction, is presented and compared.
2. Basic Balance Eqs.[3,7,19-22]
The fundamental principles used in developing the model are mass balance, energy balance, and momentum balance.
shows the nomenclature used in this paper. In
, “wall” denotes the tube wall of each heat exchanger, such as the water wall, superheater, and reheater.
- 2.1 Mass balance
Mass balance is represented in (1), which gives the rate of mass change for a heat exchanger system.
- 2.2 Energy balance
Energy balance is represented in (2), (4), and (6) for the combustion gas, tube wall, and working fluid, respectively. The dynamics of combustion gas are represented in (2), where
is the transferred heat flow from combustion gas to the tube wall. As shown in (3),
has two terms: radiative heat transfer and convective heat transfer. The temperature change of the tube wall is represented in (4), where
is the transferred heat flow from the tube wall to the internal working fluid in (5). Therefore, the combustion energy is represented as the temperature change of the tube wall using (2) and (4). Finally, the dynamics of internal working fluid energy are represented in (6).
- 2.3 Momentum balance
The exact momentum balance of fluid in the tube is difficult to describe theoretically because of the internal turbulent flow of fluid. However, in the momentum balance equation, the dynamic term can be neglected because the pressure-flow process works faster than the mass and energy balance dynamics. In addition, the inertia term can be neglected compared with the friction term. These modifications result in the following equation, which is used in
Generally, heat exchangers in boiler systems, including a water wall, superheater, reheater, and economizer, can be modelled using (1) - (7). However, the major variables in these balance equations, such as pressure (
), temperature (
), density (
), enthalpy (
), and internal energy (
), are dependent variables that are functions of thermodynamic state. The thermodynamic state of water is classified into three state regions: the compressed liquid region, saturated liquid-vapor region, and superheated region. The saturated liquid-vapor region is called the “saturation region.”
3. Development of a Water-Wall Model
- 3.1 Balance equations for water-wall systems
The detailed water-wall model considers many variables
; in this paper, several major thermodynamic variables are selected on the basis of a lumped parameter method. To describe the simple water-wall system, several assumptions are required.
The fundamental balance equations are modified according to the above assumptions. Four major variables — mass flow, enthalpy, pressure, and temperature at the outlet — are selected for both inputs and outputs. To consider the combustion energy, the mass flow of fuel (
) is included as an input variable. The selected variables are represented in
1. The pressure dynamics of the flue gas are negligible.
2. The flue gas exhibits ideal gas behavior.
3. The working fluid properties are uniform at any cross section.
4. The heat conduction in the axial direction is negligible.
5. The change in the thermodynamic properties of the internal working fluid is lumped.
6. The heat transfer from the flue gas to the wall is proportional to the combustion heat generated in the furnace.
7. The gas-wall heat transfer dynamics are sufficiently faster than the wall-fluid heat transfer dynamics.
Inputs and outputs of the water-wall model
Inputs and outputs of water-wall system
Inputs and outputs of water-wall system
Therefore, in this paper, the water wall is represented as a 5-input and 4-output system. The fundamental balance Eqs. (1)-(7), are modified as follows.
3.1.2 Mass balance
Because the working fluid enters from the economizer, the inlet of the water wall is the outlet of the economizer. Therefore, the mass balance of the working fluid in the water wall, given by (1), is modified as follows:
3.1.3 Energy balance
Regarding the energy balance of internal working fluid, (6) can be written as follows:
is the transferred heat flow from the tube wall to the fluid, which is modified from (5) as
Regarding the dynamics of
, the temperature of the tube wall, given by (4), is rewritten as follows:
is the transferred heat flow from the flue gas to the tube wall. Then,
Regarding the dynamics of
in (12), (2) is modified as
is included to consider the heat input by the fuel combustion. In (13), because the flue gas comes from an air preheater and leaves the furnace,
is the enthalpy at the air preheater outlet and
is the enthalpy at the furnace outlet. Because there is no mass flow change of the flue gas in the furnace,
in (2) are unified with
is given as follows:
is the calorific value of fuel and
is the fuel mass flow.
The energy balance Eqs. (8) - (14), can be directly used for the water-wall model. However, they require system variables from the other heat exchangers, such as the economizer, furnace, and air preheater, as well as additional system parameters such as heat transfer coefficients, the volumes of the furnace and wall, and the specific heat at constant volume of the wall and gas. Consequently, direct application of (8)-(12) results in a complicated model, which is beyond the scope of this paper.
In this study, to make the model more compact, we assume that the heat transfer from the gas to the tube wall is proportional to the combustion heat (assumption 6). Then, (12) can be expressed as follows:
where α is the ratio of
. Although α can be considered a constant, it is a function of another thermal state
. In this study, α is a function of
, which is the average temperature between two outlets. That is,
The three coefficients
can be determined using the measurement data.
Typically, heat exchange between the gas and the wall is far faster than that between the wall and the fluid (assumption 7). Therefore, the dynamics of
can be ignored in (11)
. Then, (11) is modified as a static equation with
Accordingly, (18) can be expressed using (14) and (15) as follows:
represents the ratio of
, which is equal to the product of
As a result, the energy balance of the working fluid, given by (9), is simply represented using (20) as follows:
3.1.4 Momentum Balance
In the mass and energy balance equations, given in (8) and (21), the output variable
is determined using the momentum balance equation (7). Because the outlet of the water wall is the inlet of the primary superheater, the momentum balance of the working fluid at the primary superheater is given as follows:
represent gravitational acceleration and the gravitational conversion factor, respectively, whose values are approximately 9.80665 [m/sec
] and 9.80665 [kg(mass)·m/kg(weight)·sec
], respectively. The constant 10.1772·10
is included in the denominator to change the units from [kg(weight)/m
] to [MPa].
Although the friction factor,
, in (22) is considered a constant
is proportional to
in practice. In this study,
is selected as a function of
, as follows, to obtain better accuracy of the system:
The two coefficients
can determined using the measurement data. Finally, three balance equations for the water-wall model are given by (8), (21), and (22) using (16), (20), and (23).
- 3.2 Change of state variables withPandT
The established model given by (8) and (21) explains the dynamics of density
and internal energy
of the working fluid. The state variables and the input and output variables are given as follows:
In industrial practice, the pressure
of the working fluid are directly measured and importantly managed. That is, the steam table is necessary to calculate
from measured variables. Because
are measured outputs, they can be directly compared with measured data for a real plant. Therefore, in this study, we set pressure and temperature as state variables of the waterwall system as follows:
To change the state,
in dynamic equations (8) and (21) are set as functions of
in this study. Hereafter, the subscripts of
ρ, U, P, T
are omitted for conciseness.
From the definition of enthalpy
the left side of (21) can be arranged as follows:
Accordingly, (8) and (21) can be written as follows:
Then, a steam table is used to represent
in (33) and (34) as functions of
. Because the objective system operates in the superheated region,
of the superheated vapor region of the steam table are approximated as the following simple polynomial functions of
where the coefficients are determined using the least squares method. These equations are valid only for the operation range used in the least squares method.
Using the chain rule,
and (33) and (34) are written as follows:
are determined using (39) and (41) with simple algebraic calculations as follows:
Then, we can rearrange the final water-wall equations with notations for state, input, and output as follows:
4. Simulation Results
To test the validity of the presented model, the waterwall system obtained using the APESS simulator is modeled as a target system. The presented water-wall model (46) - (55) is realized using MATLAB and a fourthorder Runge-Kutta algorithm is applied for the discrete simulation. Then, the steady-state and transient responses in superheated operation are compared.
For the simulation, three constants,
, are determined using the APESS simulator. The coefficients
are determined using off-line data from APESS in the superheated operation range. The results of interpolation using the least squares method are as follows:
shows the measurements and plot of
shows the measurements and plot of
. These two figures indicate that the interpolation is quite effective when considering real constants. In
, the measurement value of
changes from 2.8·10
according to the operating conditions, which explains why we do not use a constant
in this study.
η as a function of Tave
Fps as a function of Weco
are also determined using the least squares method with the steam table. The regions 410℃ <
< 430℃ and 25
in the steam table are selected to simulate the operation range of APESS. The results of the approximation are as follows:
Eqs. (46)-(59) form the basis for the water-wall system in APESS. To verify the performance of the system, two types of simulations are tested: steady state responses and transient responses.
- 4.1 Steady-state test
For the steady-state comparison, the APESS model is run with fixed electric power generation. Because the APESS system has internal control loops, all variables in APESS are stabilized to a steady state. Then, steady-state values of the 5 inputs and 4 outputs of the water-wall system are obtained from APESS. The same input values are applied to the presented model, and steady-state output values are compared.
shows a comparison between the APESS system and the presented model. In the table, electric power is varied from 1000 MW to 800 MW, with which the boiler operates in the supercritical region. In
are directly proportional to the electric power, whereas
is inversely proportional. Percent errors of outputs are also presented, calculated as follows: |Model- APESS|/APESS. In the table,
exhibit relatively small errors compared with
. The maximum error is 0.27% (1.14 ℃) for
with power generation of 900 MW. The average of all steady state errors is calculated to be 0.05%. Although there is no absolute criteria to determine modeling mismatch, we believe that these results are sufficient for predicting the steady state of the water-wall system.
Steady-state values of APESS and model
Steady-state values of APESS and model
- 4.2 Transient response test
For the comparison of transient responses, the electric load demand of APESS is increased and decreased in steps. The load demand signal is adjusted as follows: 800 MW → 900 MW → 1000 MW → 900 MW → 800 MW. Each step is maintained for 20 minutes to attain a new steady state.
shows graphs of the 5 inputs of the water-wall system obtained using APESS. The 5 inputs shown in
are applied to the presented model.
Five input signals for transient responses.
show a comparison between the APESS model and the presented model. According to these figures, the responses of the four outputs are similar to those of a firstorder system. Considering that (46) and (47) are very complicated, we find that the major dynamics of the waterwall system are quite simple.
Mass flow (Wwwo) of APESS and model
Enthalpy (Hwwo) graphs of APESS and model
Pressure (Pwwo) graphs of APESS and model.
Temperature (Twwo) graphs of APESS and model.
, the responses of
are almost identical, as suggested by the steady-state responses. In
, the initial value of enthalpy is not identical to the APESS data because the enthalpy is calculated using the pressure and temperature obtained by the approximated equation (58). Although the responses
exhibit different steady states, they have similar patterns with similar rising times.
We present a lumped model for the water-wall systems of supercritical once-through boilers. The model has two states, 5 inputs, and 4 outputs determined using a lumped parameter method. A steam table is approximated and used in the model equations to change the state
A water-wall system obtained using the APESS simulator is modeled as a target system. Comparison results consider both steady states and transient responses. In both simulations, the mass flow and pressure exhibited similar results, and enthalpy and temperature exhibited small errors
Although the presented model is quite complex, its dynamics are similar to those of a first-order system. We believe that this model is useful for designing an automatic controller and for analysis of water-wall systems.
This research was supported by the Chung-Ang University Excellent Student Scholarship and by the Chung-Ang University Research Grant in 2013, and by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (Grant Number: NRF-2010-0025555) .
Geon Go received his B.S. degree in Electrical and Electronics Engineering from Chung-Ang University, Seoul, Korea, in 2013. He is currently an M.S Candidate in Electrical and Electronics Engineering at Chung-Ang University. His research interests involve the operation and modeling of fossil power plants and power system analysis.
Un-Chul Moon received his B.S., M.S., and Ph.D. degrees from Seoul National University, Korea, in 1991, 1993, and 1996, respectively, all in Electrical Engineering. In 2000, he joined Woo-Seok University, Korea, and in 2002, he joined Chung-Ang University, Korea, where he is currently an Associate Professor of Electrical Engineering. His current research interests are power system analysis, computational intelligence, and automation.
“An Overview of Modeling and Simulation of Thermal Power Plant”
Proc.of IEEE International Conference on Advanced Mechatronic Systems
“A New Approach Applied to Steam Turbine Controller in Thermal Power Plant”
2nd International Conference on Control, Instrumentation and Automation
“Dynamic Modelling of Heat Transfer Processes in a Supercritical Steam Power Plant”, M.S. thesis
Dept. Energy and Environment, Chalmers University of Technology
“Estimating the Parameters of Dynamic Model of Drum Type Boilers Using Heat Balance Data as an Educational Procedure”
IEEE Trans. on Power Systems
DOI : 10.1109/TPWRS.2010.2061879
“Mathematical Model of Once-Through Boiler Dynamics”
IEEE Trans. on Power Systems
DOI : 10.1109/TPAS.1965.4766165
“Modelbased Once-through Boiler Start-up Water Wall Steam Temperature Control”
Proc. of IEEE International Conference on Control Applications
Anchorage, Alaska, U.S.A.
Van Wylen G.J.
Fundamentals of Thermodynamics
John Wiley & Sons, Inc.
Åström K. J.
Dynamic models for boiler-turbine-alternator units: Data logs and parameter estimation for a 160 MW unit
Lund Institute of Technology
“A Nonlinear Dynamic Model of a Once-Through Subcritical Steam Generator”
ASME Trans. on Dynamic Systems, Measurement, and Control
DOI : 10.1115/1.3427046
“Moving Boundary Models for Dynamic Simulations of Two-Phase Flows”
Proc. Of the 2nd Int. Modelica Conference
“Distributed parameters modeling for evaporation system in a once-through coal-fired twin-furnace boiler”
International Journal of Thermal Sciences
DOI : 10.1016/j.ijthermalsci.2011.07.010
“Mathematical modelling and design of an advanced once-through heat recovery steam generator”
Computers & Chemical Engineering
DOI : 10.1016/j.compchemeng.2004.02.034
“Moving Boundary Modeling Study on Supercritical Boiler Evaporator: By Using Enthalpy to Track Moving Boundary Location”
Power and Energy Engineering Conference
“Mathematical modeling and thermalhydraulic analysis of vertical water wall in an ultra supercritical boiler”
Applied Thermal Engineering
“A simulated model for a once-through boiler by parameter adjustment based on genetic algorithms”
Simulation Modelling Practice and Theory
DOI : 10.1016/j.simpat.2007.06.004
“Neural Network-Based Modeling for A Large- Scale Power Plant”
Power Engineering Society General Meeting
“Modeling of a 1000MW power plant ultra supercritical boiler system using fuzzy-neural network method”
Energy Conversion and Management
DOI : 10.1016/j.enconman.2012.07.028
Van Sickel J.H.
“Controller Design for a Large-Scale Ultrasupercritical Once-Through Boiler Power Plant”
IEEE Trans. on Energy Conversion
DOI : 10.1109/TEC.2010.2060488
Usoro Patrick Benedict
“Modeling and Sumulation of a Drum Boiler-Turbine Power Plant under Emergency State Control”, M.S. thesis
Dept. Mechanical Engineering, Massachusetts Institute of Technology
Cambridge, United States of America
“A Simplified Model of a Supercritical Power Plant”
University of Michigan
“A lumped parameter dynamic model of the helical coiled once-through steam generator with movable boundaries”
Nuclear Engineering and Design
DOI : 10.1016/j.nucengdes.2008.01.009
“Drum Boiler Reduced Model: a Singular Perturbation Method”
20th International Conference on Electronics, Control and Instrumentation