Abstrak — Kegiatan  pemilahan pada aktivitas produksi merupakan salah satu kegiatan yang banyak  dilakukan pada dunia industri.  Penggunaan gravity settling chamber merupakan teknologi pemilahan partikel yang paling sedehana dan murah. Maka dalam perancangan suatu gravity settling chamber  perlu diketahui parameter yang mempengaruhi kinerja dari alat tersebut.    Analisa karakteristik aliran pada gravity settling chamber dilakukan dengan metode simulasi numerik menggunakan persamaan Eularian-Lagrangian.    Sebuah aliran multifase udara  dan partikel abu terbang batu bara (fly ash) berdiameter  100  μm  mengalir  melalui  sebuah    gravity settling chamber  dengan variasi diameter hidrolis saluran inlet  sebesar 1/3 1/5 dan 1/7  kali dari diameter hidrolis  ruang pengendap dengan laju kapasitas aliran yang sama.    Selain divariasikan kecepatan inlet aliran, juga dilakukan variasi distribusi diameter dengan ukuran partikel kurang dari 50 μm, ukuran 50 -100 μm dan ukuran 100- 200 μm.  Hasil yang didapatkan dari penelitian ini adalah semakin kecil ukuran diameter partikel maka akan efisiensi partikel yang ditangkap oleh gravity settling chambersemakin rendah.    Ukuran partikel tidak berpengaruh terhadap posisi jatuh dari partikel tersebut.    Semakin besar kecepatan  inlet  tidak memiliki hubungan dengan efisiensi penangkapan partikel.  Pressure drop yang terjadi akan semakin besar sebanding dengan rasio perbandingan diameter hidrolis inlet dan ruang pengendapan, serta jumlah dan posisi secondary flow berperan signifikan pada efisiensi penangkapan partikel.
Kata Kunci — Aliran  multifase  gas-solid,  gravity settling chamber,  diameter partikel.

## I. PENDAHULUAN
EKNOLOGI  pemilahan merupakan suatu kegiatan yang banyak dilakukan dalam dunia industri.    Tujuan dari pemilahan ini pun bermacam-macam mulai dari dilakukan dalam rangka melakukan sortir suatu produk sebelum dikemas maupun dilepas ke lingkungan  seperti pada perangkat pengendali polusi.   Tipe alat pemisah debu dari gas ( gas  – solid separator ) yang paling sederhana adalah tipe gravity settling chamber . Gravity settling chamber merupakan peralatan pemisah gas-solid dengan menggunakan gayagravitasi  sebagai mekanisme pemisahan utamanya. Setiap partikel atau titik massa memiliki kecepatan terminal yakni kecepatan dimana titik massa tersebut akan mencapai kondisi setimbang antara gaya gravitasi, gaya gesek udara dan gaya ke atas dari partikel tersebut.    Pada tipe gravity settling
chamber kecepatan aliran gas dari boiler akan diturunkan hingga mencapai kecepatan terminalnya.    Penurunan kecepatan tersebut terjadi akibat dari perubahan luas penampang aliran secara mendadak pada suatu ruang. Geometri ruang yang tepat akan mampu memisahkan partikel debu dari gas dengan lebih efektif dan efisien.
Aliran multi fase adalah aliran yang terjadi dimana lebih dari satu fase material yang mengalir.  Aliran tersebut bisa terdiri dari campuran fase likuid-gas, gas-solid, likuid-solid maupun gabungan ketiga fase tersebut [1].    Aliran yang terdiri dari partikel berfase gas dan solid didalamnya memiliki tipe aliran yang berbeda.    Aliran yang memiliki fase gas dan solid yang mengalir diklasifikasikan atas aliran bertipe dense phase dan dilute phase. Dalam penentuan dan pendefinisian tipe aliran tersebut masih terdapat beberapa penafsiran khususnya untuk pendefinisian tipe dense phase .   Tapi secara umum penentuan tipe dari aliran gas-solid yang terjadi sebagai acuan dalam penentuan analisa aliran bersifat dilute, maka solid loading ratio (φ) pada persamaan (1) dari partikel yang terbawa fase gas  terhadap massa fase gas tersebut  haruslah kurang dari 15 [2].
𝜑𝜑 = 𝑚𝑚̇𝑝𝑝(𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡 𝑗𝑗𝑗𝑗𝑚𝑚 ⁄ ) 3,6 𝑚𝑚̇𝑗𝑗(𝑘𝑘𝑘𝑘 𝑡𝑡 ⁄ ) aliran maka aliran tersebut bisa dikategorikan pada tipe dense phase .
Penelitian sebelumnya tentang gravity settling chamber dilakukan oleh Kolaitis  dan Founti [3]  melakukan penelitian terhadap gravitysettling  chamber secara numeric dengan memvariasikan distribusi partikel dan kecepatan inlet fase gas (besar kapasitas aliran berbeda).    Dari penelitian tersebut didapatkan hasil bahwa semakin kecil diameter partikel yang dipisahkan dan semakin besar kecepatan inlet aliran    maka semakin kecil pula efisiensi penangkapan gravity settling chamber .    Selain itu Ko [4] melakukan penelitian mengenai sudut diffuserinlet sebesar 5 ͦ , 9 ͦ dan  15 ͦ  terhadap karakter aliran multiphase secara numerik.    Dari variasi sudut yang dilakukan bahwasannya semakin kecil sudut partikel maka separasi aliran akan terjadi semakin awal.   Kemudian Shaban [5] melakukan modifikasi saluran inlet dari gravity settling chamber dimana dari analisa numerik yang dilakukan menunjukkan bahwa menambah panjang saluran inlet maka efisiensi penangkapan partikel akan semakin meningkat.

## Studi Numerik Karakterisasi Aliran 3 Dimensi Multifase (Gas-Solid) pada Gravity settling chamber
dengan Variasi Kecepatan Inlet dan Diameter Partikel pada Aliran Dilute Phase
Adi Mochammad Isa’i dan Wawan Aries Widodo Jurusan Teknik Mesin, Fakultas Teknologi Industri, Institut Teknologi Sepuluh Nopember (ITS) Jl.  Arief Rahman Hakim, Surabaya 60111
E-mail : wawanaries@me. its. ac. id
T
(1)

## II. METODE PENELITIAN
Gravity settling chamber dalam penelitian ini didesain untuk memisahkan partikel abu terbang ( fly ash ) berdiameter minimal 100 µm dengan densitas partikel sebesar 650 kg/m 3 . Pada penelitian dilakukan 2 variasi yaitu distribusi diameter partikel yang dipisahkan dan geometri saluran inlet pada kapasitas laju aliran yang sama sebesar 5 m3/s.    Variasi diameter partikel (d p )  yang dipisahkan adalah sebesar kurang dari 50μm,  50-100  μm  dan 100-200  μm.    Sedangkan untuk variasigeometri saluran inlet dilakukan dengan acuan besar diametr hidrolis saluran inlet terhadap diameter hidrolis settlingchamber (Dh) sebesar 1/3 Dh (Dh1), 1/5 Dh (Dh 2) dan 1/7 Dh (Dh 3).  Secara  keseluruhan  desain yang dibuat dengan metode mengacu pada Theodore [6] ditunjukkan pada Gambar 1.

## A. Metode Numerik
Penelitian numerik dilakukan dengan menggunakan metode Computational Fluid Dynamics ( CFD )  dalam penyusunan model  persamaan  dan melakukan diskritisasi ( meshing ) pada geometri gravity settling chamber yang sudah dibuat. Prosedur yang dilakukan pada penelitian numerik terbagi menjadi dua tahap yaitu:
1. Pre-processing
Tahap pre-processing merupakan tahap permodelan geometri dan dimensi dari gravity settling chamber kemudian  dilanjutkan dengan proses meshing dengan kerapatan mesh tidak homogen  untuk mendapatkan nilai aspect ratio yang baik.    Kondisi batas yang digunakan adalah inlet adalah velocity inlet .  Kondisi batas semua outlet adalah outflow/escape ,  dan  kondisi batas untuk dinding adalah wall /reflect seperti pada Gambar 2 .

## 2. Processing
Tahap simulasi numerik dilakukan dengan persamaan Eularian-Lagrangian untuk menentukan trayektoori partikel dan  model turbulensi Realizable  k- ε .  Dalam penelitian ini gaya-gaya lain selain gaya gravitasi diabaikan seperti gaya Shaffman dan gaya interaksi partikel yang lain. Sehingga persamaan yang digunakan adalah seperti pada persamaan (2).
𝑑𝑑𝑢𝑢𝑝𝑝 𝑑𝑑𝑡𝑡 = 𝐹𝐹 𝑑𝑑 �𝑣𝑣 − 𝑣𝑣 𝑝𝑝 � + 𝑘𝑘 �𝜌𝜌𝑝𝑝−𝜌𝜌� 𝜌𝜌
Simulasi  diasumsikan dilakukanpada tekanan 1 atm dan pengaruh gravitasi ke arah sumbu y negatif sebesar 9,8 m/s 2 .  Proses pemberian nilai dari boundary condition yaitu pada daerah inlet dengan nilai kecepatan ditentukan seperti pada Tabel 1.  Sedangkan outlet adalah outflow/escape yang diasumsikan langsung bersinggungan dengan kondisi atmosfer. Wall /reflect  merupakan batasan semua dinding gravity settling chamber .
Pada Tabel 2 dijelaskan parameter partikel yang mengalir kedalam gravity settling chamber, dimana dengan menggunakan persamaan (1) didapatkan nilai solid loading ratio bernilai 6.   Sehingga aliran bisa dianalisa sebagai dilute phase flow.
Gambar 1. Geometri gravity settling chamber.
Gambar 2. Boundary conditions permodelan.
Tabel 1. Kecepatan inlet udara.
Parameter laju aliran inlet Kecepatan (m/s) Laju aliran udara 1 ( Dh 1) 8 Laju aliran udara 2 (Dh 2) 20,4 Laju aliran udara 3 (Dh 3) 40 Tabel 2. Parameter aliran partikel solid. Aliran Partikel Diameter velocity (m/s) mass flow rate (kg/s) 1- 50 μm 5 30 50- 100 μm 5 30 100- 200 μm 5 30
Solusi yang digunakan dalan simulasi ini adalah second order upwind untuk momentum dan turbulence dan first order untuk  tekanan , sedangkan SIMPLE coupling solution dipilih untuk simulasi ini.  Langkah berikutnya adalah initialize yang merupakan  perhitungan awal untuk memudahkan dalam mendapatkan hasil yang konvergen pada tahap iterasi dihitung dari inlet .  Kriteria konvergensi ditetapkan sebesar 10 -3 , artinya
(2)
proses iterasi dinyatakan telah konvergen setelah residualnya mencapai harga lebih kecil daripada 10 -3 atau besar mass flux inlet dan outlet kurang dari 10 -3 .

## III. HASIL DAN PEMBAHASAN

## A. Performa Gravity settling chamber Terhadap Variasi Distribusi Partikel.
Dari  jumlah partikel yang melewati maasing- masing outlet beserta distribusi diameter yang melewatinya ditampilkan pada histogram pada  Gambar 2.    Pada  Gambar 2  terlihat bahwa pada partikel dengan distribusi diameter antara 100- 200 μm partikel paling banyak keluar melewati outlet 1 dan paling sedikit keluar melalui main outlet .    hasil yang menunjukkan bahwa partikel dengan distribusi diameter tersebut telah lebih banyak jatuh ke ketiga penampung ( outlet ) dan sedikit sangat keluar melalui main  outlet mengindikasikan efisiensi penangkapan partikel sangat baik yakni 99,7 %.  Dari simulasi distribusi partikel fly ash berdiameter kurang dari 100 μm menunjukkan  tren efisiensi penangkapan partikel menurun seiring bertambah kecilnya partikel yang dipisahkan.
B. Performa Gravity settling chamber TerhadapVariasi Diameter Hidrolis.
1. Variasi Diameter Hidrolis Terhadap Distribusi Partikel
Yang Melewati Tiap Outlet .
Dengan kapasitas laju aliran yang sama, maka perbedaan ukuran geometri saluran inlet tersebut akan menyebabkan laju aliran pada masing-masing saluran inlet akan menjadi berbeda, sehingga variasi geometri saluran inlet ini merupakan juga variasi dari kecepatan inlet aliran fluida. Pada  Gambar  4  ditampilkan data kuantitatif efisiensi penangkapan partikel dari tiap variasi saluran inlet .
Pada Gambar  4  dari kesemua  variasi diameter hidrolis saluran inlet mempunyai tren yang sama bahwa semakin besar ukuran diameter  maka efisiensi penangkapan partikel akan semakin besar.    Kesimpulan tersebut diketahui dari prosentase partikel yang melalui saluran outlet utama ( main outlet) .  semakin kecil jumlah partikel yang melewati saluran  main outlet maka semakin tinggi efisiensi penangkapan partikel oleh gravity settling chamber .    Dari hasil pada Gambar  4  di atas dalam memisahkan partikel dengan  ukuran  lebih  dari  100  μm  ,  diameter  hidrolis  1 memiliki efisiensi penangkapan sebesar 99,7 %, kemudian untuk diameter hidrolis 3 dan 2 maing-masing sebesar 96,98 % dan 99,54%.    Pada variasi diameter hidrolis 2 terjadi fenomena semakin kecil diameter partikel tidak berhubungan lurus dengan posisi jatuh partikelnya. Fenomena  tersebut  dimungkinkan terjadi karena adanya secondary flow di dalam ruang pengendap dimana hubunganbesar diameter partikel terhadap posisi jatuhnya tidak  dapat diperoleh  dari perhitungan analitis dengan persamaan yang ada.    Dari fenomena ini dapat juga disimpulkan bahwa posisi jatuh dari partikel dengan diameter tertentu tidak dapat dipastikan dengan hanya dari perhitungan analitis biasa namun perlu dilakukan uji simulasi numerik maupun eksperimen.
Gambar. 3 . Gravik ditribusi diameter melalui tiap outlet.
Gambar. 4.  Distribusi diameter partikel tiap variasi saluran inlet .
Gambar 5. Gravik efisiensi penangkapan partikel tiap variasi diameter hidrolis.

## 2. Variasi Diameter Hidrolis Terhadap Performa Penangkapan Partikel Dan Penurunan Tekanan.
Proses penangkapan partikel pada ruang pengendap dilakukan perlambatan tiba-tiba pada gerakan partikel. Perlambatan tiba-tiba tersebut berpengaruh terhadap terjadinya pressure drop pada aliran.
Dari  Gambar  5  didapatkan hasil bahwa efisiensi penangkapan terbesar diperoleh dari variasi diameter hidrolis 1 , kemudian diameter hidrolis 3 dan  terakhir diameter hidrolis 2.    Variasi diameter hidrolis yang dimaksudkan sebagai parameter variasi kecepatan pada kapasitas aliran gas-solid, menunjukkan hasil bahwa hubungan antara kecepatan inlet dengan efisiensi penangkapan tidak menunjukkan hubungan yang signifikan.
88,3583,787,55 96 90,0294,93 99,7 95,9899,54 75 95 Dh 1Dh 2Dh 3 ens i P ena ng ka pa n % Variasi  Diameter Hidrolis Inlet Effisiensi vs Variasi Diameter Hidrolis Inlet d<50 50<d<10
0
Gambar. 6. Grafikk penurunan tekanan pada tiap variasi.
Menurut  Kolaitis & Founti [2],terjadinya  resirkulai aliran ( secondary flow ) akan membantu dalam peningkatan penangkapan partikel pada gravity settling  chamber . Semakin banyak secondary flow yang terbentuk akan memerangkap partikel hingga partikel tersebut kehilangan sejumlah energi kemudian partikel tersebut akan jatuh karena efek dari gaya gravitasi.    Namun perlu diketahui juga semakin banyak terbentuk secondary flow pada suatu aliran akan meningkatkan pressure drop .
Pada  Gambar6  ditampilkan data kuantitatif besar penurunan tekanan pada masing-masing variasi diameter hidrolis inlet .  Pada gravik yang ditampilkan pada Gambar5 menunjukkan bahwa semakin kecil diameter hidrolis dari saluran inlet maka semakin besar penurunan tekanan yang terjadi. Dengan membandingkan hasil pada Gambar 5 dan 6 diperoleh fenomena bahwa peningkatan efisiensi penangkapan partikel tidak berbanding lurus dengan penurunan tekanan yang terjadi.
Adanya ketidaklinieran hubungan antara efisiensi penangkapan partikel dengan penurunan tekanan di jelaskan pada Gambar 7.  Dari visualisasi aliran pada gambar 7 terlihat adanya perbedaan jumlah secondary flow pada tiap variasi diameter hidrolis.  Dari penelitian Kolaitis & Founti [1]  menyebutkan bahwa banyaknya secondary flow yang terbentuk akan menambah efisiensi partikel , tetapi dari visualisasi pada gambar 5 , 6 dan 7 didapatkan bahwa meskipun secondary flow pada variasi diameter hidrolis 1 terlihat tidak sebanyak pada diameter hidrolis 2 dan 3 namun efisiensi penangkapan partikel paling tinggi. Kemudian pada diameter hidrolis 2 terlihat bahwa posisi dari secondary flow menyebabkan aliran seakan’diantar’ keluar main  outlet sehingga menyebabkan  efisiensi penangkapan partikel paling rendah.  Dari penelitian ini dapat memberikan koreksi bahwa pernyataan hasil penelitian dari Kolaitis & Founti [3] perlu ditambah dengan pernyataan bahwasannya posisi dari secondary  flow pada tiap daerah pada gravity settling chamber akan menentukan bentuk streamline aliran dimana konfigurasi titik secondary flow yang timbul akan menentukan apakah streamline aliran cenderung ’diantar’ keluar menuju main outlet ataukah akan menyebabkan  terjadinya backflow sehingga partikel akan terperangkap di dalam gravity settling chamber dan akibat penurunan energi yang dialami partikel ,partikel akan jatuh ke dalam hopper karena efek gaya gravitasi sehingga efisiensi penangkapan partikel dapat lebih optimal.
(A) Visualisasi pathline aliran penampang midspan (B) Visualisasi pathline aliran penampang midspan (C)
Gambar 7.  (A) Visualisasi pathline aliranpada variasi diameter hidrolis 1. (B) Visualisasi pathline aliranpada variasi diameter hidrolis 2.  (C) Visualisasi pathline aliranpada variasi diameter hidrolis 3.

## IV. KESIMPULAN
Berdasarkan hasil pembahasan pada penilitian variasi diameter partikel dan kecepatan inlet kepada unjuk kerja gravity settling chamber didapatkan kesimpulan bahwa semakin besar ukuran dari diameter partikel maka efisiensi penangkapan oleh gravity settling chamber akan semakin rendah.  Posisi jatuh partikel dengan diameter tertentu tidak dapat ditentukan secara analitis namun harus dilakukan uji simulasi maupun ekspermen untuk mengetahui disribusi posisi jatuh partikel terhadap ukuran diameter partikel. Kecepatan inlet tidak berhubungan langsung terhadap efisiensi penangkapan partikel pada kapasitas aliran yang
39 7346,949 25671,149 0 10000 20000 30000 Dh 1 Dh 2 Dh 3 Pasc al Variasi Pressure Drop vs variasi Diameter Hidrolis inlet Pressure Drop vs variasi
Diameter idrolis inlet
sama.  Semakin besar rasio penampang inlet dan settling chamber maka pressure drop akan semakin besar. Efisiensi penangkaan partikel di dalam gravity settling chamber selain dipengaruhi efek banyaknya secondary flow yang terjadi pada aliran juga dipengaruhi oleh  posisi terjadinya secondary flow .    Posisi secondary flow tersebut berpengaruh besar terhadap bentuk streamline aliran  di dalam gravity settling chamber .
DAFTAR SIMBOL φ : solids loading ratio 𝑚𝑚̇ 𝑝𝑝 : Laju aliran massa partikel solid (ton/jam)
ṁ a : Laju alir massa  fluida (kg/s)
F d : Gaya drag partikel (N) g : percepatan gravitasi (m/s 2 ) v t : kecepatan terminal partikel (m/s)
d p : Diameter partikel (m) ρ : Densitas fluida (kg/m 3 ) ρ p : Densitas partikel (kg/m 3 ) A p : Luas proyeksi partikel (m 2 ) m p : Massa partikel (kg)
v
: Kecepatan aliran fluida di dalam gravity settling chamber (m/s) u p : Kecepatan Partikel (m/s)

## DAFTAR PUSTAKA
[1] Perry,Robert H and Green, Don w.  1999 “Chemical Engineer’s Haandbook “ . New York: Mc Graw Hill Inc.
[2] Mills,David.  2004. ”Pneumatic Conveying Design Guide” .  Oxford : Elsevier Butterworth-Heinemann
[3] Dionysios I.    Kolaitis, and Maria A.    Founti.  Modelling of gas-particle flow in industrial classification chambers for design optimization. Powder Technology 125 (2002) 298 – 305.
[4] T.  H, Ko.  Numerical Study on three-dimensional isothermal solid-gas flow and deposition process in a deposition chamber.    International communications in heat and mass transfer 33 (2006) 601-609.
[5] S.  Shaaban.  Numerical Optimization and Expermental investigation of the aerodynamic performance of a three stage gas solid separator, Chemical Engineering Research and Design 89(2011) 29-38.
[6] Theodore, Louis.    2008.    “ Air Polution Control Equipments and Calculations” .  New York: John Willey & Sons, Inc.