Design of Heat Recovery Steam Generator
Design of Heat Recovery Steam Generator
Abstract—The gas turbine combined cycle (GTCC) power designs requires general limitations. The boundary
plant system combination a gas power plant and a steam requirement is that the dirt factor (RD) of the two fluids
power plant using Brayton and Rankine cycles. In GTCC flowing in the HE should be greater or equal to the
specification, the heat-recovery steam generator (HRSG) is provisions of the second impurity factor, with pressure
employed as a heat exchanger to produce superheated steam.
drops obtained from piping and heat exchangers [5]. The
The utilization of waste heat is a need for sustainable energy
use. This paper reviews research on the recovery of waste calculation of the impurity factor depends on the model
heat by designing and simulating an HRSG construction. value of the heat transfer coefficient, which is in progress
Therefore, this study aims to create an HRSG with high (56 to predict the behavior of post-boiling transitions that may
bar) and low (6 bar) pressure levels at a temperature of 500°C occur during the anticipated operational events for boiling
and 137.3 kg/s of gas turbines using simulation. The use of [6]. Furthermore, it considers the effect of the most crucial
equations and design simulations can be applied to HRSG parameters, such as surface waves, evaporation, and flow
components with appropriate high and low pressure. configurations [7].
The value of each fluid’s heat transfer is used to
Index Terms—Economizer, evaporator, gas turbine
calculate the overall coefficient (UC), which determines
combined cycle, heat-recovery steam generator, superheater
𝑈𝐷 to obtain RD and provides geometrical design data from
I. INTRODUCTION HRSG [8]. This paper summarizes the detailed HRSG
design scheme and produces its height, length, and width
A Gas turbine combined cycle (GTCC) is a power-plant dimensions. The feasibility analysis of the HRSG design is
system that generates more efficient, significant, and the thickness of the pipes used, which is different from that
improved electricity than that produced by a conventional in previous studies on the Heat Exchanger. The research
cycle, such as a gas power plant and steam power plant. on the system of HRSG for various Na liquids produced
One of the benefits of GTCC is its ability to consume a differences in the heat transfer coefficient ( ℎ𝑁𝑎 ) [9].
limited amount of cooling water and shorter operational Research on the simple modification of the internal
time than those of the two conventional plants. Industries configuration of HRSG also makes the gas combustion
can improve, and high efficiency in GTCC (up to 58%) can flow more uniform [10].
be achieved by improving the heat-recovery steam Furthermore, recovering waste heat from a 60 MW gas
generator (HRSG) design [1]. We focus on introducing turbine exhaust generates additional power of 35.14 MW
multiple pressure levels by reheating steam in the steam [11]. Waste hot air from HRSG also indicated that the air
cycle on HRSG to recover more exhaust gas energy. return strategy offers significant performance
Furthermore, the HRSG has broad application prospects in improvements [12].
increasing electricity capacity and reducing its generation In this work, we propose the design and simulation of
costs. This design aims to further improve system HRSG with high and low pressure by heating the exhaust
performance above 60%. Power plant modifications, heat at a temperature of 500oC. In designing this HRSG,
configurations, and operating strategies can improve the consider the mass-energy balance equation, HRSG, and
unit’s cycling capability [2]. Also, the combined cycle boiler identical in modeling [13]. Apart from that, they also
power plants are increasingly subjected to cyclic include low and high pressures which are required to
operations [3]. GTCC employed HRSG, a Heat Exchanger maximize heat recovery from exhaust gases. The
(HE), on the exhaust gas turbine. Researched include advantage and novelty of this design are that it successfully
economizers, evaporators, and superheaters. Other combines steam and gas power with a capacity of 105 MW
research also focused on enhancement technologies by and an efficiency of 69.48%.
offering more design options for increasing thermal
efficiency [4]. Obtaining optimal performance in HE
II. SIMULATION AND DESIGN TABLE III. DATA FROM THE SIMULATION SUPERHEATER RESULTS TO
START DESIGNING
A. Simulation
This simulation’s focus is to design an HRSG with high
Low High
(56 bar) and low (6.2 bar) pressure levels at a capacity of Pressure Pressure
105 MW. This parameter is the main component of the No Parameter Unit
HRSG, which utilizes the gas turbine’s flue to convert 6 bars 58 bars
water into steam at a temperature of 500 °C. Capacity and
temperature are the initial values for simulating the initial Inlet water
1 ⁰C 140.25 269.06
temperature
design parameters. The simulation results are as shown in Water
Tables I, II, and III. 2 ⁰C 303.8 509.95
temperature out
Specific heat of
3 kJ/kg ⁰ C 4.23 4.25
TABLE I. DATA FROM THE SIMULATION SUPERHEATER RESULTS TO the water
START DESIGNING Water mass flow
4 kg/hr 27898.5 38438.3
rate
Exhaust gas
Low High
5 temperature ⁰C 372.12 650.8
No Parameter Unit Pressure Pressure
entered
6 bars 58 bars
Exhaust gas
6 ⁰C 334.92 572.53
temperature
Inlet water
1 ⁰C 151 265.5 7
Specific heat
kJ/kg ⁰ C 1.06 1.094
temperature exhaust gases
Vapor Exhaust mass
2 ⁰C 160 271 8 kg/hr 494319 494319
temperature out flow rate
Steam mass flow
3 kg/hr 27899 38438.3
rate
Exchaust gas B. Design
4 temperature ⁰C 300 464.4 The design uses the exhaust gas’s base temperature that
entered
comes out of the turbine, which is between 450 °C and
Specific heat kJ/kg
5
exhaust gases ⁰C
0.908 1.139 630°C [14]. Furthermore, the method uses a flue gas
temperature at 500oC and a mass flow rate of 137.3 kg/s
Exchaust a mass
6
flow rate
kg/hr 494319 494319 with the initial dimensions, as shown in Table IV. The
HRSG design uses a flue gas temperature at 500oC and a
Enthalpy of
7 kJ/kg 636.56 1162.49 mass flow rate of 137.3 kg/s and sets the initial design
water enters
Enthalpy of
dimensions as in Table IV.
8 kJ/kg 2756.8 2788.87 The design of the HRSG must go through the stages of
steam out
designing HRSG, which is as follows:
TABLE II. DATA FROM THE SIMULATION EVAPORATOR RESULTS TO
START DESIGNING C. Calculation of Heat Transfer Rate
The equation for calculating the heat transfer rate,
energy balance, and LMTD [15] are:
Low High Heat transfer rate:
No Parameter Unit Pressure Pressure
6 bars 58 bars q U D A tLMTD (1)
1
Inlet vapor
⁰C 160 270.5
q is the rate of heat transfer (W), A is heat transfer area
temperature (m), UD is the design of overall heat transfer coefficient
2
Vapor
⁰C 213 450 (W/(m2) (°C)), and ∆t LMTD is logarithmic mean
temperature out
temperature difference (°C).
Specific heat of kJ/kg
3
steam ⁰C
2.281 2.905 Energy Balance:
Q M C p T1 T2 m c p t2 t1
Steam mass
4 kg/hr 27758.9 38438.3
flow rate (2)
Exhaust gas
5 temperature ⁰C 350 500.5 Q ̇is energy balance (W), M is exhaust mass flow rate
entered (kg/s), m is water mass flow rate (kg/s), C is the specific
Exhaust gas
6 ⁰C 343 464.4 heat of exhaust gas (J/kg °C), and c is the specific heat of
temperature
Specific heat kJ/kg
feed water (J/kg °C). T1 is the temperature of hot fluid
7 1.101 1.135 entering (°C), T2 is the temperature of hot fluid coming out
exhaust gases ⁰C
Exhaust mass (°C), t1 is the temperature of cold fluid entering (°C), and t2
8 kg/hr 494319 494319
flow rate is the temperature cold fluid coming out (°C).
The LMTD can be calculated by equation, following: OD is the outside diameter of the tube (cm), and 𝑁𝑓 is
number of fins (per tube).
t LMTD
t2 t1 Then to calculate the projection perimeter (P) using the
t (3) following equation:
ln 2
t1
P 21 2?N f 2 2 1 N f Y (7)
LMTD is Log Mean Temperature difference ( C), ∆t1, o
∆t2 is temperature difference at the cold and hot terminals, where P is projected perimeter (cm/m), Y is fin thickness
respectively (oC) (cm), and L is fin height (cm)
2 A Ao
of the number of tubes in HRSG in one bundle uses the
f (4) equation:
de
projected perimeter Y
Nt (8)
With Af
4
OD ' OD 2 N
2 2
f
12
ST
N is the number of pipes per bundle (dimensionless), Y
(5)
is duct width (m), and ST is the vertical distance between
Ao 1 N f . y OD 12
central points of lines (cm).
(6) With ST OD 2 1 ?Dis tan cebetweenthetips
? of the fins
𝐴𝑓 is the surface of thin fins (both sides) (cm2), 𝐴𝑜 is (9)
bare surface on the outside of finned tube (cm2), OD’ outer F. Calculating the Heat Transfer Coefficient in Pipes (hi)
diameter (outer diameter of the pipe plus fin height) (cm),
Calculate Flow area
h ' fi Af Ao
hf (19)
w (12)
Gt Ai
t
Determining Reynold Numbers in Pipes (Ret) where Ω is Fin effectiveness (dimensionless)
Determination of the Reynold number in the pipe (Ret),
Calculating the Overall Heat Transfer Coefficient (UD)
using the following equation [17]
D Gt
h ' fi h 'i (20)
Ret (13) UD
h ' fi h 'i
D is tube diameter, μ and is the water’s viscosity (kg/
(m. s). The value of UD is the overall heat transfer coefficient
(W/(m2) (oC)).
G. Determine the Heat Transfer Coefficient in Pipes (hi) Calculating Heat Transfer Area (A) [18]
Q
In determining the heat transfer coefficient in the pipe A (21)
(hi), you can use the following equation: U Di tLMTD
1
cp 3
k The value of A is heat transfer area, (m2)
h jH (14)
D k Count the Number of Pipe Bundles (n)
A (22)
h is heat transfer coefficient in the pipe (W/(m2°C)), kw n
is heat conductivity of feed water at average temperatures Ai perbundle
(W/ (m °C)), μ is the viscosity of water at average
temperature (kg/ (m. s), 𝑗𝐻 is the Sieder-Tate and Colburn where n is the number of pipe bundles (dimensionless),
heat-transfer factor in the pipe (dimensionless), and c is Ai per bundle is Ai. Nt. L (m2)
specific heat (J/ (kg °C). Calculates Actual Overall Heat Transfer Coefficient
(UD. act.)
H. Calculating the Exhaust Gas Pass Area Q
U D act (23)
Ai t
s XY Nt OD.48 Nt 2 y.L.N f .48 (15)
With:
where αs is exhaust area passage area, (m2), X is duct height,
Ai n Nt a " L (24)
(m), and Y is duct width, (m).
where n is the number of pipe bundles; a” is surface per
I. Calculating Mass Gas Flue Mass Speed (Gs) linear m, m; L is pipe length (m)
W (16) Calculating Dirt Factor (Rd’)
Gs
as
Rd ' Rd Excess fouling factor Adding tothe outside fouling factor
where Gs is exhaust gas mass velocity, (kg/ (m2.s).
(25)
J. Determine Exhaust Gas Reynolds Numbers (Res) where Rd’ is the Combined dirt factor calculated
(m2 °C)/W), Rd is Combined dirt factor provisions
De Gs
Res (17) (m2 °C)/W).
s With excess fouling factor can be calculated using the
equation:
∆Ps is pressure drop on the duct side (bar), f is flue gas A. Design Results
friction factor (m2/cm2), and s is specific gravity. ST is the We note the parameters of the HRSG design as follows:
vertical distance between the center points of the pipe (cm), Pipe thickness (Table V) HRSG HP components of a
SL is the transverse distance between pipe center points standard pipe for typical HRSG models. This condition
(cm), ∅s is the exhaust gas viscosity ratio, and Lp is helps in increasing pressure dynamics and variation.
exhaust gas path length (m). The size of the path (Lp) is HRSG has many tubes (Table VI) to improve the heat
calculated using the following equation: exchange area with exhaust gas flow. Since heat transfer
in HRSG is conducted by convection, HRSG requires as
n vs
Lp (32) much surface area as possible.
12 The HP HRSG downcomer tube area is equal to 0.865
N. Calculating Pressure Reduction in Pipes (∆Pt) m2 which is higher than HRSG (0.6323 m2) [21]. This
condition is possible to increase the water head in the
ft Gt 2 L n evaporator’s downcomer, and it helps enter the natural
Pt (33) circulation process and avoid backflow problems in the
5.22 1010 D s t
evaporator.
where ∆Pt is pressure drop on the pipe side (bar), f is water This difference is due to the significant pressure drop in
friction factor (m2/cm2), s is specific gravity, ∅t is water the number of tubes in the evaporator.
viscosity ratio, and L is pipe length (m).
Outer diameter of the pipe cm 1.905 1.905 2.54 2.54 2.54 2.54
Inside diameter of the pipe cm 1.4224 1.4224 2.0574 2.1742 2.0574 1.859
Pipe length per line m 8.9977 8.9977 8.9977 8.9977 8.9977 8.9977
Pressure-drop (water for bar 0.0291 0.0345 0.02199 0.32336 0.22295 0.12728
ECO in and out; EVAP
in/steam for EVAP out; SUP
in and out)
Dirt factor 0.00306 0.00355 0.00309 0.00305 0.00349 0.00346
B. Specification of HRSG Dual Pressure Level TABLE VII. COMPARISON OF MINIMUM PIPE THICKNESSES WITH
DESIGNED PIPES
The HRSG analyzed in this study supplies two levels of
pressure, namely high and low pressures, to a steam Design pipe
Minimum pipe
Component thickness
turbine. Gas turbine exhaust gas flows through HRSG, thickness [mm]
[mm]
thus providing the heat energy used for steam production. Superheater Low
HRSG contains two superheaters (one HP SUP and one LP 0.147 3.40
Pressure
SUP), two economizers (one HP ECON and one LP Evaporator Low
0.196 1.83
ECON, and two evaporators (one HP EVAP and one LP Pressure
Economizer Low
EVAP). Pressure
0.196 3.40
Superheater High
TABLE VI. SPECIFICATION OF HRSG DUAL PRESSURE LEVEL 0.104 2.41
Pressure
Evaporator High
0.734 2.41
Parameter Unit Value Pressure
Economizer High
0.734 2.41
Pressure
Output MW 35
C. Feasibility Analysis of the Design Dimension
Flue Gas Flow kg/s 137.3
It is crucial to determine the pipe’s outer diameter
before assessing the thickness using the BWG. It is also
Flue Gas Input Temperature C 500.43
important to calculate the pressure drop and impurity
factors of HRSG components such as economists,
Flue Gas Output Temperature C 146.8 evaporators, and superheaters at high and low pressures.
The reviewed operation’s feasibility is the thickness of the
Feed Water Flow kg/s 18.43 feed pipe.
Feed Water Input Temperature C 60.14 D. Effects of Changes in Flue Gas Temperature on the
Property of Water
Flue Gas From Natural Gas
Fig. 1 illustrates an increase in flue gas mass speed due
Efficiency % 69.48 to a temperature rise. Changes in the mass rate of the
exhaust gas increase the amount of water entering the
HRSG. This condition states that the amount of energy
Table VI indicates that the HRSG design can increase from the flue gas is directly proportional to the evaporated
efficiency by up to 68.48 %, which exceeds the value feed water’s mass. HRSG is a significant component of the
indicated in the design made by Wahyu et al, who designed GTCC that utilizes the temperature of the exhaust gas
the propulsion system using HRG, namely 48.49% [22]. produced by the gas turbine to provide steam. The gas
The optimization of the double pressure steam cycle (SDC) turbine’s combustion temperature is not always constant
for waste heat recovery by Liya et al. is and 19.75 [23]. during operation [8]. This condition changes the exhaust
[22] Wahyu Nirbito, Muhammad Arif Budiyanto, and Robby Muliadi,” Sri Wuryanti was born in Indonesia on March
Performance Analysis of Combined Cycle with Air Breathing 26th, 1965. She received the Ir. degree in
Derivative Gas Turbine, Heat Recovery Steam Generator, and Chemical Engineering from Institut Teknologi
Steam Turbine as LNG Tanker Main Engine Propulsion System,” Sepuluh November (ITS), Surabaya, Indonesia
J. Mar. Sci. Eng., vol. 8, no 726, pp 1-15, 2020. in 1990, and M.Sc. degree in Chemical
[23] Recovery Liya Ren and Huaixin Wang, “Parametric Optimization Engineering from Institut Teknologi Bandung
and Thermodynamic Performance Comparison of Organic Trans- (ITB), Bandung, Indonesia in 1994. She
Critical Cycle, Steam Flash Cycle, and Steam Dual-Pressure Cycle received her doctoral degree at the Department
for Waste Heat, “Energies, vol. 12, pp. 1-22, 2019 of Physics, Universitas Indonesia (UI),
Indonesia, in 2015. From 1990 until now, she has been a Lecturer in
Energy Conversion Engineering, Politeknik Negeri Bandung, Bandung,
Indonesia. Email: sriwuryanti.lamda@gmail.com
Copyright © 2021 by the authors. This is an open access article
distributed under the Creative Commons Attribution License (CC BY-
Rakhadian Dwi Jadmiko was born in
NC-ND 4.0), which permits use, distribution and reproduction in any
Indonesia on June 20th, 1995. He received the
medium, provided that the article is properly cited, the use is non-
SST degree in Power Generation Technology
commercial and no modifications or adaptations are made.
from Politeknik Negeri Bandung, Bandung,
Indonesia, 2017. Email:
rakhadiandj@gmail.com