Ethyleneglycol Design 2520of 2520equipments

Download as pdf or txt
Download as pdf or txt
You are on page 1of 41

DESIGN OF EQUIPMENTS

Process Design of Distillation Column


The detailed process design of the Monoethylene glycol column is given below.
The pictorial representation of the column is given in fig 6.1. The feed to the column is a
mixture of Monoethylene glycol(MEG), Diethylene glycol(DEG) and Triethylene
glycol(TEG). The compositions of the components are shown in the figure. The distillate
is the required product consisting of mainly MEG.

I. Thermodynamics:
The primary requirement while designing a mass transfer contact equipment is the
thermodynamic equilibrium data. The data required is in the Vapor-Liquid Equilibrium
(VLE) data for the MEG(1)-DEG(2) system. The quantity of TEG is small enough to be
neglected. The X-Y curve is shown in the fig 6.2. To develop the VLE data, a model was
used.
yi pt = γi xi Pisat

Where,
yi = mole fraction of component i in vapor
pt = total system pressure
γi = activity coefficient of component i in liquid
xi = mole fraction of component i in liquid
Pisat = saturation vapor pressure of component i

The equilibrium vapor pressure was evaluated using correlations given in


literature. The correlation was based on the critical properties of the components. Since
the two components MEG and DEG are highly polar, they form a highly non-ideal
system. To accommodate this non-ideality, an activity coefficient term was used for the
liquid phase. The activity coefficient was evaluated using the UNIFAC model. Since the
evaluation of the VLE data is highly iterative, an algorithm was developed which was
solved using a computer program. The gas phase was assumed to be ideal. This is a valid
assumption since the column is at a very low pressure (100 mmHg, abs). The high boiling
points of the two components requires the column to be operated under vacuum. The
operating pressure was chosen to be 100 mmHg(abs)

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 MEG in liquid
xD = mole fraction of MEG in distillate
xW = mole fraction of MEG 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, K
TV = temperature of vapor, K
II. Preliminary calculations

F = 74.88 kmol/hr, xF = 0.89, MF = 66.84 kg/kmol


D = 66.96 kmol/hr, xD = 0.99, MD = 62.44 kg/kmol
W = 7.92 kmol/hr, xW = 0.045, MW = 104.02 kg/kmol

From the graph (Fig 6.3 and Fig 6.4)


xD/(Rm+1) = 0.5, Rm= 0.98
R = 1.5Rm = 1.5×0.98 = 1.47
xD/(R+1) = 0.4008

L = RD = 98.4312 kmol/hr
G = L+D = (R+1)D = 165.3912 kmol/hr
L = L+qF = 98.4312 kmol/hr
G = G-(1-q)F = 90.5112 kmol/hr

III. List of parameters used in calculation.

Enriching Section Stripping Section

PARAMETER TOP BOTTOM TOP BOTTOM


Liq, kmol/hr 98.4312 98.4312 98.432 98.4312
Vap, kmol/hr 165.3912 165.3912 90.5112 90.5112
X 0.99 0.83 0.83 0.045
Y 0.99 0.89 0.89 0.045
TL, K 409.33 411.36 411.36 433.83
TV, K 409.44 411.78 411.78 435.51
MLIQ, kg/kmol 62.44 69.48 69.48 104.02
MVAP, kg/kmol 62.44 66.84 66.84 104.02
Liq, kg/hr 6146.0441 6838.9998 6838.9998 10238.8134
Vap, kg/hr 10327.0265 11054.7478 6049.7686 9414.9750
ρL, kg/m3 1052.2884 1054.2612 1054.2612 1062.5007
ρV, kg/m3 0.2445 0.2603 0.2603 0.3830
µL, cP 0.7903 0.7259 0.7259 0.4300
qL, m3/s 1.6224×10-3 1.8020×10-3 1.8020×10-3 2.6768×10-3
qV, m3/s 11.7326 11.7970 6.4560 6.8284

IV. Design of Enriching Section

Tray Hydraulics
The design of a sieve plate tower is described below. The equations and
correlations are borrowed form the 7th edition of Perry’s Chemical Engineers’ Handbook.
The procedure for the evaluation of the tray parameters is iterative in nature. Several
iterations were performed to optimize the design. The final iteration is presented here.

1.Tray Spacing (TS)


Choose tray spacing = 9” = 228.6mm
2.Hole Diameter (dh)
Choose hole diameter = 5mm
3.Hole Pitch (lp)
Choose hole pitch = 15mm, ∆ pitch
4.Tray thickness (tT)
Choose tray thickness = 3mm
5.Ratio of hole area to perforated area (Ah/Ap) (fig 6.3)
= ½ (π/4 . dH2)/[(√3/4). lp2]
= 0.1008

6.Plate Diameter (Dc)


The plate diameter is calculated based on the flooding considerations
L/G{ρG/ρL}0.5 = 0.0207
From the flooding curve, for a tray spacing of 228.6mm,
Flooding parameter, CSB,F = 0.0582m/s
Unf = Csb x { σ / 20 ) 0.2 [ (ρL - ρG) / ρG]0.5
Unf = 2.8102 m/s

Actual velocity un= 0.8un,f = 0.8×4.4213 = 3.5370m/s


Net area available for gas flow (An)
Net area = Column cross sectional area – Downcomer area.
An = Ac – Ad

Choose weir length (Lw) = 0.60(Column diameter, Dc)

From the figure (fig 6.3),


Lw / Dc = 0.60
sin ( θc/2) = [( Lw/2) / (Dc/2) ] = 0.60
⇒ θc = 73.740
Ac = (π/4) Dc2 = 0.785Dc2
Ad = (π/4) Dc2 ( θc/3600) - (Lw/2) ( Dc/2). cos(θc/2)
= 0.1608 Dc2 - Dc2 x 0.1200
= 0.0408 Dc2

Since An = Ac -Ad
3.3352 = 0.785 Dc2 - 0.0408 Dc2
⇒ Dc = 2.1164 m

Ac = (π/4) Dc2 = 3.5179 m2


Ad = 0.1828 m2
Lw = 1.2698 m
7.Perforated plate area (Ap)
Aa = Ac - 2Ad
= 3.1523m2
Perforated Area
Ap = Ac -2Ad - Acz - Awz
Ap = 2.6246 m2

Total Hole Area


Ah / Ap = 0.1008
⇒ Ah = 0.2646 m2
8.Total Hole Area
Ah = 0.1008Ap = 0.1008×2.6246 = 0.2646m2
Number of holes (nh)
Number of holes Nh = 13476
9. Choose the weir height (hW) as 12mm
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 downcomer apron, mm of liquid
Ada = Area under the downcomer apron, m2
c = Downcomer clearance, m
hdc = Downcomer backup, mm of liquid

Calculations:

Head loss through dry hole


hd = head loss across the dry hole
= k1 + k2 (ρg/ρl) Uh2
where Uh =gas velocity through hole area
k1,k2 are constants

For sieve plates


k1 = 0 and k2 = 50.8 / ( Cv)2
Cv = 0.7419
⇒ k2 = 50.8 / 0.74192 = 92.77

Velocity through the hole area = Uh = 44.3409 m/s

⇒ hd = k2 [ρg/ρl] Uh2
= 42.16 mm

Height of Liquid Crest over Weir


how = 664Fw [(q/Lw)2/3]
q = liquid flow rate at top = 1.6224×10-3 m3/s
Fw= correction factor=1.01
Lw = weir length = 1.2698 m
⇒ how = 7.90 mm clear liquid
Head Loss Due to Bubble Formation

hσ = 409 [ σ / ( ρL.dh) ]
where σ =surface tension(mN/m)
dh =Hole dia
hσ = 409 [ 48.5 / ( 1052.2884 x 5)]
hσ = 3.77 mm clear liquid

( hd + hσ) = 42.16 + 3.77 = 45.93 mm , this is the design value


( hw + how) = 12 + 2.676 = 7.90 mm

Also, Ah/Aa = 0.0839


The minimum value of ( hd + hσ ) required is calculated from a graph given in
Perry ,plotted against Ah/Aa. The minimum value as found is 8 mm. Since the design
value is greater than the minimum value, there is no problem of weeping.

Downcomer Backup

hdc = ht + hw + how + hda +hhg

ht = total pressure drop across the plate( mm liquid)


= hd + hl`

Hydraulic Gradient
The hydraulic gradient was evaluated through an iterative procedure involving
flow parameters of the liquid on the tray. The iteration yielded a value of 17.5mm. The
large hydraulic gradient is a characteristic feature of vacuum operated towers where the
additional hydraulic gradient is required to push the liquid hold up over the plate.

Head loss over downcomer apron


hda = 165.2 {q/ Ada}2
Take clearance, C = 0.5”
hap = hds - C = 28.9 – 25.4/2 = mm

Ada = Lw x hap = 0.0206 m2


∴ hda = 165.2[1.802x 10-3/ 0.0206] 2
= 1.24 mm
ht = hd + hl`
= 42.16+8.60
= 50.76 mm

hdc = ht+ hw + how + hda + hhg


=50.76+12+7.90+1.24+17.5
= 89.40 mm
h`dc = hdc / φ , where φ is the froth density.

= 89.4/ 0.5 = 178.8 mm


which is less than the tray spacing of 228 mm.
Hence no flooding in the enriching section.

Column Efficiency

Point Efficiency
Average Vapor rate = 10690.8872 kg/hr
Average Vapor Density = 0.2524 kg/m3
Active Area = 3.1523 m2
Df = (Dc + Lw )/2 = 1.6931 m

Average Liquid rate = 6492.5220 kg/hr


Average Liquid Density =1053.2748 kg/m3
q = 1.7123 x 10 -3 m3/s

Tl = 410.34K and Tg = 410.61K

∴ ( µM)L = 0.7851 cp
( µM)G = 0.0143 cp

Diffusivity of the gas is calculated and is found out to be = 6.3003 x 10-5m2/s


Similarly the liquid diffusivity is calculated and found out to be
= 3.2411 x 10 -9 m2/s

Number of gas phase transfer units

NG = KG.a.θG
KG.a = 316 Dg 0.5(1030f+867f2 )/ h0.5 L
= 1184.4776 /s
θG = (1-φ)hLAa/(1000φQ)
= 3.6652×10-3 s
NG = 4.3414

Number of liquid phase transfer units


NL = KL.a.θL
KL.a = (3.875×108DL)0.5(0.4Uaρ0.5G+0.17)
= 1.0311 /s
θL = (1-ε)hfAa/(1000q) = 15.8263 s
⇒ NL = 16.3185

Slope of equilibrium Curve

(m) top = 0.4933


(m) bottom = .6059
GM/LM = 165.3912/98.3412

λt = mt GM/LM = 0.8289
λb = mb GM/LM = 1.0181 ⇒ λ = 0.9235

∴ NOG = 1 = 3.4851
1/NG + λ / NL
EOG = 1- e -NOG = 0.9694

Murphree Plate Efficiency

θL = 15.8263 sec
ZL = DC. Cos ( θc/2) = 1.6931m

DE = 6.675 x 10 -3( Ua)1.44 + 0.922 x 10 -4 x hL - 0.00562


= 0.0397 m2/s

Pecklet Number Npe = ( zL)2/ ( DE.θL) = 4.5624


λ EOG = 0.9235 x 0.9694 =0.8952

EMV/ EOG = 1.32

⇒ EMV = 1.32 x 0.9694 = 1.2796

Overall Efficiency ( EOC)

L/G{ρG/ρL}0.5 = 0.0094
at 80 % of the flooding value we have ψ = 0.40
⇒ E α / E MV = 1
1 + EMV [ψ/(1- ψ)]
⇒ Eα = 0.6905

The overall efficiency is given by the equation :


EOC = log[ 1 + Eα ( λ - 1) ]

log λ
EOC = 0.6819
Hence the actual number of trays can be calculated as :
{ Theoretical number of trays overall column efficiency }
= 11 0.6819 ≈ 16

Height of the enriching section can be calculated as


{ Tray Spacing x Actual number of trays }
= 16 x 228 = 3657.6 mm.

V. Design of Stripping Section

Tray Spacing =ts= 152.4mm=6 ”


Hole Diameter =dh= 5 mm
Hole Pitch = lp=15 mm
Tray Thickness = tt =3 mm

Ah/Ap = Total Hole Area


Perforated Area

= ½ (π/4 . dH2)
(√3/4). lp2
= 0.1008

Plate Diameter ( Dc)


Based on entrainment flooding, L/G{ρG/ρL}0.5 has to be evaluated at a pt where
the group is maximum
L/G{ρG/ρL}0.5 = 0.0207
Unf = Csb x { σ / 20 ) 0.2 [ (ρL - ρG) / ρG]0.5
Unf = 2.8102 m/s

U net vel. = Un = 0.7286 Unf = 2.0474 m /s

The odd value of the % flooding is taken so as to make the diameter of the
stripping section equal to that of the enriching section. This ensures a considerable saving
in mechanical design and fabrication without affecting the efficiency substantially.
Maximum flow rate of vapor = Vapor flow rate( max. at bottom ) / vapor density
= 6.8284 m/s

∴ An = 6.8284 / 2.0474 = 3.3352 m2


Lw / Dc = 0.60
sin ( θc/2) = [( Lw/2) / (Dc/2) ] = 0.60
⇒ θc = 73.740
Ac = (π/4) Dc2 = 0.785Dc2
Ad = (π/4) Dc2 ( θc/3600) - (Lw/2) ( Dc/2). cos(θc/2)
= 0.1608 Dc2 - Dc2 x 0.1200
= 0.0408 Dc2

Since An = Ac -Ad
3.3352 = 0.785 Dc2 - 0.0408 Dc2
⇒ Dc = 2.1164 m
Ac = (π/4) Dc2 = 3.5179 m2
Ad = 0.1828 m2
Lw = 1.2698 m

Active Area ( Aa)


Aa = Ac - 2Ad
= 3.1523m2
Perforated Area
Ap = Ac -2Ad - Acz - Awz
Ap = 2.6246 m2

Total Hole Area


Ah / Ap = 0.1008
⇒ Ah = 0.2646 m2
Number of holes Nh = 13476

Weir Height
Take a weir height of 12 mm.

Weeping check
hd = head loss across the dry hole
= k1 + k2 (ρg/ρl) Uh2
For sieve plates
k1 = 0 and k2 = 50.8 / ( Cv)2
Cv = 0.7476
⇒ k2 = 50.8 / 0.74762 = 90.89
Weeping check is done at a point where the gas velocity is the least. Here the
velocity of the vapor is minimum at the top of the enriching section.

Hence volumetric flow rate of the gas = 6.456 m3/s


∴ velocity through the hole area = Uh = 6.456 / hole are(Ah)
= 24.3991 m/s
⇒ hd = k2 [ρg/ρl] Uh2
= 12.77 mm

Height of Liquid Crest over Weir


how = 664 [(q/Lw)2/3]
q = liquid flow rate = 1.802 x 10-3 m3/s
Lw = weir length = 1.2698 m
⇒ how = 8.4691 mm clear liquid

Head Loss Due to Bubble Formation


hσ = 409 [ σ / ( ρL.dH) ]
= 409 [ 48.50 / ( 1054.2612 x 5)]
hσ = 3.76 mm clear liquid
( hd + hσ) = 16.52 mm , this is the design value
( hw + how) = 20.47 mm

Also, Ah/Aa = 0.0839


The minimum value of ( hd + hσ ) required is calculated from a graph plotted
against Ah/Aa. The minimum value as found is 8 mm. Since the design value is greater
than the minimum value we proceed using the same value for further calculations.

Downcomer Backup
hdc = ht + hw + how + hda +hhg
ht = total pressure drop across the plate( mm liquid)
= hd + hl`
The hydraulic gradient was iteratively evaluated and found to be equal to 10.5
mm. The large gradient is again prominent as described before.
Loss under Downcomer
hda = 165.2 {q/ Ada}2
Take clearance = 0.5” = 12.70 mm
hap = hds - c = 25.47-12.70 = 12.77 mm
Ada = Lw x hap = 0.0162m2
∴ hda = 4.51mm
ht = hd + hl`
=12.77 + 10.83
= 23.60 mm

hdc = ht+ hw + how + hda + hhg


= 23.60+12+8.47+4.51+10.50
= 59.08mm
h`dc = hdc / φ , where φ is the froth density.
= 59.08 / 0.5 = 118.2 mm
which is less than the tray spacing of 152.4 mm. Thus flooding check is fulfilled
satisfactorily for the enriching section.

Column Efficiency
Point Efficiency

Number of gas phase transfer units:


NG = KG.a.θG
KG.a = 316 Dg 0.5(1030f+867f2 )/ h0.5 L
= 578.3924 /s
θG = (1-φ)hLAa/(1000φQ)
= 4.8271×10-3 s
NG = 2.8180

Number of liquid phase transfer units:


NL = KL.a.θL
KL.a = (3.875×108DL)0.5(0.4Uaρ0.5G+0.17)
= 0.9852 /s
θL = (1-ε)hfAa/(1000q) = 15.2361 s
⇒ NL = 15.0161

Slope of equilibrium Curve


(m) top = 0.6509
(m) bottom = 2.1875

λt = mt GM/LM = 0.5572
λb = mb GM/LM = 2.4905 ⇒ λ = 1.5238

∴ NOG = 1 = 2.1914
1/NG + λ / NL

EOG = 1- e -NOG = 0.8882


Murphree Plate Efficiency
θL = 15.2361 sec
ZL = DC. Cos ( θc/2) = 1.6931m
DE = 6.675 x 10 -3( Ua)1.44 + 0.922 x 10 -4 x hL - 0.00562
= 0.0151 m2/s
Pecklet Number Npe = ( zL)2/ ( DE.θL) = 12.4599
λ EOG = 1.6314
EMV/ EOG = 1.60
⇒ EMV = 1.4211

Overall Efficiency ( EOC)


L/G{ρG/ρL}0.5 = 0.0193
at 72.58 % of the flooding value we have ψ = 0.35
⇒ E α / E MV = 1
1 + EMV [ψ/(1- ψ)]
⇒ Eα = 0.8051

The overall efficiency is given by the equation :


EOC = log[ 1 + Eα ( λ - 1) ]
log λ
EOC = 0.8354
Hence the actual number of trays can be calculated as :
{ Theoretical number of trays overall column efficiency }
= 8 0.8354 ≈ 10
Height of the enriching section can be calculated as
( Tray Spacing x Actual number of trays )
= 10 x 152.4 = 1524 mm.

Summary of the Distillation Column


Enriching section
Tray spacing = 228.6mm
Column diameter = 2.1164m
Weir length = 1.2698m
Weir height = 12mm
Hole diameter = 5mm
Hole pitch = 15mm, triangular
Tray thickness = 3mm
Number of holes = 13476
Flooding % = 80

Stripping section
Tray spacing = 152.4mm
Column diameter = 2.1164m
Weir length = 1.2698m
Weir height = 12mm
Hole diameter = 5mm
Hole pitch = 15mm, triangular
Tray thickness = 3mm
Number of holes = 13476
Flooding % = 72.58

VI. Mechanical Design of Distillation Column


Diameter of the tower Di = 2.1164 m
Working pressure = 660 mmHg(abs)
Design pressure pe = 0.1 N/ mm2 (Total Vacuum)
Working temperature = 162.36 oC
Design temperature = 178.6 oC
Shell material – IS: 2002-1962 Grade I Plain Carbon steel
Permissible tensile stress (ft)= 93.195 MN/m2
Elastic Modulus (E) = 1.88×105 MN/m2
Insulation thickness = 75mm
Density of insulation = 5.64 kN/m3
Tray spacing:
Enriching section: 9”
Stripping section: 6”
Top disengaging space = 0.5m
Bottom separator space = 1.0m
Weir height = 12mm
Downcomer clearance = 0.5”
Height of support = 2m

1. Shell minimum thickness


Considering a torispherical head, head dish:
hi= 0.4101m
Tangent to tangent length = 6.4350+2/3(0.4101)
= 6.7048m
D/L = 0.3155 => K = 0.232, = 2.46
P(allowable) = KE(t/D)2.46
0.1 = (0.232)(1.88×105)(t/2.1164)2.46
t = 10.79mm, take 12mm (standard)

Check for Plastic deformation


P = 2f(t/D)(1+1.5U(1-0.2D/L))/(100t/D)
U = 1.5% (for new equipment)
Substituting the values, we get P(allowable) = 0.224MN/m2
The allowable pressure is greater than the design pressure. Hence, the thickness
is satisfactory with respect to plastic deformation.

2. Head Design
A torispherical head of the following parameters is chosen:
Rc – crown radius, Rc = Di = 2.1164m
Rk - knuckle radius, Rk = 10% of Rc
= 0.21164 m
t = thickness of the head =12mm
Pressure at which elastic deformation occurs
P(elastic) = 0.366E(t/ Rc)2
= 2.2121 MN/ m2
The pressure required for elastic deformation, P(elastic)>3(Design Pressure)
Hence, the thickness is satisfactory. The thickness of the shell and the head are made
equal for ease of fabrication.

3. Shell thickness at different heights


At a distance ‘X’ m from the top of the shell the stress are;
3.1. Axial Stress: (compressive)
fap = pi Di___
4(ts – C)
= 4.4 N/m 2

3.2. Compressive stress due to weight of shell up to a distance ‘X’


fds = π/4 * ( Do2 – Di2 ) ρs X
π/4 * ( Do2 – Di2 )
= ρs X
= 0.0077 X N/mm2
3.3. Compressive stress due to weight of insulation
fd(ins) = π Dins tin ρins
π Dm (ts – C )
fd(ins) = 2290.4 * 75 * 0.0564X
2174.9 * (12)
= 0.0358 X N/mm2
3.4. Stress due to the weight of the liquid supported
fl = Wl_____
π ts (Di + ts)
Enriching section
F(liq) = ((X-top space)/TS+1) (πd2/4)ρL
= 14157.8X – 3842.4
fd(liq) = 0.1727X – 0.0469

Stripping section
F(liq) = ((X-top space)/TS+1) (πd2/4)ρL+ 14157.8X – 3842.4
= 35394.6X – 84145.2
fd(liq) = 0.4383X – 1.0250

3.5. Stress due to the weight of the attachments


The total weight of the attachments
Wa = (26700+1400X) N
fa = ___Wa_____
π ts (Di + ts)
= 0.3256 + 0.0171X N/mm2
The total dead load stress, fw acting along the axial direction of shell at
point is given by:
fw = fds + fd + fl + fa
Enriching section: fw = 0.2787+0.2333X
Stripping section: fw = -0.6994+0.4989X

3.6. wind loads


Stress due to wind,
fwx = 1.4 pw X2__
π Do (ts – C)
where pw is the wind pressure, which is ≅ 1300 N/mm2
fwx = __ 1.4 * 1300 X2_____
π * 21404 *(1.2 )
= 0.0226 X2 N/mm2
Taking joint efficiency as 0.8
(Total force, tensile) = jf(max) = jf(allowable)

Enriching section
0.0226 X2 + 4.4 –(0.2787+0.2333X) = ft max= 0.8(93.195)
X = 64m
Hence, the entire tower can be constructed keeping the same 12mm thickness
Stripping section
0.0226X2 – 4.4-(-0.6994+0.4989X) = 0.8(93.195)
X = 68.1m
Hence, the 12mm thickness is sufficient for the stripping section also.

Design of Support
Skirt Support
D = 2.1164 m
Minimum weight of the vessel with attachments:
W(min) = π(D+t)t(H)γs+2(shell weight)
= 55.9773 kN
W(max) = W(shell) + W(insulation) + W(water during test) + W(attachments)
= 41.1913+18.4986+226.1696+29.7526
= 316kN
Period of vibration at minimum dead weight is:
T(min) = 6.35×10-5(H/D)1.5(W(min)/t)0.5
= 0.0362 < 0.5s
∴K = 1.0
Period of vibration at maximum dead weight load:
T(max) = 6.35×10-5(H/D)1.5(W(max)/t)0.5
= 0.0860 < 0.5
∴K = 1.0

Skirt
Stresses due to wind load

P(wind)=kKp(wind)HD
p(wind) = 1 kN/m2
k = 0.7 for cylindrical surface.
For minimum weight condition, D = 2.1164m
P(wind), min = (0.7)(1)(1000)(8.7048)(2.1164)
= 12896 N
For maximum weight condition, D = 2.1164+2(0.075)
= 2.2664m
P(wind),max = (0.7)(1)(1000)(8.7048)(2.2664)
= 13810 N
Minimum wind moment
M(wind), min = P(wind),min × H/2
= 12896(8.7048/2)
= 56.13 kJ
Maximum wind load
M(wind),max = P(wind),max× H/2
= 13810(8.7048/2)
= 60.11 kJ
Assuming a small skirt thickness, D(in) = D(out) = 2.1164m
σzwm,min = 4(M(wind),min)/( πD2t)
= 4(56.13×10-3)/( π×2.11642×t)
= 0.0160 / t MN/m2
σzwm, max = 4(M(wind),max)/( πD2t)
= 4(60.11×10-3)/( π×2.11642×t)
= 0.0171 / t MN/m2

Minimum and maximum dead load stresses


σzw,min = W(min)/( 3.14×d×t)
= 55.9773×10-3/(3.14×2.1164×t)
= 0.0084/t MN/m2
σzw,max = W(max)/( 3.14×d×t)
= 316×10-3/(3.14×2.1164×t)
= 0.0475/t MN/m2

Maximum tensile stress without any eccentric load


σz (tensile) = σzwm, min - σzw,min
= 0.016/t – 0.0084/t
= 0.076/t MN/ m2
Taking a joint efficiency of 0.7 (Double welded butt joint, class 3 construction)
σz (tensile) = fJ
= 93.195×0.7
= 65.2365 MN/m2
65.2365 = 0.0076(1000/t)
t = 0.11maintenance
Maximum compressive load
σz (compressive) = σzwm, max + σzw,max
= 0.0171/t + 0.0475/t
= 0.0646/t
substituting,
σz (compressive) = 0.125E(t/D)
= 0.125(1.88×105)(t/2.1164)
= 1.1104×104t
equating,
t2 = 0.0646/(1.1104×104)
t = 2.4mm
As per IS 2825-1969, minimum corroded skirt thickness = 7mm

Design of skirt bearing plate


Maximum compressive stress between bearing plate and foundation:
σC = W(max)/A + M(wind),max/Z
W(max) =319 kN
A = 3.14(D-l)l
D =outer diameter of skirt = 2.1164m
L = outer radius of bearing plate – outer radius of skirt
M(wind),max = 60.11 kJ
Z = 3.14R2l
R = (D-l)/2
In calculating A, it is assumed that the bearing plate ID is no much less than the ID os the
skirt.
σC = 0.316/(3.14l(2.1164-l)) + 0.06011×4/(l(2.1164-l)2)
Allowable compressive strength of the concrete varies from 5.5 to 9.5 MN/m2. Taking σC
=5.5, and equating,
l = 11.8 mm
Since this is very small, a standard length of l=80mm is chosen.

Thickness of the bearing plate


t(bp) = l(3σC/f)0.5
l = 80mm
σC = Maximum compressive load calculated for l = 80mm
= 0.848 MN/m2
f = allowable stress = 93.195 MN/m2
t(bp) = 13.2mm.
Take a standard thickness of 14mm

As the bearing plate thickness is less than 20mm, gussets are not required.
Rolled angle bearing plate of 14mm thickness is used. (80×80×14)

σ(min) = W(min)/A – M(wind)/Z


= 0.056/(3.14(2.1164-0.08)(0.08) – 0.05613/(3.14(2.1164-0.08)2(0.08))
= 0.1094-0.0539
= 0.0555 MN/m2

R = 0.42D(bp), where D(bp) is the outer diameter of the bearing plate.

j factor = W(min)R/(M(wind),min)
= (0.056)(0.42)(2.1164+2×0.08)/(0.05613)
= 0.9539
j < 1.5
hence, anchor bolts are required.
Design of anchor bolts
P(bolt).N = σ(min)A
= (0.0555)(3.14×(2.1164-0.08)×0.08)
= 0.0284 MN
Hot rolled plain carbon steel is selected for bolts.(f = 53.5 MN/m2)
(area of bolts)f = P(bolt)N = 0.0284 MN
area of bolts = 0.0284/53.5 = 5.3084×10-4 m2
Choose M16×1.5 bolts. Area = 133 mm2
Number of bolts = 4.

Trays
The trays are standard sieve plates throughout the column. The plates have 13476
holes of 5mm dia arranged on a 15mm triangular pitch. The trays are supported on
purloins. The details of the trays are shown in fig 6.3

Nozzles
Nozzles are required for compensation where a hole is made in the shell. The
following nozzles are required:
1.Vapor discharge
Nozzle diameter = 0.1D = 211.64mm.
t = 10mm
2.Reflux Inlet
Nozzle diameter = 126mm
t = 10mm
3.Feed Inlet
Nozzle diameter = 96mm
t = 10mm
4.Reboiled vapor Inlet
Nozzle diameter = 116mm
t = 10mm
5.Liquid Bottoms Outlet
Nozzle diameter = 126mm
t = 10mm
All nozzles are provided with a standard compensation pad of 30mm thickness.
This small compensation is sufficient as the design pressure is low (0.1 N/mm2)

Process Design of Condenser


The following is the detailed design of the total condenser for the Distillation
column. The condenser is operated at the same pressure as that of the column. The vapor
from the column is condensed and sent as reflux and product. Any changes in
compensation are neglected. The optimization of the condenser is done in several
iterative steps. The final trial is presented here. The design methods used here are from
Chemical Engineering by Coulson and Richardson, vol.6

Shell side
Feed = 165.3912 kmol/hr
Average molecular weight Mf = 62.44
T = 409.44 K = 136.29 °C
Mass flow rate = 10327.0265 kg/hr
Heat of vaporization = ∆HV = 53191.7011 kJ/kmol
Heat Load, Q = 2443.7331 kW

Tube Side
Cooling water at 20 °C

1. Heat Balance
Heat load = Q = 2443.7331 = m(4.187)(30-20)
m = 58.3648 kg/s, where m is the cooling water flow rate.

2. LMTD
LMTD = ((136.29-20)-(136.29-30))/ln((136.29-20)/(136.29-30)) = 111.22 °Constraint

3. Heat transfer area.


Choose the overall heat transfer coefficient (U) as 500 W/(m2K)

Q = UA(LMTD)
∴ A = 2443.7331/(500×111.22)
= 43.9442 m2

4. Tubes
Choose 1” OD, 16 BWG tubes of 6ft length laid in a 1.25” Square pitch
Tube OD = 25.4 mm
Tube ID = 22.098 mm
Flow Cross sectional area = 3.8353×10-4 m2
Surface area/tube = 0.1459 m2
Number of tubes (Nt) = 43.94415/0.1459 = 301 tubes
Bundle diameter Db = OD(Nt/k)1/n

For a TEMA 1-2 exchanger with tube pitch = 1.25 (OD),


K = 0.156; n = 2.291
∴ Bundle diameter Db = 690.1 mm
Number of tubes in the central row (Nr) = 2/3(Db/pt) = 2/3(690.1/(1.25×25.4)) = 15

5. Shell side film transfer coefficient


The film temperature Tf is evaluated by an iterative procedure by first assuming a
film coefficient and recalculating the film coefficient using the film temperature. The
film temperature was calculated to be 71.62 °C. The corresponding wall temperature is
118.24 °C
The shell side film transfer coefficient is calculated by using a modified Nusselt’s
equation:
Mass flow rate per unit length (Γh) = W/(L.Nt)
= 10327.0265/(3600×1.83×301)
= 5.2112×10-3 kg/(m.s)

Film transfer coefficient (hC)


hC = 0.95kL[ρL(ρL-ρV)g / (µΓh)]1/3Nr-1/6
= 0.95×0.1935[1074.5432(1074.5432-0.2428)9.81/(4.4×10-3×5.2112×10-3)]1/315-1/6
= 919.0 W/(m2K)

6. Tube Side Heat Transfer Coefficient


Tube velocity (ut) = 58.3648/(995×3.8353×10-4×301/2)
= 1.0162 m/s

Tube side coefficient, hi


hi = 4200(1.35+0.02Tav)ut0.8/di0.2
= 4200(1.35 + 0.02×25)(1.0162)0.8/(22.098)0.2
= 4237.7348 W/(m2K)

7. Overall Heat Transfer Coefficient


Assuming dirt coefficient as 6000 W/(m2K) for both sides and k(wall) = 50 W/mK

1/U = 1/919 + 1/6000 + 25.4×10-3ln(25,4/22.098)/100 + 25.4/(22.098×6000)


+ 25.4/(22.098×4237.7348)
U = 570.4902 W/(m2K)

8. Pressure drop in tube side


Reynold’s number
Re = 22.098×10-3×1.0162×995/(0.9×10-3)
= 24826
Friction factor, jH = 3.8×10-3
Pressure drop, ∆Pt
∆Pt = NP[8 jH(L/D)+2.5]ρut2/2
= 2[8×3.8×10-3×1.83/(22.098×10-3)+2.5](995×1.01622/2)
= 5.1538 kPa
Pressure drop is acceptable.

9. Shell side pressure drop


Select baffle spacing (lb) = shell diameter, baffles 45% cut
For TEMA Pull-Through Floating Head Heat Exchanger,
Clearance = 93 maintenance
Shell diameter = Ds = 783.1 mm
Cross Flow area:
As = (pt-Do)Dslb/pt
= (1-1/1.25)(0.7831)(0.7831)
= 0.1227 m2
Equivalent diameter:
d = 1.27/Do(pt2- 0.785Do2)
= 1.27[(1.25×25.4)2-0.785(25.4)2]
= 25.0806 mm
Mass flux:
GS = w/As = 10327.0265/ (3600×0.1227) = 23.3791 kg/(m2s)
Reynold’s number:
Re = dGS/ (Asµ)
= 25.0806×23.3791/(0.1227×9.4484×10-4)
= 62107
From the correlation, friction factor
jf = 2.5×10-3
Vapor velocity:
us = GS/ρV
= 23.3791/0.2428
= 96.2895 m/s
Pressure drop:
∆PS = 8jf(DS/d)(L/lb)ρuS2/2
= 8×2.5×10-3(783.1/25.0806) (1.8288/0.7831)(0.2428×96.28952/2)
= 0.82 kPa
The maximum allowable pressure drop on the shell side for medium vacuum
operation is 10% of the absolute pressure which is 1 kPa. Hence, the condenser pressure
drop is within the limits.

Summary of Condenser design

Shell outer diameter = 783.1 mm


Bundle diameter = 690.1 maintenance
Number of tubes = 301
Tube OD = 1”
Pitch = 1.25”
Tube length = 6ft
Shell side pressure drop = 0.82 kPa
Tube side pressure drop = 5.15 kPa
Condenser type: TEMA Pull-Through Floating Head 1-2 Heat Exchanger

Mechanical Design of Condenser

Design Temperature = T = 150 °C


Design pressure = 0.1 N/mm2 (external)
Number of tubes = 301
Bundle diameter = 690 mm
Shell diameter = 783 mm
Baffles: Spacing=783 mm, 45% cut

1.Shell Thickness
Material: IS 2825-1969 Grade I plain Carbon steel.
Assume shell thickness = tS = 8mm
L/D = 2.36, D/tS = 97.875
B factor = 2250 (Obtained from graphs in IS 2825)
Maximum allowable pressure = B/(14.22D/t) = 0.16 N/mm2
Hence, the thickness is sufficient.

1. Nozzles
Take inlet and outlet nozzles as 100mm diameter.
Vent nozzle = 25mm diameter
Drain nozzle = 25mm diameter
Relief Valve = 50 mm diameter.
Only the inlet and outlet nozzles need compensation. The compensation required is
minimum and is given by pads of 30mm thickness.

2. Head
Torispherical heads are taken for both ends.
R (Crown radius) = 783 mm
R (knuckle radius) = 78.3 mm
Head thickness = shell thickness = 8mm

3. Transverse Baffles
Baffle spacing = 783mm
Baffle cut = 45%
Baffle thickness = 6mm (standard)

6. Tie rods and spacers


Diameter of tie rods = 10mm
Diameter of Spacers = 8mm

7.Flange Design
Flange is ring type with plain face.
Design pressure = P = 0.1 MN/m2 (external)
Flange material: IS 2004-1962 Class 2 Carbon Steel
Bolting steel: 5% Chromium, Molybdenum Steel
Gasket Material: Asbestos composition
Shell OD = 0.791m = B
Shell Thickness = 0.008m = g
Shell ID = 0.783m
Allowable stress for flange material = 100 MN/m2
Allowable stress of bolting material = 138 MN/m2

(a) Determination of gasket width


dO/di = [(y-Pm)/(y-P(m+1))]0.5
Assume a gasket thickness of 1.6mm
y = minimum design yield seating stress = 25.5 MN/m2
m = gasket factor = 2.75
dO/di = 1.002m

Let di = B+10 = 0.801m


Minimum gasket width = 0.801(1.002-1)/2 = 0.0008m = N
Choose N =0.04. do = 0.809m
Basic gasket seating width = 4/2 = 2mm =b
Diameter at location of gasket load reaction G = di + N = 0.805m

(b) Estimation of bolt loads


Load due to design pressure
H = πG2P/4
= 0.5090 MN
where P is the design pressure
Load to keep joint tight under operation:
Hp = πG(2b)mp
= π(0.805)(0.004)(2.75)(0.1)
= 0.0028 MN
Total Operating Load Wo = H+HT = 0.5118 MN
Load to seat the gasket under bolting condition:
Wg = πGby
= 0.1290 MN
Wo > Wg Hence, the controlling load is Wo = 0.5118 MN

(c) Calculation of Minimum bolting area:


Am = Ao = W/S = 0.5118/S
So = allowable stress for bolting material
Am = Ao = 0.5118/138 = 0.00371 m

Calculation of optimum bolt size.


g1 = g/0.707 = 1.415g
Choose M16×1.5 Bolts
Root area = 133 mm
Minimum number of bolts = 28
Radial clearance from bolt circle to point of connection of hub or nozzle and back of
flange = R = 0.025 m
Bs = 0.075m (Bolt spacing)
C = nBs/π = 0.6690
C =ID + 2(1.415g + R)
= 0.783 +2[(1.415)(0.008)+0.025]
= 0.8566
Choose C = 0.86m
Bolt circle diameter = 0.86m

(d) Flange outside diameter (A)


A = C +bolt dia + 0.02
= 0.896m
select A = 0.90m
(e) Check for gasket width
AbSG / (πGN) = 50.80 < 2y, where SG is the Allowable stress for the gasket material

(f) Flange Moment Calculations


Mo = W1(a1-a3) + W2(a2-a3)
a1 = (C-B)/2 = 0.0385m
a3 = (C-G)/2 = 0.0275m
a2 = (a1+a3)/2 = 0.033m
Mo = 3.0646×10-3 MJ
Mg = Wa3
W =(Ab+Ag)Sg/2
Ab = 28(1.33×10-4) = 3.724×10-4
Ag = Wg/Sg = 0.129/138 = 9.35×10-4
W= 0.3215 MN/m2
Mg = 8.8413×10-3 MJ
Mg > Mo
Hence, Mg is controlling.

(g) Calculation of flange thickness


t2 = M CF Y / (B SF), SF is the allowable stress for the flange material
K =A/B = 0.9/0.783 = 1.15
For K = 1.15, Y = 13
Assuming CF =1
t2 = 1.47×10-3
t = 0.038m
Actual bolt spacing BS = πC/n = (3.14)(0.86)/(28) = 0.0965m
Bolt Pitch Correction Factor
CF = [Bs / (2d+t)]0.5
= 1.1741
t(act) = t(CF) = 0.0446m
Select 45mm thick flange. Both flanges have the same thickness.

8. Saddle Support Design


Material : Carbon Steel
Shell diameter = 783mm
R = D/2
l = 200mm
Torispherical Head: Crown radius = D, knuckle radius = 0.1D
Total Head Depth = 152mm = H
Shell Thickness = Head Thickness = 8mm
ft = 95 MN/m2
Weight of the shell and its contents = 12100 N = W
Distance of saddle center line from shell end = A = 200mm

Longitudinal Bending Moment


M1 = QA[1-(1-A/L+(R2-H2)/(2AL))/(1+4H/(3L))]
Q = W/2(L+4H/3) = 49626 Nm
M1 = 49626(1-0.9649) = 1741.8726 Nm
M2 = QL/4[(1+2(R2-H2)/L)/(1+4H/(3L))-4A/L]
= 24813(0.5671)
= 14071.4523 Nm

Stresses in shell at the saddle


f1 =M1/(π R2 t) = 0.4522 N/mm
f2 = 0.4522 N/mm
f3 =M2/(π R2 t) = 3.65 N/mm
All stresses are within allowable limits. Hence, the given parameters can be
considered for design.
Process Design of Reboiler
The following is the design of a kettle reboiler for the MEG column. The reboiler
is operated as the same pressure as that of the column. Any enhancement of the vapor in
the reboiler is neglected.

Feed = 10238.8134 kg/hr


Quantity of liquid to be vaporized = 9414.9750 kg/hr
Feed Temperature = 433.83 K
Vaporization Temperature = 435.51 K
Average temperature = 434.67 K = 161.52 °C

Properties at the average temperature


Latent heat of vaporization = 550.2221 kJ/kg
Cp = 2.3298 kJ/(kgK)
Critical pressure = 4.66 MPa

(a) Heat Loads


Maximum sensible heat = 2.3298(435.51-433.83) = 3.9141 kJ/kg
Total Heat Load = 3.9141(10328.8134/3600) + 9414.9750/3699(550.2211)
= 1450.1094 kW
Adding 5% losses, Maximum heat load (duty) = 1.05(1450.1094)
= 1522.6148 kW

(b)Heat Balance and number of tubes


Assume U =1000 W/(m2K)
Choose steam at 15.55 bar, T(sat) = 200 °C, λ = 1938.6 kJ/kg
LMTD = [(200.160.68)-(200-162.36)]/ln[(200-160.68)-(200-162.36)]
= 38.47 °C
Area required = (1522.6148×1000)/(38.47×1000) = 39.5798 m2
Select 25mm ID, 30mm OD plain carbon U-tubes, l = 4.8m
No. of tubes = 88
Use square pitch, pt = 1.5(30) = 45mm
Minimum bend radius = 3(OD) = 90mm
From the tube layout, for 88 tubes, outer diameter limit = 502mm

(c) Boiling Heat Transfer Coefficient


q = 1522.6148/39.5793 = 38.47 kW/m2
hnb = 0.104 (PC)0.69 (q)0.7 [1.8(P/PC)0.17 + 4(P/PC)1.2 + 10(P/PC)10]
= 1595.8582 W/(m2 K)
This is acceptable.

(e) Maximum allowable heat flux


σ = 48.4×10-3 N/m
ρL = 1062.5007 kg/m3
ρV = 0.3830 kg/m3
N = 88
K = 0.44
q(max) = k (pt/d) (λ/N0.5) [σg (ρL-ρV) ρV2 ]0.25
= 113.53 kW/m2
Applying a factor of 0.7,maimum flux should not exceed 0.7×113.53 = 79.471 kW/m2
q < q(max). Hence, the design is within permissible limits.

(f). Velocity Check


Bundle diameter = 502mm
Take shell diameter = 2(bundle diameter) ~1010mm
Taking liquid level as 600mm from base (see fig)
Free board = DS – 600 = 1010-600 = 410mm
Width at liquid level = 984.7maintenance
Surface area of liquid = A = 984.7×10-3(4.8/2) = 2.3632 m2
Vapor velocity at the surface = 9414.9750/(3600×A×0.3038) = 3.6427 m/s

Maximum allowable velocity


u(max) = 0.2[(ρL-ρV)/ρV]0.5
= 11.826 m/s
Hence, the velocity is permissible.

Summary of Reboiler design


Type: Kettle Reboiler (TEMA AKT)
Number of Tubes = 88
Tubes: U- Tube, 30mm OD laid on 45mm square pitch.
Shell diameter = 1010mm
Free board = 410mm

Process Design of Reactor


The reactor for the hydrolysis of ethylene oxide is a high pressure adiabatic Plug
Flow Reactor. The rectors initially used were stirred tank reactors or distillation column
reactors. But literature shows that these are rapidly being replaced by Tubular Plug Flow
Reactors. Since there is a large heat of reaction, no catalyst is required and the reaction
goes to completion. The following assumptions are made for the design of the reactor.

(a) Since the molar ratio of water to other components is large, all properties are taken
for pure water.
(b) The heat of reaction if assumed to remain constant over the temperature range.
However, the variation of specific heat is considered.
(c) The pressure in the reactor is maintained such that there is no flashing of ethylene
oxide. The reactor is thus a liquid-liquid reactor.
(d) The operation is purely adiabatic.

Reactor conditions:
Feed: 9.3182 kg/s of reactant mixture containing:
Ethylene oxide: 10.8905 wt%
Water: 89.1095 wt%
Temperature of Feed = 200 °C
Operating pressure = 20atm
Heat of reaction = 21.8 kcal/mol of ethylene oxide
The reaction follows first order kinetics. The kinetics of the reaction are
represented in fig 3.2.

Reactor Design
Design equation:
γXA(-∆HR) = ∫ (CPdT)
Where, γ is the molar ratio of ethylene oxide
XA is the conversion of ethylene oxide
FAo dXA/dV = -rA(XA,T)
Where, -rA(XA,T) = kCA= [koe-E/(RT)][CAo(1-XA)]
Solving the two equations,
V/ FAo = ∫[dXA/{k(T)CAo(1-XA)]
From the kinetic data, the value of the kinetic constants are as follows:
ko = 6.3856×109 min
Activation energy = E = 8.1668×104 J/mol

Solving the integral using Simpson’s 1/3rd rule, we get the reactor volume as:
V =3.0034 m3
For Plug Flow, Re>10000 and L/d > 100
Take L/d = 150.
L = 44.15m, d= 0.2943m
Take L = 45m and hence, d = 0.2915m
The reactor can be arranged in 9 tube lengths of 5 meters each. The line representation of
the PFR is shown in fig 6.4.

You might also like