DESIGN OF EQUIPMENT

5.1 PROCESS DESIGN OF DISTILLATION COLUMN:

5.1.1 Glossary of notations used:

F = molar flow rate of Feed, kmol/hr.
D = molar flow rate of Distillate, kmol/hr.
W = molar flow rate of Residue, kmol/hr.
xF = mole fraction of Acetaldehyde in liquid/Feed.
yD = mole fraction of Acetaldehyde in Distillate.
xW = mole fraction of Acetaldehyde in Residue.
MF = Average Molecular weight of Feed, kg/kmol
MD = Average Molecular weight of Distillate, kg/kmol
MW = Average Molecular weight of Residue, kg/kmol
Rm = Minimum Reflux ratio
R = Actual Reflux ratio
L = Molar flow rate of Liquid in the Enriching Section, kmol/hr.
G = Molar flow rate of Vapor in the Enriching Section, kmol/hr.
L = Molar flow rate of Liquid in Stripping Section, kmol/hr.
G = Molar flow rate of Vapor in Stripping Section, kmol/hr.
q = Thermal condition of Feed
L = Density of Liquid, kg/m3.
V = Density of Vapor, kg/m3.
qL = Volumetric flow rate of Liquid, m3/s
qV = Volumetric flow rate of Vapor, m3/s
L = Viscosity of Liquid, cP.
TL = Temperature of Liquid, 0K.
TV = Temperature of Vapor, 0K

T x- y data:

T 0C

98.5

89.9

80

71

60.5

50

39

0.000

0.069

0.164

0.286

0.445

0.664

1.000

0.000

0.317

0.578

0.761

0.879

0.954

1.000

Table 5.1 T-x-y data.

5.1.2 Preliminary calculations:

F = 152.798 kmol/hr, xF = 0.938, MF = 44.123 kg/kmol.
D = 144.7 kmol/hr, xD = 0.99, MD = 44.04 kg/kmol.
W = 8.0931 kmol/hr, xW = 0.177, MW = 45.64 kg/kmol.
Distillation column temperature = 400 C.
Distillation column pressure = 2.08 atm. = 1586.41 mm Hg.
Basis: One-hour operation.
From the graph,
xD / (Rm+1) = 0.94
Thus, Rm = 0.0476
Let, R= 1.5*Rm
Therefore, R= 1.5*0.0476= 0.0714
Thus, xD/ (R+1) = 0.99/ (0.0714+1)
i.e., xD/ (R+1) = 0.924
Number of Ideal trays = 4 (including the reboiler).
Reboiler is the last tray.
Number of Ideal trays in Enriching Section = 2
Number of Ideal trays in Stripping Section = 2
Now, we know that,
R = Lo/ D
=> Lo = R*D
i.e., Lo= 0.0714*144.7
i.e., Lo =10.33 kmol/hr.
Therefore, Lo = 10.33 kmol/hr.
L= Liquid flow rate on the Top tray = 10.33 kmol/hr.

Since feed is Liquid, entering at bubble point,

q= (HV-HF) / (HV-HL) = 1
Now,
Slope of q-line = q/ (q-1)
= 1/ (1-1) = 1/0 =
Now we know that,
(L -L) = q = 1
F
(L - L) = F
L=F+L
i.e., L = 10.33 + 152.798
i.e., L = 163.128 kmol/hr.
Therefore, liquid flow rate in the Stripping Section = 163.128 kmol/hr.
Also, we know that,
G = [(q-1) F] + G
i.e., G = [(1-1) F] + G
i.e., G = [0F] +G
i.e., G = 0 +G
G=G
Now, we know that,
G=L+D
i.e.,

G = Lo +D

i.e.,

G= 10.33 + 144.7

i.e.,

G= 155.03 kmol/hr.

Thus, the flow rate of Vapor in the Enriching Section = 155.03 kmol/hr.

Since G =G
G = G = 155.03 kmol/hr.
Therefore, the flow rate of Vapor in the Stripping Section = 155.03 kmol/hr.

5.1.3 List of parameters used in calculation:

SECTION

ENRICHING SECTION

STRIPPING SECTION

PROPERTY

TOP

BOTTOM

TOP

BOTTOM

0.99

0.95

0.95

0.177

0.99

0.97

0.97

0.177

L 10.33

10.33

163.128

163.128

G 155.03

155.03

155.03

155.03

Liquid,
kmol/hr.
Vapor,
kmol/hr.
T liquid, 0C

39.07

39.77

39.77

80.05

T vapor, 0C

53.00

54.01

54.01

94.13

liquid

44.02

44.1

44.1

45.646

Vapor

44.02

44.06

44.06

45.646

Liquid, L kg/hr.

454.726

455.55

7193,9

7446.14

Vapor, G kg/hr

6824.42

6830.02

6830.6

7076.5

!l 784.69

784.50

784.50

747.87

!g 3.4376

3.425

3.425

3.361

0.004

0.06958

0.0705

Mavg.
kg/kmol
Mavg.
kmol/hr

'HQVLW\
3

kg/m

'HQVLW\

kg/m3
/* !g 

!l) 0.5

0.0039

Table 5.2 Parameters used in calculations.

5.1.4 Design Specification:

a) Design of Enriching Section:

Tray Hydraulics,
The design of a sieve plate tower is described below. The equations and correlations are
borrowed from the 6th and 7th editions of Perrys Chemical Engineers Handbook.
1. Tray Spacing, (ts) :
Let ts = 18 = 457 mm. (range 0.15 1.0 m).
2. Hole Diameter, (dh):
Let dh = 5 mm. (range 2.5 12 mm).
3. Hole Pitch (lp):
Let lp = 3* dh (range 2.5 to 4.0 times dh).
i.e., lp = 3*5 = 15 mm.
4.

Tray thickness (tT):

Let tT = 0.6* dh (range 0.4 to 0.7 times dh).
i.e., tT = 0.6*5 = 3 mm.

5.

Ratio of hole area to perforated area (Ah/Ap):

Refer fig 3
Now, for a triangular pitch, we know that,
Ratio of hole area to perforated area (Ah/Ap) = (/4*dh2)/ [(3/4) *lp2]
i.e., (Ah/Ap) = 0.90* (dh/lp)2
i.e., (Ah/Ap) = 0.90* (5/15)2
i.e., (Ah/Ap) = 0.1
Thus,
(Ah/Ap) = 0.1

6.

Plate Diameter (Dc):

The plate diameter is calculated based on entrainment flooding considerations

L/G {g/l} 0.5 = 0.004

---------- (maximum value)

Now for,
L/G {g/l} 0.5 = 0.004 and for a tray spacing of 500 mm.

We have,
From the flooding curve, ---------- (fig.18.10, page 18.7, 6th edition Perry.)
Flooding parameter, Csb, flood = 0.29 ft/s .
Now,
Unf = Csb, flood * ( / 20) 0.2 [(l - g) / g]0.5
---- {eqn. 18.2, page 18.6, 6th edition Perry.}
Where,
Unf = gas velocity through the net area at flood, m/s (ft/s)
Csb, flood = capacity parameter, m/s (ft/s, as in fig.18.10)
= liquid surface tension, mN/m (dyne/cm.)
l = liquid density, kg/m3 (lb/ft3)
g = gas density, kg/m3 (lb/ft3)
Now, we have,
= 19.325 mN/m = 19.325 dyne/cm.
l = 784.5 kg/m3.
g = 3.425 kg/m3.
Therefore,
Unf = 0.29*(19.325/20)0.2*[(784.50-3.4250)/ 3.4250]0.5
i.e.,Unf = 4.349 ft/s = 1.325 m/s.
Let,
Actual velocity, Un= 0.8*Unf
i.e., Un = 0.81.325
i.e., Un = 1.06 m/s
It is desired to design with volumetric flow rate maximum (therefore, actual is less than the
maximum).
Volumetric flow rate of Vapor at the bottom of the Enriching Section
= qo = 6830.62 / (3600*3.4250) = 0.554 m3/s.
Now,
Net area available for gas flow (An)
Net area = (Column cross sectional area) - (Down comer area.)

An = Ac - Ad
Thus,
Net Active area, An = to/ Un = 0.554/ 1.06 = 0.522 m2.
Let Lw / Dc = 0.77 (range 0.6 to 0.85 times Dc ).
Where, Lw = weir length, m
Dc = Column diameter, m
Now,

,c = 2*sin-1(Lw / Dc) = 2*sin-1 (0.77) = 100.70

Now,
Ac

2
 'c =

0.785*Dc2 , m2

Ad = [(/4) * Dc2 * (c/3600)] - [(Lw/2) * (Dc/2) *cos (c/2)]

i.e., Ad = [0.785*Dc2 *(100.70/3600)]-[(1/4)* (Lw / Dc) * Dc2 * cos(100.70)]
i.e., Ad = (0.2196* Dc2) - (0.1288* Dc2)
i.e., Ad = 0.0968*Dc2, m2
Since

An = Ac -Ad
0.522 = (0.785*Dc2) - (0.0968* Dc2)

i.e., 0.6882* Dc2 = 0.522

Dc2 = 0.522/ 0.6882 = 0.7585
Dc = 0.7585
Dc = 0.87 m
Since Lw / Dc = 0.77,
Lw = 0.77* Dc = 0.77*0.87 = 0.67 m.
Therefore, Lw = 0.67 m.
Now,
Ac = 0.785*0.872 = 0.5944 m2
Ad = 0.0968*Dc2 = 0.0968*0.872 = 0.0724 m2
Aa = Ac 2* Ad
i.e., Aa = 0.5944- 2*0.0724
7.

Aa= 0.4496 m2

Perforated plate area (Ap):

Now,
Lw / Dc = 0.67/ 0.87 = 0.7701

,c = 100.73 0
.  0 - ,c
0
0
LH .
 - 100.73
.  0
Now,
Acz = 2* Lw* (thickness of distribution)
Where, Acz = area of calming zone, m2 (5 to 20% of Ac )
Acz = 2*0.67* (3010-3) = 0.0402 m2 -------- (which is 6.76% of Ac)
Also,
Awz

,c /360 0) -  (Dc 30*10-3)2 ,c /360 0)}

2
 ^ 'c 

Where, Awz = area of waste periphery, m2 (range 2 to 5% of Ac)

i.e., Awz

2
 ^  *

(100.73 0/360 0-  -30*10-3)2 * (100.73 0/360 0)}

i.e., Awz = 0.0225 m2 --------- (which is 3.8% of Ac)

Now,
Ap = Ac - (2*Ad) - Acz - Awz
i.e., Ap = 0.5944- (2*0.0724) - 0.0402 - 0.0225
Thus, Ap = 0.387 m2.

8. Total Hole Area (Ah):

Since,
Ah / Ap = 0.1
Ah = 0.1* Ap
i.e., Ah = 0.1*0.387
Ah = 0.0387 m2
Thus, Total Hole Area = 0.0387 m2
Now we know that,
Ah = nh  Gh2
Where, nh = number of holes.
nh = (4*Ah  Gh2)
i.e., nh
nh

   


Therefore, Number of holes = 1971.

9.

Weir Height (hw):

Let hw = 50 mm.

10.

Weeping Check
The static pressure below the tray should be capable enough to hold the liquid
above the tray so that no liquid sweeps through the holes.
All the pressure drops calculated in this section are represented as mm head of
liquid on the plate. This serves as a common basis for evaluating the pressure
drops.

Notations used and their units:

hd = Pressure drop through the dry plate, mm of liquid on the plate
uh = Vapor velocity based on the hole area, m/s
how = Height of liquid over weir, mm of liquid on the plate
h = Pressure drop due to bubble formation, mm of liquid
hds= Dynamic seal of liquid, mm of liquid
hl = Pressure drop due to foaming, mm of liquid
hf = Pressure drop due to foaming, actual, mm of liquid
Df = Average flow length of the liquid, m
Rh = Hydraulic radius of liquid flow, m
uf = Velocity of foam, m/s
(NRe) = Reynolds number of flow
f = Friction factor
hhg = Hydraulic gradient, mm of liquid
hda = Loss under down comer apron, mm of liquid
Ads = Area under the down comer apron, m2
c = Down comer clearance, m
hdc = Down comer backup, mm of liquid

Calculations:
Head loss through dry hole
hd = head loss across the dry hole
hd = k1 + [k2* (g/l) *Uh2] --------- (eqn. 18.6, page 18.9, 6th edition Perry)
Where, Uh =gas velocity through hole area
k1, k2 are constants
For sieve plates,
k1 = 0

and

k2 = 50.8 / (Cv)2
Where, Cv=discharge coefficient, taken from fig 18.14, page 18.9 6th edition Perry.
Now,
(Ah /Aa) = 0.0387/ 0.4496 = 0.086
Also,

tT/dh = 3/5 = 0.60

Thus for (Ah/Aa) = 0.086 and tT/dh = 0.60

We have from fig. edition 18.14, page 18.9 6th Perry.
Cv = 0.74

k2 = 50.8 / 0.742 = 92.77

Volumetric flow rate of Vapor at the top of the Enriching Section

=qt = 1.8956/ (3.4376) = 0.5514 m3/s -------- (minimum at top)
Volumetric flow rate of Vapor at the bottom of the Enriching Section
= qo = 1.897 / (3.425) = 0.554 m3/s. ---- (maximum at bottom)
Velocity through the hole area (Uh):
Now,
Velocity through the hole area at the top = Uh, top = qt /Ah
= 0.5514/0.0387= 14.25 m/s
Also, Velocity through the hole area at the bottom= Uh, bottom = qo /Ah
= 0.554/0.0387 = 14.31 m/s
Now,
hd, top = k2 [g/l] (Uh,top)2
= 92.77(3.4376/784.69) 14.252

 hd, top = 82.526 mm clear liquid. -------- (minimum at top)

Also,
hd, bottom = k2 [g/l] (Uh, bottom)2
= 92.77(3.425/784.50)14.312
hd, bottom = 82.94 mm clear liquid ----- (maximum at bottom)
Head Loss Due to Bubble Formation
h = 409 [ / ( Ldh) ]
where =surface tension, mN/m (dyne/cm) = 19.325 dyne/cm.
dh =Hole diameter, mm
l = density of liquid in the section, kg/m3
= 784.69 kg/m3
h = 409 [ 19.325/(784.69 *5)]
h = 2.014 mm clear liquid
Height of Liquid Crest over Weir:
how = 664Fw [(q/Lw)2/3]
q = liquid flow rate at top, m3/s
= 0.1263*60/ (784.69)
= 0.009 m3/min.
Thus, q = 2.377 gal/min.
Lw = weir length = 0.67 m = 2.198 ft
Now,
q/Lw2.5 = 2.377/ (2.198)2.5 = 0.3318
Now for q/Lw2.5 = 0.3318 and Lw /Dc =0.7701
We have from fig.18.16, page 18.11, 6th edition Perry
Fw= correction factor =1.03
Thus, how = 1.03664 [0.00015/0.67] 2/3
how = 2.52 mm clear liquid.
Now,

(hd + h) = 82.526 + 2.014 = 84.54 mm ------ Design value

(hw + how) = 50 + 2.52 = 52.52 mm
For, Ah/Aa = 0.086 and (hw + how) = 52.52 mm
The minimum value of (hd + h) required is calculated from a graph given in Perry,
plotted against Ah/Aa.
i.e., we have from fig. 18.11, page 18.7, 6th edition Perry
(hd + h)min = 13.0 mm ------- Theoretical value.
The minimum value as found is 13.0 mm.
Since the design value is greater than the minimum value, there is no problem of
weeping.

Down comer Flooding:

hds =hw + how + (hhg /2) ------- (eqn 18.10, page 18.10, 6th edition Perry)
Where,
hw = weir height, mm
hds = static slot seal (weir height minus height of top of slot above plate floor,
height equivalent clear liquid, mm)
how = height of crest over weir, equivalent clear liquid, mm
hhg = hydraulic gradient across the plate, height of equivalent clear liquid, mm.
Hydraulic gradient, hhg
Let hhg = 0.5 mm.
hds = hw + how + hhg/2
= 50 + 2.52 + 0.5/2 = 52.77 mm.
Now, Fga = Ua g0.5
Where Fga = gas-phase kinetic energy factor,
Ua = superficial gas velocity, m/s (ft/s),
g = gas density, kg/m3 (lb/ft3)
Here Ua is calculated at the bottom of the section.
Thus, Ua = (Gb/g)/ Aa = 1.8974/(3.425 * 0.4496) = 1.232 m/s

Thus, Ua = 4.042 ft/s

g = 3.4250 kg/m3 = 0.209 lb/ft3
Therefore, Fga = 4.042(0.209) 0.5
Fga = 1.848
Now for Fga = 1.848, we have from fig. 18.15, page 18.10 6th edition Perry
Aeration factor = = 0.6
Relative Froth Density = t = 0.2
Now hl= hds ---- (eqn. 18.8, page 18.10, 6th edition Perry)
Where, hl= pressure drop through the aerated mass over and around the disperser,
mm liquid,
hl= 0.652.77 = 31.662 mm.
Now,
hf = hl/t ------- (eqn. 18.9, page 18.10, 6th edition Perry)
hf = 31.662/ 0.2 = 158.31 mm.
Average width of liquid flow path, Df = (Dc + Lw)/2
= (0.87 + 0.67)/2 = 0.77 m.
Hydraulic radius of aerated mass Rh = hf * Df /(2*hf + 1000*Df) (from eq. 18.23, page
18.12 6th edition Perry)
Rh = 158.31*0.77/(2*158.31 + 1000*0.77)
= 0.112 m.
Velocity of aerated mass, Uf = 1000*q/ (hl * Df )
Volumetric flow rate, q = 1.6061*10-4 m3/s.
Uf = 1000* 1.6061*10-4 / (31.662* 0.77)
= 0.0066 m/s.
Reynolds modulus NRe = Rh * Uf * l / liq
= 0.112 * 0.0066 * 784.5 /(1.03 * 10-3)
= 563.012
hhg = 1000* f* Uf2 *Lf/(g * Rh)
f = 0.6 for hw = 1.97 and NRe = 563. 012
Lf = 2 * Dc FRV ,c / 2) = 0.5549 m
hhg = 1000* 0.6 *0.00662*0.5549/(9.81* 0.112)

= 0.0132 mm.

Head loss over down comer apron:

hda = 165.2 {q/ Ada}2 ----- (eqn. 18.19, page 18.10, 6th edition Perry)
Where, hda = head loss under the down comer apron, as millimeters of liquid,
q = liquid flow rate calculated at the bottom of section, m3/s
And Ada = minimum area of flow under the down comer apron, m2
Now,
q = 1.60610-4 m3/s
Take clearance, C = 1 = 25.4 mm
hap = hds - C = 52.77 - 25.4 = 27.37 mm
Ada = Lw * hap = 0.6727.3710-3 = 0.0183 m2
hda = 165.2[(1.6061* 10-4)/ (0.0183)] 2
hda = 0.0127 mm
Now,
ht = hd + hl`
Here hd and hl are calculated at bottom of the enriching section.
Now we have,
hd, bottom = 82.94 mm
hl, bottom = 31.662 mm
ht = hd + hl`
= 82.94+31.662
ht = 114.602 mm
Down comer Backup:
hdc = ht + hw + how + hda +hhg ---- (eqn 18.3, page 18.7, 6th edition Perry)
ht = total pressure drop across the plate (mm liquid)
= hd + hl`
hdc = height in down comer, mm liquid,
hw = height of weir at the plate outlet, mm liquid,
how =height of crest over the weir, mm liquid,
hda = head loss due to liquid flow under the down comer apron, mm liquid,

hhg = liquid gradient across the plate, mm liquid.

hdc = 114.602 +50 +2.52 + 0.0132 + 0.0127
hdc = 167.148 mm.
Let dc = average relative froth density (ratio of froth density to liquid density)
=0.5
h`dc = hdc / dc = 167.148/ 0.5
h`dc = 334.29 mm.
which is less than the tray spacing, ts= 457 mm.
Hence no flooding in the enriching section and hence the design calculations are
acceptable.

b). Design of Stripping Section:

Tray Hydraulics
The design of a sieve plate tower is described below. The equations and correlations are
borrowed from the 6th and 7th editions of Perrys Chemical Engineers Handbook.
1

Tray Spacing, (ts) :

Let ts = 18 = 457 mm.

Hole Diameter, (dh):

Let dh = 5 mm.

Hole Pitch (lp):

Let lp = 3*dh
i.e., lp = 3*5 = 15 mm.

Tray thickness (tT):

Let tT = 0.6* dh
i.e., tT = 0.6*5 = 3 mm.

Ratio of hole area to perforated area (Ah/Ap):

Refer fig 3
Now, for a triangular pitch, we know that,
Ratio of hole area to perforated area (Ah/Ap) = (p/4*dh2)/ [(3/4) *lp2]

i.e., (Ah/Ap) = 0.90* (dh/lp)2

i.e., (Ah/Ap) = 0.90* (5/15)2
i.e., (Ah/Ap) = 0.1
Thus,
(Ah/Ap) = 0.1
6

Plate Diameter (Dc):

The plate diameter is calculated based on the flooding considerations
L/G {g/l}0.5 = 0.0705

---------- (maximum value)

Now for,
L/G {g/l}0.5 = 0.0705 and for a tray spacing of 457 mm.
We have,
From the flooding curve, ---------- (fig.18.10, page 18.7, 6th edition
Perry.)
Flooding parameter, Csb, flood = 0.27 ft/s.
Now,
Unf = Csb, flood * ( / 20) 0.2 [(l - g) / g]0.5
---- {eqn. 18.2, page 18.6, 6th edition Perry.}
where,
Unf = gas velocity through the net area at flood, m/s (ft/s)
Csb, flood = capacity parameter, m/s (ft/s, as in fig.18.10)
= liquid surface tension, mN/m (dyne/cm.)
l = liquid density, kg/m3 (lb/ft3)
g = gas density, kg/m3 (lb/ft3)
Now, we have,
= 18.330 mN/m = 18.330 dyne/cm.
l = 747.87 kg/m3.
g = 3.361 kg/m3.
Therefore,
Unf = 0.27* (18.33/20) 0.2 [(747.87-3.361)/ 3.361] 0.5
i.e., Unf = 3.949 ft/s

Let,
Actual velocity, Un= 0.8*Unf
i.e., Un = 0.83.949
i.e., Un = 3.159 ft/s
Un = 0.9628 m/s
Now,
Volumetric flow rate of Vapor at the bottom of the Stripping Section
= qo =1.9657/ (3.361) = 0.5848 m3/s.
Now,
Net area available for gas flow (An)
Net area = (Column cross sectional area) - (Down comer area.)
An = Ac - Ad
Thus,
Net Active area, An = qo/ Un = 0.5848/ 0.9628 = 0.6074 m2.
Let Lw / Dc = 0.77
Where, Lw = weir length, m
Dc = Column diameter, m
Now,

,c = 2*sin-1(Lw / Dc) = 2*sin-1 (0.77) = 100.70

Now,
Ac

2
 'c =

0.785*Dc2, m2

Ad = [(/4) * Dc2 * (c/3600)] - [(Lw/2) * (Dc/2) *cos (c/2)]

i.e., Ad = [0.7854* Dc2 * (100.70/3600)]-[(1/4) * (Lw / Dc) * Dc2 * cos (100.70/2)]
i.e., Ad = (0.2196* Dc2) - (0.1288* Dc2)
i.e., Ad = 0.0968*Dc2, m2
Since An = Ac -Ad
0.6882 = (0.785*Dc2) - (0.0968* Dc2)
i.e., 0.6882* Dc2 = 0.6074
Dc2 = 0.6074/ 0.6882 = 0.8826
Dc = 0.8826
Dc = 0.94 m

Therefore, Dc = 0.94 m
Since Lw / Dc = 0.77
Lw = 0.77* Dc = 0.77*0.94 = 0.724 m.
Therefore, Lw = 0.724 m.
Now,
Ac = 0.785*0.942 = 0.694 m2
Ad = 0.09688*Dc2 = 0.0968*0.942 = 0.0866 m2
An = Ac - Ad
i.e., An = 0.694 - 0.0866

An = 0.6074 m2

7 Perforated plate area (Ap):

Aa = Ac - (2*Ad)
i.e., Aa = 0.694- (2*0.0866)
Aa = 0.5208 m2
Now,
Lw / Dc = 0.724/ 0.94 = 0.7702

,c = 100.746 0
.  0 - ,c
0
0
LH .
 - 100.746
.  0
Now,
Acz = 2* Lw* (thickness of distribution)
Where, Acz = area of calming zone, m2
Acz = 2*0.724* (30*10-3) = 0.04344 m2 -------- (which is 6.26% of Ac)
Also,
Awz

,c /360 0) -  'c -0.03) 2 ,F  0)}

2
 ^ 'c 

Where, Awz = area of waste periphery, m2

i.e., Awz

2
 ^  *

(100.746 0/360 0)

-  -0.03) 2 * (100.746 0/360 0)}

i.e., Awz = 0.0244 m2 --------- (which is 3.515% of Ac)

Now,
Ap = Ac - (2*Ad) - Acz - Awz
i.e., Ap = 0.694- (2*0.0866) - 0.04344 - 0.0244
Thus, Ap = 0.453 m2
8 Total Hole Area (Ah):
Since,
Ah / Ap = 0.1
Ah = 0.1* Ap
i.e., Ah = 0.1*0.453
Ah = 0.0453 m2
Thus, Total Hole Area = 0.04147 m2
Now we know that,
Ah = nh  Gh2
Where nh = number of holes.
nh = (4*Ah   Gh2 )
i.e., nh

   

nh = 2307.21 
Therefore, Number of holes = 2308.

Weir Height (hw):

Let, hw = 50 mm.

10

Weeping Check
All the pressure drops calculated in this section are represented as mm head of
liquid on the plate. This serves as a common basis for evaluating the pressure
drops.
Notations used and their units:
hd = Pressure drop through the dry plate, mm of liquid on the plate
uh = Vapor velocity based on the hole area, m/s

how = Height of liquid over weir, mm of liquid on the plate

h = Pressure drop due to bubble formation, mm of liquid
hds= Dynamic seal of liquid, mm of liquid
hl = Pressure drop due to foaming, mm of liquid
hf = Pressure drop due to foaming, actual, mm of liquid
Df = Average flow length of the liquid, m
Rh = Hydraulic radius of liquid flow, m
Uf = Velocity of foam, m/s
(NRe) = Reynolds number of flow
f = Friction factor
hhg = Hydraulic gradient, mm of liquid
hda = Loss under down comer apron, mm of liquid
Ada = Area under the down comer apron, m2
C = Down comer clearance, m
hdc = Down comer backup, mm of liquid
Calculations:
Head loss through dry hole
hd = head loss across the dry hole
hd = k1 + [k2* (g/l) *Uh2] --------- (eqn. 18.6, page 18.9, 6th edition Perry)
where Uh =gas velocity through hole area
k1, k2 are constants
For sieve plates
k1 = 0

and

k2 = 50.8 / (Cv)2
where Cv = discharge coefficient, taken from fig. edition 18.14, page 18.9 6th Perry).
Now,
(Ah/Aa) = 0.0453/ 0.5208 = 0.087
also tT/dh = 3/5 = 0.60

Thus for (Ah/Aa) = 0.087 and tT/dh = 0.60

We have from fig. edition 18.14, page 18.9 6th Perry.
Cv = 0.73

k2 = 50.8 / 0.732 = 95.327

Volumetric flow rate of Vapor at the top of the Stripping Section

=qt =1.8974/ (3.425) = 0.554m3/s -------- (minimum at top)
Volumetric flow rate of Vapor at the bottom of the Stripping Section
= qo = 1.9657 / (3.361) = 0.5848 m3/s. ------- (maximum at bottom).
Velocity through the hole area (Uh):
Now,
Velocity through the hole area at the top = Uh, top = qt /Ah
= 0.554/0.0453= 12.23 m/s
also, Velocity through the hole area at the bottom= Uh, bottom = qo /Ah
= 0.5848/0.0453 = 12.91 m/s
Now,
hd, top = k2 [g/l] (Uh,top)2
= 95.327(3.425/784.50) 12.232
hd, top = 62.25 mm clear liquid. -------- (minimum at top)
also
hd, bottom = k2 [g/l] (Uh, bottom)2
= 95.327(3.361/747.87)12.912
hd, bottom = 71.4 mm clear liquid ----- (maximum at bottom)
Head Loss Due to Bubble Formation
h = 409 [ / ( Ldh) ]
Where =surface tension, mN/m (dyne/cm)
dh = Hole diameter, mm
l = average density of liquid in the section, kg/m3
l = 784.5 kg/m3

h = 409 [18.33 / ( 784.5 * 5)]

h = 1.911 mm clear liquid.
Height of Liquid Crest over Weir
how = 664Fw [(q/Lw)2/3]
q = liquid flow rate at top, m3/s
= 0.0035 m3/s.
q = 1.998 * 60 / 7193.9 = 0.0166 m3/min = 4.384 gal/min.
Thus, q = 4.384 gal/min.
Lw = weir length = 0.724 m = 2.3753 ft
Now,
q/Lw2.5 = 4.384/ (2.375)2.5 = 0.504
Now for q/Lw2.5 = 0.504 and Lw /Dc =0.7702
We have from fig.18.16, page 18.11, 6th edition Perry
Fw = correction factor =1.02
Thus, how = 1.02664 [(0.00035)/0.724] 2/3
how = 4.17 mm clear liquid.
Now,
(hd + h) = 62.25 + 1.911 = 64.161 mm ------ Design value
(hw + how) = 50 + 4.17 = 54.17mm
Also, Ah/Aa = 0.087 and (hw + how) =50 +4.17 = 54.17 mm
The minimum value of (hd + h ) required is calculated from a graph given in Perry,
plotted against Ah/Aa.
i.e., we have from fig. 18.11, page 18.7, 6th edition Perry
(hd + h)min = 12.0 mm ------- Theoretical value.
The minimum value as found is 12.0 mm.
Since the design value is greater than the minimum value, there is no problem of
weeping.
Down comer Flooding:
hds =hw + how + (hhg /2) ------- (eqn 18.10, page 18.10, 6th edition Perry)
Where,

hw = weir height, mm
hds = static slot seal (weir height minus height of top of slot above plate
floor, height

equivalent clear liquid, mm)

how = height of crest over weir, equivalent clear liquid, mm

hhg = hydraulic gradient across the plate, height of equivalent clear liquid, mm.
Hydraulic gradient, hhg
Let hhg = 0.5 mm.
hds = hw + how + hhg/2
= 50 + 4.17 + 0.5/2 = 54.42 mm.
Now, Fga = Ua g0.5
Where Fga = gas-phase kinetic energy factor,
Ua = superficial gas velocity, m/s (ft/s),
g = gas density, kg/m3 (lb/ft3)
Here Ua is calculated at the bottom of the section.
Thus, Ua = (Gb/g)/ Aa = 1.9657 / (0.52083.361) = 1.123 m/s
Thus, Ua = 3.684 ft/s
g = 3.361 kg/m3 = 0.205 lb/ft3
Therefore, Fga = 3.684 (0.205) 0.5
Fga = 1.668
Now for Fga = 1.668, we have from fig. 18.15, page 18.10 6th edition Perry)
Aeration factor = = 0.61
Relative Froth Density = t = 0.21
Now hl= hds ---- (eqn. 18.8, page 18.10, 6th edition Perry)
Where, hl= pressure drop through the aerated mass over and around the disperser,
mm liquid,
hl= 0.61 54.42 = 33.1962 mm.
Now,
hf = hl/t ------- (eqn. 18.9, page 18.10, 6th edition Perry)
hf = 33.1962/ 0.21 = 158.07 mm.

Average width of liquid flow path, Df = (Dc + Lw)/2

= (0.94 + 0.724)/2 = 0.832 m.
Hydraulic radius of aerated mass Rh = hf * Df /(2*hf + 1000*Df) (from eq. 18.23, page
18.12 6th edition Perry)
Rh = 158.07*0.832/(2*158.07 + 1000*0.832)
= 0.1145 m.
Velocity of aerated mass, Uf = 1000*q/ (hl * Df )
Volumetric flow rate, q = 2.068/747.87 =0.00276 m3/s.
Uf = 1000* 0.00276 / (33.1962* 0.832)
= 0.0999 m/s.
Reynolds modulus NRe = Rh * Uf * l / liq
= 0.1145 * 0.0999 * 747.87 /(0.924 * 10-3)
= 9257.17
hhg = 1000* f* Uf2 *Lf/(g * Rh)
f = 0.18 for hw =1.97 and NRe = 9257.17
Lf = 2 * Dc FRV ,c / 2) = 0.5995 m.
hhg = 1000* 0.18 *0.09992*0.5995/(9.81* 0.1145)
= 0.958 mm.

Head loss over down comer apron:

hda = 165.2 {q/ Ada}2 ----- (eqn. 18.19, page 18.10, 6th edition Perry)
Where, hda = head loss under the down comer apron, as millimeters of liquid,
q = liquid flow rate calculated at the bottom of section, m3/s
And Ada = minimum area of flow under the down comer apron, m2
Now,
q = 0.00276 m3/s
Take clearance, C = 1 = 25.4 mm
hap = hds - C = 54.42 - 25.4 = 29.02 mm
Ada = Lw x hap = 0.724 29.0310-3 = 0.021 m2
hda = 165.2[(0.00276)/ (0.021)] 2
hda = 2.85 mm

Now
ht = hd + hl`
Here hd and hl are calculated at bottom of the Stripping section.
Now we have,
hd, bottom = 71.4 mm
hl, bottom = 33.1962 mm
ht = hd + hl`
= 71.4+33.1962
ht = 104.6 mm
Down comer Backup:
hdc = ht + hw + how + hda +hhg ---- (eqn 18.3, page 18.7, 6th edition Perry)
ht = total pressure drop across the plate (mm liquid)
= hd + hl`
hdc = height in down comer, mm liquid,
hw = height of weir at the plate outlet, mm liquid,
how =height of crest over the weir, mm liquid,
hda = head loss due to liquid flow under the down comer apron, mm liquid,
hhg = liquid gradient across the plate, mm liquid.
hdc = 104.6 +50 +4.17 + 0.958 + 2.85
hdc = 162.58 mm.
Let dc = average relative froth density (ratio of froth density to liquid density) =
0.5
h`dc = hdc / dc = 162.58/ 0.5
h`dc = 325.16 mm.
Which is less than the tray spacing, ts= 457 mm.
Hence no flooding in the Stripping section and hence the design calculations are
acceptable.

Formulas used in calculation of properties:

1

VISCOSITY:

(i). Average Liquid Viscosity:

( liq)1/3 = [x1 ( 1)1/3] + [x2 ( 2)1/3]
2

DIFFUSIVITIES:
(i). Liquid Phase Diffusivity:
For the case of Organic solutes diffusing in Organic solvents
DAB = (1.173*10-13*(, 00.5 7  >B (VA)0.6] (Richardson coulson vol.6)
Where,

FRQVWDQW

M = molecular weight.
T = absolute temperature, 0K,

B = viscosity of solvent B, cP,

VA =molar volume of solute A at its normal boiling temperature, cm3/g-mol.
DAB =mutual diffusivity coefficient of solute A at very low concentration in
solvent B, cm2/s
(ii). Gas Phase Diffusivity:
DAB = 1.013*10-7T1.75 [(MA+MB)/ (MAMB)]1/2}/{P[(YA)1/3+ (YB)1/3]2
------ (Richardson coulson vol.6 ).

Where P = Pressure in atmospheres,

T = Temperature in 0K
DAB = Diffusivity, cm2/s
YA

and

YB

= summation of atomic diffusion volumes for

components A and B respectively.

MA and MB = Molecular weights of components A and B respectively.
3. SURFACE TENSION:

>3ch

(!l - !g)/M]4 10-12 ----- (eqn. 8.23, page 293, Coulson and Richardson

vol.6)
Where,

VXUIDFH WHQVLRQ G\QHFP

Pch =Sugdens Parachor,

!l = liquid density, kg/m3

!g = density of saturated vapor, kg/m3
M = Molecular weight

1 !l DQG !g are evaluated at system temperature.

1mix = [i 1i) where i=1,2,3,n.
4.

LIQUID DENSITY:
= Pc/ ( R * Tc * Zc[ 1 + ( 1 Tr)2/7] )

(Coulson and Richardson vol.6)

Where,
Pc = critical pressure = M/(0.34 + ( x3 2 )
M = Molecular weight.
Tc = Critical temperature = Tb / ( 0.567 + x 7 ( x 72 )
Tb = Normal boiling temperature 0K.
Zc = Pc * Vc / (R * Tc)
Vc = critical volume
R = universal gas constant.
5. GAS DENSITY:
= P * M /( R * T )
P = pressure
M = Molecular weight.
R = universal gas constant.
T = temperature.

Enriching section:
Column efficiency ( AIChe method )
1.

Point Efficiency, (Eog):

Eog = 1-e-Nog = 1-exp (-Nog) ----- (eqn. 18.33, page 18.15, 6th edition Perry)
Where Nog = Overall transfer units
Nog = 1/ [(1/Ng  1l)] ---- (eqn. 18.34, page 18.15, 6th edition Perry)
Where Nl = Liquid phase transfer units,
Ng = Gas phase transfer units,

P *m)/

Lm = Stripping factor,

m = slope of Equilibrium Curve,

Gm = Gas flow rate, mol/s
Lm = Liquid flow rate, mol/s
Ng= (0.776 + (0.0045*hw) - (0.238*Ua !g0.5) + (105*W))/ (NSc, g)0.5
----- (eqn. 18., page 18., 6th edition Perry)--- *
Where,
hw = weir height = 50.00 mm
Ua = Gas velocity through active area, m/s
= 1.232 m/s.
Ua = 1.232 m/s
Df = (Lw + Dc)/2 = (0.87 + 0.67)/2 = 0.77 m
q = 161.30 * 10 -6 m3/s
W = Liquid flow rate, m3/ (s.m) of width of flow path on the plate,
= q/Df = 161.3*10-6/0.77 = 209.48*10-6 m3/ (s.m)
NSc, g = Schmidt number = g !g*Dg) = 0.6256
Dg = Diffusivity = 4.433 * 10-6 m2/s.
Now,
Number of gas phase transfer units,
Ng=(0.776+(0.0045*50)-(0.238*1.232*3.4250.5)+(105*209.48*10-6))/ (0.6256)0.5
10 Ng = 0.6073
Also,

Number of liquid phase transfer units,

Nl = kl* a*l ----- (eqn 18.36a, page 18.15, 6th edition Perry)
Where kl = Liquid phase transfer coefficient kmol/ (sm2 kmol/m3) or m/s
a = effective interfacial area for mass transfer m2/m3 froth or spray on the plate,
l = residence time of liquid in the froth or spray, s
l = (hl*Aa)/ (1000*q) ---- (eqn. 18.38, page 18.16, 6th edition Perry)
Now,

q = liquid flow rate, m3/s

q = 161.30*10 -6 m3/s
hl = hl = 31.662 mm
Aa = 0.4496 m2
l = 31.662*0.4496/ (1000*161.3*10-6) = 88.25 s
kl *a = (3.875*108*DL)0.5* ((0.40*Ua !g0.5) + 0.17)
--- (eqn. 18.40a, page 18.16, 6th edition Perry)
DL= liquid phase diffusion coefficient, m2/s
kl *a = (3.875*108*2.002*10-9)0.5* ((0.40*1.232*3.4250.5) + 0.17)
kl *a = 0.933 m/s
Nl = kl* a*ql

i.e., Nl = 0..933*88.25
m = mm * Gm/Lm
b = 0.5990
t = 0.3
= 0.4495
Nog = 1/ [(1/Ng  1l)]
= 1/ [(1/1.093) + (0.4495/82.33)]
Nog = 1.0865
Eog = 1-e-Nog = 1-exp (-Nog)
= 1-e-1.0865 = 1-exp (-1.0865)
Eog = 0.6626
Point Efficiency = Eog = 0.6626

Murphree Plate Efficiency (Emv):

Now, Peclet number =NPe = Zl2 / (DE* ql)
Zl = length of liquid travel, m
DE = (6.675 * 10 3* (Ua) 1.44) + (0.922 * 10 4* hl) - 0.00562
----- (eqn. 18.45, page 18.17, 6th edition Perry)
Where DE = Eddy diffusion coefficient, m2/s
DE = (6.675 * 10 3* (1.232) 1.44) + (0.922 * 10 4* 31.662) - 0.00562
DE = 0.0063 m2/s
Also,
Zl = Dc FRV c/2) = 0.87* cos (100.73 0/2) = 0.555 m
NPe = Zl2 / (DE* l)
= 0.5552 / (0.0063 * 88.25)
NPe = 0.554

 (og = 0.4495 * 0.6626 = 0.2978

1RZ IRU

 (og = 0.2978 and NPe = 0.554

We have from fig.18.29a, page 18.18, 6th edition Perry

Emv/ Eog = 1.09
Emv = 1.09* Eog = 1.09*0.6626 = 0.722
Murphee Plate Efficiency = Emv = 0.722

3 Overall Efficiency ( EOC):

Overall Efficiency = EOC = log [1 + Ea ( - 1)]
log
----- (eqn. 18.46, page 18.17, 6th edition Perry)
Where, E /Emv= 1/ (1

EMV [/ (1- )])

----- (eqn. 18.27, page 18.13, 6th edition Perry)

Emv = Murphee Vapor efficiency,
E. = Murphee Vapor efficiency, corrected for recycle effect of liquid
entrainment.

(L/G)*{g/l}0.5 = 0.004
Thus, for (L/G)*{g/l}0.5 = 0.004 and at 80 % of the flooding value,
We have from fig.18.22, page 18.14, 6th edition Perry
= fractional entrainment, moles/mole gross down flow = 0.095
E / Emv = 1 / (1

Emv [/ (1- )]

E = Emv / ( 1 + Emv [/ (1- )] )

= 0.722/ (1+0.722[0.095/ (1-0.095)])
E = 0.6711
Overall Efficiency = EOC = log [1 + E ( - 1)]
log
EOC = log [1+ 0.6711(0.4495-1)]/ log 0.4495
Overall Efficiency = EOC = 0.5767
Actual trays = Nact = NT/EOC = (ideal trays)/ (overall efficiency)
Where NT = Theoretical plates,
Nact = actual trays
Nact = 2/0.5767 = 3.47 
Thus, Actual trays in the Enriching Section = 4
Thus 4th tray is the feed tray.
Total Height of Enriching section = 4*ts = 4*457 = 1828 mm = 1.828 m  P

B) Stripping Section:

Point Efficiency, (Eog):

Eog = 1-e-Nog = 1-exp (-Nog) ----- (eqn. 18.33, page 18.15, 6th edition Perry)
Where Nog = Overall transfer units
Nog = 1/ [(1/Ng  1l)] ---- (eqn. 18.34, page 18.15, 6th edition Perry)
Where Nl = Liquid phase transfer units,
Ng = Gas phase transfer units,

P *m)/

Lm = Stripping factor,

m = slope of Equilibrium Curve,

Gm = Gas flow rate, mol/s

Lm = Liquid flow rate, mol/s

Ng= (0.776 + (0.00457*hw) - (0.238*Ua !g0.5) + (104.6*W))/ (NSc, g)0.5
----- (eqn. 18., page 18., 6th edition Perry)--- *
where hw = weir height = 50.00 mm
Ua = Gas velocity through active area, m/s
= ( vapor flow rate in kg/hr)/ ( vapor density active area)
= 1.123 m/s.
Ua = 1.123 m/s
Df = (Lw + Dc)/2 = (0.724 + 0.94)/2 = 0.832 m
q = 0.00276 m3/s
W = Liquid flow rate, m3/ (s.m) of width of flow path on the plate,
= q/Df = 0.00276/0.832 =0.0033 m3/ (s.m)
NSc, g = Schmidt number = g !g*Dg) = 0.0095*10-3/(3.361*4.433*10-6)
= 0.6776
Now,
Number of gas phase transfer units,
Ng= (0.776 + (0.00457*50) - (0.238*1.232*3.3610.5) + (105*0.0033))/ (0.6776)0.5
Ng = 1.046
Also,
Number of liquid phase transfer units,
Nl = kl* a*l ----- (eqn 18.36a, page 18.15, 6th edition Perry)
Where, kl = Liquid phase transfer coefficient kmol/ (sm2 kmol/m3) or m/s
a = effective interfacial area for mass transfer m2/m3 froth or spray on the plate,
l = residence time of liquid in the froth or spray, s
l = (hl*Aa)/ (1000*q) ---- (eqn. 18.38, page 18.16, 6th edition Perry)
now, q = liquid flow rate, m3/s
q = 0.00276 m3/s
hl = hl = 33.1962 mm
Aa = 0.5208 m2
l = 33.1962*0.5208/ (1000*0.00276) = 6.264 s

kl *a = (3.875*108*DL)0.5* ((0.40*Ua !g0.5) + 0.17)

--- (eqn. 18.40a, page 18.16, 6th edition Perry)
DL= liquid phase diffusion coefficient, m2/s
kl *a = (3.875*108*2.002*10-9)0.5* ((0.40*1.232*3.3610.5) + 0.17)
kl *a = 0.875 m/s
Nl = kl* a*l
i.e., Nl = 0.875*6.264 =5.481 m
Slope of equilibrium Curve
mtop = 0.2
mbottom = 0.3
t = mt *Gm/Lm = 2.85
b = mb*Gm/Lm = 0.19

= 1.52

Nog = 1/ [(1/Ng  1l)]

= 1/ [(1/1.046) + (1.52/5.481)]
Nog = 0.8108
Eog = 1-e-Nog = 1-exp (-Nog)
= 1-e-0.8108 = 1-exp (-0.8108)
Eog = 0.5555
Point Efficiency = Eog = 0.5555
2

Murphee Plate Efficiency (Emv):

Now, Pelect number =NPe = Zl / (DE* ql)
Zl = length of liquid travel, m
DE = (6.675 * 10 3* (Ua)1.44) + (0.922 * 10 4* hl) - 0.00562
----- (eqn. 18.45, page 18.17, 6th edition Perry)
Where, DE = Eddy diffusion coefficient, m2/s
DE = (6.675 * 10 3* (1.123) 1.44) + (0.922 * 10 4* 33.1962) - 0.00562
DE = 0.0053 m2/s
Also,
Zl = Dc FRV c/2) = 0.94* cos (100.746 0/2) = 0.5995 m

NPe = Zl2/ (DE* l)

= 0.59952 / (0.0053 * 6.264)
NPe = 10.82

 (og = 1.52 * 0.5555 = 0.844

1RZ IRU

 (og = 0.844 and NPe = 10.82

We have from fig.18.29a, page 18.18, 6th edition Perry

Emv/ Eog = 1.49
Emv = 1.49* Eog = 1.49*0.5555 = 0.8276
Murphree Plate Efficiency = Emv = 0.8276
3

Overall Efficiency ( EOC):

Overall Efficiency = EOC = log [1 + Ea ( - 1)]
log
----- (eqn. 18.46, page 18.17, 6th edition Perry)
where E /Emv= 1/(1

Emv [/ (1- )])

----- (eqn. 18.27, page 18.13, 6th edition Perry)

Emv = Murphee Vapor efficiency,

E. = Murphee Vapor efficiency, corrected for recycle effect of liquid
entrainment.
(L/G)* {g/l}0.5 = 0.0705
thus, for (L/G)*{g/l}0.5 = 0.0705 and at 80 % of the flooding value,
we have from fig.18.22, page 18.14, 6th edition Perry
= fractional entrainment, moles/mole gross down flow = 0.04
E /Emv = 1/ ( 1 + Emv [/ (1- )] )

E = Emv/1 + Emv [/ (1- )]

= 0.8276/ (1+0.8276[0.04/ (1-0.04)])

E = 0.8
Overall Efficiency = EOC = log [1 + E ( - 1)]
log

EOC = log [1+ 0.8(1.52-1)]/ log 1.52

Overall Efficiency = EOC = 0.83
Actual trays = Nact = NT/EOC = (ideal trays)/ (overall efficiency)
Where NT = Theoretical plates,
Nact = actual trays
Nact = 2/0.83 = 2.41 
Thus, Actual trays in the Stripping Section = 3
Total Height of Stripping section = 3*ts = 3*457 = 1371 mm = 1.371 m
Total Height of Column =HC = Height of Enriching section + Height of Stripping section
= 2+ 1.371= 3.371 m  P

SUMMARY OF THE DISTILLATION COLUMN:

A) Enriching section
Tray spacing = 457 mm
Column diameter = 870 mm = 0.87 m
Weir length = 0.67 m
Weir height = 50 mm
Hole diameter = 5 mm
Hole pitch = 15 mm, triangular
Tray thickness = 3 mm
Number of holes = 1971
Flooding % = 80%
B) Stripping section
Tray spacing = 457 mm
Column diameter = 940 mm = 0.94 m
Weir length = 0.724 m
Weir height = 50 mm
Hole diameter = 5 mm
Hole pitch = 15 mm, triangular
Tray thickness = 3 mm
Number of holes = 2308, Flooding % = 80%

5.2 MECHANICAL DESIGN OF DISTILLATION COLUMN:

a) Shell:
Diameter of the tower =Di = 940 mm =0.940 m
Working/Operating Pressure = 2.087 atmosphere = 2.1558 kg/cm2
Design pressure = 1.1*Operating Pressure = 1.1*2.1558 = 2.37138 kg/cm2
Working temperature = 95 0C = 368 0K
Design temperature = 104.5 0C = 377.5 0K
Shell material - IS: 2002-1962 Carbon steel (specific gravity 7.7)
Permissible tensile stress (ft) = 95 MN/m2 = 970 kg/cm2
Insulation material - asbestos
Insulation thickness = 2= 50.8 mm
Density of insulation = 2700 kg/m3
Top disengaging space = 0.3 m
Bottom separator space = 0.4 m
Weir height = 50 mm
Down comer clearance = 1 = 25.4 mm
b) Head - torispherical dished head:
Material - IS: 2002-1962 Carbon steel
Allowable tensile stress = 95 MN/m2 = 970 kg/cm2
c) Support skirt:
Height of support = 1000 mm = 1.0 m
Material - Carbon Steel
d) Trays-sieve type:
Number of trays = 7
Hole Diameter = 5 mm
Number of holes:
Enriching section = 1971
Stripping section = 2308
Tray spacing:
Enriching section: 18 = 457 mm
Stripping section: 18 = 500 mm

Thickness = 3 mm
e) Support for tray:
Purlins - Channels and Angles
Material - Carbon Steel
Permissible Stress = 127.5 MN/m2 =1299.7 k gf/cm2

1. Shell minimum thickness:

Considering the vessel as an internal pressure vessel.
ts = ((P*Di)/ ((2*ft*J)- P)) + C
where ts = thickness of shell, mm
P = design pressure, kg/cm2
Di = diameter of shell, mm
ft = permissible/allowable tensile stress, kg/cm2
C = Corrosion allowance, mm
J = Joint factor
Considering double welded butt joint with backing strip
J= 85% = 0.85
Thus, ts = ((2.37138*940)/ ((2*970*0.85)- 2.1558)) + 3 = 4.35 mm
Taking the thickness of the shell = 6 mm (standard)

2. Head Design- Shallow dished and Torispherical head:

Thickness of head = th = (P*Rc*W)/ (2*f*J)
P =internal design pressure, kg/cm2
Rc = crown radius = diameter of shell, mm
W= stress intensification factor or stress concentration factor for torispherical
head,
W= * (3 + (Rc/Rk)0.5)
Rk = knuckle radius, which is at least 6% of crown radius, mm
Now, Rc = 940 mm
Rk = 6% of Rc = 0.06*940 = 56.4 mm
W= * (3 + (Rc/Rk)0.5) = * (3 + (940/56.4)0.5) = 1.7706 mm

th = (2.37138*940*1.7706)/ (2*970*0.85) = 2.39 mm

Including corrosion allowance take the thickness of head = 6 mm

Weight of Head:
Diameter = O.D + (O.D/24) + (2*sf) + (2*icr/3) --- (eqn. 5.12 Brownell and Young)
Where O.D. = Outer diameter of the dish, inch
icr = inside cover radius, inch
sf = straight flange length, inch
From table 5.7 and 5.8 of Brownell and Young
sf =1.5
icr = 2.31
Also, O.D.= 940 mm = 37
Diameter = 37+ (37/24) + (2*1.5)+(2/3*2.31)
d = 43.08 = 1094.23 mm.
:HLJKW RI +HDG

2
 G W !

 

*0.2362)/4) * (7700/1728) = 1534.15 lb

= 695.87 kg

3. Shell thickness at different heights

At a distance Xm from the top of the shell the stresses are:
3.1 Axial Tensile Stress due to Pressure:
fap =

P*Di_ =

2.37138*940_ = 185.758 kgf/cm2 .

4(ts -c)

4(6 - 3)

This is the same through out the column height.

3.2 Circumferential stress
2 * fap = 2*185.758 = 371.516 k gf/cm2

3.3 Compressive stress due Dead Loads:

3.3.1 Compressive stress due to Weight of shell up to a distance X meter from top.
fds = weight of shell/cross-section of shell
2
 'o -

Di2 !s ;  'o2- Di2)

ZHLJKW RI VKHOO SHU XQLW KHLJKW ;  'm

* (ts- c))

Where Do and Di are external and internal diameter of shell.

!s = density of shell material, kg/m3

Dm = mean diameter of shell,
ts = thickness of shell,
c = corrosion allowance
1RZ

!s = 7700 kg/m3 =0.0077 kg/cm3

fds !s* X = (7700*X) kg/m2 = (0.77*X) kg/cm2

The vessel contains manholes, nozzles etc., additional weight may be estimated 20% of
the weight of the shell.
fT,ds = 1.2 * 7700*X = 0.924* (X) kg/cm2

3.3.2 Compressive stress due to weight of insulation at a height X meter:

fd(ins) = *Dins* tins* ins *X = weight of insulation per unit height (X)
*Dm* (ts - c)

*Dm* (ts - c)

where Dins, tins, ins are diameter, thickness and density of insulation respectively.
Dm = (Dc+ (Dc+2ts))/2
Assuming asbestos is to be used as insulation material.

!ins = 2700 kg/m3

tins = 2 = 5.08 cm.
Dins =Dc+2ts+2tins = 94+ (2*0.6) + (2*5.08) = 105.36 cm.
Dm = (94+ (94+ (2*0.6)))/2 = 94.60 cm.
fd(ins) = *105.36* 5.08*2700*X = 50920.28 *X kg/m2
*94.6* (0.6 - 0.3)
= 5.092028*X kg/cm3

3.3.3 Stress due to the weight of the liquid and tray in the column up to a height X
meter.
fd, liq. =

weight

of liquid and tray per unit height X

*Dm* (ts - c)

The top chamber height is 0.3 m and it does not contain any liquid or tray. Tray

spacing is 457 mm.

Average liquid density = 775.45 kg/m3
Liquid and tray weight for X meter
Fliq-tray = [(X-   @  'i2/4) !l
= [(X-   @  2/4) *775.4
= [2X + 0.4] * 538.11 kg
fd (liq) = Fliq-tray *10/ (*Dm* (ts - c))
= [2X + 0.4] * 538.11 *10/ (*946* (6 - 3))
= [2X + 0.4] * 0.6035
= 1.207*X + 0.2414 kg/cm2
3.3.4 Compressive stress due to attachments such as internals, top head, platforms
and ladder up to height X meter.
fd (attch.) = weight of attachments per unit height X
*Dm* (ts - c)
Now total weight up to height X meter = weight of top head + pipes +ladder, etc.,
Taking the weight of pipes, ladder and platforms as 25 kg/m = 0.25 kg/cm
Total weight up to height X meter = (695.87+25X) kg
fd (attch.) = (695.87+25X) * 10/ *946* (6 - 3) = 0.7805 + 0.028X kg/cm2
Total compressive dead weight stress:
fdx = fds + fins +fd (liq) + fd (attch)
= 0.924X + 5.092X + [1.207X+0.2414] + [0.7805 +0.028X]
fdx = 7.251X + 1.0219 kg/cm2

4. Tensile stress due to wind load in self supporting vessels:

fwx = Mw /Z
Where, Mw = bending moment due to wind load = (wind load* distance)/2
= 0.7*Pw*D*X2/2
Z = modulus for the section for the area of shell 'm2* (ts-c)/4
Thus, fwx =1.4*Pw*X2  'm* (ts-c))
Now Pw = 25 lb/ft2

--- (from table 9.1 Brownell and Young)

= 37.204 kg/m2

Bending moment due to wind load

Mwx = 0.7*37.204*0.94*X2/2 = 12.24(X2) kg-m
fwx= 1.4*37.204*X2   -3)*10-3) = 0.58792(X2) kg/cm2
5. Stresses due to Seismic load:
fsx = Msx  'm2* (ts-c) / 4)
Where, bending moment Msx at a distance X meter is given by
Msx = [C*W*X2/3] * [(3H-X)/H2]
Where, C = seismic coefficient,
W= total weight of column, kg
H = height of column
Total weight of column = W= Cv !m*Dm*g* (Hv+ (0.8*Dm))*ts*10-3
----- (eqn. 13.75, page 743, Coulson and Richardson 6th volume)
Where W = total weight of column, excluding the internal fittings like plates, N
Cv = a factor to account for the weight of nozzles, man ways, internal
supports, etc.
= 1.5 for distillation column with several man ways, and with plate
support rings or equivalent fittings
Hv = height or length between tangent lines (length of cylindrical section)
g = gravitational acceleration = 9.81 m/s2
t = wall thickness

!m = density of vessel material, kg/m3

Dm = mean diameter of vessel = Di + (t *10-3) = 0.94+ (6 *10-3) = 0.946 m
:

      )*6*10

Weight of plates:
3ODWH DUHD

-3

=7590.341 N=773.73 kg.

------- (Coulson and Richardson 6th volume)

2/4 = 0.694 m2

Weight of each plate = 1.2*0.694 = 0.8328 kN

Weight of 7 plates = 7*0.8328 = 5.8296 kN = 594.25 kg.
Total weight of column = 773.73 + 594.25 = 1367.98 kg.
Let, C = seismic coefficient = 0.08
Msx = [0.08*1367.98*X2/3] * [((3*3.4)-X)/3.42]

= 36.48X2 * [0.8823-0.086X] kg-m

fsx = Msx*103 'm2* (ts-c)/4
=36.48X2 * [0.8823-0.086X * 103  2* (6-3)/4)
= [1.526X2- 0.14878X3], kg/cm2
On the up wind side:
ft,max = (fwx or fsx) + fap -fdx
Since the chances of, stresses due to wind load and seismic load, to occur together is rare
hence it is assumed that the stresses due to wind load and earthquake load will not occur
simultaneously and hence the maximum value of either is therefore accepted and
considered for evaluation of combined stresses.
Thus,
ft,max = 0.58792X2 + 168.871- [7.215X + 1.0129]
i.e., 0.58792X2- 7.251X + 168.871 - 1.0129- 824.5 = 0
0.58792X2- 7.251X - 656.64 =0
=>

X = 40.15 m

On the down side:

fc,max = (fwx or fsx) - fap +fdx
3.075X2 - 86.1618+ [7.3580X + 0.6701] = fc,max
The column height is 3.4 m, for which the maximum value is
fc,max = 0.58792(3.4)2 - 168.871+ [7.251(3.4) + 1.0129]
= -136.408 kg/cm2
this shows that the stress on the down wind side is tensile.
ft,max = 85% of allowable tensile stress.
ft,max = 970 * 0.85 = 824.5 kg/cm2.
ft,max = 0.58792(X)2 168.871 + [7.251(X) + 1.0129] = 824.5
Therefore, X = 35.38 m.
Hence we see that the design value of the column height is more than 3.4 m, which is the
actual column height. So we conclude that the design is safe and thus the design
calculations are acceptable.
Hence a thickness of 6 mm is taken throughout the length of shell.
Height of the head = Dc/4 = 0.94/4 = 0.235 m

Skirt support Height = 1.0 m

Total actual height = 3.4 + 1 + 0.235 = 4.635 m

5.2.1 Design of Support:

a) Skirt Support:
The cylindrical shell of the skirt is designed for the combination of stresses due
to vessel dead weight, wind load and seismic load. The thickness of skirt is uniform and
is designed to withstand maximum values of tensile or compressive stresses.
Data available:
(i)

Diameter = 940 mm.

(ii)

Height = 3400 mm = 3.40 m

(iii)

Weight of vessel, attachment = 2148.85 kg.

(iv)

Diameter of skirt (straight) = 940 mm

(v)

Height of skirt = 1.0 m

(vi)

Wind pressure = 37.204 kg/m2

1. Stresses due to dead Weight:

fd = :  'ok* tsk)
fd = stress,
:

GHDG ZHLJKW RI YHVVHO FRQWHQWV DQG DWWDFKPHQWV

Dok = outside diameter of skirt,

tsk = thickness of skirt,
fd

    Wsk)

= 7.1848 / tsk kg/cm2

2. Stress due to wind load:

pw = k * p1* h1* Do
p1 = wind pressure for the lower part of vessel,
k = coefficient depending on the shape factor
= 0.7 for cylindrical vessel.
Do = outside diameter of vessel,
The bending moment due to wind at the base of the vessel is given by
Mw = pw * H/2

fwb = Mw/Z = 4 * Mw 'ok)2 * tsk )

Z- Modulus of section of skirt cross-section
pw = 0.7* 37.204*1.0*0.9 = 120.785 kg
Mw = pw *H/2 = 120.79510/2 = 603.975 kg-m
Substituting the values we get,
fwb = 708.4737/tsk kg/cm2

3. Stress due to seismic load:

Load = C*W
C = seismic coefficient,
W= total weight of column.
Stress at base, fsb

 & + :  5ok)

* tsk)

C=0.08
fsb

   

(95.2/2)2 * tsk = 0.5474/ tsk kg/cm2

Maximum tensile stress:

ft, max = (8.9458/ tsk) - (7.1848/ tsk) = (1.761/ tsk) kg/cm2
Permissible tensile stress = 925 kg/cm2
Thus, 925 = (1.761/ tsk)
=>

tsk = 1.761/925 = 0.0019 cm = 0.019 mm

Maximum compressive stress:

fc, max = (8.9458/ tsk) + (7.1848/ tsk) = (16.1306/ tsk) kg/cm2
Now,
fc, (permissible) <= ( \LHOG SRLQW
= 1500/ 3 = 500 kg/cm2
Thus, tsk = 16.1306/500 = 0.03 cm = 0.3 mm
As per IS 2825-1969, minimum corroded skirt thickness = 7 mm
Thus use a thickness of 7 mm for the skirt.

Design of skirt bearing plate:

Assume both circle diameter = skirt diameter + 32.5 = 94+ 32.5 = 126.5 cm
Compressive stress between Bearing plate and concrete foundation:
fc = (:$  0w/Z)
:

GHDG ZHLJKW RI YHVVHO FRQWHQWV DQG DWWDFKPHQWV

A = area of contact between the bearing plate and foundation,

Z = Section Modulus of area,
Mw = the bending moment due to wind,
fc =   2- 942))+(0.7*37.204*3*42.32  4 -944)/(32*126.5))
= 0.0954 + 0.506
fc = 0.6014 kg/cm2
Which is less than the permissible value for concrete.
Maximum bending moment in bearing plate
Mmax = (0.6014*16.252/2) = 79.4 kg-cm
Stress, f = (6*0.6014* 16.252)/ (2 *tB2) = 476.42/ tB2
Permissible stress in bending is 1000 kg/cm2
Thus, tB2 = 476.42/1000 => tB = 0.6902 cm = 6.902 mm
Therefore, a bolted chair has to be used.

Anchor Bolts:
Minimum weight of Vessel = Wmin = 1400 kg. ------ (assumed value)
fc,min = ( Wmin/A) - (Mw/Z)
= [(4*1400)/ 2-942))]-(0.7*37.204*3*42.32  4-944)/(32*126.5))
= 0.2487 0.5059 = - 0.2572 kg/cm2
Since fc is negative, the vessel skirt must be anchored to the concrete foundation by

anchor bolts.
Assuming there are 24 bolts,
Pbolts = (0.25  2 - 942))/4) = 19.199 kg
Trays:
The trays are standard sieve plates throughout the column. The plates have 1971
holes in Enriching section and 2308 holes in the Stripping section of 5mm diameter
arranged on a 15mm triangular pitch. The trays are supported on purloins.

5.2.2 Nozzle Design:

Nozzles are required for compensation where a hole is made in the shell. The following
nozzles are required:

1. Feed Nozzle:
Liquid Velocity = VL= 2 m/s
Area of Nozzle = (Mass of liquid in)/ !L * VL)
Mass of liquid in = 6741.976 kg/hr.
= 1.87277 kg/s
Thus, Area of Nozzle = (1.87277)/ (784.50 * 2) = 1.1936 10-3 m2
1RZ $UHD RI 1R]]OH

GN2/4 = 1.1936 *10-3 m2

dN2 = (4*1.1936 *10-3

dN = 0.03898 m = 38.98 mm.

2. Nozzle for distillate:

Gas Velocity = VG= 25.0 m/s
$UHD RI 1R]]OH

0DVV RI OLTXLG LQ !G *

VG)

Mass of vapor in = 6372.56 kg/hr.

= 1.77 kg/s
Thus, Area of Nozzle = (1.77)/ (3.4376 * 25) = 0.0206 m2
1RZ $UHD RI 1R]]OH

dN2

GN2/4 = 0.0206 m2

 

dN = 0.1619 m = 16.19 cm.

3. Nozzle for residue:
Liquid Velocity = VL= 1.0 m/s
$UHD RI 1R]]OH

0DVV RI OLTXLG LQ !L *

VL)

Mass of liquid in = 369.416 kg/hr.

= 0.1026 kg/s
Thus, Area of Nozzle = (0.1026)/ (784.87 * 1) = 1.3072 *10-4 m2
1RZ $UHD RI 1R]]OH

GN2/4 = 1.3072 m2

dN2 = (4*1.3072*10-4
dN = 0.0129 m = 12.9 mm.

5.3 Process Design of Heat exchanger

Heat exchanger used is shell and tube.
The ethanol entering from vaporizer must be heated from 1000C to 2000C using ethanol,
acetaldehyde and hydrogen mixture available at 3100C.

Shell side:
Feed (mh)=2.008 kg/sec
Inlet temperature (T1)= 1000C
Outlet temperature (T2)= 2000C
Tube side:
Inlet temperature (t1)= 3100C
Outlet temperature (t2)= 232.6900C

1) Heat balance
Qh=mh Cp (T2-T1)
= 2.008*1.97*(200-100)
= 395.576 KW

2) LMTD

LMTD=120.990C
FT=LMTD correction factor.
R=0.7731 & S=0.476
From graph of FT Vs S
FT =0.91
LMTD (corrected)= 110.10090 C.

3) Heat transfer area:

Choose overall heat transfer coefficient= 120 W/(m2K)
Q = UA(LMTD)
A=395576 / (120*120.99*0.91)
A=29.94m2

4) Tube selection:
in OD ,10 BWG Tubes
OD=3/4 in=19.05 mm
ID=0.685 in=17.399 mm
Length of tube =L=16ft=4.88m
Heat transfer area per tube =0.292 m2
Number of tubes= 29.94/0.292=102.53
TEMA P or S, Floating head type:
Nearest tube count from tube count table
NT= 102
2 tube passes and 1 shell pass
in tubes arranged in triangular pitch
Shell ID (Df)=305mm=12in
Corrected heat transfer area=0.292*102=29.784 m2
Corrected over all heat transfer coefficient (U)=120.63 W/(m2K)

5) Average properties of fluids

a) Shell side (ammoniated brine) at 1500C

=3.98 kg/m3
=1300*10-8 mNs/m2
Cp=1.97KJ/kg.K
k=0.0256 w/m.k
b) Tube side (water) at 250C
=2.965 kg/m3
=4.7577*10-5 mNs/m2
Cp=1.7117 KJ/kg.K
k=0.081w/m.k

6) Tube side velocity

Number of passes NP=2
Flow area =(*ID2/4)*NT/NP
=(3.14*0.0173992/4)*102/2
Aa=0.012 m2
Vt=mc/ (Aa )
=2.008/(0.012*2.965)
=56.43 m/s.
Velocity is with the range (for vapor
7) Shell side velocity
Sm=[(Pl-Do)Ls]Ds/ Pl

cross flow area at center of shell.

N Number of baffles.
LTube length.
(P -D )*L Flow area between two adjacent tube rows.
D /P Number of tube rows.
Sm

Sm =[(25.4-19.05)*244] 305/25.4
=0.018605 m2.
Vs =mh/( Sm)
=2.008/(3.8*0.018605)
=28.4 m/s

P1 =25.4 mm.
LS = 0.8 * DS
= 0.244 m.

Nb+1=L/LS
=4.88/0.244
Nb=19 baffles
8) Shell side heat transfer coefficient:
NNU=jH Nre(NPr)1/3
NRe=VsDo/

NNu=nusselt number
NRe=Reynolds number

=28.40*19.05*10-3*3.98/(1300 * 10-8 ).
=165635
jh= 3*10-3
NPr=Cp/k
=1300*10-8 *1.97/(2.855 * 10-4) = 0.09
NNU=3*10-3 *165635 * 0.090.33 =222.68
ho=222.68* 0.0256 / 0.01905 = 299.244 W / m2 K.

9) Tube side heat transfer coefficient:

NNu=0.023(NRe)0.8 (NPr)0.3
NRe=61187.4
NPr=0.796
NNu=0.023(61187.4)0.8 (0.796)0.3 =2 26.82
hi=1055.9 w/m2.K

10) Overall heat transfer coefficient:

Dirt coefficient =3.522*10-4 w/m2.K
1/U=1/ho+(Do/Di)(1/hi)+Doln(Do/Di)/(2*KW)+dirt coefficient
1/U=1/299.24+(19.05/17.399)(1/1055.9)+0.01905*ln(19.05/17.399)/(2*50)+
+3.522*10-4
U=210.608 w/m2.K
Designed value is greater than the assumed value.
11) Pressure drop calculation:
11a) Tube side pressure drop:

Tube side Reynolds number=NRe= 61187.4

Friction factor=f=0.079(NRe)-1/4 = 0.079(61187.4) -1/4 = 5.023*10-3
PL= (4fLvt2/2gDi)*tg
= (4*5.023*10-3*4.88*56.432/2*9.8*17.399*10-3)*2.965*9.8
= 20603.08 N/m2
PE= 2.5(t vt2/2)
= 2.5(2.965*56.432/2)
= 11802 N/m2
(P)T = Np(PL+PE)
= 2*(20603.08 +11802) = 64810 N/m2 = 64.810 kPa.

11b) Shell side pressure drop (Bells method):

Shell side Reynolds number=NRe=165635
fk=0.1
Pressure drop for cross flow zones
PC = (bfkw2NC/fSm2)(w/f)
Nc= number of tube rows crossed in one cross flow section.
Nc=Ds[1-2(LC/Ds)]/PP
Where, Lc baffle cut,25% of Ds
PP=((3)/2)PI
Nc=0.305*[1-2*0.5]/0.022
Nc= 7
PC = (2*10-3*0.1*2.0082*7)/(3.98*0.0186052)
PC = 0.076 K Pa
Pressure drop in end zones:
PE= PC(1+Ncw/Nc)
Ncw=0.8LC/PP, number of cross flow rows in each window.
Ncw= 3
PE= 0.076*(1+3/7)
PE= 0.10857 kPa.

Pressure drop in window zones:

Pw= bw2(2+0.6Ncw)/(SmSw )
Sw=Swg- Swt
Sw=area for flow through window zone.
Swg= gross window area
Swt= area occupied by tubes
Swg= 25 in2 =0.01613 m2, for DS=12in & LC/DS=0.25
Swt= (NT/8)(1- FC) DO2
FC =0.63 for LC/DS=0.25
Swt= (102/8)(1-0.63) * *0.019052
Swt= 5.378*10-3 m2.
SW = (0.01613-5.378*10-3) = 0.010752 m2
PW = 5*10-5*2.0082*(2+0.6*8)
0.018605*0.010752 *3.98
PW = 0.962 kPa
(PS)T = 2PE + (Nb-1)PC + Nb Pw
(PS)T = 2*1.69 + (8-1)*1.19 + 8*1.127
(PS)T = 21.8174 kPa

5.4 Mechanical design of Heat Exchanger:

(a) Shell side details:
Material: carbon steel
Number of shell passes: one
Working pressure: 0.3N/mm2
Design pressure: 0.33N/mm2
Inlet temperature: 1000C
Out let temperature:2000C
Permissible stress for carbon steel: 95N/mm2
(b) Tube side details:
Number tubes: 102
Number of passes: 2

Outside diameter: 19.05mm

Inside diameter: 17.399 mm.
Length: 4.88m
Pitch triangular:1 inch
Working pressure: 0.3 N/mm2
Design pressure: 0.33N/mm2
Inlet temperature: 3100C
Outlet temperature: 232.690C

Shell side:
(1) Shell thickness:
ts= PD/(2fJ+P)
= 0.33*305/(2*95*0.85+0.33) = 0.57
Minimum thickness of shell must be=6.0 mm
Including corrosion allowance shell thickness is 8mm
(2) Head thickness:
Shallow dished and torispherical
ts = PRcW/2fJ
= 0.33*305*1.77/(2*95*0.85)
= 1.103 mm.
Minimum shell thickness should be 10mm including corrosion allowance.
(3) Transverse Baffles:
Baffle spacing =0.8*Dc = 244mm
Number of baffles,
Nb+1=L/LS=4.88/0.244=20
Nb=19
Thickness of baffles, tb=6mm
(4) Tie Rods and spacers
For shell diameter, 300-500mm
Diameter of Rod = 9mm
Number of rods=4

(5) Flanges
Design pressure=0.33 N/mm2
Flange material IS: 2004-1962,class 2
Bolting steel: 5% Cr-Mo steel
Gasket material: asbestos composition
Shell thickness: 8mm=go
Outside diameter of shell: 305 mm
Allowable stress of flange material: 100MN/m2
Allowable stress of bolting material = 138 MN/m2
Shell thickness = 10 mm.
Outside diameter = 325 mm.

Determination of gasket width:

dO/di = [(y-Pm)/(y-P(m+1))]0.5
Assume a gasket thickness of 10 mm
y = minimum design yield seating stress = 25.5 MN/m2
m = gasket factor = 2.75
dO/di = [(25.5-0.33*2.75)/(25.5-0.33(2.75+1))]0.5
dO/di = 1.0067
Let, di of gasket equal 335mm
do= 1.002*di
do= 0.33724 m
Minimum gasket width = (337.24-335)/2 = 1.12mm =0.00112 m.
Taking gasket width of N= 0.010m
do=0.35924 m.
Basic gasket seating width, bo=5mm
Diameter of location of gasket load reaction is
G= di + N = 0.335 + 0.01 = 0.345 m

Estimation of Bolt loads:

Load due to design pressure

H = G2P/4 = 3.14*0.3452*0.33/4 = 0.03085MN
Load to keep joint tight under operation
b = 2.5 (b0)0.5 = 5.59 mm.
Hp=*G*(2b)*m*p = 3.14*0.345*(2*0.00559)*2.75*0.33 = 0.011 MN
Total operating load,
Wo= H + Hp = 0.03085 + 0.011 = 0.04185 MN.
Load to seat gasket under bolting condition
Wg = *G*b*y = 3.14*0.345*5.59*10-3*25.5 = 0.1545 MN.
Wg>Wo, controlling load=0.1545 MN

Calculation of optimum bolting area:

Am = Ag = Wg/Sg = 0.1545 /138 = 1.12*10-3 m2
Calculation of optimum bolt size:
Bolt size, M18 X 2
Actual number of bolts =20
Radial clearance from bolt circle to point of connection of hub or nozzle and back
of flange = R = 0.027 m
C =ID + 2(1.415g + R) = 325 +2[11.315+0.027*103] = 401.63mm = 0.40163
Bolt circle diameter = 0.40163 m.
Calculation of flange outside diameter
Let, bolt diameter = 18 mm.
A=C+ bolt diameter +0.02 = 0.40163 +0.018+0.02 = 0.43963 m.
Check for gasket width,
AbSG / (GN) = 1.54*10-4*20*138/(3.14*0.345*10-2) = 39.21 < 2*y.
Where, SG is the Allowable stress for the gasket material.
Flange moment computation:
(a) For operating condition
Wo=W1+W2+W3

W1=*B2*P/4 = *0.3252*0.33/4 = 0.027 MN

W2 = H-W1= 0.03085-0.027 = 0.00385 MN.
W3= Wo-H = Hp= 0.011 MN.
Mo=Total flange moment
Mo=W1a1 + W2a2 + W3a3
a1=(C-B)/2=(0.40163-0.325)/2
a1=0.038315 m
a3=(C-G)/2=(0.40163-0.345)/2
a3=0.028315 m
a2=(a1 + a3)/2= (0.038315 +0.028315)/2=0.033315 m
Mo=0.027 *0.038315 +0.00385 *0.033315 +0.011 *0.028315
Mo= 1.474*10-3 MN-m

(b) For bolting condition

Mg=Wa3
W=(Am+Ab)*Sg/2
Ab=20*1.54*10-4 =3.08*10-3 m2
Am= 1.12*10-03 m2
W=(1.12*10-03 +3.08*10-3)*138/2
W= 0.2898 MN
Mg= 0.2898 *0.028315 = 8.205*10-3 MN-m
Mg>Mo ,Hence moment under operating condition Mg is controlling, Mg=M
Calculation of flange thickness
t2 = M CF Y / (B SF), SF is the allowable stress for the flange material
K =A/B = 0.43963/0.325 = 1.3527
For K = 1.3527, Y = 10
Assuming CF =1
t2 = 8.205*10-3 *1*10(0.325*100)
t= 0.0502 m=50.2 mm
Actual bolt spacing BS = *C/n = (3.14*0.40613)/(20) = 0.063m

Bolt Pitch Correction Factor

CF = [Bs / (2d+t)]0.5 = (0.063/(2*0.018+0.0502)1/2 = 0.855
CF=0.9246
Actual flange thickness = CF*t = 0.9246*0.063 = 0.04713 m
= 0.0464 m.
Standard flange thickness available is 50 mm

Channel and channel Cover

th=Gc(K*P/f) = 0.345*(0.3*0.33/95) = 0.01114m =11.14mm
th=14mm including corrosion allowance

Tube sheet thickness

tts=F*G(0.25*P/f) = 1*0.345(0.25*0.33/95) = 0.0101m=10.1 mm
tts=13 mm including corrosion allowance.
Nozzle design:
1. Tube side Nozzle:
Velocity = VG = 28 m/s
$UHD RI 1R]]OH

0DVV RI YDSRU LQ !G *

VG)

Mass of liquid in = 2.008 kg/s

Thus, Area of Nozzle = (2.008)/ (2.965 * 28) = 0.024187 m2
1RZ $UHD RI 1R]]OH

dN 2

GN2/4 = 0.024187 m2

 

dN = 0.17548 m = 17.548 cm.

2. Shell side Nozzle:
Velocity = VG = 27 m/s
$UHD RI 1R]]OH

0DVV RI YDSRU LQ !G *

VG)

Mass of liquid in = 2.008 kg/s

Thus, Area of Nozzle = (2.008)/ (3.8 * 27) = 0.01957 m2
1RZ $UHD RI 1R]]OH

GN2/4 = 0.01957 m2

dN 2

 

dN = 0.1578 m = 15.78 cm.

Saddle support
Material: low carbon steel
Total length of shell: 4.88 m
Diameter of shell: 325 mm
Knuckle radius: 18.3 mm
Total depth of head (H)= (Doro/2) = (325*18.3/2) = 54.53 mm
Density of the steel = 7600 kg/m3.
Weight of steel vessel = 3707. 21 kg.
R=D/2=162.5 mm
Distance of saddle center line from shell end = A =0.5R=81.25 mm

Longitudinal Bending Moment

M1 = QA[1-(1-A/L+(R2-H2)/(2AL))/(1+4H/(3L))]
Q = W/2(L+4H/3) = 3707.21/2*(5.88 +4*0.03085/3) = 10975.44 kg m
M1 = 18.6 kg-m
Bending moment at center of the span
M2 = QL/4[(1+2(R2-H2)/L)/(1+4H/(3L))-4A/L]
M2 =15706.74 kg-m.
Stresses in shell at the saddle
(a) At the top most fiber of the cross section
f1 = M1/(k1 R2 t)

k1=k2=1

= 18.6/(3.14*0.16252*0.01) = 0.02242 kg/mm2

Stress in the shell at mid point
f2 =M2/(k2 R2 t)
= 9.9656 kg/mm2
Axial stress in the shell due to internal pressure

fp= PD/4t
= 3.4089*940/(4*8) = 100.136 kg/cm2 = 1.00136 kg/mm2
f2 + fp = 10.96696 kg/mm2
The sum f2 and fp is well within the permissible values.