Evaporators PDF
Evaporators PDF
Evaporators PDF
(3.9)
where,
1 1 2 2 3 3
( ) ( ) ( )
b w b w b w
BPE T T T T T T = + +
(3.10)
NPTEL Chemical Engineering Chemical Engineering Design - II
Joint initiative of IITs and IISc Funded by MHRD Page 22 of 31
2. Calculate the total amount of solvent vaporized from the feed and product
concentration and feed flow rate. It is assumed that heat transfer rate in each effect
is roughly equal. This signifies that the rates of vaporization in each effect are also
roughly equal. Calculate the approximate vaporization rate in each effect (it is
one-third of total amount of solvent vaporized in one effect in case triple effect
system). Calculate the concentration in each effect and find out the BPE in each
effect. Then calculate the overall temperature drop (
tot
T A ).
3. Redistribute the overall temperature drop (
tot
T A ) among all the effects. Since the
areas are the same (
1 2 3
A A A = = ), the temperature difference in each effect is
roughly proportional to the overall heat transfer coefficients.
1 1 2 2 3 3
1 1
2 1 3 1
2 3
1 1
1 2 3 1
2 3
,
1
D D D
D D
D D
D D
tot
D D
U T U T U T
U U
T T T T
U U
U U
T T T T T
U U
A = A = A
| | | |
A = A A = A
| |
\ . \ .
| |
A = A + A + A = A + +
|
\ .
(3.11)
Thus, calculate
1
T A ,
2
T A and
3
T A .
4. Use the calculated value of
1
T A ,
2
T A and
3
T A ; and the composition estimated to
calculate the enthalpy values. The same reference temperatures for enthalpy must
be used for all streams, including those taken from steam tables, etc. Solve the
enthalpy balance equation sequentially to find out
s
m ,
1 s
m and
2 s
m .
5. Use heat transfer equations to calculate the heat transfer area for each effect.
1 1 2 2
1 2 1
1 1 2 2 3 3
, ,
s s s s s s
D D D
m m m
A A A
U T U T U T
= = =
A A A
(3.12)
6. Compare the areas calculated. If they are not equal, repeat the calculation. Using
the areas obtained to revise the temperature estimates. The recommended
approach is to use the ratio of the calculated heat transfer area for an effect to the
arithmetic mean of the calculated areas as shown:
1
1, 1 new
mean
A
T T
A
| |
A = A
|
\ .
(3.13)
Repeat the calculations until the area of each effect is equal.
NPTEL Chemical Engineering Chemical Engineering Design - II
Joint initiative of IITs and IISc Funded by MHRD Page 23 of 31
The procedure above discussed is applicable for forward feed evaporators. The energy
balance equation can developed accordingly for backward feed system. Here also the
area for each effect is considered to be equal.
The heat transfer area required to obtain the product purity from a specified feed is
illustrated in the next example for forward feed evaporator system.
NPTEL Chemical Engineering Chemical Engineering Design - II
Joint initiative of IITs and IISc Funded by MHRD Page 24 of 31
Lecture 4: Solved Example
Design problem
A 5% aqueous solution of a high molecular weight solute has to be concentrated to
40% in a forward-feed double effect evaporator at the rate of 8000 kg.h
-1
. The feed
temperature is 40C. Saturated steam at 3.5 kg.cm
-2
is available for heating. A
vacuum of 600 mm Hg is maintained in the second effect. Calculate the area
requirements, if calandria of equal area are used. The overall heat transfer coefficients
are 550 and 370 kcal.h
-1
m
-2
C
-1
in the first and the last effect respectively. The
specific heat of the concentrated liquor is 0.87 kcal.kg
-1
C
-1
.
SOLUTION
Part 1. Thermal design
Pressure in effect I to be decided.
Pressure in effect II = 760 - 600 = 160 mm Hg
Boiling point at this pressure = 60 C (from steam table)
(high molecular wt. solute, BPE is neglected)
Latent heat vapor generated in effect II at 160 mm Hg (0.2133bar) = 563 kcal.kg
-1
(
2 s
)
Heating steam is at 3.5 kg.cm
-2
gauge; temp (
s
T )=148 C; Latent heat (
s
)= 506
kcal.kg
-1
Feed rate = 8000 kg.h
-1
, Solute content = 5%
Final concentration = 40%
Solid in = 8000 0. 05 = 400 kg.h
-1
, water in = 8000-400= 7600 kg.h
-1
Product out (40% solid) = 400/0.40 kg.h
-1
= 1000 kg.h
-1
Water out with the product = 1000 (1-0.40) kg.h
-1
= 600 kg.h
-1
Total evaporation rate in two effects (
1 2 s s
m m + ) = 7600- 600 = 7000 kg.h
-1
(3.A)
Allow equal areas to two effects, i.e.,
1 1 2 2 D D
U T U T A = A
1 2 tot
T T T A = A +A = 148-60 = 88 C
NPTEL Chemical Engineering Chemical Engineering Design - II
Joint initiative of IITs and IISc Funded by MHRD Page 25 of 31
U
D1
= 550, U
D2
= 370 kcal.h
-1
m
-2
C
-1
1
T A = 35.4 C &
2
T A = 52.6 C (3.11)
Temperature of the vapor leaving effect I (T
b1
) =148 35.4 = 112.6 C
Latent heat vapor generated in effect I at 112.6 C = 531 kcal.kg
-1
(
1 s
)
Energy balance for effect I:
1 1 1 1
( )
f f s s f s s s
m H m m m H m H + = +
(3.2)
Enthalpy values: reference temperature = 112.6 C (temp of solution leaving effect I)
f
H = (40 - 112.6)(1 kcal/ kg C) = -72.6 kcal/ kg
1
H = 0 kcal/ kg (w.r.t. the reference temperature of 112.6 C)
s
= 506 kcal/ kg;
1 s
= 531 kcal/ kg
(8000)(-72.6) + m
s
(506) = (8000-m
s1
)(0) + m
s1
(531)
m
s
= 1.05m
s1
+ 1148 (3.B)
Energy balance for effect II:
1 1 1 1 1 2 2 2 2
( ) ( )
f s s s f s s s s
m m H m m m m H m H + = + (3.
3)
Enthalpy values: reference temperature = 60 C (temp of solution leaving effect II)
1
H = (112.6 - 60)(0.94 kcal/ kg C)= 49.4 kcal/ kg
(Specific heat of the solution leaving the 1
st
effect is taken as the mean value of the sp.
heat of feed and the concentrated liquor, i.e.,
1 0.87
2
+
=0.94)
2
H = 0 kcal/ kg (w.r.t. the reference temperature of 60 C)
1 s
= 531 kcal/ kg;
2 s
= 563 kcal/ kg
(8000 - m
s1
)(49.4) + m
s1
(531) = (8000-m
s1
- m
s2
)(0) + m
s2
(563)
m
s2
= 0.855m
s1
+ 702.6
(3.C)
Solving Eqs. 3.A, 3.B and 3.C for m
s,
m
s1
and m
s2
:
m
s
= 4713;
m
s1
= 3395 and m
s2
= 3605 kg/ h
Areas:
2
1
1 1
(4713)(506)
122.5 m
(550)(35.4)
s s
D
m
A
U T
= = =
A
;
2 1 1
2
2 2
(3395)(531)
92.6 m
(370)(52.6)
s s
D
m
A
U T
= = =
A
The areas in the two effects are not equal. Revised calculation is required.
NPTEL Chemical Engineering Chemical Engineering Design - II
Joint initiative of IITs and IISc Funded by MHRD Page 26 of 31
The revised temperature difference in 1
st
effect,
o o
1 2
40 C and 48 C T T A = A = taken for
the calculation (you may also continue the calculation with this revised
value
o
1,
38.84 C
new
T A = (3.13)
Temperature of the saturated vapor from the first effect = 148 -60= 108 C
Corresponding evaporator drum pressure = 1.2116 bar
Latent heat of the vapor leaving the 1
st
effect at 108 C= 530 kcal/ kg (
1 s
)
Corresponding pressure in the vapor drum = 1.3317 bar
Revised calculation: Energy balance for effect I (reference temperature 108 C)
(8000)(-72.6) + m
s
(506) = (8000-m
s1
)(0) + m
s1
(530)
(3.2)
m
s
= 1.047m
s1
+ 1075
(3.D)
Revised calculation: Energy balance for effect II (reference temperature 60 C)
(8000 - m
s1
)(0.94)(108-60) + m
s1
(530) = (8000-m
s1
- m
s2
)(0) + m
s2
(563)
(3.3)
m
s2
= 0.8612m
s1
+641
(3.E)
Solving 3.A, 3.D and 3.E, m
s
= 4652;
m
s1
= 3417 and m
s2
= 3583 kg/ h
Areas:
2
1
1 1
(4652)(506)
107 m
(550)(40)
s s
D
m
A
U T
= = =
A
;
2 1 1
2
2 2
(3417)(530)
102 m
(370)(48)
s s
D
m
A
U T
= = =
A
These areas are fairly close, use A
1
= A
2
= 104 m
2
plus overdesign. About 10% excess
area will provide a reasonable overdesign.
Tube details:
Let us select 1 inch nominal diameter, 80 schedule, brass tubes of 12 ft in length
Outer tube diameter (d
o
) = 42.16 mm
Inner tube diameter (d
i
) = 32.46 mm
Tube length (L) = 12 ft = 3.6576 m
Surface area of each tube (a) = d
o
L = 42.16410
-3
3.6576 = 0.4845 m
2
Number of tubes required providing 10% overdesign (N
t
) = A /a = (115/0.4845) ~238
Tube pitch (triangular), P
T
= 1.25 d
o
= 1.25 42.164 = 52.71 53 ~ mm
Total area occupied by tubes = N
t
(1/2) P
T
P
T
sin (where = 60)
= 238 0.5 (5310
-3
)
2
0.866
NPTEL Chemical Engineering Chemical Engineering Design - II
Joint initiative of IITs and IISc Funded by MHRD Page 27 of 31
= 0.2894 m
2
This area is generally divided by a factor which varies from 0.8 to 1 to find out the
actual area. This allows for position adjustment of peripheral tubes as those cant be
too close to tube sheet edge.
Actual area required = 0.2894/ 0.9 (0.9 is selected)
= 0.3216 m
2
The central downcomer area is generally taken as 40 to 70% of the total cross
sectional area of tubes. Consider 50% of the total tube cross sectional area.
Therefore, downcomer area = 0.5 [N
t
(/4) d
o
2
]
= 0.5 [238 (/4) (0.04216)
2
]
= 0.1661 m
2
Downcomer diameter = [(4 0.1661)/]
= 0.460 m
Total area of tube sheet in evaporator = downcomer area + area occupied by tubes
= 0.1661+ 0.3216 m
2
= 0.4877 m
2
Tube sheet diameter = [(4 0.4877)/ ]
= 0.788 m
NPTEL Chemical Engineering Chemical Engineering Design - II
Joint initiative of IITs and IISc Funded by MHRD Page 28 of 31
Lecture 5: Mechanical Design
Part 2. Mechanical design
A few basic parts of mechanical design of evaporator are shown below. The
systematic approach of mechanical design of shell and tube heat exchanger is
discussed in module #1.
Process design parameters and materials of construction:
Consider a standard vertical short tube evaporator (calandria type) for this service
Evaporator drum operated at 1.3317 bar pressure
Amount of water to be evaporated = 7000 kg/hr
Heating surface required A =115 m
2
Steam is available to first effect at 3.433 bar pressure
Density of 5% feed liquid (
l
)= 1050 kg/m
3
(assumed)
Density of water vapor (
v
)=
PM
RT
at saturated vapor temperature in the first effect
(108C).
= (1.331718) / (8.314 10
-5
381)
= 0.757 kg/m
3
Design pressure (P) = 5% more than the maximum working pressure
= 1.053.433
= 3.605 bar
= 3.6749 kgf/cm
2
Volumetric flow rate of water vapor (V ): 7000/(0.7573500)= 2.57 m
3
/s
Evaporator shell: low carbon steel (IS-2062)
Tube material: brass
Permissible stress for low carbon steel = 980 kg/cm
2
Modulus of elasticity for low carbon steel = 1910
5
kg/cm
2
Modulus of elasticity for brass = 9.510
5
kg/cm
2
Conical head at bottom: cone angle =120
Conical head at top: cone angle= 120
NPTEL Chemical Engineering Chemical Engineering Design - II
Joint initiative of IITs and IISc Funded by MHRD Page 29 of 31
1. Check for tube thickness
The tube thickness is given by
2
i
t
Pd
t
fJ P
=
The permissible stress of brass ( f ) = 538 kg/cm
2
; Welding or joint efficiency, J = 1 is
used for seamless tube
Therefore, t
t
= (3.6749 32.46) / [(2 538 1) 3.6749]
= 0.111 mm
The specified thickness is 4.85 mm. Therefore, the selected tube is suitable for this
service.
2. Calandria sheet thickness calculation
Thickness is given by:
2
o
s
PD
t
fJ P
=
+
= (3.6749 785) /[(29800.85) +3.6749]
= 1.73 mm
Normally, the corrosion allowance of 3 mm is used for carbon steel.
It may be taken as t
s
= 10 mm
3. Tube sheet thickness
Please refer module #2 (section 2.3.5) for the calculation of tube sheet thickness
according to the TEMA code and the minimum tube sheet thickness as per IS:4503
specification also needs to be checked.
4. Evaporator drum diameter determination
The following equation helps to estimate the drum diameter. The diameter of the
drum may be same as that of the calandria. However, it is necessary to check the size
from the point of satisfactory entrainment separation ([5] page 253].
( )
/ 0.0172 ( ) /
d l v v
V
R
A
| |
=
|
\ .
where, V=volumetric flow rate of vapor [m
3
/s]
A=cross sectional area of drum
For drums having wire mesh as entrainment separator device, R
d
may be taken as 1.3.
( )
0.0172 ( ) /
d l v v
V
A
R
| |
|
=
|
\ .
= (2.57) / [1.30.0172 {(1050 0.757)/0.757}]
= 3.087 m
2
NPTEL Chemical Engineering Chemical Engineering Design - II
Joint initiative of IITs and IISc Funded by MHRD Page 30 of 31
The drum diameter= {(4 3.087)/3.14}
= 1.98 m
Drum height is usually taken as 2 to 5 times of tube sheet diameter
Thus, the drum height = 40.788 = 3.152 m
The tube length is in between the border line of short tube (4 to 10) and long tube
( >12). Therefore higher head space is required to arrest entrainment by mis-tallying
deflector plates.
5. Drum thickness
Drum is operating at 1.3317 bar pressure
Design pressure = 1.398 bar
= 1.4253 kgf/cm
2
Drum thickness:
2
o
d
PD
t
fJ P
=
+
= (1.4253 785) /[(29800.85) +1.4253]
= 0.671 mm
Therefore the same thickness of 10 mm for both the drum and calandria sheet may be
used including the tolerance for corrosion.
Note: The drum thickness is to be calculated based on the external pressure if the
drum is under vacuum.
Practice problem:
A triple-effect forward feed evaporator of the long-tube vertical type is to be used to
concentrate 4,000 kg/h of a 9.5% solution of caustic soda available at 40C to 50%
solution. Saturated steam at 3.5 kg.cm
-2
is available. A vacuum of 700 mm Hg is
maintained in the last effect. The overall heat transfer coefficients are 5800, 3300 and
2400 Wm
-2
C
-1
, corrected for BPEs may be used for the 1
st
, 2
nd
and 3
rd
effects
respectively. Calculate the heat transfer area required (assume equal areas in all three
effects), steam economy and rate of steam consumption.
NPTEL Chemical Engineering Chemical Engineering Design - II
Joint initiative of IITs and IISc Funded by MHRD Page 31 of 31
Reference books
[1]. McCabe W.L., Smith J.C. and Harriott P. Unit Operations of Chemical
Engineering, McGraw-Hill, 5
th
Ed., (1993), pp. 463-494.
[2]. Dutta B.K. Heat Transfer-Principles and Applications, PHI Pvt. Ltd., New
Delhi, 1
st
Ed. (2006), pp. 361-420.
[3]. Parker N.H. How to specify evaporators, Chem. Eng., May 27, 1963, pp.
135-140.
[4]. Minton P.E. Hand Book of Evaporation Technology, Noyes Publications,
New Jersey, 1986.
[5]. Joshi M.V. and Mahajani V.V. Process Equipment Design, McMillan
Publishers India Ltd., 3
rd
ed., 1996, New Delhi.