CPD3297 Appendix
CPD3297 Appendix
Table of Content
Number Title Page
APPENDIX A A-1
Appendix A.1 Stoichiometry and Reaction Kinetics A-2
Appendix A.2 Thermodynamic data A-6
Appendix A.3 Input output streams and properties A-8
APPENDIX B B-1
Appendix B Binary system comparison B-2
APPENDIX C C-1
Appendix C.1 Selection of operating temperature C-2
Appendix C.2 Criteria and Selections for CO2 Removal C-4
Appendix C.3 Heat Integration C-10
Appendix C.4 Options and selections of Heat exchanger C-18
Appendix C.5 Propane and propylene separation C-23
Appendix C.6 Comparison of the Tray and Packed column properties C-29
Appendix C.7 Recommendations for treatment of light gas C-30
APPENDIX D D-1
Appendix D. Process Stream Summary D-2
APPENDIX E E-1
Appendix E.1 Aspen Plus simulation results E-2
Appendix E.2 Simulation results of distillation columns T301 T302 and
T303 E-7
Appendix E.3 Reactor Design E-34
Appendix E.4 Shell and Tube Heat Exchanger Design E-45
Appendix E.5 T302 Column sizing report E-59
Appendix E.6 Calculation of CO2 removal equipment E-60
Appendix E.7 Gas-liquid separators calculation E-69
Appendix E.8 Equipment Summary & Specification Sheets E-71
APPENDIX F F-1
Appendix F Process safety F-2
APPENDIX G G-1
Appendix G Economics G-2
APPENDIX H H-1
Appendix H.1 PFS for process with heat integration H-2
Appendix H.2 Summary of utilities for the process with Heat integration H-5
Appendix H.3 Heat and Mass balance for the process with Heat
integration H-9
Appendix H.4 Economy H-20
i
APPENDIX A
A-1
Appendix A.1 Stoichiometry and Reaction Kinetics
a) Propane oxidative dehydrogenation, exothermic
1
C3 H 8( g ) + O2( g ) → C3 H 6( g ) + H 2O( g ) (A.1.1)
2
C3 H 8 → CO, CO2 , CH 4 and H 2O
Mechanism:
C 3 H 8 + [O ] → C 3 H 7 * + [OH ] (A.1.2)
The radical mechanism of reaction can be indicated as below
Initiation:
O 2 → O 2∗ + O 2∗
CH 3 CH 2 CH 3 → CH 3 CH 2 CH 2 ∗ + H ∗
CH 3 CH 2 CH 3 → C H 3 CH 2 ∗ + CH 3 ∗
Propagation:
∗
O 2∗ + CH 3CH 2CH 3 → OH * + CH 3CH 2CH 2
∗
OH * + CH 3CH 2CH 3 → H 2O + CH 3CH 2CH 2
∗
CH 3CH 2CH 2 → H 3CHC = CH 2 + H ∗
∗ ∗
CH 3CH 2CH 2 + CH 3CH 2CH 3 → H 3CHC = CH 2 + CH 3CH 2CH 2 + H 2
∗
CH 3CH 2 → H 2C = CH 2 + H ∗
∗ ∗
CH 3CH 2 + CH 3CH 2CH 3 → H 2C = CH 2 + CH 3CH 2CH 2 + H 2
∗
CH 3 + CH 3CH 2CH 3 → CH 4 + H 2C = CH 2
∗
H ∗ + CH 3CH 2CH 3 → CH 3CH 2CH 2 + H 2
etc.
A-2
Termination:
H ∗ + OH * → H 2O
OH * + CH 3CH 2CH 2∗ → H 2O + H 3CHC = CH 2
OH * + CH 3CH 2∗ → H 2O + H 2C = CH 2
H ∗ + H ∗ → H2
H ∗ + CH 3CH 2CH 2∗ → H 2 + H 3CHC = CH 2
H ∗ + CH 3CH 2∗ → H 2 + H 2C = CH 2
H ∗ + CH 3∗ → CH 4
CH 3CH 2∗ + CH 3∗ → CH 4 + H 2C = CH 2
CH 3CH 2∗ + → H 2C = CH 2 + CH 3CH 2CH 3
CH 3∗ + CH 3CH 2CH 2∗ → CH 4 + H 3CHC = CH 2
etc.
2 [OH ] + 0.5O2 → H 2O + [O ] (A.1.3)
1
C3 H 8( g ) + O2( g ) → C3 H 6( g ) + H 2O( g ) ; ∆H = −28.3kcal / mol (A.1.6)
2
1 1
C3 H 8( g ) + O2( g ) → C2 H 4( g ) + H 2O( g ) ; ∆H = −33.5kcal / mol (A.2.9)
2 2
1
CO( g ) + 2 H 2( g ) → C2 H 4( g ) + H 2O( g ) ; ∆H = −25.1kcal / mol (A.2.10)
2
A-3
60
50
C3H8 Selectivity (mol%)
40
C2H4
30 C3H6
CH4
20
COx
10
0
650 700 750 800 850 900
T (oC)
Figure A.1.1. Effects of temperature on the COx, CH4, C2H4, and C3H6 selectivity in oxidative
conversion of propane (VH Rane et al, 2003)
From Figure A.1.1, the favorable temperature for reaction is at 850oC See also
Appendix C.1. At this condition, the selectivity of the outlet stream for this
experiment is shown as 50%C2H4, 25%C3H8, 15%CH4, and 10% COX respectively.
Then, we can back calculate to obtain the stoichiometry of oxidative dehydrogenation
as eq.A.2.11. This reaction will take place in the tube side of this design project.
C3 H 8( g ) + 0.4759O2 ( g ) → 0.6946C2 H 4 ( g ) + 0.3053CH 4 ( g ) + 0.3607C3 H 6 ( g ) + 0.2236 H 2 ( g )
+ 0.6945H 2O( g ) + 0.1894CO( g ) + 0.0339CO2 ( g ) (A.2.11)
; ∆H = −20 kcal / mol
The experimental data for reaction rate is supported by kinetic modeling using
Langmuir Huinshelwood (LHHW).
k LH K O2 K C3 H8 PO2 PC3 H8
r= (A.2.12)
(1 + K O2 PO2 )(1 + K C3 H8 PC3 H8 )
Stochiometry
For the shell side, the reaction takes place only propane dehydrogenation, that we can
combine the reactions in eq. A.2.13 and eq.A.2.14. and result in eq. A.2.15.
A-4
C3 H 8( g ) → 0.002C2 H 4( g ) + 0.002CH 4( g ) + 0.998C3 H 6( g ) + 0.998H 2( g ) (A.2.15)
Dehydrogenation reaction: C3 H 8( g ) ⇔ C3 H 6( g ) + H 2( g )
⎛ PC H PH ⎞
k1 ⎜ PC3 H8 − 3 6 2 ⎟
(− r1 ) = ⎝
K1 ⎠
(A.2.16)
1 + PC3 H 6 K 2
⎛ − Ea1 ⎛ 1 1 ⎞ ⎞
k1 = k01 exp ⎜ ⎜ − ⎟⎟ (A.2.17)
⎝ R ⎝ T Tm ⎠ ⎠
⎛ −∆H ⎛ 1 1 ⎞ ⎞
K 2 = K 02 exp ⎜ ⎜ − ⎟⎟ (A.2.18)
⎝ R ⎝ T Tm ⎠ ⎠
Cracking reaction: C3 H 8( g ) ⇔ CH 4( g ) + C2 H 4( g )
(− r2 ) = k4 PC3 H8 (A.2.19)
⎛ − Ea4 ⎛ 1 1 ⎞ ⎞
k4 = k04 exp ⎜ ⎜ − ⎟⎟ (A.2.20)
⎝ R ⎝ T Tm ⎠ ⎠
A-5
Appendix A.2 Thermodynamic data
Table A.2a. List of thermodynamics heat data-- Gibbs Energy
A-6
Table A.2b. List of thermodynamics heat data—Saturated liquid density
Design
Formula CAS-Nr. [g/ml]
Systematic
A B n Tc
Propane C3H8 74-98-6 0.22151 0.27744 0.28700 369.82
Propylene C3H6 115-07-1 0.23314 0.27517 0.30246 364.76
Ethylene C2H4 74-85-1 0.21428 0.28061 0.28571 282.36
Methane CH4 74-82-8 0.15998 0.28810 0.27700 190.58
Hydrogen H2 1333-74-0 0.03125 0.3473 0.27560 33.18
Oxygen O2 7784-44-7 0.43533 0.28772 0.29240 154.58
Water H2O 7732-18-5 0.34710 0.27400 0.28571 647.13
Carbon dioxide CO2 124-38-9 0.46382 0.26160 0.29030 304.19
Carbon monoxide CO 630-08-0 0.29818 0.27655 0.29053 132.92
Methyl CH3N(CH2
CH2OH)2 105-59-9 n/a n/a n/a n/a
Diethanolamine
A,B,n=regression coefficient for chemical components
T=Temperature, K
Tc=critical temperature, K
A-7
Appendix A.3 Input output streams and properties
Feedstocks
a. Propane
b. Oxygen
Products
a. Ethylene
A-8
Heavy Ends %wt <5 0.0
Total 100.0
Process Conditions and Price
Temp. K 303
Press. Bara 15
Phase V/L/S L
Price US$/ton 450 (7)
b. Propylene
d. Carbon dioxide
A-9
Available Design note numbers)
Carbon dioxide %wt >95 100.0
Heavy Ends %wt <5 0.0
Total 100.0
Process Conditions and Price
Temp. K 298
Press. Bara 4.5
Phase V/L/S L
Price US$/ton 6.5
A-10
APPENDIX B
B-1
Appendix B Binary system comparison
Keep the temperature constant, to find different pressure according to difference liquid
mole fraction. Compare the result from experimental data1 and simulation by Aspen. The
results are listed as bellow.
Table B.1 The comparison between experiment data and results from ASPEN
Experimental data ASPEN
Number T
x1 P x1 p
1 0.2053 3.447 0.2 3.356
260.93
2 0.5557 3.731 0.55 3.634
3 0.472 4.826 0.475 4.725
269.54
4 0.857 5.185 0.85 5.106
5 0.6554 6.322 0.65 6.253
277.59
6 0.8454 6.536 0.85 6.505
7 0.799 12.169 0.8 12.171
301.32
8 0.9 12.348 0.9 12.395
9 0.1794 13.604 0.175 13.622
310.93
10 0.7904 15.241 0.8 15.309
11 0.756 23.194 0.75 23.152
330.32
12 0.874 23.594 0.875 23.654
13 0.1769 27.49 0.175 27.614
344.26
14 0.7901 30.569 0.8 30.736
35
30
Pressure [bar]
25
20
15
10
5
0
1 2 3 4 5 6 7 8 9 10 11 12 13 14
liquid fraction literature P of literature [bar] liquid fraction ASPEN P of ASPEN [bar]
Figure B.1 The results from ASPEN compared with experimental data
1
Reference: H.Knapp, Vapor-liquid Equilibria for mixtures of low boiling substances, Chemistry data
series.
B-2
Figure B.1 shows that at the same liquid mole fraction, the vapor pressure of the binary
system from Aspen is almost the same with experimental data at the same temperature.
These verify the correct of Aspen simulation prediction.
B-3
APPENDIX C
E-1
Appendix C.1 Selection of operating temperature
In the design stage, operating temperature of reaction is imperative to be determined,
because it influences to the product selectivity, and then the margins of the process. In
view of exergy loss reduction, operating temperature also impact the loss work or loss
energy from the temperature difference, thereby the operating temperature should be
estimated to optimum both economic point of view and exergy loss reduction.
In this project, the reactions of both, endothermic and exothermic, take place in shell
and tube reactor respectively. The temperature of tube side should be higher than that
of the shell side. From the shell side, the reaction takes place at the temperature 540oC
in order to achieve the high selectivity and high conversion. So, we will base on this
data for the shell side and figure the suitable operating temperature in tube side
following data.
From the data in Appendix A.1, Figure A.1.1, show the relationship between product
selectivity and temperature. When temperature is changed, the production selectivity
will also change. Consequently, the results of reactants and products in chemical
reaction estimated will change also. Therefore, to estimate the best operating
temperature, the calculations of 3 temperature levels, which are 750oC, 800oC, and
850oC indicated as the equation below.
In fact, the process can operate with the recycle stream of unconverted propane, C3H8.
So, to compare the whole results of these three operating temperature. The iterative
calculation is necessary to figure it out. The results to achieve the design product and
compatible with the data from literature show in the Table C.1.1.
Table C.1.1 The results of iterative calculation for different operating temperature in tube side
E-2
Operating Input Output
Temperature (kton/a) (kton/a)
(oC) C3H8 O2 C2H4 C3H6 CH4 H2 H2O CO CO2
750 240.87 63.78 81.51 118.5 20.48 3.35 52.36 22.1 6.22
800 236.98 102.27 41.26 158.74 3.61 3.76 80.34 42.31 8.97
850 222.18 40.6 5.72 194.29 1.36 7.46 28.3 20.24 5.22
And work loss that related to exergy loss mentioned in Chapter 8, eq.8.3.1.15 can be
figured out. Consequently, the margins and work loss of different operating
temperature for tube side can be estimated as the Table C.1.3
Table C.1.3. The margins and Wlost of difference operating temperature.
The re lation of margins and Wlost with diffe re nt ope rating te mperature
Margin (US$/a)
Wlost (kcal)
35000 50
45
30000
40
25000 35
Margins
30
20000
Wlost
25
15000
20
10000 15
10
5000
5
0 0
740 760 780 800 820 840 860
E-3
Appendix C.2 Criteria and Selections for CO2 Removal
Carbon dioxide is green house gas, and will cause global warming. It is not allowed
that High concentration of carbon dioxide present in product. So CO2 produced in
process should be removed.
Solvent chosen
In practice, physical and chemical absorption are both used to remove CO2.
According to higher separation efficiency of chemical absorption, it is applied in
design. CO2 is acid gas and should use amine to absorb it. Normally, methyl
ethanolamine (MEA), diethanolamine (DEA) and methyl diethanolamine (MDEA)
are chosen. Compare these three amines, MDEA has higher energy efficiency, greater
acid gas removal capacity, higher resistance to degradation, smaller equipment size
for the new plants and above all much less corrosivity as compared to primary and
secondary amines. As a result of following advantages, MDEA is the best choice to
remove CO2 in this process.
E-4
experience and laboratory data, relative corrosivity of amines are ranked as
follows: MEA >> DEA >> MDEA. Table C.2.4 generates corrosion data for
various amine-based solutions.
• MDEA is tertiary amine and therefore carbamate formation with CO2 does not
take place in MDEA based system. MEA and DEA form carbamates with
CO2. Therefore operation with MDEA is far more stable with no spurious
shutdowns over longer periods.
CONCENTRATION % 15 30 35-50
*Selectivity is defined as ratio of (mole percent of H2S removed to mole percent of H2S in feed gas)
to (mole percent of CO2 removed to mole percent of CO2 in feed gas)
30% Wt MEA 32
E-5
50% Wt DEA 25
15% Wt MEA 13
20% Wt DEA 8
50% Wt MDEA 3
E-6
Packed column chosen
Due to one liquid input (MDEA solution) and one vapor input (raw product stream) a
vertical absorption column is preferred which can get high exchange efficiency and
need less volume. Followed the criteria below, packed column absorber is chosen and
material of equipment is carbon steel.
Here, MDEA is corrosive for equipment and much pressure drop is not good for the
separations followed. As a result, packed column is chosen. Figure C.2.1 shows the
scheme of packed absorption column. Actually structured packing has low HETP
(typically less than 0.5) and low-pressure drop (around 100Pa/m); however, the cost
of structured packing per cubic meter is significantly higher than that of random
E-7
packing. So, random packed column is preferred. Assume a 99 per cent recovery of
the carbon dioxide is reached. In order to improve the liquid distribution
characteristics the type of packing is Pall ring, which increases the free area. Ring
packings are available in a variety of materials: ceramics, metals, plastics and carbon.
Metal and plastics (polypropylene) rings are more efficient than ceramic rings, as it is
possible to male the walls thinner. Due to MDEA a kind of amine that is corrosive to
plastics, metal (carbon steel) is used in design.
In general, the largest size of packing that is suitable for the size of column should be
used, up to 50mm. The reason is that small sizes are appreciably more expensive than
the larger sized. However, above 50mm the lower cost pre cubic meter does not
normally compensate for the lower mass transfer efficiency. If packing size is too
large in a small column it can cause poor liquid distribution. Recommended size
ranges are: [Couson, volume 6]
According to the stream flow rate into the unit 700m3/h approximate, the large range
of column diameter and related range of use packing size are chosen.
And Table C.2.6 [Distillation principles and practices] lists typical values of specific area and
porosity for several random column packings. Herewith related data is shown in Table
C.2.6.
E-8
Particle Diameter dn, Specific Area a Porosity ε
mm m2/m3
Pall rings 10 515 0.92
20 360 0.93
25 215 0.94
35 145 0.94
50 105 0.95
80 78 0.96
In order to use the largest packing size suitable for the size of column and get higher
specific area, 50mm Pall rings are used. And Onda’s method is based on a large
amount of data on gas absorption and distillation; with a variety of packings, which
included Pall rings. Accordingly, it is used for the calculation of column design.
E-9
Appendix C.3 Heat Integration
The method which is used in MER targeting is the Temperature Interval (TI) method.
The temperature-interval method was applied according to Linnhoff (1987).
After MER targeting, HEN is designed with a unit-by-unit method beginning at the
closest-approach temperature difference (the pinch analysis) (Linnhoff et al, 1987).
For the minimum utility requirements aver all possible HENs, the minimum approach
Step by steps of MER targeting for heat integration, which has been followed
according the Temperature Interval (TI) method (Linnhoff, 1987). And Seader (2003)
are presented shortly. First, “Hot and Cold Composite” Curve is constructed in one
graph. Then, “Hot and Composite” Curve plus ‘the pinch point’ is constructed. After
that, calculation for Tpinch by making ‘Cascade table’ is described as follows.
The first step in pinch analysis is to determine the pinch point by using Temperature
Interval (TI) method. An interval is an imaginary boundary to make our works easier.
We defined a 10oC for the minimum temperature difference, where interval is
1 1
∆Tmin (5 oC ) ∆Tmin (5 oC )
2 below the hot stream and 2 above cold stream.
The procedure of TI method is usually applied for the streams with constant the heat-
capacity flowrate (FCp). Meanwhile, in the alkenes plant; there are some phase
changes involved. Phase changes (latent heats) are counted into TI method formalism
simply by assuming 1oC temperature change at the temperature of the phase change;
FC p (1) = F ∆H v
E-10
Thus, phase changes simply increase the number of temperature interval considered
(Douglas, 1988).
In the Alkenes plant, heat integration is applied for pre-condition and reaction section.
The four hot steams and three cold streams in precondition and reaction systems can
be found in Table C.2.1.
The heating and cooling requirement in separation section is not included in heat
integration task. The temperature of condenser in separation section is very low due to
cryogenic distillation. Therefore, special design methods to reduce energy loss and
utilities cost is applied in the separation section as described below:
1. The light gas column T301.
The Coolant used in condenser E301B is hydrogen expanded, which obtain
from the F-T process, from 20 bar to 5 bar. This is because the overhead gas
temperature is very cold, around –130 C. and expanded H2 can be used as
refrigerant. The temperature of H2 expanded is around –250 C.
2. The ethylene column; T302.
Overhead stream as <306> is condensed by E302B, and reflux stream as
stream<307> can be used for purity adjust. The condenser of this column uses
part of ethylene product as refrigerant itself, by passing through expansion
valve. This is because at the top column, the temperature is around –40 C.
Ethylene after expansion valve will be reach to –70 C that can be used to
condense the overhead gas.
3. The last distillation column, T303.
This overhead gas is condensed by E303B, which use part of propylene
product as heat pump process referring to Appendix C.5. After heat pump,
propylene will be used as the heating media for reboiler, E303A.
E-11
Table C.3.1 The stream conditions and properties for precondition amd reaction section
Cp (J/kg- Temp(in)
No Stream Hot/Cold F (kg/hr) F (kg/s) K) (oC) Temp(out)(oC)FCp (kJ/oC/s)
1 102 Hot 39022.39 10.8 2637.2 540 25 28.59
3 205 Hot 32755.24 9.1 2761.2 850 25 25.12
4 104 Hot 39022.39 10.8 2736.9 227 55 29.67
5 208 Hot 27875.32 7.7 2307.3 308 30 17.87
13 001 Cold 29644 8.2 3509.6 25 540 28.90
6 312 Cold 35392.96 9.8 2805.9 43 540 27.59
2 203 Cold 32756.43 9.1 3254.7 540 850 29.61
After making the cascade of temperature intervals, which is shown in Figure C.3.1
and Enthalpy Differences for Temperature Intervals in Table C.3.2, the pinch point is
acquired quickly.
855oC
845oC
850oC 850oC
545oC
540oC
308oC
o
303 C
227oC
222oC
50oC
55oC
48oC
o
43oC
o o o
25oC 25 C 25 C 30 C 25 C
20oC
E-12
Table C.3.2 Enthalpy Differences for Temperature Intervals
From
hot utilities
Interval
855oC i
Ti − Ti −1 ∑ FC P,hot − ∑ FCP,cold ∆H (kJ / s ) + or - Accumulate
2602
E-13
From Hot Utility Cumulative From Hot Utility Cumulative Hot
Utility
855oC 0 kW +2602 kW
∆H = 0
∆H = −296
845oC - +2305kW
∆H = −1347 ∆H = −1347
o
535 C -1957kW +654kW
∆H = −644 ∆H = −644
o
303 C -2602kW 0 kW
∆H = +1222 ∆H = +1222
o
222 C -1379kW +1222kW
∆H = +7698 ∆H = +7698
o
50 C +6319kW +8920kW
∆H = +30 ∆H = +30
48oC +6349kW +8950kW
∆H = +768 ∆H = +768
o
25 C +7117kW +9718kW
∆H = +269 ∆H = +269
+9987kW
o
20 C +7385kW
∆H = 0 ∆H = 0
+7385kW +9987kW
To cold Utility
To cold Utility
From the Cascade table above, we found that Tpinch is at Tinterval = 303oC. In the other
word, Tpinch = 308oC for Hot Streams (5oC above Tinterval)
Tpinch = 298oC for Cold Streams (5oC below Tinterval).
E-14
z The Utility Consumption for Heating and Cooling
E-15
Figure C.3.3 Pinch Design (Heat Exchanger Network)
The list of heat exchangers after heat integration can be found in Chapter 11,
economic part and in equipment specification (Appendix E.8).
E-16
is applied, the usage of hot and cold utilities decreases significantly. However, the
heat integration has an impact in the investment cost. Most of liquids in the Alkenes
plant are in gaseous phase, which have very low heat transfer coefficient (h). Typical
overall heat transfer coefficient (U) for gas-gas heat exchanger is 10-50 (W/m2.oC).
Consequently, the heat exchange area becomes very large and will results the very
expensive heat exchanger investment costs.
After conducting Heat integration procedure and making the new flowsheet, the CPD
project goes to economic assessment. The economic evaluation or profitability
analysis shows that the investment for the process with heat integration has significant
increasing due to the cost of heat exchanger with large heat exchange area. Some
large heat exchangers require very expensive purchased cost and result a high
investment cost.
In this CPD report, two flowsheets are presented. One flowsheet is the process ‘with’
heat integration and the other flowsheet is the process ‘without’ heat integration. Both
flowsheets; the process ‘with’ and ‘without’ heat integration are presented together
with mass and heat balance and the economic evaluation by profitably analysis.
Finally, after the economic evaluation, the profitability analysis shows that the
alkenes process ‘without’ heat integration is more profitable and more robust than the
alkenes process ‘with’ heat integration due to the high investment in heat exchanger
equipments.
E-17
Appendix C.4 Options and selections of Heat exchanger
Selection is the process in which the designer selects a particular type of heat
exchanger for a given application from a variety of heat exchangers. There are a
number of alternatives for selecting heat transfer equipment, but only one among
them is the best for the given set of conditions.
Selection Criteria
Selection criteria are many, but primary criteria are type of fluids to be handled,
operating pressures and temperatures, heat duty, and cost. Fluids involved in heat
transfer can be characterized by temperature, pressure, phase, physical properties,
toxicity, corrosivity, and fouling tendency. Operating conditions for heat exchangers
vary over a very wide range, and a broad spectrum of demands is imposed for their
design and performance. All of these must be considered when assessing the type of
unit to be used. When selecting a heat exchanger for a given duty, the following points
must be considered:
• Materials of construction
• Operating pressure and temperature, temperature program, and temperature
driving force
• Flow rates
• Flow arrangements
• Performance parameters i.e. thermal effectiveness and pressure drops
• Fouling tendency
• Types and phases of fluids
• Maintenance, inspection, cleaning, extension, and repair possibilities
• Overall economy
• Fabrication techniques
• Intended applications
Materials of Construction
For reliable and continuous use, the construction materials for pressure vessels and
heat exchangers should have a well-defined corrosion rate in the service environments.
Furthermore, the material should exhibit strength to withstand the operating
temperature and pressure. Shell and tube heat exchangers can be manufactured in
virtually any materials that may be required for corrosion resistance, e.g., from
nonmetals like glass, Teflon, and graphite to exotic metals like titanium, zirconium,
tantalum, etc. compact heat exchangers with extended surfaces are mostly
manufactured from any metal that has drawability, formability, and malleability. Heat
exchanger types like plate heat exchangers normally require a material that can be
pressed or welded.
E-18
The design pressure is important to determine the thickness of the pressure retaining
components. The higher the pressure, the greater will be the required thickness of the
pressure-retaining membranes and the more advantage there is to placing the high-
pressure fluid on the tubeside. The pressure level of the fluids has a significant effect
on the type of unit selected.
• At low pressures, the vapor-phase volumetric flow rate is high and low
allowable pressure drops may require a design that maximizes the area
available for flow, such as crossflow or split flow with multiple nozzles.
• At high pressures, the vapor-phase volumetric flow rates are lower and
allowable pressure drops are greater. These lead to more compact units.
• In general, higher heat transfer rates are obtained by placing the low-pressure
gas on the outside of tubular surfaces.
• Operating pressures of the gasketed plate heat exchangers and spiral plate heat
exchangers are limited because of the difficulty in pressing the required plate
thickness, and by the gasketed materials in the case of PHEs. The floating
nature of floating-head shell and tube heat exchangers and lamella heat
exchangers limits the operating pressure.
Temperature
Design Temperature
This parameter is important as it indicates whether the material at the desired
temperature can withstand the operating pressure and various loads imposed on the
components. For low-temperature and cryogenic applications toughness is prime
requirement, and for high-temperature applications the material has to exhibit creep
resistance.
Temperature Program
Temperature program in both a single pass and multipass shell and tube heat
exchanger decides (1) the mean metal temperatures of various components like shell,
tube bundle, and tubesheet, and (2) the possibility of temperature cross. The mean
metal temperatures affect the integrity and capability of heat exchangers and thermal
stresses induced in various components.
Temperature Driving Force
The effective temperature driving force is a measure of the actual potential for heat
transfer that exists at the design conditions. With a counterflow arrangement, the
effective temperature difference is defined by the log mean temperature difference
(LMTD). For flow arrangements other than counterflow arrangement, the LMTD
must be corrected by a correction factor, F. The F factor can be determined
analytically for each flow arrangement but is usually presented graphically in terms of
the thermal effectiveness P and the heat capacity ratio R for each flow arrangement.
E-19
Shell and tube heat exchanger units can be designed for almost any combination of
pressure and temperature. In extreme cases, high pressure may impose limitations by
fabrication problems associated with material thickness, and by the weight of the
finished unit. Differential thermal expansion under steady conditions can induce
severe thermal stresses either in the tube bundle or in the shell. Damage due to flow-
induced vibration on the shellside is well known. In heat exchanger applications
where high heat transfer effectiveness is required, the standard shell and tube design
may require a very large amount of heat transfer surface. Depeding o the fluids and
operating conditions, other types of heat exchanger design should be investigated.
Compact Heat Exchanger
Compact heat exchangers are constructed from thinner materials, which are
manufactured by mechanical bonding, soldering, brazing, welding, and etc. Therefore,
they are limited in operating pressures and temperatures.
Gasketed Plate Heat Exchanger and Spiral Exchanger
Gasketed plate heat exchanger and spiral exchanger are limited by pressure and
temperature, wherein the limitations are imposed by the capability of the gaskets.
Flow Rate
Flow rate determines the flow area: the higher the flow rate, the higher will be the
crossflow area. Higher flow area is required to limit the flow velocity through he
conduits and flow passages, and the higher velocity is limited by pressure drop,
impingement, erosion, and, in the case of shell and tube exchanger, by shell-side floe-
induced vibration. Sometimes a minimum flow velocity is necessary to improve the
heat transfer, to eliminate stagnant area, and to minimize fouling.
Flow Arrangement
As defined before, the choice of a particular flow arrangement is dependent upon the
required exchanger effectiveness, exchanger construction type, upstream and
downstream ducting, package envelope, and other design criteria.
E-20
Fouling Tendencies
Fouling is defined as the formation on heat exchanger surfaces of undesirable deposits
that impede the heat transfer and increase the resistance to fluid flow, resulting in
higher pressure drop. The growth of these deposits causes the thermohydraulic
performance of heat exchanger to decline with time. Fouling affects the energy
consumption of industrial processes, and it also decides the amount of extra material
required to provide extra heat transfer surface to compensate for the effects of fouling.
Compact heat exchangers are generally preferred for nonfouling applications. In a
shell and tube unit the fluid with more fouling tendencies should be put on the tube
side for ease of cleaning. On the shellside with cross baffles, it is sometimes difficult
to achieve a good flow distribution if the baffle cut is either too high or too low.
Stagnation in any region of low velocity behind the baffles is difficult to avoid if he
baffles are cut more than about 20-25%. Plate heat exchangers and spiral plate
exchangers are better chosen for fouling services. The flow pattern in plate heat
exchanger induces turbulence even at comparable low velocities; in the spiral units,
the scrubbing action of fluids on the curved surfaces minimizes fouling.
E-21
Overall Economy
There are two major costs to consider in designing a heat exchanger: the
manufacturing cost and the operating costs, including maintenance cost. In general,
the less the heat transfer surface area and the less the complexity of he design, the
lower is the manufacturing cost. The operating cost is the pumping cost due to
pumping devices such as fans, blowers, pumps, and etc. The maintenance costs
include costs of spares that require frequent renewal due to corrosion, and costs due to
corrosion/fouling prevention and control. Therefore, the heat exchanger design
requires a proper balance between thermal sizing and pressure drop.
Fabrication Techniques
Fabrication techniques are likely to be the determining factor o the selection of a heat
transfer surface matrix or core. They are the major factors in the initial cost and to a
large extent influence the integrity, service life, and ease of maintenance of the
finished heat exchanger. For example, shell and tube units are mostly fabricated by
welding, plate-fin heat exchanger and automobile aluminum radiators by brazing,
copper-brass radiators by soldering, most of the circular tube-fin exchangers by
mechanical assembling, etc.
E-22
Appendix C.5 Propane and propylene separation
Propane and propylene have similar boiling points (propane: -42.1 °C, propylene, -
47.70 °C) and as a result separation of these compounds requires highly complicated
units. Distillation is by far the most commonly used separation process in the
chemical industry today. The variants that are in use are:
1. Single-Column Process
2. Double-Column Process
3. Heat Pump Process
Single-Column Process: This process requires a large number of trays (150 – 200),
resulting in units of about 100 meters. The reflux can be condensed with cooling
water (column pressure 16 – 19 bar) or in air coolers (column pressure 21 – 26 bar).
Double-Column Process: For the large throughputs that are common today, the
double column process is preferred over the single-column process, since it does
requires smaller columns with smaller column diameters, which makes transportation
of these units easier. A schematic of the double-column process is given in Figure
C.5.1. Only the reflux from the second column is condensed with cooling water. The
pressure of the first column is sufficiently high (ca. 25 bar) that the overhead vapours
(ca. 59 °C) can be condensed in the reboiler of the second column and serve as the
heat carrier. Heating the first column with warm water is still possible. Both columns
provide approximately half the propene product. Since the reboiler for the second
column also serves as the condenser of the first column, the first column does not
require any cooling water. As a result the cooling water requirements are about half
that of the single column process.
propane/propene
I II
propane
propene
E-23
Figure C.5.1. Schematic of double column process [adapted from Ullman]
Heat pump process: In the aforementioned processes, the heat for the reboiler is
usually available as waste heat from e.g. the steam cracker, and is essentially cost free.
If this heat is not available, a heat pump can be used. A schematic of the heat pump
process is given in Figure C.5.2. The overhead vapors are heated slightly in the reflux
sub cooler, which enables these vapors to be compressed and cooled in the condenser-
reboiler.
reflux subcooler
compressor
propane/propene
I
after-cooler
condensor-reboiler
propene
propane
Figure C.5.2 Schematic of heat pump process [adapted from Ullman]
Thermodynamic analysis
E-24
Q out
rD
L
D
D /B = 1
V Q in
Per mole of feed F, the distillate D amounts to ½ mole and so does the bottom
product B. With a reflux ratio r = L/D, in which L is the number of moles of liquid
re-introduced at the top of the column, L = rD = ½r moles above the feedpoint and
(½r + 1) below the feedpoint. The vapour flow V = ½(r+1) throughout the column.
The heat introduced at the bottom of the column is therefore
Qin = 1 (1 + r )∆ v H (C.5.1)
2
We assume that the heat of vaporisation is roughly the same for both components and
that the temperature dependence is negligible over the range of the column. From this,
it follows that
Qin = Qout (C.5.2)
in which Qout is the cooling duty of the condenser. Now, it is interesting to note that
the overall separation does not require any energy! The number of Joules entering the
column equals the number of Joules leaving the column. However, the “quality” of
these heat streams, or equivalently, the exergy of these heat streams is not equal, due
to the Carnot factor. In an ideal column, that is a column operating under reversible
conditions, the heat is stripped of its quality and pays for the separation of the liquid
mixture into its constituents in the liquid state. The minimum work required to
separate the liquid mixture into its constituents is given below (see Figure. C.5.4):
ideal
Wsep = − RT0 ∑ xi lnxi (C.5.3)
i
where the assumption is made that the mixture behaves in an ideal fashion and is
close to the temperature of the surroundings, which, for the propane-propylene
mixture is a fair assumption. If we further insert the assumption that the mixture is
equimolar, equation (C.4.3) reduces to:
E-25
ideal
Wsep = RT0 ln 2 (C.5.4)
which is the minimum amount of exergy that needs to be introduced into the column.
And from
⎛
min ⎜ 1 1 ⎞⎟
Winmin = Qin T0 − (C.5.5)
⎜ Ttop Tbottom ⎟
⎝ ⎠
which has to equal the minimum amount of work that has to be spent on the
separation:
Winmin = Qsep
ideal
(C.5.6)
This also defines the minimum reflux ratio, r, according to
min 1 min
Qin = r
2
(
+1 ∆v H ) (C.5.7)
Qout
rD
L
D
F
Winmin
D/B = 1
V Qin
We stress that this equation dictates the minimum reflux ratio based purely on
E-26
The previous analysis begs the following question. What will the efficiency be of a
“real” propane-propylene distillation column? To answer this question, we must
realize that heat cannot be transferred into the column without a temperature
difference. In a “real” column with less stringent product quality constraints, the heat
is supplied at 377 K, the bottom temperature of the column is 331 K, the temperature
at the top of the column is 320 K and is transferred to the surroundings at 298 K
[Seader]. The minimum heat required for separation is, according to equation (C.4.7):
min 1 min
Qin = r
2
(
+1 ∆v H ) (C.5.9)
with rmin = 9.64 from the data in [Seader]. The separation inside the column does not
take place according to thermodynamic ideal processes, and the real heat is larger:
real
Qin
2
( )
= 1 r real + 1 ∆ v H (C.5.10)
The heat has to be transferred over a temperature difference of 377 – 331 = 46 K and
the resulting lost work can easily be calculated. Then the heat flows from 331 K to
320 K inside the column and is used to perform the separation. Finally, the heat is
discarded at the top of the column at 320 K to the surroundings at 298 K. The overall
thermodynamic efficiency of the column can be computed as follows:
min ⎛ 1 ⎞
Qin ∆⎜ ⎟ T0
⎝ T ⎠ column
ηoverall = (C.5.11)
real
⎡ ⎛ 1 ⎞ ⎛ 1 ⎞ ⎛ 1 ⎞ ⎤
Qin ⎢ ∆⎜ ⎟ + ∆⎜ ⎟ + ∆⎜ ⎟ ⎥T0
⎢⎣ ⎝ T ⎠bottom ⎝ T ⎠ column ⎝ T ⎠top ⎥⎦
which yields ηoverall = 0.093:
Closer scrutiny of equation (C.5.11) reveals that the main sources of inefficiency are
min
Qin
real
the temperature driving forces in the condenser and reboiler (the ratio Qin equals
0.63 and quantifies losses inside the column). The only way of improving the
efficiency is to reduce these temperature-driving forces. A noteworthy point is,
however, that the heat for the reboiler should be supplied at 377 K, which is often
available as waste heat in a chemical plant and integration with other heat sources
should therefore be contemplated. For the purpose of this example, however, we will
not do so. Another way of improving the single-column process is to use a membrane
to split the feed into two different feed streams. See Figure C.4.5.
E-27
Qout
rD
D
xD
xR
F
D/B = 1
xF xP
Qin
B
xB
E-28
Appendix C.6 Comparison of the Tray and Packed column
properties
Tray column
1. Generally employed in large diameter (larger than 1m)
2. Several down comers necessary
3. Gas load in tray columns must be kept within a relatively narrow range
4. Only valve trays allow greater operational flexibility
5. Liquid load can be varied over a very wide range
6. Allow heat to be added or removed easily
Disadvantages:
1. Relatively high pressure drop (7mbar/equilibrium stage)
2. Decomposition of thermally unstable substances
Packed column
1. Small diameter (smaller than 0.7m), development allows large diameter also
2. Extremely flexible as far as gas load is concerned, but require a minimum
liquid load
3. Small pressure loss (0.5mbar/equilibrium stage), more than 1 order magnitude
lower than in tray columns
4. Liquid hold up very small, thermally unstable substance is less
5. Countercurrent of gas and liquid, efficiency of mass transfer
6. Ceramic packing. Less corrosion
Disadvantages
1. Not suitable for the treatment of liquids that obtain particulate contaminants
or tent to crystallize
E-29
Appendix C.7 Recommendations for treatment of light gas
Membrane Technology
Membrane is a novel technology in separation field. One can select a qualified
material for the specified purpose of separation. Hydrogen and carbon monoxide are
the composition of synthesis gas and methane is natural gas. Therefore, to separate
methane from the stock is the desired way. Unfortunately, after searching the huge
amount of literature, we found that due to the fact that the molecular size of methane
and carbon monoxide are very close (Diameter of CO: 3.76, Diameter of CH4: 3.80),
it is not applicable by separate only by physical sense.
By physical properties membrane separation, only hydrogen can be separated from
the stock.
Membrane Reactor
The wide application of membrane reactor is methane steam reforming. In the
membrane reactor, two reactions will happen in parallel:
CH 4 steam reforming: CH 4 + H 2O → CO + 3H 2
Water shifting: CO + H 2O → CO2 + H 2
Methane will convert to CO and H2, which is valuable syngas can be sold to Fischer-
Tropsch process. However, this treatment is quite time consuming in calculation of
the reactor parameter design and very costly. And this process will bring the by-
product CO2 that requires further treatment for the sustainability point of view.
Absorption
In order to separate methane from the stock, one can use the chemical way such as
using a right solvent. It is required to cost money for the solvent and to design an
extra equipment.
Chromatography
Gas Chromatography (GC) is used for the qualitative and quantitative analysis of
complex mixtures of gases, liquids, and sometimes solids. A sample is vaporized and
transported by an inert carrier (usually He) gas through a column of sufficient length
to provide the separation. There are many different columns with both different
E-30
mobile and solid phases used with different carrier gases. The right combination has
to be chosen such that the desired separation takes place with reasonable retention
times. The different components of the vapor mixture are separated as a result of their
different vapor pressures and relative affinities for the bonded liquid phase. As the
components of the mixture are separated and elude from the column, they enter a
detector, where a signal proportional to the concentration of the component is
amplified and displayed. Identification is done by the retention time.
The Chromatography is more commonly used in analysis way and not applicable in
the huge amount separation unit. Also the problem is it is very costly.
Conclusion:
In this design, we choose to sale whole products mixture as our decision. The other
alternatives can serve as the candidates if there are specified requirements for the
treatment, such as purity and H2/CO ratio.
E-31
APPENDIX D
E-1
Appendix D. Process Stream Summary
E-2
Table D.1.a Process stream summary (Con’t)
STREAM Nr. : A203 A204 A205 A206 A207
Tube gas prod. After
Name : Tube Feed furnace Feed to tube reactor Tube gas Prod. cooler Process water
COMP MW kg/s kmol/s kg/s kmol/s kg/s kmol/s kg/s kmol/s kg/s kmol/s
Propane 44 6.54 0.1487 6.54 0.1487 1.36 0.0309 1.36 0.0309 0.00 0.0000
Oxygen 32 1.79 0.0559 1.79 0.0559 0.00 0.0000 0.00 0.0000 0.00 0.0000
Ethylene 28 0.00 0.0000 0.00 0.0000 2.29 0.0818 2.29 0.0818 0.00 0.0000
Propylene 42 0.68 0.0163 0.68 0.0163 2.47 0.0588 2.47 0.0588 0.00 0.0000
Carbonmon-oxide 28 0.00 0.0000 0.00 0.0000 0.62 0.0223 0.62 0.0223 0.00 0.0000
Carbondi-oxide 44 0.00 0.0000 0.00 0.0000 0.18 0.0040 0.18 0.0040 0.00 0.0000
Methane 16 0.00 0.0000 0.00 0.0000 0.58 0.0360 0.58 0.0360 0.00 0.0000
Hydrogen 2 0.00 0.0000 0.00 0.0000 0.05 0.0265 0.05 0.0265 0.00 0.0000
Water 18 0.00 0.0000 0.00 0.0000 1.47 0.0817 0.12 0.0064 1.36 0.0753
Nitrogen 28 0.08 0.0029 0.08 0.0029 0.08 0.0029 0.08 0.0029 0.00 0.0000
Total 9.10 0.2239 9.10 0.2239 9.10 0.3449 7.74 0.2695 1.36 0.0753
Enthalpy kW -4308 4957 -4243 -6063 -21804
Phase V V V V L
Press. Bara 1 1 1 1 1
Temp oC 539.978075 850 850 25 25
E-3
Table D.1.a Process stream summary (Con’t)
STREAM Nr. : 301 302 303 304 305
Name : Gas Prod. to LG column Overhead T301 Light gas Product T301 Reflux Feed to T302
COMP MW kg/s kmol/s kg/s kmol/s kg/s kmol/s kg/s kmol/s kg/s kmol/s
Propane 44 6.44 0.1464 0.00 0.0000 0.00 0.0000 0.00 0.0000 6.44 0.1464
Oxygen 32 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Ethylene 28 2.29 0.0819 7.90 0.2821 0.03 0.0010 7.87 0.2811 2.26 0.0809
Propylene 42 4.69 0.1116 0.00 0.0000 0.00 0.0000 0.00 0.0000 4.69 0.1116
Carbonmon-oxide 28 0.62 0.0223 1.19 0.0426 0.62 0.0223 0.57 0.0203 0.00 0.0000
Carbondi-oxide 44 0.01 0.0003 0.01 0.0002 0.00 0.0000 0.01 0.0002 0.01 0.0003
Methane 16 0.58 0.0361 4.79 0.2994 0.58 0.0361 4.21 0.2633 0.00 0.0000
Hydrogen 2 0.15 0.0769 0.16 0.0797 0.15 0.0769 0.01 0.0027 0.00 0.0000
Water 18 0.01 0.0006 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.01 0.0006
Nitrogen 28 0.08 0.0029 0.14 0.0050 0.08 0.0029 0.06 0.0021 0.00 0.0000
Total 14.88 0.4790 14.19 0.7089 1.47 0.1392 12.72 0.5697 13.42 0.3398
Enthalpy kW -14585 -5787 -5786 -5787 -14573
Phase V V V V L
Press. Bara 30 15 15 15 15
Temp oC 42 -131 -131 -131 10
E-4
Table D.1.a Process stream summary (Con’t)
STREAM
Nr. : 311 312 313 314 314-W
T303
Name : Overhead T303 Reflux Propylene Prduct Propane recycle1 Moisture removed
COMP MW kg/s kmol/s kg/s kmol/s kg/s kmol/s kg/s kmol/s kg/s kmol/s
Propane 44 0.05 0.0011 0.05 0.0011 0.00 0.0001 8.91 0.2025 0.00 0.0000
Oxygen 32 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Ethylene 28 0.04 0.0015 0.04 0.0014 0.00 0.0001 0.00 0.0000 0.00 0.0000
Propylene 42 105.15 2.5036 100.44 2.3914 4.71 0.1121 0.92 0.0219 0.00 0.0000
Carbonmon-oxide 28 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Carbondi-oxide 44 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Methane 16 0.00 0.0002 0.00 0.0002 0.00 0.0000 0.00 0.0000 0.00 0.0000
Hydrogen 2 0.05 0.0242 0.05 0.0231 0.00 0.0011 0.00 0.0000 0.00 0.0000
Water 18 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.01 0.0006
Nitrogen 28 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Total 105.29 2.5306 100.58 2.4172 4.72 0.1133 9.83 0.2244 0.01 0.0006
Enthalpy kW 0 0 2093 -24058 -168
Phase V L V L L
Press. Bara 15 15 15 15 15
Temp oC 35 35 35 43 43
E-5
APPENDIX E
E-1
Appendix E.1 Aspen Plus simulation results
E-2
303
D301
302
T301 304
106
E101
D101
306
103
C101
M301
105 T302
305 310
E102
301
104
M302
307
T303
V301 309
312
C201
M001
001 313 E301
RX-ENDO
D201
E203
101
206 208 209
201
RX-EXO
M201
204
E202
E201 210
202 205 212
203
207
314
E-3
Table E.1.1 The results of Aspen Plus simulation
Heat and Material Balance Table
Stream ID 001 101 102 103 104 105 106 107
Temperature C 540.0 540.0 540.0 25.0 215.1 55.0 54.7 54.7
Pressure bar 1.000 1.000 1.000 1.000 30.000 30.000 30.000 30.000
Vapor Frac 1.000 1.000 1.000 1.000 1.000 0.730 1.000 0.000
Mole Flow kmol/hr 675.348 888.943 1073.215 1073.215 1073.215 1073.215 785.355 287.860
Mass Flow kg/hr 29644.355 39022.384 39022.384 39022.384 39022.384 39022.384 26657.857 12364.527
Volume Flow cum/hr 45666.790 60110.020 72578.525 26334.839 1379.850 587.287 547.674 28.420
Enthalpy MMkcal/hr -5.121 -6.772 -1.069 -14.167 -10.581 -14.979 -9.506 -5.473
Mass Flow kg/hr
C3H8 26802.447 35329.320 27203.576 27203.576 27203.576 27203.576 18291.666 8911.911
O2
N2
C2H4 10.339 10.339 10.339 10.339 8.799 1.540
C3H6 2841.908 3693.064 11431.830 11431.830 11431.830 11431.830 7988.958 3442.872
CH4 5.912 5.912 5.912 5.912 5.497 0.416
H2 370.727 370.727 370.727 370.727 362.939 7.788
H2O
CO
CO2
MDEA+
Mole Flow kmol/hr
C3H8 607.813 801.182 616.910 616.910 616.910 616.910 414.810 202.100
O2
N2
C2H4 0.369 0.369 0.369 0.369 0.314 0.055
C3H6 67.535 87.762 271.665 271.665 271.665 271.665 189.849 81.816
CH4 0.369 0.369 0.369 0.369 0.343 0.026
H2 183.903 183.903 183.903 183.903 180.040 3.863
H2O
CO
CO2
MDEA+
Mole Frac
C3H8 0.900 0.901 0.575 0.575 0.575 0.575 0.528 0.702
O2
N2
C2H4 343 PPM 343 PPM 343 PPM 343 PPM 399 PPM 191 PPM
C3H6 0.100 0.099 0.253 0.253 0.253 0.253 0.242 0.284
CH4 343 PPM 343 PPM 343 PPM 343 PPM 436 PPM 90 PPM
H2 0.171 0.171 0.171 0.171 0.229 0.013
H2O
CO
CO2
MDEA+
E-4
Table E.1.1 The results of Aspen Plus simulation (Con’t)
Heat and Material Balance Table
Stream ID 201 202 203 204 205 206 207 208 209 210 211 212
Temperature C 540.0 540.0 540.0 850.0 850.0 25.0 25.0 288.6 30.0 30.0 30.0 30.0
Pressure bar 1.000 1.000 1.000 1.000 1.000 1.000 1.000 30.000 30.000 30.000 30.000 30.000
Vapor Frac 1.000 1.000 1.000 1.000 1.000 1.000 0.000 1.000 1.000 0.000 1.000 1.000
Mole Flow kmol/hr 592.629 212.000 804.629 804.629 1238.872 967.996 270.876 967.996 947.059 20.937 933.860 13.198
Mass Flow kg/hr 26014.923 6741.501 32756.424 32756.424 32755.239 27875.317 4879.922 27875.317 27498.122 377.194 26917.260 580.862
Volume Flow cum/hr 40073.347 14337.397 54415.038 75172.427 115724.474 23865.336 6.460 1501.437 667.676 0.378 658.313 9.362
Enthalpy MMkcal/hr -4.514 0.822 -3.692 4.249 -3.637 -5.197 -18.689 -1.647 -4.239 -1.427 -2.995 -1.244
Mass Flow kg/hr
C3H8 23552.880 23552.880 23552.880 4891.314 4891.314 < 0.001 4891.314 4891.314 4891.314
O2 6444.558 6444.558 6444.558
N2 296.943 296.943 296.943 296.943 296.943 trace 296.943 296.943 296.943
C2H4 8246.499 8246.497 0.002 8246.497 8246.497 8246.497
C3H6 2462.043 2462.043 2462.043 8885.550 8885.549 0.001 8885.549 8885.549 8885.549
CH4 2072.763 2072.763 < 0.001 2072.763 2072.763 2072.763
H2 190.757 190.757 trace 190.757 190.757 190.757
H2O 5294.891 414.984 4879.907 414.984 37.789 377.194 37.789
CO 2245.138 2245.138 < 0.001 2245.138 2245.138 2245.138
CO2 631.383 631.372 0.011 631.372 631.372 50.510 580.862
MDEA+
Mole Flow kmol/hr
C3H8 534.121 534.121 534.121 110.923 110.923 trace 110.923 110.923 110.923
O2 201.400 201.400 201.400
N2 10.600 10.600 10.600 10.600 10.600 trace 10.600 10.600 10.600
C2H4 293.953 293.953 < 0.001 293.953 293.953 293.953
C3H6 58.508 58.508 58.508 211.155 211.155 < 0.001 211.155 211.155 211.155
CH4 129.202 129.202 < 0.001 129.202 129.202 129.202
H2 94.627 94.627 trace 94.627 94.627 94.627
H2O 293.911 23.035 270.876 23.035 2.098 20.937 2.098
CO 80.154 80.154 trace 80.154 80.154 80.154
CO2 14.346 14.346 < 0.001 14.346 14.346 1.148 13.198
MDEA+
Mole Frac
C3H8 0.901 0.664 0.664 0.090 0.115 10 PPB 0.115 0.117 0.119
O2 0.950 0.250 0.250
N2 0.050 0.013 0.013 0.009 0.011 trace 0.011 0.011 0.011
C2H4 0.237 0.304 318 PPB 0.304 0.310 0.315
C3H6 0.099 0.073 0.073 0.170 0.218 77 PPB 0.218 0.223 0.226
CH4 0.104 0.133 66 PPB 0.133 0.136 0.138
H2 0.076 0.098 7 PPB 0.098 0.100 0.101
H2O 0.237 0.024 1.000 0.024 0.002 1.000 0.002
CO 0.065 0.083 4 PPB 0.083 0.085 0.086
CO2 0.012 0.015 929 PPB 0.015 0.015 0.001 1.000
MDEA+
E-5
Table E.1.1 The results of Aspen Plus simulation (Con’t)
Heat and Material Balance Table
Stream ID 301 302 303 304 305 306 307 308 309 310 311 312 313 314
Temperature C 42.5 -131.4 -131.4 -131.4 9.5 -39.0 39.9 37.4 37.9 34.8 42.7 42.7 42.7 540.0
Pressure bar 30.000 15.000 15.000 15.000 15.000 15.000 15.000 15.000 15.000 15.000 15.000 15.000 15.000 1.000
Vapor Frac 1.000 1.000 1.000 0.000 0.000 0.000 0.000 0.200 0.060 1.000 0.000 1.000 0.000 1.000
Mole Flow kmol/hr 1719.215 498.572 494.958 3.614 1220.642 292.954 927.688 287.860 1215.549 407.227 808.322 806.224 2.098 806.224
Mass Flow kg/hr 53575.117 5275.333 5173.947 101.386 48299.784 8255.705 40044.079 12364.527 52408.606 16977.864 35430.742 35392.953 37.789 35392.953
Volume Flow cum/hr 1225.285 371.054 369.247 0.173 105.967 19.365 93.977 102.049 214.237 542.593 84.936 1078.658 0.051 54516.577
Enthalpy MMkcal/hr -12.501 -4.960 -4.999 0.028 -12.491 2.566 -14.503 -5.473 -19.975 1.794 -20.621 -17.924 -0.144 -6.165
Mass Flow kg/hr
C3H8 23182.980 trace trace 23182.980 7.022 23175.958 8911.911 32087.869 8.116 32079.753 32079.753 32079.753
O2
N2 296.943 296.943 296.943 trace trace trace trace
C2H4 8255.295 101.386 101.386 8153.909 8148.609 5.300 1.540 6.841 6.841 trace
C3H6 16874.507 trace trace 16874.507 49.528 16824.979 3442.872 20267.851 16954.651 3313.199 3313.199 3313.199
CH4 2078.259 2078.124 2078.124 0.136 0.136 trace 0.416 0.416 0.416 trace
H2 553.696 553.696 553.696 trace trace trace 7.788 7.788 7.788 trace
H2O 37.789 trace 37.789 trace 37.789 37.789 trace 37.789 37.789
CO 2245.138 2245.138 2245.138 trace trace trace trace
CO2 50.510 0.047 0.047 50.462 50.410 0.052 0.052 0.052 trace
MDEA+
Mole Flow kmol/hr
C3H8 525.733 trace trace 525.733 0.159 525.573 202.100 727.673 0.184 727.489 727.489 727.489
O2
N2 10.600 10.600 10.600 trace trace trace trace
C2H4 294.267 3.614 3.614 290.653 290.464 0.189 0.055 0.244 0.244 trace
C3H6 401.004 trace trace 401.004 1.177 399.827 81.816 481.643 402.909 78.735 78.735 78.735
CH4 129.545 129.537 129.537 0.008 0.008 trace 0.026 0.026 0.026 trace
H2 274.667 274.667 274.667 trace trace trace 3.863 3.863 3.863 trace
H2O 2.098 trace 2.098 trace 2.098 2.098 trace 2.098 2.098
CO 80.154 80.154 80.154 trace trace trace trace
CO2 1.148 0.001 0.001 1.147 1.145 0.001 0.001 0.001 trace
MDEA+
Mole Frac
C3H8 0.306 trace trace 0.431 544 PPM 0.567 0.702 0.599 452 PPM 0.900 0.902 0.902
O2
N2 0.006 0.021 0.021 trace trace trace trace
C2H4 0.171 0.007 1.000 0.238 0.992 204 PPM 191 PPM 201 PPM 599 PPM trace
C3H6 0.233 trace 3 PPB 0.329 0.004 0.431 0.284 0.396 0.989 0.097 0.098 0.098
CH4 0.075 0.260 0.262 7 PPM 29 PPM trace 90 PPM 21 PPM 64 PPM trace
H2 0.160 0.551 0.555 trace trace trace 0.013 0.003 0.009 trace
H2O 0.001 trace 0.002 trace 0.002 0.002 trace 0.003 1.000
CO 0.047 0.161 0.162 trace trace trace trace
CO2 668 PPM 2 PPM 2 PPM 939 PPM 0.004 1 PPM 976 PPB 3 PPM trace
MDEA+
E-6
Appendix E.2 Simulation results of distillation columns T301
T302 and T303
**********************
**** INPUT DATA ****
**********************
NUMBER OF STAGES 20
ALGORITHM OPTION STANDARD
ABSORBER OPTION NO
INITIALIZATION OPTION STANDARD
HYDRAULIC PARAMETER CALCULATIONS NO
INSIDE LOOP CONVERGENCE METHOD BROYDEN
DESIGN SPECIFICATION METHOD NESTED
MAXIMUM NO. OF OUTSIDE LOOP ITERATIONS 200
MAXIMUM NO. OF INSIDE LOOP ITERATIONS 10
MAXIMUM NUMBER OF FLASH ITERATIONS 50
FLASH TOLERANCE 0.000100000
OUTSIDE LOOP CONVERGENCE TOLERANCE 0.000100000
*******************
**** RESULTS ****
*******************
OUTLET STREAMS
--------------
302 305
COMPONENT:
C3H8 .25189E-11 1.0000
E-7
N2 1.0000 .24305E-09
C2H4 .12281E-01 .98772
C3H6 .27710E-10 1.0000
CH4 .99993 .65298E-04
H2 1.0000 .13345E-13
H2O .64603E-20 1.0000
CO 1.0000 .72575E-09
CO2 .93762E-03 .99906
**NOTE** REPORTED VALUES FOR STAGE LIQUID AND VAPOR RATES ARE THE FLOWS
FROM THE STAGE EXCLUDING ANY SIDE PRODUCT. FOR THE FIRST STAGE,
THE REPORTED VAPOR FLOW IS THE VAPOR DISTILLATE FLOW. FOR THE
LAST STAGE, THE REPORTED LIQUID FLOW IS THE LIQUID BOTTOMS FLOW.
ENTHALPY
STAGE TEMPERATURE PRESSURE KCAL/MOL HEAT DUTY
C BAR LIQUID VAPOR MMKCAL/H
E-8
5 2374. 2870.
6 2351. 2872.
7 2212. 2850.
8 1826. 2711.
9 1474. 2325. 1719.2147
10 1560. 253.8
11 1578. 339.4
19 1587. 368.8
20 1221. 366.2 1220.6425
E-9
1 0.26562E-11 0.21261E-01 0.72487E-02 0.22288E-10 0.25981
2 0.13375E-07 0.71006E-02 0.39828 0.68259E-07 0.42238
4 0.32873E-05 0.39461E-02 0.79407 0.99057E-05 0.72777E-01
5 0.29753E-04 0.38918E-02 0.80976 0.70509E-04 0.58272E-01
6 0.26148E-03 0.38855E-02 0.81175 0.48798E-03 0.55537E-01
7 0.22681E-02 0.39130E-02 0.80582 0.33353E-02 0.55328E-01
8 0.18587E-01 0.40962E-02 0.76265 0.21595E-01 0.57192E-01
9 0.11009 0.47186E-02 0.56189 0.10257 0.63605E-01
10 0.13253 0.10914E-02 0.67435 0.12272 0.46920E-01
11 0.13771 0.21258E-03 0.70232 0.12737 0.27161E-01
19 0.16177 0.30406E-09 0.68600 0.14961 0.15557E-03
20 0.22127 0.49524E-10 0.57808 0.19830 0.59509E-04
E-10
2 0.63496E-06 0.48763E-03 0.92091 0.23369E-05 0.72358E-01
4 0.57036E-04 0.24099E-03 0.98602 0.12899E-03 0.91573E-02
5 0.50016E-03 0.23689E-03 0.98628 0.89073E-03 0.72632E-02
6 0.43300E-02 0.23435E-03 0.97714 0.60762E-02 0.68553E-02
7 0.35062E-01 0.22234E-03 0.91367 0.38874E-01 0.64525E-02
8 0.19310 0.17655E-03 0.62526 0.17169 0.50274E-02
9 0.43899 0.13808E-03 0.23056 0.32365 0.34019E-02
10 0.42736 0.34212E-04 0.25124 0.31647 0.25057E-02
11 0.42498 0.67522E-05 0.25637 0.31503 0.14489E-02
19 0.43933 0.95268E-11 0.23140 0.32725 0.79686E-05
20 0.47998 0.14942E-11 0.16882 0.34937 0.28097E-05
********************************
E-11
***** HYDRAULIC PARAMETERS *****
********************************
TEMPERATURE
C
STAGE LIQUID FROM VAPOR TO
1 -131.44 -67.400
2 -67.400 -51.158
4 -47.815 -47.138
5 -47.138 -46.810
6 -46.810 -45.298
7 -45.298 -36.418
8 -36.418 -11.228
9 -11.228 23.802
10 -4.4281 -2.9161
11 -2.9161 -2.3618
19 1.4536 9.5163
20 9.5163 9.5163
E-12
6 442.20 23.499 0.86141E-01 0.87204E-02 6.6716
7 447.94 22.556 0.88349E-01 0.89623E-02 6.8618
8 470.10 21.540 0.98510E-01 0.95095E-02 7.7388
9 481.76 21.257 0.10657 0.10080E-01 8.4143
10 466.85 25.580 0.98158E-01 0.94588E-02 7.3843
11 463.79 25.777 0.96570E-01 0.94570E-02 7.1710
19 460.49 26.824 0.95084E-01 0.95824E-02 6.9703
20 455.80 26.824 0.93374E-01 0.95824E-02 6.8167
***************************************
***** PACKING SIZING CALCULATIONS *****
***************************************
*******************
*** SECTION 1 ***
*******************
DESIGN PARAMETERS
-----------------
OVERDESIGN FACTOR 1.00000
SYSTEM FOAMING FACTOR 1.00000
FRAC. APP. TO MAXIMUM CAPACITY 1.00000
MAXIMUM CAPACITY FACTOR M/SEC MISSING
DESIGN CAPACITY FACTOR M/SEC MISSING
PRESSURE DROP FOR THE SECTION BAR MISSING
PRESSURE DROP PER UNIT HEIGHT MM-WATER/M MISSING
PACKING SPECIFICATIONS
----------------------
PACKING TYPE PALL-RING
PACKING MATERIAL METAL
PACKING SIZE 50-MM
VENDOR RASCHIG
PACKING SURFACE AREA SQCM/CC 1.05000
E-13
PACKING VOID FRACTION 0.96000
HETP METER 0.85000
PACKING HEIGHT METER 15.3000
HEIGHT
FROM TOP FRACTIONAL PRESSURE PRESSURE LIQUID
STAGE OF SECTION CAPACITY DROP DROP/HEIGHT HOLDUP HETP
METER BAR MM-WATER/M CUM METER
2 0.0000E+00 0.8974 0.55373E-03 6.6429 0.1292 0.8500
3 0.8500 0.9770 0.66984E-03 8.0358 0.1393 0.8500
4 1.700 0.9967 0.70300E-03 8.4337 0.1419 0.8500
5 2.550 1.000 0.70908E-03 8.5066 0.1423 0.8500
6 3.400 0.9959 0.70378E-03 8.4430 0.1416 0.8500
7 4.250 0.9666 0.66496E-03 7.9773 0.1369 0.8500
8 5.100 0.8836 0.56537E-03 6.7826 0.1265 0.8500
9 5.950 0.8380 0.52178E-03 6.2597 0.1223 0.8500
10 6.800 0.3483 0.15217E-04 0.18255 0.1250 0.8500
11 7.650 0.3592 0.17070E-04 0.20479 0.1263 0.8500
12 8.500 0.3630 0.17755E-04 0.21300 0.1268 0.8500
13 9.350 0.3648 0.18088E-04 0.21700 0.1270 0.8500
14 10.20 0.3658 0.18266E-04 0.21913 0.1271 0.8500
15 11.05 0.3663 0.18366E-04 0.22032 0.1272 0.8500
16 11.90 0.3666 0.18436E-04 0.22117 0.1272 0.8500
17 12.75 0.3671 0.18538E-04 0.22239 0.1273 0.8500
18 13.60 0.3686 0.18846E-04 0.22609 0.1275 0.8500
19 14.45 0.3743 0.20015E-04 0.24011 0.1284 0.8500
**********************
**** INPUT DATA ****
E-14
**********************
NUMBER OF STAGES 24
ALGORITHM OPTION STANDARD
ABSORBER OPTION NO
INITIALIZATION OPTION STANDARD
HYDRAULIC PARAMETER CALCULATIONS NO
INSIDE LOOP CONVERGENCE METHOD BROYDEN
DESIGN SPECIFICATION METHOD NESTED
MAXIMUM NO. OF OUTSIDE LOOP ITERATIONS 25
MAXIMUM NO. OF INSIDE LOOP ITERATIONS 10
MAXIMUM NUMBER OF FLASH ITERATIONS 50
FLASH TOLERANCE 0.000100000
OUTSIDE LOOP CONVERGENCE TOLERANCE 0.000100000
*******************
**** RESULTS ****
*******************
OUTLET STREAMS
--------------
306 307
COMPONENT:
C3H8 .30289E-03 .99970
N2 1.0000 .50283E-12
C2H4 .99935 .64987E-03
C3H6 .29351E-02 .99706
CH4 1.0000 .61080E-08
H2 1.0000 .94838E-15
H2O .14964E-11 1.0000
CO 1.0000 .10715E-11
CO2 .99897 .10339E-02
E-15
**** MANIPULATED VARIABLES ****
BOUNDS CALCULATED
LOWER UPPER VALUE
MOLAR REFLUX RATIO 0.10000 15.000 3.7639
**NOTE** REPORTED VALUES FOR STAGE LIQUID AND VAPOR RATES ARE THE FLOWS
FROM THE STAGE EXCLUDING ANY SIDE PRODUCT. FOR THE FIRST STAGE,
THE REPORTED VAPOR FLOW IS THE VAPOR DISTILLATE FLOW. FOR THE
LAST STAGE, THE REPORTED LIQUID FLOW IS THE LIQUID BOTTOMS FLOW.
ENTHALPY
STAGE TEMPERATURE PRESSURE KCAL/MOL HEAT DUTY
C BAR LIQUID VAPOR MMKCAL/H
E-16
LIQUID VAPOR LIQUID VAPOR MIXED LIQUID VAPOR
1 1103. 0.0000E+00 292.9542
2 1070. 1396.
3 970.3 1363.
5 762.3 1121.
6 761.5 1055.
7 768.5 1054.
12 772.3 1066.
13 771.5 1065. .40890-01
14 2023. 1064. 1220.6016
15 2066. 1095.
22 2162. 1233.
23 2162. 1234.
24 927.7 1235. 927.6883
E-17
15 0.00000E+00 0.11846E-02 0.23231E-12 0.36177E-03
22 0.00000E+00 0.12154E-02 0.00000E+00 0.58507E-05
23 0.00000E+00 0.13754E-02 0.00000E+00 0.28948E-05
24 0.00000E+00 0.22611E-02 0.00000E+00 0.12779E-05
E-18
12 44.762 0.18154 21.431 2.6251
13 44.568 0.18355 21.387 2.6344
14 44.460 0.18529 21.378 2.6447
15 41.770 0.22677 21.022 2.8638
22 37.036 0.31073 20.042 3.2033
23 36.952 0.31228 20.019 3.2082
24 36.905 0.31395 20.018 3.2158
E-19
5 0.20857E-15 0.98237E-10 0.47136E-10 0.51554E-02
6 0.20401E-15 0.86941E-09 0.45898E-10 0.33127E-02
7 0.19587E-15 0.48720E-08 0.44019E-10 0.23395E-02
12 0.18812E-15 0.66813E-05 0.42240E-10 0.17494E-02
13 0.18797E-15 0.26646E-04 0.42207E-10 0.17442E-02
14 0.18760E-15 0.10510E-03 0.42162E-10 0.17403E-02
15 0.00000E+00 0.12263E-03 0.34662E-11 0.11553E-02
22 0.00000E+00 0.15807E-03 0.00000E+00 0.19163E-04
23 0.00000E+00 0.17958E-03 0.00000E+00 0.94859E-05
24 0.00000E+00 0.29650E-03 0.00000E+00 0.41932E-05
********************************
***** HYDRAULIC PARAMETERS *****
********************************
TEMPERATURE
C
STAGE LIQUID FROM VAPOR TO
1 -39.006 -38.142
2 -38.142 -34.431
3 -34.431 -22.714
5 -3.8103 10.247
6 10.247 16.704
7 16.704 19.256
12 21.462 21.739
13 21.739 22.033
14 22.034 28.582
15 28.582 33.141
22 39.487 39.654
23 39.654 39.881
24 39.881 39.881
E-20
7 30964. 39220. 68.141 1389.1 40.291 36.948
12 31700. 39956. 70.855 1401.6 41.047 37.509
13 31719. 39974. 71.024 1401.8 41.113 37.554
14 83268. 43224. 186.73 1454.6 41.162 39.464
15 86596. 46552. 197.06 1516.7 41.920 40.905
22 93233. 53189. 218.43 1647.8 43.121 43.087
23 93303. 53259. 218.77 1648.9 43.146 43.132
24 40044. 53259. 93.977 1648.9 43.165 43.132
***************************************
***** PACKING SIZING CALCULATIONS *****
***************************************
*******************
*** SECTION 1 ***
*******************
E-21
PRESSURE PROFILE UPDATED NO
DESIGN PARAMETERS
-----------------
OVERDESIGN FACTOR 1.00000
SYSTEM FOAMING FACTOR 1.00000
FRAC. APP. TO MAXIMUM CAPACITY 1.00000
MAXIMUM CAPACITY FACTOR M/SEC MISSING
DESIGN CAPACITY FACTOR M/SEC MISSING
PRESSURE DROP FOR THE SECTION BAR MISSING
PRESSURE DROP PER UNIT HEIGHT MM-WATER/M MISSING
PACKING SPECIFICATIONS
----------------------
PACKING TYPE PALL-RING
PACKING MATERIAL METAL
PACKING SIZE 50-MM
VENDOR RASCHIG
PACKING SURFACE AREA SQCM/CC 1.05000
PACKING VOID FRACTION 0.96000
HETP METER 0.85000
PACKING HEIGHT METER 18.7000
HEIGHT
FROM TOP FRACTIONAL PRESSURE PRESSURE LIQUID
STAGE OF SECTION CAPACITY DROP DROP/HEIGHT HOLDUP HETP
METER BAR MM-WATER/M CUM METER
2 0.0000E+00 0.6166 0.21921E-03 2.6298 0.7964E-01 0.8500
3 0.8500 0.5882 0.20521E-03 2.4618 0.7595E-01 0.8500
4 1.700 0.5503 0.18697E-03 2.2430 0.7207E-01 0.8500
5 2.550 0.5506 0.18890E-03 2.2662 0.7317E-01 0.8500
6 3.400 0.5731 0.20308E-03 2.4362 0.7639E-01 0.8500
7 4.250 0.5881 0.21265E-03 2.5511 0.7840E-01 0.8500
8 5.100 0.5950 0.21703E-03 2.6037 0.7932E-01 0.8500
9 5.950 0.5980 0.21890E-03 2.6260 0.7974E-01 0.8500
10 6.800 0.5997 0.21980E-03 2.6368 0.7999E-01 0.8500
11 7.650 0.6008 0.22036E-03 2.6436 0.8017E-01 0.8500
12 8.500 0.6017 0.22078E-03 2.6487 0.8032E-01 0.8500
13 9.350 0.6024 0.22107E-03 2.6521 0.8045E-01 0.8500
14 10.20 0.8500 0.39393E-03 4.7258 0.1499 0.8500
15 11.05 0.8986 0.46044E-03 5.5238 0.1552 0.8500
16 11.90 0.9373 0.51993E-03 6.2374 0.1595 0.8500
17 12.75 0.9634 0.56414E-03 6.7678 0.1626 0.8500
18 13.60 0.9793 0.59291E-03 7.1129 0.1645 0.8500
19 14.45 0.9884 0.61012E-03 7.3194 0.1657 0.8500
20 15.30 0.9936 0.62001E-03 7.4381 0.1663 0.8500
21 16.15 0.9966 0.62577E-03 7.5071 0.1667 0.8500
22 17.00 0.9986 0.62939E-03 7.5505 0.1670 0.8500
E-22
23 17.85 1.000 0.63178E-03 7.5792 0.1672 0.8500
**********************
**** INPUT DATA ****
**********************
*******************
**** RESULTS ****
*******************
OUTLET STREAMS
--------------
310 311
COMPONENT:
C3H8 .25292E-03 .99975
E-23
C2H4 1.0000 .76830E-32
C3H6 .83653 .16347
CH4 1.0000 .21693E-60
H2 1.0000 .10998E-98
H2O .55068E-51 1.0000
CO2 1.0000 .64856E-31
BOUNDS CALCULATED
LOWER UPPER VALUE
DISTILLATE TO FEED RATIO 0.10000 1.0000 0.33501
**NOTE** REPORTED VALUES FOR STAGE LIQUID AND VAPOR RATES ARE THE FLOWS
FROM THE STAGE EXCLUDING ANY SIDE PRODUCT. FOR THE FIRST STAGE,
THE REPORTED VAPOR FLOW IS THE VAPOR DISTILLATE FLOW. FOR THE
LAST STAGE, THE REPORTED LIQUID FLOW IS THE LIQUID BOTTOMS FLOW.
E-24
ENTHALPY
STAGE TEMPERATURE PRESSURE KCAL/MOL HEAT DUTY
C BAR LIQUID VAPOR MMKCAL/H
E-25
6 0.95917E-03 0.12138E-04 0.99902 0.33499E-06 0.10823E-04
7 0.10724E-02 0.11868E-04 0.99890 0.33493E-06 0.10823E-04
8 0.11961E-02 0.11789E-04 0.99878 0.33493E-06 0.10824E-04
91 0.52205 0.11406E-04 0.47789 0.34330E-06 0.11847E-04
92 0.53178 0.11398E-04 0.46809 0.34341E-06 0.11865E-04
93 0.54044 0.10050E-04 0.45925 0.22906E-06 0.36656E-05
94 0.54098 0.32977E-05 0.45872 0.26867E-07 0.10835E-06
152 0.85689 0.00000E+00 0.14277 0.00000E+00 0.00000E+00
153 0.86835 0.00000E+00 0.13127 0.00000E+00 0.00000E+00
154 0.87959 0.00000E+00 0.11987 0.00000E+00 0.00000E+00
155 0.89037 0.00000E+00 0.10859 0.00000E+00 0.00000E+00
156 0.90000 0.00000E+00 0.97405E-01 0.00000E+00 0.00000E+00
E-26
155 0.35940E-03 0.00000E+00
156 0.89904E-03 0.00000E+00
E-27
91 0.13873E-04 0.59558E-07
92 0.42508E-04 0.59472E-07
93 0.12514E-03 0.52750E-07
94 0.12514E-03 0.17934E-07
152 0.13783E-03 0.00000E+00
153 0.15846E-03 0.00000E+00
154 0.22323E-03 0.00000E+00
155 0.42668E-03 0.00000E+00
156 0.10666E-02 0.00000E+00
********************************
***** HYDRAULIC PARAMETERS *****
********************************
E-28
ML IS THE MASS FLOW OF LIQUID FROM THE STAGE
MV IS THE MASS FLOW OF VAPOR TO THE STAGE
RHOL IS THE MASS DENSITY OF LIQUID FROM THE STAGE
RHOV IS THE MASS DENSITY OF VAPOR TO THE STAGE
QV IS THE VOLUMETRIC FLOW RATE OF VAPOR TO THE STAGE
TEMPERATURE
C
STAGE LIQUID FROM VAPOR TO
1 34.806 35.336
2 35.336 35.354
6 35.358 35.359
7 35.359 35.360
8 35.360 35.361
91 39.418 39.499
92 39.499 39.583
93 39.597 39.610
94 39.610 39.615
152 42.228 42.328
153 42.328 42.430
154 42.430 42.546
155 42.546 42.703
156 42.703 42.703
E-29
156 417.14 32.758 0.80855E-01 0.96326E-02 5.0282
***************************************
***** PACKING SIZING CALCULATIONS *****
***************************************
*******************
*** SECTION 1 ***
*******************
DESIGN PARAMETERS
-----------------
OVERDESIGN FACTOR 1.00000
SYSTEM FOAMING FACTOR 1.00000
FRAC. APP. TO MAXIMUM CAPACITY 1.00000
MAXIMUM CAPACITY FACTOR M/SEC MISSING
DESIGN CAPACITY FACTOR M/SEC MISSING
PRESSURE DROP FOR THE SECTION BAR MISSING
PRESSURE DROP PER UNIT HEIGHT MM-WATER/M MISSING
PACKING SPECIFICATIONS
----------------------
PACKING TYPE PALL-RING
PACKING MATERIAL METAL
PACKING SIZE 80-MM
VENDOR RASCHIG
PACKING SURFACE AREA SQCM/CC 0.78000
PACKING VOID FRACTION 0.96000
HETP METER 0.85000
PACKING HEIGHT METER 130.900
E-30
***** SIZING RESULTS *****
HEIGHT
FROM TOP FRACTIONAL PRESSURE PRESSURE LIQUID
STAGE OF SECTION CAPACITY DROP DROP/HEIGHT HOLDUP HETP
METER BAR MM-WATER/M CUM METER
2 0.0000E+00 0.9252 0.68429E-03 8.2092 0.6172 0.8500
3 0.8500 0.9254 0.68460E-03 8.2129 0.6173 0.8500
4 1.700 0.9254 0.68463E-03 8.2132 0.6173 0.8500
5 2.550 0.9254 0.68464E-03 8.2134 0.6173 0.8500
6 3.400 0.9254 0.68465E-03 8.2135 0.6173 0.8500
7 4.250 0.9254 0.68466E-03 8.2137 0.6173 0.8500
8 5.100 0.9254 0.68467E-03 8.2138 0.6173 0.8500
9 5.950 0.9254 0.68469E-03 8.2140 0.6173 0.8500
10 6.800 0.9254 0.68470E-03 8.2141 0.6174 0.8500
11 7.650 0.9254 0.68472E-03 8.2143 0.6174 0.8500
12 8.500 0.9254 0.68473E-03 8.2145 0.6174 0.8500
13 9.350 0.9254 0.68475E-03 8.2147 0.6174 0.8500
14 10.20 0.9255 0.68477E-03 8.2149 0.6174 0.8500
15 11.05 0.9255 0.68479E-03 8.2152 0.6174 0.8500
16 11.90 0.9255 0.68481E-03 8.2155 0.6174 0.8500
17 12.75 0.9255 0.68484E-03 8.2158 0.6174 0.8500
18 13.60 0.9255 0.68487E-03 8.2161 0.6174 0.8500
19 14.45 0.9255 0.68490E-03 8.2165 0.6175 0.8500
20 15.30 0.9255 0.68493E-03 8.2168 0.6175 0.8500
21 16.15 0.9256 0.68496E-03 8.2173 0.6175 0.8500
22 17.00 0.9256 0.68500E-03 8.2178 0.6175 0.8500
23 17.85 0.9256 0.68505E-03 8.2183 0.6175 0.8500
24 18.70 0.9256 0.68509E-03 8.2188 0.6175 0.8500
25 19.55 0.9257 0.68514E-03 8.2194 0.6176 0.8500
26 20.40 0.9257 0.68520E-03 8.2201 0.6176 0.8500
27 21.25 0.9257 0.68526E-03 8.2208 0.6176 0.8500
28 22.10 0.9258 0.68533E-03 8.2216 0.6177 0.8500
29 22.95 0.9258 0.68540E-03 8.2225 0.6177 0.8500
30 23.80 0.9259 0.68548E-03 8.2235 0.6177 0.8500
31 24.65 0.9259 0.68557E-03 8.2245 0.6178 0.8500
32 25.50 0.9260 0.68566E-03 8.2256 0.6178 0.8500
33 26.35 0.9260 0.68576E-03 8.2269 0.6179 0.8500
34 27.20 0.9261 0.68588E-03 8.2282 0.6179 0.8500
35 28.05 0.9262 0.68600E-03 8.2297 0.6180 0.8500
36 28.90 0.9263 0.68613E-03 8.2313 0.6181 0.8500
37 29.75 0.9264 0.68628E-03 8.2331 0.6181 0.8500
38 30.60 0.9264 0.68644E-03 8.2350 0.6182 0.8500
39 31.45 0.9265 0.68661E-03 8.2370 0.6183 0.8500
40 32.30 0.9267 0.68680E-03 8.2393 0.6184 0.8500
41 33.15 0.9268 0.68700E-03 8.2417 0.6185 0.8500
42 34.00 0.9269 0.68722E-03 8.2444 0.6186 0.8500
43 34.85 0.9271 0.68746E-03 8.2473 0.6187 0.8500
44 35.70 0.9272 0.68773E-03 8.2504 0.6188 0.8500
45 36.55 0.9274 0.68801E-03 8.2538 0.6190 0.8500
46 37.40 0.9276 0.68832E-03 8.2575 0.6191 0.8500
E-31
47 38.25 0.9278 0.68865E-03 8.2615 0.6193 0.8500
48 39.10 0.9280 0.68901E-03 8.2658 0.6194 0.8500
49 39.95 0.9282 0.68940E-03 8.2705 0.6196 0.8500
50 40.80 0.9284 0.68982E-03 8.2756 0.6198 0.8500
51 41.65 0.9287 0.69028E-03 8.2810 0.6200 0.8500
52 42.50 0.9290 0.69077E-03 8.2869 0.6203 0.8500
53 43.35 0.9293 0.69129E-03 8.2932 0.6205 0.8500
54 44.20 0.9296 0.69186E-03 8.3000 0.6208 0.8500
55 45.05 0.9300 0.69247E-03 8.3073 0.6211 0.8500
56 45.90 0.9304 0.69312E-03 8.3151 0.6214 0.8500
57 46.75 0.9308 0.69382E-03 8.3235 0.6217 0.8500
58 47.60 0.9312 0.69456E-03 8.3324 0.6221 0.8500
59 48.45 0.9317 0.69535E-03 8.3419 0.6224 0.8500
60 49.30 0.9322 0.69619E-03 8.3520 0.6228 0.8500
61 50.15 0.9327 0.69709E-03 8.3627 0.6233 0.8500
62 51.00 0.9332 0.69804E-03 8.3741 0.6237 0.8500
63 51.85 0.9338 0.69903E-03 8.3861 0.6242 0.8500
64 52.70 0.9344 0.70009E-03 8.3987 0.6247 0.8500
65 53.55 0.9351 0.70119E-03 8.4120 0.6252 0.8500
66 54.40 0.9357 0.70235E-03 8.4259 0.6257 0.8500
67 55.25 0.9364 0.70356E-03 8.4404 0.6263 0.8500
68 56.10 0.9372 0.70482E-03 8.4555 0.6269 0.8500
69 56.95 0.9379 0.70613E-03 8.4712 0.6275 0.8500
70 57.80 0.9387 0.70749E-03 8.4875 0.6281 0.8500
71 58.65 0.9395 0.70888E-03 8.5042 0.6287 0.8500
72 59.50 0.9403 0.71031E-03 8.5214 0.6294 0.8500
73 60.35 0.9411 0.71178E-03 8.5390 0.6300 0.8500
74 61.20 0.9419 0.71328E-03 8.5569 0.6307 0.8500
75 62.05 0.9428 0.71480E-03 8.5752 0.6314 0.8500
76 62.90 0.9437 0.71633E-03 8.5936 0.6321 0.8500
77 63.75 0.9445 0.71788E-03 8.6122 0.6328 0.8500
78 64.60 0.9454 0.71944E-03 8.6309 0.6335 0.8500
79 65.45 0.9462 0.72100E-03 8.6496 0.6341 0.8500
80 66.30 0.9471 0.72255E-03 8.6682 0.6348 0.8500
81 67.15 0.9479 0.72409E-03 8.6867 0.6355 0.8500
82 68.00 0.9488 0.72562E-03 8.7050 0.6362 0.8500
83 68.85 0.9496 0.72712E-03 8.7230 0.6368 0.8500
84 69.70 0.9504 0.72859E-03 8.7407 0.6375 0.8500
85 70.55 0.9512 0.73003E-03 8.7579 0.6381 0.8500
86 71.40 0.9520 0.73144E-03 8.7748 0.6387 0.8500
87 72.25 0.9527 0.73280E-03 8.7911 0.6393 0.8500
88 73.10 0.9534 0.73412E-03 8.8069 0.6398 0.8500
89 73.95 0.9541 0.73538E-03 8.8221 0.6404 0.8500
90 74.80 0.9547 0.73658E-03 8.8365 0.6409 0.8500
91 75.65 0.9553 0.73768E-03 8.8497 0.6414 0.8500
92 76.50 0.9558 0.73857E-03 8.8604 0.6418 0.8500
93 77.35 0.9767 0.81731E-03 9.8050 0.6968 0.8500
94 78.20 0.9768 0.81764E-03 9.8089 0.6969 0.8500
95 79.05 0.9768 0.81773E-03 9.8101 0.6969 0.8500
96 79.90 0.9769 0.81783E-03 9.8112 0.6970 0.8500
97 80.75 0.9769 0.81793E-03 9.8124 0.6970 0.8500
98 81.60 0.9770 0.81803E-03 9.8137 0.6971 0.8500
99 82.45 0.9770 0.81815E-03 9.8151 0.6971 0.8500
100 83.30 0.9771 0.81827E-03 9.8165 0.6971 0.8500
101 84.15 0.9772 0.81841E-03 9.8182 0.6972 0.8500
102 85.00 0.9772 0.81855E-03 9.8199 0.6972 0.8500
103 85.85 0.9773 0.81871E-03 9.8218 0.6973 0.8500
104 86.70 0.9774 0.81888E-03 9.8238 0.6974 0.8500
105 87.55 0.9774 0.81906E-03 9.8260 0.6974 0.8500
106 88.40 0.9775 0.81926E-03 9.8284 0.6975 0.8500
107 89.25 0.9776 0.81947E-03 9.8309 0.6976 0.8500
E-32
108 90.10 0.9777 0.81970E-03 9.8337 0.6977 0.8500
109 90.95 0.9778 0.81995E-03 9.8366 0.6978 0.8500
110 91.80 0.9780 0.82021E-03 9.8398 0.6978 0.8500
111 92.65 0.9781 0.82050E-03 9.8432 0.6980 0.8500
112 93.50 0.9782 0.82081E-03 9.8469 0.6981 0.8500
113 94.35 0.9784 0.82114E-03 9.8509 0.6982 0.8500
114 95.20 0.9785 0.82149E-03 9.8552 0.6983 0.8500
115 96.05 0.9787 0.82187E-03 9.8597 0.6985 0.8500
116 96.90 0.9789 0.82228E-03 9.8647 0.6986 0.8500
117 97.75 0.9791 0.82273E-03 9.8700 0.6988 0.8500
118 98.60 0.9793 0.82320E-03 9.8756 0.6989 0.8500
119 99.45 0.9795 0.82370E-03 9.8817 0.6991 0.8500
120 100.3 0.9798 0.82425E-03 9.8882 0.6993 0.8500
121 101.1 0.9800 0.82483E-03 9.8952 0.6995 0.8500
122 102.0 0.9803 0.82545E-03 9.9026 0.6997 0.8500
123 102.8 0.9806 0.82611E-03 9.9105 0.7000 0.8500
124 103.7 0.9809 0.82681E-03 9.9190 0.7002 0.8500
125 104.5 0.9813 0.82756E-03 9.9280 0.7005 0.8500
126 105.4 0.9816 0.82836E-03 9.9376 0.7008 0.8500
127 106.2 0.9820 0.82921E-03 9.9478 0.7011 0.8500
128 107.1 0.9824 0.83012E-03 9.9586 0.7014 0.8500
129 107.9 0.9828 0.83107E-03 9.9701 0.7018 0.8500
130 108.8 0.9833 0.83208E-03 9.9822 0.7021 0.8500
131 109.6 0.9838 0.83315E-03 9.9950 0.7025 0.8500
132 110.5 0.9843 0.83428E-03 10.009 0.7029 0.8500
133 111.3 0.9848 0.83546E-03 10.023 0.7033 0.8500
134 112.2 0.9853 0.83671E-03 10.038 0.7038 0.8500
135 113.0 0.9859 0.83802E-03 10.053 0.7042 0.8500
136 113.9 0.9865 0.83939E-03 10.070 0.7047 0.8500
137 114.7 0.9871 0.84082E-03 10.087 0.7052 0.8500
138 115.6 0.9878 0.84230E-03 10.105 0.7057 0.8500
139 116.4 0.9885 0.84385E-03 10.123 0.7063 0.8500
140 117.3 0.9892 0.84546E-03 10.143 0.7068 0.8500
141 118.1 0.9899 0.84712E-03 10.163 0.7074 0.8500
142 119.0 0.9906 0.84883E-03 10.183 0.7080 0.8500
143 119.8 0.9914 0.85059E-03 10.204 0.7086 0.8500
144 120.7 0.9922 0.85240E-03 10.226 0.7092 0.8500
145 121.5 0.9929 0.85426E-03 10.248 0.7099 0.8500
146 122.4 0.9938 0.85614E-03 10.271 0.7105 0.8500
147 123.2 0.9946 0.85807E-03 10.294 0.7112 0.8500
148 124.1 0.9954 0.86002E-03 10.317 0.7118 0.8500
149 124.9 0.9962 0.86199E-03 10.341 0.7125 0.8500
150 125.8 0.9971 0.86397E-03 10.365 0.7132 0.8500
151 126.6 0.9979 0.86594E-03 10.388 0.7138 0.8500
152 127.5 0.9987 0.86788E-03 10.412 0.7145 0.8500
153 128.3 0.9995 0.86965E-03 10.433 0.7151 0.8500
154 129.2 1.000 0.87093E-03 10.448 0.7156 0.8500
155 130.0 0.9999 0.87068E-03 10.445 0.7158 0.8500
E-33
Appendix E.3 Reactor Design
For this project, we propose to use Shell and Tube reactor. The endothermic reaction,
dehydrogenation of propane, takes place in the shell side and exothermic reaction, oxidative
dehydrogenation of propane, takes place in the tube side. In this chapter, it shows the idea
how to estimate the volume of reactor, volume and type of catalyst used and the material of
reactor.
There are many reactions take place in tube side. The rate expressions of reactions are
showed
3C 3 H 8 ( g ) + 9 O 2 → C 3 H 6 ( g ) + 3C O + 3C O 2 + 9 H 2 O ( g ) (E.3.1.1)
⎡ ⎛ k O 2 PO 2 ⎞
2 ⎤
k O 2 PO 2 k O2 PO2
( − r1 ) = k O2 PO2 ⎢1+ − ⎜ ⎟ + ⎥ (E.3.1.2)
⎢ 2 ⋅ k C H PC H ⎜ 2 ⋅ k C H PC H ⎟ k C H PC H ⎥
⎢⎣ 3 8 3 8 ⎝ 3 8 3 8 ⎠ 3 8 3 8
⎥⎦
2. Propylene oxidation to CO
2 C 3 H 6 ( g ) + 7 .5 O 2 → 3 C O ( g ) + 3 C O 2 ( g ) + 6 H 2 O ( g ) (E .3 .1 .3 )
( − r2 ) = k 2 P C m n
PO (E .3 .1 .4 )
3H 6 2
4. Cracking of propane
C3 H 8( g ) → CH 4( g ) + C2 H 4( g ) (E.3.1.7)
PC3H6 PH2
k 4 (PC3H8 − )
K1
(− r4 ) = (E.3.1.8)
1 + PC3H6 K 2
There are many reactions take place in shell side. The rate expressions of reactions are
showed
1. Cracking of propane
E-34
C3 H 8( g ) → C2 H 4( g ) + CH 4( g ) (E.3.2.1)
(A) (B) (C)
PCS3H6 PHS2
k S4 (PCS3H8 − )
K1
(− r5 ) = (E.3.2.2)
1 + PC3H6 K 2
S
2. Propane dehydrogenation
C3 H 8( g ) → C3 H 6( g ) + H 2( g ) (E.3.2.3)
(A) (D) (E)
(− r6 ) = k 6 PCO,SH (E.3.2.4)
3 8
We consider that the Tube side of this reactor is plug flow reactor system both sides, and then
make the mass balance and energy balance at the steady state condition.
For Propane
In = O u t + C o n v e rs io n + h o ld u p
dC C
FT C C 3 H 8 = FT C C 3 H 8 + ( − rC ) d VT + VT 3H 8
V =V o V =V o+dV 3H 8
dt
d C C3H 8 ⎛ ( − r1 ) + ( − r4 ) ⎞
− = ⎜ ⎟ (E .3 .3 .1 )
d VT ⎝ FT ⎠
For oxygen
dCO2 ⎛ 3( −r1 ) + 3.75( −r2 ) + 0.5( −r3 ) ⎞
− = ⎜ ⎟ (E.3.3.2)
dVT ⎝ FT ⎠
For propylene
dC C H ⎛ − 1 / 3 ( − r1 ) + ( − r2 ) ⎞
− 3 6
= ⎜ ⎟ ( E .3 .3 .3 )
d VT ⎝ FT ⎠
E-35
dC CO2 ⎛ − ( − r1 ) − 1.5 ( − r2 ) − ( − r3 ) ⎞
− =⎜ ⎟ (E.3.3.5)
dVT ⎝ FT ⎠
For water
d C H 2O ⎛ − 3 ( − r1 ) − 3 ( − r2 ) ⎞
− =⎜ ⎟ (E .3 .3 .6 )
d VT ⎝ FT ⎠
For methane
dC CH 4 ⎛ − ( − r4 )⎞
− = ⎜ ⎟ ( E .3 .3 .7 )
d VT ⎝ FT ⎠
For ethylene
d C C2H 4 ⎛ − ( − r4 )⎞
− = ⎜ ⎟ ( E .3 .3 .8 )
d VT ⎝ FT ⎠
Energy balance:
In = O ut + h e a t tra n s fe r + h e a t g e n e ra te d + h e a t a c c u m u la te
ρ T , m ix FT C P , m ix T T i = ρ T , m ix FT C P , m ix T T + U ⋅ d A ⋅ ( T T − T S ) +
+ ⎡⎣ ∆ H 1 ( − r1 ) + ∆ H 2 ( − r2 ) + ∆ H 3 ( − r3 ) + ∆ H 4 ( − r4 ) ⎤⎦ d V T
4U
ρ T , m ix FT C P , m ix T T i = ρ T , m ix FT C P , m ix T T + (TT − T S ) d V T +
D ia ., T u b e
+ ⎡⎣ ∆ H 1 ( − r1 ) + ∆ H 2 ( − r2 ) + ∆ H 3 ( − r3 ) + ∆ H 4 ( − r4 ) ⎤⎦ d V T
⎛ 4U
⎜ D ( T T − T S ) + ⎡⎣ ∆ H 1 ( − r1 ) + ∆ H 2 ( − r2 ) + ∆ H 3 ( − r3 ) + ∆ H 4 ( − r4 ) ⎤⎦ ⎞⎟
d TT
− =⎜ ia ., T u b e
⎟ (E .3 .3 .9 )
d VT ⎜ ρ T , m ix FT C P , m ix ⎟
⎜ ⎟
⎝ ⎠
For the Shell side, we also consider like plug flow reactor system, which there is temperature
and concentration profile. The volume of shell reactor is temperature and conversion
dependence. So, taking mass balance and energy balance at the steady state condition is to
model this system as following.
E-36
In = O u t + C o n v e rs io n + h o ld u p
dC C
FS C C3H 8 = FS C C3H 8 + ( − rC )dVS + VS 3H 8
V=Vo V =V o+dV 3H 8
dt
F o r p ro p a n e
d C C3H 8 ⎛ ( − r5 ) + ( − r6 ) ⎞
− = ⎜ ⎟ (E .3 .4 .1 )
dVS ⎝ FS ⎠
For ethylene
d C C2H 4 − ( − r5 )
− = (E .3 .4 .2 )
dVS FS
For methane
dC CH 4 − ( − r5 )
− = (E .3 .4 .3 )
dVS FS
For propylene
d C C3H 6 − ( − r6 )
− = (E .3 .4 .4 )
dVS FS
For hydrogen
dC H2 − ( − r6 )
− = (E .3 .4 .5 )
dVS FS
Energy balance:
In = O u t + h e a t tra n s fe r + h e a t g e n e ra te d + h e a t a c c .
ρ S F S C P , S T S i = ρ S F S C P , S T S − U d A (TT − T S ) +
+ ⎡⎣ ∆ H 5 ( − r5 ) + ∆ H 6 ( − r6 ) ⎤⎦ d V S
4U
ρ S FS C P ,S TSi = ρ S FS C P ,S TS − (TT − T S ) d V T +
D ia , T u b e
+ ⎡⎣ ∆ H 5 ( − r5 ) + ∆ H 6 ( − r6 ) ⎤⎦ d V s
⎛ − 4U d VT ⎞
⎜ (TT − T S ) ⋅ + ⎡⎣ ∆ H 5 ( − r5 ) + ∆ H 6 ( − r6 ) ⎤⎦ ⎟
− d TS D d V
= ⎜ ia , T u b e S
⎟ (E .3 .4 .6 )
dVS ⎜ ρ S FS C P ,S ⎟
⎜ ⎟
⎝ ⎠
From the mathematic model above, they are the series of ordinary differential equations that
consume time to solve precisely. In this conceptual design, the reactor will be simplified by
considering the main reaction, and the reactor both sides take place in gas phase
heterogeneously with fixed bed catalyst, that heat capacity of gas phase is usually less than
the liquid phase. Then, the heat exchange between shell and tube could be by heat conduction
of each side in order to retain its temperature. Therefore the reactions of both sides are
assumed as isothermal reactions. In addition to this, catalyst deactivation and coke formation
E-37
will not be considered during design as well. The following section shows the simplified
calculation of reactor both sides.
For the reactions in tube side, oxidation dehydrogenation of propane will be considered in the
main reaction, the kinetic rate expression and catalyst data show below
Catalyst data:
Type: V2O5
Specific surface: 95m2/g-cat
Void of bed: 0.5
Particle density: 1500kg/m3
Effective diffusivity of propane to V2O5 is less, and then assumes the catalyst effectiveness is
~1.
We can rewrite the rate reaction constant, which is catalyst mass basis into the volume of
fluid basis.
From the differential mass balance in terms of the void volume, the general expression shows
as
FA0 dX A = ( − rA'' ) dVP
( − rA'' )dVP = ( − rA ) ⋅ void fraction. ⋅ dVRe actor
so,
FA0 dX A = ( − rA ) ⋅ void fraction ⋅ dVRe actor
E-38
FC3 H8,0 dX C3 H8 = k1' ⋅ PC0.6
3 H8
⋅ void fraction ⋅ dVRe actor
Assume the gas in reaction perform like ideal gas, hence the term PC0.6
3 H8
can be changed in
= ( CC3 H8 RT )
0.6
term of concentration PC0.6
3 H8
Then,
( )
0.6
FC3 H8,0 dX C3 H8 = k1' ⋅ CC3 H8 RT ⋅ void fraction ⋅ dVRe actor
( )
0.6
= k1' ⋅ CC3 H8,0 (1 − X C3 H8 ) RT ⋅ void fraction ⋅ dVRe actor
X = 0.8 V
k1'
1
( )
0.6
In order to design as the base case, from ASPEN result, we can define
FC3 H8 ,0 = 534.121 kmol/hr
0.664 /1.01 kmol
CC3 H8 ,0 = ( )
R ⋅T m3
⎛ 0.657 ⎞ kmol
= ⎜ ⎟
⎝ R ⋅T ⎠ m
3
therefore,
X = 0.8 VReactor
1 50 ⋅ 3600
∫ ∫ ⋅ ( 0.657 ) ⋅ (1 − 0.5) ⋅ dVRe actor
0.6
dX C3 H8 =
(1 − X ) 534.121 ⋅1000
0.6
X =0 C3 H 8 V=0
E-39
⎛ -Ea i ⎛ 1 1 ⎞ ⎞
ki = k 0i exp ⎜⎜ ⎜ − ⎟ ⎟⎟ (E.3.6.5)
⎝ R ⎝ T T0 ⎠ ⎠
⎛ -∆H ⎛ 1 1 ⎞ ⎞
Ki = K 0i exp ⎜⎜ ⎜ − ' ⎟ ⎟⎟ (E.3.6.6)
⎝ R ⎝ T T0 ⎠ ⎠
The data from the literature [ Applied catalyst A :General 248(2003) 105-106, J.Gascon]
shows in Table E.3.1 below.
Catalyst specification:
Shape : Pellet
Diameter : 5.0mm
High : 4.9 mm
Specific surface : 214m2 /g-cat
Void fraction : 0.5
Catalyst density : 1,500kg/m3
E-40
PC3H6 PH2
0.0033(PC3H8 − )
(− r2 ) = 0.089 (E.3.6.8)
1 + PC3H6 6.8 ⋅10−7
The design equation for the plug flow is
x A = x AF
V dX A
= ∫
F xA =0
(-rA )
If the mass balances are based on propane, then the two equations are
V dX
=∫ 1 (E.3.6.9)
F (-r1 )
V dX
=∫ 2 (E.3.6.10)
F (-r2 )
where the conversion X1 is the gram moles of propane disappearing by the reaction 1 per
gram mole of the feed, and the conversion X2 is the gram moles of propane disappearing by
the reaction 2 per gram mole of the feed.
Based on 1.0 mole of feed of entering propane, the moles of each component at conversions
X1 and X2 are:
Component mole at conversions X1 ans X2
C3H8 1-X1-X2
C3H6 X2
H2 X2
CH4 X1
C2H4 X1
Total 1+X1+X2
Since the total mole equals to 1+X1+X2, the mole fractions of each components are also
given by these quantities. If the components are assumed to behave as ideal gases, then the
partial pressures are:
Component Partial pressure pi
C3H8 1-X1-X2 / (1+X1+X2)
C3H6 X2/ (1+X1+X2)
H2 X2/ (1+X1+X2)
CH4 X1 / (1+X1+X2)
C2H4 X1 / (1+X1+X2)
Substituting the partial pressure of the components of the components in eq. E.3.6.11 and eq.
E.3.6.12 gives
⎛ 1 − X1 − X 2 ⎞
(− r1 ) = 17.62 ⎜ ⎟ (E.3.6.11)
⎝ 1 + X1 + X 2 ⎠
⎛ ⎛ 1− X − X ⎞ X22 ⎞
0.0033 ⎜ ⎜ 1 2
⎟ − 2 ⎟
⎜ ⎝ 1 + X1 + X 2 ⎠ 0.089 (1 + X + X ) ⎟
(− r2 ) = ⎝ 1 2 ⎠ (E.3.6.12)
⎛ X2 ⎞
⎜1 + 6.8 ⋅10−7 ⎟
⎝ (1 + X1 + X 2 ) ⎠
E-41
Substituting eq. E.3.6.11 and eq. E.3.6.12 in the eq. E.3.6.13 and eq. E.3.6.14, and the
composition of the components are computed for various values of V/F.
dX1 ⎛ 1 − X1 − X 2 ⎞
= r1 (X1 ,X 2 )= 17.62 ⎜ ⎟ (E.3.6.13)
d(V/F) ⎝ 1 + X1 + X 2 ⎠
⎛ ⎛ 1− X − X ⎞ X22 ⎞
0.0033 ⎜ ⎜ 1 2
− ⎟
dX 2 ⎜ ⎝ 1 + X1 + X 2 ⎟⎠ 0.089 (1 + X + X )2 ⎟
= r2 (X1 ,X 2 ) = ⎝ 1 2 ⎠ (E.3.6.14)
d(V/F) ⎛ X2 ⎞
⎜1 + 6.8 ⋅10−7 ⎟
⎝ (1 + X1 + X 2 ) ⎠
Calculating by function ODE45, MATLAB, resulting in Table E.3.2 and Figure E.3.1 shows
the plot of the rates of each reaction as a function of V/F. In both instance, the rates decrease
toward zero as V/F increases.
Table E.3.2 Conversion versus V/F for the dehydrogenation of propane
V/F X1 X2 XT rate1 rate2
0.0000 0.0000 0.0000 0.0000 1.76E+01 0.0033
0.0018 0.0307 5.8E-06 0.0307 1.66E+01 0.0031
0.0036 0.0597 1.12E-05 0.0597 1.56E+01 0.0028
0.0054 0.087 1.63E-05 0.0870 1.48E+01 0.0025
0.0072 0.113 2.12E-05 0.1130 1.40E+01 0.0022
0.0126 0.1835 3.44E-05 0.1835 1.22E+01 0.0014
0.0180 0.245 4.59E-05 0.2450 1.07E+01 0.0006
0.0234 0.2994 5.61E-05 0.2995 9.50E+00 -0.0002
The reactor volume required to process from Aspen Plus Simulation, 889kmol/hr of propane
is estimated from Table E.3.2. For the total conversion of 23 %,
V/F = 0.018
V = 0.018F
= 0.018*889*1000/3600
= 4.45m3
Catalyst particle void in the reactor is 0.5, so the volume of catalyst in shell side is equal to
2.225 m3 and catalyst used 3.3 tons.
E-42
Rates of reaction for propane dehydrogenation
2.00E+01
1.80E+01
1.60E+01
1.40E+01
Rates of reaction ( mole/m3 s)
1.20E+01
1.00E+01
Rate1
Rate2
8.00E+00
6.00E+00
4.00E+00
2.00E+00
0.00E+00
0.0000 0.0018 0.0036 0.0054 0.0072 0.0126 0.0180 0.0234
-2.00E+00
V/F (m3/mol s)
Since, in this design stage, we use the simplified model, the heat transfer coefficient of shell
and tube is assumed to a general figure approximately. The heat transfer coefficient used in
further calculation like cost evaluation is assumed as 500 W/m2 oC. The heat exchange
between shell and tube is roughly equal to the heat reaction of exothermic, which will be
removed, as mention is Appendix A.2. The amount of heat is 20kcal/molC3H8 ( 84 kJ/molC3H8 )
and feed propane through the exothermic reaction is 804.63 kmol/hr based on the required
product 200 kton annual. Therefore the removed heat is 18.7 x 103 kW. Therefore the total
area of the tube can be estimated.
Q = UA∆TLM
Q 18.7 ⋅103 kW
A= = = 121m 2
U ∆TLM 500W / m C ⋅ (850 − 540)
2 D
Then, the ratio of area and volume results in the diameter of the tube side 0.3 m.
The number and length of tube can be calculated listed in the Table E.3.3. We choose 20
tubes with the 6.41-meter long.
Table E.3.3. Number and length of tube reacto
Dia. Tube Numbe of tube Height of Tube Area Volume
(0.3 m) (m) (m2) (m3)
0.3 10 12.82 121 9.06
0.3 20 6.41 121 9.06
0.3 30 4.27 121 9.06
0.3 40 3.20 121 9.06
0.3 50 2.56 121 9.06
0.3 60 2.14 121 9.06
0.3 70 1.83 121 9.06
0.3 80 1.60 121 9.06
0.3 90 1.42 121 9.06
0.3 100 1.28 121 9.06
E-43
0.3 110 1.17 121 9.06
0.3 120 1.07 121 9.06
And for the shell side, the total volume can estimate the diameter of reactor.
Total Volume = 4.45 + 9.06 = 13.51m3
4 ⋅ Volume total
Dia.ofReactor = = 1.63 m.
π ⋅ height
E-44
Appendix E.4 Shell and Tube Heat Exchanger Design
The design of shell and tube heat exchanger involves the determination of the heat transfer
coefficient and pressure drop on both the tube side and the shell side. A large number of
methods are available for determining the shellside performance. Before the design procedure,
some guidelines for shellside design and points to be raised while specifying a heat exchanger
are listed, followed by preliminary sizing of a shell and tube heat exchanger.
Guideline for Shell-Side Design
Recommend guidelines for shellside design include:
1. Accept TEMA fabrication clearances and tolerances, and enforce these standards during
fabrication.
2. For segmental baffles employ 20%battle cuts.
3. Employ no-tubes-in-window design to eliminate the damage from flow-induced vibration.
4. Evaluate the heat transfer in the clean condition and pressure drop in the maximum fouled
condition.
5. Employ sealing devices to minimize bypassing between the bundle and shell for pull
through floating heat exchanger, and through pass partition lanes.
6. Ration of baffle spacing to shell diameter may be restricted between 0.2 and 1.0. Baffle
spacing much greater than the shell diameter must be carefully evaluated.
7. Avoid shell longitudinal baffles that are not welded to the shell; all other sealing methods
are inadequate.
E-45
Design Considerations for a Shell and Tube Heat Exchanger
The basic criterion that a given or designed heat exchanger should satisfy is that it should
perform the given heat duty within the allowable pressure drop. The design is also to satisfy
additional criteria such as:
1. Withstand operating conditions, startup, shutdown, and upset conditions that influence the
thermal and mechanical design.
2. Maintenance and servicing.
3. Multiple shell arrangement.
4. Cost.
5. Size limitations.
In terms of five factors just mentioned, multishell arrangement needs some comments on it.
Consider the advantages of a multishell arrangement to allow one unit to be taken out of
service for maintenance without severely upsetting the rest of the plant. For part load
operation, multiple shells result into an economical operation. Shipping and handling may
dictate restrictions on the overall size or weight of the unit, resulting in multiple shells for an
application.
E-46
Problem Identification
Acceptable
Figure E.4.1 Basic logic structure for process heat exchanger design
Example of Design
Here present the design of one shell and tube heat exchanger in dehydrogenation process. For
the physical properties and conditions of the components in shell and tube is shown as Table
E.4.1.
Table E.4.1. Physical properties and conditions of components
Shell- side Inlet Mean Outlet
Stream No. 201
Phase Gas Gas Gas
Temperature oC 40 59 78
Specific heat kJ/kg oC 2.01 2.05 2.09
Thermal conductivity W/m oC 0.135 0.134 0.133
Density kg/m3 840 820 800
Viscosity mNsm-2 4.3 3.2 2.4
Tube-side Inlet Mean Outlet
Stream No. 202
Phase Gas Gas Gas
E-47
Temperature oC 200 145 90
Specific heat kJ/kg oC 2.72 2.47 2.26
Thermal conductivity W/m oC 0.130 0.132 0.135
Density kg/m3 690 730 770
Viscosity mNsm-2 0.22 0.43 0.80
E-48
⎡ 1− S ⎤
(R 2
+ 1) ln ⎢
⎣1 − RS ⎥⎦
F= = 0.876
(
⎡ 2 − S R + 1 − R2 + 1 ⎤
( R − 1) ln ⎢⎢ ⎥ )
(
⎢⎣ 2 − S R + 1 + R + 1 ⎥⎦
2 ⎥
)
Method to Determine Number of Shells
Quickly check the limits:
2Tho ≥ Tci + Tco hot fluid on the shell side
2Tco ≤ Thi + Tho cold fluid on the shell side
If the limits are approached, it is necessary to use multiple 1-2N shells in series. There is a
rapid graphical technique for estimating a sufficient number of 1-2N shells in series. The
terminal temperatures of the two streams are plotted on the ordinates of an arithmetic graph
paper sheet as shown as Figure E.4.2:
Temperature
Heat Transfer
Figure E.4.2 Heat transfer efficiency v.s. temperature
The distance between the ordinates is arbitrary. Starting with the cold fluid outlet temperature,
a horizontal line is laid off until it intercepts the hot fluid line. From that point a vertical line
is draw to the cold line at or below the cold fluid inlet temperature. The number of horizontal
lines (including the one that intersects the right-hand ordinate) is equal to the number of
shells in series that is clearly sufficient to perform the duty. Following this procedure will
usually result in a number of shells having an overall F close to 0.8.
Estimation of U
First, need to guess a number of U. The assumption can be chosen from Table E.4.2. [Coulson &
Richardson’s Volume 6]
. Here start with 300 W/m2oC
Table E.4.2 Typical overall coefficients
E-49
E-50
Fouling Resistance
Fouling resistance values may be chosen from TEMA Table E.4.3 [Coulson & Richardson’s Volume 6].
Table E.4.3 Fouling factors (coefficients), typical values
Calculation of Ao
Once q, U, LMTD, and F are known, the total outside heat transfer area Ao can be calculated:
q
Ao = (E.4.3)
UF ( LMTD )
1509.4 ×103
Ao = = 71.15 m 2
300 × 0.876 × 80.7
E-51
Determination of Shell Size and Tube Length from Heat Transfer Area, Ao
The objective is to find the right number of tubes of diameter, Do, and the shell diameter, Ds,
to accommodate the number of tubes, Nt, with given tube length, L.
Ao = π Do N t L (E.4.4)
The problem now arises of how to interpret the value of value Ao in terms of tube length and
shell diameter, when both values are not known.
For estimation purpose, the tube number Nt is given by
⎛ CTP ⎞ Ds2
N t = 0.785 ⎜ ⎟ (E.4.5)
⎝ CL ⎠ ( PR ) Do2
2
From equation (27), an expression for the shell diameter in terms of main constructional
diameters can be obtained as:
1
CL ⎡ Ao ( PR ) Do ⎤
2 2
Ds = 0.637 ⎢ ⎥ (E.4.7)
CTP ⎢⎣ L ⎥⎦
Heat transfer surface Ao can be obtained by various combinations of parameter L and Do for
any given tube layout pattern. An initial assumption of these values is necessary.
Here use 19.05 mm outside diameter, 14.83 mm inside diameter, and 5 m long tubes on a
square pitch arrangement.
Area of one tube:
A1 = π × 0.01905 × 5 = 0.2992 m 2
Number of tubes = 71.15/0.2992 = 237, say 240
So, for 2 passes, tubes per pass, Nt = 120
1
CL ⎡ Ao ( PR )2 Do ⎤ 2
1 ⎡ 71.15 × (1.25 ) × 0.01905 ⎤
2
E-52
hi d e
Nu = Nusselt number =
kf
ρ ut d e
Re = Reynolds number =
µ
Cpµ
Pr = Prandtl number =
kf
C = 0.021 for gases
C = 0.023 for non-viscous liquids
C = 0.027 for viscous liquids
and
hi = inside cofficient, W/m°C,
d e = equivalent diameter, m = Di
ut = fluid velocity, m/s,
k f = fluid thermal conductivity, W/m°C,
Gt = mass velocity, mass flow per unit area, kg/m 2s,
µ = fluid viscosity at the bulk fluid temperature, Ns/m 2 ,
µ w = fluid viscosity at the wall,
C p = fluid specific heat, heat capacity, J/kg°C.
Calculation of ut:
π
( 0.01483)
2
tube cross-sectional area = = 0.0001727 m 2
4
so area per pass = 120 × 0.0001727 = 0.02073 m 2
19.44 3
volumetric flow rate = = 0.0237 m
820 s
0.0237
tube-side velocity, ut = = 1.14 m
0.02073 s
It is often convenient to correlate heat transfer data in terms of a heat transfer j factor. The
heat transfer factor is defined as:
−0.14
⎛ µ ⎞
jh = St Pr 0.67
⎜ ⎟
⎝ µw ⎠ (E.4.9)
St = Stanton number = 0.0225exp ⎡ −0.0225 ( ln Pr ) ⎤ Re−0.205 Pr −0.505
2
⎣ ⎦
eq. E.4.8 can be rearranged as:
0.14
hi Di ⎛ µ ⎞
= jh Re Pr 0.33 ⎜ ⎟ (E.4.10)
kf ⎝ µw ⎠
After having Re and L/Di, from Figure E.4.3 [Coulson & Richardson’s volume 6], it can get the jh value.
By eq. E.4.10, hi can be calculated.
Calculation of hi:
E-53
820 × 1.15 × 0.01483
Re = = 4388
3.2 × 10−3
2.996 × 103 × 3.2 × 10−3
Pr = = 71.57
0.134
L 5
= = 337
Di 0.01483
From Figure E.4.3,
jh = 3.3 ×10−3
N u = 3.3 ×10−3 ( 4388 )( 71.57 )
0.33
= 59.27
⎛ 0.134 ⎞ W 2o
hi = 59.27 × ⎜ ⎟ = 536
⎝ 0.01483 ⎠ m C
E-54
FigureE.4.3 Tube side heat transfer factor
E-55
Pt = tube pitch
Do = tube outside diameter
Ds = shell inside diameter
lB = baffle spacing
The term (Pt-Do)/Pt is the ratio of the clearance between tubes and the total distance between
tube centers. The baffle spacing used arrange from 0.2 to 1.0 shell diameter. A close baffle
spacing will give higher heat transfer coefficient but at the expense of higher pressure drop.
The optimum spacing will usually be between 0.3 to 0.5 times the shell diameter. The baffle
cut is the height of the segment removed to form the baffle, expressed as a percentage of the
baffle disc diameter. Generally, a baffle cut of 20 to 25% will be the optimum, giving good
heat transfer rate, without excessive pressure drop.
(1.25 × 0.01905 − 0.01905) 0.437
As = × 0.437 × = 0.009548 m 2
1.25 × 0.01905 5
Calculate the shell-side mass velocity Gs and the linear velocity us:
W
Gs = s
As
Gs
us =
ρ
where
Ws = fluid flow rate on the shell-side, kg/s
ρ = shell-side fluid density, kg/m3
Calculate the shell-side equivalent diameter (hydraulic diameter). For a square pitch
arrangement:
1.27 2
de =
Do
( Pt − 0.785 Do2 ) (E.4.12)
E-56
de =
1.27 2
Do
( Pt − 0.785Do2 ) =
1.27
0.01905
(
(1.25 × 0.01905) − 0.785 ( 0.01905)o = 0.01881 m
2 2
)
5.556 3
volumetric flow rate = = 0.00761 m
730 s
0.0761
shell-side velocity = = 0.797 m From
0.009548 s
730 × 0.797 × 0.01881
Re = = 25453
0.43 × 10−3
3.515 ×103 × 0.43 ×10−3
Pr = = 11.45
0.132
Figure E.4.4,
jh = 6 × 10−3
⎛ 0.132 ⎞
⎟ × 6 ×10 × 25453 × (11.45 ) = 2396
ho = ⎜ −3 0.33
W 2o
⎝ 0.0188 ⎠ m C
E-57
can be modified as:
1
Uo =
1 ⎛ Do ⎞ 1 ⎛ Do ⎞ Do ⎛ Do ⎞ 1 1
⎜ ⎟+ ⎜ ⎟+ ln ⎜ ⎟+ +
hdi ⎝ Di ⎠ hi ⎝ Di ⎠ 2km ⎝ Di ⎠ ho hdo
1
=
1 ⎛ 0.01905 ⎞ 1 ⎛ 0.01905 ⎞ 0.01905 ⎛ 0.01905 ⎞ 1 1
⎜ ⎟+ ⎜ ⎟+ ln ⎜ ⎟+ +
0.0002 ⎝ 0.01483 ⎠ 536 ⎝ 0.01483 ⎠ 2 × 45 ⎝ 0.01483 ⎠ 2396 0.0002
= 300.6921 W 2 o
m C
Check the calculation:
U − U o ,ass
If 0 < o ,calc < 30% , the initial assumption of Uo is acceptable.
U o , ass
U o ,calc − U o ,ass 300.6921 − 300
×100% = = 0.23%
U o ,ass 300
Therefore, the initial assumption of Uo is acceptable.
E-58
Appendix E.5 T302 Column sizing report
Table E.5.1 T302 column sizing report
Item Sizing Data Comments
Parameter
MG 10.9 kg/s Gas phase flow rate
RhoG 21.7 kg/m3 Gas phase density
ML 22.05 kg/s Liquid phase flow rate
RhoL 426 kg/m3 Liquid phase density
ViscL 1.00E-04 Pa s Liquid phase dynamic viscoscity
Capacity Factor Correlation
FLG 0.455 [-] ( ML /MG )( RhoG/RhoL)0.5
Design Pressure drop 4.1 mbar/m distillation range 4-8 mbar/m
Crp 0.65 [-] Read from the Figure 8.2.3.1
Frp 66 1/m use dia.50 mm pall ring in order to saving cost
FG,oper 1.1874 [-] [Crp*(RhoL-RhoG)/(13*Frp(viscL/RhoL)0.1]0.5
Dimension of Column
Dia 1.586 m [4*MG/(pi*FGOper*Rho0.5)]
This valid for Rasching and pall ring in distillation
HETP 0.85 m with moderate value of surface tension .
low viscoscity and pressure drop below 6 mbar/m
No. of ideal stages 24 [-]
Height of packing 20.4 m
Height of top 1 m
Height of bottom 1 m
Height of skirt 1.5 m
Height of column 23.9 m
E-59
Appendix E.6 Calculation of CO2 removal equipment
T401 CO2 absorber
From CO2 solubility data are given as blow2:
2
http://www.cape.canterbury.ac.nz/Apcche_Proceedings/ APCChE/Data/802rev.pdf
E-60
90
80
70
CO2 pressure, kPa 60
25oC
50
50oC
40
75oC
30
20
10
0
0.0000 0.1000 0.2000 0.3000 0.4000 0.5000 0.6000
CO2_loading, mol CO2/mol MDEA
Due to the volatility and related high concentration of MDEA and followed low temperature
separation process, 25oC absorption temperature is designed.
Partial pressure of CO2 in the feed = 1.5%*101.325 = 1.52[kPa]
Partial pressure of CO2 in the exit gas at 99 percent recovery =1.52*0.01= 0.0152 [kPa]
Assume the equilibrium line at 25oC in Figure E.6.1 is straight line so Figure E.6.2 [volume 6 p597]
can be used to estimate the number of stages needed.
Figure E.6.2 Number of transfer units NOG as a function of y1/y2 with mGm / Lm as parameter
From the data above: partial pressure at 0.2988 [mol CO2/mol MDEA] = 10.34 [kPa].
E-61
10.34
Mol. fraction in vapor = = 0.1020
101.325
0.2988
Mol. fraction in liquid = = 0.0377
0.2988 + 1 + 1*119.2 /18
0.1020
The slope of operation line m = = 2.71
0.0377
Using Figure E.6.2 the number of stages required at different MDEA solution flow rates will
be determined and the “optimum” rate chosen:
y1 p1 14.54
= = = 100
y2 p2 0.1454
where y1 and y2 = the mol fractions of the solute in the gas at the bottom and top of the
column, respectively
Gm
m 0.1 0.2 0.3 0.4 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9
Lm
N OG 5.0 5.4 6.1 6.7 7.8 8.4 9.2 10.2 11.3 13.0 15.1 18.3 23.6
where Gm = molar gas flow-rate per unit cross-sectional area
Lm = molar liquid flow-rate per unit cross-sectional area
N OG = the number of overall gas-phase transfer units
Below 0.55 there is only a small decrease in the number of stages required with increasing
liquid rate; above 0.7 the number of stages increases rapidly with decreasing liquid rate. It
can be seen that the “optimum” will be between mGm / Lm =0.55 to 0.7, as would be expected.
Check the liquid outlet composition at 0.55 and 0.7:
Material balance Lm x1 = Gm ( y1 − y2 )
G m Gm
So x1 = m (1.5% *0.99) = (0.0148)
Lm 2.71 Lm
G
at m m = 0.55, x1 = 0.003 mol fraction
Lm
G
at m m = 0.7, x1 = 0.0038 mol fraction
Lm
Use 0.7, as the higher concentration will favor the stripper design and operation, without
significantly increasing the number of stages needed in the absorber.
N OG = 11.3 , =12 (say)
Column diameter
From the data of mass balance,
Gas flow-rate (stream 209) = 7.64kg/s, = 0.263kmol/s
Carbon dioxide mole fraction in gas inlet is 0.015.
Liquid flow-rate = 2.71/0.7*0.263 = 1.02kmol/s, = 31.91kg/s
Select 50mm metal Pall rings.
From Table E.6.2 [volume 6], Fp = 66m-1
Gas density ρV = 27498.122/667.676 = 41.18kg/m3
E-62
Assume: Liquid density ρ L = 1000kg/m3
Liquid viscosity µL = 10-3Ns/m2
L* ρV 69.97 41.18
FLV = W* = = 1.86
VW ρ L 7.64 1000
where L*W = liquid mass flow-rate per unit column cross-sectional area, kg/m2s
VW* = gas mass flow-rate per unit column cross-section area, kg/m2s
FLV = flow factor
Design for a pressure drop of 20mm H2O/m packing
From Figure E.6.3 [volume 6]
K4 = 0.22
At flooding, K4 = 0.42
Here K4 is the function [volume 6]:
13.1(VW* ) 2 FP ( µ L / ρ L )
0.1
K4 =
ρV ( ρ L − ρV )
where FP = packing factor, characteristic of the size and type of packing, m-1
2
µL = liquid viscosity, Ns/m
3
ρL, ρV = liquid and vapor densities, kg/m
E-63
Figure E.6.3Generalized pressure drop correlation, adapt from a figure by Norton Co. with permission
0.22
Percentage flooding = = 0.72 , satisfactory.
0.42
From equation above
0.5
⎡ K ρ (ρ − ρ ) ⎤
V =⎢ 4 V L
* V
0.1 ⎥
⎢⎣13.1FP ( µ L / ρ L ) ⎥⎦
W
0.5
⎡ ⎤
⎢ 0.22* 41.18(1000 − 41.18) ⎥
=
⎢ 13.1*66 (10-3 /1000 )0.1 ⎥
⎣ ⎦
= 6.32kg/m 2s
Column height
Onda’s method
The equation for the effective area is:
E-64
⎡ 0.75 0.1
⎛ L*W 2 a ⎞
−0.05
⎛ L*W 2 ⎞ ⎤
0.2
aw ⎛ σ ⎞ ⎛ L*
⎞
= 1 − exp ⎢ −1.45 ⎜ C
⎟ ⎜
W
⎟ ⎜ ⎟ ⎜ ⎟ ⎥
a ⎢ ⎝ σ L ⎠ ⎝ aµ L ⎠ ⎝⎜ ρ L 2 g ⎠⎟ ⎝⎜ ρ Lσ L a ⎠⎟ ⎥
⎣ ⎦
and for the mass coefficients:
1/ 3 2/3 −1/ 2
⎛ ρ ⎞ ⎛ L* ⎞ ⎛ µL ⎞
kL ⎜ L ⎟ = 0.0051⎜ W ⎟ ⎜ ⎟ (ad P )0.4
⎝ µL g ⎠ ⎝ aw µ L ⎠ ⎝ ρ L DL ⎠
0.7 1/ 3
kG RT ⎛ V* ⎞ ⎛ µ ⎞
= K 5 ⎜ W ⎟ ⎜ V ⎟ (ad P ) −2.0
a Dv ⎝ a µV ⎠ ⎝ ρV Dv ⎠
where K5 = 5.23 for packing sizes above 15mm, and 2.00 for sized below 15mm
aw = effective interfacial area of packing per unit volume, m2/m3
a = actual area of packing per unit volume m2/m3
dP = packing size, m
σC = critical surface tension for the particular packing material, N/m [data
see volume 6]
σL = liquid surface tension, N/m
kG = gas film mass transfer coefficient, kmol/m2s atm or kmol/m2s bar
kL = Liquid film mass transfer coefficient, kmol/m2s (kmol/m3)=m/s
R = gas constant, bar m3/kmol K
T = operation temperature, K
g = gravity, 9.81m/s2
Dv, Dl= Mass diffusivity in vapor or liquid phase, m2/s
Gas and liquid diffusivities can be calculated using given equation as below [volume 6,
chapter 8, p291]
1.013e-7T 1.75 (1/ M a + 1/ M b )1/ 2
Dv =
1/ 3 2
⎡⎛ ⎞
1/ 3
⎛ ⎞ ⎤
P ⎢⎜ ∑ vi ⎟ + ⎜ ∑ vi ⎟ ⎥
⎣⎢⎝ a ⎠ ⎝ b ⎠ ⎦⎥
where Ma,Mb = molecular weights of components a and b
∑a vi , ∑b vi = the summation of the special diffusion volume coefficients for
components a and b, given in Table E.6.3 [volume 6]
1.173e-13(φ M )0.5 T
DL =
µVM 0.6
where φ = an association factor for the solvent
2.6 for water
1.9 for methanol
1.5 for ethanol
1.0 for unassociated solvents
M = Molecular weight of solvent
µ = viscosity of solvent, mNs/m2
Vm = molar volume of the solute at its boiling point, m3/kmol. This can be
estimated from the group contributions given in Table E.6.4 [volume 6]
E-65
Table E.6.4 structural contributions to molar volumes, m3/mol
1.173e-13*(1.0*31.28)0.5 *303.13
DL =
1.002*0.03400.6
= 1.51e-9m 2 /s
The film transfer unit heights are given by:
Gm
HG =
kG a w P
Lm
HL =
k L awCt
where P = column operating pressure, bar
E-66
Ct = total concentration, kmol/m3
HL, HG= overall height of a transfer unit and the individual film transfer units
Assume the flow of gas and liquid is essentially constant throughout the column, the height of
packing required, Z, is given by:
Gm y1 dy
Z= ∫
K G aCt y 2 y − ye
Z = H OG N OG
G
H OG = H G + m m H L
Lm
where Z = height of packing required, m
ye = concentration in the gas that would be in equilibrium with the liquid
concentration at any point
HOG = height of an overall gas-phase transfer unit
L*W = 31.91/1.33 = 23.99kg/m 2s
R= 0.08314bar m3/kmol K
Assume σL equals water = 70e-3N/m
g = 9.81m/s2
dP = 50e-3m
From Table E.6.2, for 50mm Pall rings
a =102m2/m3
For metal materials,
σC = 75e-3N/m
aw ⎡ 0.75 0.1
⎛ 75e-3 ⎞ ⎛ 23.99 ⎞ ⎛ 23.99 *102 ⎞
2 −0.05
⎛ 23.992 ⎞ ⎤
0.2
= 1 − exp ⎢ −1.45 ⎜ ⎟ ⎜ ⎟ ⎜ 2 ⎟ ⎜ ⎟ ⎥
a ⎢⎣ ⎝ 70e-3 ⎠ ⎝ 102*1e-3 ⎠ ⎝ 1000 *9.81 ⎠ ⎝ 1000*70e-3*102 ⎠ ⎥⎦
= 0.87
aw = 0.87 *102 = 88.74m 2 /m3
1/ 3 2/3 −1/ 2
⎛ 1000 ⎞ ⎛ 52.61 ⎞ ⎛ 10−3 ⎞
k L ⎜ −3 ⎟ = 0.0051⎜ −3 ⎟ ⎜ ⎟ (102*50e-3)0.4
⎝ 10 *9.81 ⎠ ⎝ 88.74*10 ⎠ ⎝ 1000*1.51e-9 ⎠
k L = 5.88e-4m/s
VW* on actual column diameter = 7.64/1.33 = 5.74kg/m2s
0.7 1/ 3
kG 0.08314*303.15 ⎛ 5.74 ⎞ ⎛ 0.018e-3 ⎞ −2.0
= 5.23 ⎜ ⎟ ⎜ ⎟ (102*50e-3)
102 5.45e-7 ⎝ 102*0.018 e -3 ⎠ ⎝ 41.18*5.45 e -7 ⎠
kG = 1.15e-4kmol/sm bar 2
5.74
Gm = = 0.198kmol/m 2s
27498.122 / 947.059
52.61
Lm = = 1.68kmol/m 2s
0.1312*119.2 + (1 − 0.1312) *18
0.198
HG = = 0.65m
1.15e-4*88.74*30
1.68
HL = = 1.08m
5.46e-4*88.74*(1000 / 31.28)
H OG = 0.65 + 0.7 *1.08 = 1.4m
Z = 1.4*12 = 16.8m
E-67
Round up packed bed height to 17m.
1.32
Then the column volume is deduced: V = π *17=22.6m3
4
After the calculation:
Column diameter: 1.30m
Column height: 17m
Column volume: 22.6m3
E-68
Appendix E.7 Gas-liquid separators calculation
E-69
2. D-301 T301 Reflux Drum
From data of Aspan simulation, at 15bar, -131.44oC, liquid density 502.50kg/m3, vapor
density 14.22kg/m3.
ut = 0.07 ⎡⎣( 502.5 − 14.22 ) /14.22 ⎤⎦
1/ 2
= 0.4102 m/s
Vapor volumetric flow-rate is 5275kg/hr =261.9m3/hr
Take hv = Dv / 2 and Lv / Dv = 3
π Dv2
Cross-sectional area for vapor flow = *0.5 = 0.393Dv2
4
Vapor velocity, uv = 261.9 / 3600 /(0.393Dv2 ) = 0.185Dv−2
Vapor residence time required for the droplets to settle to liquid surface
= hv / ut = 0.5Dv / 0.4102 = 1.22 Dv
Actual residence time = vessel length/vapor velocity
= Lv / uv = 3Dv / 0.185Dv−2 = 16.2 Dv3
For satisfactory separation required residence time = actual.
So, 1.22 Dv = 16.2 Dv3
Dv = 0.42 m , say 0.45m.
Liquid hold-up time,
Liquid volumetric flow-rate is 45810kg/hr =91.2m3/hr.
liquid cross-sectional area = π *(0.45)2 / 4*0.5 = 0.0795 m 2
Length, Lv = 3*0.45 = 1.35 m
Hold-up volume = 0.0795*1.5 = 0.1073 m3
Hold-up time = liquid volume/liquid flow-rate
= 0.1073 /(91.2 / 3600) = 4.23 s << 10 minutes
As same as D101 design, increase the vessel diameter to reach higher residence time. After
several trying, new Dv is given.
New Dv = 2.4 m
Recalculate:
Liquid hold-up time,
Liquid volumetric flow-rate is 91.2m3/hr.
liquid cross-sectional area = π *(2.4) 2 / 4*0.5 = 2.26 m 2
Length, Lv = 3* 2.4 = 7.2 m ,
Hold-up volume = 2.26*7.2 = 16.272 m3
Hold-up time = liquid volume/liquid flow-rate
= 16.272 /(91.2 / 3600) = 642.3 s = 10.7 minutes
This is satisfactory.
D2 (2.4) 2
The vessel volume is π v Lv = π * *7.2 = 32.6 m3
4 4
Result:
Vessel diameter: 2.4m.
Vessel length: 7.2m.
Vessel volume: 32.6m3.
E-70
Appendix E.8 Equipment Summary & Specification Sheets
E-71
REACTORS, COLUMNS & VESSELS – SUMMARY
E-72
REACTORS, COLUMNS & VESSELS – SUMMARY
E-73
REACTORS, COLUMNS & VESSELS – SUMMARY
E-74
HEAT EXCHANGERS & FURNACES – SUMMARY
E-75
HEAT EXCHANGERS & FURNACES – SUMMARY
E-76
HEAT EXCHANGERS & FURNACES – SUMMARY
E-77
HEAT EXCHANGERS & FURNACES – SUMMARY
E-78
HEAT EXCHANGERS & FURNACES – SUMMARY
E-79
HEAT EXCHANGERS & FURNACES – SUMMARY
E-80
PUMPS, BLOWERS & COMPRESSORS – SUMMARY
E-81
PUMPS, BLOWERS & COMPRESSORS – SUMMARY
E-82
PUMPS, BLOWERS & COMPRESSORS – SUMMARY
E-83
Specification Sheets
E-84
SHELL & TUBE REACTOR – SPECIFICATION SHEET
Remarks:
(1) Novel reactor type
E-85
VESSEL – SPECIFICATION SHEET
Type : - Drum
Position : - Horizontal
- Vertical
Internals : - Demister / Plate / Coil / _________
Heating/Cooling medium : - none / Open / Closed / External Hxgr /________
- Type : n.a.
- Quantity [kg/s] : n.a.
- Press./Temp.’s [bara/oC] : n.a.
Vessel Diameter (ID) [m] : 1.5
Vessel Height [m] : 6
Vessel Tot. Volume [m3] : 10.6
Vessel Material : Carbon steel
Other :
Process Conditions
Stream Data Feed Top Bottom
E-86
VESSEL – SPECIFICATION SHEET
Type : - Drum
Position : - Horizontal
- Vertical
Internals : - Demister / Plate / Coil / _________
Heating/Cooling medium : - none / Open / Closed / External Hxgr /________
- Type : n.a.
- Quantity [kg/s] : n.a.
- Press./Temp.’s [bara/oC] : n.a.
Vessel Diameter (ID) [m] : 3
Vessel Height [m] : 6
Vessel Tot. Volume [m3] : 42.4
Vessel Material : Carbon steel
Other :
Process Conditions
Stream Data Feed 1 Feed 2 Out
E-87
VESSEL – SPECIFICATION SHEET
Type : - Drum
Position : - Horizontal
- Vertical
Internals : - Demister / Plate / Coil / _________
Heating/Cooling medium : - none / Open / Closed / External Hxgr /________
- Type : n.a.
- Quantity [kg/s] : n.a.
- Press./Temp.’s [bara/oC] : n.a.
Vessel Diameter (ID) [m] : 3.2
Vessel Height [m] : 6.4
Vessel Tot. Volume [m3] : 51.5
Vessel Material : Carbon steel
Other :
Process Conditions
Stream Data Feed Outlet
Temperature [oC] 25 25
Pressure [bara] 1 1
Density [kg/m3] 1.2 1.2
Mass Flow [kg/s] 7.7 7.7
Composition mol% wt% mol% wt%
C3H8 0.115 0.18 0.115 0.18
N2 0.011 0.01 0.011 0.01
C2H4 0.304 0.30 0.304 0.30
C3H6 0.218 0.32 0.218 0.32
CH4 0.133 0.07 0.133 0.07
H2 0.098 0.01 0.098 0.01
H2O 0.024 0.01 0.024 0.01
CO 0.083 0.08 0.083 0.08
CO2 0.015 0.02 0.015 0.02
Remarks:
Type : Drum
Position : - Horizontal
- Vertical
Internals : - Demister / Plate / Coil / _________
Heating/Cooling medium : - none / Open / Closed / External Hxgr /________
- Type : n.a.
- Quantity [kg/s] : n.a.
- Press./Temp.’s [bara/oC] : n.a.
Vessel Diameter (ID) [m] : 2.4
Vessel Height [m] : 7.2
Vessel Tot. Volume [m3] : 32.6
Vessel Material : Carbon steel
Other :
Process Conditions
Stream Data Feed Outlet Reflux
E-89
VESSEL – SPECIFICATION SHEET
Type : - Drum
Position : - Horizontal
- Vertical
Internals : - Demister / Plate / Coil / _________
Heating/Cooling medium : - none / Open / Closed / External Hxgr /________
- Type : n.a.
- Quantity [kg/s] : n.a.
- Press./Temp.’s [bara/oC] : n.a.
Vessel Diameter (ID) [m] : 1.3
Vessel Height [m] : 2.6
Vessel Tot. Volume [m3] : 3.45
Vessel Material : Carbon steel
Other :
Process Conditions
Stream Data Feed Outlet Reflux
E-90
VESSEL – SPECIFICATION SHEET
Type : -
Position : - Horizontal
- Vertical
Internals : - Demister / Plate / Coil / _________
Heating/Cooling medium : - none / Open / Closed / External Hxgr /________
- Type : n.a.
- Quantity [kg/s] : n.a.
- Press./Temp.’s [bara/oC] : n.a.
Vessel Diameter (ID) [m] : 1.3
Vessel Height [m] : 2.6
Vessel Tot. Volume [m3] : 3.45
Vessel Material : Carbon steel
Other :
Process Conditions
Stream Data Feed Bottom
Temperature [oC] 30 30
Pressure [bara] 30 30
3
Density [kg/m ] 656.45 656.45
Mass Flow [kg/s] 70.0 70.0
E-91
DISTILLATION COLUMN – SPECIFICATION SHEET
E-92
DISTILLATION COLUMN – SPECIFICATION SHEET
E-93
DISTILLATION COLUMN – SPECIFICATION SHEET
E-94
ABSORBER COLUMN – SPECIFICATION SHEET
E-95
DISTILLATION COLUMN – SPECIFICATION SHEET
Column Internals
Trays Number of Packing Type : pall ring
caps / sieve holes / : … Material : Plastic
Active Tray Area [m2] : … Volume [m3] :
Weir Length [mm] : … Length [m] :
Diameter of Width [m] :
chute pipe / hole / [mm] : … Height [m] : 4.25
Remarks:
E-96
HEAT EXCHANGER – SPECIFICATION SHEET
E-97
HEAT EXCHANGER – SPECIFICATION SHEET
E-98
HEAT EXCHANGER – SPECIFICATION SHEET
(1) Use during start up and used for temperature controllability during normal operation, target
temperature is 540 C
E-99
HEAT EXCHANGER – SPECIFICATION SHEET
E-101
HEAT EXCHANGER – SPECIFICATION SHEET
E-102
HEAT EXCHANGER – SPECIFICATION SHEET
E-103
HEAT EXCHANGER – SPECIFICATION SHEET
E-105
HEAT EXCHANGER – SPECIFICATION SHEET
E-106
HEAT EXCHANGER – SPECIFICATION SHEET
E-109
HEAT EXCHANGER – SPECIFICATION SHEET
E-110
HEAT EXCHANGER – SPECIFICATION SHEET
E-113
HEAT EXCHANGER – SPECIFICATION SHEET
E-115
HEAT EXCHANGER – SPECIFICATION SHEET
E-117
HEAT EXCHANGER – SPECIFICATION SHEET
E-118
HEAT EXCHANGER – SPECIFICATION SHEET
E-121
HEAT EXCHANGER – SPECIFICATION SHEET
E-125
CENTRIFUGAL PUMP – SPECIFICATION SHEET
E-127
CENTRIFUGAL PUMP – SPECIFICATION SHEET
E-128
CENTRIFUGAL PUMP – SPECIFICATION SHEET
E-129
CENTRIFUGAL PUMP – SPECIFICATION SHEET
E-130
CENTRIFUGAL PUMP – SPECIFICATION SHEET
E-131
CENTRIFUGAL PUMP – SPECIFICATION SHEET
E-132
CENTRIFUGAL PUMP – SPECIFICATION SHEET
E-133
APPENDIX F
F-1
Appendix F Process safety
F-2
F.2. Dow Fire and Explosion Index (F & EI)
1. The procedure for calculating unit hazard factor F & EI shown here is referred from
Figure F.2. [Guide]
Figure F.2 The procedure for calculating unit hazard factor F & EI
2. To determine general process hazards and special process hazards, some figures
from Guide are used in calculation.
For relief pressure item, Figure A10.3 determines the relation of pressure penalty and
set pressure.
F-3
Figure F.3 Pressure penalty for flammable and combustible liquids
For quantity of flammable and unstable materials item, flammability vs. penalty is
given in Figure F.4
3. In order to check if shell and tube reactor is the Pertinent Process Unit, other three
distillation columns’ F & EI are also calculated. The results list as below:
F-4
Table F.2 Dow Fire and Explosion Index Form of light ends distillation column
PREPARED BY CPD 3297 APPROVED BY DATE 03-Dec-03
SITE Light ends distillation column LOCATION Grey area
MATERIALS IN PROCESS UNIT
Propane, nitrogen, ethylene, propylene, methane, hydrogen, water, carbon monoxide
STATE OF OPERATION BASIC MATERIAL(S) FOR MATERIAL FACTOR
Ethylene
_X_DESIGN NORMAL OP __SHUTDOWN
MATERIAL FACTOR (Table I or Appendices A or B) Note requirements when unit temp over 60oC 24
F-5
Table F.3 Dow Fire and Explosion Index Form of ethylene distillation column
PREPARED BY CPD 3297 APPROVED BY DATE 03-Dec-03
SITE Ethylene distillation column LOCATION Grey area
MATERIALS IN PROCESS UNIT
Propane, ethylene, propylene, water
STATE OF OPERATION BASIC MATERIAL(S) FOR MATERIAL FACTOR
Ethylene
_X_DESIGN NORMAL OP __SHUTDOWN
MATERIAL FACTOR (Table I or Appendices A or B) Note requirements when unit temp over 60oC 24
Penalty Factor Penalty Factor
1. GENERAL PROCESS HAZARDS
Range Used1
BASE FACTOR…………………………………………………………………………………………. 1.00 1.00
A. EXOTHERMIC CHEMICAL REACTIONS 0.30 to 1.25 0.00
B. ENDOTHERMIC PROCESSES 0.20 to 0.40 0.00
C. MATERIAL HANDLING AND TRANSFER 0.25 to 1.05 0.50
D. ENCLOSED OR INDOOR PROCESS UNITS 0.25 to 0.90 0.00
E. ACCESS 0.20 to 0.35 0.00
F. DRAINAGE AND SPILL CONTROL 0.00
GENERAL PROCESS HAZARDS FACTOR (F1)………………………………………………… 1.50
2. SPECIAL PROCESS HAZARDS
BASE FACTOR…………………………………………………………………………………………. 1.00 1.00
A. TOXIC MATERIAL(S) 0.20 to 0.80 0.40
B. SUB ATMOSPHERIC PRESSURE (>500 mm Hg) 0.50 0.00
C. OPERATION IN OR NEAR FLAMMABLE RANGE _Inerted _X_ Not Inerted
1. TANK FARMS STORAGE FLAMMABLE LIQUIDS 0.50
2. PROCESS UPSET OR PURGE FAILURE 0.30
3. ALWAYS IN FLAMMABLE RANGE 0.80 0.80
D. DUST EXPLOSION (See Table 3) 0.25 to 2.00 0.00
E. PRESSURE (See Fig 2) Operating Pressure 217.6 psig 0.45
F. LOW TEMPERATURE 0.20 to 0.30 0.00
G. QUALITY OF FLAMMABLE / UNSTABLE MATERIAL : Quantity 17747 lb / kg
HC = 21.5*103 BTU/lb
1. LIQUIDS OR GASES IN PROCESS (See Fig 3) 0.80
2. LIQUIDS OR GASES IN STORAGE (See Fig 4)
3. COMBUSTABLE SOLIDS IN STORAGE, DUST IN PROCESS (See Fig 5)
H. CORROSION AND EROSION 0.10 to 0.75 0.00
I. LEAKAGE - JOINTS AND PACKING 0.10 to 1.50 0.00
J. USE OF FIRED EQUIPMENT (See Fig 6) 0.00
K. HOT OIL EXCHANGE SYSTEMS (See Table 5) 0.15 to 1.15 0.00
L. ROTATING EQUIPMENT 0.50 0.00
F-6
Table F.4 Dow Fire and Explosion Index Form of propylene distillation column
PREPARED BY CPD 3297 APPROVED BY DATE 03-Dec-03
SITE Propylene distillation column LOCATION Grey area
MATERIALS IN PROCESS UNIT
Propane, propylene, water
STATE OF OPERATION BASIC MATERIAL(S) FOR MATERIAL FACTOR
Proylene
_X_DESIGN NORMAL OP __SHUTDOWN
MATERIAL FACTOR (Table I or Appendices A or B) Note requirements when unit temp over 60oC 21
Penalty Factor Penalty
1. GENERAL PROCESS HAZARDS
Range Factor Used1
BASE FACTOR………………………………………………………………………………………… 1.00 1.00
A. EXOTHERMIC CHEMICAL REACTIONS 0.30 to 1.25 0.00
B. ENDOTHERMIC PROCESSES 0.20 to 0.40 0.00
C. MATERIAL HANDLING AND TRANSFER 0.25 to 1.05 0.50
D. ENCLOSED OR INDOOR PROCESS UNITS 0.25 to 0.90 0.00
E. ACCESS 0.20 to 0.35 0.00
F. DRAINAGE AND SPILL CONTROL 0.00
GENERAL PROCESS HAZARDS FACTOR (F1)………………………………………………… 1.50
2. SPECIAL PROCESS HAZARDS
BASE FACTOR…………………………………………………………………………………………. 1.00 1.00
A. TOXIC MATERIAL(S) 0.20 to 0.80 0.40
B. SUB ATMOSPHERIC PRESSURE (>500 mm Hg) 0.50 0.00
C. OPERATION IN OR NEAR FLAMMABLE RANGE __Inerted_X_NotInerted
1. TANK FARMS STORAGE FLAMMABLE LIQUIDS 0.50
2. PROCESS UPSET OR PURGE FAILURE 0.30
3. ALWAYS IN FLAMMABLE RANGE 0.80 0.80
D. DUST EXPLOSION (See Table 3) 0.25 to 2.00 0.00
E. PRESSURE (See Fig 2) Operating Pressure 217.6psig 0.45
F. LOW TEMPERATURE 0.20 to 0.30 0.00
G. QUALITY OF FLAMMABLE / UNSTABLE MATERIAL : Quantity 19257 lb / kg
HC = 21.5*103 BTU/lb
1. LIQUIDS OR GASES IN PROCESS (See Fig 3) 0.84
2. LIQUIDS OR GASES IN STORAGE (See Fig 4)
3. COMBUSTABLE SOLIDS IN STORAGE, DUST IN PROCESS (See Fig 5)
H. CORROSION AND EROSION 0.10 to 0.75 0.00
I. LEAKAGE - JOINTS AND PACKING 0.10 to 1.50 0.00
J. USE OF FIRED EQUIPMENT (See Fig 6) 0.00
K. HOT OIL EXCHANGE SYSTEMS (See Table 5) 0.15 to 1.15 0.00
L. ROTATING EQUIPMENT 0.50 0.00
SPECIAL PROCESS HAZARDS (F2)…………………………………………………………………………………………… 3.49
F-7
APPENDIX G
G-1
Appendix G Economics
Purchased Equipment Costs (PCE) with Lang method, can be found in the Table G.1.
Table G.1 Purchased Equipment Costs
Purchased Equipment Costs (PCE) in December 2003
Type Name Equipment Cost
(US$)
Reactor RX001 Shell & Tube 113,643
RX002 Shell & Tube 113,643
Total 227,286
Drum D101 Vessel 12,354
D102 Vessel 23,488
D201 Vessel 26,382
D301 Vessel 24,217
D302 Vessel 6,298
D401 Vessel 5,212
Total 97,951
Column T301 Packed Column 115,198
T302 Packed Column 128,018
T303 Packed Column 3,669,718
T401 Packed Column 38,049
T402 Packed Column 38,049
Total 3,989,032
Heat exchanger E202 Shell & Tube 117,196
E203 Shell & Tube 74,385
E101 Shell & Tube 96,306
E102 Shell & Tube 26,667
E301A Shell & Tube 24,795
E301B Shell & Tube 164,194
E302A Shell & Tube 113,555
E302B Shell & Tube 82,650
E303A Shell & Tube 166,315
E303B Shell & Tube 77,971
E303C Shell & Tube 22,456
E402 Shell & Tube 61,987
E401 Shell & Tube 148,149
Total 1,028,474
Furnace E201 Process, Cylindrical 1,118,956
Compressor C101 Reciprocating 935,706
C201 Reciprocating 998,811
C303 Reciprocating 209,010
Total 2,143,526
Total Purchased cost (US $) 8,605,226
From the Table 6.1 [Coulson& Richardson, Volumn6], we can calculate Direct Capital Cost and
Fixed Capital Cost in the Table G.2.
G-2
Table G.2 Capital costs estimation
Process type
Item
Fluids
1. Major equipment as total purchased cost
f1 : Equipment erection 0.40
f2 : Piping 0.70
f3 : Instrumentation 0.20
f4 : Electrical 0.10
f5 : Buildings,process 0.15
f6 : Utilities 0.20
f7 : Storages 0.15
f8 : Site development 0.05
f9 : Ancillary buildings 0.15
∑ (f1+f2+f3+…+f9) 2.10
2. Total physical plant cost (PPC)
PPC = PCE(1+ f1 + f2 + …+ f9)
PPC = Direct Cost 26,676,201
3. Indirect cost
f10 : Design and engineering 0.20
f11 : Contractor's fee 0.05
f12 : Contingency 0.05
∑ (f10+f11+f12) 0.30
Indirect Capital Cost 8,002,860
Fixed Capital = PPC (1 + f10 + f11 +f12)
Fixed Capital = Direct + Indirect
cost 34,679,062
2003
Costs
US$
1. Direct Capital Cost 26,676,201
2. Indirect Capital Cost 8,002,860
3. Fixed Capital Cost 34,679,062
Raw material and utilities costs, which is used in the process, determined in Table G.3
G-3
Table G.3 Raw material and utilities costs
Stream hrs/annum = 8040
Cost
t/t Cost US$/a Million
Raw Materials Str.No. kg/s m3/a ton/hrs ton/a Price US$/unit Unit
Alkenes @2003 US$/a
@2003
Propane 001 8.23 - 29.64 238,341 1.1748 160 Ton 38,194,083 38.194
Oxygen 008 1.87 - 6.74 54,194 0.2671 143 Ton 7,771,432 7.771
Catalyst1
(V2O5/CeO2/SA5205) 6.8 0.00003 106,549 Ton 724,535 0.725
Catalyst2 (Pt on MFI zeolite) 3.3 0.00002 181,939 Ton 600,399 0.600
Total catalyst cost Ton 1,324,934 1.325
Total Raw Material cost (IN) 47,290,450 47.290
Income
t/t Income US$/a Million
Product Str.No. kg/s m3/a ton/hrs ton/a Price US$/ton Unit
Alkenes @2003 US$/a
@2003
Ethylene 308 2.29 - 8.26 66,376 0.3272 518 ton 34,382,700 34.383
Propylene 313 4.72 16.98 136,502 0.6728 408 ton 55,692,827 55.693
Light gas for syngas plant 303 1.47 5.28 42,414 0.2091 64 ton 2,718,717 2.719
CO2 for EOR 402 0.16 0.58 4,670 0.0230 6.5 ton 30,356 0.030
Water for EOR 208-5 1.46 5.26 42,267 0.2083 0.01 ton 423 0.000
Total Income (OUT) 92,825,022 92.825
Water (t/a)
11,188,806 0.01 Ton 111,888 0.112
Electric (kWh/a)
342,596,748 0.04 kWh 14,469,791 14.470
Total Utilities cost 14,581,679 14.582
G-4
Economic Criteria
a) Net Cash Flow (NCF) can be calculated from Gross Income and Production
Costs.
Net Cash Flow annual = ∑ (Gross Income annual - Production Costs annual )
⎛ Price ⎞
Gross Income annual = Σ ⎜ Products annual x ⎟
⎝ Unit ⎠
= 92,825,022 US$/a
= 92.83 US$ million
Production Cost annual = 79.10 US$ million
Net Cash Flow annual = 92.83- 79.10 US$ million
= 13.72 US$ m illion
b) Rate of Return (ROR) and Pay Out (Back) Time (POT or PBP) from the total
Investment and NCF along ref(1) approach.
G-5
In calculating cash flow (NCF), the project is usually considered as an isolated system,
and taxes on profits and the effect of depreciation of the investment are not considered,
since tax rates are not constant and depend on government policy as well as the rates
of depreciation. Depreciation rates also depend on the accounting practices of the
particular company. Therefore during evaluating projects, the effect of government
policy must be taken into account at some stages particularly when considering
projects in different country.
200
150
100
Cash Flow (mil $US)
NFV (i=0%)
50
NPV (i=8%)
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
-50
-100
Year
G-6
APPENDIX H
H-1
Appendix H.1 PFS for process with heat integration
H-2
Figure H.1.1 Process flow scheme Part 1.
H-3
Figure H.1.2 Process flow scheme Part 2
H-4
Appendix H.2 Summary of utilities for the process with
Heat integration
H-5
Table H.2. Summary of utilities for process with Heat integration
SUMMARY OF UTILITIES
EQUIPMENT UTILITIES
H-6
Table H.2. Summary of utilities for process with Heat integration (con’t)
SUMMARY OF UTILITIES
EQUIPMENT UTILITIES
H-7
Table H.2. Summary of utilities for process with Heat integration (con’t)
SUMMARY OF UTILITIES
EQUIPMENT UTILITIES
H-8
Appendix H.3 Heat and Mass balance for the process with Heat integration
Table H.3.1a Mass and heat balance for the process with heat integration
STREAM Nr. : 001 IN 002 003 004 005=004+318
Name : Propane Feed In Preheated Propane1 Preheated Propane2 Preheated Propane3 Total Propane Feed
COMP MW kg/s kmol/s kg/s kmol/s kg/s kmol/s kg/s kmol/s kg/s kmol/s
Propane 44 7.45 0.1692 7.45 0.1692 7.45 0.1692 7.45 0.1692 16.36 0.3717
Oxygen 32 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Ethylene 28 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Propylene 42 0.79 0.0188 0.79 0.0188 0.79 0.0188 0.79 0.0188 1.71 0.0407
Carbonmon-oxide 28 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Carbondi-oxide 44 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Methane 16 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Hydrogen 2 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Water 18 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Nitrogen 28 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Total 8.23 0.1880 8.23 0.1880 8.23 0.1880 8.23 0.1880 18.07 0.4124
Enthalpy kW -17379 -12289 -6334 -5975 -13167
Phase V V V V V
Press. Bara 1 1 1 1 1
Temp oC 25 298 527.5 540 540
H-9
Table H.3.1a Mass and heat balance for the process with heat integration (Con’t)
STREAM Nr. : 006 007 008 009 010
Name : Propane to Shell Rx Propane splitted to Tube Oxygen feed Material to Tube Feed to Tube
COMP MW kg/s kmol/s kg/s kmol/s kg/s kmol/s kg/s kmol/s kg/s kmol/s
Propane 44 9.81 0.2230 6.54 0.1487 0.00 0.0000 6.54 0.1487 6.54 0.1487
Oxygen 32 0.00 0.0000 0.00 0.0000 1.79 0.0559 1.79 0.0559 1.79 0.0559
Ethylene 28 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Propylene 42 1.03 0.0244 0.68 0.0163 0.00 0.0000 0.68 0.0163 0.68 0.0163
Carbonmon-oxide 28 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Carbondi-oxide 44 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Methane 16 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Hydrogen 2 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Water 18 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Nitrogen 28 0.00 0.0000 0.00 0.0000 0.08 0.0029 0.08 0.0029 0.08 0.0029
Total 10.84 0.2475 7.23 0.1650 1.87 0.0589 9.10 0.2239 9.10 0.2239
Enthalpy kW -7900 -5267 959 -4308 4957
Phase V V V V V
Press. Bara 1 1 1 1 1
Temp oC 540 540 25 795 850
H-10
Table H.3.1a Mass and heat balance for the process with heat integration (Con’t)
STREAM Nr. : 101 102 103 104 105
Name : Shell gas product1 Shell gas product2 Shell gas product3 Shell gas product4 Shell gas compressed
COMP MW kg/s kmol/s kg/s kmol/s kg/s kmol/s kg/s kmol/s kg/s kmol/s
Propane 44 7.56 0.1717 7.56 0.1717 7.56 0.1717 7.56 0.1717 7.56 0.1717
Oxygen 32 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Ethylene 28 0.003 0.0001 0.003 0.0001 0.003 0.0001 0.003 0.0001 0.003 0.0001
Propylene 42 3.18 0.0756 3.18 0.0756 3.18 0.0756 3.18 0.0756 3.18 0.0756
Carbonmon-oxide 28 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Carbondi-oxide 44 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Methane 16 0.002 0.0001 0.002 0.0001 0.002 0.0001 0.002 0.0001 0.002 0.0001
Hydrogen 2 0.10 0.0515 0.10 0.0515 0.10 0.0515 0.10 0.0515 0.10 0.0515
Water 18 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Nitrogen 28 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Total 10.84 0.2990 10.84 0.2990 10.84 0.2990 10.84 0.2990 10.84 0.2990
Enthalpy kW -1247 -9325 -16392 -16528 -12344
Phase V V V V V
Press. Bara 1 1 1 1 30
Temp oC 540 308 32 25 227
H-11
Table H.3.1a Mass and heat balance for the process with heat integration (Con’t)
STREAM Nr. : 106 107 108 109 201
Name : Shell gas compressed2 Shell gas comp.after cooler Shell gas to sep.Unit Shell liq to sep. Unit Tube gas product1
COMP MW kg/s kmol/s kg/s kmol/s kg/s kmol/s kg/s kmol/s kg/s kmol/s
Propane 44 7.56 0.1717 7.56 0.1717 5.08 0.1155 2.48 0.0563 1.36 0.0309
Oxygen 32 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Ethylene 28 0.003 0.0001 0.003 0.0001 0.002 0.0001 0.000 0.0000 2.291 0.0818
Propylene 42 3.18 0.0756 3.18 0.0722 2.22 0.0528 0.96 0.0228 2.47 0.0588
Carbonmon-oxide 28 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.62 0.0223
Carbondi-oxide 44 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.18 0.0040
Methane 16 0.002 0.0001 0.002 0.0000 0.002 0.0001 0.000 0.0000 0.576 0.0360
Hydrogen 2 0.10 0.0515 0.10 0.0023 0.10 0.0504 0.00 0.0011 0.05 0.0265
Water 18 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 1.47 0.0817
Nitrogen 28 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.08 0.0029
Total 10.84 0.2990 10.84 0.2464 7.40 0.2189 3.43 0.0801 9.10 0.3449
Enthalpy kW -16373 -17475 -11090 -6385 -4243
Phase V L/V V L V
Press. Bara 30 30 30 30 1
Temp oC 65 55 55 55 850
H-12
Table H.3.1a Mass and heat balance for the process with heat integration (Con’t)
STREAM Nr. : 202 203 204 205 206
Name : Tube gas product2 Tube gas product3 Tube gas product4 Tube gas product5 Tube gas compressed
COMP MW kg/s kmol/s kg/s kmol/s kg/s kmol/s kg/s kmol/s kg/s kmol/s
Propane 44 1.36 0.0309 1.36 0.0309 1.36 0.0309 1.36 0.0309 1.36 0.0309
Oxygen 32 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Ethylene 28 2.291 0.0818 2.291 0.0818 0.082 0.0029 2.291 0.0818 2.291 0.0521
Propylene 42 2.47 0.0588 2.47 0.0588 2.29 0.0545 2.47 0.0588 2.47 0.0561
Carbonmon-oxide 28 0.62 0.0223 0.62 0.0223 2.47 0.0882 0.62 0.0223 0.62 0.0142
Carbondi-oxide 44 0.18 0.0040 0.18 0.0040 0.58 0.0131 0.18 0.0040 0.18 0.0040
Methane 16 0.576 0.0360 0.576 0.0360 0.053 0.0033 0.576 0.0360 0.576 0.0131
Hydrogen 2 0.05 0.0265 0.05 0.0265 0.12 0.0576 0.05 0.0265 0.05 0.0012
Water 18 1.02 0.0566 0.57 0.0315 0.62 0.0346 0.12 0.0064 0.12 0.0026
Nitrogen 28 0.08 0.0029 0.08 0.0029 0.18 0.0063 0.08 0.0029 0.08 0.0019
Total 8.65 0.3198 8.19 0.2946 7.74 0.2915 7.74 0.2695 7.74 0.1760
Enthalpy kW -12982 -18961 -3102 -6063 -6063
Phase V V V V V
Press. Bara 1 1 1 1 1
Temp oC 550 308 219 25 25
H-13
Table H.3.1a Mass and heat balance for the process with heat integration (Con’t)
STREAM Nr. : 207 208(208-5) =208-1 to -4 209 301 302
Name : Tube gas compressed2 Water discharged Tube gas to CO2 removal Propane to Tube Overhead T301
COMP MW kg/s kmol/s kg/s kmol/s kg/s kmol/s kg/s kmol/s kg/s kmol/s
Propane 44 1.36 0.0309 0.00 0.0000 1.36 0.0309 6.44 0.1464 0.00 0.0000
Oxygen 32 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Ethylene 28 0.082 0.0029 0.00 0.0000 2.29 0.0818 2.29 0.0819 7.90 0.2821
Propylene 42 2.29 0.0545 0.00 0.0000 2.47 0.0588 4.69 0.1116 0.00 0.0000
Carbonmon-oxide 28 2.47 0.0882 0.00 0.0000 0.62 0.0223 0.62 0.0223 1.19 0.0426
Carbondi-oxide 44 0.58 0.0131 0.00 0.0000 0.18 0.0040 0.01 0.0003 0.01 0.0002
Methane 16 0.053 0.0033 0.00 0.0000 0.58 0.0360 0.58 0.0361 4.79 0.2994
Hydrogen 2 0.12 0.0576 0.00 0.0000 0.05 0.0265 0.15 0.0769 0.16 0.0797
Water 18 0.62 0.0346 1.46 0.0332 0.01 0.0006 0.01 0.0006 0.00 0.0000
Nitrogen 28 0.18 0.0063 0.00 0.0000 0.08 0.0029 0.08 0.0029 0.14 0.0050
Total 7.74 0.2915 1.46 0.0332 7.64 0.2637 14.88 0.4790 14.19 0.7089
Enthalpy kW -1921 -20359 -4946 -14585 -5787
Phase V L V V V
Press. Bara 30 3 30 30 15
Temp oC 289 25 30 42 -131
H-14
Table H.3.1a Mass and heat balance for the process with heat integration (Con’t)
STREAM Nr. : 303 304 305 306 307
Name : Light gas Product Reflux T301 Feed to T302 Overhead T302 Reflux T302
COMP MW kg/s kmol/s kg/s kmol/s kg/s kmol/s kg/s kmol/s kg/s kmol/s
Propane 44 0.00 0.0000 0.00 0.0000 6.44 0.1464 0.01 0.0002 0.01 0.0002
Oxygen 32 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Ethylene 28 0.03 0.0010 7.87 0.2811 2.26 0.0809 2.26 0.0808 0.00 0.0000
Propylene 42 0.00 0.0000 0.00 0.0000 4.69 0.1116 8.53 0.2032 8.52 0.2029
Carbonmon-oxide 28 0.62 0.0223 0.57 0.0203 0.00 0.0000 0.05 0.0018 0.05 0.0018
Carbondi-oxide 44 0.00 0.0000 0.01 0.0002 0.01 0.0003 0.01 0.0003 0.00 0.0000
Methane 16 0.58 0.0361 4.21 0.2633 0.00 0.0000 0.00 0.0000 0.00 0.0000
Hydrogen 2 0.15 0.0769 0.01 0.0027 0.00 0.0000 0.00 0.0000 0.00 0.0000
Water 18 0.00 0.0000 0.00 0.0000 0.01 0.0006 0.00 0.0000 0.00 0.0000
Nitrogen 28 0.08 0.0029 0.06 0.0021 0.00 0.0000 0.05 0.0019 0.05 0.0019
Total 1.47 0.1392 12.72 0.5697 13.42 0.3398 10.93 0.2883 8.63 0.2068
Enthalpy kW -5786 -5787 -14573 -14264 11270
Phase V V L V L
Press. Bara 15 15 15 15 15
Temp oC -131 -131 10 -39 -39
H-15
Table H.3.1a Mass and heat balance for the process with heat integration (Con’t)
STREAM Nr. : 308 309 310 311 312
Name : Ethylene Prduct Bottom product T303 Feed to T303 Overhead T303 Reflux T302
COMP MW kg/s kmol/s kg/s kmol/s kg/s kmol/s kg/s kmol/s kg/s kmol/s
Propane 44 0.00 0.0000 6.44 0.1463 8.91 0.2026 0.05 0.0011 0.05 0.0010927
Oxygen 32 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0
Ethylene 28 2.26 0.0808 0.00 0.0001 0.00 0.0001 0.04 0.0015 0.04 0.0014474
Propylene 42 0.01 0.0003 4.67 0.1113 5.63 0.1340 105.15 2.5036 100.44 2.391449
Carbonmon-oxide 28 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0
Carbondi-oxide 44 0.01 0.0003 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 7.027E-06
Methane 16 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0002 0.00 0.000154
Hydrogen 2 0.00 0.0000 0.00 0.0000 0.00 0.0011 0.05 0.0242 0.05 0.0230689
Water 18 0.00 0.0000 0.01 0.0006 0.01 0.0006 0.00 0.0000 0.00 6.849E-39
Nitrogen 28 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0
Total 2.29 0.0815 11.12 0.2582 14.56 0.3384 105.29 2.5306 100.58 2.4172
Enthalpy kW 2994 -16920 -23305 0 0
Phase L L L V L
Press. Bara 15 15 15 15 15
Temp oC -39 40 38 35 35
H-16
Table H.3.1a Mass and heat balance for the process with heat integration (Con’t)
STREAM Nr. : 313 314 314-W 315 316
Name : Propylene Prduct Propane recycle1 Moisture removed Propane recycle2 Propane recycle3
COMP MW kg/s kmol/s kg/s kmol/s kg/s kmol/s kg/s kmol/s kg/s kmol/s
Propane 44 0.00 0.0001 8.91 0.2025 0.00 0.0000 8.91 0.2025 8.91 0.2025
Oxygen 32 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Ethylene 28 0.00 0.0001 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Propylene 42 4.71 0.1121 0.92 0.0219 0.00 0.0000 0.92 0.0219 0.92 0.0219
Carbonmon-oxide 28 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Carbondi-oxide 44 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Methane 16 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Hydrogen 2 0.00 0.0011 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Water 18 0.00 0.0000 0.00 0.0000 0.01 0.0006 0.00 0.0000 0.00 0.0000
Nitrogen 28 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Total 4.72 0.1133 9.83 0.2244 0.01 0.0006 9.83 0.2244 9.83 0.2244
Enthalpy kW 2093 -24058 -168 -17002 -14849
Phase V L L V V
Press. Bara 15 15 15 15 15
Temp oC 35 43 43 217 298
H-17
Table H.3.1a Mass and heat balance for the process with heat integration (Con’t)
STREAM Nr. : 317 318 401 402
Name : Propane recycle4 Propane recycle5 CO2 less gas product CO2 By product
COMP MW kg/s kmol/s kg/s kmol/s kg/s kmol/s kg/s kmol/s
Propane 44 8.91 0.2025 8.91 0.2025 1.36 0.0309 0.00 0.0000
Oxygen 32 0.00 0.0000 0.00 0.0000 0.00 0.0000 0.00 0.0000
Ethylene 28 0.00 0.0000 0.00 0.0000 2.29 0.0818 0.00 0.0000
Propylene 42 0.92 0.0219 0.92 0.0219 2.47 0.0588 0.00 0.0000
Carbonmon-oxide 28 0.00 0.0000 0.00 0.0000 0.62 0.0223 0.00 0.0000
Carbondi-oxide 44 0.00 0.0000 0.00 0.0000 0.01 0.0003 0.16 0.0037
Methane 16 0.00 0.0000 0.00 0.0000 0.58 0.0360 0.00 0.0000
Hydrogen 2 0.00 0.0000 0.00 0.0000 0.05 0.0265 0.00 0.0000
Water 18 0.00 0.0000 0.00 0.0000 0.01 0.0006 0.00 0.0000
Nitrogen 28 0.00 0.0000 0.00 0.0000 0.08 0.0029 0.00 0.0000
Total 9.83 0.2244 9.83 0.2244 7.48 0.2601 0.16 0.0037
Enthalpy kW -9407 -7193 -3494 -1451
Phase V V V V
Press. Bara 15 1 30 1
Temp oC 476 540 30 101
H-18
Table H.3.1a Mass and heat balance for the process with heat integration (Con’t)
STREAM Nr. : 403 404
Name : MDEA sol feed Spent MDEA sol
COMP MW kg/s kmol/s kg/s kmol/s
Propane 44 0.00 0.0000 0.00 0.0000
Oxygen 32 0.00 0.0000 0.00 0.0000
Ethylene 28 0.00 0.0000 0.00 0.0000
Propylene 42 0.00 0.0000 0.00 0.0000
Carbonmon-oxide 28 0.00 0.0000 0.00 0.0000
Carbondi-oxide 44 0.00 0.0000 0.16 0.0037
Methane 16 0.00 0.0000 0.00 0.0000
Hydrogen 2 0.00 0.0000 0.00 0.0000
Water 18 15.96 0.8864 15.96 0.8864
Nitrogen 28 0.00 0.0000 0.00 0.0000
MDEA 119.2 15.96 0.1339 15.96 0.1339
Total 31.91 0.8864 32.07 0.8901
Press. Bara 30 30
Temp oC 30 30
Table H.3.1bOverall mass and heat balance for the process with heat integration (Con’t)
Overall Component
Mass Balance &
Stream Heat balance
STREAM Nr. : 001+008 303+308+313+402+208+314W OUT-IN
Name : Total Plant IN Total Plant OUT Total Plant
COMP MW kg/s kmol/s kg/s kmol/s kg/s kmol/s
Propane 44 7.45 0.1692 0.00 0.0001 -7.44 -0.1691
Oxygen 32 1.79 0.0559 0.00 0.0000 -1.79 -0.0559
Ethylene 28 0.00 0.0000 2.29 0.0819 2.29 0.0819
Propylene 42 0.79 0.0188 4.72 0.1125 3.93 0.0937
Carbonmon-oxide 28 0.00 0.0000 0.62 0.0223 0.62 0.0223
Carbondi-oxide 44 0.00 0.0000 0.18 0.0040 0.18 0.0040
Methane 16 0.00 0.0000 0.58 0.0361 0.58 0.0361
Hydrogen 2 0.00 0.0000 0.16 0.0780 0.16 0.0780
Water 18 0.00 0.0000 1.47 0.0817 1.47 0.0817
Nitrogen 28 0.08 0.0029 0.08 0.0029 0.00 0.0000
Total 10.11 0.2439 10.11 0.4165 0.00
H-19
Appendix H.4 Economy
Profitability analysis is required to evaluate the economic aspect of a project design.
Some crucial measures such as Purchased Equipment (PCE), Total Investment cost,
production cost, profit margin, return on investment (ROI), payback period (PBP). The
cash flow, net present value (NPV) and the investor’s rate of return (IRR) (also know as
the discounted cash-flow rate of return (DCFRR) to count the time value of money.
Margin is defined as
follows:
Margin=Total Value (Products, Wastes OUT ) - Total Value (Feedstock's, Process Chemicals, IN )
In order to prevent confusion, in this conceptual stage, profitability analysis would not
include some local regulation such as, local taxes, depreciation, subsidy, grant etc.
Thereby, they will be taken on board in further stage.
In this chapter, it presents the economic indexes such as investment cost, margin, and
economic criteria in order to view the feasibility of proposed design.
Investment
Lang method is used in order to estimate the investment of the designed process.
Investment is considered in term of direct capital cost, indirect capital cost, fixed capital
cost, license cost, and working capital cost.
Some commodity chemicals can be obtained from the Chemical Marketing Reporter and
Chemical week magazine.
Direct Capital Cost: The summation of the material costs required to build the complete
Process that are incurred in the construction of a plan, in addition to the cost of
equipment are
1. Equipment erection, including foundations and minor structural work.
2. Piping, including insulation and painting
3. Electrical, power and lighting.
4. Instruments, local and control room
5. Process building and structures
6. Ancillary buildings, offices, laboratory buildings, workshops.
7. Storages, raw materials and finished products.
8. Utilities, provision of plant for steam, water, air, firefighting services
9. Site and site preparation
Indirect Capital Cost: The costs that are first, for Design and engineering, which cover
the cost of design and the cost of engineering, purchasing, procurement and construction
supervision. Secondly, for Contractor’s fees and for Contingency allowance.
Fixed Capital Cost is the summation of Direction Capital Cost and Indirection Capital
Cost. Purchased Equipment Costs (PSE) with Lang method and capital cost can be found
in the Table H.4.1. The summary of the capital investment cost presents in Table H.4.2.
H-20
Purchased Equipment Costs (PCE) with Lang method, can be found in the Table H.4.1.
Table H.4.1 Purchased Equipment Costs
Purchased Equipment Costs (PCE) in December 2003
Type Name Equipment Cost
(US$)
Reactor RX001 Shell & Tube 113,643
RX002 Shell & Tube 113,643
Total 227,286
Drum D101 Vessel 12,354
D102 Vessel 23,488
D201 Vessel 26,382
D301 Vessel 24,217
D302 Vessel 6,298
D401 Vessel 5,212
Total 97,951
Column T301 Packed Column 115,198
T302 Packed Column 128,018
T303 Packed Column 3,669,718
T401 Packed Column 38,049
T402 Packed Column 38,049
Total 3,989,032
Heat exchanger E-001 Shell & Tube 3,868,958
E-002 Shell & Tube 2,446,242
E-004 Shell & Tube 1,173,004
E-101 Shell & Tube 35,126
E-102 Shell & Tube 643,563
E-103 Shell & Tube 30,994
E-201 Shell & Tube 636,338
E-202 Shell & Tube 1,961,731
E-203 Shell & Tube 65,577
E-204 Shell & Tube 61,987
E-301A Shell & Tube 24,795
E-301B Shell & Tube 164,194
E-302A Shell & Tube 113,555
E-302B Shell & Tube 82,650
E-303A Shell & Tube 166,315
E-303B Shell & Tube 77,971
E-303C Shell & Tube 22,456
E-402 Shell & Tube 61,987
E-401 Shell & Tube 148,149
Total 11,785,592
Furnace E-003 Process, Cylindrical 34,504
E-005 Process, Cylindrical 111,404
E-006 Process, Cylindrical 117,354
Total 263,262
Compressor C101 Reciprocating 935,706
C201 Reciprocating 998,811
C303 Reciprocating 209,010
Total 2,143,526
H-21
From the table 6.1 (Reference: Coulson& Richardson, Volumn6), we can calculate Direct
Capital Cost and Fixed Capital Cost in the Table H.4.2.
Table H.4.2 Capital costs estimation
Process type
Item
Fluids
1. Major equipment as total purchased cost
f1 : Equipment erection 0.40
f2 : Piping 0.70
f3 : Instrumentation 0.20
f4 : Electrical 0.10
f5 : Buildings,process 0.15
f6 : Utilities 0.20
f7 : Storages 0.15
f8 : Site development 0.05
f9 : Ancillary buildings 0.15
∑ (f1+f2+f3+…+f9) 2.10
2. Total physical plant cost (PPC)
PPC = PCE(1+ f1 + f2 + …+ f9)
PPC = Direct Cost 57,370,614
3. Indirect cost
f10 : Design and engineering 0.20
f11 : Contractor's fee 0.05
f12 : Contingency 0.05
∑ (f10+f11+f12) 0.30
Indirect Capital Cost 17,211,184
Fixed Capital = PPC (1 + f10 + f11 +f12)
Fixed Capital = Direct + Indirect cost 74,581,798
2003
Costs
US$
1. Direct Capital Cost 57,370,614
2. Indirect Capital Cost 17,211,184
3. Fixed Capital Cost 74,581,798
To estimate the total investment cost of this process, the total investment is fixed capital cost and working capital.
Therefore the total capital cost is summarized in Table H.4.3.
H-22
The raw material and utilities are presented in Table H.4.4. Meanwhile, Production cost is shown in Table H.4.5.
Table H.4.4 The raw material and utilities cost
Stream hrs/annum = 8040
Cost
t/t Price Cost US$/a Million
Raw Materials Str.No. kg/s m3/a ton/hrs ton/a Unit
Alkenes US$/unit @2003 US$/a
@2003
Propane 1 8.23 - 29.64 238,341 1.1748 160 Ton 38,194,083 38.194
Oxygen 8 1.87 - 6.74 54,194 0.2671 143 Ton 7,771,432 7.771
Catalyst1 (V2O5/CeO2/SA5205 6.8 0.00003 106,549 Ton 724,535 0.725
Catalyst2 (Pt on MFI zeolite) 3.3 0.00002 181,939 Ton 600,399 0.6
Total catalyst cost Ton 1,324,934 1.325
Total Raw Material cost (IN) 47,290,450 47.29
Income
Income
t/t Price Million
Product Str.No. kg/s m3/a ton/hrs ton/a Unit US$/a
Alkenes US$/ton US$/a
@2003
@2003
Ethylene 308 2.29 - 8.26 66,376 0.3272 518 ton 34,382,700 34.383
Propylene 313 4.72 16.98 136,502 0.6728 408 ton 55,692,827 55.693
Light gas for syngas plant 303 1.47 5.28 42,414 0.2091 64 ton 2,718,717 2.719
CO2 for EOR 402 0.16 0.58 4,670 0.023 6.5 ton 30,356 0.03
Water for EOR 208-5 1.46 5.26 42,267 0.2083 0.01 ton 423 0
Total Income (OUT) 92,825,022 92.825
H-23
Table H.4.4 The raw material and utilities cost (Con’t)
Cost
Cost Cost US$/a Million
Utilities Load Unit
US$/unit @2003 US$
@2003
Water (t/a)
11,281,221 0.010 Ton 112,812 0.113
Steam (t/a)
12.000 Ton 0.000
Electric (kWh/a)
106,097,231 0.04 kWh 4,481,084 4.481
Fuel (t/a) 0
Total Utilities cost 4,593,896 4.594
H-24
Economic criteria
In order to determine the economic, Process Cash Flow, Rate of Return (ROR) and
Pay Out (Back) Time (POT) of the investment are summarized in Table E.4.6.
Table E.4.6 Economic Criteria
Economic criteria
Cash Flow 15.77
M US$
Rate of 11.17
Return
(ROR) [%]
Pay Out 10
Time yrs
DCFROR
% 15.31
Economic Criteria
c) Net Cash Flow (NCF) from Gross Income and Production Costs.
Net Cash Flow annual = ∑ (Gross Income annual - Production Costs annual )
⎛ Price ⎞
Gross Income annual = Σ ⎜ Products annual x ⎟
⎝ Unit ⎠
= 92,825,022 US$/a
= 93.83 US$ million
Production Cost annual = 77.057 US$ million
Net Cash Flow annual = 93.83- 77.05 US$ million
= 15.77 US$ million
d) Rate of Return (ROR) and Pay Out (Back) Time (POT or PBP) from the total
Investment and NCF along approach.
Accu. Cash Flow
ROR =
(Project life*Tot.Investment)
H-25
Project life time = Construction time + Plant life time + Salvage
= 2 + 15 +1
= 18 years
So,
161.94 US$ million
ROR =
18 years * 80.55 US$ million
= 11.17 %
and,
POT = 10 years
200
150
100
Cash Flow (mil $US)
NFV
50
NPV (i=8%)
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
-50
-100
Year
H-26