Nama : Gissa Navira Sevie
NIM : 23017001
Tugas 5A –TK 5201
Tanggal tugas : 1-2-2018
Pengumpulan : 21-2-2018
Problem 1. CSTR Unit Tunggal
Tinjau reaksi fasa gas ideal tak-reversible : A + B C
Reaksi dilangsungkan dalam sebuah CSTR isothermal. Umpan adalah campuran 50%-mol A
dan 50%-mol B, dengan konsentrasi masing – masing 50 mol/m3. Waktu ruang adalah 500
detik dan temperatur reaktor adalah 25oC.
𝑚3
Laju reaksi A : (−𝑟𝐴 ) = 𝑘 𝐶𝐴 𝐶𝐵 ,dimana k = 3,5 x10-5 𝑚𝑜𝑙.𝑑𝑒𝑡𝑖𝑘 .
Tekanan umpan sama dengan tekanan [Link] konversi A dan waktu tinggal rata – rata
Jawab:
reaksi : (𝑟𝐴 ) = 𝑘 𝐶𝐴 𝐶𝐵
Data diketahui
𝑚𝑜𝑙
CAin = 50 𝑚3
𝑚𝑜𝑙
CBin = 50 𝑚3
CCin =0
nAin = 0.5 mol
nBin = 0.5 mol
nCin =0
t = 500 detik
𝑚3
parameter kinetika : k = 3,5 x10-5
𝑚𝑜𝑙.𝑑𝑒𝑡𝑖𝑘
Pada kondisi awal saat t = 0
CA = CAin
CB = CBin
CC = CCin
Persamaan pendukung
𝑛 +𝑛
V=𝐶 Ain+ 𝐶Bin
Ain Bin
𝑉
Qin= 𝑡
Q1=Qin
𝐶𝐴
𝑋𝐴 = 1 −
𝐶𝐴𝑖𝑛
Persamaan rancangan:
Komponen A:
Input – output + bereaksi = akumulasi
𝑑(𝑉. 𝐶𝐴 )
𝐹𝐴𝑖𝑛 − 𝐹𝐴𝑜𝑢𝑡 − (𝑟𝐴 . 𝑉) =
𝑑𝑡
𝑑(𝑉. 𝐶𝐴 )
𝑄𝑖𝑛 . 𝐶𝐴𝑖𝑛 − 𝑄1 . 𝐶𝐴 − (𝑟𝐴 . 𝑉) =
𝑑𝑡
𝑄𝑖𝑛 . 𝐶𝐴𝑖𝑛 − 𝑄1 . 𝐶𝐴 𝑑(𝐶𝐴 )
− ( 𝑟𝐴 ) =
𝑉 𝑑𝑡
Komponen B:
Input – output + bereaksi = akumulasi
𝑑(𝑉. 𝐶𝐵 )
𝐹𝐵𝑖𝑛 − 𝐹𝐵𝑜𝑢𝑡 − (𝑟𝐴 . 𝑉) =
𝑑𝑡
𝑑(𝑉. 𝐶𝐵 )
𝑄𝑖𝑛 . 𝐶𝐵𝑖𝑛 − 𝑄1 . 𝐶𝐵 − (𝑟𝐴 . 𝑉) =
𝑑𝑡
𝑄𝑖𝑛 . 𝐶𝐵𝑖𝑛 − 𝑄1 . 𝐶𝐵 𝑑(𝐶𝐵 )
− ( 𝑟𝐴 ) =
𝑉 𝑑𝑡
Komponen C:
Input – output + bereaksi = akumulasi
𝑑(𝑉. 𝐶𝐶 )
𝐹𝐶𝑖𝑛 − 𝐹𝐶𝑜𝑢𝑡 + (𝑟𝐴 . 𝑉) =
𝑑𝑡
𝑑(𝑉. 𝐶𝐶 )
𝑄𝑖𝑛 . 𝐶𝐶𝑖𝑛 − 𝑄1 . 𝐶𝐶 + (𝑟𝐴 . 𝑉) =
𝑑𝑡
𝑄𝑖𝑛 . 𝐶𝐶𝑖𝑛 − 𝑄1 . 𝐶𝐶 𝑑(𝐶𝐶 )
+ ( 𝑟𝐴 ) =
𝑉 𝑑𝑡
Grafik 1. Profil konsentrasi keluar reaktor (CA, CB, CC,) terhadap waktu selama 2000 detik
Grafik 2. Profil konversi A (XA) terhadap waktu selama 2000 detik
a) Didapatkan XA = 0.359183
b) waktu tinggal rata – rata = waktu ruang =500 detik
dikarenakan densitas gas tidak berubah sehingga nilai Q tetap
Skrip Flex PDE
TITLE CSTR Tunggal Gas' { the problem identification }
COORDINATES cartesian2 { coordinate system, 1D,2D,3D, etc }
SELECT
ngrid =1
errlim = 1e-3
penwidth =4
VARIABLES
Ca
Cb
Cc (0.00000001)
DEFINITIONS
!A+BC
Cain = 50 !mol/m3
Cbin =50 !mol/m3
Ccin =0.0000001
nain = 0.5 !mol/m3
nbin =0.5 !mol/m3
ncin =0.0000001
tao = 500 !detik
k = 0.000035
V =(nain+ nbin)/( Cain+ Cbin)
Qin=V/tao
Q1=Qin
ra = k*Ca*Cb
Xa = 1 -Ca/Cain
INITIAL VALUES
Ca = Cain
Cb=Cbin
Cc=Ccin
EQUATIONS
Ca: (Cain*Qin-Ca*Q1)/V-ra = dt(Ca)
Cb: (Cbin*Qin-Cb*Q1)/V-ra = dt(Cb)
Cc: (Ccin*Qin-Cc*Q1)/V+ra = dt(Cc)
BOUNDARIES
REGION 1
START(0,0)
LINE TO (1,0) TO (1,1) TO (0,1) TO CLOSE
TIME 0 TO 2000 { if time dependent }
PLOTS { save result displays }
!for time=endtime
!transfer(Ts) file='[Link]'
for cycle =1
!elevation(NO) from (0,Ly/2) to (Lx,Ly/2)
HISTORIES
history(Ca,Cb,Cc) at (0.5,0.5) as 'konsentrasi'
export file ='[Link]' format 't#r#i'
history(Xa1) at (0.5,0.5) as 'konversiA'
export file ='[Link]' format 't#r#i'
END