FYP PROJECT Dheeraj & Kashyap
FYP PROJECT Dheeraj & Kashyap
FYP PROJECT Dheeraj & Kashyap
Project
On
2-METHOXY 2-METHYLE HEPTANE
of
B.Tech.
in
Chemical Engineering
By
Kashyap Parmar (U20CH048)
Dheeraj Salvi (U20CH036)
Certificate
This is to certify that the B.Tech. seminar, entitled “Production of 2-Methoxy 2-Methyl
Heptane” submitted by Kashyap Parmar (U20CH048) & Dheeraj salvi (U20CH036) in partial
fulfillment for the award of the degree of “B.Tech. in Chemical Engineering” of Sardar
Vallabhbhai National Institute of Technology, Surat, is a record of his own work carried out under
my supervision and guidance.
Supervisor:
Dr.G.C.JADEJA
(Assistant Professor)
Date: 08/12/2023
Examiners:
1. Dr. Meghal.A.desai
(Associate Professor)
2. Dr.Smita Gupta
(Associate Professor)
3. Dipesh kachhadiya
(Teaching assistant)
Date: 08/12/2023
ABSTRACT
This project report is an attempt to understand the basic principles and technicalities of
manufacturing of a product at an industrial scale. Here we have selected 2-Methoxy- 2-
Methyl Heptane as our product. This report describes about introduction of 2-
Methoxy-2-Methyl Heptane, its various properties including physical, chemical as well
as thermodynamic properties. Also, various methods of industrial production practiced
today are introduced and a comparison between them is made and the best and optimized
method of production is chosen according to the feasibility and context of the project.
For the method chosen we have presented an attempt to the material balance and energy
balance for the entire process as well as detailed balance at each stage of the process.
We have done equipment’s design as well as plant location and layout.
Thermodynamics and kinetics of the reactions are also represented. We also made a
cost estimation of 2-Methoxy-2-Methyl Heptane and safety issues and ETP of plant.
Thus, this project represents an idea of industrial production of 2-Methoxy-2-
Methyl Heptane.
(U20CH048) ( U20CH036)
Contents
1 Introduction
1.1 Properties.....................................................................................................................
3 Process Description
3.1 Chemistry Behind the Process...............................................................................
3.2 Process Steps................................................................................................................
3.3 Vapour-Liquid Equilibrium Graphs....................................................................
4 Material Balance
4.1 Column-3......................................................................................................................
4.2 Column-2......................................................................................................................
4.3 Column-1......................................................................................................................
4.4 Reactor
5 Energy Balance
5.1 Reactor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
5.2 Column-1 .. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
5.3 Column-2 .. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
5.4 Column-3 .. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
6 Thermodynamics and Kinetics of the Process
6.1 Thermodynamics of the Process.................................................................................
6.2 Thermodynamics of Desired Reaction...................................................................
6.2.1 Calculating Change in Enthalpy, ∆H.................................................
6.2.2 Calculating Change in Entropy, ∆S...................................................
6.2.3 Calculating Change in Gibbs Free Energy, ∆G......................................
6.3 Kinetics of the Process.................................................................................................
6.3.1 Finding the Volume of Reactor..................................................................
9 Plant Layout
9.1 Factors in Planning Layout...........................................................................
9.2 Facilities to be Included in a Layout................................................................
9.3 Methods of Layout Planning................................................................................
9.4 Unit Area Diagram of the Plant Layout.......................................................
Cost Estimation
9.5 Purchased Equipment Cost..............................................................................
9.6 Estimation of Total Capital Investment...........................................................
9.6.1 Direct Costs...........................................................................................
9.6.2 Indirect Costs........................................................................................
9.7 Estimation of Manufacturing Cost...................................................................
9.7.1 Direct Production Cost.........................................................................
9.7.2 Fixed Charges........................................................................................
9.8 General Expenses..................................................................................................
9.9 Gross Earnings and Net Profit.........................................................................
9.10 Rate of Return, Payback Period and Turnover Ratio......................................
10 Bibliography
Chapter 1
Introduction
2-Methoxy-2-Methylheptane (MMH) is a flammable and irritation causing
substance with a molecular weight 144.25 g/mol.
The trade name of the following chemical is Daikon Ether. The chemical 2-methoxy- 2-
methylheptane (MMH) has been proposed as a gasoline additive to replace methyl tert-
butyl ether (MTBE) in order to avoid groundwater contamination.
MMH Structure
Current advancements are going on in this field with different methods like reactive dis-
tillation being tested to get better and better quality MMH.
Currently not a lot of companies are involved in producing MMH. Some prominent in-
dustries like BASF, Harmacore Co. Ltd, Alfa Aesar, etc.
Properties
g
Molecular Weight 144.25
mol
Density 875.8 Kg
m3
◦
Melting Point -77.4 C
Boiling Point 179◦ C
Vapour Pressure 9.48 mmHg
Water Solubility 3.49 × 10−3 mol
l
Hazard Flammable
Flash Point 37.1◦ C
kJ
Heat Capacity (STP) 1.87522
kg·K
Viscosity (STP) 0.0004189 Pa· s
Thermal Conductivity 119 mW
m·K
Property Table
Chapter 2
• The major key players which drive the demand and supply of MMH are Alfa
Aesar, Harmacore Co., Hangzhou J&H Chemicals, BASF.
• The global market is classified on the basis of product type which is namely
Phar- maceutical Grade and Chemical Grade.
• If seen on the basis of use or application MMH is classified into Analytical Reagents,
Pharmaceutical Intermediates, Others.
Applications
Demand Data
mg
Product Type Daily Systematic Exposure
bw
kgday
Cosmetic Products
Body Lotion 0.0611
Face Cream 0.0120
Hand Cream 0.0169
Fine Fragrances 0.0586
Deodorant 0.0117
Shampoo 0.0008
Conditioner 0.0003
Shower Gel 0.0015
Hand Wash Soap 0.0016
Hair Styling Products 0.0031
Sub Total 0.1676
House Hold Products
Laundry Liquid 0.00170
Fabric Softener 0.0007
Dish Washing Liquid 0.0001
All Purpose Cleaner 0.0011
Sub Total 0.0036
Aerosol Exposure
Hair Spray 0.0016
Total 0.1728
Demand Data
Production Trend
Production Trend
Consumption Trend
Chapter 3
Process Description
The chemistry involves the liquid phase reversible reaction of methanol with 2- methyl-
1-heptene (MH) to form MMH. However, methanol and MH also undergo an un-
desirable reaction to form dimethyl ether (DME) and 2-methyl-2-heptanol (MHOH).
The activation energies of the two competing reactions favour high reactor temper-
ature, which is limited by catalyst activity to 400 K. Yield of the desired product is in-
creased by either increasing reactor size or keeping the concentration of methanol
in the reactor low by using a large excess of MH in the reactor. The latter strategy
results in a large recycle flow rate of MH, which increases separation costs. Thus, this
process demonstrates the classic design trade-off between reactor costs and separation
costs.
• Reactions Involved
1. Desired Reaction
CH3OH + C8H16 − → C9H20O R1 = k1F · xMeOH · xMH − k1R · x MMH
2. Undesired Reaction
2 CH3OH + C8H16 − → C2H6O + C8H18O R2 = k2 · (xMeOH)2
Different R1 R2 Units
Components Forward Backward
kmol
k (Overall Reaction Rate) 6.7 × 107 2.1 × 10−6 1.3 × 109 s·kgcat
kJ
E (Activation Energy) 90000 900 105900 kmol
Concentration xMeOH · xMH x MMH x2MeOH mole fraction
Kinetics Table
Process Steps
MH Recycle
DME
C1 - DME Column
MMH
MH Feed
C2 - MH Recovery Column
CSTR
MeOH Feed
C3 - MMH Column
MHOH
1. Reactor:
The reactor is a 12 m3 continuous stirred tank reactor (CSTR) operating at 400
kmol
K and 15 atm. A fresh methanol stream (MeOH) is fed at 49.98 h
. Fresh
kmol
2-methyl-1-heptene (MH) is fed to the reactor at 49.43 h
and combined with a
kmol
recycle stream from a downstream distillation column to give a total of 130.43
h
of MH fed to the reactor. Concentrations of components in the reactor are 0.40
mol % dimethyl ether (DME), 2.47 mol % MeOH, 58.80 mol % MH, 37.30 mol %
MMH, and 0.40 mol % MHOH. The low methanol concentration is needed to
kmol
achieve the desired 99% yield of MMH. The reactor output is 131.25 . The per-
h
pass conversion of methanol is 93.9%, and the per-pass conversion of MH is
39.2%.
2. Column 1:
The reactor effluent is fed on stage 8 of a 12-stage distillation column. The column
operates at 10 atm, which gives a reflux-drum temperature of 318 K and permits
the use of cooling water in the condenser. The distillate product is a small stream
of 0.52kmol
h
of 99.9 mol % DME. The separation between DME and methanol is
easy, so the reflux ratio is small (RR 0.8) and few stages are required. The base
temperature is 509 K. Most of the reboiler heat goes to increase the sensible heat
of the feed from 400 K up to the base temperature of 509 K.
3. Column 2:
The bottoms from the first column is fed to column C2 separates the product MMH
kmol
from the recycle MH. The distillate is 81
h
with a composition of 4.1 mol %
MeOH, and 95.82 mol % MH. The bottoms is a mixture of the product MMH and
kmol
the impurity MHOH. The bottoms stream flow rate is 49.98
h
and composition
of the constituents is 98 mol% MMH and 1.04 mol% MHOH. The column has 42
stages and is fed on stage 23. It operates under vacuum conditions (0.4 atm) with
a reflux-drum temperature of 358 K. A reflux ratio of 2.198 is required to achieve
the specified separation. The base temperature is 410 K, so medium-pressure
steam can be used in the reboiler (11 bar, 457 K).
4. Column 3:
The third column has 22 stages and is fed on stage 14. The distillate is the MMH
kmol kmol
product 48.944 h
with a purity of 98.9 mol % MMH. The bottoms is 0.4953 h
of
byproduct MHOH with a purity of 99.8 mol %. The column operates under vacuum
(0.1 atm). The reflux drum temperature is 352 K, and the base temperature is
426 K using medium-pressure steam. The reflux ratio is 0.427.
Chapter 4
Material Balance
G (MH Recycle)
E (DME)
A I (MMH)
C1 - DME Column
(MH Feed)
D
C2 - MH Recovery Column
CSTR
C (MeOH Feed)
C3 - MMH Column
F
J (MHOH)
Material Balance Equation:
(Materialin)−(Materialout)+(Materialgenerated)−(Materialdisappearance) = Accumalation
kmol
Basis of Calculation: 48.94 based on amount of MMH i.e stream I
h
Given/Assumptions
Calculations
⬦ MeOH Balance
MeOHout - MeOHin = -ϵ1 - 2ϵ2
kmol
MeOHin = ϵ1 + 2ϵ2 = 49.98 h
kmol
C = 49.98 h
Column-3
Column-3
Column-2
Column-1
Column-1
Stream Species Fraction
D MeOH 0.0247
MMH 0.373
MH 0.588
DME 0.004
MHOH 0.004
E DME 0.999
MeOH 0.001
F MH 0.589
MMH 0.374
MeOH 0.0254
DME 0.0008
MHOH 0.004
∗ As F - G = 49.98
∗ Now D = E + F
kmol
D = 131.25 h
∗ At Mixing
Point, A + G =
B
49.43 + 81 = B
kmol
B = 130.43 h
Reactor
Reactor
Stream Species Fraction
D MeOH 0.0247
MMH 0.373
MH 0.588
DME 0.004
MHOH 0.004
B MH 0.974
MeOH 0.0254
C MeOH 1.0
Energy Balance
Basic Energy Balance Equation:
(Energyin) − (Energyout) + (Energygenerated) − (Energylost) = Accumalation
Reactor
Reactor
Every data at the reference temperature 298 K.
∗ Stream C:
393 K → 298 K
kJ
Species Flow Rate ( kmol ) Specific Heat ( ) ∆ T (K) ∆H( kJ
)
h kmol·K s
Stream C Energy
∗ Stream B:
353 K → 298 K
kJ
Species Flow Rate ( kmol ) Specific Heat ( ) ∆T (K) ∆H ( kJ
)
h kmol·K s
Stream B Energy
kJ
Species ∆Hf ( )
kmol
MeOH -7200677.4
MH -96990.3
MMH -9790.6
DME -183826.58
MHOH -0.001237
⬦ Desired Reaction
Σ
kJ
∆Hr1 = niHi = 287877.1
kmol · · · with respect to MeOH
(MeOH)B+(MeOH)c−(MeOH)D (MeOH)B+
Now XMeOH = (MeOH)c
kJ
∆Hr1 = (0.918)×(287877.1)×(52.3) = 3907.84 s
⬦ Undesired Reaction
Σ
kJ
∆Hr2 = niHi = 157259.26 · · · with respect to MeOH
kmol
kJ
∆Hr2 = (0.01)×(157259.26)×(52.3) = 45.8673
s
kJ
∆Hr = ∆Hr1 + ∆Hr2 = 3953.7073 s
∗ Stream D
298 K → 400 K
kJ
Species Flow Rate ( kmol ) Specific Heat ( ) ∆T (K) ∆H ( kJ
)
h kmol·K s
Stream D Energy
Column-1
Stream E Energy
⬦ Stream F (Bottom Product)
298 K → 509 K
kJ
Species Flow Rate ( kmol ) Specific Heat ( ) ∆T (K) ∆H ( kJ
)
h kmol·K s
Stream F Energy
∗ Balance on Condenser
Where L is the flow rate of reflux stream, D the flow rate of distillate and V the
flow rate of vapour coming into the condenser.
∗ Balance on Reboiler
Feed 1050.29 -
Distillate - 0.18099
Condenser - 8.4514
Reboiler 1281.97 -
Residue - 2323.484
Total 2332.207 2332.116
kJ
Energy Loss = Energy In - Energy Out = 0.091 s
Column-2
⬦ Stream G (Recycle)
298 K → 358 K
kJ
Species Flow Rate ( kmol ) Specific Heat ( ) ∆T (K) ∆H ( kJ
)
h kmol·K s
Stream G Energy 42
⬦ Stream H (Bottom Product)
298 K → 410 K
kJ
Species Flow Rate ( kmol ) Specific Heat ( ) ∆T (K) ∆H ( kJ
)
h kmol·K s
Stream H Energy
∗ Balance on Condenser
Where L is the flow rate of reflux stream, D the flow rate of distillate and V the
flow rate of vapour coming into the condenser.
∗ Balance on Reboiler
Feed 2323.484 -
Distillate - 329.935
Condenser - 3368.59
Reboiler 2282.661 -
Residue - 724.175
Total 4606.145 4422.7
kJ
Energy Loss = Energy In - Energy Out = 183.445 s
Column-3
Stream I Energy
⬦ Stream J (Bottom Product)
298 K → 426 K
kJ
Species Flow Rate ( kmol ) Specific Heat ( ) ∆T (K) ∆H ( kJ
)
h kmol·K s
Stream J Energy
∗ Balance on Condenser
Where L is the flow rate of reflux stream, D the flow rate of distillate and V the
flow rate of vapour coming into the condenser.
∗ Balance on Reboiler
Feed 843.317 -
Distillate - 233.354
Condenser - 962.584
Reboiler 221.692 -
Residue - 6.0868
Total 1065.01 1202.0248
kJ
Energy Loss = Energy In - Energy Out = -137.015 s
Thermodynamics of the reaction tells us whether the proposed route or reaction is pos-
sible or not.
It also tells us if the reaction is possible and is reversible what can be done in order to
favour a forward reaction (Le Chatelier’s Principle)
Now there are three important thermodynamic properties;
• Change in Enthalpy, ∆H
• Change in Entropy, ∆S
The change in Enthalpy, ∆H tells us about the nature of reaction i.e whether it is
exothermic or endothermic (∆H < 0 or ∆H > 0) which becomes very important
in order to guide the reaction in correct direction.
The change in Entropy, ∆S tells us regarding spontaneity as well as reversibility
of a reaction. For a completely reversible process ∆S = 0.
The change in Gibbs Free Energy, ∆G tells us whether reaction is possible or
not at given conditions or not. If ∆G < 0 means the reaction is possible.
Thermodynamics of Desired Reaction
where;
∆HP = Enthalpy Change of Product w.r.t ref temperature
∆HR = Enthalpy Change of Reactant w.r.t ref temperature
In order to calculate ∆H we need the values of change of enthalpy of products and
reactants which we can get from table below;
Component kJ
∆H mol
MeOH -228.3
MH -108
MMH -20
The properties in the above table are at reaction temperature i.e 400 K w.r.t reference
temperature i.e 298.15 K
∆H = ∆HMMH − ∆HMeOH − ∆HMH
kJ
∆H = -20 + 108 + 228.3 = 316.3 mo
l
As we see ∆H is positive as a result our reaction is endothermic in nature and ∴
we have provided a steam jacket for our CSTR so that the reaction can proceed in
forward direction.
Calculating Change in Entropy, ∆S
where;
∆SP = Entropy Change of Product w.r.t ref temperature
∆SR = Entropy Change of Reactant w.r.t ref temperature
In order to calculate ∆S we need the values of change of entropy of products and reactants
which we can get from table below;
Component kJ
∆S molK
MeOH -0.215
MH -0.66
MMH 0.025
The properties in the above table are at reaction temperature i.e 400 K w.r.t reference
temperature i.e 298.15 K
∆S = ∆SMMH − ∆SMeOH − ∆SMH
kJ
∆S = -0.025 + 0.66 + 0.215 = 0.85mol
K
where;
∆GP = Gibbs Free Energy Change of Product w.r.t ref temperature
∆GR = Gibbs Free Energy Change of Reactant w.r.t ref temperature
In order to calculate ∆G we need the values of change of gibbs free energy of products
and reactants which we can get from table below;
Component kJ
∆G mol
MeOH -143
MH 156
MMH -10
The properties in the above table are at reaction temperature i.e 400 K w.r.t reference
temperature i.e 298.15 K
∆G = ∆GMMH − ∆GMeOH − ∆GMH
kJ
∆G = -10 + 143 - 156 = -23 mo
l
Now we can also find ∆G by using the formula;
∆G = ∆H − T × ∆S
l
Either ways we are getting the same values of ∆G which is negative and that means
that the reaction is feasible thermodynamically.
The chemistry involves the liquid phase reversible reaction of methanol with 2- methyl-
1-heptene (MH) to form MMH. However, methanol and MH also undergo an un-
desirable reaction to form dimethyl ether (DME) and 2-methyl-2-heptanol (MHOH).
The activation energies of the two competing reactions favour high reactor temper-
ature, which is limited by catalyst activity to 400 K. Yield of the desired product is in-
creased by either increasing reactor size or keeping the concentration of methanol
in the reactor low by using a large excess of MH in the reactor. The latter strategy
results in a large recycle flow rate of MH, which increases separation costs. Thus, this
process demonstrates the classic design trade-off between reactor costs and separation
costs.
• Reactions Involved
1. Desired Reaction
CH3OH + C8H16 − → C9H20O R1 = k1F · xMeOH · xMH − k1R · x MMH
2. Undesired Reaction
2 CH3OH + C8H16 − → C2H6O + C8H18O R2 = k2 · (xMeOH)2
Different R1 R2 Units
Components Forward Backward
kmol
k (Overall Reaction Rate) 6.7 × 107 2.1 × 10−6 1.3 × 109 s·kgcat
kJ
E (Activation Energy) 90000 900 105900 kmol
Concentration xMeOH · xMH x MMH x2MeOH mole fraction
Kinetics Table
W XA 1
= = × )exit
FA −rA −r
exi
A
exi (XA
o t t
where;
W = Catalyst Loading (kg)
XA = Conversion of A
kmol
FAo = Molar Flow Rate of A ( s
)
rA = Net Rate of Reaction of A (kgkmol )
cat
XMeOH = 0.939
kmol
FMeOH = 53.24 h
∴ W = 10422.57 kg
Density of the catalyst = 900 mkg3
10422.57
So, the volume of the catalyst = 3
900 = 11.58 m 11.58
Assuming Void fraction to be 0.45, the volume of fluidised bed will be = = 20.91m3
0.55
Considering Volume of reactor = 1.1 × Volume of bed = 23 m 3
Chapter 7
• Simplification of Representation
• Economic Evaluation
• Environmental Concerns
• Controlability Analysis
• Safety Analysis
• Flexibility
Based on process requirement and above factors the best suited column type here will
be Plate Column.
1. Sieve Trays
3. Valve Trays
Here we would be following Fenske, Underwood, Gilliland Method for calculations
√
In the above table αmean = 2
α1 × α2
Now on the basis of relative volatility calculated;
Nmin = 10 trays
i=1
where,
λ = Latent Heat of Vaporisation of
Feed. Tb = Bubble Point of Feed.
Tf = Feed Temperature.
q = fraction of feed in liquid phase.
zi,F = mole fraction of ith component in feed.
xi,D = mole fraction of ith component in distillate.
αi,HK = relative volatility of ith component w.r.t HK
q = 1.1 (Through latent heat and sensible heat calculations)
From equation 6.1, θ = 2.792
From equation 6.2, Rmin = 1.809
Erbar-Maddox Correlation
R Rmin
Now from above graph, as = 0.68 and = 0.645
R+ Rmin + 1
1
Nmi
= 0.24 ∴N = 41.666 ≈ 42 trays
n
N
• Feed Tray Location
Feed tray location can be calculated using Kirkbride empirical equation
" x
2
Nr
B xf,HK × # b,LK
log10 = 0.206 × log10 xd,HK
Ns D × xf,LK
where,
Nr = number of stages above the feed
Ns = number of stages below the feed.
B = flow rate of bottom product.
D = flow rate of distillate.
xf,HK = mole fraction of HK in feed.
xf,LK = mole fraction of LK in feed.
xd,HK = mole fraction of HK in
distillate.
xb,LK = mole fraction of LK in bottom product.
Nr
= 0.621 (7.4)
Ns
Nr + Ns = 42 (7.5)
From above two equations Ns ≈ 18 ∴ feed tray location will be the 24th tray
uf = K1 r ρl − ρv
× ρv
For calculating K1 we need a graph of K1 vs FLV
K1 vs FLV
V = L + 2.4916 (7.6)
L
= 0.687 (7.7)
V
3
kg
From above two equations V = 7.9606 s = 5.76 ms
For Bottom of the Column:
L = V + 1.98 (7.8)
L
= 1.632 (7.9)
V
From above two equations V = 3.133 kg 3
s = 2.191 ms
∗ Calculating Diameter at top and bottom
V olumetric Flow Rate
Net Area =
V elocity
We will be designing the column at 85% of flooding velocity.
For Top of the Column:
ufo
= 1.589m ∴
s
Net Area = 4.03 m2(Considering 10% downcomer area) & D = 2.26 m
For Bottom of the Column:
ufo
= 1.67m ∴
s
Net Area = 1.457 m2(Considering 10% downcomer area) & D = 1.3624 m
how = 750 × L w
3
ρLlw
where;
Lw = Liquid Flow
Rate
lw = Weir Length
ρL = Density of
Liquid
howmax = 21.24mm
If we see the actual minimum vapour velocity is greater than minimum vapour
velocity (12.12 > 10.33) ∴ our design of column is verified.
• Calculating Column Wall Thickness
Using Carbon Steel as MOC for our column ;
ρMOC kg
= 7850
m3
ft = 6.2 ×107 N
m2
Operating Pressure P = 0.4 atm
External Pressure Pext = 1 atm
Assuming design pressure 10 % more than the operating pressure (IS2825)
PD i
t= + C.A
2ftJ − P
where;
P is design pressure (External or Internal)
J = Joint Efficiency
ft = Permissible Tensile Stress
C.A = Corrosion Allowance
J assumed to be 0.85 assuming spot radiography and C.A as 2mm
Based on Internal Pressure : t = 2.973
mm Based on External Pressure : t =
6.1 mm
Considering thickness based on external pressure as it is greater in magnitude.
PD i
fa =
4t c
where tc is corroded thickness
N
fa = 15.62 m 2
m
∗ Compressive Stress Due to Weight of Shell upto height H
Weight of shell
f1 =
Cross Section of shell
f1 = ρshell × H × g
N
f1 = 0.077H m 2
m
∗ Compressive Stress Due to Liquid, Tray etc
Σ
Liquid Weight Per Unit Height × H
f2 = π Dm t c
where Dm is mean diameter of shell = 2302.77 mm
f2 = H × Ac × Liquid Hold Up
Tray Spacing πD t
m c
Σ
Weight of Attachments Per Unit Height × H
f3 = π Dm t c
Weight of Attachments per unit height = 100mkg
N
f3 = 0.033H m 2
m
∗ Stress Due to Wind
Mw 1.4PwH2
f4 = = πDotc
Z
Where Pw = Wind Pressure = 1300 N
m2
f4 = 0.06136H2
To Find the Height H;
ft = f1 + f4 − f2 − f3
H = 30.27 m
The Original Height of the Column is 28.2 m < 30.27 m
∴ t = 6.1 mm satisfies our calculations.
As per Standards for a column with 2.2987 m diameter the minimum wall
thick- ness should be 19.7 mm which would definitely satisfy the above
calculations as 6.1 mm satisfies all the calculations.
where;
n
Pd n 2ftJ − P
′
t =
N
ft = Permissible Tensile Stress = 62 2
m
tn′ =m0.34mm
Based on market availability tn = 4mm
Now applying area for area method in order to see if any compensation is required
or not;
Reinforcing Pad Limits for Nozzle
A′ = Ah + Ao + Ai = 3151.82 mm2
Now the area for which the compensation is required is;
A = dn × ts′ = 1402.207 mm2
As A < A′ no compensation is required.
Ao 2
p Wop = 8551.787 mm
=
fop
2 π × Ya × G ×
A =
N = 60415.38 mm2 (Ab > & Aatm)
b fatm Aop
G(in cm)
Number of Bolts, n = = 92 bolts
2.54
Further using M16 bolts the area per bolt, Y = 0.51(X)2.09 = 167.56 mm2
In the above equation X = Nominal Diameter in mm = 16mm
Total Bolt Area, A = n × Y = 15415.89 mm2
Here if we see Aop < A < Ab and ∴ it meets the design norms.
Bolt Circle Diameter, B = Go + 2X + 12 = 2382mm
π×
Bolt Pitch, P =
B = 80 mm
n
Outer Diameter of Flange = B + 2X = 2414mm
Flange Plate Thickness Calculation;
P
s
tf = G
Kf
where; ×
• Design of Support
We will be using Skirt Support for the column.
In order to find the skirt thickness
Column Specs
Reactor Design
Reactor becomes an important part of our process as it is the equipment where chemical
reaction will take place.
16. Calculate weight of catalyst using design equation for Continuous Stirred
Tank Re- actor.
19. Calculate the Diameter and length of Reactor using L/D Ratio.
Reactor Design
Different R1 R2 Units
Components Forward Backward
kmol
k (Overall Reaction Rate) 6.7 × 107 2.1 × 10−6 1.3 × 109 s·kgcat
kJ
E (Activation Energy) 90000 900 105900 kmol
Concentration xMeOH · xMH x MMH x2MeOH mole fraction
Kinetics Table
W XA 1
= = × )exit
FAo −rA −r
exi
A
exi (XA
t t
where;
W = Catalyst Loading (kg)
XA = Conversion of A
kmol
FAo = Molar Flow Rate of A ( s
)
rA = Net Rate of Reaction of A (kgkmol )
cat
XMeOH = 0.939
kmol
FMeOH = 53.24 h
∴ W = 10422.57 kg
Density of the catalyst = 900 mkg3
10422.57
So, the volume of the catalyst =
900 = 11.58 m3
Assuming Void fraction to be 0.45
11.58
The volume of fluidised bed = = 20.91m3
0.55
Considering Volume of reactor = 1.1 × Volume of bed = 23 m3
• Mechanical Design of the Reactor Shell
H π
Assuming = 3 & through equation, V = × D2 × H ⇒ H = 6.5 m and
D 4
D=2m
Here we will be using Carbon Steel as our MOC
∴ ρMOC kg
= 7500
m3
N
ft = 100 2
m
m
Density of different components involved in the process;
Components Density ( kg )
m3
MMH 877.392
MH 717.899
MeOH 789.219
DME 657.674
MHOH 866.67
Density Table
Thickness of shell;
PD i
t= + C.A
2ftJ − P
where;
P is design pressure (External or Internal)
J = Joint Efficiency
ft = Permissible Tensile Stress
C.A = Corrosion Allowance
Assuming design pressure 10 % more than the operating pressure (IS2825)
Operating Pressure = 15 atm
J assumed to be 0.85 assuming spot radiography and C.A as 2mmt = 21.86
mm ⇒ 25 mm (Based on Market Availability)
tc = Corroded Thickness = 23 mm
P (Di + tc)
f c= 2tc
N
fc = 73.525 m 2
m
∗ Axial Stresses Calculation
PD i
f 1 = 4tc
N
f1 = 36.345 m 2
m
⬦ Compressive Stress Due to Self Weight and Contents Weight
f2 = W
πt (D + t )
c i c
Where W includes weight of all the things in the reactor along with the self
weight of the reactor.
π 2 2
W1 = ρMOC × Vshell & Vshell = 4
× (do − di )
Now fa = f1 − f2 + f3
N
fa = 36.345 - 2.49 = 33.847 m 2
m
∗ Offset Piping Stress
As compared to other stresses involved this stress becomes negligible ∴ we are
N
neglecting it. fs = 0 m
2
m
Now after calculating all the stresses
q we calculate the resultant stress.
N
fr = 63.74 m 2
m
Conditions to be checked;
1. fc < ft → (73.525 m
N
2 < 100 m
N
2 )
m m
2. fa < ft → (33.847 N
m 2 < 100 m
N
2 )
m m
3. f2 < fcper → (2.49 N
m 2 < 39.74 m
N
2 )
m m
For the last condition we need to find fcper which is given by;
1 √ E t
fcper = × ×
Do
12 3(1 − µ2) 2
where;
fcper = Permissible Compressive Stress
E = Modulus of Elasticity = 35 GPa
µ = Poissons Ratio = 0.31
N
fcper = 39.744 m 2
m
∴ The thickness of 25mm is ideal for this reactor shell.
Pjdi
tj =
2ftJ − Pj
where;
Pj = Jacket fluid Pressure(Considering 10% extra as per IS2825)
di = Inside diameter of jacket = 2.25m
J = Joint efficiency = 0.85 (Assuming spot radiography)
tj = 2.95mm
As theoretically jacket thickness comes out to be 2.95mm we will take a jacket with
thickness of 4mm depending on market availability.
Pd n
tn ′ = 2ftJ − P
where;
P = Nozzle design
pressure
J = Joint Efficiency = 1
N
ft = Permissible Tensile Stress = 100 2
m
tn′ = 1.686mm
m
A′ = Ah + Ao + Ai = 887.4 mm2
Now the area for which the compensation is required is;
A = dn × th′ = 3344 mm2
As A > A′ ∴ compensation is required;
∴ now H1 = 2.5tn = 10mm
As a result now Ao = 46.28 mm2 and A′ = 702.28 mm2
As a result compensation to be given = A − A′ = 2641.72 mm2
4
Aat 2
m Watm = 24138.368 mm
=
fatm
Ao 2
p Wop = 98154.4 mm = A max
=
fop
2 π × Ya × G ×
A
N = = 122219.58 mm2 (Ab > & Aatm)
b fatm Aop
G(in cm)
Number of Bolts, n = = 80 bolts
2.54
Now keeping in mind all the factors we take 124 bolts in order to satisfy all the
design conditions
Further using M34 bolts the area per bolt, Y = 0.51(X)2.09 = 809.77 mm2
In the above equation X = Nominal Diameter in mm = 34mm
Total Bolt Area, A = n × Y = 100411.348 mm2
Here if we see Aop < A < Ab and ∴ it meets the design norms.
Bolt Circle Diameter, B = Go + 2X + 12 = 2120mm
π×
Bolt Pitch, P =
B = 53.71 mm
n
Outer Diameter of Flange = B + 2X = 2188mm
where;
G = Mean Diameter of Gasket Load Reaction = 2020mm
P = Design Pressure = 1.67 mN 2
m
f = Permissible Tensile Stress = 100mN 2
m
1
K= 1.5WmaxhG
0.3 + HG
where,
Wmax = Wop = 5692958.6 N
B−
h = Radial Distance from Gasket Load Reaction to Bolt Circle = = 50mm
GG
π 2
2
H = Total Hydrostatic End Force = G P = 5357882.398 N
4
K = 2.946
Now after finding K we will find tf ;
tf = 152mm
Based on market availability 160 mm should be used.
• Design of Support
We will be using Bracket Support for the reactor;
The material of construction will be structural steel IS-800
Permissible Stresses for Structural Steel IS-800:-
Tension = 140 mN 2
m
Compression = 123.3 mN 2
m
N
Bending = 157.5 m 2
m
Number of Brackets, n = 8
Figure 7.9: Vessel Support on Brackets
B2 4
a4 4
a +B
where; 1
fb = 0.7 × Pav × 2
×
T
fb = Permissible Bending Stress
a, B = Base Plate Dimensions
T1 = Thickness of Base Plate
Now, PFT
=
Pav a×B
For finding Pavg we need PFT ;
Σ
4 × fw × (H − F ) W
PF = n × Db + n
T
where;
PFT = Total Force
F = Bottom Clearance = 1000mm
H = Height of Vessel above the Foundation = 7500mm
Db = Diameter of Bolt Circle = 1650mm
Σ
W = Weight of Vessel with its Contents = 37255.57 kg
n = Number of Brackets = 8
fw = Force due to Wind Load = 0.7 × p × h × Do = 12125.75 N
PFT = 70453.5 N
Now for finding base plate thickness we consider base plate dimensions as
a = 140 mm & B = 150 mm
∴ 140907.0311
= = 3.355 mN 2
Pav 140 × 150 m
Now as fb = 157.5 and other information known, T1 = 12.02 mm
As per standards and market availability we should use 14 mm thick base plate.
Cost Estimation
A plant design must present a process that is capable of operating under conditions
which will yield a profit. Since net profit equals income minus all expenditure,
it is essential that chemical engineers should be aware of various costs involved in the
manufacturing process. Capital must be allocated for direct plant expenses, such as
those for raw materials, labour, and equipment. Besides direct expenses, many other
indirect expenses are incurred, and these must be included if a complete analysis of the
total cost is to be obtained. Some examples of these indirect expenses are administra-
tive salaries, product-distribution costs, and costs for inter-plant communications.
A capital investment is required for any industrial process, and determination of the
necessary investment is an important part of a plant-design project. The total
investment for any process consists of fixed-capital investment for physical equipment
and facili- ties in the plant plus working capital which must be available to pay salaries,
keep raw materials and products on hand, and handle other special items requiring a
direct cash outlay. Thus, in an analysis of costs in industrial processes, capital-
investment costs, manufacturing costs, and general expenses including income taxes
must be taken into consideration.
Purchased Equipment Cost consists of the cost of all the major and minor equip-
ments that are needed for the smooth functioning of the plant like reactor, distil-
lation column, heat exchanger, pumps etc.
1. Cost of Reactor
In 1987 the cost of a similar reactor of capacity 2000 gallons was 32000$ Capacity
of Reactor to be purchased in 2023 is 3170 gallons.
Chemical Engineering Plant Cost Index in 1987 = 320
Chemical Engineering Plant Cost Index in 2023 = 607.5
Now;
0.53
Reactor Cost = Cost in 1987 × Capacity in 2023 Cost Index in 2023
×
Capacity in 1987 Cost Index in 1987
Reactor Cost = 77546$ = 59,38,666 Rs
2. Cost of Column-1
• Tower Cost
In 1987 the cost of distillation tower of 6 ft diameter was 900 $f
t
Diameter of Column-1 = 3.54 ft
Chemical Engineering Plant Cost Index in 1987 = 320
Chemical Engineering Plant Cost Index in 2023 = 607.5
Now;
0.78
Tower Cost = Cost in 1987× Diameter in 2023 × Cost Index in 2023
Diameter in 1987 Cost Index in 1987
Tower Cost = 1132.14 f$
t
Height of Tower = 10.05
m
∴ Tower Cost = 37328.94 $
• Tray Cost
$
In 1987 the cost of sieve trays with column diameter 3.54 ft was 450 tra
y
Chemical Engineering Plant Cost Index in 1987 = 320
Chemical Engineering Plant Cost Index in 2023 = 607.5
Number of trays in Column-1, n = 12
∗ Condenser Cost
kW
Heat Transfer Coefficient = 0.852 m 2K
Differential Temperature = 8 K
Duty = 4.4 kW
∴ Area = 0.6455 m2
Cost Index in 2023
Capital Cost = 729.6 A0.65
Cost Index in 2010
× ×
Capital Cost = 606.3 $
∗ Reboiler Cost
Heat Transfer Coefficient = 0.568 kW
m 2K
• Tower Cost
In 1987 the cost of distillation tower of 6 ft diameter was 900 $f
t
Diameter of Column-2 = 7.54 ft
Chemical Engineering Plant Cost Index in 1987 = 320
Chemical Engineering Plant Cost Index in 2023 = 607.5
Now;
• Tray Cost
$
In 1987 the cost of sieve trays with column diameter 7.54 ft was 900 tra
y
Chemical Engineering Plant Cost Index in 1987 = 320
Chemical Engineering Plant Cost Index in 2023 = 607.5
Number of trays in Column-2, n = 42
∗ Condenser Cost
kW
Heat Transfer Coefficient = 0.852 m 2K
∗ Reboiler Cost
kW
Heat Transfer Coefficient = 0.568 m 2K
4. Cost of Column-3
• Tower Cost
In 1987 the cost of distillation tower of 6 ft diameter was 900 $f
t
Diameter of Column-3 = 5.21 ft
Chemical Engineering Plant Cost Index in 1987 = 320
Chemical Engineering Plant Cost Index in 2023= 607.5
Now;
0.78
Tower Cost = Cost in 1987× Diameter in 2023 × Cost Index in 2023
Diameter in 1987 Cost Index in 1987
Tower Cost = 1530.43 f$
t
Height of Tower = 16.1 m
∴ Tower Cost = 80838.66 $
• Tray Cost
$
In 1987 the cost of sieve trays with column diameter 5.21 ft was 700 tra
y
Chemical Engineering Plant Cost Index in 1987 = 320
Chemical Engineering Plant Cost Index in 2023= 607.5
Number of trays in Column-3, n = 22
∗ Condenser Cost
kW
Heat Transfer Coefficient = 0.852 m 2K
Here we assume that Purchased Equipment Cost is 30% of Fixed Capital Invest-
ment.
∴ Fixed Capital Investment = 40 Crores.
Now Fixed Capital Investment comprises of Direct Cost and Indirect Cost.
Also here we consider Working Capital Investment as 15% of Total Capital
Investment.
All the below calculations will be on annual basis.
Direct Costs
Direct Costs are the cost of all the materials and labour involved in actual installation
of complete facility.
Breakdown of the Direct Costs is as follows;
Now if we see;
Fixed Capital Investment = Direct Cost + Indirect Cost
Fixed Capital Investment = 31.6 + 8.4 = 40 Crores
Now,
Total Capital Investment = Fixed Capital Investment + Working Capital Investment
Total Capital Investment = 47 Crores
Direct Production Costs are the expenses directly associated with the manufacturing
operation. This type of cost involves expenditures for raw materials (including trans-
portation, unloading, etc.,); direct operating labor; supervisory and clerical labor directly
connected with the manufacturing operation; plant maintenance and repairs; operating
supplies; power; utilities; royalties; and catalysts.
Here we assume that the Total Product Cost = x Crores
3. Maintenance Cost
We assume the Maintenance Cost to be 6% of the Fixed Capital Investment.
As Fixed Capital Investment = 40 Crores
∴ Maintenance Cost = 2.4 Crores
5. Utility Cost
We assume the Utility Cost to be 10% of the Total Product Cost
∴ Utility Cost = 0.1x Crores
6. Direct Supervision Cost
We assume the Direct Supervision Cost to be 17% of Operating Labour Cost.
As Operating Labour Cost = 0.13x Crores
∴ Direct Supervision Cost = 0.0221x Crores
7. Laboratory Charges
We assume the Laboratory Charges to be 15% of the Operating Labour Cost.
As Operating Labour Cost = 0.13x Crores
∴ Laboratory Charges = 0.0195x Crores
Certain expenses are always present in an industrial plant whether or not the manufac-
turing process is in operation. Costs that are invariant with the amount of production
are designated as fixed costs or fixed charges. These include costs for depreciation, local
property taxes, insurance, and rent.
1. Depreciation Charges
We assume the Depreciation Charges to be 10% of the Fixed Capital Investment.
As Fixed Capital Investment = 40 Crores
∴ Depreciation Charges = 4 Crores
2. Local Taxes
We assume the Local Taxes to be 2.5% of the Fixed Capital Investment
Fixed Capital Investment = 40 Crores
∴ Local Taxes = 1 Crore
3. Insurance Charges
We assume the Insurance Charges to be the 0.7% of the Fixed Capital Investment
Fixed Capital Investment = 40 Crores
∴ Insurance Charges = 0.28 Crores
The expenditures required for routine plant services are included in plant overhead
costs. Non-manufacturing machinery, equipment, and buildings are necessary for many
of the general plant services, and the fixed charges and direct costs for these items are
part of the plant-overhead costs. Expenses connected with the following comprise the
bulk of the charges for plant overhead:
• Safety services
• Control laboratories
• Distribution of utilities
• Warehouses
• Lighting
For Plant Overhead Charges we assume the charges to be 60% of the sum of the
Operating Labour Charges, Direct Supervision Charges and Maintenance Charges
As;
Operating Labour Charges = 0.13x Crores
Direct Supervision Charges= 0.0221x Crores
Maintenance Charges = 2.4 Crores
∴ Plant Overhead Charges = 0.09126x + 1.44 Crores
As a result;
Total Manufacturing Cost = 0.59286x + 9.48 Crores
General Expenses
In addition to the manufacturing costs, other general expenses are involved in any com-
pany’s operations. These general expenses may be classified as administrative
expenses, distribution and marketing expenses, research and development expenses,
financing ex- penses.
1. Administrative Expenses
We assume the Administrative Expenses to be 15% of the sum of the Operating
Labour Charges, Direct Supervision Charges and Maintenance Charges
Operating Labour Charges = 0.13x Crores
Direct Supervision Charges= 0.0221x Crores
Maintenance Charges = 2.4 Crores
∴ Administrative Expenses = 0.022815x + 0.36 Crores
4. Financing Charges
We assume Financing Charges to be 5% of the Total Product Cost
∴ Financing Charges = 0.05x Crores
General Expenses in a Nutshell
Now;
The total income minus the total production cost gives the gross earnings
made by the particular production operation, which can then be treated mathematically
by any of several methods to measure the profitably of the proposed venture or project.
Because of income-tax demands, the final net-profit is often much less than the gross
earnings. Income-tax rates are based on the gross earnings received from all the company
interests. Consequently, the magnitude of these costs varies widely from one company
to another.
The Price of MMH = 11235 kRs
g
kg
The Annual Production of MMH = 51000 year
∴ Total Income = 58 Crores
Gross Earning = Total Income − Total Product Cost
Crores
These 3 factors help us know how profitable our plant can be and in how much time
can we recover our capital invested back. These factors help us measure whether our
investment would prove to be good or not.
Turnover Ratio
Turnover ratio is defined as the ratio of gross annual sales to the fixed-capital in-
vestment, where the product of the annual production rate and the average selling
price of the commodities is the gross annual sales figures. The reciprocal of the
turnover ratio is sometimes defined as the capital ratio or the investment ratio.
Turnover ratios of up to 5 are common for some business establishments and some are
as low as 0.2.
Rate of Return
Net Profit
Rate of Return = 100
Fixed Capital Investment
×
10.56
40 × 100 = 26.4%
Rate of Return =
Payback Period
Refreence
1. William L. Luyben
Design and Control of the Methoxy-Methyl-Heptane Process; American Chemical
Society, 2010
4. https://www.marketwatch.com/press-release/global-2-methoxy-2-methylheptane-cas-
76589-16-7-market-size-2020-emerging-trends-industry-share-future-demands-market-
potential-traders-regional-overview-and-swot-analysis-till-2025-2020-04-14
5. http://m.molbase.com/moldata/2132090.html
6. http://www.molbase.com/en/synthesis 76589-16-7-moldata-2132090.html
7. https://pubchem.ncbi.nlm.nih.gov/compound/2-Methoxy-2-methylheptane
8. https://m.sigmaaldrich.com/IN/en/search?term=methoxy+heptane
9. http://www.chemspider.com/Chemical-Structure.28698534.html