FYP Report Final PDF
FYP Report Final PDF
FYP Report Final PDF
RECOVERY BOILER
Reference#: 15/2017
It is to certify that the following students have completed their project “Optimization
of Waste Heat Recovery Boiler” satisfactorily.
Internal Advisor
Mohammad Ehteshamul Haque
Assistant Professor
Department of Mechanical Engineering
NED University of Engg. & Tech.
Projects’ Coordinator
Mr. Masood Ahmed Khan
Department of Mechanical Engineering
NED University of Engg. & Tech.
We would like to take this opportunity to thank Almighty Allah to whom we owe
everything. It is by the will of Almighty that we are finally able to complete this
project.
We would like to express our gratitude, appreciation and special thanks to Mr Masood
Ahmed Khan (project coordinator) for his time to time co-operation.
We would also like to acknowledge with much appreciation the crucial role of the
staff of Mechanical Engineering Department’s computer programming laboratory and
the Digital Library for their extended support, who gave the permission to use all
required software’s and the necessary material to complete the project.
II
ABSTRACT
Waste heat recovery is relatively better technology to improve the process in the plant
by extracting waste energy and utilizing it in the processing and power generation.
Waste heat boilers are the equipment that extract these energies and used for steam
generation which is used in various processes. In cement plant the waste heat is
recovered from Air quenched cooler (AQC) and Suspension Preheater (SP).
The primary of objective is to increase the efficiency of the Waste heat recovery
boilers which includes Boiler analysis, optimization of parameters, custom and design
of boiler.
Boiler analysis includes the energy balance for the all the equipment of the boiler at
actual data and heat transfer across the tubes.
Optimization of parameters includes the effect of Inlet temperature and flow rate of
exhaust gases on the boiler efficiency and life. Custom design of boiler includes the
changing tube materials, and tube area.
Estimation of boiler tube life is also analyzed and optimized through erosion theory
because the erosion is the most important parameter in defining the boiler tube life
III
TABLE OF CONTENTS
1. INTRODUCTION........................................................................................ 9
1.1 WASTE HEAT RECOVERY: .................................................................. 9
1.2 CYCLES FOLLOWED IN WASTE HEAT RECOVERY SYSTEMS .... 9
1.2.1 Simple Rankine cycle: ......................................................................... 10
1.2.2 Organic Rankine cycle: ........................................................................ 10
1.2.3 Kalina cycle: ........................................................................................ 11
1.3 CONCERNED INDUSTRY .................................................................... 11
1.4 BASIC OPERATIONS OF CEMENT PLANT: ..................................... 12
1.4.1 Sources of waste heat in a cement plant: ............................................. 13
1.5 PROJECT & BACKGROUND: .............................................................. 14
1.5.1 Pre-requisites: ...................................................................................... 14
1.6 OBJECTIVES OF PROJECT: ................................................................. 15
1.7 SCOPE OF PROJECT ............................................................................. 15
4
3.2.1 Estimating the overall heat transfer co-efficient: ................................. 32
3.3 EVAPORATOR DESIGN ....................................................................... 36
3.3.1 Design and operational criterion: ......................................................... 36
3.4 EVAPORATOR DESIGN CALCULATIONS: ...................................... 37
3.5 SUPER HEATER DESIGN..................................................................... 44
3.5.1 Design consideration:........................................................................... 44
3.5.2 Design aspects:..................................................................................... 44
3.5.3 Design assumptions: ............................................................................ 45
3.6 SUPERHEATER DESIGN CALCULATIONS: ..................................... 45
3.6.1 OVERALL HEAT TRANSFER COEFFICIENT OF PLAIN TUBES
50
3.6.2 TOTAL SURFACE AREA FOR HEAT TRANSFER ........................ 50
3.6.3 NUMBER OF TUBES WIDE ............................................................. 51
3.6.4 TOTAL AREA FOR MASS FLOW.................................................... 51
3.6.5 TRANSVERSE PITCH ....................................................................... 51
3.6.6 WIDTH OF SUPERHEATER ............................................................. 51
3.7 INDIRECT EFFICIENCY BOILER OLD CONDITION ....................... 52
5
5.2.2 TOTAL VOLUME .............................................................................. 77
5.2.3 MASS REQUIRED BY DIFFERENT GRADES ............................... 77
5.2.4 COST FOR EACH STEEL GRADE ................................................... 78
Conclusion ……………………………………………………………….. 90
References ………………………………………………………………. 91
6
LIST OF FIGURES
7
LIST OF TABLES
Table 1-1: High Temperature Source for Waste Heat Recovery ................................... 9
Table 2-1: Weight Fraction and Molecular Weight of Flue gas constituents ............. 17
Table 2-2: Viscosity of Flue gas constituents at different temperatures ..................... 17
Table 3-1: Design specifications for Economizer ........................................................ 26
Table 3-0-3: Standard dimension for steel tubes ......................................................... 28
Table 3-4: Maximum Velocity Ranges ......................................................................... 30
Table 3-5: Tube Specifications for Economizer........................................................... 30
Table 3-7: Design specifications for Evaporator ........................................................ 37
Table 3-10: Design specifications for Super-heater .................................................... 46
Table 4-3: Heat Loss across the Outside Wall Vs Inlet Gas Temperature .................. 57
Table 4-4: Efficiency VS Temperature ......................................................................... 57
Table 5-1: Chemical composition of ASTM A36 ......................................................... 63
Table 5-2: Mechanical properties of ASTM A36 ......................................................... 63
Table 5-3: Chemical composition of AISI 1040 ........................................................... 64
Table 5-4: Mechanical properties of AISI 1040 .......................................................... 64
Table 5-5: Chemical composition of AISI 1010 ........................................................... 64
Table 5-6: Mechanical properties of AISI 1010 .......................................................... 65
Table 5-7: Chemical composition of Stainless Steel Grade 410.................................. 65
Table 5-8: Mechanical properties of Stainless Steel Grade 410 ................................. 65
Table 5-9: Chemical composition of Stainless Steel 316 ............................................. 66
Table 5-10: Mechanical properties of Stainless Steel 316 .......................................... 66
Table 5-11: Thermal conductivities of different steel grades ...................................... 67
Table 5-12: Overall Heat Transfer coefficient and Surface Area of economizer for
different materials ........................................................................................................ 69
Table 5-13: Overall Heat Transfer coefficient and Surface Area of evaporator for
different material ......................................................................................................... 70
Table 5-14: Overall Heat Transfer coefficient and Surfaces Area of Super-heater for
different material ......................................................................................................... 72
Table 5-15: Ranges of Price per Ton for different Steel Grades ................................. 75
Table 5-16: Mass required by boiler tubes for different material ............................... 78
Table 5-17: Cost of each steel grade ........................................................................... 78
8
CHAPTER #1
1. INTRODUCTION
9
1.2.1 Simple Rankine cycle:
There are four components:
10
1.2.3 Kalina cycle:
A recent development in power generation technology is the Kalina cycle, which
basically follows the Rankine cycle concept except that the working fluid is 70%
ammonia–water mixture. It has the potential to be 10–15% more efficient than the
Rankine cycle and uses conventional materials of construction, making the
technology viable.
So the waste heat recovery boilers (WHRB) are efficient but not useful for every
industry.
Industries like:
Chemical industry
Cement industry
Sulphuric acid industry
Coke dry quenching industry etc.
11
We have chosen Cement Industry as our primary focus, reasons are:
Among energy-intensive industries, cement is the one where waste heat
recovery (WHR) has been most developed.
Eligible for thorough learning.
Ease of approach to learn.
1. The raw material required to manufacture cement are limestones and clay. Rocks
extracted from the quarry are routed to cement plant nearby.
2. The minerals from the quarry are routed to grinding plant where they undergo
initial milling before being reduced to fine powder. The raw material (80%
limestone and 20% clay) stored in pre homogenization pile. This mixture is called
the “raw meal”.
3. The raw mix is fed into pre-heating tower at 800 °C before returning to inclined
rotary kiln where it is heated to about 1450 °C. Combustion process causes de-
carbonization of limestone. The fired material take the form of hard granules
called “clinkers”.
4. Following the cooling, the clinker is stored in the silos then transformed into
cement according to the production requirement. At final stage, for increasing the
settling time gypsum is added to the clinker in proportion of 3-5% and mixture is
finely ground.
5. The cement is stored in silos before being delivered in bulk using trucks or packed
into for the shipment.
Blasting
Crushing
Grinding
Pre-heating (suspension pre-heater)
Combustion
Cooling (Air Quenched Cooler)
Fine grinding
Packaging
12
Figure 1-4: Schematic of Cement Plant
Afterwards, the clinker is send to grate cooler, where clinker is spread on perforated
moving bed where air is blown from the bottom to cool the clinker. This process is air
quenching, also called Air Quenched Cooling.
Therefore, the exhaust gases leaving the suspension pre heater and the hot air by the
quenching are the two major sources of heat recovery.
13
1.5 PROJECT & BACKGROUND:
The energy crisis is the concern that the world & demand on the limited natural
resources that are used to power the industrial sector as resources are limited and
demand is continuously increasing.
Natural resources may take hundred/ thousands of years to emerge out again to
replenish the stores.
Causes of energy crisis are;
Poor infrastructure of equipment.
Over – population.
Delay in commissioning of power plants.
Wastage of energy.
1.5.1 Pre-requisites:
To accomplish the project in a right manner, there are several things to be fulfilled.
14
1.6 OBJECTIVES OF PROJECT:
Performance evaluation and optimization of the waste heat recovery boiler by varying
different parameters.
15
CHAPTER # 2
It can also be seen in the design calculations that we are interested in finding the
heating surface area of the tubes. Using formula:
And the formula for both Reynold and Nusselt we will be using are:
And,
Finally to calculate :
16
- Find we need (viscosity)
- Find , we need (Prandtl Number)
- Find , we need (Thermal Conductivity)
MOLECULAR
CONSTITUENTS WEIGHT FRACTION
WEIGHT
CO2 0.22 44
H2 O 0.06 18
SO2 0.015 64
O2 0.02 32
NO2 0.67 46
Table 2-1: Weight Fraction and Molecular Weight of Flue gas constituents
2.2 CALCULATIONS:
2.2.1 Viscosity:
Temperature (℃)
μ of
380 300 245.5 219
CO2 0.1079 0.0953 0.0883 0.848
H2 O 0.087 0.074 0.067 0.0641
SO2 0.1002 0.0868 0.0796 0.076
O2 0.133 0.1207 0.1135 0.11
NO2 0.1264 0.1138 0.0891 0.1
Table 2-2: Viscosity of Flue gas constituents at different temperatures
17
- FOR TEMPERATURE = 380 ˚C
Gases . . .
(weight
(each)
frac. each)
Gases . . .
(weight
(each)
frac. each)
18
- FOR TEMPERATURE = 245.5 ˚C
Gases . . .
(weight (each)
frac. each)
CO2 0.22 0.0883 6.633249 0.12886 1.45931478
H 2O 0.06 0.067 4.24264 0.01706 0.2545584
SO2 0.015 0.0796 8 0.00955 0.12
O2 0.02 0.1135 5.65685 0.01284 0.113137
NO2 0.67 0.0891 6.78233 0.40488 4.5441611
Formula:
Cp =
Cp=
19
FOR TEMP = 380 ˚C
Cp * weight
GASES Cp@380 kcal/kg ˚C weight fraction
fraction
CO2 0.2347 0.22 0.051634
H2O 0.4652 0.06 0.027912
SO2 0.1689 0.02 0.003378
O2 0.2301 0.03 0.006903
NO2 0.236 0.67 0.15812
20
FOR TEMP = 219 ˚C
Cp * weight
GASES Cp@219 kcal/kg˚C weight fraction
fraction
CO2 0.22 0.22 0.0484
H2O 0.4545 0.06 0.02727
SO2 0.1599 0.02 0.003198
O2 0.224 0.03 0.00672
NO2 0.2312 0.67 0.154904
Formula:
K=
K=
21
FOR TEMP = 245.5 ˚C
Gases k@245.5 k.cal/m.hr.˚C
CO2 0.0287
H2O 0.0333
SO2 0.0198
O2 0.0337
NO2 0.0321
K=
FOR TEMP = 219 ˚C
Gases k@219 k.cal/m.hr.˚C
CO2 0.028
H2O 0.0315
SO2 0.085
O2 0.0339
NO2 0.035
K=
PRANDTL NUMBER =
22
FOR TEMP = 300 ˚C
(Sp.Heat (Thermal (Prandtl
(Viscosity)
Cap.) Cond.) Number)
1024.0000 2.994 E – 05 0.043 0.713
23
CHAPTER # 3
Traditionally, the heat exchanger performance analysis and simulation are performed
using steady-state energy balance across the heat exchanger. The energy balance on
the hot and cold fluids together with the heat-transfer equation constitutes the model
of heat exchangers. A simplified model generally uses an average driving force such
as log mean temperature difference (LMTD) and assumes uniform properties of the
fluids along the length of the heat exchanger to determine the overall heat-transfer
coefficient. Under the assumption that there is no heat loss to the surroundings, the
heat lost by the hot fluid stream shall be equal to the heat gained by the cold fluid
stream.
Waste heat recovered in cement plant is gases with high dust content, which can be as
high as 150 g/m3. Dust on the pre-heater side is sticky, and dust at the clinker cooler
exhaust is abrasive. This aspect has an impact on the design of heat exchangers used
for waste heat recovery. In addition, sulfur contained in those gases can condensate,
which limits the amount of waste heat that can be recovered.
Figure 3-1:
Rough
Schematic
of Boiler
24
1. Its function is to preheat the feed water before it enters in the evaporator. Only
sensible heating take place, no phase change process otherwise steaming can
cause blockage of tubes and may disturbs the circulation.
2. Therefore required heat transfer in economizer is lower than the saturation
temperature of water at given pressure.
3. But the decision is also dependent on mass flow rate of water. Large quantity of
water can cause steaming too; therefore we kept a margin of 15 to 20 degrees
below saturation temperature.
25
3.1.4 Considerations and assumptions:
1. The operating pressure of the boiler and the degree of superheat of steam
should be fixed.
2. Normally the pressure of steam is high for power generation.
3. The operation pressure should be around 13 bars.
4. The pinch and approach point should be selected; these greatly affect the area
for heat exchange and stability of operation.
5. Suitable temperature of feed water is set at the inlet of economizer.
6. The log-mean temperature difference should be evaluated.
7. The gas-mass velocity should be fixed (depending on the nature of gas i-e
clean or dirty).
8. Dew point should be estimated.
9. Transport properties should be evaluated.
10. The area for the flow or the number of tubes required for a fixed size tube and
velocity can be determined by the equation ṁ= . Here A is the water side
flow area.
11. The gas side area is determined by the equation:
26
The outlet temperature of water is taken short of the saturation temperature at
operating pressure of 13 bars, in order to prevent steaming.
The values of enthalpies are taken from book of “Heat Transfer” by Yunus A. Cengel.
= 2.222(741.5 – 167.57)
= 1274.124 kW
1274.124 kW = ṁ ∆T
= 306.83 °C
= = 314.385°C
27
∆T LMTD =
Where
∆ Ta = –
∆T b =
321.94°C
306.83°C
40°C
175°C
∆T LMTD =
∆T LMTD = 200.96 °C
28
compensate fouling and corrosion and also reduces pressure drop.
Thickness: t =3.175 mm =0.003175 m ( )
Cross-sectionalqareaqforqoneqtube:
Theqcrossqsectionalqareaqofqtheqtube would be used in the formula of mass flow
rate to find out the number of tubes width and hence width of the boiler.
A= = (0.01905)2
A = 2.85× m2
Number of tubes
Continuity equation tells us that in how many tubes does a mass flow rate of water
(ṁ) of 2.22 kg/sec would be divided.
Nw=
Where,
Since,
= = = 937 kg/m3
ṁ=2.22 kg/sec
A = 2.85× m2
For velocity of the fluid iterations are required within the suitable range as given in
the table below
29
Table 3-3: Maximum Velocity Ranges
= 16.626 ≈ 17 tubes
Thus 7 tubes are quite reasonable for lesser frontal area. This would curtail the length
of the tubes as given by the following relation =
Tubes specs for economizer:
thickness 3.175 mm
number of tubes 17
area m2
Table 3-4: Tube Specifications for Economizer
Fixing the gas mass velocity:
Gas mass velocity is basically the mass flux (amount of exhaust gas passing through
per unit area). The range for the gas flux typically ranges from 6-16 kg/m2sec,
depending upon the type of gas which ascends from dirtiest to clean.
ṁ = ρAV
= ρV =G
30
Where is the frontal area, re-arranging the equation and putting values we get
= 9.278 m2
Then;
ST = + 0.0254
31
3.2.1 Estimating the overall heat transfer co-efficient:
Overall-heat-transfer coefficient is useful for finding the surfaceqareaqforqheat
transfer. Since convective heat-transfer-coefficient of gas side has much lower value
than water side, therefore gas side coefficient is the governing parameter. Even
though they don’t have much impact, but for higher accuracy we have considered it.
Also we have to consider the fouling factor on gas side and corrosion factor on water
side. Therefore the equation for overall coefficient becomes
32
For water side:
= 0.0254m
( ) H₂O = 107.5°C = 225.5°F
Reynolds number:
Re = =
Re =
Re = 27546.06
Nusselt number:
As we know that critical Reynold number for internal flows is Re 2300, here
Reynold number is greater than this, hence the flow is turbulent.
Now using Dittus-Boelter correlation for finding Nusselt number (Book: Incropera)
NU = 0.023
0.8
NU = 0.023 ( (2.02)0.4
NU = 108.6215
hi =
h = 3.843 kW/m2K
33
For Gas side:
Properties of gases evaluated at the film temperature
= 314.385
= 1.010 kJ/kg
Nusselt number:
Since Reynold number for flow over tube bank is less than the critical Reynold
numbers. For outside heat transfer coefficient using correlation
Nu = 0.33 .
This relation is known as GRIMSON MODEL and is used for flow over tube banks
when
Re 4
Pr
and m are constants. For values of these, two tables are given below
34
Using tube material to be low Carbon steel thermal conductivity k = 50 W/m.K
Now;
U = 90.36635 W/m2K
Flow area required:
Now calculating the surface area:
A = 70.1605 m2
Number of tubes high required:
= ×
= ×
= 13.61 13
35
3.3 EVAPORATOR DESIGN
36
16. In case of clean surface, gas film resistance is important and considerable as
compared to steam or water film resistance or tube wall conductance.
Prandtl Number Pr
Table 3-5: Design specifications for Evaporator
37
Amount of heat to be generated:
This heat will be equal to the heat gained by water:
=
= (2.22 2022)
= 4488 kW
Applying Energy Balance:
Energy gained by the water = Energy lost by flue gas
= 321 C
Where,
=
After Rearranging,
38
Where,
Here,
For and :
We will follow the following steps, first we will calculate Nusselt Number, then
Reynolds and then Convective Heat Transfer Coefficient.
39
For gas side:
:
Here,
……. [Gas mass velocity is nothing but the mass
flux ]
----------- (1)
Nu = 1.13 .
Nu = 1.13 .
According to the calculated Reynold Number Figure 3-7: Constants for Nusselt
– C = 0.229 and m = 0.632 correlation
Hence, Nu = 102.3
Now,
40
For water side:
:
Here,
Solving, we get,
Nu= 0.023
Nu= 1529.0108
Now,
We get,
41
For gas side:
Also,
Here,
……… (B)
Also we have,
A = Total required surface area taking all the tubes.
A= ----------------- (C)
Putting from eq (B):
42
For water side:
Now finding the velocity of fluid (water) in evaporator
Again using equation (1)
=
Where,
After rearranging,
-------------------- (D)
Where,
= Velocity of fluid in evaporator
--------------- (E)
Here,
Also, C. R = 4,
Substituting, we get
= 244.59 kg/
= 0.0312
43
3.5 SUPER HEATER DESIGN
1. The production of steam at higherqtemperatureqthanqtheqsaturationqtemperature
is called superheating.
2. The temperature added to the saturationqtemperatureqisqcalledqtheqdegreeqof
superheat.
3. A superheater is employed in the steam boiler to add additional energy into steam
and raise its temperature.
4. Superheater is in demanding condition, and their failure could mean shut down
for few days resulting in large financial losses.
5. The function of superheater is to increase the capacity of the steam boiler,
eliminates erosion of the steam turbine, and reduces steam consumption of the
steam turbine.
1. Uniform distribution of steam and gas across all the sections to minimize
imbalance of flows.
2. Provision of thermal expansion of headers, tubes, spacers and supports.
3. Accessibility of cleaning, examination and removal of elements
4. Optimally high steam velocity in all the tubes to keep the metal temperatures as
low as possible
5. Minimum steam pressure losses.
44
3.5.3 Design assumptions:
Velocity:
Superheaters transfer heat from flue gas to steam. Heat transfer between two gases is
not very effective compared to heat transfer from gas to fluid. For that reason, steam
must flow fast enough (10-20 m/s) in order to give the superheater tubes enough
cooling. Lower steam pressure weakens the heat transfer rate, so with lower pressures,
steam must have a greater velocity (15-40 m/s).
Spacing:
Superheater of boiler consists of banks of tubes. A system of tubes is located in the
path of the furnace gases in the top of furnace. A superheater must be built so that it
superheats approximately the same amount of steam from low to high loads.
Changing tube lengths between passes can control temperature differences. The
outermost tube that receives the most radiative flux should be shorter than the rest of
the tubes.
Tubes:
Tubes in superheaters can be arranged according to inline or staggered arrangement.
Inlineqtubeqarrangementqis preferred forqfoulingqboilers, andqrecovery. Staggered
arrangementqisqpreferredqforqoil, gasqandqheatqrecoveryqsteam generator. As free
space with staggered arrangement is much smaller than with inline arrangement the
reason for decreased fouling with inline is evident.
45
Outlet steam temperature Tso 294.28°C
Thermalqconductivity K 0.0364qW/m.K
E gain = 457.32 kW
Temperature for Transport property:
Average temperature of steam:
Where
0.95 is for fouling effects
0.98 is for radiation losses
374.34
46
Average temperature of flue gas:
377.17
Where,
∆ = out,– sin,
∆ = , – s,
, 380
, 373.4°C
s, 205.7°C
s, 294.28°C
Area =
Area = 12.568 m2
Assuming volumetric flow rate of flue gases = = 166.3 m/s
47
Finding the velocity of flue gases:
Vf =
Vf = 13.23 m/s
Reyonld number =
At
ρ = 0.5 kg/m3
V = Vf = 13.23 m/s
D = Do = 0.0508 m
Re = 10000
Finding Reynold number at particular gas mass velocity
At G = 12 kg/m2 sec
48
Figure 3-9: Table of constant for Nusselt correlation
C2 selected on the basis of number of tubes C2 (8 rows) = 0.96
49
Specific heat of steam s 3.37 kJ/kg K
= 1300 W/m.K
Here,
=
=
= (steam ) = .0002 m2.oC/W
= ( ) = .001 m2.oC/W
When,
ho = 82.98 W/m.K
= 1300 W/m.K
50
TRANSFER
Assuming Nh 4
Nw = 22.06
A = 0.032 m2
From the definition of continuity equation, velocity of steam is given by
V= 12.327 m/s
W = 2.95 m
51
3.7 INDIRECT EFFICIENCY BOILER OLD CONDITION
In order to find the current condition of boiler, we use indirect method of efficiency
calculation. The efficiency can be measured easily by measuring all the losses
occurring in the boilers using the principles to be described. The disadvantages of the
direct method can be overcome by this method, which calculates the various heat
losses associated with boiler. The efficiency can be arrived at, by subtracting the heat
loss fractions from 100.An important advantage of this method is that the errors in
measurement do not make significant change in efficiency.
52
Dry flue gas loss:
Radiation losses
2%
Unaccountable loses
2%
Manufacturer’s margin
Therefore
53
CHAPTER # 4
There are two ways to increase the efficiency of boiler which include either to change
the design of the boiler or its auxiliaries and second one is to optimize the boiler
parameters to obtain the best efficiency of the boiler. There are various parameters of
the boiler which can be controlled. Following are the variables that may change
without altering the design of plant and boiler to obtain the optimization.
1. Inlet gas temperature
2. Inlet mass flow rate of gas
3. Feed water flow
54
And the effectiveness of each component of the boiler is shown below the following
graph showing the increase in effectiveness of the boiler as a function of inlet gas
temperature.
55
But due to high temperature of gases there will be an overall increase in efficiency of
the boiler as the heat transfer enhance more than the losses. The graph between
efficiency and the temperature clearly showing the efficiency so make the temperature
about 385oC with constant mass flow of the gases.
So, these are the advantages of the increase in inlet gas temperature.
Further if we increase the temperature of the inlet gas the feed water flow is adjusted
accordingly and it will increase the overall power output of the turbine but there are
several disadvantages of temperature increment on the gas side which includes
1. Reduction of the tubes life due to increase in velocity of flue gas which will
enhance the erosion rate discuss in the next chapters.
2. Heat transfer across the outside wall will also increase significantly led to
increment of the Heat loss from 1.2% to 1.3%. Following graph shows the loss of
the energy as the inlet gas temperature increases.
3. Increase in creep which is another factor of life reduction.
4. Increase in boiler temperature also increases the irreversibility’s of the cycle and
boiler by increasing the pinch point of the system and decrease the efficiency of
the boiler.
Pinch and approach points are important parameters which determine directly the
steam generation rate, the temperature profiles, and the heat transfer surface area
required in the WHRB.
PINCH POINT: The pinch point is defined as the minimum temperature difference
between two streams, in this case the exhaust gas and the water or steam. In practice,
the pinch point is usually the difference between saturation temperature and
temperature of the exhaust gas at the evaporator outlet in waste heat recovery steam
boilers. After the exhaust gases have transmitted a certain amount of heat and cooled
to a certain temperature, there is a point below where the further cooling of the gases
requires disproportionately large heat transfer surface area resulting in expensive
WHRB design. As the pinch point decreases, the increase of the heat transfer surface
area of the evaporator is exponential, whereas the increase of recovered heat is only
linear. In other words, the smaller the pinch point, more efficient the WHRB, but also
more expensive the design. The commonly used pinch point in WHRB designs is
usually between 8 to 15 °C.
56
Table 4-1: Heat Loss across the Outside Wall Vs Inlet Gas Temperature
57
Figure 4-3: Mass flow rate of gas Vs Overall Heat Transfer
coefficient for Economizer, Evaporator and Superheater
58
increase in temperature which includes No heat transfer increment across the wall of
the boiler and No growth in creeping phenomenon of the tubes.
But the only disadvantage of flow rate increment is the increase in velocity of the
gases and gases has significant amount of kiln dust content which will enhance the
erosion of the wall exponentialy .The complete relation of the velocity and erosion
and discuss in the chapter of erosion. So it should be 45 to 55 kg/s when operating at
about 350oC to 380 oC.
Itqisqimportantqtoqnoteqthatqincreaseqin mass flow will increase the pressure losses
across the tubes which will lead to the requirement of larger Induced Draft (ID) fans
operating at its maximum power for obtaining the maximum efficiency of the boiler.
59
4.3 Constant Inlet energy:
If the inlet energy is constant and the inlet temperature of the gases is increased then
their will be reduction of mass-flow-rate of the gas.The following graph is showing
the detail of the decrease in effectiveness with increase in inlet temperature of the
gases at constant energy because of the mass flow rate reduction which will reduce
the effectiveness and so the efficiency.
From the graph it is clear that the effectiveness is reduced and it is recommended that
one should increase the mass flow rate and reduce the temperature if constant amount
of energy is transferred to the boiler. However the erosion rate will also increase with
the increase in mass flow significantly and reduce the boiler life.
Figure 4-5: Mass flow rate Vs Effectiveness & Inlet Temperature Vs Effectiveness for
Economizer
Figure 4-6: Mass flow rate Vs Effectiveness & Inlet Gas Temperature Vs
Effectiveness for Evaporator
60
Figure 4-7: Mass flow rate Vs Effectiveness & Inlet Temperature Vs Effectiveness
From the graph it is clear that the effectiveness is reduced and it is recommended that
one should increase the mass flow rate and reduce the temperature if constant amount
of energy is transferred to the boiler. However the erosion rate will also increase with
the increase in mass flow significantly and reduce the boiler life.
4.4 Conclusion:
So from the results obtain in the form of the graphs showing above it can be
concluded that:
1. Temperature should be greater than 370oC and less than 450oC provided that
the mass flow rate in constant to obtain the maximum efficiency.
2. Flow rate increase the efficiency and effectiveness of boiler greater than the
inlet temperature of the gases.
3. All the efficiency increment is at the cost of boiler life.
61
CHAPTER # 5
5. MATERIAL OPTIMIZATION:
A general procedure that could be used for identifying the most appropriate material
for a specific heat exchanger application would consist of the following steps:
1. Defineqtheqheatqexchangerqrequirements.
2. Establishqaqstrategyqforqevaluatingqcandidateqmaterials.
3. Identifyqcandidateqmaterials.
4. Evaluateqmaterialsqinqdepth.
5. Selectqtheqoptimumqmaterial.
62
Mechanical properties:
1. Good tensile and creep properties (High creep rupture strength at the highest
temperature of operation and adequate creep ductility to accommodate localized
strain at notches are important).
2. Good fatigue, corrosion fatigue and creep-fatigue behavior.
3. High fracture toughness and impact strength to avoid fast fracture.
Corrosion resistance:
1. Low corrosion rate to minimize the corrosion allowance.
2. Resistance to corrosion from off normal chemistry resulting from leak in upstream
heat exchanger or failure in the chemistry control
3. Tolerance to chemistry resulting from mix up of shell and tube fluids.
ASTM A36:
ASTMqA36qisqtheqmostqcommonlyqusedqmildqandqhotqrolledqsteel.qItqhas
excellentqweldingqpropertiesqandqisqsuitableqforqgrinding,qpunching,qandqtapping,
drilling and machining processes. Chemical composition is given below.
Carbon 0.25-0.29%
Copper 0.20%
Iron 98.0%
Manganese 1.03%
Phosphorus 0.04%
Silicon 0.28%
Sulfur 0.05%
Table 5-1: Chemical composition of ASTM A36
Mechanical properties:
Ultimate tensileqstrength 400-550qMPa
Yield strength 250qMPa
Elasticqmodulus 200qGPa
Bulkqmodulus 140qGPa
Shearqmodulus 79.3qGPa
Poissonqratio 0.26
Table 5-2: Mechanical properties of ASTM A36
Machining:
The machinabilityqrateqofqASTMqA36qisqestimatedqtoqbeq72%, andqtheqaverage
surfaceqcuttingqfeedqofqASTMqA36qisq120 ft. /min. Machining of ASTM A36 steel
is not as easy as that of AISI 1018 steel.
63
Welding:
ASTM A36 steel is easy to weld using any type of welding methods, and the welds
and joints so formed are of excellent quality.
AISI-1040
AISIq1040qcarbonqsteelqhasqhighqcarbonqcontentqandqcanqbeqhardenedqbyqheat
treatment followed by quenching and tempering to achieve 150 to 250 ksi tensile
strength. Chemical composition is given below.
Carbon 0.37-0.44%
Iron 98.6-99.0%
Manganese 0.60-0.90%
Phosphorus ≤ 0.04%
Sulfur ≤ 0.05%
Table 5-3: Chemical composition of AISI 1040
Mechanical properties:
Ultimateqtensileqstrength 620qMPa
Yieldqstrength 415qMPa
Elasticqmodulus 190-210qGPa
Bulkqmodulus 140qGPa
Shearqmodulus 80qGPa
Poissonqratio 0.27-0.30
Table 5-4: Mechanical properties of AISI 1040
Machining:
The machinability rating of AISI 1040 carbon steel is 60.
Weldability:
AISI 1040 carbon steel can be welded using all welding techniques. It can be
preheated at 149 to 260°C (300 to 500°F) and post heated at 594 to 649°C (1100 to
1200°F) due to its high carbon content.
AISI-1010
AISI 1010qcarbonqsteelqisqaqplainqcarbonqsteelqwithq0.10%qcarbonqcontent. This
steel has relatively low strength but it can be quenched and tempered to increase
strength. Chemical composition is given as:
Carbon 0.080%
Manganese 0.30%
Phosphorus 0.05%
Sulfur 0.04%
Iron 99.18%
Table 5-5: Chemical composition of AISI 1010
64
Mechanical properties:
Machinability:
TheqmachinabilityqofqAISIq1010qcarbonqsteel, especiallyqinqtheqcoldqdrawnqor
cold worked state, is considered as fairly good.
Welding:
AISI 1010 carbon steel can be welded using all the conventional welding techniques.
STAINLESS STEEL GRADE-410:
Gradeq410qstainlessqsteelsqareqgeneralqpurposeqmartensiticqstainlessqsteels
containingq11.5%qchromium,qwhichqprovideqgoodqcorrosionqresistanceqproperties
. However, the corrosion resistance of grade 410 steels can be further enhanced by a
series of processes such as hardening, tempering and polishing. Quenching and
tempering can harden grade 410 steels. Chemical composition is given as:
Carbon 0.15%
Manganese 1.00%
Silicon 1.00%
Phosphorus 0.04%
Sulfur 0.03%
Chromium 13.00%
Nickel 0.75%
Iron 84.03%
Table 5-7: Chemical composition of Stainless Steel Grade 410
Mechanical properties:
Gradeq410qsteelsqcanqbeqweldedqusingqallqconventionalqweldingqtechniques,qbut
theqmaterialsqshouldqpreheatedqat 150 to 260°C followedqbyqpostweldqannealing
treatment,qtoqmitigateqcracking. According to AS 1554.6 standards, grade 309
electrodes or rods are preferred for welding 410 steels.
65
Machining:
Grade 410 steels can be easily machined in highly tempered or annealed conditions.
However, it is hard to machine grade 410 steels if they are hardened above 30HRC.
Free machining grade 416 is the best alternative.
Excellent weldability by all standard fusion methods, both with and without filler
metals. AS 1554.6 prequalifies welding of 316 with Grade 316 and 316L with Grade
316L rods or electrodes (or their high silicon equivalents). Heavy welded sections in
Grade 316 require postweld annealing for maximum corrosion resistance. This is not
required for 316L. Grade 316Ti may also be used as an alternative to 316 for heavy
section welding.
Machinability:
66
5.1.2 TREND OF THERMAL CONDUCTIVITY ON ‘U’ VALUE:
Different grades of AISI and ASTM standards are considered whose composition and
properties are mentioned in the above section. Now thermal conductivities of each
grade are considered. The table below shows thermal conductivities:
Where,
5.1.4.1 ECONOMIZER:
In order to calculate the overall-heat-transfer-coefficient as the function of thermal
conductivities, all other values are needed to be fixed. From the design section the
values are given below:
67
K/W
K/W
W/ .K
W/ .K
Thermal conductivity of mild steel is 42 W/m.K at elevated temperatures. Now,
+ + In +0.001+0.0002
4.64×
1.135
U 89.07
= 71.16
In a similar manner as above, the values of surface area and overall heat transfer
coefficient is calculated. This work is done directly on Microsoft Excel sheet. The
values are as follows:
Thermal
Material U A
Cond (k)
68
Economizer
5.1.4.2 EVAPORATOR:
In order to calculate the overall-heat-transfer-coefficient as the function of thermal
conductivities, all other values are needed to be fixed. From the design section the
values are given below:
K/W
K/W
W/ .K
W/ .K
Thermal conductivity of mild steel is 42 W/m.K at elevated temperatures. Now,
+ + In +0.001+0.0002
2.36×
69
1.2983
U 75.0236
= 369.55
In a similar manner as above, the values of surface area and overall heat transfer
coefficient is calculated. This work is done directly on Microsoft Excel sheet. The
values are as follows:
Thermal
Material U A
cond. (k)
Evaporator
70
5.1.4.3 SUPERHEATER:
In order to calculate the overall-heat-transfer-coefficient as the function of thermal
conductivities, all other values are needed to be fixed. From the design section the
values are given below:
K/W
K/W
W/ .K
W/ .K
Thermal conductivity of mild steel is 42 W/m.K at elevated temperatures. Now,
+ + In +0.001+0.0002
2.31×
1.511
U 66.6
= 64.57
71
U-values of other grades:
In a similar manner as above, the values of surface area and overall heat transfer
coefficient is calculated. This work is done directly on Microsoft Excel sheet. The
values are as follows:
Thermal
Material U A
cond. (k)
Superheater
5.1.5 CONCLUSIONS:
From the above table you can observe the trend that thermal conductivity has on the
overall-heat-transfer-coefficient and surface area required for heat transfer. As
thermal conductivity increases, overall-heat-transfer-coefficient increases, therefore
lesser area will be required for same amount of energy transfer.
Conversely we can also say that if material of greater thermal conductivity is used,
then for same surface area amount of energy transfer would be maximized.
We will be showing this result through calculations. For this purpose we have to
redesign the boiler.
72
5.1.6 EFFECT ON FINAL TEMPERATURES OF COLD FLUID:
For the same surface area and new overall-heat-transfer-coefficient the amount of
power ( ) would also be maximized. This will give us a rise in the outlet
temperatures of water at every section.
After material shifted from mild steel to high thermal conductivity material AISI
1010, then U values of every section becomes
= 92.92 W/ .K
= 77.61 W/ .K
= 77.81 W/ .K
Now,
=
= 1290.87
= 4525.02
= 469.50
ECONOMIZER
For the changed amount of heat rate the outlet temperature in economizer is given as
138.35
= 178.35
73
EVAPORATOR
For the previous temperature difference and phase change, evaporator requires 4488
KW of heat rate. Now here are two changes:
1. Inlet temperature of evaporator changes.
2. Material change will further increase the heat rate and outlet temperature.
Now
At 178.35 =C = 440 = 784.74
At 191.5 = 2792
After material changed, the heat rate available has been calculated and is found to be
greater than the previous one.
= 4525.02
Therefore extra amount of energy will be utilized in increasing the outlet temperature
of saturated steam.
= -
68.9
This extra energy will increase the outlet temperature of saturated steam coming out
of evaporator section.
= = 31.04
= = 2824
Now from steam table at pressure 13 bars superheated section across enthalpy, final
outlet temperature is found through interpolation.
= 205.7
74
SUPERHEATER
For the changed amount of heat rate the outlet temperature in economizer is given as:
=
469.5 = 2.22
246.38
= = 3071
Now from steam table at pressure 13 bars superheated section across enthalpy, final
outlet temperature is found through interpolation.
= 313
75
5.2.1 VOLUME OF MATERIAL USED IN HEAT EXCHANGERS
PIPING:
5.2.1.1 ECONOMIZER
Tube specification of economizer are:
5.2.1.2 EVAPORATOR
Tube specification of evaporator are:
76
5.2.1.3 SUPERHEATER
Tube specification of super-heater are:
Also
Since each steel grade has different density, therefore the mass required for each
grade will be different. This is a good way to truly justify our material.
77
MATERIAL DENSITY TOTAL TOTAL MASS
VOLUME
78
CHAPTER # 6
6. MATERIAL INSPECTION:
production of sulfuric acid may be probable. Thus testing our materials in acidic
medium of sulfuric acid will give us right magnitudes of corrosion.
0.5 molar solution of sulfuric acid without corrosion inhibitor conditions were taken
from research articles for mild steel and AISI 1010.
FOR MILD STEEL
For mild steel directly measured corrosion rate is given:
Where,
79
Therefore surface area is
The specimen is tested for 3 days dipped in the solution, therefore time in hours
would be 72 hours
Now,
It is found that the corrosion rate of AISI 1010 is approximately 1.5 times the
corrosion rate of mild steel and comparing the tubes on same dimensions (thickness)
AISI 1010 needs to be changed earlier than mild steel.
On the other hand by using AISI 1010 sufficient increase in the degree of superheat is
observed i.e. 20 degrees increment. So you cannot let it go by seeing only one minor
difference of corrosion estimation.
80
Figure 6-1: T-s diagram of Rankine
Cycle
While our system utilizing mild steel as tubing material the degree of superheat was
When we upgrade our material from mild steel to AISI 1010, our degree of superheat
is found to be:
Just because we neglected the pump work, enthalpy at point 1 and point 2 are same.
For mild steel @ 295
Now,
81
WORK OUTPUT
The difference in superheat in terms of work output is given as:
FOR MILD STEEL
This means that upgrading the material is improving our work output up-to 13.46%.
Therefore changing the material is better.
82
6.2 EROSION IN BOILER TUBES
6.2.1 THEORY
EROSION:
Erosion-is-a-form-of-wear. Wear can be defined as progressive loss of a surface
material due to mechanical action involving impingement of abrasive particles.
Unlike corrosion, which is a chemical or electrochemical action, wear is purely
mechanical.
Abrasion and erosion are two types of wear.
1. Abrasionqcanqbeqlikenedqtoqsand-paperingqinqwhichqsolidqparticlesqmoveqin
contactqwithqaqparallelqsurface. Abrasionqaffectsqtheqhighqspotsqofqtheqsurface
withoutqmuchqeffectqonqtheqmainqbody. The resulting loss of material is smaller in
comparison with erosion. Abrasion resistance can be built by a boundary layer of high
and preferably hard spots.
2. Erosion is the impingement of hard particles at an inclination, and it has more
energy and destructive power than abrasion. The impinging particles cut through the
boundary layer and destroy the main matrix. Hence, abrasion-resistant material cannot
withstand erosion.
Erosionqofqhotqpartsqdueqtoqash,qparticularlyqtheqtubes,qisqaqseriousqproblem
affecting the life of boiler tubes. The hard constituents moving with flue gas at high
velocities impinge on tubes, refractory, and other parts in the gas path, causing
erosion.
Erosion is influenced by the following:
Erosion is an inseparable. It can be minimized but not eliminated altogether. The aim
is to minimize and predict the erosion rate so that the intervals between the downtimes
83
can be extended to coincide with the planned outages, thus improving the unit
availability. The erosion prevention and protection measures essentially fall into two
categories:
1. Design stage measures
2. Provision of sacrificial protective material
84
= Gas velocity in the narrowest section of tubes
= strike efficiency
85
According to BeckmannqandqGotzmann, the ratio isqaquniversalqparameterqused
Figure: 6-2
Source: Solid Particle Erosion: Occurrence, Prediction & Control by Ilmar Kleis, Priit Kulu
86
= Modulus of Elasticity of Target material
6.2.3 CALCULATION
EROSION RATE AT DIFFERENT IMPINGMENT ANGLE
6.2.3.1 OBJECTIVE:
FORMULA USED:
Erosion by plastic contact:
The reason of using this formula is that it is more general and does not depend on the
type of fuel used, while other approaches are restricted to the erosion caused by fly
ash of coal using as a fuel in fired tube boiler. In our case, the boiler is waste heat
boiler no firing happens only the cause of dust content is the Cement kiln dust (CKD).
Assumptions:
The particles are spherical
Particlesqcausingqerosion is homogeneousqandqelasticallyqdeformable
The velocity of particle before immediate impact is constant
The particles have no rotational energy
Homogenous dust content distribution
Flow rate and dust content are taken at nominal operating point. Changes due
to temperature is not incorporated
DATA:
Density of particle = = 2799.92
87
Hardness of target material (AISI 1010) = 60 (Rockwell B)
Hardness of the Particle (CKD) = 113 (Rockwell B)
ModulusqofqElasticityqof Target material = 190GPa
ModulusqofqElasticityqof Particle = 20GPa
Poisson’s ratio of target material = = 0.27
= 0.07 (from the graph presented by to Beckmann and Gotzmann given in literature
review heading)
1 03
5 20
10 15
20 20
30 16
40 10
50 06
60 03
>60 07
88
= 18.7096
R = 9.3548
= 18.6031028GPa
= 2.47963e-07 m
= 5.97204e-4
=0
= 3.5644e-09
89
6.2.3.4 WEAR RATE ( AT DIFFERENT ANGLE OF IMPINGEMENT:
The velocity is kept constant as well as the hardness, density and viscosity the
resultant effect will be in following figure:
6.2.3.5 RESULT:
Angle of impingement =
90
Poisson’s ratio of particle = = 0.15 (approx.)
It is given the nominal flow rate and dust inlet content. To find the volume of material
removed from the boiler tubes we have just find the dust content mass flow rate and
then simply by multiplying the mass flow rate by the wear rate
The angle of impingement is kept constant i.e. 90 degrees
Angle of impingement =
91
At different wear rate calculated using above formulae at different velocities the graph
between volumes removed per hour with velocity will be:
92
CONCLUSION
We have a boiler operating at degraded performance. Through heat transfer analysis
boiler was defined in terms of number of tubes, heat transfer area and overall heat
transfer coefficient. Our work was on the material of tubes in boiler and optimizing
boiler gas side parameters.
In our boiler material was mild steel that has a thermal conductivity of 42 W/m.K
through this the boiler was giving superheated steam at 295 C. while optimizing gas
sides parameters we considered:
1. Inlet gas temperature effect on overall heat transfer coefficient and effectiveness.
2. Mass flow rate of gas effect on overall heat transfer coefficient.
3. Mass flow rate of gas effect on effectiveness of heat exchanger.
4. Inlet gas temperature effect when constant energy input.
5. Mass flow rate of gas effect when constant energy input.
For inlet gas temperature effects at constant mass flow rate of gas, it is
recommended to use temperature of 380 to 400 degrees not greater than this
because this will lead to lesser tube life due to creep loading. At temperature 385
degrees overall efficiency found to be 85.11%.
For mass flow rate of gas effects at constant inlet temperature, it is recommended
to use 80 to 90 kg/sec otherwise inlet dust content and increase in velocity will
increase, thus increase erosion rate of tubes.
For constant energy inlet, increasing mass flow rate increases the effectiveness,
while increasing the inlet temperature reduces it. The optimum point is achieved
by increasing the mass flow rate and decreasing the inlet temperature. Keeping
erosion rate in mind.
Also by material upgrading to AISI 1010 that has high thermal conductivity, it is
found that degree of superheat enhanced and work output increased by approximately
13%.
Overall optimization can be concluded by considering the life of boiler. Factors like
corrosion, erosion can reduce lifespan and thus increase frequency of overhauling.
Final conclusion is increase the mass flow rate up-to 90 kg/sec for fixed hea0t transfer
ra0te, th0us reducing the high inlet temperature requirement. This will increase the
tubes life by reducing creep load and also increase the overall efficiency of boiler.
93
REFERENCES