Ethyleneglycol Design 2520of 2520equipments
Ethyleneglycol Design 2520of 2520equipments
Ethyleneglycol Design 2520of 2520equipments
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
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
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.
Since An = Ac -Ad
3.3352 = 0.785 Dc2 - 0.0408 Dc2
⇒ Dc = 2.1164 m
Calculations:
⇒ hd = k2 [ρg/ρl] Uh2
= 42.16 mm
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
Downcomer Backup
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.
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
∴ ( µM)L = 0.7851 cp
( µM)G = 0.0143 cp
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
λ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
θL = 15.8263 sec
ZL = DC. Cos ( θc/2) = 1.6931m
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
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
= ½ (π/4 . dH2)
(√3/4). lp2
= 0.1008
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
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
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.
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
Column Efficiency
Point Efficiency
λt = mt GM/LM = 0.5572
λb = mb GM/LM = 2.4905 ⇒ λ = 1.5238
∴ NOG = 1 = 2.1914
1/NG + λ / NL
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
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.
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
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
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)
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)
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
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
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)
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) 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.