Production of Mono Ethylene Glycol
Production of Mono Ethylene Glycol
Production of Mono Ethylene Glycol
Session: 2002-2005
Project Advisor
Prof.Dr.Shahid Naveed
Project Members
Fiaz Ahmed Tahir 2002-POLY-1062
M.Rehan Maqbool 2002-POLY-1071
Mohsin Abbas 2002-POLY-1052
Hafiz Abdul Samad 2002-POLY-1063
Abbas Ali 2002-POLY-1067
Bachelor’s Degree
In
CHEMICAL ENGINEERING
With
SPECIALIZATION IN POLYMER ENGINEERING
Name:______________
Name:______________
Our
Beloved Parents,
Respected Teachers,
And Sincere Friends
ACKNOWLEDGEMENT
Chapter # 1
INTORDUCTION
1.0: History 1
1.1: Uses of Mono Ethylene Glycol 2
Chapter # 2
PROCESS SELECTION AND DESCRIPTION
2.1: Process Selection 4
2.2: Chemistry of Reaction 6
2.3: Process Description 8
2.4: Product Specification 10
Chapter # 3
MATERIAL BALANCE
3.1: Material Balance around Pressure Vessel 13
3.2: Material Balance around Reactor 14
3.3: Material Balance around Preheater 14
3.4: Material Balance around Evaporator 15
3.5: Material Balance around Light End Column 15
3.6: Material Balance around EG Refiner 16
Chapter # 4
ENERGY BALANCE
4.1: Energy Balance around Pressure Vessel 17
4.2: Energy Balance around Preheater 18
4.3: Energy Balance around Reactor 18
4.4: Energy Balance around Evaporator 19
4.5: Energy Balance around Light End Column 24
4.6: Energy Balance around EG Refiner 25
Chapter # 5
EQUIPMENT DESIGN
5.1: Pressure Vessel Design 27
5.2: Mechanical Design of Pressure Vessel 28
5.3: Preheater Design 38
5.4: Reactor Design 44
5.5: Evaporator Design 52
5.6: Light End Column Design 71
5.7: EG Refiner Design 105
Chapter # 6
INSTRUMENTATION AND CONTROL
6.1: Objectives 123
6.2: Temperature Measurement Control 124
6.3: Pressure Measurement and Control 124
6.4: Flow Measurement and Control 125
6.5: Control Loops
125
6.6: Control of Distillation Column 127
6.6.1: Flow rate & composition control 127
6.6.2: Temperature control 129
6.6.3: Pressure control 133
Chapter # 7
HAZOP STUDY
7.1: Objectives 134
7.2: HAZOP Study of storage tank for Ethylene Oxide 136
Chapter # 8
MATERIAL SAFETY DATA SHEET
8.1: Chemical Identification 139
8.2: Composition 139
8.3: Hazards Identification 140
8.4: First Aid Measures 141
8.5: Firefighting Measures 141
8.6: Accidental Release Measures 142
8.7: Handling and Storage 142
8.8: Exposure Controls and Personal Protection 143
8.9: Physical and Chemical Properties 144
8.10: Stability and Reactivity 144
8.11: Toxicological Information 145
8.12: Ecological Information 146
8.13: Disposal Information 146
8.14: Transportation Information 146
Chapter # 9
COST ESTIMATION
9.1: Cost index 147
9.2: Total purchase cost major equipments 147
9.3: Fixed Capital Cost 150
9.4: Operating Costs 151
9.5: Production cost
REFERENCES
Reference books 153
Graphs & charts 155
CHAPTER -1
INTRODUCTION
Ethylene glycol was first prepared in 1859 by the French Chemist Charles Wurtz. It was
produced on a small scale during World War I as a coolant and as an ingredient in
explosives. Widespread industrial production began in 1937 when its precursor ethylene
oxide became cheaply available.
When first introduced it created a minor revolution in aircraft design because when used
in place of water as a radiator coolant, its higher boiling point allowed for smaller
radiators operating at higher temperatures. Prior to the widespread availability of ethylene
glycol, many aircraft manufacturers tried to use evaporative cooling systems which used
water at high pressure. Invariably these proved to be rather unreliable and easily damaged
in combat because they took up large amounts of room on the plane, where they were
easily hit by gunfire.
Ethylene glycol is a colorless, practically odorless, low-volatility, low-viscosity,
hygroscopic liquid. It is completely miscible with water and many organic liquids. The
hydroxyl groups on glycols undergo the usual alcohol chemistry, giving a wide variety of
possible derivatives. Hydroxyls can be converted to aldehydes, alkyl halides, amines,
azides, carboxylic acids, ethers, mercaptans, nitrate esters, nitriles, nitrite esters, organic
esters, peroxides, phosphate esters and sulfate esters. This chemistry permits ethylene
glycol to act as an intermediate in a wide range of reactions. Especially significant is
resin formation, including the condensation with dimethyl terephthalate or terephthalic
acid resulting in a polyester resins. The reactivity and solubility of ethylene glycol
provide the basis for many applications. The widespread use of ethylene glycol as an
antifreeze is based on its ability to lower the freezing point when mixed with water. The
physical properties of ethylene glycol-water mixtures are therefore extremely important.
1.1: USES OF MONOETHYLENE GLYCOL
The following is a summary of the major uses of Ethylene Glycol.
1.1.1: Antifreeze
A major use of ethylene glycol is as antifreeze for internal combustion engines. Solutions
containing ethylene glycol have excellent heat transfer properties and higher boiling
points than pure water. Accordingly, there is an increasing tendency to use glycol
solutions as a year-round coolant. Ethylene Glycol solutions are also used as industrial
heat transfer agents.
Mixtures of ethylene glycol and propylene glycol are used for defrosting and
deicing aircraft and preventing the formation of frost and ice on wings and fuselages of
aircraft while on the ground. Ethylene glycol based formulations are also used to de-ice
airport runways and taxiways.
Asphalt-emulsion paints are protected by the addition of ethylene glycol against
freezing, which would break the emulsion. Carbon dioxide pressurized fire extinguishers
and sprinkler systems often contain ethylene glycol to prevent freezing.
Ordinary dynamite will freeze at low temperatures and cannot then be detonated.
Ethylene glycol dinitrate, which is an explosive itself, is mixed with dynamite to depress
its freezing point and make it safer to handle in cold weather. Mixtures of glycerol and
ethylene glycol are nitrated in the presence of sulfuric acid to form solutions of
nitroglycerin in ethylene glycol dinitrate, which are added to dynamite in amounts
ranging from 25 % to 50 %.
1.1.2: Polyester Fibers
The use of ethylene glycol for fibres is becoming the most important consumer of glycol
worldwide. These fibres, marketed commercially under various trade names like Dacron,
Fortel, Kodel, Terylene etc are made by the polymerization of ethylene glycol with
bishydroxyethyl terephthalate (Bis HET).
1.1.3:Resins
Polyester resins made from maleic and phthalic anhydrides, ethylene glycol, and vinyl
type monomers have important applications in the low-pressure lamination of glass fibre,
asbestos, cloth and paper. Polyester fibre glass laminates are used in the manufacture of
furniture, automobile bodies, boat hulls, suitcases and aircraft parts.
1.1.4: Hydraulic Fluids
Ethylene glycol is used in hydraulic, brake and shock absorber fluids to help dissolve
inhibitors, prevent swelling of rubber, and inhibit foam formation.
1.1.5: Capacitors
Ethylene glycol is used as a solvent and suspending medium for ammonium perborate,
which is the conductor in almost all electrolytic capacitors. Ethylene glycol, which is of
high purity, is used because it has a low vapor pressure, is non-corrosive to aluminum and
has excellent electrical properties.
1.1.6: Glyoxal
An important use for ethylene glycol is as the intermediate for the manufacture of
glyoxal, the corresponding dialdeyde. Glyoxal is used to treat polyester fabrics to make
them “permanent press”.
1.1.7: Other uses
Ethylene glycol is used to stabilize water dispersions of urea-formaldehyde and
melamine-formaldehyde from gel formation and viscosity changes. It is used as a
humectant for textile fibres, paper, leather and adhesives and helps make the products
softer, more pliable and durable.
CHAPTER 2
PROCESS SELECTION AND
DESCRIPTION
The formation of these higher glycols is inevitable because ethylene oxide reacts
faster with ethylene glycols faster than with water. The most important variable is the
water-to-oxide ration, and in commercial plants the production of diethylene glycol
(DEG) and triethylene glycol (TEG) can de reduced by using a large excess of water.
The reactor product distribution is essentially unaffected by temperature and
pressure over the ranges (90-200 °C, 1-30 bars) which are normally of commercial
interest. The product distribution is substantially identical for the uncatalyzed and acid
catalyzed reactions, whereas with basic catalysts, the production of higher glycols is
substantially increased.
The rate of the hydration reaction is greatly affected by temperature and by the
catalytic action of acids. Bases are about 1/100 as effective as equivalent concentrations
of acids in promoting the hydration reaction, though they more strongly catalyze the
glycols, as noted above.
The use of acid catalysts in the hydration reaction permits it to be operated at lower
temperature and pressure than the uncatalyzed reaction, but makes the reaction solution
highly corrosive and requires acid removal. Sulfuric acid will partly combine with the
glycol product and is best removed using anion-exchange resin. The potential use of other
acids such as oxalic and trichloroacetic has been proposed to reduce the neutralization
problem. In 1959, the acid-catalyzed reaction reportedly was the one most commonly
used, but by 1968 it was concluded that there was no readily available evidence to show
that was still the case. Currently, it is believed that most plants use uncatalyzed hydration
reactions.
An interesting approach to avoiding the problems associated with the use of free
acids involved ion-exchange resins in heterogeneous solid-liquid systems. Though acidic
resins were reported to be only one-tenth as active as the same number of equivalents of
sulfuric acid, their use was bed of resin was proposed. There is no indication from the
literature that such processes have been employed commercially.
CHAPTER-3
MATERIAL BALANCE
(200,000 TPA MEG)
If we use 12 fold molar excess of water, then from literature selectivity of different
products is as under:
MEG = 85 %, DEG = 14 %, TEG = 1 %
Reactions:
C2H4O+H2O OHCH2CH2OH (MEG)
Suppose MEG is concentrated from 19.73 % to 60 % in this long tube vertical falling
film evaporator.
CHAPTER-4
ENERGY BALANCE
[Rate of Accumulation of Energy within system =Rate of Energy entering the system –
Rate of energy leaving the system + Rate of Energy generation]
For steady state system no accumulation of mass or energy with in system, so by
modifying above equation, and applying across each equipment, the energy
balance of each equipment is as under.
For case of energy balance across each equipment to determine the enthalpy of streams
we used reference temperature equal to 25 degC
Datum Temperature = 25 0C
Heat of feed entering = Q = mCp∆T
Heat of Formation
EG -386 kJ/mole
EO -129 kJ/mole
Water -242 kJ/mole
DEG -571 kJ/mole
Reaction 1:
H2O + EO EG
∆Hr1 = [(-386)] – [(-129)-1(-242)]
= -15 kJ/mole
= -15,000 kJ/kmole
Reaction 2:
H2O + 2EO DEG
∆Hr2 = [(-571)] – [2(-129)+(-242)]
= -71 kJ/mole
= -71,000 kJ/kmole
Total heat of reaction = ∆Hro = [(-15000)*411.4] + [(-71,000)*(40)]
- ∆Hro = 9011000 kJ
Total Heat In = 84542103 + 9011000
= 93553103 kJ
Heat of product leaving = 129286 * 4.24 * (468-298)
= 93189349 kJ
Difference in heat = 93553103-93189349
= 363754 kJ
Temperature change = 363754/129286*4.24
= 1oC
This temperature change is very negligible so cooling required.
EQUIPMENT DESIGN
In this chapters introduction to different equipments of plant along with their designing
procedure and specification sheets are presented.
Weight of Vessel
Wv = 240 * Cv * Dm (Hv + 0.8 Dm) t
Cv = 1.08 for vessel
Hv = Height of vessel = 6.40 m
g = 9.81 m/sec2
t = wall thickness = 13.25 * 10-3 m
Dm = mean vessel diameter = (Di + t * 10-3) m
Dm = (2.56 + 13.25 * 10-3)
= 2.57325 m
Wv = 74753 N
Fluid Weight
Wf = 129324 * 10/60 = 21554 kg = 21554 * 9.8
= 211229 N
Fittings,plateforms weight
= 360 + 150 + 1.7 * 1000 * 5.147
= 9260 N
Total Weights = 9260 + 211229 + 74753
= 295242 N
Primary Stresses
δw = W / π(Di + t) * t
= 295242 / π (2.56 *103 + 13.25) * 13.25
= 2.756 N/mm2
c) Bending Stresses
Wind Loads
Mx = wx2/2
Assume wind pressure = 1280 N/m2
Dm = 2.57325 m
Loading per linear meter = 1280 * 2.57325
= 3294 N/m
Bending moment at bottom = Mx = 3294/2 * 6.42
= 67456 Nm
Bending Stresses
δb = ± M / Iv (Di/2 + t)
Iv = second moment of area
Iv = π / 64 (Do4 – Di4)
Do = 2586.5 mm
Di = 2560 mm
Iv = π / 64 ((2586.5)4 – (2560)4)
= 8.866 * 1010
δb = {67456 Nm / 8.866 * 1010}* (2560/2 + 13.25) * 103 mm/1m
= ± 0.984 N/mm2
d) Torsional Shear Stresses
τ = I/Ip (Di/2 + t)
For preliminary vessel design. These are neglected.
Principal Stresses
δ1 – δ2 = 86-84 = 2 N/mm2
δ1 – δ3 = 86 – (-0.446) = 86.446 N/mm2
δ2 – δ3 = 84 – (-0.446) = 84.446 N/mm2
Maximum Stress Intensity = δ1 – δ3 = 86.446 N/mm2
Design Stress = f = 120 N/mm2
Maximum Stress Intensity < Design Stress
Check Elastic Stability (Bubling)
Critical Bubling Stress = δc = 2 * 104 (t/Do)
δc = 2 * 104 (13.25/2586.5)
= 102 N/mm2
Maximum Compressive Stress is well below critical bubling stress, so design is
satisfactory.
Vessel Supports
Brackets and lugs are used.
As pressure vessel is vertical and also bending moment is negligible, so we will use
brackets.
Bolt Flanged Joints
Welding,neck flanges is used for process vessels or process equipments.
Gaskets
Pressure = 8 atm
T = 128 oC
Corrosiveness = negligible
So Teflon or compressed asbestos gaskets are selected
.
Flange Types
Flange Faces
(a)Full face (b) Gasket within bolt circle (c) Spigot and socket (d)Ring type joint
Standard Flange Design
T1 - t 2 - T2 - t1
LMTD T - t
ln 1 2
T2 - t1
LMTD = 43 K
Estimated UD
UD= 1448 W/m2.K
A= 195 m2
Tube sheet layout
Tube length = 16 feet
Tube Dia = 1 inch
Pitch = 1 ¼ inch triangular
BWG = 14
Surface area per lin ft = 0.2618 ft2
Surface Area of 1 tube = 0.2618 x 16 = 4.19 ft2 =0.3893 m2
Number of tubes = Nt = Total Area /Surface area of 1 tube
Nt = 501 , Nearest tube count = 522
Number of passes = 2
5.3.1:SPECIFICATION SHEET
PREHEATER
Identification: Item Preheater
No. Required = 1
Function: Preheating of the feed to the required temperature of the Reactor
5.4:DESIGN OF REACTOR
5.4.1:REACTOR SELECTION
There are three types of Reactors.
PLUG FLOW REACTOR
CONTINUOUS STIRRED TANK REACTOR
BATCH REACTOR
Reactor is selected on the basis of following parameters.
1- Conversion
2- selectivity Reaction
rate
3- productivity
4- safety Volume
5- economics Volume flow
of reactor
6- availability rate
7- flexibility
8- compatibility with processing
Model
9- energy utilization
Outlet
10- feasibility concentra
11- investment operating cost tion
Inlet
12- heat exchange and mixing Concentr
ation
Contacting
pattern
rA kCA n
CA CA 0 (1 XA)
rA kCA 0 (1 XA)
By integration,
1
V / FA 0 ln(1 XA)
kCA 0
We need C A0 and F ;
A0
(2.4870 E 01 )
02
6.4713 E
1.836
Density of Ethylene Oxide ( EO ) = 468 0.2696
0.26024[1(1 469.15 ) ]
= 22.5 m3
Thus
Volume of PFR= 22.5 m3
VOLUME OF CSTR
As we know the formula for calculating the volume of CSTR
i.e.
XA
V / FA 0 =
rA
FA 0 XA
V
kCA 0 (1 XA)
Assume 8 passes:-
Length of each pipe = 104/8
= 13 m (bend is also included)
Length of pipe = 13 m
Total no. of pipes = 112 * 8
= 896 pipes
Thus Total pipe passes = 8
Bundle Diameter:
Selection: U tube and Fixed.
Db = Do(nT/K1)1/n1
n1 = 2.675
K1 = 0.0365
Triangular Pitch
Tube Pitch = PT = 1.25 * Do
Do = Outside Tube Diameter in “mm”.
Do = 2.38 in = 60.452 mm
nT = 896
Db = 60.452(896/0.0365)1/2.675
Db = 2646 mm
Minimum Wall Thickness of the Pipe:
Min. wall thickness = {P (inside radius) / SEj – 0.6P} + Cc
Where
P = Max. allowable internal pressure, kPa
S = Max. allowable working stress, kPa
Ej = Efficiency of joints expressed as a fraction
Cc = Allowance for corrosion, m
Check for 2”, Schedule No. 80
Putting the values of P, inside radius and other values from Table 12.10;
Min. wall thickness = t = {1519.875(1.939/2*12*3.281)}/{94,500(0.85)-0.6(1519.875)}
+ 0.0032
By solving the above Equation, we get
Min. wall thickness = 0.00367 m = 0.1445”
Actual Wall Thickness = 2.38-1.939/2
= 0.2205
Safety Factor = .2205-.1445/.1445
= 50%
Nominal Pipe Size 2” and Schedule No. 80 is selected.
5.4.4:SPECIFICATION SHEET
Identification
Item Reactor
Item No. R-1
No. required 1
Function: Production of Mono ethylene Glycol by Hydration of ethylene oxide.
Operation: Continuous
Type: Non-Catalytic
Multi tube
Chemical Reaction:
5.5:EVAPORATOR DESIGN
As we use excess of water to increase the selectivity of Monoethylene glycol ,so after
completion of reaction and obtaining required selectivity in the reactor excess water is
removed. About 90% of water is removed in evaporator by using triple effect evaporator.
EVAPORATION
In the evaporation process, concentration of a product is accomplished by boiling out a
solvent, generally water. It is a unit operation that is used extensively in processing foods,
chemicals, pharmaceuticals, fruit juices, dairy products, paper and pulp, and both malt
and grain beverages. Also it is a unit operation which, with the possible exception of
distillation, is the most energy intensive.
While the design criteria for evaporators are the same regardless of the industry involved,
two questions always exist: is this equipment best suited for the duty, and is the
equipment arranged for the most efficient and economical use? As a result, many types
of evaporators and many variations in processing techniques have been developed to take
into account different product characteristics and operating parameters.
MODE OF EVAPORATION
The user needs to select one or more of the various types of evaporator modes
that were described in the previous section. To perform this selection, there are a
number of ‘rules of thumb’ which can be applied.
1) Falling film evaporation:
either plate or tubular, provides the highest heat transfer coefficients.
is usually the mode chosen if the product permits.
will usually be the most economic.
is not suitable for the evaporation of products with viscosities over 300cp.
is not suitable for products that foul heavily on heat transfer surfaces
during boiling.
2) Forced circulation evaporators:
can be operated up to viscosities of over 5,000cp.
will significantly reduce fouling.
are expensive; both capital and operating costs are high.
3) Paravap evaporators:
are suitable for viscosities up to 10,000cp for low fouling duties.
are suitable for very high viscosities, i.e., over 20,000cp, usually the
only suitable evaporation modes are the wiped film and thin film systems.
FILM EVAPORATORS–PLATE OR TUBULAR
1)Plate evaporators:
Provide a gentle type of evaporation with low residence times and are
often the choice for duties where thermal degradation of product can occur.
Often provide enhanced quality of food products.
Require low headroom and less expensive building and installation costs.
2)Tubular evaporators:
Are usually the choice for very large evaporators.
Are usually the choice for evaporators operating above 25 psia (1.7 bar).
Are better at handling large suspended solids.
Require less floor space than plate evaporators.
Have fewer gasket limitations.
Fig. (Evaporator Selection)
.
From above discussion and fig. Most suitable and economical evaporator for our
application is long tube vertical (LTV) falling film Evaporator
5.5.3:CALENDRIA DESIGN
1ST Effect
Selection: Long Tube Vertical Falling Film Evaporator
Mode of feed supply is Forward Feed.
Design:
If uD1 , uD2 , uD3 are the corresponding overall heat transfer coefficients and A1, A2
and A3 are the areas of heat transfer required, then we may write
Effect 1: Q1 = ws λs = UD1A1(Ts-Tb1)
Effect 2: Q2 = ws1 λs1 = UD2A2(Tw1-Tb2)
Effect 3: Q3 = ws2 λs2 = UD3A3(Tw2-Tb3)
Assumptions:
1. Equal heat transfer rates in all effects. From this rate of evaporation is same in all
effects.
2. Equal areas of heat transfer for all effects.
From assumption (1);
Q1 = Q2 = Q3
UD1A1ΔT1 = UD2A2ΔT2 = UD3A3ΔT3
A1 = A2 = A3
UD1 ΔT1 = UD2 ΔT2 = UD3 ΔT3
Assume values of UD1 and UD2 and UD3
1
T 1 uD1 T
1 1 1
uD1 uD 2 uD 3
uD 1
T 2 T 1
uD 2
uD 1
T 3 T 1
uD 3
From these ΔT1 , ΔT2 and ΔT3 ; select the pressures of three effects. Then calculate values
of UD1 ,UD2 and UD3
To select the pressure of evaporators is an iterative procedure. The final iteration results
are given here;
Pressure ΔP (kPa) Steam or λ (kJ/kg)
(kPa) vapor
temperature
(K)
Steam chest, 2026.5 ------- ts = 486 λ s = 1886
1st effect
Steam chest, 962 1064.5 t1 = 451 λ 1 = 2020
2nd effect
Steam chest, 304 658 t2 = 401 λ 2 = 2179
3rd effect
Vapor to 35.46 268.54 T3 = 346 λ 3 = 2326
condensor
By Heat Balance
1st Effect:
ws λs + wF CF (tF-ts1) = w1 λ1 --------------------(A)
2nd Effect:
w1 λ1 + (wF - w1) C1 (ts1-ts2) = w2 λ2 --------------------(B)
3rd Effect:
w2 λ2 + (wF - w1 – w2) C2 (ts2-ts3) = w3 λ3 --------------------(C)
ts1 ,ts2 and ts3 are saturation temperatures of liquor in effects1,2,3 respectively.
CF ,C1 and C2 are heat capacities of feed, liquor in effects (1) and (2).
These are approximated for 1st time by using equal evaporation in each effect and
calculating weight fractions of liquors in different effects. After which it is iterated and
revised. Final iteration results will be discussed here.
wT = w1 + w2 + w3 ---------------------(D)
by solving simultaneously equations (A), (B), (C) and (D), following relations are
obtained.
w1
w
t 3 wf 2C 2(ts 2 ts 3) wf 2C 1(ts1 ts 2 ) wfC 1C 2(ts1 ts 2)(ts 2 ts 3) wf 3C 1(ts1 ts 2 )
2
w2 = 22194 kg/hr
w3 = 33927 kg/hr
ws = 27931 kg/hr
233
ua =0.2 Ta =448 oK
Lv = latent heat of vaporization = *1000
Estimation of Properties (Vapour):
As in all these effects, water is evaporated. So for vapor properties are properties of pure
water are used and calculated as under.
4.60*104 N * M .W .0.5 * Pc 2/ 3
Viscosity: v *10
3
Tc 1/ 6
And N 0.0003400* Tr 0.94
A = 3.470 B = 1.450*10-3 C=0 D = 0.121*105
8.314*1000
Cpv ( A BT CT 2 D / T 2 ) * J/kg.oK
18
kv = A BT CT 2
PM
v
RT
Calculation of UD1:
Heat transfer coefficients inside the tubes”
by using above correlation at
T = 470.5 oK
Following properties are obtained
Liquid Properties:
l 1.8623*104 Pa sec
Cpl 4.2057 *103
kl .840W / m.K
l 916.4kg / m3
Lv 2020000 J / kg
Cpl l
Pr l
kl
l 0.0813
Vapour Properties:
v 4.4293kg / m3
Cpv 1943 J / kgK
v 1.274*105 Pa sec
kv 0.1161
g 9.8m / sec 2
Csf 0.013
n 1( water )
Te Ts ts1 486 455 31K
putting all these values in Equation (1)
Heat Flux = qs 3.133*107 W / m 2
And from Equation (2)
Critical Heat Flux = q max = 3.2773*106 W/m2
As
qs > qmax
So by using Equation (3)
d = ¾ “ = 0.019 m
hb = 1255 W/m2K
hio = 1255 W/m2K
Heat Transfer Coefficient inside Shell:
As we are using steam for heating purpose, the condensing coefficient of steam
ho = 5700 W/m2K
Rd = 0.001/5.6783
= 0.000173 m2 K/W
hio * ho 1255*5700
UcI
(hio ho ) (1255 5700)
UcI 1028
1 1
Rd
UdI UcI
= (1/1028)+0.000173
UdI = 873 W/m2K
UdI = 3142 kJ/m2. K.hr
2nd Effect:
Calculation of UD2:
For calculation of UD first we will estimate liquid and vapor properties by using the
correlations discussed above
Liquid Properties (at average temperature)
t 1 ts 2
T = 451+405/2 = 428 oK
2
all properties are calculated by correlations given before.
Weight fraction of MEG = x = 0.3337
l 965kg / m3
Cpl 3988.3J / kgK
l 3.2428*104 Pa sec
kl 0.4835W / m.K
l 0.0872
Vapor Properties:
v 1.5378kg / m3
Cpv 1919 J / kgK
v 1.1422*105 Pa sec
kv 0.1122W / m.K
d = 0.019 m n=1 (water)
Csf 0.013
by using equation (1) and (2);
qs = 8.278*106 W/ m2
qmax = 1.8922*106 W/ m2
qs > qmax
so film boiling is there. By applying equation (3)
hio = hb = 888 W/ m2 K
Shell Side H.T Coefficient
Considering coefficient for steam = 5700 W/ m2 K
ho = 5700 W/ m2 K
uc = (hio * ho) / (hio + ho)
= (5700*888)/(5700+888)
= 768 W/ m2 K
1/ uD 1/ uc Rd
= 1/768 + 0.000173
uD = 678 W/ m2 K
uD2 = 2441 kJ/hr m2 K
Calculation of uD3:
For calculation of UD first we will estimate liquid and vapor properties by using the
correlations discussed above
Liquid Properties:
t 2 ts 3
T
2
= 401+350/2
= 375.5 oK
weight fraction of MEG = x = 0.60
l 6.6795*104 Pa sec
l 1032kg / m3
Cpl 3483J / kgK
kl 0.3728W / mK
0.0553
Vapor Properties:
(At T=ts3)
v 9.9577 *104 Pa sec
v 0.2044kg / m3
Cpv 1894 J / kgK
kv 0.1076W / mK
d 0.0190
Csf 0.0130
n 1
qs = 6.94*105 W/ m2
qmax =6.59*105 W/ m2
qs > qmax
so
hio = hb = 585 W/ m2 K
uc = 530 W/ m2 K
uD = 485 W/ m2 K
uD3 = 485 W/ m2 K
uD3 = 1746 kJ/hrm2K
Area Calculation:
A1 = (27931*1886)/3142*(486-455)
= 541 m2
A2 = (30652*2020)/2441*(451-405)
= 551 m2
A3 = (22194*2179)/1746*(401-350)
= 543 m2
A1 A2 A3
A = A1 + A2 + A3 /3
= 541+551+543/3
A = 545 m2
Thus area of each effect of triple effect evaporator will be 545 m2
¾ “O.D, 14 B.W.G tubes, L = 9.1435 m (30 ft)
Surface of one tube = 0.547 m2
nt = 545/0.547 = 996 corrected nt = 1102 m2 corrected area= 603 m2
Safety factor = 603-545/545
= 11 % (satisfactory)
5.5.4: CONDENSER DESIGN
Detailed selection criterion and design of condenser is given in design of condenser for
light end column.
Specification sheet is given here.
SPECIFICATION SHEET
CONDENSER
Operation: Continuous
Type: 1-1 Horizontal Condenser
Shell side condensation
Heat Duty = 78812421 KJ/ hr
Tube Side: Tubes:0.01905 m diam. 14 BWG
Fluid handled :Cold water No of Tubes =730
Flow rate = 942056 Kg /hr Length of Each Tube =5m
Pressure = 101325 N/m2 (1 Atm) No of Passes =1
Temperature = 25oC to 45 oC Pitch :0.0238 m triangular
pressure drop = 0.072 psi
Shell Side: Shell: 0.736 m . diam.
Fluid handled = Water vapours No of Pass = 1
Flow rate = 8475.99 lb/hr Baffles spacing = 0.61 m .
Temperature 70 oC Pressure drop =0.18 psi
Utilities: Cold water
Ud assumed = 10875 KJ /hr m2 K Ud calculated =10691 KJ /hr m2
Remaining water from evaporator is pumped to the light end column where it is
separated, to obtain our required purity. In this column about 99.9% of water is removed.
To obtain such high purity we have selected distillation column.
Designing Steps
Plate columns are designed to handle wide range of liquid flow rates without
flooding than packed column.
For large column heights, weight of the packed column is more than plate column.
Periodic cleaning is easy in plate columns as compared to packed columns.
When temperature change is involved, packing may be damaged.
If the liquid causes fouling or contains solids, it is easier to make provision for
cleaning in plate column
Selection Of Plate Column
More positive contact by repeated mixing and separation.
For the liquid causes fouling. So man ways are installed on plate for cleaning
purposes.
For the large temperature changes.
Plates are mostly used for large diameter more than 0.6m.
Easy to maintain than that of packed column.
Suitable for very low liquid rates.
Selection Criteria of Trays
Cost
Cost of plate depends upon material of construction used.
For mild steel, the ratio of cost between plates is
Sieve plate : valve plate : bubble-cap plate
3 : 1.5 : 1.0 .
Capacity
Sieve tray > valve tray > bubble-cap tray
Operating Range:-
It is the ranges of liquid and vapor flow rates which must be above the weeping
conditions and below the flooding conditions. Operating range flexibility comparison is
bubble cape tray > valve tray > sieve tray
For good design, sieve plate gives satisfactory operating range.
Pressure drop
Bubble-cap tray > valve tray > sieve tray
Choice of Tray Type (Sieve Tray)
Sieve plates are lighter in weight and less expensive. It is easier and cheaper to
install.
Pressure drop is low as compared to bubble cap trays.
Maintenance cost is reduced due to ease of cleaning.
If properly designed, sieve tray gives desired separation
5.6.2:COLUMN DESIGN
Material Balance:
Feed
Wt Wt Wt Moles Moles Mole
(kg/hr (kg/sec) (%) (Kgmoles/hr) (kgmoles/sec) s
) (%)
Water 12504 3.4733 29.41 694.6667 0.19296 60.55
MEG 25508 7.0856 60.00 410.9554 0.11415 35.82
DEG 4201 1.1669 9.68 39.5872 0.01099 3.45
TEG 300 0.0833 0.71 1.9977 0.00055 0.174
TOTAL 42513 11.80917 100.00 1147.2070 0.31867 100.0
0
Average Molecular weight of the Feed =37.0578
Distillate
Wt Wt Wt Moles Moles Moles
(kg/hr) (kg/sec) (%) (Kgmoles/hr) (kgmoles/sec (%)
Bottom
Wt Wt Wt Moles Moles Moles
(kg/hr) (kg/sec) (%) (Kgmoles/hr) (kgmoles/s (%)
ec
Water 10 0.0278 0.03 0.55556 0.000154 0.123
MEG 25408 7.0578 84.92 409.3443 0.113707 90.67
DEG 4201 1.1669 14.04 39.5872 0.010996 8.768
TEG 300 0.0833 1.00 1.9977 0.000555 0.44248
TOTAL 29919 8.3108 100.00 451.4848 0.125412 100.00
Average Molecular weight of the Bottom = 66.2680
NUMBER OF PLATES
Fenske equation
X X
ln LK HK
X HK d X LK b
Nm
ln LK
0.997684 0.90666
ln
0.002316 d 0.00123 b
ln(62)
=3.069 plates with Reboiler
q = 1 (feed at its bubble point)
α LK xf LK α HK xf HK
1 q
α LK θ α HK θ
(62 * 0.60553) (1 * 0.358)
1 1
(62 θ) (1 θ)
α LK xD LK α HK xD HK
Rm 1
α LK θ α HK θ
Rm = 0.019698 R = 0.21
.0.566
N Nm R Rm
N ! 0. 75 1
R 1
.0.566
N 3.069 0.21 0.019698
0.75 1
N ! 0.21 1
N= 7
2
N B x xbLK
log D 0.206 log fHK
D x fLK x dHK
NB
ND
0.6325 N D 0.6325 N B
NB
N total N D N B 7 0.6325 1 N B
N D 3
NB 4
ND=3 / 0 .60=5
Column Diameter
The plate diameter is calculated based on the flooding considerations
Maximum liquid flow rate in rectifying section
Ln = D Rm = (0.21 0.19326) = 0.0405846 kmol/sec
Maximum vapor flow rate in rectifying section
Vn = Ln + D = (0.0405846 + 0.19326) =0.2338446 kmol/sec
Maximum liquid flow rate in stripping section
Lm = Ln + F = ( 0.0405846+ 0.3186682) = 0.3592528 kmol/sec
Maximum vapor flow rate in stripping section
Vm = Lm – B = (0.3592528– 0.125412456)) =0.2338446 kmol/sec
Density of liquid in the stripping section
ρLm = 1052kg/ m3
Density of the vapour in the stripping section
ρVn = 0.557 kg/ m3
Density of the liquid in the rectifying section
ρVm = 976 kg/ m3
Density of the vapour in the rectifying section
ρLn = 0.190694 kg/ m3
0.358665903 55.46
Bottom 0.002 m3/sec
1052
Calculate FLV = liquid vapor flow factor
LW V
FW
VW L
0.0405846 0.19694
FWTop
0.2338446 976
= 0.002489427
Bottom Calculation
0.352528 0.557
FLV Bottom
0.2338446 1052
= 0.035389951
and for a tray spacing of Ts=609 mm.
Csb=0.0105+(8.127/10000)*(Ts^0.755)*EXP(-1.463*FLG^0.842)
Flooding parameter
Csb, flood = 0.104801008 m/s.
Now,
Unf = Csb, flood × (σ/ 20) ^0..2 [( l - g) / g] 0.5
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)
Unf=Csb,flood×(42.2/20)^0.2[(1052 - 0.557)/0.557] 0.5
= 5.283950816 m/sec
Flooding = 80 %
Let
Actual velocity, Un= 0.8×Unf
i.e., Un = 4.227 m/s
Net area available for gas flow (An)
Net area = (Column cross sectional area)-(Down comer area.)
An = Ac - Ad
An = (0.235133*55.46 / (0.557))/4.4227 = 5.53 m2
Let Ad = 15 % of the Ac (cross sectional are) Ac= An / 0.85
Ac = 5.53/0.85 Ac = 6.51m2
D = √(4 * Ac/pi)
= √(4 * 6.51/(22/7)
= 2.88 m
Calculation at the Top
0.0405846 0.19694
FWTop = 0.002489427
0.2338446 976
=√(4 * 6.51/(22/7)
= 2.5598 m
Column Area Ac d2
4
Ac = 6.51 m2
Net area An = Ac – Ad
Weir length
Check Weeping
Maximum liquid rate=12.72 kg/sec
Minimum liquid rate at 70% turn down 0.7 12.72 = 8.90 kg/sec
2/3
12.72
how 750
1052 0. 9765
= 40.14 mm liquid
2/3
8.90
Minimum how 750 = 31.63 mm liquid
1052 0.9765
K 0.9 25.4 a n
K2 = 30.3 U min 2
v 1 / 2
30.3 0.9 25.4 3.175
U min = 13.79 m/s
0.557 1 / 2
minimum vapou vol rate
Actual minimum vapour velocity
Ah
0.7 23.22
= 35.67 m/s
0.4557
2
50.95 0.557
hd 51 = 41.046 mm liquid
0.93 1052
Residual Head
12.5 10 3
hr 11 .88 mm liquid
1052
ht = hd + (hw + how) + hr
ht = 114.556 mm liquid
Check Entrainment
FLV(base) = 0.035
No of Holes
Height of Column
# Of plates = 18
Height of column = (#of plates - 1) × space between each plate) + (space for
disengagement on top and bottom)
Identification:
Item Distillation column
Item No.
No. required 1
Tray type Sieve tray
Function: Separation of Water for the Recovery of MEG.
Operation: Continuous
Material handled
Feed Top Bottom
Quantity F =1147.2070 D= 695.72 B =451.48
kg.moles/hr kg.moles/hr
kg.moles/hr
Design Data
As we are operating distillation column at vacuum , high rate of vaporization are required
and also viscosity is low .therefore; the best choice is kettle reboiler
Calculate Heat Duty – Q
Shell Side: (cold) Tube Side: (hot)
Water and MEG Steam
Inlet Temp of cold fluid = t1 = 76 oC
Outlet Temp of cold fluid= t2 = 90 oC
Pressure at shell side =0.300 atm=228 mm of Hg
Inlet Temp of hot fluid = T1 = 125oC
Outlet Temp of hot fluid = T2 = 125oC
LMTD = ((125-76)-(125-90)) / Ln ((125-76)/(125-90))
=35 (41.60)
Vapour flow rate (to be evaporated) = Vm = 0.2338446 kgmoles / sec
Latent heat of Water at 76 oC = λ =38865.6 KJ/kg mole
Latent heat of MEG at 76 oC = λ = 61010.488 KJ/kgmole
For designing we use latent heat of MEG
Heat to be provided Q = Vm * λ
= 0.2338446 * 61010.488
= 14266973.32 J/sec
= 5136113.95 KJ/hr
Latent Heat of Steam at 125 oC =2187 KJ/Kg
Mass Flow rate of the Steam = 5136113.95 / 2187
= 23484.72 kg/hr
Over Design heat transfer Coefficient
Assumed Ud=752.3353 W /m2 K
Estimated Area
A=Q/(Ud*LMTD)*1000/3600
A = 541.8167 m2
Tube side Calculations
Tubes: ¾ in. dia . 16 BWG
Effective Length of the tubes L=9m
O.D of the tubes=0.75 in = 0.0625 ft = 0.01905 m
I.D of the tubes =0.584 in = 0.04867 ft = 0.01484 m
Surface area per linear ft of tube
=0.1963 ft2/ft= 2.3556 in2/in =0.0598m2/m
Surface area =at = 9*0.0598 = 0.547 m2
No of tubes
nt=A/at =541.8167/0.547= 990.5241
corrected number of tubes nt =1068
No of Tubes pass
No of tube passes tp = 2
Corresponding Shell ID = 35 in =0.889 m
Corrected A
A=nt*a t= 1068*0.547= 584.1960m2
Cross sectional Area
Act= nt * atc / np
Atc = 1068 *(0.268*2.542/1002 ) /2
= 0.0923 m2
Mass Velocity
Gt = w/ Atc =23484.72/0.0923=254439
Corrected Ud
Ud = Q/(Acor*LMTD)*1000/3600
Ud = 697.7586 W/m2 K
Inside tube heat transfer coeffience based on the outer dia
hio=8521.812252 W/m2 K
Shell Side Calculations
Shell Diameter= 35 in = 0.889 m
Pitch =Pt=15/16 in =0.0238 m
B=ID=0.889 m
Number of baffles = nb=L/B = 9/0.889=10.12
Corrected baffle number = 10
O.D of the tubes=0.75 in = 0.0625 ft = 0.01905 m
Clearance C=pt-OD=0.0238 - 0.01905 = 0.00475 m= 0.1875 in
B=L/nb
B= 9/10= 0.900= 35.43
Shell Cross sectional Area
as =((pi*ID2/4)-(nt*pi*OD2/4))
as =((3.14*0.8892/4)-(1068*3.14*.019052/4))
as = 0.3163 m
Mass Flow rate =wf=53391.179 kg/hr
Mass Velocity
Gs=wf/as;
Gs = 53391.179 / 0.3163
Gs =168799.1748 kg/m2hr
Equivalent Diameter
De= 4*as / /(nt*pi*OD)
De= 4*0.3163/(1068*3.14*0.01905)
De=0.0198 m
Bubble point Temperature ta=90 0C
μ=0.0011 kg/m sec
Reynold Number
Re.s = De Gs /
Re.s = (0.0198*168799.1748)/(0.001*3600)
Re.s = 843.7998
JH Factor
jH = 15
ho = jH *k/D (c / k) 1/3 ( / w)0.14
( / w)0.14 = ?
Dirt Factor
Rd= 0.001/5.6812;
Rd = 0.000176 m2 K/W
Ud =697.7586 W/m2 K
Estimated Over all Clean Coefficient
Uc=Ud/(1-Ud*Rd)
Uc = 697.75/(1-697.75*0.000176)
Uc = 795.4556 W/m2k
hio= 8521.81 W/m2 K
Outside film Coefficient
Estimation
ho=Uc*hio/(hio-Uc)
ho=795.4556*8521.81/ (8521.81-795.4556)
ho =877.3506 W/m2 K
Wall temperature
Bubble point Temperature =ta=90 0C
Steam Temperature =Ta=125 oC
μ=1.1 Cp = 0.0011 kg/m sec
tw=ta+hio/(hio+ho)*(Ta-ta)
tw = 90+8521.81/(8521.81+877.3506)*(125-90)
tw =121 oC
μw=0.9 Cp = 0.00090 kg/m sec
k=0.645 W/m k
c = 4153 J /kg k
jH = 15
Corrected Out Side Film Coefficient
ho=Jh*k/De*(c* μ / k)(1/3) * ( μ / μw)0.14
ho=15*0.645/0.0198 *(4153* 0.0011/ 0.645 /)(1/3) * 0.0011 / 0.0009)0.14
ho= 965.4099 W/m2 K
Corrected Overall Clean Coefficient
hio= 8521.81 W/m2 K
Uc=hio*ho/(hio+ho)
Uc=8521.8*965.499/(8521.81+965.4099)
Uc = 867.1708 W/m2 K
Rd = 0.000176 m2 K/W
Corrected Over all Dirt Coefficient
Ud=Uc / (1+Uc*Rd)
Ud=867.171 / (1+867.171*0.000176)
Ud= 752.3353
Bundle Dia
Db = do (Nt / K1)1/n1
For Triangular Pitch: K1 = 0.249 N1 = 1068
Db = ¾* (1068 / 0.249)1/2.207 = 33.182 in = 0.8428 m
Shell-bundle clearance :
C = shell inside dia – bundle dia = 0.889- 0.8428 = 0.0462m
Heat Flux
Heat Load = 51361104 kJ / hr = 14266973.33 J / sec
Area = 542 m2
Heat Flux = 14266973.33/ 542 =26322.8 W/ m2
Shell / Bundle dia =1.4 ( for the range of Heat Flux (25,000 to 40,000W/m2 )
Shell Diameter and Disengagement space
Shell Dia = 1.4 * 0.8428= 1.1799 m
Space for Disengagement = 1.1799 - 0.889=0.2909 m
For single phase flow
At T = 76 oC
ap =0.001 kg /m sec = 0.36kg / m hr
De = 0.01397 m
Res = DeGs/= 0.0198 168799.1748 / (0.001 *3600)=843.8
For Res = 843.8
Friction factor for shell side f = 0.0033= 0.4752 ft2/ ft2
No. of crosses, N + 1 = 11
Shell Side Pressure Drop
f Gs 2 Ds (N + 1)
Δ Ps = ½ ( )
2.54x 1011 De S s
Re.s = 843.8
f =0.0033
De = 0.0198m.
s = 0.002169
Gs =168799.1748 kg / hr m2
0.003*(168799.1748) 2 x0.889 x 11
Δ Ps = ( )
2 x 2.54 x 1011 x 0.0198 x 0.002169
= 38.31 kg/m2=0.054 psi
Pressure Drop (Tube Side)
Reynolds No.
Ret = DG/= 0.01484 254439/(0.00017 3600 ) = 6169.73
For Ret = 6169.7 Friction factor for tube side f = 0.4464 ft2/ft2
fG 2t Ln
Pt
2.54 1011 Dti st
0.4464 254439 9 2
2
=
2.54 1011 0.01484 0.87
= 160 kg / m2 =0.227psi
Pr = (4n / s)(v2/2g) =
=4*2/0.87*0.0005=0.00459 psi
=4*2/0.87*0.3518 =3.2357kg / m2
PT Pt Pr =163 kg /m2
= 0.231 psi
SPECIFICATION SHEET
REBOILER
Identification: Item Reboiler
No. Required = 1
Function: Boils the solution of Water & MEG
Operation: Continuous
Type: 1-2 Horizontal Kettle type
Shell side: Boiling
Heat Duty: 5136113.95 KJ/hr
Tube Side Tubes 0.01905 m dia. 16 BWG
Fluid handled: Saturated Steam No of Tubes=1068
Flow rate = 23484.72 Kg /hr Length = 9m
Pressure = 30397.5 N/m2 (0.3 Atm) No of Passes = 2
Temperature = 125oC Pitch=0.0238 m triangular
pressure drop =163 kg /m2=0.231psi
Shell Side: Shell: 0.889 m diam.
Fluid handled = MEG & DEG No of Pass =1
Flow rate = 53391.179 kg/hr Baffles spacing = 0.818 m .
Temperature = 76 to 90 oC Pressure drop = 38.31 kg/m2=0.054
Pressure = 30397.5 N/m2 (0.3 Atm) psi
Bundle Dia: 0.8428 m
Disengagement Space = 0.2909 m
Introduction
A condenser is a two-phase flow heat exchanger in which heat is generated from
the conversion of vapor into liquid (condensation) and the heat generated is removed
from the system by a coolant.
Condenser Types
Condensers may be classified into two main types:
Direct Contact
Those in which the coolant and condensing vapor are brought into
direct contact.
A vapor which is bubbled into a pool of liquid
A liquid which is sprayed into a vapor
A packed-column in which the liquid flows downwards as a film over
a packing material against the upward flow of vapor.
Indirect Contact
Those in which the coolant and condensate stream are separated by a
solid surface, usually a tube wall,
Condensers in which the streams are separated
Air-cooled
In the air-cooled type, condensation occurs inside tubes with
cooling provided by air blown or sucked across the tubes. Fins with
large surface areas are usually provided on the air side to compensate
for the low air-side heat transfer coefficients
Shell-and-tube
In shell-and-tube condensers, the condensation may occur inside
or outside the tubes. The orientation of the unit may be vertical
or horizontal
Coiled
Evaporator and condenser coils are used when the second fluid is
air because
of the low heat transfer coefficient on the air side
Selection of Indirect Contact
In direct contact , the two streams which are mixed together to remove the latent heat of
one stream, are to be separated again .So another equipment is required so the capital and
operating cost will be high so usually they are not recommended
Best choice is the indirect contact
Four Condenser Configuration are Possible
1) Horizontal with condensation is shell side and cooling medium in the tubes.
2) Horizontal with condensation in tube side cooling medium in shell side.
3) Vertical with condensation in the shell.
4) Vertical with condensation in the tubes.
Horizontal shell side and vertical tube side are the most commonly used types of
condensers.
Selection of Horizontal Condenser
When a condenser is employed on a distillation column, several specific
factors must be taken into consideration.
For the same height , horizontal condensers are more suited for the
Reflux to send back by gravity
Vertical Condenser has to be elevated considerably above the column, which in
many cases is very tall by itself. Maintenance and structural support may be costly , are
pump is installed for the reflux to send back which is again costly
Most suited choice is Horizontal Indirect contact type Condenser
In this process we have used the normal mechanism for heat transfer in
commercial condenser which film wise condensation.
Since vapor-liquid heat transfer changes usually occur at constant or really
constant pressure in industry, the vaporization or condensation of a single compared
normally occurs isothermally.
If a mixture of vapors instead of a pure vapor is condensed at constant pressure,
the change does not take place isothermally in most instances.
=(3.337*107)/(1*(113-77)) m =926944.44 lb / hr
LMTD = ((161.6-77)-(161.6-113)) / Ln ((161.6-77)/(161.6-113))
LMTD = 64.9 oF
LMTD = ((72-25)-(72-45)) / Ln ((72-25)/(72-45))
LMTD =36 oC
Caloric Temperature
Tc and tc Tc=T=161.6 oF =72 oC
tc=(t1+t2)/2 t1=77 oF = 25 oC t2 =113oF=45 oC
tc=(25+45)/2= 35 oC
Overall Design heat transfer Coefficient
Assumed Ud= 513 BTU/hr/sq.feet/F = 2914.4597W /m2 K
Estimated Area
=Q/(Ud*LMTD)
= 3.52127 *107 / ( 2914.459 * 36 )*1000 / 3600
A = 93.2 m2
Tube side Calculations
Tubes: ¾ in. dia . 14 BWG
Effective Length of the tubes L=5m
O.D of the tubes=0.75 in = 0.0625 ft = 0.01905 m
I.D of the tubes =0.584 in = 0.04867 ft = 0.01484 m
Surface area per linear ft of tube
=0.1963 ft2/ft= 2.3556 in2/in =0.0598m2/m
Surface area =at = 5*0.0598 = 0.292 m2
= 3.1408 ft2
No of tubes
nt=A/at =93.2/0.292 = 318.752
Density of the water ρ=60.2 lb/ft3=964 kg / m3
Volumetric Flow Rate of water
V=w/ ρ V= 420310.55 /(964*3600) V=0.121m3/sec
velocity of water
v = 4 ft/sec =1.2 m /sec
Total Cross Sectional Area
Ax=V/v Ax =0.121/1.2 Ax = 0.10083 m2
Cross Sectional Area of a tube x=0.546 in2 =0.00035226 m2
Tubes per Pass
tp=Ax/x tp = 0.10083 / 0.00035226 = 286.2
No of tube pass
np=nt / tp np= 318.752/286.2 = 1.1137
Corrected number of tube pass
np =1
shell and tube specification
Corrected value of number of tubes
Correct nt=361
Corresponding Shell Dia
D=21.25 in =0.53975 m
Corrected A
A=nt*a t= 361*0.292 A = 105.412 m2
Corrected Ud
Ud = Q/(Acor*LMTD)
Ud = 3.52127 *107 /(105.412*36)*1000/36000
= 2570.8 W/m2 K=452.96 BTU/ hr ft2 oF
Tube Cross sectional Area
ac=0.268 in2=0.0001729 m2
Total Cross sectional Area
Act=at*nt /np
At=0.0001729*361 / 1 At=0.0624 m2
Mass velocity
Gt=w/At Gt=420310.55 /0.0624 Gt=6733810.058 kg/ hr m2
Velocity
v=Gt/Rhoc
v=6733810.058 /(964*3600)
v=1.940 m/sec (assumption is true)
Caloric Temperature
tc= 35 oC
Inside tube heat transfer coefficient based on the outer diameter
hi= 1440 BTU/hr/sq.feet/F =8180.928 W/m2 K
hio = hi*id/OD
hio =1440*0.584 / 0.75
hio =1121.3 BTU/hr/sq.feet/F = 6370.32 W/m2 K
Shell Side Calculations
Shell Diameter= 21.25 in = 0.53975 m
Pitch =Pt=15/16 in =0.0238 m
B=ID=0.5397 m
Number of baffles = nb=5/0.5397 =9.26
Corrected baffle number
Nb =10 O.D of the tubes= 0.75 in = 0.0625 ft = 0.01905 m
I.D of the tubes = 0.584 in =0.04866 ft = 0.0148 m
Clearance
C=pt-OD=0.0238-0.01905 = 0.00475 m= 0.1875 in
B=L/nb*12
B= 5/10= 0.5 m =19.68 in
Shell Cross sectional Area
as=ID*c*B/(pt)
as = 0.53975 * 0.00475*0.625/(0.0238)
as=0.0657968 m2 =0.7083 ft2
Mass Flow rate =wf= 15128 kg/hr
Mass Velocity
Gs= wf / as;
Gs = 15128 /0.0657968=229919.996
Gs =229919.996 kg/m2hr
at= 0.292m2 np= 1
Condensing Load
Gl=w/(L*nt2/3)
Gl=15128/(5*3612/3)
Gl=59 kg / m hr =39 lb / hr ft
Dirt Factor
Rd= 0.0005/5.6812 = 0.000088 m2 K/W
Overall Clean Coefficient Uc=Ud/(1-Ud*Rd) Uc=Ud/(1-Ud*Rd)
Uc=452.96 /(1-452.96 *0.0005)
Uc=585 BTU/hr/sq.feet/F = 3323.5 W/m2 K
hio = 1121.3 BTU/hr/sq.feet/F= 6770.329 W/m2 K = 22933.186 KJ /m2 hr
ho=Uc*hio/(hio-Uc) ho= 1225.7 BTU/hr/sq.feet/F
ho =6963.44 W/m2 K
Wall Temperature
tw=tc+hio/(hio+ho)*(Tc-tc)
tw=35+24195.06 /(24195.06 +3852.72)*(72-35) tw = 66
tw = 54.3 oC= 129.7818 oF tf=(Tc+tw)/2
tf = (72 + 54.3)/2 = 63.15 oC= 145 .6909 oC
kf at tf' = 0.3860 Btu ft / hr /sq.ft / F
sg at tf= 0.9750
μf at tf=0.45 Cp ho=1800 Btu/hr ft2 oF
ho = 10119.24 W/m2 K
Over all Clean Coefficient
Uc=hio*ho/(hio+ho) =326 Btu / hr ft F =1852.073863 W/m2 K
Over all Desing Coefficient
Ud=Uc/(1+Uc*Rd) = 280 Btu ft / hr ft F= 1590.73 W/m2 K
=5723 KJ /hr m2
Pressure Drop (Shell Side)
At T = 72 oC
ap =0.00042 kg /m sec =1.512 kg / m hr
De = 0.01397 m
Res= DeGs/= 0.01397 229919.996 /0.00042 *3600
= 2124.32
For Res = 2124.32
Friction factor for shell side f = 0.4032 ft2/ ft2
No. of crosses, N + 1 = 11
Shell Side Pressure Drop
f Gs 2 Ds (N + 1)
Δ Ps = ½( )
2.54x 1011 De S s
Re.s = 2124.32 =0.0028
De = 0.01396m. s = 0.002169
Gs =229919.996 kg / hr m2
0.028*(229919.996) 2 x0.53975 x 11
Δ Ps =
2 x 2.54 x 1011 x 0.01396x 0.002169
= 571.kg/m2=0.812 psi
Pressure Drop (Tube Side)
Reynolds No.
Ret = DG/ = 0.01484 6733810/(0.00054 3600 ) = 50469.57
for Ret = 50469.57
Friction factor for tube side f = 0.02448
fG 2t Ln
Pt
2.54 1011 Dtit
0.02448 6733810 5 1
2
=
2.54 1011 0.01484 0.964
= 1527 kg / m2 =2.17psi
Pr = (4n / s)(v2/2g) =
=4*1/0.964*0.21=0.87 psi
=4*1/0.964*147.79 =613.2kg / m2
PT Pt Pr =2140 kg /m2
= 3.041 psi
SPECIFICATION SHEET
CONDENSER
Identification: Item condenser
No. Required = 1
Function: Condense vapors by removing the latent heat of vaporization
Operation: Continuous
Type: 1-2 Horizontal Condenser
Shell side: condensation
Heat Duty = 35212712 KJ/ hr
Tube Side: Tubes :0.01905 m diam. 14 BWG
Fluid handled: Cold water No of Tubes=361
Flow rate = 420310.55 Kg /hr Length of Each Tube = 5 m
Temperature = 25oC to 45 oC No of passes =1
Pitch =0.0238 m triangular
pressure drop =2140 kg /m2 =3.041
psi
Shell Side: Shell: 0.53975 m diam.
Fluid handled = Water & MEG 1 pass
Flow rate = 15128 kg/hr Baffles spacing = 0.5 m
Temperature =72 oC Pressure drop =571.kg/m2 =0.812psi
Bottom product from Light end column contains Monoethylene glycol and higher glycols
DEG and TEG and also a negligible amount of water .MEG is separated from higher
glycols and impurities. For polyester grade MEG 99.9% pure MEG is required. For this
purpose we selected tray column. Selection criteria has been discussed with light end
column
Material Balance
Feed
Wt Wt Wt Moles Moles Moles
(kg/hr) (kg/sec) (%) (Kmole/hr) (kmole/sec) (%)
Distillate
Wt Wt Wt Moles Moles Moles
(kg/hr) (kg/sec) (%) (Kgmoles/hr) (kgmoles/sec (%)
Bottom
MEG DEG
F 0.125412456 kmoles/sec Xf 0.906662 Xf 0.0876
D 0.1131962 kmoles/sec XD 0.9983 XD 0.0002
B 0.012216256 kmoles/sec XB 0.05678 XB 0.8977
NUMBER OF PLATES
Fenske equation
X X
ln LK HK
X HK d X LK b =6.71 plates with reboiler
Nm
ln LK
q =1+CpL*(Tb-Tf) / λ=1+0.44*(157-95) / 214.77= 1.157
By hit and Trial Θ = 1.0675 Rm = 0.2538
R = 1.2xRm R = 0.3046 N= 20
Column efficiency
using O’connell method
Molar avg. liquid viscosity =0.3899mNs/m2
Average relative volatility of the light key =5.179
So, Actual no. of plates = 20 / 0.42 = 47 plates
NB=5 / 0.42= 12
ND=15 / 0 .42=35
Maximum liquid flow rate in rectifying section = Ln = D Rm= 0.304598x0.113196
= 0.0345479kmole/sec
= 0.0345479+0.1131962
= 0.147744 kmole/sec
0.034479+(1.157)(0.1254)=0.17958 kmole/sec
Weir length
Ad / Ac = 1.563 / 10.42 = 0.15 Lw / dc = 0.81 Lw =2.95 m
Check Weeping
Maximum liquid rate = 0.17958 × 105.62
= 18.98 kg/sec
K2 = 30.7
0.1674 105.62
Û h = Volumetric Flow Rate / Hole Area Uˆ h =37.56 m/s
0.64159 0.72
Ah Ah
For plate thickness/hole dia = 5/5 = 1 and 0.1
Ap Aa
lo = 0.84
2
Uˆ h V 2
37.56 0.64519
hd 51 hd 51
lo L 0.84 944.23
= 69.67mm liquid
Residual Head
12.5 103
hr 13.24 mm liquid
l
ht = 161.6 mm liquid
hb = 0.249 liq. m
0.249 < ½ (Tray spacing + weir height) So tray spacing is acceptable
Check Entrainment
Uv = Maximum Volumetric Flow Rate of vapors/Net Area
FLV = 0.028
No of Holes
0.7294
No. of holes
1.964 105
= 37138
Height of Column:
Number Of plates = 47
Height of column = (#of plates × space between each plate) + (space for
So
SPECIFICATION SHEET
Identification:
Item Distillation column
Item No.
No. required 1
Tray type Sieve tray
Function: Recovery of MEG
Operation: Continuous
Material handled
Feed Top Bottom
Quantity F=451.48 D=407.42 B =43.98
kg.moles/hr kg.moles/hr
kg.moles/hr
Design Data
Hole area/Active area = 0.10
No. of trays = 47 Weir length = 2.95 m
Height of column = 24.5 m Reflux ratio = 0.3046
Diameter of column = 3.64 m Hole size = 5 mm
Tray spacing = 0.500 m Down comer area = 1.563 m2
Tray thickness = 5 mm Hole area = 0.7294 m2
Flooding = 80 % Active area = 7.294 m2
6.7.4:REBOILER DESIGN
Calculate Heat Duty – Q
Shell Side: (cold) Tube Side: (hot)
MEG & DEG Steam
Inlet Temp of cold fluid = t1 = 90 oC
Outlet Temp of cold fluid= t2 = 157 oC
ta = 123.5 oC
Inlet Temp of hot fluid = T1 = Ta= 212oC
Outlet Temp of hot fluid = T2 =Ta = 212oC
Preheat = Q1=mCp∆t = 9238925 KJ/hr
Latent Heat = Q2 = m λ = 38649042 kJ/ hr
Total Heat =Q = 47887967 kJ/ hr
LMTD
LMTD 1= ((Ta-t1)-(Ta-t2)) / Ln ((Ta-t1)/(Ta-t2))= 84.0982 oC
LMTD 2 = 54 oC
Weighted LMTD = Qt / ((Q1/LMTD1) + (Q2/LMTD2))
Weighted LMTD = 58.9341 oC
DesignOverall heat transfer Coefficient
Assumed Ud=1.3639e+003 W /m2 K
Ud=1.3639e+003 W /m2 K Ud1=374.1021 W /m2 K
Ud2=481.5266 W /m2 K
Estimated Area
A=Q/(Ud*LMTD)*1000/3600 Ac1=Q1*1000/(Ud1*lmtd1*3600)
Ac1 = 81.57 m2
Ac2=Q2*1000/(Ud2*lmtd2*3600)
A=Qtotal/(Ud*wlmtd)*1000/3600
A = 165.49m2
5.6.2:SPECIFICATION SHEET
REBOILER
Operation: Continuous
Type: 1-2 Horizontal Kettle type
Shell side Boiling
Heat Duty = 47887967 KJ/hr
Tube Side: Tubes 0.01905 m diam. 14 BWG
Fluid handled Saturated Steam No of Tubes =506 tubes
Flow rate = 23484.72 Kg /hr Length of Each Tube = 5m
Temperature = 212oC No of Passes =2
Pitch=0.0238 m triangular pitch
5.6.3:CONDENSER DESIGN
Complete designing and selection have been discussed in detail in designing of condenser
for light end column. The design procedure is same. so here is given the specification
sheet for condenser.
SPECIFICATION SHEET
CONDENSER
Operation: Continuous
Type: 1-2 Horizontal Condenser
Shell side condensation
Heat Duty = 30292120 KJ/ hr
Tube Side: Tubes:0.01905 m diam. 14 BWG
Fluid handled :Cold water No of Tubes =224
Flow rate = 120695 Kg /hr Length of Each Tube =5m
Pressure = 101325 N/m2 (1 Atm) No of Passes =2
Temperature = 25oC to 85 oC Pitch :0.0238 m triangular
pressure drop = 5.078 psi
Shell Side: Shell: 0.43815 m . diam.
Fluid handled = MEG & DEG No of Pass = 1
Flow rate = 8475.99 lb/hr Baffles spacing = 0.4167 m .
Pressure = 14.7 psia Pressure drop =0.18 psi
Temperature 154 oC =128 kg/m2
Utilities: Cold water
Ud assumed = 5723 KJ /hr m2 K Ud calculated =5723 KJ /hr m2
Uc calculated = 1522.56 W/m2 K Rd = 0.000088 m2 K/W
In next chapter we will discuss about instrumentation and control.
CHAPTER-6
INSTRUMENTATION AND
CONTROL
6.1:OBJECTIVES
The primary objective of the designer when specifying instrumentation and control
schemes are:
6.1.1: Safe Plant Operation
This is used to control the temperature of outlet and inlet streams in heat
exchanger, condenser, heater and column, Most temperature measurement in the
industrial processes is made by means of thermocouples to facilitate bringing
measurements to a controlised location. Usually where high measurement accuracy
is required resistance thermometer are used. For the control of temperature,
pneumatic feedback control is used. Outlet temperature of the system is measured,
with the help of thermocouple. This temperature is fed to a comparator. The
difference between the set point and the input signal is called error, is fed to the
controller, which controls diaphragm valve to adjust the flow of hot or cooling
medium.
The temperature and pressure are valuable indications of material state and
composition. In fact, these two measurements considered together are the evaluating
devices of industrial materials. In compressed gas systems etc., pressure is of
primary importance; pumping equipments are furnished with pressure measuring
devices. Pressure measurement becomes an indication of the increase or decrease.
Most pressure measurement in the industry is made through elastic devices, either
directly connected for local use or transmission type to a controlised location.
Most flow measurements in the industry are done by variable head devices.
To lesser extent variable area and displacement types are used.
6.5:CONTROL LOOPS
For instrumentation and control of different sections and equipments of plants,
following control loops are most often used.
1. Feed back control loop
2. Feed forward control loop
3. Ratio control loop
4. Auctioneering control loop
5. Split range control loop
6. Cascade control loop
Here is given a short outline of these control schemes, so that to justify out
selection of a control loop for specified equipment.
6.5.1:Feed Back Control Loop
A method of control in which a measured value of a process variable is
compared with the desired value of the process variable and any necessary action is
taken. Feed back control is considered as the basic control loops system. Its
disadvantage lies in its operational procedure. For example if a certain quantity is
entering in a process, then a monitor will be there at the process to note its value.
Any changes from the set point will be sent to the final control element through the
controller so that to adjust the incoming quantity according to desired value (set
point). But in fact changes have already occurred and only corrective action can be
taken while using feed back control system.
6.5.2:Feed Forward Control Loop
A method of control in which the value of disturbance is measured than action is
taken to prevent the disturbance by changing the value of a process variable. This is
a control method designed to prevent errors from occurring in a process variable.
This control system is better than feed back control because it anticipates the change
in the process variable before it enters the process and takes the preventive action.
While in feed back control system action is taken after the change has occurred.
6.5.3:Ratio Control
A control loop in which, the controlling element maintains a predetermined ratio of
one variable to another. Usually this control loop is attached to such as system where
two different systems enter a vessel for reaction that may be of any kind. To maintain
the stoichiometric quantities of different streams, this loop is used so that to ensure
proper process going on in the process vessel.
6.5.4:Auctioneering Control Loop
This kind of control loop is normally used for a huge vessel where readings of a
single variable may be different at different locations. This type of control loop
ensures safe operation because it employs all the readings of different locations
simultaneously, and compares them with the set point, then the controller sends
appropriate signal to final control element.
6.5.5:Split Range Loop
In this loop controller is preset with different values corresponding to different
actions to be taken at different conditions. The advantage of this loop is to maintain
the proper conditions and avoid abnormalities at very differential levels.
6.5.6:Cascade Control Loop
This is a control in which two or more control loops are arranged so that the output of
one controlling element adjusts the set point of another controlling element. This
control loop is used where proper and quick control is difficult by simple feed forward
or feed backward control. Normally first loop is a feed back control loop.
6.6.2.1:Transfer Function
For this example, we'll assume that the process is second-order with the following
transfer function:
Gp = 1/(s2 + 10s + 20)
For simplicity's sake we are assuming no gains in the instrumentation. The goal for this
example is to control the distillation column temperature to 100 °F, minimize response
time, and minimize overshoot.
The open loop response tells us what we need to work on with our controller. As most
steam control valves fail to a closed position, without any control, the temperature will
only rise 5 °F.
6.6.2.3:P Control
Proportional control is used to improve the rise time. It does bring the temperature close
to the set point, but there is still a steady-state error. The following graph has a
proportional gain, Kc = 300.
6.6.2.4:PI Control
Proportional-derivative control improves the overshoot, but fails to improve the steady-
state error. The following graph displays the response of PD control with Kc = 300 and D
= 0.033.
6.6.2.5:PID Control
Tuning proves the best response with Kc = 350 and I = 1.17 and D = 0.143. As shown
below, the closed loop response for the combined proportional-integral-derivative control
gives the quickest response with no overshoot and converges to the desired 100 °F.
Thus, we can see that PID control is the optimal technique for temperature control of a
distillation column.
One thing to note is that in none of them is a valve simply placed on the vapour line. This
would lead to the use of a large expensive control valve. Instead the pressure is controlled
indirectly involving the use of the condenser and/or reflux drum
CHAPTER-7
HAZOP STUDY
A HAZOP survey is one of the most common and widely accepted methods of
systematic qualitative hazard analysis. It is used for both new or existing facilities and
can be applied to a whole plant, a production unit, or a piece of equipment It uses as its
database the usual sort of plant and process information and relies on the judgment of
engineering and safety experts in the areas with which they are most familiar. The end
result is, therefore reliable in terms of engineering and operational expectations, but it is
not quantitative and may not consider the consequences of complex sequences of human
errors.
The objectives of a HAZOP study can be summarized as follows:
1) To identify (areas of the design that may possess a significant hazard
potential.
2) To identify and study features of the design that influence the probability of a
hazardous incident occurring.
6) To provide a mechanism for feedback to the client of the study team's detailed
comments.
MATERIAL SAFETY DATA
SHEET
8.1:CHEMICAL IDENTIFICATION
8.2:COMPOSITION
Eye Immediately flush eyes with water for at least 15 minutes. Hold
Contact: eyes open while flushing out with water. Seek medical attention
immediately.
Skin Immediately remove contaminated clothing and shoes. Flush
Contact: skin with water for at least 15 minutes. Use soap if available or
follow by washing with soap and water. Do not reuse
contaminated clothing without laundering. If irritation persists,
seek medical attention.
Inhalation: Remove victim to fresh air. If breathing is difficult, give
oxygen. If not breathing, administer artificial respiration. Seek
medical attention immediately.
Ingestion: Induce vomiting immediately as directed by medical personnel.
Seek medical attention immediately.
Unusual Fire and Do not use direct stream of water to fight fire.
Explosion Information: Containers can build up pressure if subjected to
heat of the fire and may explode. Flashback hazard
– vapors are heavier than air and can collect in low
areas forming an explosive MEG and air mixture.
Environmental Note: Prevent product from getting into sewers or surface
waters.
Isolate the hazard area and deny entry to nonessential personnel. Emergency responders
and/or clean-up personnel should wear appropriate protective clothing and equipment
when responding. Remove all ignition sources. Prevent from entering sewers or surface
waters. Collect liquid in containers and seal shut. Absorb remaining material with a
noncombustible absorbent such as earth, sand, or vermiculite and collect for disposal.
Keep away from heat, sparks, and open flames. Keep containers tightly closed. Store
away from incompatible materials in a cool, dry place. Use adequate explosion-proof
ventilation to prevent accumulation of static charge. When pouring or transferring
materials, containers must be bonded and grounded.
Do not weld, heat, or drill on or near full or empty containers. Empty containers can
contain explosive vapors.
Do not breathe vapors or mist. Minimize skin contact. Wash with soap and water before
eating, drinking, smoking, or using toilet facilities. Launder contaminated clothing
before reuse. Properly dispose of contaminated leather articles, including shoes that
cannot be decontaminated.
8.11:TOXICOLOGICAL INFORMATION
8.12:ECOLOGICAL INFORMATION
8.13:DISPOSAL INFORMATION
Place in a city, state, or federally permitted disposal facility. Handle in accordance with
all applicable regulations.
8.14:TRANSPORTATION INFORMATION
COST ESTIMATION
Before the plant to be operated, specified money must be supplied to purchase and install
the equipment. The capital needed to supply the necessary plant facilities is called fixed
capital investment while that for the operation of the plant is called the working principal
and sum of two capitals is called total capital investment.
9.1:COST INDEX
Is and index value for a given point in time showing the cost at that time
relative to certain base time. if the cost at any time in past is know , the equivalent cost at
the present time can be determined by multiplying the original cost by the ratio of present
of index value to the index value applicable when the original cost was obtained. There
are many types of cost indexes and the most commune type are marshall and swift all
industry and process industry equipment indexes and engineering News Record
construction index.
Present cost =
Original cost *cost index value at present time/cost index value at past.
Fixed costs
Maintenance cost = ₤3600684
Operating labor = ₤5*10^5
Laboratory cost =₤1*10^5
Supervision =₤1*10^5
Plant overheads = ₤2.5*10^5
Capital charges = ₤5401026
Local taxes = ₤720137
Insurance = ₤360068
10) Kuppan, T., “Heat Exchanger Design Hand Book, Marcel Dekker
Inc., New York, 2000.
15) www.meglobal.com
16) www.nist.com
17) www.uspto.gov
REFERENCE GRAPHS AND TABLES
9.4:OPERATING COSTS
Operating costs are dependent on the amount of product produced, so these costs are
variable and changes with production rate.
Raw materials and product cost
= 0.30 ton/kg
Maintenance cost = 0.10 * fixed capital cost
= 0.10 * 3934522
= ₤393452.2
Miscellaneous = 0.10 * maintenance cost
= 0.10 * 393452.2
=₤39345
Operating labor = ₤5*10^5
Laboratory cost = 0.20* labor cost=₤1*10^5
Supervision=0.20* labor cost=₤1*10^5
Plant overheads = 0.50 * labor cost =₤2.5*10^5
Capital charges = 0.15 * fixed capital cost =590178
Local taxes = 0.02 * fixed capital cost =78690
Insurance = 0.01 * fixed capital cost=39345
Regulations and license fees = 2 % of sales price
9.5:PRODUCTION COST
PRODUCTION COST IS CALCULATED BY
Variable cost
Raw materials = 0.30 * 200,000,000
= ₤60,000,000
Miscellaneous = ₤39345
Shipping and packing negligible
Sub total = A = ₤60039345
Fixed costs
Maintenance cost = ₤393452.2
Operating labor = ₤5*10^5
Laboratory cost =₤1*10^5
Supervision =₤1*10^5
Plant overheads = ₤2.5*10^5
Capital charges = ₤590178
Local taxes = ₤78690
Insurance = ₤39345