Perancangan Menara Distilasi
Perancangan Menara Distilasi
Perancangan Menara Distilasi
QC
D, XD
MD-301
F, XF
QB
B, XB
Keterangan :
F = umpan masuk
B = hasil bawah
D = hasil atas
2
Dalam perancangan menara distilasi ini dipilih jenis Tray dengan pertimbangan
diameter kolom lebih dari 3 ft (0,91 m) (Walas, 1990). Sedangkan jenis tray
yang digunakan adalah sieve tray dengan pertimbangan:
1.) Pressure drop rendah dan efesiensi tinggi (tab. 9.22, ludwig, 1980).
2.) Lebih ringan, murah karena pembuatannya lebih mudah.
3.) Biaya perawatan murah karena mudah dibersihkan.
C. Kondisi Operasi
Langkah-langkah perhitungannya adalah sebagai berikut :
1.) Menghitung kondisi operasi atas dan bawah menara distilasi.
2.) Menentukan Volatilitas Rata-rata
3
3.) Mengecek pemilihan light key component (LK) dan heavy key component
(HK) dengan persamaan Shira’s et. Al. pada Rm.
x j .D .D 1x
j LK . D .D LK
j x HK .D .D
x j .F .F LK 1x LK .F .F LK 1x HK .F .F
( Treybal,1981pers.9.164)(F.1)
Batasan:
Komponen j tidak terdistribusi jika:
x j , D .D x j , D .D
0,01 atau 1,01 , atau (F.2)
x j , F .F x j , F .F
N B x x LK , B
2
log r 0,206. log HK , F .
Ns D x LK , F x
HK , D
(Coulson, 1983, pers. 11.62)(F.7)
4
33406,48 kg/jam
Umpan dalam kondisi cair jenuh. Untuk menentukan temperatur umpan maka
perlu ditrial temperatur bubble point feed pada tekanan 9 atm. Tekanan uap tiap
komponen dihitung dengan menggunakan persamaan:
log10(P) = A + B/T + C LogT + DT + ET2 (Yaws, 1996)
keterangan:
A, B, C,D, E = konstanta
P = tekanan uap komponen i (mmHg)
T = temperatur (K)
Komponen A B C D E
C2H4Cl2 92,355 -6920,4 -10,651 9,1426E-06 2
C2H3Cl 91,432 -5141,7 -10,981 0,000014318 2
HCl 105,16 -3748,4 -15,214 0,031737 1
Sumber : (Chemical Properties Handbook ; Carl L Yaws)
Dengan cara trial T pada tekanan, 9 atm hingga yi = 1 maka akan diperoleh
Massa Masuk
Komponen Pi (atm) K=Pi/P yi= Ki x xi
kg/jam kmol/jam fraksi (xi)
C2H4Cl2 13341,5 134,82 0,2494 0,1162 0,0058 0,0014
C2H3Cl 12651,6 202,43 0,3745 4,1846 0,2093 0,0784
HCl 7413,42 203,32 0,3761 48,9158 2,4465 0,9202
Jumlah 33.406,48 540,57 13341,5 53,2166 2,6617 1,0000
P = 9 atm
Trial temperatur digunakan metode goal seek pada program Ms. Excel, dengan
Pada keadaan dew point, xi = (yi/Ki) = 1. Dengan cara trial T pada tekanan
tabel berikut.
6
Massa
Komponen Pi (atm) K=Pi/P Xi= yi/Ki
kg/jam kmol/jam fraksi (yi)
C2H4Cl2 0 0 0,0000 0,0027 0,0003 0,0000
C2H3Cl 12,65 0,202 0,0010 0,3978 0,0442 0,0225
HCl 7406,00 203,12 0,9990 9,2039 1,0227 0,9775
Jumlah 7418,65 203,32 1,0000 9,6044 1,0672 1,0000
P = 9 atm
Trial temperatur digunakan metode goal seek pada program Ms. Excel, dengan
Pada keadaan bubble point, yi = (Ki x xi) = 1. Dengan cara trial T pada
sebagai berikut:
Massa Masuk
Komponen Pi (atm) K=Pi/P yi= Ki x xi
kg/jam kmol/jam fraksi (xi)
C2H4Cl2 13341,5 134,82 0,3998 0,8404 0,0934 0,0373
C2H3Cl 12638,9 202,23 0,5996 14,4464 1,6052 0,9625
HCl 7,41 0,203 0,0006 138,9918 15,4435 0,0000
Jumlah 25987,83 37,246 1 154,2786 17,1421 1,0000
P = 9 atm
Trial temperatur digunakan metode goal seek pada program Ms. Excel, dengan
x j,D D
Komponen i terdistribusi jika: -0,01 ≤ ( ) ≤ 1,01
x j,F F
8
x iD .D x .D
Komponen i tak terdistribusi jika: < -0,01 atau iD >1,01
z iF .F z iF .F
Pengambilan light key dan heavy key tepat karena dari hasil perhitungan di atas
C2H3Cl dan HCl terdistribusi.
𝑥 𝐿𝐾 𝑥
log ∙ 𝐻𝐾
𝑥 𝐻𝐾 𝐷 𝑥 𝐿𝐾 𝐵
𝑁𝑚 = Coulson ,1983,pers .11.58) (F.17)
log 𝛼 𝑎𝑣𝑒 ,𝐿𝐾
Nm = 28 Plate
Keterangan:
Karena feed yang masuk adalah liquid pada boiling point, maka q = 1. Subsitusi
persamaan (12) menjadi:
𝛼𝑖 ×𝑥𝑖,𝐹
=0 (F.19)
𝛼𝑖 −𝜃
Untuk menghitung nilai refluks minimum dicari dengan cara trial nilai sampai
diperoleh nilai persamaan diatas sama dengan nol.
Keterangan :
Rm = Refluks minimum
R = Refluks
i xi , F
Nilai ditrial hingga 0. Nilai harus berada di antara nilai
i
Maka :
i xi , D
Rm + 1
i
1,1925 = Rm + 1
Rm = 0,1925
R operasi = 0,231
R op 0,231
= = 0,187
R op +1 0,231 + 1
Rm 0,1925
= 0,1925+1 = 0,161
R m +1
11
Nm
0,2
N
Nm = 5,893 Plate
N = 5,893 Plate
0,2
x B
top avg = x = 0,013 cp
B
bottom avg =
x B
= 0,0879 cp
x
B
Keterangan :
(Xlk, Xhk)F : Fraksi mol light key dan heavy key di Feed
Diketahui :
D (kmol/jam) = 203,324
B (kmol/jam) = 210,855
XHK F = 0,3745 XHK d = 0,001
XLK F = 0,3761 XLK b = 0,00096
28
Ns =
1 0,9968
Ns = 17,52 plate plate ke 18 (tidak termasuk reboiler)
Suhu distilat dicari secara trial dengan cara dew point calculation seperti
persamaan (1). Suhu uap dari puncak menara, V1, dicari dengan menggunakan
dew point calculation seperti persamaan (2).
Tref = 25 oC (298,15 K)
Tref = 25 oC (298,15 K
Tref = 25 oC (298,15 K)
Q = n Cp dT
Vapor juga memiliki energi laten (perubahan fase liquid menjadi vapor).
Perhitungan Energi laten vapor ditabulasikan pada tabel di bawah ini
V. Hv = L. HL + D. HD + ∆H cond.
∆H cond = V. Hv – L. HL – D. HD
= 4.156.712,8170 – [(-19.6054,6991) + (-140.705,8874)]
= 4.493.473,4035 kJ/jam
Suhu bottom dan komposisi uap yang berada pada kesetimbangan dengan bottom,
Vn+1, dicari secara trial dengan bubble point calculation seperti pada persamaan
(1).
Tref = 25 oC (298,15 K)
∆H in = ∆H out
∆H umpan + ∆Hreboiler = ∆H bottom + ∆H distilat + ∆H condensor
17
ρ liquid =
y D
y
D
= 478,933 kg/m3
Densitas uap
P = 9 atm = 911925 Pa
R = 8.314,34 m3.Pa/kmol.K
ρ uap = (BMxP)/RT
18
ρ uap = y
D
y
D
= 16,778 kg/m3
ρ liquid = x
B
x
B
= 351,206 kg/m3
Densitas uap:
P = 9 atm = 911925 Pa
R = 8314,34 m3.Pa/kmol.K
ρ uap = (BMxP)/RT
19
ρ uap = x B
x
B
= 3,5155 kg/m3
𝑃𝑐ℎ 𝜌 𝐿 −𝜌 𝑣 4
𝜎= 𝑥 10−12 (Coulson ,1983,pers .8.23.hal .258) (F.22)
𝑀
Keterangan :
Cl 54,3000
double bond, '=' 23,2000
rangkap tiga 46,6000
Tinggi plate spacing pada umumnya antara 0,3 – 0,6 m (Coulson, 1983, hal.448).
Diambil plate spacing = 0,30 m.
L’ = L + qF = 35.120,51 kg/jam
V’ = L’ – B = 9.132,68 kg/jam
L’-V’ = 25.987,83 kg/jam
L’/V’ = 1,351
1.714,03 kg / h 16,778kg / m 3
FLV , top
9.132,68 kg / h 478,933 kg / m 3
= 0,0352
Untuk tray spacing = 0,3 m dan FLV = 0,0352 maka K1 = 0,065 (Coulson,
1983, fig.11.27).
top
0, 2
K '
1 top K1
0,02
K1’top= 0,0288
22
1.714,03 kg / h 23,877kg / m 3
FLV , bottom
9.132,68 kg / h 351,21 kg / m 3
= 0,0489
Untuk tray spacing = 0,3 m dan FLV = 0,048 maka K1 = 0,06 (Coulson, 1983,
fig.11.27).
Koreksi nilai K1 bottom:
0, 2
K '
1 bottom K1 bottom
0,02
K1’bottom= 0,0251
K1 = Konstanta
u f , top 0,0288
478,93kg / m 16,78kg / m 3
3
16,78kg / m 3
= 0,1511 m/s
23
u f , bottom 0,0251
351,21kg / m 23,877 kg / m 3
3
23,877 kg / m 3
= 0,0928 m/s
9.132,68 kg / h
Qv , bottom
23,877 kg / m 3 x 3600 sec/ h
= 0,106 m3/s
Menentukan Luas Area Netto Untuk Kontak Uap-Cair
Qv
An
uv
Keterangan :
An = Luas Area Netto (m2)
Qv = Laju alir volumetrik (m3/s)
uv = Kecepatan uap (m/s)
= 1,25 m2
Luas area netto bagian bawah (bottom) :
0,106 m 3 / s
An , bottom
0,0742m / s
= 1,145 m2
1,25 m 2
=
(1 0,2)
= 1,563 m2
25
An , bottom
Ac , bottom
1 Ad
1,145 m 2
=
(1 0,2)
= 1,432 m2
4. Ac
Menghitung diameter menara : Dc
Diameter menara bagian atas (top) :
4 1,563 m 2
Dc , top
3,14
= 1,262 m
4 1,431 m 2
Dc , bottom
3,14
= 1,35 m
35.120,51 kg / jam
QL,B =
351,21 kg / m 3 3600 sec/ jam
= 0,028 m3/s
Keterangan:
QL.bot = laju alir volumetrik bagian bottom (m3/s)
Lw = laju alir massa cairan bagian bottom (kg/s)
ρL = densitas cairan bagian bottom (kg/m3)
26
dari fig. 11.28 (Coulson, 1983) untuk QL,B = 0,0070 m3/s maka jenis
alirannya adalah cross flow (single pass).
Perancangan Tray
Diameter menara , Dc = 1,35 m
Luas menara, Ac (/4xDc2) = 1,495 m2
Luas downcomer,Ad = 0,12 Ac = 0,179 m2
Luas netto, An = Ac - Ad = 1,315 m2
Luas aktif, Aa = Ac – 2,Ad = 1,136 m2
Luas hole, Ah = 0,03,Aa = 0,034 m2
Dari fig,11.31 (Coulson, 1983), untuk Ad/Ac = 0,12 maka :
lw/Dc = 0,72
Panjang weir,lw = 0,72 x Dc
= 0,72 x 1,379 m
= 0,994 m
27
Tebal Tray
Untuk bahan carbon steel tebal plate yang digunakan adalah 5 mm (3/16 in),
sedangkan untuk bahan stainless steel tebal plate yang digunakan adalah 3
mm.
Untuk menara distilasi ini digunakan bahan stainless steel, sehingga tebal
plate yang digunakan = 3 mm. (Coulson vol 6 1ed p465, 1983).
28
= 9,361 mm liquid
0,381 kg / s
how min = 750
478,93 kg / m 0,994 m
3
= 8,07 mm liquid
Keterangan:
u h = kecepatan uap minimum desain, m/s
K2 = konstanta
dh = diameter hole, mm
v = densitas uap, kg/m3
uh
30,4 0,9025,4 5 mm
= 1,866 m/s
16,78 kg / m 3 1/ 2
uh
30,5 0,9025,4 5 mm
23,877kg / m 3 1/ 2
= 1,544 m/s
= 4,633 m/s
Dari figure 11.34, untuk ketebalan plate/diameter lubang = 0,6
(Ah/Aa) x 100 = 10
Didapatkan nilai orifice coefficient (Co) = 0,74
𝑈ℎ 2 𝜌 𝑣
ℎ𝑑 = 51 (F.33)
𝐶𝑜 𝜌𝐿
ℎ𝑑 = 29,583 𝑚𝑚 𝑙𝑖𝑞𝑢𝑖𝑑
Keterangan :
uh = Kecepatan uap maksimum melewati hole (m/s)
Ah = Luas hole (m2)
Qv.bottom = Laju alir volumetrik uap maksimum bagian bottom (m3/s)
ρv = densitas uap bagian Bottom (kg/m3)
ρL = densitas liquid bagian Bottom (kg/m3)
Co = orifice coefficient
hd = Dry plate drop (mm)
31
Residual Head
12,5 𝑥 10 3
ℎ𝑟 = (F.34)
𝜌𝐿
ℎ𝑟 = 26,099 mm liquid
Keterangan :
hr = residual head (mm)
ρL = densitas liquid bagian bottom (kg/m3)
Keterangan :
hr = residual head (mm)
hd = dry plate drop (mm)
hw = tinggi weir (mm)
how = tinggi weir liquid (mm cairan)
hT = total plate pressure drop (mm liquid)
uh = 3,255 m/s
hd = 116,273 mm cairan
hr = 35,592 mm cairan
ht = 211,225 mm cairan
Pt 9,81 x 10 3 ht L
ΔPt = 727,74 Pa = 0,0072 atm
32
hb <
1
lt hw , maka pemilihan plate spacing 0,3 m sudah tepat
2
Menara bagian bawah
hap = 43 mm cairan
Aap = hap x lw = 0,0481 m2
hdc = 5,382 mm cairan
hb = 173,968 mm cairan
hb = 0,174 m
1
lt hw 0,175 m
2
hb <
1
lt hw , maka pemilihan plate spacing 0,3 m sudah tepat
2
Check Entrainment
Actual Percentage Flooding For Design Area
Layout Tray
50 mm
Dc =
0,8309
50 mm
c
1,35
Lw =
Gambar F.9. Layout plate.
( ) = 180 - c
= 180 - 90 = 90o
lh/Dc = 0,2
Lav = Dc 0,05
180 (F.38)
= 2,088 m
(Acz) = 0,089 m2
= 0,9222 m2
Ah
0,0369
Ap
Ip
= 3,8
dh
Hole pitch
Ip
(Ip) = dh (F.57)
dh
= 3,8 5 mm = 19 mm (0,019 m)
Luas 1 lubang = d h2 (F.58)
4
= 0,785 (5 mm)2
Ah
Jumlah lubang = (F.59)
luas 1 lubang
0,0341 m 2
=
1,9625.10 -5 m 2
= 1.736,53 buah
= 1737 buah
37
Spesifikasi Tray
OD
b = tinngi
icr dish
OA
B A
sf
ID t
a
r
Keterangan :
th = Tebal head (in)
icr = Inside corner radius ( in)
r = Radius of dish( in)
sf = Straight flange (in)
OD = Diameter luar (in)
ID = Diameter dalam (in)
b = Depth of dish (in)
OA = Tinggi head (in)
P.ri
t c
f .E 0,6.P (Brownell & Young,1959, pers. 13.11)
158,717 psi26,582 in
11500 psi0,8 (0,6 158,717 ) 0,125 in
=
= 0,513 in
1 r
w . 3 c = 4,382 in
4 icr
P.rc .w
th c = 1,65 in
2 f 0,2 P
t head standar = 1,75 in
maka tebal yang digunakan :
t head = 1,75 in
Untuk tebal head 1,75 in, dari tabel 5,8 Brownell and Young maka
sf = 1,5 – 4,5 in.
Diambil sf = 3 in
b = 9,814 in
= 0,249 m
Dari Tabel 5.8 Brownell & Young, untuk ketebalan head = 1,75 in,
sf = 1,5 – 4,5 in dan dipilih sf = 3 in.
OA = th + b + sf (Brownell and Young,1959)(F.62)
= (1,75 + 0,249+ 3) in
= 4,999 in (0,127 m)
AB = ri – icr
= 26,58 – 3,25
= 23,332 ( 0,593 m)
BC = rc – icr
= 54 – 3,25
= 44,75 in (1,137 m)
AC = BC 2 AB 2
= 44,75 2 23,332 2
= 38,156 in (0,9699 m)
Tinggi Menara
Data perhitungan :
Diameter kolom (Dc) = 1,35 m
Luas kolom (Ac) = 1,495 m2
Volume head = 0,000049 Di3
= 0,000049 (1,38)3
= 0,000121 m3
Volume head pada sf = π/4 ID2 sf
= 3,14/4 (53,165)2 (3)
= 6.656,45 in3
= 0,109 m3
Volume total head = Vhead tanpa sf+ Vhead pada sf
= 0,000121 m3 + 0,109 m3
= 0,1092 m3
41
a) Pipa feed
1
ρ mix
wi
Σ ρi
= 399,1204 kg/m3
x B
x
µmix =
B
= 1,057E-05 Pa.s
= 1,057E-02 cp
0,52
Di ,optimum = 352,8 G .μ0,03. ρ-0,37 (Coulson 5.13 Vol. 6., 1983)(F.66)
= 86,816 mm = 3,418 in
ID = 4,026 in = 0,1023 m
OD = 4,5 in
ρ = (BM .P)/(R.T)
P = 9 atm
R = 0,08026 m3.atm/kmol.K
Tabel F.27. Densitas dan Viskositas Gas Keluar dari Puncak Menara
Vapor µ
Komponen wi ρkg/m3 wi/ρ Wi/µ
(kg/jam) (Pa.s)
C2H4Cl2 0.0000 0.0000 40.9540 0.0000 0.00001 0.0000
C2H3Cl 12.6500 0.0017 25.8648 0.0001 0.00001 197.07
HCl 7,406.0000 0.9983 15.0892 0.0662 0.00001 76,845.15
Total 7.418,65 1,0000 0,0662 77.042,223
1
ρ mix (Coulson, 1983, pers.8.2) (F.67)
wi
Σ ρi
1
=
0,0662
= 15,099 kg/m3
Keterangan :
P = Tekanan (atm)
T = Temperatur (K).
= 1,298E-05 cP
Diameter Optimum :
0,52
Di ,optimum = 352,8 G .μ0,03. ρ-0,37 (Coulson 5.13 Vol. 6., 1983) (F.68)
= 134,28 mm = 5,287 in
ID = 6,065 in = 0,154 m
OD = 6,625 in
A = 28,9 in2
45
1
ρ mix
wi
Σ ρi
1
= = 478,933 kg/m3
0,0021
= 0,013 cP
0,52
Di ,optimum = 352,8 G .μ0,03. ρ-0,37
= 6,807 mm = 0,268 in
ID = 0,269 in
OD = 0,405 in
A = 0,058 in2
46
Tabel F.29. Densitas dan Viskositas Cairan Keluar dari Dasar Menara
Bottom µ
Komponen wi ρkg/m3 wi/ρ Wi/µ
(kg/jam) (Pa.s)
C2H4Cl2 13.341,51 0,5134 344,71 0,0015 1,259E-05 40.757,69
C2H3Cl 12.638,91 0,4863 358,31 0,0014 1,298E-05 37.444,81
HCl 7,41 0,0003 406,33 0,0000 1,927E-05 14,804
Total 25.987,83 1,0000 0,0029 78.217,24
1
ρ mix
wi
Σ ρi
1
=
0,0029
= 351,21 kg/m3
= 0,0128 cP
0,52
Di ,optimum = 352,8 G .μ0,03. ρ-0,37
ID = 0,493 in
47
OD = 0,675 in
A = 0,192 in2
ρ = (BM .P)/(R.T)
P = 9 atm
R = 0,08026 m3.atm/kmol.K
Temperatur = 393,31 K
1
ρ mix
wi
Σ ρi
1
=
0,0983
= 10,175 kg/m3
= 1,925E-02 cP
48
0,52
Di ,optimum = 352,8 G .μ0,03. ρ-0,37
= 175,19 mm = 6,897 in
ID = 7,981 in
OD = 8,625 in
A = 50 in2
Densitas ( ) = 18 lb/ft3
r11 r2 r3
T1
T2
r1 T3
r2 Tu
r3
Perpindahan panas melalui tiap lapis tahanan dihitung dengan hukum Fourier
dan A = 2πrL, diperoleh:
2L (T1 Tu )
Q (Holfman, 1997, pers.2-9) (F.69)
ln r2 ln 3
r
1 r2
r
k1 k2
Jika perpindahan panas disertai konveksi dan radiasi, maka persamaan di atas
dapat dituliskan:
2L (T1 Tu )
Q (Holman, 1997, pers.2-12) (F.70)
ln r2 ln 3
r
1 r2
r 1
k1 k2 h c h r r3
Jika diaplikasikan dalam perhitungan perancangan tangki maka diperoleh:
50
2L (T1 Tu )
Q= (F.71)
r x is
ln r2 ln 2
1
r r2 1
kp k is h c h r (r2 x is )
Keterangan :
Untuk menghitung perpindahan panas dari luar ke dalam shell, harus dihitung
terlebih dahulu temperatur kesetimbangan radiasi pada permukaan dinding luar
yang terkena sinar matahari pada temperatur udara lingkungan sekitar shell.
q
sun low temp . T Tsurr
4 4
(J P Holman, 1979, 6th ed)(F.72)
A sun
Keterangan ;
51
q
= Fluks radiasi matahari (W/m2)
A sun
W
σ = Konstanta Boltzman = 5,669 x 10-8
m K
2 4
Data perhitungan :
r1 = 54 in (2,25 ft)
kp = 25,7305 Btu/hr.ft2 oF
Isolasi yang digunakan akan di lapisi dengan cat (pigmen) berwarna putih.
Berdasarkan Tabel 8.3 (Holman,1979), diperolah data :
q
= 500 W/m2
A sun
surya = 0,18
W
σ = 5,669 10-8
m K
2 4
52
W
500 2 0,18 0,8 5,669 10
8 W 4
[T 303,15 ] K
4 4
m m K
2
T3 = 323,8620 K
= 50,7120 oC
Ti 4 Tu4
hr
T1 T2
= 3,9429 W/m2.K
= 0,6944 Btu/hr.ft2 oF
Keterangan :
= ½ (323,8620 + 308,15)
= 316,0060 K
53
ρf = 1,1201 kg/m3
µf = 1,923x10-5 kg/m.s
kf = 0,0274 W/m K
β = 3,1714x10-5 1
K
L3ρ 2 gβΔT
N Gr ( SI ) (Geankoplis, 1993, Pers.4.7-4)(F.75)
μ2
=
16,189 m3 1,1201 kg / m 3 2 9,8067 m / s 3,1714.10 5 K 1 323,8620 308,15o K
0,00002 kg / m.s
2
= 7,031E+10
C p .
N Pr (Geankoplis, Pers.4.7-4.1993)(F.76)
k
=
1,0056
kJ / kg.o K 0,00002kg / m.s
0,0274 W / m.o K
= 0,70446
N Ra N Gr N Pr (Geankoplis, Pers.4.7-4.1993)(F.77)
= 7,031E+10 0,70446
= 4,953E+10
Berdasarkan Tabel 4.7-2 (Geankoplis,1993, hal. 256), untuk silinder vertikal dan NRa = >
9
10 , maka koefisien perpindahan panas konveksi dirumuskan sebagai berikut :
hc 1,24 . T 1 / 3
(F.78)
= 1,24T3 Tu
1/ 3
= 1,24323,8620 308,151 / 3
54
= 1,2413 Btu/hr.ft2.oF
= 11.710 r3……………………………………………………..(1)
2L (T1 Tu )
qc
ln r2 ln 3
r
1 r2
r 1
k1 k2 h c h r r3
2 53,113 ft (123,2816 95) o F
ln 2,854 ft ln r 3 ft
2,75 ft 2,854 1
25,7305 btu / hr. ft . F 0,0530 btu / hr. ft . F 1,2413 btu / hr. ft 2 .o F r3
o o
11.109,6590
ln 3
r
2,0208 1
0,0004
0,0530 1,2413 r3 ………………………………….(2)
qr = qc
11.109,6590
13783,8236 r3
ln 3
r
2,0208 1
0,0004
0,0530 1,2413 r3
Dengan substitusi pers. (1) ke (2) maka diperoleh nilai jari - jari isolator (r3 )
adalah 2,9479 ft
55
xis = r3 – r2
= 2,1146 ft – 2,0208 ft
= 17.260,05 Btu/jam
A = Surface of head
= 17,235 ft2
= 605,068 Btu/jam
= 18.470,2 Btu/jam
56
Perhitungan awal tebal shell dan head menara telah dilakukan. Menara cukup
tinggi sehingga perlu dicek pengaruh angin dan gempa,
Spesifikasi menara:
OD shell = 54 in = 4,5 ft
OD
Diameter,d = OD + + 2 sf + 2 icr (F.80)
24 3
54
= 54 + + 2(3) + 2 (3,25)
24 3
= 64,42 in
d 2 t
Beban head =
4 1728 (F.81)
= 577,296 lb
fwx
atau fap fap
fsx
fwx
atau
fdx fdx
fsx
di = 1,379 in
ts = 0,25 in
Pd
f ap (Pers. 3.13, Brownell, 1959)(F.82)
4 t s c
158,72 66
f ap =4.285,365 psi
4 0,625 0,125
keterangan :
c = corrosion allowance, in
Shell
Diketahui :
Wshell .( Do2 Di2 ). s . X (Pers. 9.1, Brownell, 1959)(F.83)
4
58
Isolator
Diketahui :
= 0,0937 ft = 1,1250 in
Wins 2
.Dins . X .tins .ins (Pers. 9.2, Brownell, 1959)(F.85)
12
ins . X .tins
fdead ins. = (Pers. 9.4a, Brownell, 1959)(F.86)
144(t s c)
Attachment
53,11
= 3,14 x (4,6042 2 – 4,5 2) x
4
= 39,54 lb/ft
Dm = diameter shell
ts = 0,62500 in
Σ Weight of attachment s
fdead wt attachment . = (F.88)
12. .D m .(t s c)
175,5575 + 70,0045 X
=
12 4,0417 (0,625 0,125)
= 6,6553 + 0,964 X
X 4 X
n 1 1 (F.89)
2 2
X .Dm
1 x 25 x
(liquid trays ) wt 2 4
fdead wt (liquid + trays) = (F.90)
12. .Dm .t s c 12. .Dm. t s c
X
= 1,0417 1
2
= 0,5208 X – 1,0417
15,89 d eff . X 2
fwx = (Pers. 9.20, Brownell, 1959) (F.91)
d o2 t s c
= 69,5 in
fwx = 0,4884 X2
f = 17000 psi
fallowable =f x E
= 10120 psi
fallowable = ft(max)
X2 = a = 0,4884
X = b = -5,1661
c = -4.889,1
x1 = 105,484 ft
x2 = -94,91 ft
61
fc = 1,5 x 106 (t/r) < 1/3 y.p (Pers. 2,25, Brownell, 1959)(F.94)
keterangan :
t = ketebalan shell = 0, 25 in
karena fc lebih kecil dari 1/3 y.p, maka digunakan fc = 13.814,36 psi
fc = fc(max)
X2 = a = 0,4884
X = b = 5,1661
c = -19.045,26
x1 = 192,26 ft
x2 = -202,84 ft
Kondisi ereksi yaitu kondisi tower kosong, tanpa tray, tanpa insulasi, tanpa
tekanan, pipa uap, dan hanya dipengaruhi oleh beban angin.
62
Upwind side,
= 3,907 X + 7,796
deff = 69,5 in
15,89 d eff . X 2
fwx = = 0,4884 X2
d o2 t s c
Upwind side
f = 11650
E = 0,8000
fallowable =fxE
= 10.120 psi
63
X2 = a = 0,4884
X = b = – 0,5073
c = – 10.127,796
X1 = 144,53 ft
X2 = - 143,488 ft
Untuk ketinggian total menara (vessel + skirt ) 63,11 ft, berat menara plus
attachment, liquids, dan lainnya dapat dihitung dengan mengalikan compressive
stress total terhadap berat dengan luas permukaan penampang menara
= 1328,6119 lb
ΣW = fdw(total) π d ts
= 465,427 lb
5. Vibration
Vibrasi ditemui pada menara tinggi. Perioda dari vibrasi pada menara tinggi harus
dibatasi, karena vibrasi yang berlangsung dalam perioda yang cukup lama akan
menimbulkan suatu kerusakan pada menara.
periode vibrasi:
2 1/ 2
H wD
T = 2,65 10 5 (Pers. 9.68. Brownell, 1959)(F.99)
D t
= 0,0863 s
Keterangan :
T = periode vibrasi, s
dari tabel 9.3 Brownell and Young untuk zone 1 & T < 0,0312 s
diperoleh, C = 0,05
65
keterangan :
Msx = 4 C WX 2 3 H X
H2
Msx = 256.345,233 lb
M sx
π r ts c
fsx = 2 (Pers. 9.72. Brownell, 1959)(F.101)
= 886,31 psi
fwx = 2,7426 X2
= 0,5595 x (53,11)2 ft
= 1377,55 psi
fwx > fsx, maka fwx yang mengontrol dan perhitungan pengecekan tinggi menara
benar.
4B
Pa (hal 32, Megyesy, 1983) (F.102)
3(Do / t )
Keterangan :
= 640,659 in
L / Do = 9,4562
Do / t = 38,714
bottom)
Karena nilai dari maximum allowable design pressure lebih besar dari design
shut down), maka pada bagian atas menara dipasang valve yang dihubungkan
Skirt adalah penyangga yang digunakan dan paling aman untuk menyangga
vertikal vessel. Skirt disatukan dengan vessel menggunakan pengelasan
kontinyu (continous welding), ukuran pengelasan ditentukan berdasarkan
ketebalan skirt. Ketebalan dari skirt harus mampu untuk menahan berat mati
dan bending moment dari vessel. Ketebalan skirt harus lebih dari 6 mm.
M Pw .Dis .H .hl
(Megesy, 1983)(F.103)
Keterangan :
= 199.063,056 ft.lb
hT = ketinggian skirt = 10 ft
= 131.157,796 ft.lb
12 M T W
t (Megesy, 1983)(F.105)
R π S E D π S E
2
S = Nilai stress dari head atau material skirt menggunakan bahan stainless steel
= 15.000 psi
Butt Weld
tebal skirt
Vertikal vessel harus merekat erat pada concrete fondation, skirt atau yang lain
dengan anchor bolt dan base (bearing) ring. Jumlah anchor bolt harus 4 atau
kelipatannya untuk setiap vertikal vessel, pada vessel yang tinggi sebaiknya
menggunakan 8 buah anchor bolt. Agar merekat kuat pada concrete fondation,
anchor bolt sebaiknya tidak dipasang terlampau dekat, yakni tidak kurang dari 18
in. Pada vessel diameter kecil agar jarak minimal dari anchor bolt terpenuhi,
sebaiknya menggunakan conical skirt atau wider base ring with gussets, atau
anchor bolt chair.
12M W
T (Megesy, 1983)(F.106)
AB CB
70
keterangan :
= 1.384,7400 in2
= 131,8800 in
T.CB
B4 (Megesy, 1983) (F.107)
SB .N
keterangan :
dari tabel A (Megesy, 1983) untuk area bolt seluas = 2,0490 in2
2,0395 in2
l3 = 2,3750 in
l2 = 1,8750 in
Pada menara distilasi, anchor bolt didesain dengan menggunakan chair agar
lebih kuat dan mampu untuk menahan menara bermuatan berat, digunakan bolt
72
dengan ukuran 1,875 in maka dari tabel standar chair anchor bolt, tabel
berdasarkan Scheiman A.D. Shorts Cuts to anchor Bolting and Base Ring Sizing,
Petroleum Refiner, June 1963. (Megesy hal 76, 1983)
A = 2,625 in E = 1,5 in
B = 5 in F = 2,125 in
C = 3,5 in G = 2,375 in
D = 0,75 in
2
3
1,5
1
1
/4
1
/4
2,5 1/2
1,75
T.CB
SB (Megesy, 1983)(F.108)
B4 .N
= 14.015,0659
shell
skirt
chair
anchor
bolt
bearing
plate
Beban yang ditopang pada skirt, dilanjutkan ke pondasi menara melalui base
ring. Base ring harus cukup lebar agar bisa mendistribusikan beban ke pondasi
secara merata, sehingga cukup kuat untuk menahan beban menara.
12M W
Pc (Megesy, 1983)(F.113)
As Cs
keterangan :
12 x199.063,06 13.286,12
Pc
1.384,7400 131,8800
= 1.825,81 lb/lin-in
Pc
l (Megesy, 1983)(F.114)
fb
keterangan :
= 1.825,81 lb/lin-in
1.825,81
l
750
l = 2,4344 in
Dari tabel A (Megesy, 1983 hal 69) digunakan bolt dengan ukuran 1,8750 in.
l2 = 1,8750 in
l3 = 2,3750 in
75
li = 12 + 13 = 4,25 in
tB = 0,32 x 4,25 in
Data Perancangan :
t h
Gasket
W
hG
hT R hD C
go
HG
HT g1
g1/2
do y pm
= 1,0052 (Pers. 12.2 Brownell, 1959)(F.116)
di y p(m 1)
keterangan :
assumsi :
digunakan gasket dengan tebal 1/16 in, dari fig 12.11 B & Y, diperoleh :
assumsi :
diamater dalam gasket = diameter luar shell, do yaitu = 1,0204 in, sehingga:
do = 1,0204 66 = 67,349 in
G = 67,3489 in + 0,1500 in
= 67,4989 in
N
bo (F.117)
2
77
maka bo = b = 0,0750 in
Wm1 = Hy = b π G y (F.118)
= 58.815,174 lb
Hp =2bπGmp (F.119)
= 11.566,72 lb
π G2
H = p (F.120)
4
2
3,14 x 67,4989
= 132,3
4
= 473.176,327 lb
Wm2 = H + Hp (F.121)
= 473.176,327 lb + 11.566,72 lb
= 484.743,05 lb
Wm2 lebih besar dari Wm1 sehingga Wm2 sebagai beban pengontrol.
Perhitungan luas baut minimum (minimum bolting area) – pers 12.92 B & Y
78
Wm2
Am1 = (F.122)
fb
keterangan :
484.743,05
Am1 =
20000
= 24,237 in2
Dari tabel 10.4, Brownell & Young Digunakan baut berukuran 2 in sebanyak 8
baut. Bolt circle diameter yang digunakan 59,684 in.
C = 59,684 in.
E
d
= 63,434 in = 1,6112 m
79
Ab actual = 2,0490 5
= 9,2205 in2
A b actual f allow
Nmin = (F.124)
2yπG
9,2205 x 17000
=
2 x 3700 x 3,14 x 54,4989
Perhitungan momen
= 284.390,04 lb
= ½ ( 59,6838 – 54,4989) in
= 2,215 in
= 284.390,04 lb x 2,215 in
= 481.002,032 lb in
HD = 0,785 x B2 x p (F.128)
HD = 0,785 x 54 2 in x 37,2530
HD = 302.842,638 lb
hD = ½ (C – B) (F.129)
= ½ (59,6838 - 54)
= 2,8419 in
MD = HD x h D (F.130)
= 302.842,638 lb x 2,8419 in
= 860.640,92 lb in
HG = W – H = Wm1 – H (F.131)
= 9.468,353 lb
= 2,215 in
MG = H G x hG (F.133)
= 20.972,86 lb in
HT = H - HD (F.134)
= 14.224,69 lb
Hubungan lever arm adalah dengan persamaan 12.102 Brownell & Young.
= 2,5285 in
MT = HT x h T (F.136)
= 35.966,58 lb in
= 917.580,37 lb in
t = Y M max (F.138)
fa B
Diketahui:
K = A/B = 1,1747
Bolt
t = tebal flange
Gasket
d = diameter baut
Gambar F.17. Detail untuk flange dan bolt pada head menara.
Setiap pressure vessel yang dalam operasinya melibatkan liquid ataupun vessel
yang di dalamnya terdapat alat lain seperti impeler, packing, plate dan lainnya
sebaiknya dilengkapi dengan manhole yang tujuannya untuk pemeriksaan dan
perbaikan. Untuk vessel dengan diameter antara 48 in sampai 96 in, digunakan
manhole dengan diameter dalam minimal 15 in. (Megyesy, 1983). Manhole dipasang
dengan tujuan sebagai tempat untuk perbaikan plate.
42,25"
12,25"
42,25"
27"
9"
26.25"
1,25"
28,75"
3"
6"
5"
1,25"
2,0764"
20'’
24,5"
Diameter vessel = 54 in
83
Diameter manhole = 20 in
Kontrol Temperatur
Sangatlah penting halnya untuk melakukan pengontrolan dengan ketat pada distilasi
kolom dengan tujuan untuk mengoptimasikan pemisahan, mencegah flooding,
meminimalkan biaya steam, dan memaksimalkan keuntungan.
Bila temperatur top kolom meningkat, maka alat kontrol akan menghasilkan
ukuran temperatur yang kemudian dikirim ke DCS, kemudian DCS akan
memproses data, lalu memerintahkan control valve untuk menambah laju
alir cooling water.
Bila temperatur top kolom menurun, maka yang terjadi adalah sebaliknya,
DCS akan memerintahkan control valve untuk mengurangi laju alir cooling
water.
Kontrol tekanan
Bila tekanan dalam kolom menurun, hingga pada keadaan vacuum (biasanya
terjadi saat pengosongan kolom) maka kontrol tekanan akan mengirimkan ukuran
tekanan ke DCS, selanjutnya DCS akan memerintahkan control valve untuk
terbuka, sehingga tekanan di dalam kolom akan normal.
85
Bila terjadi perubahan ekstrim pada kolom hingga berada diluar batas tekanan
desain kolom, maka kontrol tekanan akan mengirimkan sinyalnya ke DCS,
kemudian DCS akan memerintahkan control valve untuk terbuka.
Kontrol komposisi
Pada menara distilasi, bottom yang dihasilkan berupa EDC 3,9 % , VCM 96 %
dengan sedikit HCl. Sedangkan produk atas yang dihasilkan yaitu VCM 0,01 % dan
HCl 99,9 %. Oleh karena itu perlu dilakukan pengontrolan komposisi pada hasil atas.
Bila HCl pada distilat tidak mencapai mencapai 99,9 %, maka alat kontrol komposisi
akan mengirimkan sinyal ke DCS, kemudian DCS akan memerintahkan control valve
pada refluk untuk menambah laju dari refluk yang masuk ke kolom.
Untuk kontrol temperatur digunakan PID Controller, karena PID memberikan respon
yang paling cepat terhadap penyimpangan keadaan (temperatur) yang terjadi di
kolom.
Dari grafik dapat dilihat bahwa PID Control memberikan cara yang baik untuk
55,75 in
6 in
3
2
1
0,1667 ft
0,75 in
4
35
11
34
10
33
32 9
20
19 8
53,11 ft 4 in
18
5
5
4 14
3
2
1
8 in
12
13
120 in 7
15 0,375 in
1. Shell
2. Flange dan bolt
3. Lubang keluar uap atas
4. Lubang masuk refluks
5. Lubang masuk umpan
6. Lubang masuk dari reboiler
7. Lubang keluar bottom
8. Plate area
9. Sieve tray
10. Apron
11. Weir
12. Manhole
13. Skirt
14. Downcomer
15. Chair anchor bolt
OD = 53,414 in
OA = 4,9999 in
b = 9,814 in
icr = 3,25 in
sf = 3 in
B A
ID = 52,165 in t= 1,75 in
rc = 54 in
3 mm 50 mm
300 mm
5 mm
50 mm