Chemical and Petroleum Engineering Department
Chemical and Petroleum Engineering Department
Submitted by:
I.D:
3115204
supervised :
Dr. Mohammed Edali
1
Fall (2019-2020)
Problem statement
It is required to design 30 m3 a storage tank to
store mixture ( i-Butane 30% , n-Butane 70%)
under 5.5 bar And 48.5 °C.
Information table:
Description Symbo Valu Uni
l e t
Volume V 30 m3
o
Temperature T 48.5 C
2
Page
Chapters Contents
number
1.1. Introduction 7
32
4
References 52
Page
Chapters List of figures numbe
r
Page
Chapters List of tables
number
6
1.1. Introduction
Chemical engineering involves the application of sciences to the process
industries, which are primarily concerned, with the conversion of one material
into another by chemical or physical means. These processes require the
handling or storing of large quantities of materials in containers of varied
constructions, depending upon the existing state of the material, it's physical and
chemical properties and the required operations, which are to be performed. For
handling such liquids and gases, a container or vessel is used. It is called a
pressure vessel, when they are containers for fluids subjected to pressure. They
are leak proof containers. They may be of any shape ranging from types of
processing equipment. Most process equipment units may be considered as
vessels with various modifications necessary to enable the units to perform
certain required functions, e.g. an autoclave may be considered as high pressure
vessel equipped with agitation and heating sources.
7
concentration which may develop when the end closure of a high-pressure vessel
is attached to a conically shaped nozzle. The main reason for this occurrence is
that the conical nozzle must be connected separately. This process would result
in geometrical discontinuities between the pressure vessel and the nozzle at the
point of attachment . The solution for the value of stress at the connection of
cylindrical nozzle to an ellipsoidal shape pressure vessel. The stress calculations
were carried out using finite element method, and a parametric model was
developed . Tests were carried out for various positions of the cylindrical nozzle
on the pressure vessel along with multiple orientations of the connection instead
of a central radial hole . The tests carried out confirmed that the actual stress
value of the maximum allowable internal pressure is different for various
orientations of the attachment and states that further investigation is necessary
to determine the most efficient connection.
According to the intended use of the pressure vessel, they can be divided into
8
storage containers and process vessels.
The first classes are only used for storing fluids under pressure, and in
accordance with the service are known as storage tanks. Process pressure vessels
have multiple and varied uses, among them we can mention heat exchangers,
reactors, fractionating towers, distillation towers, etc.
9
Figure 1.3: Vertical pressure vessel
11
2.1. Calculation of vessel dimensions:-
4∗V D
D i= 3
√ L
D
∗π
… … … … … … … … … … … … … … … … …(Eq .2 .1)
Where :
V : Volume operating = 30 m3
D: Vessel diameter, m
L: Length of vessel, m
L/D = (2.5 – 5)
12
By substitution :
The vessel will be horizontal (because : volume design, VD > 5000 gal)
4∗33
√
D i= 3
( 2.5 )∗(π)
= 2.5615 m or 2561.5 mm
L
L= ∗D i … … … … … … … … … … … … … … … … … … … … … … … ( Eq.2 .2)
D
By substitution :
P D∗Di
t v= +C … … … … … … … … … … … … … … …(Eq .2 .3)
2∗J∗f −PD
Where :
13
f: Design stress, N/mm², from table 2.1
By substitution :
P D=0.605 N /mm 2.
14 J=0.85
Di= 2561.9 mm
0.605∗2561.5
t v= +2=9.313568 mm
2∗125∗0.85−0.605
D o =D i+ 2∗( t v ∗10−3 ) … … … … … … … … … … … … … … ( Eq .2 .3 )
By substitution :
Di + Do
Dm = … … … … … … … … … … … … … … … … … … … … … . ( Eq: 2.1.5 )
2
By substitution :
15
Where :
By substitution :
W v =¿11213. 9N
L/Di Di L tv Do Dm Wv
- m m mm m m N
2.5 2.5615 6.4038 9.313568762 2.58013 2.570815 11213.9
3 2.4105 7.2314 8.882331705 2.42823 2.419348 10064.3
3.5 2.2897 8.0141 8.53762497 2.30681 2.298273 9189.35
4 2.1901 8.7602 8.253013334 2.20656 2.198306 8496.45
4.5 2.1057 9.4758 8.012270529 2.12176 2.113748 7931.16
5 2.0331 10.165 7.804783373 2.04867 2.04087 7459.26
From the table , the less wight when the ratio L/D = 5
Then:
16
Length of vessel is 10.165327 m
The minimum thickness from calculations is less than the minimum thickness
from table, so that is acceptable
17
Chapter three:
Physical properties
18
3.1. Physical properties of material :
The material that is used is a carbon steel, the physical properties of carbon steel are shown
that in table 3.1.
Density kg/m3 24
19
3.3. Physical properties of air:-
Physical properties of air taken at Tfilm= 36°C (309K) . By interpolation between
temperature (300 K and 350 K). this properties shown in table 3.3.
Prandtl 0.69652
20 Prandtl 3.67856
Chapter four:
Mechanical design
21
4.1. Choice of closure and head thickness
Flat plates are used as covers for manways, and as the channel covers of heat
exchangers. Formed flat ends, known as “flange-only” ends, are manufactured by
turning over a flange with a small radius on a flat plate, Figure 4.1.a. The corner
radius reduces the abrupt
Figure .4.1. Domed heads (a) Hemispherical (b) Ellipsoidal (c) Torispherical
change of shape, at the junction with the cylindrical section; which reduces the
local stresses to some extent: “Flange-only” heads are the cheapest type of formed
head to manufacture, but their use is limited to low-pressure and small-diameter
vessels. Standard torispherical heads (dished ends) are the most commonly used
end closure for vessels up to operating pressures of 15 bar. They can be used for
higher pressures, but above 10 bar their cost should be compared with that of an
equivalent ellipsoidal head. Above 15 bar an ellipsoidal head will usually prove to
be the most economical closure to use. A hemispherical head is the strongest
22
shape; capable of resisting about twice the pressure of a torispherical head of the
same thickness. The cost of forming a hemispherical head will, however, be
higher than that for a shallow torispherical head. Hemispherical heads are used
for high pressures.
4.1.1. For section Head thickness (Torispherical) :
Because the Design pressure is less then 15 bar. The thickness of the upper and
lower head is given :
PD∗Rc∗103∗Cs
t h= +C …………………………………………(Eq.4.1)
2∗f∗J + PD (Cs−0.2)
Where:
1 R
4 √
C s= * (3+ 2 c )………………………………………...(Eq.4.2)
Rk
Rc :Crow radius=D i , m.
Rk : Knuckle radius=0.06∗Rc ,m .
By substitution:
Rc = 2.0330654 m
23 1
C s= * (3+ 2 2.0330654 )
4 √0.121983925
Cs = 1.770621
0.605∗2.0330654∗103∗1.770621
t h= +2
2∗0.85∗125+0.605∗( 1.770621−0.2 )
th = 12.20318 mm
. Low cost.
. Low density.
∆T
q= …………………………………(Eq.4.3)
∑R
Where:
∑ R=Rm + R¿ + Rf + ¿ R a ¿………………………..(Eq.4.4)
24
Where:
T f −T w
q without Insulation=
r2
ln ( ) ……………………………..
1 r1
+
hiAi 2 π∗km∗L
(Eq.4.5)
Where:
Assume:-
53.35−30
boundary layer ¿ = 11.7 °C
2
T w =53.35+11.7=41.7 ° C
T w + T air
T film= ………………………………………………(Eq.4.8)
2
41.7+30
T film ¿ = 36°C or 309 K at properties of air from table.3.3.
2
Nu∗k a
h0 = … … … … … … … … … … … … … … … … … … … … … … … … … … .(Eq .4 .9)
l
1 1
0.62 R e 2 p r 3 pr > 0.2 ………( Eq .4 .10)
Nu=0.3+
26
¿¿
V Do
ℜ= ……………………………………………….(Eq.4.11)
v
(Table:6-8: J.P. Holman)
Where:
By substitution :
10.2778∗2.048674976
ℜ= −5
=1.276927∗106 Ra<10 7
1.6603∗10
1 1
0.62 ( 1.276927∗106 ) 2 0.6965183 3
Nu=0.3+ ¿¿
Nu=1507.712882
1507.712882∗0.0274687 W
h0 =
27 2.048674976
=20.2155 2
m .K
Calculate of hi heat transfer coefficient of mixture :
Nu∗k
hi = … … … … … … … … … … … … … … … … … … … …(Eq .4 .12)
Di
1
β= … … … … … … … … … … … … … … … … … … … … …(Eq .4 .16)
T film
T w+ T D
T film= … … … … … … … … … … … … … … … … … … … .( Eq .4 .17)
2
28
β: Thermal expansion coefficient for mixture
By substitution :
( 41.7+273.15)+(53.35+273.15)
T flim= =320.6625 K
2
1
β= =¿0.003118544 K-1
320.6625
Gr=¿ 5.079*1013
Ra=Gr * Pr = 1.868*1014
29 Nu=0.13 ( 1.868∗1014 )
0.33
Nu =7431.676
7431.676∗0.092427172
hi =
2.0330654
W
hi=337.8587
m2 . K
D o 2.048675
r2 = = =1.02434 m , Ao =
2 2
π*Do*l
Ao = 3.14*2.048675*10.16533 =65.3919 m2
Di 2.0330654
r1 = = =1.01653 m , Ai = π *Di*l
2 2
Ai = 3.14*2.0330654*10.16533 = 64.9266 m2
T f −T w
q without Insulation=
r2
ln ( )
1 r1
+
h i Ai 2 π∗k m∗L
30
By substitution :
35.35−41.7
q with out Insu .=
1.02434
ln ( )
1 1.01653
+
337.8587∗64.9266 2 π∗54∗10.16532705
T f −T air
q with Insu .= … … … … … … … … … … … .(Eq .4 .18)
r2 r3
ln ( ) ln ( )
1 r1 r2 1
+ + +
hi A i 2 π∗k m∗L 2 π∗k ¿∗L ho A o
Where:
By substitution :
31
The maximum heat loss should be not exceed 10% then :
q with Insu .=qwith out Insu .∗0.1
53.35−30
24422.27=
1.02434 r
ln (
) ln ( 3 )
1 1.01653 r2 1
+ + +
337.8587∗64.9266 2 π∗54∗10.16533 2 π∗0.038∗10.16533 20.2155∗2 π∗r 3∗10
By solved
r3 = 1.02472 m
Thickness insulation = r3 – r2
32
4.3. Weight of loads :
4.3.1. Weight of ladders (Wl):
By substitution :
contour=2∗3.14∗ ( 2.048675
2 )=6.4361m
L=( 0.25∗6.4361 ) +1.28=2.88903m
W l =150∗2.88903=433.354 N
Where:
33
V ins : surface area of insulation of length m³
= 0.02481 m3
Where:
D0=¿ 80.612 in
34
a=2
W H =488.037 ¿
W H = 1290.62 N
By substitution :
WD = 10479.7N
35
4.3.6. Life weight ( Wf ) :
W f = ρf ∗V f ∗g … … … … … … … … … … … … … … … … … … … … …(Eq .4 .25)
Where:
π π
V f : ∗Di2∗l = ∗( 2.0334 )2∗10.1670=33 m3
4 4
By substitution :
W f =530.16∗33∗9.81=171621 N
W T =W f +W D ……………………………………….. ( Eq .4 .26)
36
By substitution :
W T =171621+10479.7
W T =182101 N
4.4. The calculation of wind :
Where:
m
uw :Velocity of wind (10.2778 )
sec
By substitution :
F w∗X 2
M x= ………………………………………………( Eq .4 .29 )
2
Where:
By substitution :
10.7791∗10.1653 2
M x= = 556.925 N.m
2
P D Di
σ h= ……………………………………………… ( Eq .4 .30 )
2t v
By substitution :
0.605∗2.0330654
σ h= = 78.7981 N/mm2
2∗7.8047833/1000
38
4.5.2. Circumferential stress (σ l , N /mm2):-
P D∗Di
σl= …………………………………………
4tv
( Eq .4 .31 )
By substitution :
0.605∗2.0330654
σl= = 39.3991 N/mm2
4∗7.8047833/1000
W
σ w = ¿ ¿D …………………………………………………( Eq .4 .32 )
By substitution :
M x Di
σ b=± ( +t ) ………………………………………………( Eq .4 .33 )
Iv 2 v
39 I v=
π
( D4 + Di4 )…………………………………………………( Eq .4 .34 )
64 0
Where:
I v :Second moment of area of the vessel about plant of bending.
By substitution :
π
I v= ( 2048.6754 +2033.0654 4 ) =1.7033*1012mm4
64
556.925∗1000 2033.0654
σ b=± ( +7.8047834)
1.7033∗1012 2
By substitution :
At the up wind ( + σ b ) .
By substitution :
σ max=78.7981−(39.189)=39.6088 N /mm ²
Check :
2∗104∗t v
σ c= ………………………………………………………( Eq .4 .37 )
41 Do
2∗104∗7.8047833 2
σ c= =76.1934759 N / mm
2048.675
Check :
42
Figure 4.2. Horizontal cylindrical vessel on saddle supports
Figure 4.3 : Typical skirt-support designs. (a) Straight skirt. (b) Conical skirt.
Figure
4.4:
43
Bracket
supports. (a) Supported on legs. (b) Supported from steelwork.
4.7. Shirt thickness :
4∗Mxs
σ bs = … … … … … … … … … … … … … … … … … … … … ( Eq.4 .38 )
π∗(Ds+ts)∗ts∗Ds
Fw
Mxs= ∗( L+ Ls )2 … … … … … … … … … … … … … … … … … … … … … … … ( Eq .4 .39 )
2
Where :
Assume: t s =t v (¿ 7.805mm)
By substitution :
44 M xs =
10.77913
∗( 10.1653270+ 2.8873 )2= 918.2262N.m or 918226.2 N.mm
2
4∗918226.2
σ bs = =0.036102 N/mm2
π∗(2033.0654+ 7.805)∗7.805∗2033.0654
WD
σ ws ( test )= … … … … … … … … … … … … … … … … … … ( Eq .4 .40 )
π∗( D s +t s )∗t s
By substitution :
10479.7059
σ ws (test )= =¿ 209422.3 N/m2
π∗(2.0330654+.007805)∗.007805
Wf
σ ws ( operation )= … … … … … … … … … … … … … … … … … ( Eq .4 .41 )
π∗( D s +t s )∗t s
By substitution :
45 σ ws (operation)=
171621.237
π∗(2.0330654+ 0.007805)∗0.007805
=¿ 3429611 N/m2
4.7.3. The Resultant Stress on the Skirt:
By substitution :
Check :
46
Where:
By substitution :
= 96911.9 N/mm2
Check :
47
Where:
Assume : 𝐺 = 0.5 kg/sec
𝜌𝑚 : density of mixture from table 3.4 = 530.137 kg/m3
By substitution :
Figure 4.5. Flange types. (a) Welding-neck. (b) Slip-on. (c) Lap- joint. (d) Screwed.
48
Figure 4.6 Typical standard flange design (All dimensions mm).
49
Chapter five:
Cost estimation
50
5.1. Storage Tank Cost Estimation :-
Diameter of vessel=2.0330654 m
Pressure Design=6.05 ¯¿
Presser factor=1.1 ¿( Figure 5.1)
Material factor for Carbon steel=1.0
Purchase cost = (base cost from figure) * Material factor *Pressure factor……………( Eq .5.1 )
52
5.3. Total Cost Estimation:-
The total cost of design is the summation of vessel cost and insulation cost
Total cost= cost of vessel +cost of insulation………………………………………………..
( Eq .5.3 )
IndexValue(2019)
Total cost (2019) = Total cost (2004) * ……………………………………………….
Index Value(2004)
( Eq .5.4 )
Table 5.1. Inventory Index Factors for 2019. Index Code: A , Average For All
53
Index Value ( 2019 ) =1716.2
54
References
1) “Design of pressure vessel project report” Submitted by MIJO JOSEPH.
12) Index values as published by Marshall & Swift Valuation Service 01/201