Reaktor (R-01) Abi Revisi Sesuai Pak Kun
Reaktor (R-01) Abi Revisi Sesuai Pak Kun
Reaktor (R-01) Abi Revisi Sesuai Pak Kun
REAKTOR (R-01)
Tugas : Mereaksikan Etilen dengan Oksigen dan Asam Klorida dengan
katalis padat Tembaga (II) Klorida yang menghasilkan Etilen
Diklorida dan air serta produk samping berupa gas karbon
monoksida dan gas karbon dioksida
Jenis reaktor : Reaktor non isothermal non adiabatic
Tipe reaktor : Fixed bed multitube
Fase : Gas
Katalisator : CuCl2
Kondisi Operasi:
P = 3.5 atm Produk :
- N2
T = 230 °C
- CO
- O2
- CH4
- C2H4
- C2H6
- Pendingin - HCl
Keluar - CO2
- C2H4Cl2
- H2O
Pendingin
Reaktan : Masuk
- N2
- O2
- CH4
- C2H4
- C2H6
- HCl
- H2O
REAKTOR (R-01)
Dimisalkan:
1
A + 2B + 2 C D+E
A+5B E+2F+2G
Laju reaksi
𝑘 𝑥 𝐾𝑎 𝑥 𝐶𝑒 𝑥 𝐶𝑐
𝑟=
1+𝐾𝑎 𝑥 𝐶𝑒
REAKTOR (R-01)
−37.8
𝑘 = 269 exp( )
1.987×𝑇
𝐾𝑎 = 0.63
(Asian Journal of Applied Sciences Volume 06 – Issue 05, October 2018)
FA |z+ ∆z
z + ∆z
z
FA |z
Keterangan :
FA0 : Laju A masuk reaktor (kmol/jam)
XA : Konversi A
REAKTOR (R-01)
2. Neraca Panas
Asumsi :
Kondisi operasi steady state
Tidak ada gradien suhu dan konsentrasi ke arah radial
Kecepatan umpan masuk reaktor konstan
T|z+ ∆z
z + ∆z
z
T|z
Keterangan :
T : suhu cairan (K)
Tc : suhu pendingin (K)
Fi : laju masing-masing komponen (kmol/jam)
Cpi : kapasitas panas masing-masing komponen ( J/(kmol.K))
ΔHR : panas reaksi (J/kmol)
UD : koefisien transfer massa overall (J/(m2.jam.K))
[Q c z U D . .ID.Nt.(Tc - T)] - Q C z Δz
0
QC QC
z Δz z
U D . .ID.Nt.(Tc - T) 0
Δz
dQ c
lim Δz 0 U D . .ID.Nt.(T - Tc )
dz
dTc
m c .cp c . U D . .ID.Nt.(T - Tc )
dz
dTc U D . .ID.Nt.(T - Tc )
dz m c .cp c
Keterangan:
T = suhu cairan, K
Cpi = kapasitas panas masing-masing komponen, J/kmol.K
∆HR = panas reaksi, J/kmol
UD = koefisien transfer panas overall, J/m2.jam.K
mc = massa pendingin, kg/jam
cpc = kapasitas panas pendingin, J/kg.K
Tc = suhu pendingin, K
REAKTOR (R-01)
Untuk menghitung nilai UD, digunakan persamaan sebagai berikut (Kern, 1965):
UC
UD
1 U C .R d
h io .h o h .ID
di mana, U C dan h io i
h io h o OD
1
k 0,55 cp .μ
= 0,36. c .Re shell . c c
3
De kc
kc = konduktivitas panas pendingin, J/m.s.K
De = diameter ekivalen, m
Reshell = bilangan Reynolds pada shell
μc = viskositas pendingin, kg/m.jam
4. Pressure drop
Penurunan tekanan gas yang mengalir di dalam pipa yang berisi katalisator dapat
dihitung dengan persamaan Tallmagade (Bird, R.b.,”transport phenomena”,ed
2,halaman 191)
𝑑𝑃 1−𝜀 1 − 𝜀 1/6 1 − 𝜀 𝐺𝑝2
= [150 ( ) + 4,2 ( ) ].( 3 ). … … … … … … … (5)
𝑑𝑧 𝑅𝑒𝑝 𝑅𝑒𝑝 𝜀 𝜌. 𝐷𝑝
Dimana:
𝐷𝑝 = Diameter Katalis (m)
REAKTOR (R-01)
clear;clc;
// Initial Condition
dz=0.5;
z=(0:dz:10)';// Panjang pipa reaktor, meter
n=length(z);
Y(1,1)=0;// Konversi
Y(1,2)=230+273.15;// Temperatur Reaktor, Kelvin
Y(1,3)=50+273.15;// Temperatur Pendingin, Kelvin
Y(1,4)=3.5;// Pressure Drop, atm
// BM Campuran (kg/kmol)
BMT=yA*BMA+yB*BMB+yC*BMC+yD*BMD+yE*BME+yF*BMF+yG*BMG+
yH*BMH;
REAKTOR (R-01)
IDs=3.61*0.0254;// m
Pt=1.65*0.0254;// m
// Data Ukuran Katalis jurnal wachi
eps=0.68;// Porositas
Dp=2.59e-5;// m, diameter katalis
rhoK=500;// kgkat./m^3
REAKTOR (R-01)
// Cp pendingin ( kJ/kgram/K)
CpP=(143437+372.894*Tc+0.111497*Tc^2)/1000/166;// Dowtherm A
// Data Densitas (rho = A*B^(-(1-T/Tc)^n), kg/L) Carl Yaws Handbook
rhoA=P*yA*BMA/(0.082*T);// Etilen
rhoB=P*yB*BMB/(0.082*T);// HCl
rhoC=P*yC*BMC/(0.082*T);// Oksigen
REAKTOR (R-01)
// Persamaan Differensial
dxAdz=(rA*3.14*IDt^2*nt/4/FA0);
dTdz=(Ud*%pi*IDt*nt*(Tc-T)*3.6-dHr0*FA0*dxAdz)/FiCpi;
dTcdz=Ud*3.14*IDt*nt*(T-Tc)*3.6/FP/CpP;
dPdz=-f*G^2*(1-eps)/eps^3/Dp/rhoT/1000/101325/3600;
REAKTOR (R-01)
D=[dxAdz,dTdz,dTcdz,dPdz];
endfunction
// Metode RK-4
functionk=K(z, xA, T, Tc, P)
k1=F(z,xA,T,Tc,P);
k2=F(z+dz/2,xA+k1(1)*dz/2,T+k1(2)*dz/2,Tc+k1(3)*dz/2,P+k1(4)*dz/2);
k3=F(z+dz/2,xA+k2(1)*dz/2,T+k2(2)*dz/2,Tc+k2(3)*dz/2,P+k2(4)*dz/2);
k4=F(z+dz,xA+k3(1)*dz,T+k3(2)*dz,Tc+k3(3)*dz,P+k3(4)*dz);
k=k1+2*(k2+k3)+k4;
endfunction
fori=1:n-1
Y(i+1,1:4)=Y(i,1:4)+K(z(i),Y(i,1),Y(i,2),Y(i,3),Y(i,4))*dz/6;
end
disp(' z xA T Tc P')
printf('%10.2f %15.5f %15.5f %15.5f %15.5f \n',z,Y)
figure
plot(z,Y(:,1));xlabel('Panjang reaktor (m)');ylabel('Konversi');
figure
plot(z,[Y(:,2),Y(:,3)]);xlabel('Panjang reaktor (m)');ylabel('Temperatur
(K)');legend('Reaktor','Pendingin');
figure
plot(z,Y(:,4));xlabel('Panjang reaktor (m)');ylabel('Tekanan (atm)');
REAKTOR (R-01)
0.8
y = 3.81E-04x4 - 6.70E-03x3 + 1.60E-02x2 + 2.30E-01x - 7.69E-04
0.6 R² = 9.99E-01 Konversi
Poly. (Konversi)
0.4
0.2
0
0 2 4 6 8 10 12
Reaksi:
Selektivitas C2H4Cl2 = 98 %
Selektivitas CO2 =2%
XA = 0,9958
XA1 (Reaksi Utama) = 0,9958 × 0.98 = 0,97588
XA2 (Reaksi Samping) = 0,9958 x 0,02 = 0,019916
Maka Didapatkan Hasil :
1. Neraca Massa
= 2,58 kg/jam
0,0001
C2H6 masuk reaktor = 0,9997×1000 kmol/jam ×30 kg/kmol × 0.460185
= 1,29 kg/jam
HCl masuk reaktor = 2000 kmol/jam ×36,5 kg/kmol × 0.460185
= 33616.88 kg/jam
O2 masuk reaktor = 600 kmol/jam × 32 kg/kmol × 0.460185
= 8835,55 kg/jam
0,79
N2 masuk reaktor = 0,21×8835,55 kg/jam × 0.460185
= 33238,5 kg/jam
H2O masuk reaktor = 124,45 kmol/jam × 18 kg/kmol
= 2240,06 kg/jam
REAKTOR (R-01)
Neraca Panas
618.84
∆𝐻𝑟
473.15 K
∆𝐻𝑝
∆𝐻𝑟
298 K 298 K
∆𝐻𝑓𝑜 298
ΔHR pada reaksi utama = ΔHF H2O × Xa1× nH2O + ΔHF C2H4Cl2 × nC2H4Cl2
× Xa1 - ΔHF C2H4× nC2H4 × Xa1 - ΔHF HCl × 2 × nHCl × Xa1 - ΔHF O2 × nO2 ×
Xa1 = -107511029.9 kJ/jam
ΔHR pada reaksi samping = ΔHF H2O × 4 × Xa2 × nH2O + ΔHF CO2 × nCO2 × 2
× Xa2 + ΔHF CO× nCO× 2 × Xa2 - ΔHF C2H4 × nC2H4 × Xa2 - ΔHF O2 × nO2 ×
5 × Xa2= -107511029.9 kJ/jam
=-9445311.52 kJ/jam
ΔHR = -116956341.4 kJ/jam
3. Menghitung Viskositas
µ gas = A + B×T + C×T^2
Komponen
A B C
C2H4 -3.985 0.3873 -0.0001122
HCl -9.118 0.555 -0.000111
O2 44.224 0.562 -0.000113
C2H4Cl2 1.0252 0.3179 -0.000041853
H2O -36.83 0.429 -0.0000162
CO2 11.336 0.49918 -0.00010876
CO 35.086 0.50651 -0.00013314
(Yaws , 1997)
REAKTOR (R-01)
𝐼𝐷𝑡 𝑚𝑡̇
𝑅𝑒𝑡 =
𝜇𝐿 𝑛𝑡 𝑎𝑡
𝐼𝐷𝑡 𝑚𝑡̇
𝑛𝑡 =
𝜇𝐿 𝑅𝑒𝑡 𝑎𝑡
0.026645 𝑚 × 57224.4375 𝑘𝑔⁄3600 𝑑𝑒𝑡𝑖𝑘
𝑛𝑡 =
𝑘𝑔
15.969 × 10−6 ⁄𝑚 𝑑𝑒𝑡𝑖𝑘 × 15000 × 0.000557 𝑚2
𝑛𝑡 = 3174.47 ≈ 3200 buah
Dimana K1 dan n1 merupakan konstanta yang ditentukan oleh bentuk pitch dan
jumlah pass yangakan digunakan. Pitch yang digunakan adalah tipe triangular
dan jumlah pass sebanyak satu. Diperoleh nilai K1 = 0.319 dan n1 = 2.142.
1⁄
3180 2.142
𝐼𝐷𝑏 = 0.033528 𝑚 (0.319) = 2.4671 m
𝐼𝐷𝑡
ℎ𝑖𝑜 = ℎ𝑖
𝑂𝐷𝑡
Untuk mengevaluasi nilai hi, digunakan persamaan:
𝑘𝐿
ℎ𝑖 =
𝐼𝐷𝑡 0.027(𝑅𝑒𝑡 )0.8 (𝑃𝑟𝑡 )1/3
Dengan:
𝐶𝑝𝐿 𝜇𝐿
𝑃𝑟𝑡 =
𝑘𝐿
𝑚𝑡̇
𝐺𝑡 =
𝑎𝑡 × 3600
Dimana:
hio = Koefisien transfer panas di tube, kJ/m2-jam-K
kL = Konduktivitas panas campuran, kJ/m-jam-K
CpL = Kapasitas panas campuran, kJ/kmol-K
µL = Viskositas campuran, kg/m-s
Ret = Bilangan Reynold di tube
Prt = Bilangan Prandtl di tube
ODt = Diameter luar tube, m
IDt = Diameter dalam tube, m
b. Kecepatan transfer panas bagian shell (ho)
Pada bagian shellakan dialirkan pemanas berupa Dowtherm Adengan
spesifikasi sebagai berikut, dimana kapasitas panas, viskositas, dan
konduktivitas termalnya merupakan fungsi temperatur dowtherm A:
cp = (143437 + 372.894 * Tp + 0.111497 * Tp ^ 2) / 1000 / 166 (kapasitas
panas dowterm) [kJ/kg.K ]
Rd = Rd shell + Rd tube
= 0.001 + 0.001
= 0.002 jam ft2 °F/BTU
= 9.7838E-05 jam m2 K/kJ
1
𝑈𝑑 = 1
𝑅𝑑 + 𝑈
𝑐
Dimana:
Uc = koefisien transfer panas overall bersih, kJ/jam-m2-K
Ud = koefisien transfer panas overall kotor, kJ/jam-m2-K
REAKTOR (R-01)
1. PIPA (TUBES)
Bahan: Stainless Steel SA–336 Type 316
Susunan tube: TriangularPitch
Spesifikasi ukuran tube:
a. Diameter luar (OD) = 1.049 in
b. Diameter dalam (ID) = 1.32 in
c. Luas penampang = 0.864 in2
d. Luas permukaan luar per panjang = 0.344 ft2/ft
e. Luas permukaan dalam per panjang = 0.274 ft2/ft
f. Berat pipa per satuan panjang = 1.68 lb/ft
g. Pitch (Pt) = 1.65 in
h. Clearance (C’) = 0.33 in
i. Jumlah tube (Nt) = 3200 buah
2. SHELL
Bahan: Carbon Steel SA–53 (Schedule STD)
Dengan spesifikasi sebagai berikut:
Tekanan yang diizinkan (f) = 7985 psi
Efisiensi pengelasan (E) = 0.85
REAKTOR (R-01)
3. HEAD
Untuk tekanan operasi antara 15 – 200 psig, digunakan torispherical head
(flange and dished head). Bahan head diambil sama dengan shell, sehingga
nilai f dan C sama (Brownell & Young, 1979).
OD
icr b
OA
B A
sf
r
ID
t
a
C
Gambar 1. Design Head Reaktor
REAKTOR (R-01)
b. Tinggi head = th + b + sf
sf = 1½ – 4½ in (Tabel 5.4. Brownell & Young, 1979)
icr = 3.75 in
diambil sf = 3 in
a = ½ IDs
= 71.063 in
AB = a – icr
= 67.313 in
r = IDs
BC = r – icr
= 138.376 in
b = r − √BC2 − AB2
= 21.2257 in
REAKTOR (R-01)
7. TEBAL ISOLASI
Perhitungan tebal isolasi yang diperlukan menggunakan asumsi sebagai
berikut:
a. Suhu di dalam reaktor = suhu pemanas rata-rata
b. Suhu dinding luar isolator = 50°C
c. Perpindahan panas pada keadaan steady sehingga q1 = q2 = q3= q4
X1 X2
Keterangan:
R1 = jari-jari dalam shell
R1 R2 = jari-jari luar shell
R2 R3 = jari-jari luar setelah
R3 diisolasi
X1 = tebal dinding
T1 q1 Tu
X2 = tebal isolasi
T2 T3
q4 Tu = suhu udara luar (30°C)
q2 T1 = suhu dinding dalam shell
q3 T2 = suhu dinding luar shell
T3 = suhu dinding isolator
(50°C)
Dengan:
g = percepatan gravitasi bumi (9.8 m/s2)
L = panjang pipa over design reaktor (m)
9.8 × 3.1949 × 10−3 × (50 − 30) × 2.653
𝑅𝑎 = × 0.7051
(1.7601 × 10−5 )2
Ra = 4.3017×1010
Dengan Ra> 109, maka aliran udara adalah turbulen, sehingga untuk
menghitung nilai hc dapat dipakai persamaan berikut (Tabel 7.2
Holman, 1986):
hc = 1.31(∆T)0.33
= 3.5206 W/m2.K
b. Menentukan koefisien perpindahan panas radiasi asbestos-udara
Digunakan persamaan berikut (Holman, 1986):
hr (T3 – Tu) = εσ (T34 – Tu4)
REAKTOR (R-01)
dengan:
σ = Konstanta Boltzman (5.669×10-8 W/m2.K4)
ε = Emisivitas dinding isolator (0.93)
hr = ε·σ (T34 – Tu4) / (T3 – Tu)
= 6.4733 W/m2.K
c. Menghitung panas tiap lapisan
Nilai q4 dihitung dengan persamaan:
q 4 (hr hc).2 .R3 .L.(T3 T4 )
T2 = 567.98686 – R3/24.9569
q4 = q3
0.4462 (𝑇2 − 323.15)
399.756 𝑅3 = 𝑅
3
ln (1.8235 )
3 𝑅
0.4462 (244.83686 − 24.9569)
399.756 𝑅3 = 𝑅 3
ln (1.8235)
R3 = 2.079 m
Sehingga, tebal isolasi dan suhu dinding luar shell didapatkan sebesar:
tisolator = R3 – R2
= 2.079 – 1.8235
= 0.256 m
T2 = 567.98686 – R3/24.9569
= 567.90 K = 294.75 °C
REAKTOR (R-01)
1
ODS = 4,4913 m
16
2 in
4
m
m
HS = 10,769
HR = 12,894
5
9
7
10
REAKTOR (R-01)
RESUME
Jenis Reaktor Fixed bed multitube
Kode R-01
Katalis CuCl2
Pendingin Dowtherm A