Evaporator Design Calculation
Evaporator Design Calculation
Formulas and design parameters in Robert Evaporator Design with online calculation sheet
1.1 Types of Evaporators:
1.2 Steps in calculation of evaporator design:
1.2.1 Number of tubes (N ):
1.2.2 Tube plate & Downtake dia :
1.2.2.1 Area occupied for tubes in tube plate :
1.2.2.2 Dia of the down take :
1.2.2.3 The final required tube plate diameter.
1.2.3 Dia required for vapour inlet & Dia of the Calendria in radial steam/ vapour entry
1.2.3.1 Dia of the Vapour Inlet :
1.2.3.2 Evaporator Bodies Vapour Velocities Recommended by E.Hugot
1.2.3.3 Calendria dia at the entry of the steam/vapour jocket
1.2.4 4. Vapour outlet pipe dia :
1.2.5 5. Dia of the condensate line :
1.2.6 6. Noxious gases connections :
1.2.7 7. Vapour space height :
1.2.8 8. Velocity in vapour space of body (Cross checking of the system ):
1.2.9 9. Calendria shell thickness :
1.2.10 10. Vapour shell thickness :
1.2.11 11. Tube plate thickness :
1.2.11.1 Generally shell thickness will be taken as follows as (in mm).
1.2.12 12. Vapour doom dia :
1.2.13 13. Top cone Height :
1.2.14 14. Center Umbrella area :
1.3 Online Calculation Sheet for Robert Evaporator Body Design
ulation sheet
ur entry
Types of Evaporators:
1. Robert type natural circulation evaporator.
2. Long tube vertical raising film evaporator (Kestner)
3. Long tube vertical falling film evaporator (FFE)
4. Plate type evaporator.
5. Thin film evaporators, high viscous liquid
6. Inclined tube evaporators
7. Horizontal tube evaporators.
8. Flash evaporators.
9. Compact evaporators.
Steps in calculation of evaporator design:
1. Number of tubes (N ):
Mean dia of the tube ( Dm) in mtr = Tube OD- Tube Thickness ( In some designers also take ID of the tube in the place of mean dia.
Effective Length of the tube ( L ) = Tube length – 2(Tube plate thickness)-2(Tube expansion allowance)
Number of tubes = Heating surface / π x Dm x L.
The thickness of the tube generally take for juice heaters and for evaporator 18G and for pans take 16G. ( 18g = 1.22mm , 16g = 1.625
Take extra dia in percentage on area occupied for tubes in tube plate for stay roads arrangement, free withdrawal of condensate and no
purpose. while providing the multiple down design than this percentage may go higher side. Its value lies in the range 10 to 20% on ar
Tube plate area required for tubes only ( AT ) = (0.866 x P2 x N /β) x %extra
Tube plate dia required for tubes only = SQRT ( AT x 4/π )
According to peter rein down take dia consider less than 25% of the tube plate dia.
According to E. Hugot The diameter of the centre well varies from ¼ to ⅛ of the interior diameter of the vessel.
Certain manufacturers replace the centre well by a lateral well or by a series of down takes of small diameter distributed over the calen
downtakes).
From the above generally downtake dia take 20% on tube plate.
Dia of the single downtake = Tube plate dia for tubes x % of downtake on tube plate.
Dia of the central downtake in multiple down takes design = SQRT [(Area of the single downtake – Total area of peripheral down take
(Note : This dia to be maintained at the vapour entry side later it may reduced in vapour travel direction.).
4. Vapour outlet pipe dia :
Vapour volume outlet vapour in M3/sec = Heating surface X Evap. Rate x Specific volume of outlet
vapour/3600.
Vapour outlet pipe dia in mtrs = SQRT [vapour volume /(0.785 x velocity of vapour)]
Volume of the condensate in M3/sec = [Heating surface X Evap. Rate ] / [ Density of water x 3600].
Dia of the each condensate line = SQRT (Volume of the condensate each./(0.785 x velocity of condensate)).
6. Noxious gases connections :
Generally 10 m2 heating surface area required 1cm2 area for removal of noncondensable gases.
Cross section area of non condensable gases in cm2 = Heating surface in m2 /10
Dia of the each non condensable gases line = SQRT( Total area of non condensable gases /0.785*no. of points)
Generally in evaporator design vapour space dia may be take same as to calendria dia meter. But we have to
check, how much vapour velocity maintained in body. It is helpful to avoid the entertainment of the system. The
velocity of vapour leaving the liquid surface would then be approximately 10 cm/sec.
In present designs, The vapour velocities in the vessels to be maintained below 3.6m/sec and for lost body it can
go upto 4.6 m/sec . Hence the necessity of providing entrainment separators or save all at the vapour outlets
from the evaporator vessels.
In any case, it is considered that entrainment rapidly becomes excessive from the moment when the vapour
velocity in the vessel reaches more than 6m/sec.
%
Value
0
0
0
0
20
0
0
0
0
0
0
0
0
200
8
0
0
0
ntry
2
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
2
0
0
0
6
0
0
2
0
0
0
0
0
3
0
1400
0.7
1.5
0
0
2
1400
0.7
1.5
0
0
0
1.5
1400
2.7 2
2100 00 0
2.7 2
2100 00 0
1900 00 0
0
0
0
0
0
0
0
0
0
2
0
0
0
30
0
0
0
Inlet vapour temperature 1 03
outlet vapour temperature 93
Velocity of inlet vapour 30
velocity of outlet vapour 35
Velocity of condensate 0 .6
Evaporation rate of the body 25
Inlet vapour specific volume 0
outlet vapour specific volume 0
Formula
Generally for Robert type bodies taken for lost effect 2.5 times on calendria tube height and for remaining bodies take 2 tmes o
The velocity of vapour leaving the liquid surface would then be approximately 10 cm/sec
m/sec
m/sec
m/sec
Kgs/m2 /hr
M3/kg
m3/kg
ing bodies take 2 tmes on Calendria tube height
EVAPORATOR
Pressure 13 Bar
mass flow rate 2.22 Kg/Sec
Tsat @ 13 bar 191.6 deg C
Q = m x ∆h
Q= 4488 kW
Also,
Tg, out 321 Deg.C
Q = UA(LMTD)
For LMTD,
Tg, in Tg,out Tw,in
373.4 321 175
Tw,out Tg,in - Tw, out Tg,out - Tw, in
191.6 181.8 146
LMTD LN
0.219
Re^m 441.1
Pr^1/3 0.893
ho (W/m2K) 86.32
vf (m3/kg) 1.14E-03
(rho)f 714.2
vg (m3/kg) 0.15
(rho)g 6.7
Now,
Nuselt Nu 1530.1
hi 19321.7
do/di 1.181
(1/ hi)(do/di) 4.3808515056E-05
Rfi (do/di) 0.000236279069767
do/2k.ln(do/di) 8.82100930798E-05
1/ho 0.011584800158949
Rfo 0.001
1/U 0.013
U W/m2K 77.2
A=Q/U*LMTD
Area A m2 356.1
FRONTAL AREA
Aff = mg/Aff
Aff= (ST-do)Nw.L
Putting we get
mg/G = (ST-do)Nw.L
So,
Nw.L = mg/(ST-do).G
mg kg/s 49.86
ST m 0.15
do m 0.0508
G kg/m2s 9
Nw.L 55.8
Area A = pi.do.L.Nw.Nh
Nh = A/pi.do.Nw.L
Nh A/pi.do Nw.L Nw
9 963.29 107.03 10
10 963.29 96.33 11
10 963.29 96.33 12
15 963.29 64.22 20
69.4 963.29 13.87 30
50 963.29 19.27 50
20 963.29 48.16 25
25 963.29 38.53 20
20 963.29 48.16 20
25 963.29 38.53 25
30 963.29 32.11 30
35 963.29 27.52 35
50 963.29 19.27 40
SPECIFICATIONS:
Nw 20
Nh 20
V m/s 5.78 m/s
L m 9.62
ST m 0.15
do m 0.0508
di m 0.04445
Length L
10.70
8.76
8.03
3.21
0.46
0.39
1.93
1.93
2.41
1.54
1.07
0.79
0.48
Q = m x ∆h
Q= 4488 kW
Also,
Tg, out 321 Deg.C
Q = UA(LMTD)
For LMTD,
Tg, in Tg,out Tw,in Tw,out Tg,in - Tw, out Tg,out - Tw, in
0.219300560035042
Re^m 441.1263
Pr^1/3 0.893367
ho (W/m2K) 86.32
WATER SIDE CALCULATION
Inside fouling factor Ri 0.0002
Outside fousing factor Ro 0.001
K metal W/m.K 63.9
Circulation ratio 4
vf (m3/kg) 1.14E-03
(rho)f 714.2
vg (m3/kg) 0.15
(rho)g 6.67
Now,
Nuselt Nu 1693.1
hi 18838.9
do/di 1.041
(1/ hi)(do/di) 5.525708E-05
Rfi (do/di) 0.000208
do/2k.ln(do/di) 1.596584E-05
1/ho 0.0116
Rfo 0.001
1/U 0.01286422
U W/m2K 77.73
A=Q/U*LMTD
Area A m2 353.665759
FRONTAL AREA
Aff = mg/Aff
Aff= (ST-do)Nw.L
Putting we get
mg/G = (ST-do)Nw.L
So,
Nw.L = mg/(ST-do).G
mg kg/s 83.5
ST m 0.15
do m 0.0508
G kg/m2s 9
Nw.L 93.52599
Area A = pi.do.L.Nw.Nh
Nh = A/pi.do.Nw.L
SPECIFICATIONS:
Nw 20 LMTD deg C 75.153
Nh 20 A m2 614.4231
V m/s 5.78 m/s U W/m2.K 255.12
L m 9.62
ST m 0.15
do m 0.0508
di m 0.04445
Tg,out - Tw, in
EVAPORATOR
Pressure 13 Bar
mass flow rate 2.22 Kg/Sec
Tsat @ 13 bar 191.6 deg C
Q = m x ∆h
Q= 4488 kW
Also,
Tg, out 321 Deg.C
Q = UA(LMTD)
For LMTD,
Tg, in Tg,out Tw,in Tw,out Tg,in - Tw, out Tg,out - Tw, in
0.219300560035042
ho (W/m2K) 114.8
vg (m3/kg) 0.15
(rho)g 6.666667
hi 20000
do/di 1.0471142521
(1/ hi)(do/di) 5.235571E-05
Rfi (do/di) 0.0002094229
do/2k.ln(do/di) 1.601245E-05
1/ho 0.0087108014
Rfo 0.001
1/U 0.00998859
U W/m2K 100.114206
A=Q/U*LMTD
Area A m2 274.608423
FRONTAL AREA
Aff = mg/Aff
Aff= (ST-do)Nw.L
Putting we get
mg/G = (ST-do)Nw.L
So,
Nw.L = mg/(ST-do).G
mg kg/s 49.86
ST m 0.15
do m 0.0508
G kg/m2s 9
Nw.L 55.84677
Area A = pi.do.L.Nw.Nh
Nh = A/pi.do.Nw.L
Q = m x ∆h
Q= 4488 kW
Also,
Tg, out 321 Deg.C
Q = UA(LMTD)
For LMTD,
Tg, in Tg,out Tw,in Tw,out Tg,in - Tw, out
0.248
Re^m 441.13
Pr^1/3 0.89
ho (W/m2K) 86.3200043401
WATER SIDE CALCULATION
Inside fouling factor Ri 0.0002
Outside fousing factor Ro 0.001
K metal W/m.K 42
Circulation ratio 4
vf (m3/kg) 1.14E-03
(rho)f 714.2
vg (m3/kg) 0.15
(rho)g 6.67
Now,
Nuselt Nu 1530.1
hi 19321.7
do/di 1.181
(1/ hi)(do/di) 6.114332412E-05
Rfi (do/di) 0.0002362790698
do/2k.ln(do/di) 0.00010
1/ho 0.01158
Rfo 0.001
1/U 0.013
U W/m2K 77.0
A=Q/U*LMTD A=Q/U*LMTD
Area A-m2 362.0
FRONTAL AREA
Aff = mg/Aff
Aff= (ST-do)Nw.L
Putting we get
mg/G = (ST-do)Nw.L
So,
Nw.L = mg/(ST-do).G
mg kg/s 49.86
ST m 0.15
do m 0.0508
G kg/m2s 9
Nw.L 55.85
Area A = pi.do.L.Nw.Nh
Nh = A/pi.do.Nw.L
SPECIFICATIONS:
Nw 20 LMTD deg C
Nh 20 A m2
V m/s 5.78 m/s U W/m2.K
L m 9.62
ST m 0.15
do m 0.0508
di m 0.04445
Tg,out - Tw, in
141.8
75.153
614.4231
255.12
EVAPORATOR
Pressure 13 Bar
mass flow rate 2.22 Kg/Sec
Tsat @ 13 bar 191.6 deg C
Q = m x ∆h
Q= 4488 kW
Also,
Tg, out 321 Deg.C
Q = UA(LMTD)
For LMTD,
Tg, in Tg,out Tw,in Tw,out Tg,in - Tw, out Tg,out - Tw, in
0.248489567645351
Re^m 441.1263
Pr^1/3 0.893367
ho (W/m2K) 86.32
WATER SIDE CALCULATION
Inside fouling factor Ri 0.0002
Outside fousing factor Ro 0.001
K metal W/m.K 63.9
Circulation ratio 4
vf (m3/kg) 1.14E-03
(rho)f 714.2
vg (m3/kg) 0.15
(rho)g 6.666667
Now,
Nuselt Nu 1693.07396
hi 18838.9172
do/di 1.0409836066
(1/ hi)(do/di) 5.525708E-05
Rfi (do/di) 0.0002081967
do/2k.ln(do/di) 1.596584E-05
1/ho 0.0115848002
Rfo 0.001
1/U 0.01286422
U W/m2K 77.7
A=Q/U*LMTD
Area A m2 358.7
FRONTAL AREA
Aff = mg/Aff
Aff= (ST-do)Nw.L
Putting we get
mg/G = (ST-do)Nw.L
So,
Nw.L = mg/(ST-do).G
mg kg/s 83.5
ST m 0.15
do m 0.0508
G kg/m2s 9
Nw.L 93.52599
Area A = pi.do.L.Nw.Nh
Nh = A/pi.do.Nw.L
SPECIFICATIONS:
Nw 20 LMTD deg C 75.153
Nh 20 A m2 614.4231
V m/s 5.78 m/s U W/m2.K 255.12
L m 9.62
ST m 0.15
do m 0.0508
di m 0.04445
Tg,out - Tw, in
CLIMBING / RISING FILM EVAPORATOR
Cyclone seperator
Steam
Condensate
To condenser
Cyclone seperator
Feed from storage
FALLING FILM EVAPORATOR
To Condenser
Concentrate withdrawl
PORATOR
Feed to FFE
Outside dry bulb temp. (DBT) 115 F
Outside wet bulb temp. (WBT) 68 F
Inside dry bulb temp. (T) 84 F
Total sensible heat gain 36000 BTU
Media Saturation Efficiency (Eff.) 0.85