THF FinalReport
THF FinalReport
THF FinalReport
Spring 2019
Group 05
Production of Tetrahydrofuran
Report - Final Stage Evaluation
Members:
1
Contents
List of figures 5
List of tables 7
Introduction 9
1.1 Tetrahydrofuran 9
1.3 Specifications 12
4. Process Flowsheet 21
5.1 Reactor 22
5.4.1 DC101 23
5.4.2 DC102 23
6.1 Reactor 24
6.3 Pump 26
2
6.4 Compressor 27
8.1 Reactor 33
9. Sensitivity Analysis 40
9.1 DC101 40
9.2 DC102 42
10.1 Reactor 44
Cost of Vessel 44
Cost of Catalyst 46
10.2 Flash 47
Size Estimation: 51
Sizing of shell 51
Cost Estimation 54
3
Cost of shell 54
Cost of Reboiler 55
Cost of Condenser 57
10.5 Compressor 60
Cost Estimation 61
10.6 Pump 62
Cost Estimation 63
4
12.2 Control Strategy for F102 System 74
12.2.1 Cooler 74
5
List of figures
1.1 2D and 3D conformer of Tetrahydrofuran 09
1.2 Products and byproducts of Tetrahydrofuran 11
2.1 Total VOlume & value of imports in India 14
2.2 Trend of material used in sportswear 15
2.3 World consumption of THF 15
2.4 Forecasted Market size and CAGR of THF globally 16
2.5 Major importers of THF 17
2.6 Major exporters of THF 17
3.1 Main Reaction 18
3.2 First Side Reaction 19
3.3 Second Side Reaction 19
3.4 Third Side Reaction 19
3.5 Block Diagram for Flowsheet development 20
4.1 Process FLow Diagram 21
7.1 Aspen Plus Flowsheet (Part 1) 29
7.2 Aspen Plus Flowsheet (Part 2) 30
8.1 Reactor 33
8.2 Flash F101 34
8.3 Flash F102 35
8.4 Distillation Column DC101 36
8.5 Distillation Column DC102 37
9.1 Change in purity of THF in Distillate of DC-102 on changing input
feed stage in DC-101 40
9.2 Change in mass flowrate of THF in Distillate of DC102 on changing
input feedstage in DC-101 40
9.3 Change in mass flow rate of water in bottom of DC-101 on changing
input feed stage in DC-101 41
9.4 Change in mass flow rate of THF in Distillate of DC-101 on changing
input feed stage in DC-101 41
9.5 Change in mole fraction of water in bottom of DC-101 on changing
input feed stage in DC-101 42
9.6 Change in purity of THF in Bottom of DC-102 on changing input
feed stage in DC-102 42
9.7 Change in mass flow rate of THF in Bottom of DC-102 on changing
input feed stage in DC-102 43
10.1 Design Velocity Factor v/s Vapour-Liquid Separation Factor Curve 48
10.2 Graph for calculating Flooding Velocity 52
6
10.3 Cost of shell against weight for Distillation column 54
10.4 Installed cost of sieve trays against column diameter 55
10.5 Cost of shell & tube heat exchanger against area 56
10.6 Isentropic Efficiency v/s Compression ratio 60
12.1 Degrees of Freedom Analysis of Distillation Column 69
12.2 Pressure Control 71
12.3 Reboiler Control 72
12.4 Temperature Profile of DC101 73
12.5 Reflux Ratio Variation 73
12.6 Pump Control 74
12.7 Overall Control Diagram 75
14.1 Heat Exchanger Network 85
List of Tables
1.1 General Properties of Tetrahydrofuran
10 1.2 General Product specification
12 6.1 Comparison between Centrifugal & Reciprocating Pump
24
6.2 Heat Exchanger Comparison 25
6.3 Comparison between Centrifugal & Reciprocating Pump 26
6.4 Comparison between Centrifugal & Reciprocating Compressor 27
7.1 Stream Summary 31 32
8.1 Reactor stream summary 33
8.2 Flash 101 Stream Summary 34
8.3 Flash 102 Stream Summary 35
8.4 DC101 Stream Summary 36
8.5 DC 102 Stream Summary 37
10.1 Catalyst Composition 46
10.2 Catalyst Cost 47
10.3 Flash Vessel Sizing 50
10.4 Flash Costing 51
10.5 Distillation Column Information - DC101 52
10.6 DC101 and DC 102 Design Parameters 53
10.7 Areas of Reboiler and Condenser 53
10.8 Shell cost of DC-101 and DC-102 54
10.9 Total cost of Trays - DC101 and DC102 55
10.10 Reboiler cost of DC101 and DC102 56
10.11 Condenser cost of DC101 and DC102 57
7
10.12 Heat exchanger type factor 58
10.13 Heat exchanger pressure factor 58
10.14 Heat exchanger Material factor 58
10.15 Heat exchanger double pipe pressure factor 59
10.16 Heat Exchanger Costing 59
10.17 Power requirement of Compressor 61
10.18 Cost of Compressor 61
10.19 Power requirement of Pumps 62
10.20 Cost of Pumps 63
11.1 All Major Equipment Cost 64
11.2 Total Capital Investment 65
11.3 Revenue Calculation 66
11.4 Raw Material Cost 66
11.5 Total Product Cost 66
12.1 Variables in Process and Instrumentation Diagram 70
13.1 Emissions in kilograms/day 76
13.2 Emission from stream no. 31 77
13.3 Emission from stream no. 28 78
13.4 Emission from stream no. 22 79
13.5 Emission from stream no. 20 80
13.6 Emission from stream no. 14 81
14.1 Stream listing for heat exchanger network 82
14.2 Stream listing for heat exchanger network below pinch point 83
14.3 Stream listing for heat exchanger network above pinch point 84
14.4 Utilities cost after heat integration 85
8
1. Introduction
This chapter is aimed towards providing a brief overview about the subject and also the
motivation for the need of detailed study which includes restropection of trends and
analysis of future prospects of tetrahydrofuran.
1.1 Tetrahydrofuran
THF is composed of 4 carbon atoms, 8 hydrogen atoms, and 1 oxygen atom. It is most
importantly characterized as being a heterocyclic ether. The ring in the compound is
non-planer. It is an organic compound, cyclic aliphatic ether, with the chemical formula
C4H8O, with IUPAC name Oxolane.
THF is a highly volatile, chemically neutral compound with low boiling point and odor
similar to acetone. Important structural property of THF is that it is a polar molecule. Since
oxygen is more electronegative than carbon and hydrogen, the oxygen pulls electron
density towards itself in the form of the covalent bonds within the ring structure. As a
consequence of this, polarity is introduced and hence THF mixes readily with other polar
species like water, methanol, ethanol, and acetonitrile, just to name a few. On contact with
air it readily oxidizes to form peroxides, which can be inhibited to a certain extent by
adding stabilizers.[3]
1
[1.1] THF Image on Indiamart https://www.indiamart.com/proddetail/tetrahydrofuran-8120213188.html
2
[1.2] THF 3D conformer Image https://www.indiamart.com/proddetail/tetrahydrofuran-4621862848.html
9
Table 1.1: General Properties of Tetrahydrofuran [1.3]3
Boiling Point 66 °C
Flash Point -7 °C
1. Solvent - THF is most commonly used as organic solvent. It’s use as a solvent is
widespread from small scale laboratory experiments to manufacturing of PVC (Poly
Vinyl Chloride) which indeed has its use as material of plumbing and piping
equipments, healthcare devices, cable insulation etc.
3
[1.3]THF Properties - https://pubchem.ncbi.nlm.nih.gov/compound/tetrahydrofuran#section=Boiling-Point
4
[1.4] Application of THF - http://www.basf-petronas.com.my/sites/default/files/Tetrahydrofuran.pdf
5
[1.5] Application of THf - https://pubchem.ncbi.nlm.nih.gov/compound/tetrahydrofuran#section=MS-MS
10
3. Reaction medium - The property of being a powerful solvent is utilized by the
pharmaceutical industry. THF is used in many reactions like Grignard reaction,
Wittig reaction etc. Many Grignard and organometallic compounds can be obtained
using THF as the reaction medium.
6
[1.6] Products & Byproducts of THF
http://www3.hkexnews.hk/listedco/listconews/sehk/2010/1130/01887_947325/ewf114.pdf?fbclid=IwAR
2mo8KMPJsOywa-8ikyXgIjGsPiwk816KvzMd-FSJtJoCfJIRJEU6oj9Q0
11
1.3 Specifications
Tetrahydrofuran has various uses, and product specification needs to satisfied on several
grounds. On contact with air it readily oxidizes to form peroxides, which is required to be
inhibited to a certain extent by adding stabilizers. As an antioxidant, 250 ppm of BHT
(2,6-di-tert-butyl-4-methylphenol) is usually added.
Item Specification
7
[1.7] Product specification
https://www.m-chemical.co.jp/en/products/departments/mcc/c4/product/1201006_7922.html
12
2. Market Survey & Global Demand
In this chapter we will see the demand and supply scenario of THF over the years which
also includes the market survey of chemical plants actively producing THF to meet the
global demands. We will also analyse the price trends and the products influencing the
production.
2. Chinese market has expanded aggressively in all domains (includes rapid increase in
industries producing chemicals along with manufacturing industries consuming it),
hence Indian industry faces fierce competition from cheap Chinese products
3. Foreign suppliers from China, Vietnam, Korea etc. have reduced the price of Spandex
to ₹ 100/kg (2016) and import duty being reduced to 1% has the worsen the
situation for Indian suppliers affecting the profitability of business.
4. India must revise its anti dumping duty to protect the Indian manufacturers
8
[2.1] Article on Spandex production in India
http://www.indiantextilemagazine.in/corporate-news/indorama-seeks-sops-to-encourage-domestic-spande
x-production/
13
The below mentioned image gives a trend of total volume and value of imports between
February 2013 & November 2016
The total value of imports of THF market in India is about $27,605,949, total quantity is
about 12,478,750 kg as of Nov 2016. [2.3]9
China is the largest consumer of THF, refer the figure below to analyse the market
consumption
9
Trend of volume and value of THF import in India
https://www.zauba.com/importanalysis-tetra+hydro+furan-report.html
10
Trend of materials used in sportswear
http://emergingtextiles.com/?q=art&s=170321-spandex-market-price&r=free&i=samplearticle
14
Fig 2.2: Trend of material used in sportswear
The production and increase in of goods like - cosmetics, adhesives, surface coatings, inks,
PVC (Poly Vinyl Chloride) etc. also have contributed in increasing the demand of THF.
In the period of 2012-15 the global market has seen a growth of about 40% but in coming
years we can expect the growth by just 5%, with Chinese market witnessing the highest
growth. The largest producers of THF globally are- BASF and Dairen Chemical Corporation.
11
[2.5] world consumption of THF
https://ihsmarkit.com/products/tetrahydrofuran-chemical-economics-handbook.html
15
Projected growth of THF market is from USD 2.95 Billion in 2017 to USD 4.12 Billion in
2022 at CAGR of 6.9%, with China as the country with highest projected CAGR. Asia Pacific
is the largest market of THF globally. Due to boost in manufacturing sector industries like
chemical, automotive, packaging etc. countries like India have also contributed in creating
demand of THF.
Asia Pacific is dominating as global hub both in terms of volume and value.
Fig 2.4: Forecasted Market size and CAGR of THF globally [2.6]12
2. India ($37M)
3. Belgium-Luxembourg ($22.1M)
4. Japan ($13.3M)
12
[2.6] Forecasted market size and CAGR of THF globally
https://www.marketsandmarkets.com/Market-Reports/tetrahydrofuran-thf-market-180054945.html
16
Fig2.5 : Major importers of THF [2.7]13
1. Germany ($80.6M)
3. Malaysia ($21M)
5. China ($16.9M)
13
Global importers of THF https://atlas.media.mit.edu/en/profile/hs92/293211/
14
Global exporters of THF https://atlas.media.mit.edu/en/profile/hs92/293211/
17
15
3. Production Process of THF
In this process we are using, one main reaction and three side reactions are taking place.
All the reactions are happening in a single reactor. The main reaction is a 3-step reaction.
Besides the main reaction, 3 side reactions are also happening. Although, main reaction has
a very high selectivity i.e. about 98% and rest of the reactions have the conversion of the
order of 10-2.
15
[3.1 ]The hydrogenation of maleic anhydride to γ-butyrolactone using mixed metal oxide catalysts in a
batch-type reactor
https://www.sciencedirect.com/science/article/pii/S0378382005002055
18
This formed acid is then hydrogenated to form n-butanol.
Second side reaction happening is forming n-propanol as the side product. In this reaction,
first step is same as main reaction, maleic anhydride is being hydrogenated in the presence
of metal catalyst to obtain succinic anhydride. In second step, succinic anhydride is
subjected to hydration, propionic acid along with carbon dioxide is formed as a result. In
last step, Propionic acid undergoes hydrogenation, water is removed and formation of
n-propanol takes place.
Third side reaction which is happening is consuming the product of main reaction i.e. THF.
In this reaction THF is reacting with hydrogen to form butane, butanol and water as the
products.
19
Fig 3.5 Block Diagram for Flowsheet development
20
4. Process Flowsheet
21
5. Process Description of Major Equipments
5.1 Reactor
Mixture of Hydrogen gas (from feed and recycle) and Maleic Anhydride (stream number-4,
Pressure-40, Temperature-220°C, phase-liquid) are fed to the reactor which operates at
pressure 4 0 bar and t emperature 240°C . The reactor is fixed bed adiabatic type of reactor
which is producing THF (Tetrahydrofuran), along with unreacted Maleic Anhydride and
Hydrogen; side products in outlet stream are Butane, Butanol, Water and Propanol. The
outlet stream (stream number 8) is passed through a heat exchanger before doing flash
separation. Catalyst used is oxides of different materials at different percentages, base used is of
aluminium {oxides of- Copper(45%), Zinc(20%), Alumina(30%) and Chrome(5%)}
22
5.4 Distillation System
We are using pressure swing distillation as we need to separate azeotropic mixture of
water and THF.
5.4.1 DC101
The DC101 column is largely to separate water from the feed mixture. The distillation
column operates at ambient pressure. The feed stream (stream number 13) is fed at stage
number 5 . The distillate (stream number 15) contains large fraction of THF and small
fraction of Butane and Water, which is cooled before doing flash separation. The bottoms
stream (stream number 14) which contains large fraction of water along with some
amount of Butanol and Maleic Anhydride.
5.4.2 DC102
The distillation column DC102 functions at pressure 8 bar . The input stream (stream
number 19) containing majorly THF and some amount of THF and Water is fed at stage
number 4 . The purpose of this distillation column is to first produce THF of required purity
and second to recycle back some condensable containing mixture of Water and THF as part
of swing distillation system. The distillate (stream 21) contains large fraction of THF and
small fraction of Water and Butane, from which Butane is separated using separator and
purged, rest is cooled and recycled back to DC101. The bottoms (stream number 20)
contains our product THF of p urity 9
8.4%.
23
6. Process Equipment Selection
16
6.1 Reactor
Table 6.1: Comparison between Centrifugal & Reciprocating Pump
In our reactor we have three phases. Maleic Anhydride is coming as a liquid, hydrogen as
gas and the reaction is happening over a solid catalyst.
The fixed bed reactor is modelled as an PFR and Slurry Type reactor is modelled as an
CSTR. The patent state the PFR has better efficiency and conversions, hence we are using a
3 Phase Fixed bed reactor.
16
[6.1]Visual Encyclopedia of Chemical Engineering/Slurry Reactor
http://encyclopedia.che.engin.umich.edu/Pages/Reactors/Slurry/Slurry.html
[6.2] Visual Encyclopedia of Chemical Engineering/Fixed Bed Reactor
http://encyclopedia.che.engin.umich.edu/Pages/Reactors/Tricklebed/Tricklebed.html
24
6.2 Heat Exchanger
A comparison between Shell & Tube heat exchanger & Plate type heat exchanger has been
done to decide which type of heat exchanger is to be used.
Heat Exchanger
Advantages Disadvantages
Type
● Less expensive than plate type heat ● Less efficient than plate type heat
exchanger exchanger
Shell & Tube ● Pressure drop across the tube is less ● Its capacity can not be increased
Heat Exchangers ● Can be operated in systems with high ● Occupies more space as compared
temperature & pressure conditions to plate type heat exchanger
● Tube leaks are easy to locate and ● Cleaning and maintenance is
repair since pressure test is easy comparatively difficult
● It’s compact in size and no extra space ● Initial investment is high because of
is required for dismantling expensive plate material
● Higher heat transfer efficiency than ● Locating leakage is difficult since
Plate Type Heat shell & tube heat exchanger pressure test is difficult to perform
Exchangers ● Capacity can be increased by ● Bonding material between plates
introducing new plates cause hindrance to operating temp,
● Is easy to clean & maintain ● Pressure drop is higher than tube
exchanger
As we can see from the above mentioned table plate heat exchanger is more efficient than
shell & tube heat exchanger but it is also more expensive than the shell and tube heat
exchanger. Plate heat exchanger is more efficient than shell and tube exchanger and it’s
capacity can also be increased by adding more plates.
Although plate heat exchanger is better in terms space occupancy and efficiency, we will be
using she17ll and tube heat exchanger, since it is cheaper, offers lower pressure drop,
detection and rectification of leakages can be done easily, also it can be operated at high
temperature and pressure . These advantages gives it an edge over plate heat exchanger.
17
[6.3] Advantages and Disadvantages of Shell and Tube and Plate type Heat Exchangers
https://marineengineeringonline.com/advantages-disadvantages-shell-tube-plate-type-heat-exchangers/
[6.4] “General Engineering Knowledge” by H.D. McGeorge
25
6.3 Pump
A comparison between Centrifugal & Reciprocating pumps has been done to decide which
type of pump is to be used.
It is obvious from the comparison that centrifugal pumps should be used because of its
ability to discharge large quantity of fluid continuously. In addition to that, it is also
cheaper than reciprocating pumps.
18
[6.5] ADVANTAGES AND DISADVANTAGES OF CENTRIFUGAL PUMPS
http://constructionmanuals.tpub.com/14265/css/Advantages-and-Disadvantages-of-Centrifugal-Pumps-21
5.htm
[6.6] Comparison between Centrifugal pumps and Reciprocating Pumps
https://theconstructor.org/water-resources/comparison-between-centrifugal-reciprocating-pumps/4786/
26
19
6.4 Compressor
A comparison between Centrifugal & Reciprocating Compressor has been done to decide
which type of compressor is to be used.
Although reciprocating compressor is cheaper and can be used to develop high pressure in
comparison to centrifugal compressor, it has a major disadvantage that it is suitable for
continuous and high flow rates. That is why, centrifugal compressor should be used.
19
[6.7] “Process Equipment Selection” CL 415 slides
https://moodle.iitb.ac.in/pluginfile.php/206689/mod_resource/content/1/Compressors%202%20-%20Fin
al.pdf
27
7. ASPEN PLUS Flowsheet
To resemble actuality and also account for energy balance, we simulate our process plant in
ASPEN PLUS process simulator to get better and real life matching insights in the plant. We
initially describe the property method chosen and then look at the stream summary to get
an overview of products generated.
Feed Streams - Hydrogen is obtained at 25 bar in vapor phase and Maleic Anhydride
(MAH) is obtained at STP, which is in solid state and liquified using a heater in a vessel.
Reactor - We have modelled reactor as RStoic as we had conversion data. Later on we
didn’t get equilibrium data so finally we have kept it as RStoic. The H2:MAH ratio in the
reactor is 30:1 and it lies in the prescribed range mentioned in the patent.
Distillation Columns - We could not use DSTWU block in ASPEN because DSTWU fails
most of the time in case of azeotropic systems. So we did a logical hit and trial and found
out optimal number of stages and reflux ratio. Also, we have 2 distillation columns which in
combination are used to separate azeotropic components using Pressure Swing concept. By
altering pressure, we get different composition of THF and water, as composition of
azeotropes change with pressure.
20
http://jingwang.yolasite.com/resources/Senior%20Design.pdf | Synthesis of Tetrahydrofuran
28
7.3 ASPEN PLUS Flowsheet
29
Fig 7.2 Aspen Plus Flowsheet (Part 2)
30
7.3 Stream Summary
Table 7.1: Stream Summary
31
32
8. Key Equipments and Mass Balance
8.1 Reactor
Mass Fractions
MALEI-01 1 0 0.01
1-PRO-01 0 0 0
33
8.2 Flash F101
Mass Fractions
1-PRO-01 0 0 0
CARBO-01 0 0 0
BUTANE 0 0 0
34
8.3 Flash F102
Mass Fractions
MALEI-01 0 0 0
HYDRO-01 0 0 0
BUTANOL 0 0 0
1-PRO-01 0 0 0
CARBO-01 0.00 0 0
35
8.4 Distillation Column - DC101
Mass Fractions
36
8.5 Distillation Column - DC102
Mass Fractions
37
8.6 Overall Mass Balance
This chapter concerns with doing linear mass balance. To satisfy the degrees of freedom we
take an assumption of sharp separation in the product streams of Flash and Distillation
Column.
Points to note:
1. Basis chosen - H2:MAH = 100:1 such that feed to the reactor is 101 moles/hr
2. The distillate of DC-102 will contain THF and water in the azeotropic composition at the
pressure at which the DC-102 is operated. (Azeotropic Mole fraction of THF is 0.64 at
P=8 bar)
Assumptions:
1. Main reaction has 98% conversion [Chosen logically since 98% selectivity of main
product is given]
3. As per literature, we normally purge 1% and recycle 99% as the raw materials form a
major cost (around 50%)
38
Step 1:
Now, using the above facts mass balance is carried out. As mentioned above, basis is chosen
and knowing the conversion of each of the reaction, mass balance is carried out on the
reactor. (All numbers are shown in the Excel file)
Step 2:
Now, the products of the reactor are to be separated into 2 streams viz., one which will give
the products and the one which needs to be recycled. This mass balance is done using
Flash. Since recycle is done a fraction of that is purged to avoid accumulation of the inerts
which are produced by the side reaction.
Step 3:
Now the mass balance is done on the first Distillation column. We use the fact stated in the
patent and do a linear mass balance. We are able to separate water as residue and send
purer azeotropic mixture of THF Water in the distillate.
Step 4:
The DC-102 distillation is operated at a higher pressure due to Pressure Swing Distillation
method which helps us in obtaining THF. In this distillation column, we get pure THF as
bottom product because THF-Water form a minimum boiling azeotrope and hence the
azeotrope goes in the upper product part of the DC-102. Now this is recycled to DC-101 and
a part of it is purged.
Overall mass needs to be conserved and it is being followed in the mass balance carried
out.
The molar basis chosen was 1:100 for Maleic Anhydride:H2 (as per patent).
65,050 65,050
39
9. Sensitivity Analysis
9.1 DC101
40
Variation in local parameters
We find out Water flow rate in bottom stream of DC-101 increases slightly on
increment in feed stage then becomes constant; this implies not much effect
of feed stage on flow rates and purity.
41
Fig 9.5 : Change in mole fraction of water in bottom of DC-101 on
changing input feed stage in DC-101
9.2 DC102
i) Variation of Purity and Flow rate of THF via bottom of DC 102 against feed
stage
42
Fig 9.7 : Change in mass flow rate of THF in Bottom of DC-102 on
changing input feed stage in DC-102
We have Feed stage as stage number 5. So, we have done a sensitivity analysis to check
what happens to the product (THF) flow if we vary feed stage from 3 to 6. We can look from
the above 2 graphs, that there is a meagre change in the purity and flow rate of THF on
changing the feed stage in the vicinity of the chosen feed stage. Hence, we can conclude that
feed stage is insensitive towards purity of THF.
43
10. Equipment Sizing & Costing
10.1 Reactor
Cost of Vessel
We are using a fixed bed reactor but, kinetics data is not available for our reactions so we
assumed it to be a vertical column and calculations are done accordingly. For material of
construction we are taking carbon steel because constituents of the process are
non-corrosive in nature.
𝞽 = Time constant for fixed bed reactor = 13.6 sec
3
Vflow =
Feed volumetric flow rate = 0.577 m /sec
V reactor
τ= V f low
We get, L = 3.41 m
D = 1.7 m
P = Design Pressure
44
where, t = Thickness of vessel
P = Design pressure
W = 0.25*[ 3 + (Rc/Rk)0.5 ]
After calculating the thickness, we add corrosion and milling tolerance to the thickness to
calculate trecommended.
t+ε
t recommended = m
1 − 100
After obtaining the thickness required for construction, we calculate the weight of material
required for the construction of vessel.
Vessel length (L) = 3.41 m
For costing of reactor vessel we have followed the method stated in the book21
C = F m *C b + C a $
21
Chemical Process Equipment - Selection and Design 2nd edition by James R. Couper W. Roy Penney James
R. Fair
45
C a = 300 * D 0.7396 * L 0.7066 $
where, Fm = cost factor associated with material of construction
Cb = cost accounting the weight of material used
Ca = cost accounting the platforms and ladders with the vessel
C = total cost of vessel
Ca = $ 1,690.92
Cb = $ 43,372.57
C = $ 45,064 (in 2003)
C = $ 75,891 (in 2008) = INR 53,12,335
Cost of Catalyst
In the fixed bed reactor, we are using a catalyst comprising of mixture of oxides of copper,
zinc, aluminium & chromium. The composition of constituents of the catalyst was taken
from the Budge et al22. We have assumed that 80% of the total volume of reactor is
occupied by the catalyst bed, 20% has been left for the inflow & outflow of reactant &
product streams. We have assumed a void fraction of 0.4 for the catalyst bed. Data
regarding catalyst life was not available so, we have assumed that 10 % of the total catalyst
weight will be replaced every year during shutdown period.
Volume of reactor (V) = 7.85 m3
= 3.76 m3
Table 10.1: Catalyst Composition
22
VAPOR-PHASE HYDROGENATION OF MALEIC ANHYDRIDE TO TETRAHYDROFURAN AND GAMMA-BUTYRO
-LACTONE, PATENT NUMBER 4,965,378
46
Average Density = [ ∑ Densityi*Percentagei ] /100 where i = 1,...,4 (Sr no.)
= 52966 kg/m3
catalyst (WT)
= 19842.28 kg
ith catalyst
47
10.2 Flash
Calculation for sizing and costing of flash vessel will be done as stated in the book23.
Next we use below mentioned graph to find Kv, the design velocity factor. This plot is for
85% flooding
Sv
Fig 10.1: Design Velocity Factor v/s Vapour-Liquid Separation Factor Curve24
Kv = 0.42
Using Kv, we calculate (Uv)max , then using (Uv)max and Qv we calculate Amin and thus Dmin .
ρ liq − ρ vap
(U v ) max = K v [ ρ vap
] 0.5
Qv
A min = (U v ) max
23
Chemical Process Equipment Selection and Design by Stanley M. Walas
24
Chemical Process Equipment Selection and Design by Stanley M. Walas
48
4*A min
D min = [ π
] 0.5
Qv = Volumetric flow rate = 91.53 m3
We get, L = 2.25 m
P = Design Pressure
After calculating the thickness, we add corrosion and milling tolerance to the thickness to
calculate trecommended.
t+ε
t recommended = m
1 − 100
where, 𝞊 = Corrosion allowance = 1.5 mm
49
m = Milling tolerance = 12.5%
trecommended = 14.54 mm
After obtaining the thickness required for construction, we calculate the weight of material
required for the construction of vessel.
Vessel length (L) = 2.25 m
For costing of reactor vessel we have followed the method stated in the book25
C = F m *C b + C a $
Ca = Cost accounting the platforms and ladders with the vessel
25
Chemical Process Equipment - Selection and Design 2nd edition by James R. Couper W. Roy Penney James
R. Fair
50
C = Total cost of vessel
Ca = $ 803.91
Cb = $ 17,614.16
Size Estimation:
Sizing of shell
Diameter calculation: As we know that flooding velocity decides limit on the max
allowable vapour flow to the column. So using Souders Brown equation we found Dmin for
rectifying and stripping section then took maximum diameter as esign diameter.
51
Fig 10.2: Graph for calculating Flooding Velocity26
Reflux 0.440
26
Chemical Process Equipment Selection and Design by Stanley M. Walas
52
Diameter (ft) 3.106 1.719
Number of stages 16 20
Number of trays 22 28
53
Cost Estimation
Cost of shell
a. We can calculate cost per height of column given column diameter from
graph given below:
Shell cost USD Height (in ft) Shell cost USD Shell cost INR
(per ft) (1987) (2018)
27
Bookmatter Chemical Engineering Economics
https://www.scribd.com/document/392007315/1989-Bookmatter-ChemicalEngineeringEconomics
54
Cost of Sieve Trays
Given column diameter we can calculate installed cost of tray and apply cost index to scale
cost upto 2018.
Cost of Reboiler
We have used graph (shown below) mentioned in lecture of Sieder. To read from graph all
you need is total heat transfer area, which is calculated using standard heat exchanger
calculation.
28
Plant design and economics for chemical engineers by Peters, M. S., K. D. Timmerhaus
55
Fig 10.5: Cost of shell & tube heat exchanger against area29
Cost of Condenser
Fixed and floating head readings are taken from above graph corresponding to heat
transfer area as calculated above.
29
Lec 06 Equipment Sizing and Capital Cost Estimation
https://www.seas.upenn.edu/~dlewin/CACHE_Workshop/LECTURE_06_Equipment_Sizing_and_Capital_Cost
_Estimation.pdf
56
Table 10.11: Condenser cost of DC-101 and DC-102
Process stream is at 41℃ & 8 bar pressure and needs to be heated upto 120℃ & 8 bar.
We are using MP steam at 175℃ & 10 bar to heat the process stream, it will cool down to
174℃ & 10 bar.
ΔTLMTD = 88.7 ℃
Q = 377568.11 Watt
A = Q / ( A * ΔT LM T D )
A = 7.5 m2
Area obtained from Aspen comes out to be 8.49 m2. Both the areas are reasonably close.
For costing of heat exchangers we have used the book30 as the reference performed
calculations as mentioned in the book.
C = 1.218 * f d * f m * f p * C b
30
Chemical Process Equipment - Selection and Design 2nd edition by James R. Couper W. Roy Penney James
R. Fair
57
where, C = Total cost of heat exchanger
fd = Cost factor which accounts for the type of heat exchanger used
A = area in sq feet
Type fd
CS/CS 1
58
Table 10.15: Heat exchanger double pipe pressure factor
≤4 1.00
4-6 1.10
6-7 1.25
From these calculations we obtain the cost of heat exchanger in USD in the year 2002.
Using Chemical Engineering Plant Cost Index (CEPCI) for the year 2002 and 2018 we obtain
the cost in USD in the year 2018 and than in INR in the year 2018.
59
10.5 Compressor
Compressors can be of two types- positive displacement compressors and dynamic
compressors. Dynamic compressors can be further divided as Axial and Centrifugal
compressors. Analysing our plant condition i.e. high flow rates and high pressure difference
we choose centrifugal as our type of compressors.
We calculate the work done by the formula -
W = P 1 * V 1(n/(n − 1)) [(P 2/P 1) (n−1/n) − 1]
But since the value of n depends upon the design and operation of the machine (we know
this from slides of course- Process equipment selection) and as our references did not
mention anything about value of ‘n’ we assume it to be equal to (Cp/Cv). We calculate the
work done as an isentropic process but to approximate towards actual answer we used
polytropic efficiency by using relation between polytropic and isentropic efficiency.
E c = (P 2/P 1) ( −1)/Y
− 1/ (P 2/P 1) ( −1)/Y Ep
−1
Ec = Isentropic efficiency
Ep = Polytropic efficiency
31
Plant design and economics for chemical engineers by Peters, M. S., K. D. Timmerhaus
60
Work required (W’) = W/ Efficiency
[Here P1, V1 stands for pressure and volume at the initial conditions and P2, V2
corresponds to the final conditions, in theoretical calculation we assume no change in
)
temperature and
Table 10.17: Power requirement of Compressor
Cost Estimation
For cost estimation, the cost of compressor is calculated as a relation of its power [10.1] -
Now since these cost refer to the year 2002 we need to use cost index to map it to the
current year and also convert it to Rupees from Dollars.
61
10.6 Pump
Analysing our flow rates and head requirements we chose centrifugal pumps from various
types of pumps available.
E = Pump efficiency
Table 10.19: Power requirement of Pumps
62
Cost Estimation
For cost estimation, the cost of pump is calculated as a relation of its flow rate and head
[10.1]-
Cost = Cb*FM*FT
FM = from the table provided in the book as per material used (Cast steel)
ˆ ˆ ˆ
Cb = 1.39 * e (8.833−0.6019*(ln(f l*((H) 0.5) + 0.0519*(ln(f l*(H) 0.5) 2))))
ˆ ˆ ˆ
FT = e (9.8849+ (−1.6164*(ln(f l*(H) 0.5))+0.0834*(ln(f l*(H) 0.5)) 2))
Where, fl = Flow rate in gpm (1 m3/s = 15850.32314 gpm), H = Head(ft.)
Factors FT, Cb are mentioned in [10.1] as part of cost estimation of multistage centrifugal
pump.
Now since these cost refer to the year 2002 we need to use cost index to map it to the
current year.
63
11. Plant Economics
11.1 Major Equipment Costing
The equipment costing methods employed are taken from the books Peters et. al. and
Walas et. al.
The Chemical Engineering Plant Cost Index(CEPCI) is used to calculate the cost to year
2018.
Flash 40,73,442
Reactor 1,44,07,305
Pump 89,19,910
Compressor 2,92,96,254
The working capital capital is calculated as 15% of the total capital investment (TCI). The
total capital investment (TCI) and fixed capital investment (FCI) estimates including direct
and indirect plant costs are calculated below.
64
Table 11.2: Total Capital Investment
Piping 15 10 5,08,25,126
Buildings 10 6 3,38,83,418
Land 2 1 67,76,684
Contingency 10 6 3,38,83418
65
11.3 Revenue Calculation
The revenue generated from the plant is calculated on year basis in the given Table 11.3
Utilities 74,11,167
Catalyst 3 3 15,53,76,991
66
maintenance
Patents 2 2 10,35,84,661
67
11.7 Payback Period
Payback period is defined as the time after which you recover your total fixed capital
investment.
The price of our final products as well as raw materials is considered to be constant for
profit calculations. A total of 90 percent of FCI excluding the land cost is considered
depreciable over a period of 10 years. A linear depreciation model is used to estimate the
average depreciation occurring for the plant.
Bank Interest = 6%
68
12. Process Control and Instrumentation
12.1 Distillation System
We will treat DC101 as an independent unit for process control as minor disturbances in
this unit can directly affect the quality of product produced. Therefore, we ensure the
quality of product in distillate of DC101.
We will later integrate the CO103 and F102 in process control to stabilize the feed to
DC102.
For Designing the control strategy for Distillation column the book “Distillation Design and
Control” by William Luyben is referred.
The degrees of freedom of a process system are the independent variables that must be
specified in order to define the process completely. [An approach was developed by
Waller, V. (1992)]
There are five control valves, one on each of the following streams: distillate, reflux,
coolant, bottoms and heating medium. The feed stream is considered being set by the
upstream process. In this way, the column will have five degrees of freedom. Inventories in
any process must be always controlled which involve liquid levels and pressures. This
means that the liquid level in the reflux drum, the liquid level in the column base, and the
column pressure must be controlled.
32
Degree of Freedom Analysis for Distillation Column
69
To find the DOF we count the total number streams involved and separately add up extra
phases i.e. add up all occurrences of phases greater than one in all units. The number of
control degrees of freedom is the difference between these two numbers.
egrees of Freedom = 5
From, Theory, D
So the number of degrees of freedom is 5. However, a typical control strategy for such a
process would use only 4 of these - feedrate, column pressure, top and bottom composition.
This is because the column and condenser are normally maintained at the same pressure.
Same is the case in or DC101 system.
Controlled Variable: The process variable which we want to maintain at particular value
Manipulated Variable: Process variable that is adjusted to bring the controlled variable
back to set point.
Disturbance Variable: Any process variable that can cause the controlled variable to
change.
http://homepages.ed.ac.uk/jwp/control06/controlcourse/restricted/course/fourth/course/module3-
1.html
70
12.1.2 Pressure Control
This is required to provide good response to column disturbances, and to protect the
column from disturbances occurring in the heating medium.
For inlet steam controlled reboiler, the heat transfer rate is regulated by varying the steam
control valve opening, thereby changing the steam condensing pressure and temperature.
33
For Pressure Control of Distillation Colum
http://www.separationprocesses.com/Operations/POT_Chp03.htm
71
Fig. 12.3: Reboiler control [12.2]
When an additional boil-up is required, the valve opens and raises the reboiler pressure,
which increase the temperature, and in turn increases the boil-up rate. This scheme has the
disadvantage of non-linear relationship between pressure and boil-up, and is affected by
fouling in the reboiler.
The other to be controlled is the composition of the tops product. The reason is that the
final product is directly linked to the top of DC101 hence it is important to control its
composition.
But change in concentration can be related to change in temperature inside the distillation
column only when a sufficient gradient is present. And in our system such a gradient is
available between tray 4-5-6.
72
Fig 12.4 Temperature Profile of DC101
We need to make sure that there is required head in distillation column is maintained.
Here, we can manipulate it with adjusting the feed rate or rate to release of water at
bottom.
The Level sensor of distillation column is then cascaded with flow sensor in the water exit
line and hence liquid level in distillation column can be maintained.
Variations in the feed stream are very possible due to upstream processes. We need to
make sure that feed to distillation is maintained at constant rate. Here, a simple flow
controller can do this task.
73
12.2 Control Strategy for F102 System
12.2.1 Cooler
Variations in the temperature of the feed stream to this cooler can be present due to
upstream processes. We need to make sure that the outlet comes at a constant temperature
of 40. Therefore we sense the temperature of outlet stream of Heat Exchanger and this is in
turn controlled by the rate of cooling water.
Variations in flow in the outlet are okay because F102 can cope with this changes.
We need to make sure that there is required head in flash is maintained. Here, we can
manipulate it with adjusting the feed rate or rate to release of fluid at bottom.
The Level sensor of distillation column is then cascaded with flow sensor in the water exit
line and hence liquid level in distillation column can be maintained.
We need to make sure that operates at a constant pressure. The top of the flash is a purge
stream to remove the minimal hydrogen impurity on this loop.
The pressure in this vessel is adjusted by controlling the outlet of this line.
A pump is very sensitive to flow to it. If not properly maintained can shut down.
Therefore a controller to this is very important.
In this type of control scheme if flow to pump reduces the valve opens and liquid start
flowing in the pump loop. When a considerable flow is achieved it is noted by the flow
sensor and and then the valve can close as the pump is now safe from cavitation.
74
12.3 Overall Control Diagram
75
13. Environmental Impact Assessment
In theory, our plant is not emitting any harmful gases in the atmosphere. In our flowsheet,
we have 7 boundary streams. Two of these boundary streams are input streams from
which maleic anhydride & hydrogen enter the system. Remaining 5 streams include purge
streams from various equipments and stream from the bottom of distillation columns.
Following table gives an account of what is being emitted in the environment by boundary
streams.
Amount
Emissions Process
(kg/day)
76
13.1 Emissions from Boundary Streams
13.1.1 Emissions from Stream number 31
Stream number 31 is purge stream of flash F102. This stream is at 40℃ & 1 bar pressure.
34
[13.1] MSDS Butane
http://www.afrox.co.za/en/images/Butane_tcm266-27624.pdf
[13.2] MSDS Butanol
https://www.cdhfinechemical.com/images/product/msds/19_1391926836_n-Butanol-CASNO-71-36-3-MSD
S.pdf
77
13.1.2 Emissions from Stream number 28
Stream number 28 is purge stream of separator S1. This stream is at 142℃ & 25 bar
pressure.
Table 13.3: Emission from stream no. 28
Hydrogen (H2): Majority of H2 is being released from S1 separator. The temperature of this
stream is 142℃ and pressure is 25 bar. The major components are H2, THF, water, butane
and CO2. All of these will be in liquid form except H2. So, we first we will flash this stream
to separate H2 and remaining components. This separated H2 can be sent back to the input.
In mixture of remaining components, we will reduce the pressure so that condensables can
be removed and THF and water can be fed into the DC101.
79
13.1.4 Emissions from Stream number 20
Stream number 20 is bottom stream of distillation column DC102. This stream is at
140.95℃ & 8 bar pressure.
Table 13.5: Emission from stream no. 20
This stream is our main product stream. THF mass fraction is 0.98
80
13.1.5 Emissions from Stream number 14
Stream number 14 is bottom stream of DC101 . This stream is at 97.04℃ & 1 bar pressure.
Maleic anhydride: Majority of maleic anhydride emission is coming from the bottom of
DC101 which is about 2.89 tons/day, which is a big number considering the fact it is one of
the raw material and is very expensive. So, we may want to recycle it. In order to recycle it
we need to remove water & butanol which is present in this stream along with maleic
anhydride. Butanol boils at about 117℃ and water at 100℃. Butanol and water form an
azeotrope which boils at about 90 ℃, also maleic anhydride boils at 200℃ . What we can
do here is to raise the temperature of this stream upto 150-170℃ and flash it. This will
separate all the components in the stream from maleic anhydride. Maleic anhydride will
leave the flash vessel in liquid phase and rest of the components in vapour phase.
Else we can do is that the temperature of this stream will be dropped to 50 degree, as
melting point of Maleic Anhydride is 65 degree, its state will change to solid form.
This solid can be recovered by filtration and rest can be emitted out.
81
14. Heat Integration
14.1 Streams present in plant
14.1.1 Segregation to Sensible and Latent Heat part
Analysing the streams in the plant we found - 6 Cold stream and 4 Hot stream details of
which can be found in table below. FCph [W/C] and FCpc [W/C] represents ratio of change
in enthalpy to temperature over than enthalpy range.
Change in enthalpy values has been segregated into sensible and Latent heat part. To find
FCph and FCpc, we took out the latent heat part from total enthalpy given by aspen then
took ratio of remaining enthalpy with temperature change as shown in sample calculation
below.
82
5. FCpc = Sensible heat / delta_T
= 7925.11 [W/C]
While below the pinch temperature, we have integrated E107 with E101, E104, E105 as
FCph > FCpc, table is attached for below pinch point temperature, above pinch point stream
E102, E103, E106 and E107 (till pinch temperature) has to be cooled using utilities.
Table 14.2: Stream listing for heat exchanger network below pinch point
35
Table 14.3: Stream listing for heat exchanger network above pinch point
35
The green colour in Hot stream section shows that heat is available in this section and amount is mentioned
in Watts. The red colour in Cold stream shows that heat is required and amount is mentioned in Watts
83
We have used problem table method36 for pinch analysis, we found pinch temperature to
be 136 C.
Now, for E 110 (product stream) we would use cooling water which is discharged from
E109, E108 or E107 stream, in this way cost for cooling product would be negligible.
36
Harvard (18th ed.) - COULSON, J. M., RICHARDSON, J. F., & COULSON, J. M. (2002). Coulson & Richardson's
chemical engineering. Oxford, Butterworth-Heinemann.
84
a. Operating Cost before heat integration = INR 11,70,41,789
b. Operating Cost after heat integration = INR 8,16,02,985
c. Savings = INR 3,54,38,804
Blue lines represent hot streams which are cooled with the means of heat exchange while
red line represent the cold streams which are heated via same heat exchange.
HY - Heater; CY - Cooler
85