DESIGN OF A TRIPPLE EFFECT FORWARD FEED EVAPORATOR FOR THE

CONCENTRATION OF GLYCEROL SOLUTION:

Design: Glycerol solution is to be concentrated from 7Wt% to a final concentration of 80Wt% The
feed is available at 50993 Kg/Hr and 27C. Steam is supplied at 103.66Kpa (Abs.) to the first effect
and a vacuum of 74.2(Abs.) mm of Hg is maintained in the last effect. Heat losses by radiation and
by entertainment are neglected. Condensate are assumed to enter each effect at its saturation
temperature.
Glycerol in feed = 50993 0.07
= 3569.51 Kg/Hr
Feed rate WF = 50993/3600
= 14.16 Kg/Sec
Bottom outlet from third effect = 3569.51/(0.83600)
= 1.24 Kg/Sec
Total evaporation =14.16 1.24
=12.9 Kg/Sec
Assuming equal evaporation in all the three effects.
W1 Evaporation rate in first effect
W2- Evaporation rate in second effect
W3 Evaporation rate in the third effect
W1 = W2 = W3 = (12.9/3) = 4.3 Kg/sec
Outlet from first effect = WF - W1
= 14.16 4.3
= 9.86 Kg/sec
Outlet from second effect = WF W1- W2
= 14.16 4.3 4.3
= 5.56 Kg/sec
Outlet from third effect = WF W1 W2 W3
= 14.16 4.3 4.3 4.3
= 1.26 Kg/sec
Concentration of glycerol in outlet from first effect = (14.16 0.07 )/ 9.86
= 10 wt%

Concentration of glycerol in outlet from second effect = (14.16.07)/5.56

= 0.178
= 17.8 wt%
Concentration of glycerol in outlet from third effect = (14.16.07)/1.26
= 0.79
= 79 wt%
Saturation temperature of inlet steam TS = 100C
Boiling point of water in the third effect = 12C
Temperature difference =100 12
= 88C
BOILING POINT RISE FOR GLYCEROL SOULTION
1) BPR for the third effect at 79 Wt% glycerol = 10.4C
2) BPR for the second effect at 17.8 Wt% glycerol = 1.5C
3) BPR for the first effect at 10 Wt% glycerol = 0.5C
Sum of the BPR of the glycerol solution = 10.4 + 1.5 + 0.5
= 12.4C
Effective temperature difference = 88 12.4
= 75.6C
Now we assume that the amount of heat transfer from steam to the solution is the same in all the
three effects i.e. Q (heat transfer in the first effect) = Q (heat transfer in the second effect) = Q
(heat transfer in the third effect)
OVERALL HEAT TRANSFER COEFFICIENT
Since it is a vertical effect evaporator, the overall heat transfer coefficient can safely be assumed as
follows:
In the first effect U1 = 2325 W/mC
In the second effect U2 = 1275 W/mC
In the third effect U3 = 1031 W/mC
Now as we have assumed heat transfer rates to be equal , we have
Q1 = Q2 = Q3
(U1 A1 t1 ) = (U2 A2 t2 ) = (U3 A3 t3 )

We design the triple effect evaporator such that the heating area in all the three effect
is the same
A1 = A2 = A3
U1 t1 = U2 t 2 = U3 t3
t2 / t1 = U1 / U2
= 2352 / 1275
= 1.8
t3 / t 2 = U2 / U3
=1275 / 1031
= 0.6
t1 + t2 + t3 = 75.6C
0.6 t2 + t2 + 0.6 t2 = 75.6C
t2 = 32.3C
t1 = 21.54C
t3 = 21.64C
ACTUAL BOILING POINTS IN EACH EFFECT
First effect:
T1 = TS - t1
= 100 - 21.54
= 78.54C
Second effect:
T2 = T1 - (BPR)1 - t2
= 78.5 - 0.5 - 32.3
= 45.7C
Third effect:
T3 = T2 - (BPR)2 - t3
= 45.7 - 1.5 - 21.64
= 22.56C

Effect 1 ( C )

Effect 2 ( C )

Effect 3 ( C )

TS = 100

T1 = 78

T2 = 44.2

T1 = 78.5

T2 = 45.7

Condenser ( C )

T3 = 12.16

T3 =22.56

HEAT BALANCE
FIRST EFFECT:
WS S + WFHF = W1H1 + ( WF - W1 ) h1
latent heat of steam S = 2257.86 KJ/ kg
HF Enthalpy of feed at inlet temperature ( 27C) = Cpf ( Tf - 0)
= ( 0.576 4.18) 27
= 65 KJ / Kg
H1- Enthalpy of vapor leaving the first effect = H2S + (Cp)steam (BPR1)superhea
= 2640 + (1.884 0.5)
= 2487 KJ/ Kg
H2s- Enthalpy of steam at 78C = 2640 KJ / Kg
(Cp)steam at 78C = 1.884 KJ/ Kg
h1 - enthalpy of outlet from first effect at 78.5C = Cp1 ( t1 0 )
= 0.65 4.18 78.5
= 213.28 KJ/Kg
WS (2257.86) + (14.16 65) = ( W1 2487 ) + ( 14.16 W1 ) 213.28
WS (2257.86) = 2273.7 W1 + 2099.6
WS = W1 + 0.93

------------------------(1)

SECOND EFFECT:
W1 1 + ( WF W1 ) h1 = W2 H2 + ( WF W1 W2 ) h2
Latent heat of steam at 78C 1 = H1 - h2S
= 2487 - 325
= 2162 KJ /Kg
H3S Enthalpy of steam vapor at 44.2C = 2580 KJ /Kg
H2 - Enthalpy of vapor leaving the second effect = H3S + ( Cp)steam (BPR2)supreheat
= 2580 + ( 1.884 1.5 )
= 2583 KJ/Kg
h2-Enthalpy of outlet from the second effect at 45.7C = Cp2 ( t2 - 0 )
= 0.6 4.18 45.7
= 114.62 KJ /Kg
W1 2162 + ( 14.16 W1 ) 213.28 = (W2 2583 ) + ( 14.16 W1 W2 )114.62
2063.34 W1 = 1397 + 2468.38 W2
W1 + 0.667 = 1.196 W2 -------------------(2)
THIRD EFFECT:
W2 2+ (WF W1 W2 ) h2 = W3 H3 + (WF W1 - W2 - W3 ) h3
Latent heat of steam at 44.2C 2 = H2 - h3S
= 2583 - 190
= 2393 KJ / Kg
H4S -Enthalpy of steam vapor at 12 .16C = 2523 KJ/Kg
H3 -Enthalpy of vapor leaving the third effect = H4S + (Cp)steam (BPR3)ssuperheat
= 2523 + (1.884 10.4 )
= 2543 KJ/ Kg
h3 Enthalpy of outlet from third effect at 22.56C = Cp3 ( t3 - 0)

= 0.57 4.18 22.56

= 53.75 KJ /Kg
W2 2393 + (14.16 W1 W2 ) 114.62 = W3 2543 + ( 14.16 W1 W2 W3 ) 53.75
2332.13 W2 - 60.87 = 2489.25 W3 + 861.92
W2 + 0.370 = 1.067 W3 + 0.026 W1 --------------------(3)
W1 + W2 + W3 = 12.9 Kg/ Sec -------------------------(4)
Solving equations' (1),(2),(3) and (4), we get:
WS = 5.153 Kg/sec
W1 = 4.223 Kg/sec
W2 = 4.089 Kg/sec
W4 = 4.588 Kg/sec
Now, Q1 = WSS
=11634.75 KJ/sec
But Q1=U1A1t1
Therefore A1= (11634.75103)/(21.542325)
= 232.32 m2
Q2 = W11
= 9130.13 KJ/sec
But Q2 = U2A2t2
Therefore A2 = (9130.13103)/(32.31275)
= 228 m2
Q3 = W22
= 9784.98 KJ/sec
But Q3 = U3A3t3
Therefore A3 = (9784.98103)/(21.641031)
= 235 m2
Thus the obtained areas are within the acceptable range of 5% difference
Therefore the average area per effect of the evaporator is 232 m2.

Tube details:
Most generally used diameters today ranges from 1.25 to 2.00 in. outer diameter and most generally
used lengths of tubes ranges from 4 to 15 ft.
Let us choose 5/4-in. nominal diameter, 80 schedule, brass tubes of 10-ft length.
Therefore Outer diameter do = 42.164 mm
Inner diameter di = 32.46 mm
Length L = 10 ft
= 3.048 m
Tube pitch ()PT = 1.25 do
= 1.25 42.164
= 52.705 mm
Surface area of each tube a = doL
= 52.70510-3 3.048
= 0.4037 m2
Number of tubes required Nt = A /a
= 619
Area occupied by tubes = Nt (1/2) PT PT sin
= 619 0.5 (52.705 10-3)2 0.866
= 0.7445 m2
Where = 60o
But actual area is more than this. Hence this area is to be divided by factor which varies from 0.8 to
1.0.
Let us choose this factor as 0.9.
Therefore actual area required = 0.7445/ 0.9
= 0.827 m2
The central downcomer area is taken as 40 to 70% of the total cross sectional area of tubes. Let us
take it as 50%.
Therefore Downcomer area = 0.5 [Nt (/4) do2]
= 0.5 [619 (/4) (0.04216)2]
= 0.432 m2
Downcomer diameter = (4 0.432) /

= 0.742 m
Total area of tube sheet in evaporator = downcomer area + area occupied by tubes
= 0.432+ 0.827
= 1.259 m2
Thus tube sheet diameter = (4 7.1025)/
= 1.27 m

MECHANICAL DESIGN OF EVAPORATOR:

Take standard vertical short tube evaporator (calendria type)
Data:
Evaporator drum operating under vacuum at 0.4163 bar
Amount of water to be evaporated = 15480 kg/hr
Heating surface required A = 232 m2
Steam is available to first effect at pressure of 1.03 bar.
Density of liquid (10 wt% glycerol) = 1019 kg/m3
Density of water vapour = PM/RT
= (0.4163 105 18) / (8314 351)
= 0.258 kg/m3
Design pressure = 5% extra of maximum working pressure
= 1.05 1.03
= 1.082 bar
= 1.103 kgf/cm2
Material:
Evaporator low carbon steel
Tubes brass
Permissible stress for low carbon steel = 980 kg/cm2
Modulus of elasticity for low carbon steel = 19 105 kg/cm2
Modulus of elasticity for brass = 9.5 105 kg/cm2
Conical head at bottom
Cone angle 120o
Conical head at top
Cone angle 120o
1) Calendria sheet thickness:
The thickness is given by
ts = (PDo) /(2fJ + P)
= (1.1031270) /[(2 980 0.85) +1.103]
= 0.84 mm
The actual thickness must be much more so as to allow for corrosion and give rigidity to the shell.

Therefore it may be taken as ts = 10 mm

2) Tube sheet thickness:
To find tube sheet thickness
K = [Ests(Do ts)] / [EtNttt (do tt)]
Es = elastic modulus of shell
Et = elastic modulus of tube
Do - outside diameter of shell = 3 m
do - outside diameter of tube = 60.325 mm
ts - shell thickness = 12 mm
tt - tube wall thickness = 5.5 mm
Nt - number of tubes in shell = 2108
Therefore K = [19 105 12 (1270 10)] / [9.5 105 619 4.85 (42.164 4.85)]
= 0.027
F = [(2 + K) / (2 + 3 K)]
= [(2 + 0.027) / (2 + 3 0.027)]
= 0.99
The effective tube sheet thickness is given by
tts = FDo [(0.25 P) / f]
= 0.991270[(0.25 1.103) / 980]
= 21 mm
with corrosion allowance the thickness may be taken as 25 mm.
4) Check for tube thickness:
The tube thickness is given by
tt = Pdi / (2fJ P)
The permissible stress for brass = 538 kg/cm2 and J = 1
Therefore tt = (1.103 32.46) / [(2 538 1) 1.103]
= 0.033 mm
But provided thickness is 4.85 mm.Therefore chosen tubes have enough strength to withstand within
operating conditions.

5) Evaporator drum diameter:

The following equations help to determine the drum diameter. The diameter of the drum may be
same as that of the calendria. However it is necessary to check the size from the point of satisfactory
entrainment separation.
Rd = (V/A) / [0.0172 {(l - v) / v}]
Where V volumetric flow rate of vapour in m3/sec
A cross sectional area of drum
For drums having wire mesh as entrainment separator device, Rd may be taken as 1.3.
A = V / [Rd 0.0172 {(l - v) / v}]
= [15480 /(3600 0.258)] / [1.3 0.0172 {(1019 0.258)/0.258}]
= 11.87 m2
Therefore drum diameter = {(4 11.87) /3.14}
3.89 m
Which is very large and therefore taking the drum diameter same as that of calendria.i.e 1.27 m
Drum height can be taken as 2 to 5 times of tube sheet diameter.
Thus drum height = 21.27
= 3.81 m
5) Drum thickness:
Drum is operating at 0.4163 bar. Design is therefore based on an external pressure of 1.7856 kg/cm2
(Design pressure).
Assume thickness of 15 mm.
The critical external pressure is given by
Pc = [2.42E (t /Do)2.5] / [(1 - 2)3/4 {(L /Do) 0.45 (t / Do)0.5}]
= [2.42 19 105 (10 /1290)2.5] / [(1 0.322)3/4{(3810 /1290) 0.45 (10 /1290)0.5}]
= 8.96 kg/cm2
Pa = Pc / 4
= 2.24 kg/cm2
According to IS 2825 (Appendix F)
L /Do =3.81/ 1.27
=3
Do / t = 1270 /10

= 85
Therefore factor B = 10878
Pa= B / [14.22 (Do / t)]
= 10878 / [14.22 85]
= 9 kg/cm2
Here Pa is greater than design pressure.
Assume thickness as 15 mm.
6) Flange calculation:
Flange material = IS:2004 1962 class 2
Bolting steel = 5% Cr Mo steel
Gasket material = asbestos composition
Outside diameter of calendria = 1290 mm
Calendria sheet thickness = 10 mm
Inside diameter of calendria = 1270 mm
Allowable stress of flange material = 100 MN/m2
Allowable stress for bolting material = 138 MN/m2
Determination of gasket width
do / di = [(y Pm) /{y P(m + 1)}]0.5
Assuming gasket thickness of 3.2 mm
Therefore y = 11.0, m = 2, P = 0.11 MN/m2
do / di = [(11 0.11 2) /{11 0.11 (2 + 1)}]0.5
= 1.005
Let di of gasket equals 1300 mm i.e. 10 mm larger than calendria diameter, then
do= 1.3065 m
Minimum gasket width 2b = (1.3065 1.300) /2
= 3.2510-3 m
Basic gasket seating width, bo = 12 / 2
= 6.0 mm
Diameter at location of gasket load reaction is, G = di + N
= 1300 + (3.2510-3)
= 1300 mm

Estimation of bolt loads

Load due to design pressure
H = (/4)G2 P
= (/4) (1.3)20.11
= 0.15 MN
Load to keep joint tight under operation
Hp = G(2b)mP
= 1.3(3.2510-3 ) 2 0.11
= 2.92 103 MN
Total operating load Wo = H + Hp
= 0.15+2.92 103
= 0.153 MN
Load to seat gasket under bolting up condition Wg = Gby
= 3.046 (3.2510-3/2) 11
= 0.17 MN
Here Wg is larger than Wo and therefore, controlling load = 0.17 MN
Calculation of minimum bolting area:
Am = Ao = Wg/So = 0.17/138
= 1.23 10-3 m2
Consider the bolt size as M 20 2
Therefore root area = 210-4 m2
Number of bolts required = (1.23 10-3) / (2 10-4)
=8
Flange thickness:
An approximate value of flange thickness may be given by
tF = G{P/(kf)}
Where,
k = 1/ [0.3 +{(1.5WmhG) /(HG)}]
Wm - total bolt load = 1.2894 MN
hG = (BG) / 2

Where B is minimum pitch circle diameter.

B=G +12 +(220)
= 1300+12 + 40
= 1352 m
Therefore hG = (1352 1300) / 2
= 26
H = 1.25 MN
= 127421 kg
G = 1300 mm
Therefore k =1 / [0.3 + {(1.5 131437.31 26) / (127421 1300)}]
= 3.02
Hence tf = 1300{1.103/(3.02 1019)}
= 24.61 mm
Therefore use thickness of 30 mm including corrosion allowance.
7) Bracket design:
Data:
Diameter of vessel = 1270 mm
Height of vessel = 3.81 m
Clearance from vessel bottom to foundation = 1000 mm (assumed)
Density of carbon steel = 7820 kg/m3
Density of brass = 8450 kg/m3
Wind pressure =128.5 kg/m2
Number of brackets = 4
Diameter of bolt circle = 3.15 m
Height of bracket from foundation = 2.25 m
Permissible stress for structural steel (IS 800)
Tension = 1400 kg/cm2
Compression = 1233 kg/cm2
Bending = 1575 kg/cm2
To find weight of vessel with contents:
Weight of vapour drum = dLt

= 1.27 3.81 0.015 7820

W1= 1783 kg
Weight of tubes W2 = (/4)Nt(do2 di2) L
= (/4) 619 (0.0422 2 0.032462) 3.048 8450
= 9150 kg
weight of tube sheet W3 = 2 (/4)Ds2 t
= (/4) 1.272 0.025 7820 2
= 495 kg
Therefore total weight W = W1 + W2 + W3
= 10728 kg
a) Base plate:
Taking suitable base plate size, a = 140 mm, B = 150 mm
Maximum total compressive load is given by
P = [{4Pw(H f)} / (nDb)] + [W / n]
Where Pw - total force acting on vessel due to wind = kPhDo
k = 0.7
P - wind pressure = 128.5 kg/m2
h - height = 3.81 m
Thus Pw= 0.7 128.5 3.81 1.27
= 435.2 kg/m2
H - height of vessel above foundation =1.27 m
f - vessel clearance from foundation = 0
n - number of brackets = 4
Db diameter of bolt circle = 1.65 m
Hence compressive load,
P = [{4 435.2 (1.27 0)}/ (4 3.15)] + [10728/ 4]
= 2857 kg
Average pressure on the plate:
Pav = P / (a B)
= 2857 / (14 15)
=13.6 Kg/cm2

But f = 0.7 Pav (B2 / T12) {a4 / (a4 + B4)}

= 0.7 13.6 (152 / T12) {144 (144 + 154)}
= 997.2 / T12
Therefore T12 = 997.2 / 1575
= 0.63mm
T1 = 0.8 mm
Use 6 mm thick plate.
b) Web plate:
Bending moment of each plate = (19532.25 / 2) {(3.15 3) / 2} 100
= 73246 kg cm
Stress at the edge = f = (3 P C) / (T2 h2)
= (73246 / 0.707) / (T2 14 14)
1575 = 528.58 / T2
Therefore T2 = 0.3356 cm = 3.356 mm
T2 may be taken as 4 to 6 mm.
c) Column support:
It is proposed to use a channel section as column.
Size 150 75
Area of cross section = A = 20.88 cm2
Modulus of section = Zyy = 19.4 cm2
Radius of gyration = ryy = 2.21 cm
Weight of section = 16.4 kg/m
Height from foundation = l = 2.25 m
Equivalent length for fixed ends = le = l 2 = 2.25 2 = 1.125 m
Slenderness ratio = le r = (1.125 100) 2.21 = 51.0
Now if the load is acting eccentric on a short column, the maximum combined bending and direct
stress is given by
f = [W (A n)] + [(W e) (n Z)]
= [71861 (20.88 4)] + [(71861 7.5) (4 19.4)]
= 7805.73 kg/cm2

The permissible compressive stress is

f = [W (A n)] [1 + a (le r)2] + [(W e) (n Z)]
= [71861 (20.88 4)] [1 + (512 7500)] + [(71861 7.5) (4 19.4)]
= 8104.12 kg/cm2
The calculated values are less than the permissible compressive stress and hence the channel selected
is satisfactory.
d) Base plate for column:
The size of column is 150 75. It is assumed that the base plate extends 20 mm on either side of the
channel.
Width = 0.8 75 + 2 20
= 100 mm
Length = 0.95 150 + 2 20
= 182.5 mm
Bearing pressure on each plate Pb = P /(B C)
= 19532.25 / (10 18.25)
= 107.03 kg/cm2
This is less than the permissible bearing pressure for concrete.
Stress in the plate = [(107.03 / 2) (202 /102)] / (t2 / 6)
= 12.84 / t2 kg/cm2
But f = 1575 kg/cm2
Therefore t2 = (12.84 / 1575) 100 mm2
t = 0.8152 mm
It is usual to select a plate of 4 to 6 mm thickness.

CONDENSER:

Design:The condenser is a horizontal condenser designed to condense 7652 Kg/Hr of 98% vapour of
glycerine at 1 atm pr. and 190C .The coolant used is water which is supplied in the tube side at an
inlet temperature of 20C and leaves at an outlet temperature of 35C.
Heat of vapourisation of glycerine = 18170 cal./mole
= 18170 4.18 92 J/Kg
= 6987.46 KJ/Kg
Amount of heat removed from the vapour Q = m
= (7652/3600) 6987.461000
= 14852.23 KJ/sec
AMOUNT OF WATER TO BE CIRCULATED:
(m)W Cp t = (m)G
(m)W 4.187 (35-20) = 14852.23 KW
(m)W = (14852.23 1000) / (4.18715)
(m)W = 236.48 Kg/Sec
Amount of water required = 236.48 Kg/sec
LOGARITHMIC MEAN TEMPERATURE DIFFERENCE:
LMTD = ((190-20)-(190-35))/(ln ((190-20)/(190-35)))
= (170 - 155 )/ (ln (170/155))
= 162.36C
OVERALL HEAT TRANSFER COEFFICIENT:
Assume U = 600 W/m2k
Total heat transfer area A = Q/(U (t)LMTD)
= 14852.23/(600 162.38)
= 152.45 m2
Choose tubes of 5/4 O.D.,16 BWG ,length of 16 ft laid on a 25/16 square pitch.
Heat transfer area A = 152.45 m2
Number of tubes Nt = 152.45/( 1.25 0.0254 (4.88-0.05))
= 317
From tube count table,
For tube O.D. of 5/4" on 25/16" square pitch

TEMA P or S
No. of passes = 2
Nearest no. of tubes Nt = 310
I.D. of shell = 889 mm
Corrected heat transfer area = n d L
=310 3.14 0.03175 4.88
= 150 m2
Corrected U = (14852.23 1000)/(150 162.38)
= 614 W/m2k
SHELL SIDE HEAT TRANSFER COEFFICIENT
CONDENSING VAPOUR SIDE:
Temperature of vapour coming in = 190C
Average temperature of water = (20 + 35 )/2
= 27.5C
Wall temperature = (190 + 27.5 )/2
=108.75C
Film temperature = (190 + 108.75)/2
=149.38C
So the property of water are taken at 149C
viscosity = 1.6103 Ns/m2
Thermal conductivity K = 0.670 W/mC
Density =1175 Kg/m3
Specific heat Cp = 3.05 KJ/KgC
Mass flow rate per unit length = W/(Nt2/3 L)
=7652/(36003102/34.88)
=9.568103 Kg/m sec
Renolds number Nre = (49.568103)/(1.6103)
= 23.92
Outer film coefficient hc =1.51(K32g / 2)2/3(Nre1/3)
hc =1.51(0.300813806259.81/2.56106)2/3(23.921/3)
= 1808.3 W/m2k

TUBE SIDE HEAT TRANSFER COEFFICIENT:

Average temperature of water =27.5C
Physical properties of water at 27.5C
Specific heat Cp 4.18 KJ/KgC
Viscosity = 0.9103 Ns/m2
Thermal conductivity K = 0.616 W/mC
Prandtl number Pr = (Cp )/K
= 6.41
Mass velocity of water = 236.48 Kg/sec
at = (n0.0284)/(42)
= (3103.140.0284)/8
= 0.0492 m2/pass
Gt = mass velocity/at
= 236.48/0.0492
= 2806.5 Kg /m2sec
Inside diameter of the tube di = 1.12
= 0.0284 m
NRe = (di Gt)/
= (0.02844806.5437)/(0.9103)
= 15192 > 10000
Ditus boltern equation can be used
Nu = 0.023 (NRe)0.8 (Pr)0.3
= 0.023151920.86.410.3
= 561.17
hi = 4518 w/m2k
CALCULATED OVERALL HEAT TRANFER COEFFICIENT:
1/Uo = 1/ho + 1/hi(do/di) +0.000528
= 1/1808.3 + 1/4518(1.25/1.12) +0.000528
= 0.000553+0.000247+0.000528
=0.001328
Uo = 753 W/m2k

Which is greater then the corrected U

Therefore this value of U is good enough.
PRESSURE DROP CALCULATIONS:
TUBE SIDE:
Nre = 15192
f = 0.079/(Nre)1/4
= 0.004
PL=(4fLv2/2gdi)g
vt= Gt/
= 4.82m/sec
PL= (40.0044.884.822/20.0284)997
=15.920 Kpa
PC = (2.5/2)vt2
=14.476 Kpa
PT = 2(PL+PC)
= 2(15.920+14.476)
= 60.8Kpa
This is within the permissible limit of a maximum pressure drop of 70Kpa in the tube side
So this pressure drop is acceptable.
SHELL SIDE:
Mass flow rate of glycerin =7652/3600
=2.13 Kg/sec
Saturation temperature of vapour Tvap=190C
Clearance C = pitch O.D.
= 25/16-5/4
= 7.94 mm
Pitch Pt = 25/16
= 0.04 m
as = (I.D)CB/Pt
Here I.D.=B=Ds-diameter of shell = 889 mm
as = 0.8890.8897.94103/0.04

= 0.157 m2
Equivalent diameter De= 4[Pt2-(do2)/4]/( do)
= 4[0.042-(3.140.031752)/4]/(3.140.03175)
= 0.0324m
Gs = 2.13/0.157
=13.57 kg/sec m2
At 190C vapor viscosity vap =1.1105 Ns/m2
(Nre)vap =GsDe/ vap
=13.570.0324105/1.1
= 39970
f =1.87(Nre0.2)
=0.225
Density of vapour vap= 2.42 Kg/m3
Number of baffles Nb+1=L/Ds
=4.88/0.889
=5.5 i.e. 6
Therefore Nb =5
PS =2[(f(Nb+1)DsGs2g)/(gDevap)]0.5
=(0.22560.88913.572)/(0.03242.42)
=3194.4 pa
=3.194 Kpa
This is also within the permissible limit of a maximum pressure drop of 14 Kpa
So this is acceptable.

MECHANICAL DESIGN OF CONDENSER:

SHELL SIDE:

Material

: Carbon steel

No. of shells

:1

No. of passes

:2

Fluid

: 80 % Glycerine vapour

Internal diameter : 889mm

Working pressure : 0.1 N/mm2
design pressure

: 0.11 N/mm2

Inlet temperature : 190C

Outlet temperature: 190C
Allowable stress

: 950 Kg/Cm2

TUBE SIDE:
Material

: Stainless steel (IS grade 10)

No. of tubes

: 310

Outside diameter : 31.75 mm

Length

: 4.88 m

fluid

: water

pitch

: 39.69 mm (square)

Allowable stress

: 10.06 Kg/m2

working pressure : 1.033 Kg/cm2

Design pressure

: 1.55 Kg/cm2

Inlet temperature : 20C

Outlet temperature : 35C
SHELL SIDE:

SHELL THICKNESS:
ts = (Pd DS)/(2fJ-Pd)
= (0.11889)/((2950.85)-0.11)
= 0.61 mm
But minimum thickness of shell is 6 mm
Therefore with corrosion allowance of 2mm
Thickness of shell = 8 mm
NOZZLE DIAMETER
M = Mass velocity/sec
= 7652/3600
= 2.13 Kg/sec
Density =1175 kg/m3
Assume velocity to be 10 m/sec
( dn2 v)/4 = M
dn2= (2.134)/(103.141175)
dn =0.015 m
NOZZLE THICKNESS
tn = (Pddn)/(2fJ-P)
= (0.1115)/(2950.85-0.11)
= 0.10 mm
Nozzle thickness with corrosion allowance = 5 mm
HEAD THICKNESS:
th = (PdRcW)/(2fJ)
W = (1/4)(3+(Rc/RK)1/2 )
Rc - crown radius 80% of shell I.D. = 711.2 mm
RK - Knuckle radius 10% of shell I.D.= 88.9 mm
W = 1.46
th = 0.71 mm
Using same thickness as that of the shell = 8 mm

BAFFLE ARRANGEMENT:

Tansverse bafffles
Number of baffles =5
Baffle Spacing = Ds = 889 mm
Thickness of baffles = 6 mm
Height of baffle = 0.75 Ds
= 666 mm
TIE RODS AND SPACERS:
For shell diameter Ds = 889 mm
No. of tie rods = 6
Diameter of rods = 13 mm
FLANGE CALCULATION:
Flange material = IS:2004-1962 class 2
Bolting steel =5% Cr Mo steel
Gasket material =asbestos composition
Shell inside diameter =889 mm
Shell thickness ts = 8mm
Shell outside diameter = (2ts)+889
= (28) +889
= 905 mm
Allowable stress of flange material = 100 MN/m2
Allowable stress of bolting material = 138 MN/m2
GASKET WIDTH:
Go/Gi =[(y-Pdm)/(y-pd(m+1))]1/2
m-Gasket factor =2.75
y-Minimum design seating stress =25.5 MN/m2
Go/Gi = [25.5-(0.112.75)/25.5-(0.113.75)]1/2
= 1.0022
Minimum gasket width N =10 mm
Basic gasket seating width bo = N/2
= 5 mm< 6.3 mm
Let,inside diameter of gasket = inside diameter of shell = 889 mm

Gi = 889 +(28)
= 905 mm
Mean gasket width = Gi+N
= 915 mm
therefore G = 915 mm
Estimation of bolt load:
Load due to design pressure H = (G2Pd)/4
= (0.91520.11)/4
= 0.072 MN
Effective gasket sitting width b = bo= 5 mm since b<6.3 mm
Load to keep joint tight under pressure Hp = (2b)GmPd
= 3.140.010.9152.750.11
= 0.0087 MN
Total operating load Wo =H+Hp
=0.072+0.0087
= 0.0807 MN
Load to seat gasket under bolting up condition Wb = bGPd
=3.140.9150.0050.11
= 0.367 MN
Since Wb>Wo, controlling load =0.367 MN
Minimum bolting area:
Total cross sectional area of bolt under operating condition Am1=Wo/Sb
Sb-nominal bolt stress at design temperature of 190C=138 MN/m2
Therefore Am1 = 0.0807/138
= 0.000585 m2
Total cross sectional area of bolt required for gasket seating Am2 =Wb/Sa
Sa-nominal bolt stress at ambient temperature (30C)=138 MN/m2
Therefore Am2 = 0.367/138
= 0.00266 m2
Since Am2 > Am1 ,Am = Am2 = 0.00266 m2
Calculation of optimum bolt size:

C =2(R+g1)+B
Choosing bolt M-18 2
Total number of bolts = G/(18 2)
= 915/36
= 25
Actual number of bolts = 24
R-radial clearance from bolt circle to point of connection of hub and back of flange = 27 mm
B-inside diameter of flange = outside diameter of shell = 0.905 m
g1 =go/0.707,let go= 8 mm
g1 = 0.008
C = 2(0.027 + 0.008) + 0.905
= 0.975 m
Therefore bolt circle diameter = 0.975 m
Flange outside diameter:
A = C + bolt diameter + 0.02
= 0.975 + 0.018 + 0.02
= 1.013 m
Check of gasket width:
Ab-root area of bolt (m2)-1.54 10-4 m2
Sg-allowable stress for bolting material at atmospheric temperature =138 MN/ m2
Therefore, AbSg/GN = 15.86
Since 15.86 > 2y
Condition is satisfied.
Flange moment computations:
Wo = W1 + W2 + W3 (under operating condition)
W1 = (B2Pd)/4
= (3.140.90520.11)/4
= 0.071 MN
W2 = H W1
H =(G2Pd)/4
= (3.140.91520.11)/4

= 0.0732 MN
W2 = 0.0723 0.071
= 0.0013 MN
W3 = Wo H
= 0.0807 0.0723
= 0.0084 MN
Total flange moment ,
Mo = W1a1 + W2a2 + W3a3
a1 = (C B)/2
= 0.035 m
a3 = (C-G)/2
= 0.03 m
a2 = (a1+a2)/2
= 0.0325 m
Mo =0.0710.035 + 0.00130.03 + 0.00840.0325
= 2.78 10-3 MNm
Bolting up condition:
Total flange moment Mg = Wa3
W = (Am + Ab)Sg/2
= 0.00472138
= 0.651 MN
a3 = 0.03 m
Therefore Mg = 0.01953 MNm
= 19.5 10-3 MNm
Since Mg > Mo for moment under operating condition,Mg is controlling.
Therefore M = Mg = 19.510-3 MNm
Flange thickness:
t2 = (MCFY)/BST
K =A/B
= Outer diameter of flange/ inner diameter of flange
= 1.013/0.905

= 1.12
Y =14
Assume CF =1
Therefore thickness 't' = 0.055 m
Actual bolt spacing BS =C/n
= (3.140.975)/24
= 0.128 m
Bolt pitch correction factor CF =[BS/(2d+t)]1/2
=[0.128/(20.018+0.055)]1/2
=1.19
Therefore CF1/2 = 1.091
Actual flange thickness = CF1/2 t
= 1.0910.055
= 0.06 m
= 60 mm
TUBE SIDE:
TUBE THICKNESS:
tt =Pdo/2fJ+P
J =1 for seamless tube
Therefore tt =(1.5531.75)/(21006+1.55)
= 0.024 mm
No corrosion allowance since the tube is made of stainless steel
Thickness of tube = 1mm
TUBE SHEET:
tS = FG[0.25P/f]1/2
F-The value of F varies according to type of heat exchanger,for most cases it is taken as 1
G =915 mm
Therefore tS =915[(0.251.55)/1006]1/2
= 17.96 mm

CHANNEL AND CHANNEL COVER:

t = G[KP/f]1/2
K = 0.3 for ring type gasket
Material of construction is carbon steel
So allowable stress f =950 Kg/cm2
Therefore t =915[(0.31.55)/950]1/2
= 20.24 mm
With corrosion allowance t = 25 mm
NOZZLE THICKNESS:
Assume inlet and outlet diameter = 75 mm
Thickness of nozzle th = Pd/2fJ-P
= (1.5575)/(20.85950-1.55)
= 0.0612 mm
With corrosion allowance thickness = 6mm

SADDLE SUPPORT DESIGN:

Material -Low carbon steel
Vessel diameter=905 mm
Length of shell =4.88 m
Torispherical head:
Crown radius =711.2 mm
Knuckle radius =88.9 mm
Working pressure =1Atm
Shell thickness = 8 mm
Head thickness = 8mm
Corrosion allowance = 2mm
permissible stress = 950 Kg/cm2
R-Vessel radius = 425.5
Distance of saddle center line from shell end A = 0.45R
= 203.63 mm < 0.2L

Longitudinal bending moment:

The bending moment at the support is;

M1 =QA[1-{(1-A/L)+(R2-H2/2AL)/(1+4H/3L)}]
A =203.63 mm
Q =W/2[L+4H/3]
W-Weight of fluid and vessel.
Weight of shell material:
W1 =[(Do2-Di2)Lshell material] /4
shell material - 7700 Kg/m3
W1 =[3.14(0.9052-0.8892)4.887700] /4
= 847 Kg
Weight of tubes:
W2 = n[(Do2-Di2)Ltube material] /4
tube material -7800 Kg/m3
W2 =310[3.14(0.031752-0.02982)4.887800] /4
= 1112 Kg
Weight of tube sheet:
W3 = (2D2t)/4
= (23.140.90520.01807800)/4
= 180 Kg
Liquid load in the shell:
W4 =(shell volume -tube volume)liquid
=[(Ds2L)/4 - (ndo2l)/4 ]1175
=[(0.90524.88)/4 - (3100.0317524.88) /4 ]1175
=(3.14 - 1.20)1175
= 2280 Kg
Liquid load in tubes:
W5 = ndi2lliquid /4
=(3103.140.029824.88997)/4
= 1052 Kg
Therefore total weight WT = W1 + W2 + W3 + W4 + W5
= 847 + 1112 + 180 + 2280 + 1052

= 5471 Kg
Hence, Q =5471/2[4.88+(40.257)/3]
= 14287 Kgm
M1=(142870.20363)[1-{(1-0.20363/4.88)+(0.45252-0.2572/20.2034.88)/(1+40.257/34.88)}]
=2909.261[1-{(0.958+0.0698)/1.0702}]
=115.26 Kgm
The bending moment at the center of the span is given by:
M2 =(QL/4)[{1+2(R2-H2)/L2}/{1+(4H/3L)}- (4A/L)]
=(142874.88/4)[{1+2(0.45252-0.2572)/4.882}/{1+(40.257/34.88)}- (40.203/4.88)]
= 17430[1.312 - 0.1664]
= 19968 Kgm
Stress in shell at the saddle:
At the top most fiber of the cross-section,
f1 =M1/(K1R2t)
For an angle of 120 ,K1 =0.107 m
t-thickness of shell = 8 mm
f1 =115.6/(0.1073.140.452520.008)
=20.93 Kg/cm2
At the bottom most fiber of the cross-section
f2 =M1/(K2R2t)
For an angle of 120 ,K2 =0.192 m
f2 =115.6/(0.1923.140.452520.008)
=11.67 Kg/cm2
Stress in the shell at mid point:
f3 =M2/(R2t)
=19968/(3.140.452520.008)
=388 Kg/cm2
Thus the values of stresses are within the limited range
Hence the designed support is acceptable.
COST ESTIMATION:
Glycerin plant size =100 T/day

Fixed capital investment for cost index of 130 = Rs 108

Cost index for 2002 = 402
Therefore present fixed capital investment =108(402/130)
=Rs 30,92,30,769
Estimation of total investment cost:
1) Direct cost:
a) Purchased equipment cost:(15 40% of FCI )
Assume 40% of FCI
=Rs 12,36,92,307
b) Installation cost:(35 45% of PC)
Assume 45%
=Rs 13,91,53,846
c) Instrument and control installed:(6 30% of PEC)
Assume 30% of PEC
=Rs 9,27,69,230
d) Piping installation cost:(10 80% of PEC)
Assume 75%
=Rs.23,19,23,076
e) Electrical installation cost:(10 40% of PEC)
Assume 40% of PEC
=Rs 12,36,92,307
f) Building process and auxilliary:(10-70% of PEC)
Assume 65%
=Rs 20,09,99,999
g) Service facilities:(30-80% 0f PEC)
Assume 75%
=Rs 23,19,23,076
h) Yard improvement:(10-15% of PEC)
Assume 15%
=Rs 4,63,84,615
i) Land:(4-8% of PEC)

Assume 8%
=Rs 2,47,38,461
Therefore direct cost =Rs 1,21,52,76,917
2) Indirect cost:
Expenses which are not directly involved with material and labour of actual installation or complete
facility
a) Engineering and supervision:(5-30% of DC)
Assume 30%
=Rs 36,45,83,075
a) Construction expenses:(10% of DC)
=Rs 12,15,27,691
b) Contractors fee:(2-7% 0f DC)
Assume 7%
=Rs 8,50,69,584
c) Contingency:(8-20% of DC)
Assume 10%
=Rs 12,15,27,691
Therefore total indirect cost =69,27,08,042
3) Fixed capital investment:
Fixed capital investment(FCI) = DC+IC
= Rs 1,90,79,84,959
4) Working capital investment:
10 20% of FCI
Assume 15%
=Rs 28,61,97,743
5) Total capital investment:
= FCI + WC
=Rs 2,19,41,82,703

Estimation of total product cost(TPC):

Fixed charges:

a) Depreciation:(10% of FCI for machinery)

=Rs 2,19,41,827
b) Local taxes:(3-4% of FCI)
Assume 3%
=Rs 6,58,25,481
c) Insurances:(0.4-1% of FCI)
Assume 1%
=Rs 2,19,41,827
d) Rent:(8-12% of FCI)
Assume 9%
=Rs 4,57,97,372
Therefore total fixed charges =Rs 15,55,06,507
But, Fixed charges = (10-20% of TPC)
Assume 10%
Therefore Total product cost =155506507/0.1
=Rs 1,55,50,65,070
Direct production:
a) Raw material:(10-50% 0f TPC)
Assume 30%
=Rs 46,65,19,521
b) Operating labour(OL):(10-20% of TPC)
Assume 15%
=Rs 23,32,59,760
c) Direct supervisory and electric labour:(10-25% of OL)
Assume 15%
=Rs 3,49,88,964
d) Utilities:(10-20% of TPC)
Assume 15%
=Rs 23,32,59,760
e) Maintainence:(2-10% of FCI)
Assume 6%

=Rs 9,33,03,904
e) Operating supplies (OS):(10-20% of maintainence)
Assume 15%
=Rs 1,39,95,585
f) Laboratory charges:(10-20% of OL)
Assume 15%
=Rs 3,49,88,964
g) Patent and royalties:(2-6% of TPC)
Assume 4%
=Rs 6,22,02,602
Plant overhead cost:
50-70% of (OL+OS+M)
Assume 60%
=Rs 20,43,35,549
General expenses:
a) Administration cost:(40-60% of OL)
Assume 55%
=Rs 12,82,92,868
b) Distribution and selling price:(2-30% of TPC)
Assume 15%
=Rs 23,32,59,760
c) Research and development cost:(3% of TPC)
=Rs 4,66,51,952
Therefore general expenses(GE) =Rs 40,82,04,580
Therefore manufaacturing cost(MC)= Product cost+fixed charges+Plant overhead expenses
= Rs 1,91,49,07,126
Total production cost:
Total production cost =MC + GE
=Rs 2,32,31,11,706
Gross earnings and rate of return:
The plant is working for say 300 days a year

Selling price =Rs. 97/kg

Total income =100300100097
=Rs 2,91,00,00,000
Gross income =Total income total product cost(1555065070)
=Rs 1,35,49,34,930
Tax =50%
Net profit =Rs 67,74,67,465
Rate of return = net profit/total capital investment
=30.87%