Distilasi Dan Absorpsi Gas.pdf

  • Uploaded by: Maulina Aini
  • 0
  • 0
  • December 2019
  • PDF

This document was uploaded by user and they confirmed that they have the permission to share it. If you are author or own the copyright of this book, please report to us by using this DMCA report form. Report DMCA


Overview

Download & View Distilasi Dan Absorpsi Gas.pdf as PDF for free.

More details

  • Words: 5,491
  • Pages: 160
PERANCANGAN ALAT BIOPROSES Dr.rer.nat. Ir. Yuswan Muharam, M.T. Ir. Rita Arbiati, MSc.

TUJUAN β€’ Mahasiswa mampu merancang alat bioproses. β€’ Mahasiswa mampu menentukan ukuran alat bioproses.

RUJUKAN

DISTILASI DAN ABSORPSI GAS

DISTILASI

ABSORPSI

DISTILASI DAN ABSORPSI GAS

Heat

Liquid solvent

DISTILASI DAN ABSORPSI GAS Distilasi

Absorpsi

Pemisahan gas(uap)-cair

Pemisahan gas(uap)-cair

Komposisi fasa gas dan fasa cair pada kesetimbangan berbeda

Komposisi fasa gas dan fasa cair pada kesetimbangan berbeda

Fasa gas kaya akan komponen yang mudah menguap.

Fasa gas kaya akan komponen yang mudah menguap.

Menggunakan panas untuk penguapan dan pendingin untuk kondensasi.

Menggunakan cairan yang tidak mudah menguap dan mampu melarutkan komponen di fasa uap.

Perpindahan komponen dari gas ke cair dan dari cair ke gas

Perpindahan komponen dari gas ke cair

Alat dan persamaan desain sama dengan absorbsi.

Alat dan persamaan desain sama dengan distilasi.

Dilengkapi dengan reboiler dan kondenser.

Dilengkapi dengan kolom stripping dan reboiler.

DISTILASI DAN ABSORPSI GAS

DISTILASI DAN ABSORPSI GAS

DISTILASI DAN ABSORPSI GAS

DISTILASI DAN ABSORPSI GAS

Kesetimbangan Uap-Cair β€’ Hubungan antara komposisi uap dan cairan pada berbagai suhu dan tekanan. 𝑦𝑖 = 𝑓 𝑇, 𝑃, π‘₯1 , π‘₯2 , … , π‘₯𝑛 𝑦𝑖 = 𝐾𝑖 π‘₯𝑖 Ki atau vaporization equilibrium ratio (VER) = f(P, T, xi, yi) β€’ Hukum Raoult:

𝑦𝑖 = 𝑃𝑖sat ΀𝑃 π‘₯𝑖 𝐾𝑖

β€’ Pada P < 5 atm:

ideal

= 𝑃𝑖sat ΀𝑃

Volatilitas Relatif β€’ Volatilitas relatif: komposisi fasa uap dan cairan dari sistem biner pada kesetimbangan

(tidak konstan)

DISTILASI KOMPONEN JAMAK Dr.rer.nat. Ir. Yuswan Muharam, M.T. Ir. Rita Arbiati, MSc.

DISTILASI KOMPONEN JAMAK β€’ Kolom distilasi: pemisahan di antara dua produk (atas dan bawah), atau komponen-komponen murni antara komponen kunci ringan dan berat. βœ“ Bagian rectifying (di atas umpan) βœ“ Bagian stripping (di bawah umpan)

β€’ Komponen kunci ringan: komponen yang paling mudah menguap yang konsentrasinya dikendalikan di produk bawah (bottom). β€’ Komponen kunci berat: komponen yang paling sedikit menguap yang konsentrasinya dikendalikan di produk atas (overhead). β€’ Komponen terdistribusi: komponen yang volatilitasnya menengah yang terdistribusi di produk atas dan produk bawah. β€’ Jika produk > 2, misalnya n produk atas dan bawah, jumlah kolom = n - 1.

DISTILASI KOMPONEN JAMAK β€’ Produk samping diambil dari kolom dengan komposisi komponen antara. βœ“ Kolom fraksionasi minyak bumi: aliran campuran komponen dengan rentang titik didih tertentu. βœ“ Aliran dengan sedikit pengotor namun kritis yang konsentrasinya tinggi pada lokasi tertentu karena terjadi pembalikan volatilitas akibat gradien konsentrasi. ➒ 1-pentuna dengan kehadiran n-pentana di dalam campuran rengkahan C yang kaya akan isoprene. ➒ Rektifikasi alkohol fermentasi, ester dan alkohol tinggi yang tidak diinginkan menumpuk pada posisi tertentu karena kelarutannya di dalam etanol berkonsentrasi tinggi dan rendah dalam air sangat berbeda.

Normal

Anomali

RANGKAIAN KOLOM β€’ n produk atas dan bawah dengan n - 1 kolom dapat disusun dengan beberapa cara. β€’ Cara langsung: βœ“ Komponen-komponen yang paling mudah menguap dipisahkan satu per satu sebagai produk atas kolom-kolom yang disusun seri dengan produk terberat sebagai produk bawah kolom terakhir.

β€’ Jumlah cara naik tajam dengan jumlah produk βœ“ 3 produk β†’2 susunan. βœ“ 7 produk β†’ > 100 susunan.

RANGKAIAN KOLOM β€’ Panduan: β€’ Lakukan pemisahan termudah terlebih dahulu, yaitu pemisahan yang paling sedikit memerlukan talam dan refluks. Pemisahan paling sulit dilakukan terakhir. β€’ Apabila volatilitas relatif atau konsentrasi pada umpan tidak sangat variatif, singkirkan komponen-komponen itu satu per satu sebagai produk atas. β€’ Apabila volatilitas relatif pada umpan sangat variatif, susun pemisahan dalam urutan volatilitas relatif yang menurun. β€’ Apabila konsentrasi umpan sangat variatif tetapi volatilitas relatifnya tidak, urutkan pemisahan komponen-komponen umpan dalam urutan konsentrasi menurun.

JUMLAH VARIABEL BEBAS β€’ Kinerja kolom distilasi diprediksi dengan menyelesaikan persamaan matematis: βœ“ βœ“ βœ“ βœ“

Neraca massa Neraca panas Kesetimbangan uap-cair Batasan fraksi mol (spesifikasi)

di setiap talam

β€’ Pemisahan C komponen dengan n talam: βœ“ Variabel-variable pada persamaan, βœ“ C + 6 variabel.

atau komposisi produk atas atau komposisi produk bawah

TAHAP DESAIN KOLOM DISTILASI β€’ Tentukan spesifikasi distribusi komponen ringan dan berat yang diinginkan. β€’ Tentukan kondisi operasi spesifik. β€’ Hitung ukuran alat (sizing)

β€’ Hitung keekonomian atau tentukan berdasarkan pengalaman.

Parameter Desain β€’ Jumlah minimum talam teoritis, β€’ Distribusi komponen bukan kunci di produk atas dan bawah, β€’ Rasio refluks minimum, β€’ Rasio refluks operasi, β€’ Jumlah talam teoritis, β€’ Lokasi talam umpan, β€’ Efisiensi talam.

Metode Desain β€’ Metode pintas dan cepat (short cut). β€’ Metode teliti (rigorous). β€’ Menggunakan program komputer. β€’ Menggunakan nilai-nilai parameter hasil perhitungan metode pintas sebagai nilai taksiran awal (dasar). β€’ Proses iteratif yang bergantung pada konvergensi nilai taksiran awal.

Kolom Paking β€’ Kondisi berubah secara kontinyu dari atas ke bawah. β€’ Makroskopis (persamaan desain): ketinggian paking setara dengan talam teoretis (HETP) β†’ menggunakan teori kolom talam. β€’ Mikroskopis (persamaan perubahan): variasi kondisi secara kontinyu.

Metode Pintas (Metode FenskeUnderwood-Gilliland) β€’ Asumsi luapan (overflow) molal konstan di zona rectifying dan stripping β€’ Asumsi volatilitas relatif konstan β€’ Pada kondisi talam umpan. β€’ Pada kondisi rata-rata geometrik talam paling atas dan paling bawah. Karena kondisi pada talam paling atas dan paling bawah belum diketahui, maka evaluasi volatilitas relatif rata-rata adalah proses iteratif yang bisa dimulai dengan kondisi pada talam umpan.

β€’ Metode Winn: rasio ekuilibrium penguapan (VER)

Skala industri:  ο‚³ 1,05 β€’  dan  dari plot log–log untuk K β€’ atau

 = 1 jika tidak ada informasi

Jumlah Talam Minimum (Refluks Total) β€’ Underwood-Fenske:

β€’ VER:

Distribusi Komponen Bukan Kunci

Rasio Refluks Minimum β€’ Metode Underwood ( = parameter bantu)

β€’ Jika hanya satu komponen terdistribusi (dk):

β€’ q = kondisi termal umpan

Rasio Refluks Minimum β€’ Underwood:

β€’ Jika tidak ada komponen terdistribusi:

Rasio Refluks Operasi β€’ Rasio refluks operasi: βœ“ Lebih besar dari rasio refluks minimum. βœ“ Penetapannya harus mempertimbangkan keekonomian (biaya operasi dan biaya capital). βœ“ R = 1,2Rm (mendekati optimal).

Jumlah Talam Teoritis β€’ Jumlah talam teoretis: jumlah talam dimana kesetimbangan uap-cair dicapai untuk semua komponen. β€’ Gilliland (untuk q rendah):

Jumlah Talam Teoritis β€’ Molokanov:

β€’ Rusche:

Lokasi Talam Umpan β€’ Underwood-Fenske termodifikasi:

β€’ Akashah, Erbar dan Maddox:

βœ“ Kirkbride:

Contoh Soal: β€’ Suatu campuran memiliki kondisi termal q = 0,8 dan tekanan 10 atm. Campuran ini akan difraksionasi sehingga 98% komponen C dan 1% komponen E muncul di produk atas. Hitung rasio refluks, jumlah talam teoritis, dan posisi talam umpan!

Jawab: β€’ Jumlah talam minimum:

Jawab: β€’ Distribusi komponen A:

β€’ Distribusi komponen lainnya dihitung dengan cara yang sama.

Jawab: β€’ Distribusi komponen

Jawab: β€’ Refluks minimum:

Jawab: β€’ Refluks minimum:

Jawab: β€’ Refluks minimum: Substitusi 1 = 1,8817, 2 = 1,12403

Jawab: β€’ Refluks operasi dan jumlah talam: Apabila R = 1,2 Rm = 1,2(2,3077) = 2,7692

Jawab: β€’ Lokasi talam umpan:

Jawab: β€’ Lokasi talam umpan:

Lebih akurat!

Jawab: β€’ Distribusi komponen D:

 A B C D E F

f 3.1 2.6 2.2 1.3 1 0.8

0.03 0.07 0.15 0.33 0.3 0.12

d 0.0300 0.0698 0.1470 0.0930 0.0030 0.0001

b 0.0000 0.0002 0.0030 0.2370 0.2970 0.1199

Efisiensi Talam β€’ Perhitungan sejauh ini adalah untuk jumlah talam teoretis. β€’ Efisiensi talam aktual tergantung pada: β€’ Jenis sistem, β€’ Laju alir, β€’ Konstruksi talam. β€’ Rentang efisiensi talam: β€’ < 10% hingga > 100%, β€’ Ketidakpastian terbesar dalam perancangan kolom distilasi. β€’ 60% (fraksionasi hidrokarbon)

PEMISAHAN DI DALAM KOLOM PAKING

PEMISAHAN DI DALAM KOLOM PAKING β€’ Kontaktor: paking β€’ Pemisahan kontinyu β€’ Teori perpindahan massa antarfasa β€’ Fluks perpindahan massa antarfasa: 𝑁 = π‘˜G 𝑦 βˆ’ 𝑦 βˆ— = π‘˜L π‘₯ βˆ— βˆ’ π‘₯ 𝐴

β€’ Dua kasus: βœ“ Difusi dua arah sama-molar (distilasi) βœ“ Difusi searah (absorpsi dan stipping)

π‘˜G

effective

=

π‘˜G 𝑦 βˆ’ π‘¦βˆ—

log mean

βˆ— 1 βˆ’ 𝑦 βˆ’ 1 βˆ’ 𝑦 𝑦 βˆ’ π‘¦βˆ— log mean = ln 1 βˆ’ 𝑦 Ξ€ 1 βˆ’ π‘¦βˆ—

π‘¦βˆ— βˆ’ 𝑦 = ln 1 βˆ’ 𝑦 Ξ€ 1 βˆ’ π‘¦βˆ—

Koefisien Perpindahan Massa β€’ Koefisien perpindahan massa total (partikel bola):

Distilasi

Distilasi Satu Tahanan Perpindahan Massa: β€’ Laju pada elemen dengan tinggi dz: 𝑑𝑁 = π‘˜Gπ‘Ž 𝑦 βˆ’ 𝑦 βˆ— 𝑑𝑧 𝑑 𝐺𝑦 βˆ’ = π‘˜G π‘Ž 𝑦 βˆ’ 𝑦 βˆ— 𝑑𝑧 𝐺𝑑𝑦 βˆ’ = π‘˜Gπ‘Ž 𝑦 βˆ’ π‘¦βˆ— 𝑑𝑧 𝐺 𝑑𝑦 𝑑𝑧 = βˆ’ π‘˜Gπ‘Ž 𝑦 βˆ’ 𝑦 βˆ— 𝑍 𝑦2 𝐺 𝑑𝑦 ΰΆ± 𝑑𝑧 = βˆ’ ΰΆ± π‘˜Gπ‘Ž 𝑦1 𝑦 βˆ’ 𝑦 βˆ— 0 𝑦2 𝐺 𝑑𝑦 Tinggi kolom: 𝑍= ΰΆ± π‘˜G π‘Ž 𝑦1 𝑦 βˆ— βˆ’ 𝑦 βˆ’

a = luas permukaan antarmuka per satuan volume unggun.

Distilasi Satu Tahanan Perpindahan Massa: β€’ Laju pada elemen dengan tinggi dz: 𝑑𝑁 = π‘˜L π‘Ž π‘₯ βˆ— βˆ’ π‘₯ 𝑑𝑧 𝑑 𝐿π‘₯ βˆ’ = π‘˜L π‘Ž π‘₯ βˆ— βˆ’ π‘₯ 𝑑𝑧 𝐿𝑑π‘₯ βˆ’ = π‘˜L π‘Ž π‘₯ βˆ— βˆ’ π‘₯ 𝑑𝑧 𝐿 𝑑π‘₯ 𝑑𝑧 = βˆ’ π‘˜L π‘Ž π‘₯ βˆ— βˆ’ π‘₯ 𝑍 π‘₯2 𝐿 𝑑π‘₯ ΰΆ± 𝑑𝑧 = βˆ’ ΰΆ± π‘˜L π‘Ž π‘₯1 π‘₯ βˆ— βˆ’ π‘₯ 0 π‘₯2 𝐿 𝑑π‘₯ 𝑍= ΰΆ± π‘˜L π‘Ž π‘₯1 π‘₯ βˆ’ π‘₯ βˆ— βˆ’

Tinggi kolom:

Distilasi β€’ Hubungan antara (x*, y*) dan (x, y): 𝑑𝑁 βˆ’ = π‘˜Gπ‘Ž 𝑦 βˆ’ π‘¦βˆ— 𝑑𝑧

𝑦 βˆ’ 𝑦 βˆ— π‘˜L = π‘₯ βˆ— βˆ’ π‘₯ π‘˜G

𝑑𝑁 βˆ’ = π‘˜L π‘Ž π‘₯ βˆ— βˆ’ π‘₯ 𝑑𝑧

π‘¦βˆ— βˆ’ 𝑦 π‘˜L = βˆ’ π‘₯βˆ— βˆ’ π‘₯ π‘˜G

(x, y) β†’ neraca massa. (x*, y*) β†’ kesetimbangan (grafik, tabel, persamaan). π‘¦βˆ— = 𝑓 π‘₯ βˆ— βˆ— 𝛼π‘₯ π‘¦βˆ— = 1 + 𝛼 βˆ’ 1 π‘₯βˆ—

Distilasi β€’ Jumlah tray teoritis β†’ jumlah unit perpindahan (distilasi biner):

β€’ Tinggi unit perpindahan: 𝐺 π»π‘‡π‘ˆG = π‘˜Gπ‘Ž

𝐿 π»π‘‡π‘ˆL = π‘˜L π‘Ž

β€’ Tinggi unggun paking: 𝑍 = π»π‘‡π‘ˆGπ‘π‘‡π‘ˆG = π»π‘‡π‘ˆL π‘π‘‡π‘ˆL

Distilasi

Distilasi Dua tahanan perpindahan massa β†’ HTU overall: π»π‘‡π‘ˆoG = π»π‘‡π‘ˆG + πœ†π»π‘‡π‘ˆL Faktor stripping: πœ† =

𝐾𝑖 𝐺 𝐿

Koefisien perpindahan massa overall:

1 1 𝐻𝑖 = + π‘˜π‘œG π‘˜G π‘˜L

Konstanta hokum Henry: 𝐻𝑖 = 𝐾𝑖 𝑃

𝑍 = π»π‘‡π‘ˆoGπ‘π‘‡π‘ˆπ‘œG

Distilasi Campuran Multikomponen β€’ Height equivalent to a theoretical plate (HETP) komponen kunci: ln πœ† 𝐻𝐸𝑇𝑃 = π»π‘‡π‘ˆπ‘œG πœ†βˆ’1 π‘π‘‡π‘ˆπ‘œG ln πœ† πœ†βˆ’1

β€’ Jumlah tahap teoritis:

𝑁t =

β€’ Tinggi kolom:

𝑍 = 𝐻𝐸𝑇𝑃 𝑁t

β€’ Jika kL/kG >>> β†’ HETP = HTUOG

Contoh Soal β€’ Sebuah campuran equimolal pada titik didihnya dipisahkan menjadi 95% dan 5% komponen ringan di produk atas dan produk bawah. Volatilitas relatif  = 2, refluks minimum = 1,714, dan refluks operasi = 50% lebih besar dari refluks minimum. Koefisien perpindahan massa pada film gas kG = 0,012 kmol/m2.jam, sedangkan tahanan perpindahan massa pada sisi cair sangat kecil. Paking yang digunakan adalah berl saddle. Jika laju alir umpan = 1 kmol/m2.jam, hitung tinggi kolom paking dan tentukan posisi umpan!

x2 = 0,95

xF = 0,5

x1 = 0,05

Jawab β€’ Diketahui: F = 1 kmol/m2.jam xF = 0.5 q=1 x2 = 0.95 x1 = 0.05 =2 Rm = 1,714 R = 1,5Rm Berl saddle: kG = 0,012 kmol/m2.jam kL >>>> (tidak digunakan) a = 250 m2/m3

x2 = 0,95

xF= 0,5

x1 = 0,05

Jawab β€’ Hitung fluks produk atas, produk bawah, uap di bagian rectifying dan uap di bagian stripping: β–ͺ Gunakan rumus untuk distilasi biner dengan kontaktor talam! D = 0,5 kmol/m2.jam B = 0,5 kmol/m2.jam Vr = 1,786 kmol/m2.jam Vs = 1,786 kmol/m2.jam β€’ Hitung rasio refluks operasi: R = 1,5Rm = 2,571

Jawab 1

β€’ Buat kurva kesetimbangan!

0.9 0.8

𝐾𝑖 𝛼= =2 𝐾𝑗

0.7

0.6 0.5 0.4

0.3 0.2

0.1 0 0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

Jawab 1

β€’ Buat garis operasi! 𝑦r = 𝑦s =

𝑅

𝑅+1

π‘₯r +

𝑅𝐷+π‘žπΉ

1

𝑅+1

π‘₯ 𝑅𝐷+π‘žπΉβˆ’π΅ s

βˆ’

0.9 0.8

π‘₯D

0.7

𝐡

π‘₯ 𝑅𝐷+π‘žπΉβˆ’π΅ B

0.6 0.5 0.4

0.3 0.2

0.1 0 0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

Jawab 1

β€’ Buat garis-garis tahanan perpindahan antarfasa:

0.9

0.8

π‘¦βˆ— βˆ’ 𝑦 π‘˜L = βˆ’ =ο‚₯ βˆ— π‘₯ βˆ’π‘₯ π‘˜G

0.7

0.6 0.5 0.4

0.3 0.2

0.1 0 0

0.2

0.4

0.6

0.8

1

Jawab β€’ Hitung pasangan nilai y dan y*, tabulasikan! Kemudian hitung nilai y* 0.093 0.178 0.257 0.331 0.399 0.462 0.520 0.573 0.623 0.669 0.712 0.751 0.788 0.823 0.856 0.887 0.916 0.945 0.973

1/(y*-y) 23.424 15.621 12.623 11.281 10.795 10.924 11.668 13.251 16.347 23.163 20.077 18.719 18.397 18.881 20.170 22.445 26.160 32.313 43.371

π‘¦βˆ— βˆ’π‘¦

50 45 40 35

1/(y*-y)

y 0.050 0.114 0.178 0.242 0.306 0.370 0.434 0.498 0.562 0.626 0.662 0.698 0.734 0.770 0.806 0.842 0.878 0.914 0.950

1

30 25 20 15 10 5 0

0.0

0.1

0.2

0.3

0.4

0.5

y

0.6

0.7

0.8

0.9

1.0

Jawab β€’ Integrasikan! y* 0.093 0.178 0.257 0.331 0.399 0.462 0.520 0.573 0.623 0.669 0.712 0.751 0.788 0.823 0.856 0.887 0.916 0.945 0.973

1/(y*-y) 23.424 15.621 12.623 11.281 10.795 10.924 11.668 13.251 16.347 23.163 20.077 18.719 18.397 18.881 20.170 22.445 26.160 32.313 43.371

50 45 40 35

1/(y*-y)

y 0.050 0.114 0.178 0.242 0.306 0.370 0.434 0.498 0.562 0.626 0.662 0.698 0.734 0.770 0.806 0.842 0.878 0.914 0.950

30 25 20 15 10 5 0

0.0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1.0

y

Luas area di bawah kurva β†’ integrasi numeris!

Jawab β€’ Hasil integrasi numeris: βœ“ Luas area di bawah kurva = 15,627 (NTUG) β–ͺ Luas zona stripping = 8,052 (NTUGs) β–ͺ Luas zona rectifying = 7,575 (NTUGr) β€’ Hitung HTUG βœ“ π»π‘‡π‘ˆG =

𝐺 π‘˜G π‘Ž

β–ͺ Zona stripping: π»π‘‡π‘ˆGs = β–ͺ Zona rectifying: π»π‘‡π‘ˆGr =

𝑉s π‘˜G π‘Ž 𝑉r

π‘˜G π‘Ž

= 0,595 = 0,595

Jawab β€’ Hitung tinggi unggung paking: βœ“ Tinggi unggun paking zona stripping: β–ͺ 𝑍s = π»π‘‡π‘ˆGsπ‘π‘‡π‘ˆGs = 4,792 m βœ“ Tinggi unggun paking zona rectifying: β–ͺ 𝑍r = π»π‘‡π‘ˆGrπ‘π‘‡π‘ˆGr = 4,508 m β€’ Maka: β€’ Tinggi unggun paking, Z = 9,3 m β€’ Posisi umpan = 4,508 m dari atas

Absorpsi dan Stripping β€’ Absorpsi: βœ“ Gas β†’ cair β€’ Stripping βœ“ Cair β†’ gas β€’ Laju alir massa atau molar tidak konstan. β€’ Untuk 1 komponen berpindah β†’ laju alir fasa bebas zat terlarut, G’ atau L’ G’ = G(1 – y) = G1(1 - y1) L’ = L(1 – x) = L1(1 - x1)

Absorpsi dan Stripping β€’ Neraca massa di sekitar bagian bawah kolom: 𝐺𝑦 + 𝐿1 π‘₯1 = 𝐺1 𝑦1 + 𝐿π‘₯ Eliminasi G dan L: 𝐺′ 𝐿′ 𝐺′ 𝐿′ 𝑦+ π‘₯ = 𝑦 + π‘₯ 1βˆ’π‘¦ 1 βˆ’ π‘₯1 1 1 βˆ’ 𝑦1 1 1 βˆ’ π‘₯ 𝑦 𝐿′ π‘₯ 𝑦1 𝐿′ π‘₯1 = + βˆ’ 1 βˆ’ 𝑦 𝐺′ 1 βˆ’ π‘₯ 1 βˆ’ 𝑦1 𝐺′ 1 βˆ’ π‘₯1

Absorpsi dan Stripping β€’ Neraca massa di sekitar bagian atas kolom: 𝐺𝑦 + 𝐿2 π‘₯2 = 𝐺2 𝑦2 + 𝐿π‘₯ Eliminasi G dan L: 𝐺′ 𝐿′ 𝐺′ 𝐿′ 𝑦+ π‘₯ = 𝑦 + π‘₯ 1βˆ’π‘¦ 1 βˆ’ π‘₯2 2 1 βˆ’ 𝑦2 2 1 βˆ’ π‘₯ 𝑦 𝐿′ π‘₯ 𝑦2 𝐿′ π‘₯2 = + βˆ’ 1 βˆ’ 𝑦 𝐺′ 1 βˆ’ π‘₯ 1 βˆ’ 𝑦2 𝐺′ 1 βˆ’ π‘₯2

Absorpsi dan Stripping β€’ Neraca massa pada elemen dengan tinggi dz: 𝑑𝑁 = 𝑑 𝐺𝑦 = π‘˜G

eff π‘Ž

𝑑𝑁 = 𝑑 𝐿π‘₯ = π‘˜L

eff π‘Ž

𝑦 βˆ’ 𝑦 βˆ— 𝑑𝑧

π‘₯ βˆ— βˆ’ π‘₯ 𝑑𝑧

Absorpsi dan Stripping β€’ Neraca massa di sekitar bagian bawah kolom: βœ“ Fasa gas: 𝐺′ 𝑦 𝑦 β€² β€² 𝑑 𝐺𝑦 = 𝑑 𝑦 = 𝑑𝐺 𝑑𝑦 + 𝐺 𝑑 𝑑𝑦 1βˆ’π‘¦ 1βˆ’π‘¦ 1βˆ’π‘¦ β€² β€² 𝐺 𝐺 𝐺 =0+ 𝑑𝑦 = 𝑑𝑦 = 𝑑𝑦 1βˆ’π‘¦ 2 1βˆ’π‘¦ 2 1βˆ’π‘¦ βœ“ Fasa cair: 𝐿′ π‘₯ β€² π‘₯ β€² 𝑑 𝐿π‘₯ = 𝑑 π‘₯ = 𝑑𝐿 𝑑π‘₯ + 𝐿 𝑑 𝑑π‘₯ 1βˆ’π‘₯ 1βˆ’π‘₯ 1βˆ’π‘₯ 𝐿′ 𝐿′ 𝐿 =0+ 𝑑π‘₯ = 𝑑π‘₯ = 𝑑π‘₯ 1βˆ’π‘₯ 2 1βˆ’π‘₯ 2 1βˆ’π‘₯

Absorpsi dan Stripping β€’ Substitusi 𝑑 𝐺𝑦 =

𝐺

1βˆ’π‘¦

𝑑𝑦 dan π‘˜G

𝑑 𝐺𝑦 = π‘˜G

effective

mean π‘Ž

=

π‘˜G 𝑦 βˆ’ π‘¦βˆ—

log mean

𝑦 βˆ’ 𝑦 βˆ— 𝑑𝑧

Menghasilkan 𝐺 𝑍= π‘˜G π‘Ž atau

𝐺 𝑍= π‘˜G π‘Ž

𝑦2

𝑦 βˆ’ 𝑦 βˆ— log mean ΰΆ± 𝑑𝑦 1 βˆ’ 𝑦 𝑦 βˆ’ π‘¦βˆ—

mean 𝑦

1

𝑦2

𝑑𝑦 ΰΆ± 1 βˆ’ 𝑦 ln 1 βˆ’ 𝑦 Ξ€ 1 βˆ’ 𝑦 βˆ—

mean 𝑦1

ke

Absorpsi dan Stripping β€’ Tinggi unit perpindahan: 𝐺 π»π‘‡π‘ˆG = π‘˜G π‘Ž

π»π‘‡π‘ˆL =

𝐿 π‘˜L π‘Ž

mean

mean

Absorpsi dan Stripping β€’ Jumlah unit perpindahan: 𝑦2

𝑑𝑦 π‘π‘‡π‘ˆG = ΰΆ± 1 βˆ’ 𝑦 ln 1 βˆ’ 𝑦 Ξ€ 1 βˆ’ 𝑦 βˆ— 𝑦1

π‘₯2

𝑑π‘₯ π‘π‘‡π‘ˆL = ΰΆ± 1 βˆ’ π‘₯ ln 1 βˆ’ π‘₯ Ξ€ 1 βˆ’ π‘₯ βˆ— π‘₯1

Absorpsi dan Stripping β€’ Jumlah unit perpindahan keseluruhan: 1 1 πœ† = + π‘π‘‡π‘ˆπ‘œG π‘π‘‡π‘ˆG π‘π‘‡π‘ˆL

Absorpsi dan Stripping β€’ Garis operasi dan kurva kesetimbangan lurus dan umpan pelarut bebas dari solut: 1 π‘π‘‡π‘ˆG = ln π‘šπΊ 1βˆ’ 𝐿

π‘šπΊ 𝑦1 π‘šπΊ 1βˆ’ + 𝐿 𝑦2 𝐿

m = kemiringan kurva (garis kesetimbangan)

Absorpsi dan Stripping β€’ Dua tahanan lapisan film: 1 𝐾oG π‘Ž

=

1 π‘˜G π‘Ž

+

𝐻𝑖 π‘˜L π‘Ž

Hi = konstantan hukum Henry 𝑍=

atau

𝐺

𝐾oG π‘Ž

𝑦2

𝑑𝑦 ΰΆ± 1 βˆ’ 𝑦 ln 1 βˆ’ 𝑦 Ξ€ 1 βˆ’ 𝑦 βˆ—

mean 𝑦1

𝑍 = π»π‘‡π‘ˆoG π‘π‘‡π‘ˆoG

Absorpsi dan Stripping

Contoh Soal β€’ Karbondioksida (CO2) di dalam gas bumi dipisahkan di acid removal gas unit (ARGU) dengan menggunakan kolom absorpsi jenis paking. CO2 ini akan digunakan untuk pertumbuhan mikroalga di industri biofuel. Pelarut yang digunakan adalah dietanolamina (DEA) 2 M. Kadar CO2 di dalam gas umpan adalah 50%. Kadar CO2 di dalam pelarut lean DEA adalah 10%, sedangkan kadarnya di dalam pelarut rich DEA yang keluar kolom adalah 35%. Apabila laju alir gas dan cairan yang bebas CO2 adalah 2 dan 4 kmol/m2.jam, hitung tinggi kolom!

Jawab β€’ Diketahui: x2 = 0,1 (kadar CO2 di dalam larutan lean DEA ) x1 = 0,35 (kadar CO2 di dalam larutan rich DEA) y1 = 0,5 (kadar CO2 di dalam larutan umpan gas bumi) P = 2 atm (tekanan operasi) T = 323 K (temperatur operasi) L’ = 4 kmol/m2.jam (fluks larutan DEA bebas CO2) G’ = 2 kmol/m2.jam (fluks gas bumi bebas CO2) kG = 0,025 kmol/m2.jam.atm (koef perpindahan massa di film gas) -kL/kG = βˆ’0,94199 a = 250 m2/m3 (luas area spesifik) HCO2 = 0.50956 atm.m3/kmol (konstanta Henry)

Jawab β€’ Buat kurva kesetimbangan: 0.6

0.5

y

0.4

0.3

0.2

0.1

0

0

0.1

0.2

0.3 x

0.4

0.5

Jawab β€’ Buat garis operasi (neraca massa): 0.6

0.5

0.4

y

𝑦 1βˆ’π‘¦ 𝐿′ π‘₯ = 𝐺′ 1 βˆ’ π‘₯ 𝑦1 𝐿′ π‘₯1 + βˆ’ 1 βˆ’ 𝑦1 𝐺′ 1 βˆ’ π‘₯1

0.3

0.2

0.1

0

0

0.1

0.2

0.3 x

0.4

0.5

Jawab β€’ Buat garis-garis tahanan perpindahan antarfasa: π‘˜L βˆ’ = βˆ’0,94199 π‘˜G

0.6

0.5

y

0.4

0.3

0.2

0.1

0 0

0.1

0.2

0.3

x

0.4

0.5

Jawab β€’ Hitung jumlah satuan perpindahan: 0.500 0.470 0.438 0.405 0.371 0.335 0.297 0.258 0.216 0.173 0.127 Luas =

Integran

0

0.0 0.53002 0.42646 0.36638 0.32838 0.30368 0.28832 0.28103 0.28343 0.30463 0.37180 3.484

0.1

0.2

0.3

-5 1/(1βˆ’π‘¦)ln[(1βˆ’π‘¦)βˆ•(1βˆ’π‘¦βˆ— ) ]

y

-10

-15

-20

NTUoG = 3,484

-25

y

0.4

0.5

0.6

Jawab β€’ Hitung tinggi satuan perpindahan: 1 𝐾oG π‘Ž

=

1 π‘˜G π‘Ž

+

𝐻𝑖 π‘˜L π‘Ž

π»π‘‡π‘ˆoG =

= 0,330 m3 jam. atm/kmol 𝐺

𝐾oG π‘Ž mean

= 2,020 m

β€’ Hitung tinggi kolom: 𝑍 = π»π‘‡π‘ˆoGπ‘π‘‡π‘ˆoG = 7,038 m

Tugas Kelompok β€’ Senyawa A di dalam gas dipisahkan dengan menggunakan kolom absorpsi jenis paking dengan pelarut B. Kadar A di dalam gas umpan adalah 75%. Kadar A di dalam pelarut lean adalah 2,5%, sedangkan kadarnya di dalam pelarut rich adalah 40%. Apabila fluks aliran gas dan cair yang bebas A adalah 4 dan 10 kmol/m2.jam, hitung tinggi kolom! Data lain: a = 250/m kG = 0,015 kmol/m2.hr.atm kL = 0,012 m/hr H = 0.5 atm.m3/kmol

 = 2,5 π‘¦βˆ—

𝛼π‘₯ βˆ— = 1 + 𝛼 βˆ’ 1 π‘₯βˆ—

Distilasi Kolom Paking Komponen Jamak β€’ Tinggi kolom: 𝑍 = 𝐻𝐸𝑇𝑃 𝑁t Nt = jumlah talam teoritis HETP = height equivalent to a theoretical plate (untuk komponen kunci) 𝐻𝐸𝑇𝑃 =

lnπœ† π»π‘‡π‘ˆoG πœ† βˆ’1

 = rasio kemiringan garus kesetimbangan terhadap kemiringan garis operasi = m/(L/G) m sangat bervariasi apabila volatilitas relatif tinggi.

ABSORPSI FISIK 𝑑𝐢𝐢CO ,G 𝑑 2 𝐢CO2,G 2 βˆ— 𝑠G πœ€P 𝐷𝑧,G βˆ’π‘ˆ βˆ’π‘˜ π‘Ž 𝐢 βˆ’ 𝐢CO2,L = 0 G L,CO CO 2 2 ,L 𝑑𝑧 2 𝑑𝑧 𝑑 2 𝐢CO2,L 𝑑𝐢CO2 ,L βˆ— 𝑠L πœ€P 𝐷𝑧,L βˆ’π‘ˆ +π‘˜ π‘Ž 𝐢 βˆ’ 𝐢CO2,L = 0 L L,CO CO 2 2 ,L 𝑑𝑧 2 𝑑𝑧

ABSORPSI KIMIA 𝑑𝐢𝐢CO ,G 𝑑 2 𝐢CO2,G 2 βˆ— 𝑠G πœ€G 𝐷𝑧,G βˆ’π‘ˆ βˆ’π‘˜ π‘Ž 𝐢 βˆ’ 𝐢CO2,L = 0 G L,CO CO 2 2 ,L 𝑑𝑧 2 𝑑𝑧 𝑑 2 𝐢CO2,L 𝑑𝐢CO2,L βˆ— 𝑠L πœ€L 𝐷𝑧,L βˆ’π‘ˆ +π‘˜ π‘Ž 𝐢 βˆ’ 𝐢CO2,L βˆ’ 𝑅CO2 ,L = 0 L L,CO CO 2 2 ,L 𝑑𝑧 2 𝑑𝑧 𝑑 2 𝐢MEA,L 𝑑𝐢MEA,L 𝑠L πœ€L 𝐷𝑧,L βˆ’π‘ˆ + 𝑅MEA,L = 0 L 2 𝑑𝑧 𝑑𝑧

ABSORPSI KIMIA β€’ Reaksi kimia CO2-MEA: k

CO2 + 2RNH2 ⇋kβˆ’1 RNCOOβˆ’ + RH2+ 2

𝑅 = π‘˜2 CO2 MEA π‘™π‘œπ‘” π‘˜2 = 10.99 βˆ’

2152 𝑇

Kondisi Batas β€’ Fasa gas: βœ“ Pada z = 0: CCO2,G = CCO2,G,inlet βœ“ Pada z = L:

𝑑𝐢CO2,G =0 𝑑𝑧

β€’ Fasa cair: βœ“ Pada z = L β†’ CCO2,L = CCO2,L,inlet CMEA,L = CMEA,L,inlet βœ“ Pada z = 0:

𝑑𝐢CO2,L =0 𝑑𝑧 𝑑𝐢MEA, L =0 𝑑𝑧

CONTOH β€’ CO2 di dalam aliran gas bumi ingin diturunkan kadarnya dari 20% hingga 5% dengan menggunakan pelarut kimia MEA. Apabila tinggi kolom 15 m dan diameter kolom 1,5 m, apakah desain dengan kondisi di bawah sudah sesuai? Data: mol CO2/mol MEA di dalam lean solven = 0,1 Konsentrasi MEA inlet = 1000 mol/m3 Kecepatan superfisial gas = 1 m/detik Kecepatan superfisial cair = 1 m/detik Tekanan gas = 10 bar Suhu = 323 K Paking berupa silinder berongga dengan jari-jari 1 cm.

JAWAB

JAWAB

Pengaruh UG

Pengaruh UG

Pengaruh UG

Pengaruh UL

Pengaruh UL

Pengaruh P

Pengaruh P

Pengaruh T

Pengaruh T

Pengaruh CMEA,inlet

Pengaruh CMEA,inlet

Pengaruh mol CO2/mol MEA lean

Pengaruh mol CO2/mol MEA lean

Pengaruh jari-jari paking

Pengaruh jari-jari paking

Pengaruh diameter kolom

Pengaruh tinggi kolom

5m

10 m

15 m

20 m

Pengaruh diameter kolom

0,5 m

1m

1,5 m

2m

Optimalisasi

rp = 100 mm

UL = 5 m/detik

rp = 20 mm, UL = 3 m/detik

Desain β€’ CO2 di inlet gas bumi = 20% β€’ CO2 di outlet gas bumi = 5% β€’ Tinggi kolom =15 m β€’ Diameter kolom = 1,5 m β€’ Mol CO2/mol MEA di dalam lean solven = 0,1 β€’ Konsentrasi MEA inlet = 1000 mol/m 3 β€’ Kecepatan superfisial gas = 1 m/detik β€’ Kecepatan superfisial cair = 3 m/detik β€’ Tekanan gas = 10 bar β€’ Suhu = 323 K β€’ Paking berupa silinder berongga dengan jari-jari 2 cm.

Sizing Kolom Talam

Sizing Kolom Talam β€’ Uap harus mengalir hanya melalui daerah terbuka pada talam di antara downcomer-downcomer. β€’ Cairan harus mengalir hanya melalui downcomer-downcomer. β€’ Cairan tidak boleh menetes melalui lubang-lubang talam. β€’ Cairan tidak boleh terbawa naik ikut aliran uap (entrainment flooding). β€’ Uap tidak boleh turun ikut aliran cair (downcomer flooding). β€’ Uap tidak boleh menggelembung melalui downcomer.

Aliran Cairan dan Aliran Uap β€’ VG <<< cairan menetes melalui lubang. β€’ VG >>> entrainment flooding β†’ efisiensi ο‚―

Aliran Cairan dan Aliran Uap β€’ Beban gas minimum: β€’ Pada kecepatan gas rendah, gas tidak lagi mengalir secara merata melalui semua lubang talam atau cairan menetes melalui lubang talam. β€’ Kedua mode operasi itu harus dihindari karena efisiensi talam berkurang akibat tingkat pencampuran kedua fase tidak cukup. β€’ Faktor utama yang mempengaruhi tetesan adalah diameter lubang (beban gas minimum meningkat dengan meningkatnya diameter lubang).

Aliran Cairan dan Aliran Uap β€’ Beban gas maksimum: β€’ Pada kecepatan gas tinggi, gas meniup cairan dari talam dalam bentuk droplet-droplet halus (entrainment, jet flood) β€’ Cairan mengalir tidak lagi berlawanan arah dengan aliran gas. β€’ Beban gas maksimum yang layak tergantung pada sistem (kadar gas dan cairan, tegangan permukaan) serta desain talam. β€’ Entrainment flooding pada talam terjadi pada jarak talam z yang sangat besar. Untuk jarak talam yang lebih kecil tinggi buih pada talam menentukan batasan yang lebih rendah.

Aliran Cairan dan Aliran Uap β€’ Beban cair minimum: β€’ Pada beban cairan yang sangat rendah, aliran cairan tidak merata di atas talam (maldistribution) sehingga efisiensi perpindahan massa berkurang. β€’ Tinggi minimum overflow weir hwo ο‚³ 5 mm β€’ Nilai ini sama dengan beban weir cair minimum

𝑄L 𝑙w

z

lw

= 2 m3/(m/jam)

Aliran Cairan dan Aliran Uap β€’ Beban cair maksimum: β€’ Aliran cairan ke bawah melalui downcomer didesak oleh gaya gravitasi sehingga beban cairan maksimum terbatas. β€’ Empat aturan empiris untuk menentukan laju alir cairan maksimum: 1. Beban weir, 𝑄L ΀𝑙w < 60 m3/(m.jam). 2. Kecepatan cairan di downcomer ο‚£ 0,1 m/s. 3. Volume downcomer memberikan waktu tinggal cairan ο€Ύ 5 detik. 4. Tinggi cairan jernih di downcomer hld ο‚£ Β½ z.

Floading β€’ Downcomer flooding terjadi apabila cairan tergenang di atas talam karena downcomer area terlalu kecil. β€’ Entrainment flooding terjadi karena terlalu banyak cairan yang terbawa oleh aliran uap β†’ masalah besar βœ“ Kecepatan desain = 80%βˆ’85% kecepatan flooding. βœ“ Kecepatan flooding, π‘ˆf = 𝐾1 K1 = faktor kapasitas

𝜌L βˆ’πœŒG 𝜌G

Floading Faktor kapasitas: 𝐿m 𝜌G 𝐹LV = 𝐺m 𝜌L Lm = laju alir massa cairan [kg/detik] Gm = laju alir massa gas [kg/detik] β€’ Persen flooding = UG/Uf

Kecepatan Penetesan (Weeping) β€’ Kecepatan uap di lubang: 𝑄G π‘ˆh = π‘Žh QG = laju alir volumetrik gas [m3/detik] ah = luas total lubang [m2] β€’ Kecepatan weeping:

𝐾2 βˆ’ 0,90 25,4 βˆ’ 𝑑h π‘ˆh = 𝜌G dh = diameter lubang [m], K2 dari chart

Efisiensi Talam dan Paking β€’ Jumlah talam teoritis (kesetimbangan) dihitung dengan menggunakan data kesetimbangan. Namun kenyataannya kesetimbangan tidak tercapai sepenuhnya pada talam. β€’ Tinggi paking dihitung dengan anggapan bahwa pencampuran komponen di fasa gas dan fasa cair berlangsung sempurna. Namun kenyataannya, untuk kolom berdiameter besar, penyebaran komponen di dalam kolom tidak berlangsung sempurna. β€’ Jumlah talam atau tinggi paking yang dihitung dengan cara yang sudah dilakukan perlu dikoreksi sehingga jumlah talam dan tinggi paking aktual lebih besar dari yang dihitung sejauh ini.

Efisiensi Talam β€’ Jumlah talam aktual: 𝑁teoritis 𝑁aktual = 𝐸oc Eoc = efisiensi kolom keseluruhan Dipengaruhi oleh: βœ“ βœ“ βœ“ βœ“

Geometri talam, Laju alir, Rasio kesetimbangan uap-cair, Sifat fisik uap dan cair.

Efisiensi Talam β€’ Kolom distilasi (metode O’Connel):

 [cP]

Efisiensi Talam β€’ Kolom absorpsi/stripping (metode O’Connel):

H [lb mol/(ft3)(atm)], P [atm], and  [cP]

Efisiensi Talam β€’ Metode rigorous:

Jatuh Tekan β€’ Jatuh tekan total pada talam: β„Žt = β„Žd + β„ŽL [mm cair] β€’ Jatuh tekan melalui lubang (sieve tray): β„Žd = 92,8 𝜌G ΀𝜌L π‘ˆh2 Uh = kecepatan di lubang [m/detik]

Jatuh Tekan β€’ Jatuh tekan melalui campuran teraerasi (holdup cairan): β„ŽL = 𝛽 β„Žw + β„Žow

 = faktor aerasi (=0,5); hw = tinggi weir β€’ Crest di atas outlet weir: β„Žow = 664 𝑄L ΀𝑙w

2Ξ€3

QL= laju alir volumetrik cairan [m3/detik] lw = panjang weir [m]

Koefisien Perpindahan Massa β€’ Koefisien perpindahan massa volumetrik fasa uap: π‘˜G π‘Ž =

316𝐷G0,5 1030𝑓 βˆ’ 867𝑓 2

1Ξ€detik

β„Ž0,5 L

β€’ Koefisien perpindahan massa volumetrik fasa cair: π‘˜L π‘Ž = 0,40𝐹va + 0,17 3,875E8𝐷L

0,5

DG, DL = koefisien difusi fasa uap dan cair [m2/detik] hL = holdup cairan [mm] f = fraksi flooding f Va = factor-F luas aktif = π‘ˆπ‘Ž 𝜌G0,5

1Ξ€detik

Waktu Tinggal β€’ Fasa uap dalam campurannya di dalam talam: 𝑑G =

volume uap di dalam busa (π‘“π‘Ÿπ‘œπ‘‘β„Ž) 𝐴a πœ€π‘f = laju alir volumetrik uap 𝑄G

β€’ Fasa cair dalam campurannya di dalam talam: volume cairan di dalam busa (π‘“π‘Ÿπ‘œπ‘‘β„Ž) 𝐴a 1 βˆ’ πœ€ 𝑍f 𝑑L = = laju alir volumetrik cairan 𝑄L Zf = tinggi busa (froth); dari chart. Aa = luas talam aktif

ο₯ = fraksi kosong (uap) di dalam busa = 1 – hL/Zf

Waktu Tinggal

Satuan Perpindahan Massa β€’ Fasa uap: π‘π‘‡π‘ˆG = π‘˜G π‘Žπ‘‘G β€’ Fasa cair:

π‘π‘‡π‘ˆL = π‘˜L π‘Žπ‘‘L

Efisiensi Lokal (Titik) 1 1 πœ† = + π‘π‘‡π‘ˆoG π‘π‘‡π‘ˆG π‘π‘‡π‘ˆL

𝐸oG = 1 βˆ’ exp βˆ’π‘π‘‡π‘ˆoG

 = faktor stripping = m(L/G) EoG = efisiensi lokal (titik), [fraksi]

Efisiensi Talam Murphree β€’ Pencampuran sempurna di talam: β€’ Kolom kecil β€’ Kolom vakum (G >>> L) 𝐸mv = 𝐸oG

β€’ Aliran sumbat melalui talam: 𝐸mv =

1 exp πœ†πΈoG βˆ’ 1 πœ†

Efisiensi Kolom Keseluruhan β€’ Refluks total: 𝐸oc = 𝐸mv β€’ Selain refluks total: 𝐸oc =

ln 1 + 𝐸mv πœ† βˆ’ 1 lnπœ†

Koreksi Entrainment β€’ Resirkulasi cairan sebagai entrainment antara talam-talam. β€’ Mempengaruhu efisiensi. β€’ Terutama terjadi pada kolom vakum. β€’ Konversi dari efisiensi kering ke efisiensi basah (dengan entrainment) 𝐸wet

𝐸dry = 1 βˆ’ 𝐸dryπœ“Ξ€ 1 βˆ’ πœ“

Efisiensi Talam β€’ Contoh: β€’ Lihat Example 13.17 pada buku Couper dkk, 2005, Chemical Process Equipment – Selection and Design, Elsevier.

Sistem Multikomponen β€’ Menggunakan komponen kunci sebagai sistem biner semu.

Diameter Kolom β€’ Dikontrol oleh kecepatan uap. β€’ Menggunakan kecepatan superfisial uap. β€’ UG = 0,5βˆ’2,5 m/detik (3 βˆ’ 6 m/s di kolom dengan bubble cup). β€’ UL = 1x10-3βˆ’15x10-3 m/detik. β€’ Diameter kolom: 𝐺max 𝑇 𝐷c = 102,89 [m] π‘ƒπ‘ˆG,allow 𝐺max = laju alir molar uap maksimum [kmol/detik]

Diameter Kolom β€’ Kecepatan uap yang diperolehkan tergantung pada: β€’ Jenis talam dan geometrinya, β€’ Kecepatan cairan, β€’ Sifat fisik fasa uap dan cair. 𝐴a 𝐴a 𝜌L βˆ’ 𝜌G π‘ˆG,allow = 𝑓 π‘ˆ =𝑓 𝐾 𝐴Q f 𝐴Q 1 𝜌G f = fraksi flooding AQ = luas penampang lintang kolom =

πœ‹ 2 𝐷 4 c

Aa = luas aktif talam [mΒ²] 𝑀G = berat molekul rata-rata campuran gas

=

𝐺ሢ max 𝑀G 𝜌G π‘ˆG

[mΒ²]

Efisiensi Paking β€’ Tinggi paking aktual: 𝑍 = 𝐻𝐸𝑇𝑃 𝑁teoritis Z = tinggi unggun paking untuk pemisahan yang ekuivalen dengan N tahap teoritis. HETP = height equivalent to a theoretical plate βœ“ Tinggi unggun dimana satu tahap teoritis dilaksanakan. βœ“ Istilah empiris yang hanya berhubungan dengan kecepatan transfer massa setelah dikonversi ke nilai HTU. βœ“ Dipengaruhi oleh: ❑ Konfigurasi geometri paking, ❑ Laju alir, ❑ Rasio kesetimbangan uap-cair, ❑ Sifat fisik uap dan cair.

Paking Acak Cornell dkk (1960) β€’ Fasa cair π»π‘‡π‘ˆ

0,5 Ξ€ = πœ™Sc L L 𝐢 𝑍 10

0,15

β€’ Fasa gas:

πœ“ScG0,5 𝑑cΞ€12 π‘š 𝑍΀10 π»π‘‡π‘ˆ G = 𝐿m 𝑓1 𝑓2 𝑓3 0,6

1Ξ€3

Paking Acak HTUL, HTUG = tinggi satuan perpindahan, ft

 dan  = parameter ScL dan ScG = Bilangan Schmidt fasa cair dan gas (= πœ‡ ΀𝜌𝐷) C = faktor koreksi Z = tinggi bagian paking (di antara distributor-distributor), m dc = diameter kolom, in Lm = laju alir massa cair (lb/jam.ft 2) f 1 = koreksi viskositas cair (= L/1,005)0,16, L dalam cP f 2 = koreksi densitas cair (= 1/ L)0,16, L dalam g/cm3 f 3 = koreksi tegangan permukaan (= 72,8/ )0,8,  dalam dyn/cm m = eksponen = 1,24 (ring) dan 1,11 (saddle)

Paking Acak

Paking Acak

Paking Acak

C

Paking Acak Onda (1968) β€’ Fasa cair:

β€’ Fasa gas:

𝐿 π»π‘‡π‘ˆ L = π‘˜L π‘Žw𝜌L

𝐺 π»π‘‡π‘ˆ G = π‘˜G π‘Žw𝑃𝑀G

L = laju alir massa cairan G = laju alir massa uap aw = luas permukaan paking yang basah (m2/m3)

Paking Acak Onda (1968) π‘˜L

𝜌L π‘”πœ‡L

π‘˜G

1Ξ€3

𝑅𝑇 π‘Žp 𝐷G

= 0,0051

1Ξ€3

𝐿 π‘Žwπœ‡L

𝐺 =𝐢 π‘Žp πœ‡G

2Ξ€3

0,7

ScLβˆ’0,5 π‘Žp𝑑p

Ξ€ ScG1 3

π‘Žp𝑑p

DG = koefisien difusi uap [m2/detik] ap = luas permukaan spesifik paking (m2/m3) dp = ukuran nominal paking C = 5,23 (normal) atau 2 (Raschig ring atau Berl saddles)

βˆ’2

0,4

Paking Acak Onda (1968) βˆ’0,05 0,2 π‘Žw = π‘Žp 1 βˆ’ exp βˆ’1,45Re0,1 Fr We L L L

𝜎 𝜎c

βˆ’0,75

ReL = Bilangan Reynolds = (L/apL) FrL = Bilangan Froude = (π‘Žp 𝐿2 Ξ€π‘”πœŒL2) WeL = Bilangan Weber = (𝐿2 ΰ΅—π‘Žp𝜎𝜌L) Berlaku untuk 0,04 < ReL < 500; 1,2x10-8 < WeL < 0,27; 2,5x10-9 < FrL < 1,8x10-2

0,3 < (c/) < 2.

Paking Terstruktur Metode pintas: 𝐻𝐸𝑇𝑃 = 𝐴1 + 𝐡1 𝐹s,80 [m] βœ“ Berlaku pada refluks total dan  = m/(L/V) = 1,0 βœ“ Fs,80 = faktor-F superfisial pada 80% flooding

𝐻𝐸𝑇𝑃 =

2450𝐹s,80 π‘Žp0,5

[m]

Paking Terstruktur β€’ Locket: 4,82 𝜌L βˆ’ 𝜌G 𝐻𝐸𝑇𝑃 = 𝛼

𝛼 = 𝛼p

0,5 πœ‡βˆ’0,006 r

1 + 0,78𝑒 0,00058π‘Žp

πœ‡π‘Ÿ =

πœ‡L πœ‡H2O pada 20 C

ap = luas area spesifik paking (m2/m3)

𝜌G 𝜌L

0,25 2

Paking Terstruktur Model mekanistik (Bravo dkk, 1985): β€’ Fasa gas 0,333 ShG = 0,0338Re0,8 Sc G G

β€’ Fasa cair

Paking Terstruktur ShG = Bilangan Sherwood = (kGdeq/DG) ReG = Bilangan Reynolds = (deqG/G)(UG,eff + UL,eff) ScG = Bilangan Schmidt = (G/GDG) S = Panjang sisi lekukan = diameter ekuicalen [m] 𝐷L dan 𝐷G = koefisien difusi fasa cair dan fasa gas [[m2/detik] π‘ˆL,eff dan π‘ˆG,eff = laju alir efektif cairan dan gas [m/detik] π‘ˆG,eff =

π‘ˆG,s πœ€ 1 βˆ’ β„Žt sinπœƒ

π‘ˆL,eff =

π‘ˆL,s πœ€β„Žt sinπœƒ

Paking Terstruktur ht = holdup cairan total

 = sudut lekukan yand diukur dari bidang horisontal, deg UG,s dan UL,s = kecepatan superdisial gas dan cairan [m/detik]

Related Documents


More Documents from "RINA MARDYAH S"

Modul 3.pdf
December 2019 30
Cerpen.docx
June 2020 20
Bab 3.docx
December 2019 18
Gosdp.docx
April 2020 40
Sdpprint.docx
April 2020 45