UNIVERSITAS DIPONEGORO
SIMULASI PEMISAHAN FLUIDA 2 FASA DI DALAM COALESCER DENGAN MENGGUNAKAN CFD
TUGAS AKHIR
M. ANDRIANTO ADITYA PUTRA L2E 005 461
FAKULTAS TEKNIK JURUSAN TEKNIK MESIN
SEMARANG MARET 2011
i
ii
iii
iv
v
ABSTRAK
Coalescer digunakan dalam beberapa aplikasi industri. Dalam kilang, coalescer digunakan dalam pemisahan larutan dari hidrokarbon seperti dalam pemisahan dispersi selama tahap pencucian. Dalam aplikasi lepas pantai coalescer digunakan untuk memisahkan air dari minyak dan juga untuk air dan gas alam glikol dari kondensat. Fungsi coalescer yaitu untuk memisahkan minyak dan air. Proses pemisahan coalescer oil water separators adalah air dan minyak masuk separator dan menyebar horizontal, didistribusikan melalui energi dan turbulensi akan menyebar ke perangkatnya. Air dan minyak memasuki flopak media dan fungsi flopak media yaitu untuk menyaring kotoran minyak dan air dan melaminarkan aliran fluida untuk mempercepat aliran fluida. Dalam penelitian ini, Computational Fluid Dynamics (CFD) digunakan untuk mengetahui berapa kandungan minyak yang keluar dari keluaran air. Simulasi ini dilakukan dengan memvariasikan jarak sekat di dalam coalescer yaitu sekat 1 inch, 1.2 inch, dan 1.5 inch. Hasil simulasi coalescer ini mendapatkan hasil kandungan minyak dalam ppm ( part per million ) yang paling sedikit yaitu pada sekat 1 inch yaitu 27.9 ppm. Semakin kecil jarak sekat didalam coalescer maka semakin kecil hasil kandungan minyak yang didapatkan karena minyak akan menempel di flopak media dan akan menghambat aliran minyak karena pengaruh viskositas.
Kata kunci : Coalescer, flopak media, jarak sekat didalam coalescer, viskositas
vi
ABSTRACT
Coalescer is used in several industrial applications. In refineries, coalescer used in the separation of hydrocarbon solution such as the separation of the dispersion during the washing stage. In the offshore applications coalescer is used to separate water from oil and natural gas for water and glycol from the condensate. Coalescer is used to separate oil from water and all mixture. In the process of separation of oil water coalescer separators water and oil flow into the separator and spread horizontally, distributed through the energy and turbulence would spread to the device. Water and oil entering flopak media, flopak media filter out dirt oil and water and serves to accelerate the flow. In this study, Computational Fluid Dynamics (CFD) is used to determine how much oil content out of the water output. The simulation was done by varying the distance bulkhead in the coalescer is 1 inch, 1.2 inch and 1.5 inch. The simulation results of this coalescer shows that the oil content in ppm (parts per million), which is at least 1 inch bulkhead, is 27.9 ppm. The smaller the distance bulkhead inside the coalescer, the smaller the oil content is obtained because the oil will stick in the flopak media and will inhibit the flow of oil due to the influence of viscosity.
Keywords : Coalescer, flopak media, the distance bulkhead in coalescer, viscosity
vii
MOTTO
Sesungguhnya, Aku mengingatkan kepadamu supaya kamu tidak termasuk orang-orang yang tidak berpengetahuan. (QS Hud : 46) Tak ada rahasia untuk menggapai sukses. Sukses itu dapat terjadi karena persiapan, kerja keras, dan mau belajar dari kegagalan. (General Colin Powell)
PERSEMBAHAN
Kupersembahkan Tugas Sarjana ini kepada Ibu, Bapak, Adikku yang tercinta... Teman Teknik Mesin Angkatan 2005 UNDIP... Terima kasih atas berbagai dukungan dan doa yang telah diberikan...
viii
KATA PENGANTAR
Segala puji syukur senantiasa penulis panjatkan kepada Allah SWT, karena berkat rahmat-Nya, penulis dapat menyelesaikan laporan tugas akhir ini dengan judul “SIMULASI PEMISAHAN FLUIDA 2 FASA DI DALAM COALESCER DENGAN MENGGUNAKAN CFD”. Tugas akhir ini merupakan salah satu syarat yang harus dipenuhi pada program strata satu (S1) di Jurusan Teknik Mesin, Fakultas Teknik, Universitas Diponegoro Semarang. Pada kesempatan ini penulis ingin mengucapkan terima kasih atas bimbingan, bantuan, serta dukungan kepada : 1. Dr. MSK. Tony Suryo Utomo, ST, MT, selaku Dosen Pembimbing I dan Koordinator Tugas Akhir. 2. Dr.Ir.Dipl.Ing. Berkah Fajar T K, selaku Ketua Jurusan Teknik Mesin Universitas Diponegoro Semarang. Dalam penulisan tugas akhir ini penulis menyadari banyak kekurangan. Oleh karena itu segala kritik yang bersifat membangun akan diterima dengan senang hati untuk kemajuan bersama. Akhir kata penulis berharap semoga laporan tugas akhir ini dapat memberikan manfaat kepada siapa saja yang membutuhkan data maupun referensi yang ada dalam laporan ini. Terima kasih.
Semarang, 29 Maret 2011
Penulis
ix
DAFTAR ISI
JUDUL .................................................................................................................. i TUGAS SARJANA .............................................................................................. ii HALAMAN PERNYATAAN ORISINALITAS .................................................. iii HALAMAN PENGESAHAN............................................................................... iv ABSTRAK ............................................................................................................ v ABSTRACT ............................................................................................................ vi MOTTO dan PERSEMBAHAN ........................................................................... vii KATA PENGANTAR .......................................................................................... viii DAFTAR ISI ......................................................................................................... ix DAFTAR GAMBAR ............................................................................................ xii DAFTAR TABEL ................................................................................................. xv NUMENKLATUR ................................................................................................ xvi BAB I PENDAHULUAN .................................................................................. 1 1.1. Latar Belakang............................................................................... 1 1.2. Batasan Masalah ............................................................................ 4 1.3. Tujuan Penelitian ........................................................................... 4 1.4. Metodologi Penelitian ................................................................... 5 1.5. Sistematika Penulisan .................................................................... 6 BAB II DASAR TEORI...................................................................................... 7 2.1. Klasifikasi Aliran........................................................................... 7 2.1.1
Aliran Invicid dan Viscous................................................. 8
2.1.2
Aliran Laminar dan Turbulen ............................................ 9
2.1.3
Aliran Kompresibel dan Inkompresibel ............................ 10
2.1.4
Aliran Internal dan Eksternal............................................. 10
2.1.5
Aliran 2 Fasa...................................................................... 11
2.2. Persamaan Dasar Aliran Fluida ..................................................... 11 2.2.1 Persamaan Kekekalan Massa............................................. 11 2.2.2 Persamaan Kekekalan Momentum .................................... 14 2.2.3 Persamaan Navier-Stokes untuk Sebuah Fluida Newtonian ........................................................................................... 18 x
2.3. Bilangan Reynolds......................................................................... 21 2.4. Coalescer ....................................................................................... 21 2.5. Design coalescer ............................................................................ 23 BAB III METODA VOLUME HINGGA DALAM FLUENT 6.2.16 Dan PERMODELAN SIMULASI ............................................................... 27 3.1.
Fluent ............................................................................................ 27
3.2.
Skema numerik ............................................................................. 30 3.2.1 Metode Solusi Segregated ................................................. 30 3.2.2 Metode solusi Coupled ...................................................... 31
3.3.
Diskritisasi (Discretization) .......................................................... 32 3.3.1 First-Order Upwind ........................................................... 34 3.3.2 Second-Order Upwind Scheme .......................................... 34 3.3.3 Bentuk Linearisasi Persamaan Diskrit ............................... 35 3.3.4 Under-Relaxation .............................................................. 35
3.4.
Diskretisasi Coupled Solver .......................................................... 36 3.4.1 Aturan Persamaan-persamaan dalam Bentuk Vektor ........ 36 3.4.2 Preconditioning ................................................................. 37
3.5.
Model Multifasa ............................................................................ 39
3.6.
Model Turbulen (Turbulence Models) .......................................... 40 3.6.1 Permodelan k-epsilon (k-ε) ................................................ 40 3.6.1.1 Standard ............................................................................. 40 3.6.1.2 RNG ................................................................................... 41 3.6.1.3Realizable ........................................................................... 41 3.6.2 Permodelan k-omega (k-ω) ................................................ 42 3.6.2.1 Standard ............................................................................. 42 3.6.2.2 SST .................................................................................... 42
3.7.
Adapsi Grid ................................................................................................. 43 3.7.1 Proses Adapsi..................................................................... 43 3.7.2 Adapsi Y+ dan Y*............................................................... 45
xi
3.8.
Permodelan Bentuk dan Simulasi Numerik pada Fluent 6.3.26 ....................................................................................................... 47 3.8.1 Permodelan Bentuk............................................................ 47 3.8.1.1 Bentuk Coalescer Sekat 1.2 inch ....................................... 47 3.8.1.2 Bentuk Coalescer Sekat 1 inch .......................................... 48 3.8.1.3 Bentuk Coalescer Sekat 1.5 inch ....................................... 48 3.8.2 Pembentukan Mesh ............................................................ 49 3.8.3 Solver dan Model Turbulensi ............................................ 50 3.8.4 Pendefinisian Kondisi Batas .............................................. 52
BAB IV Hasil Analisa......................................................................................... 55 4.1.1 Hasil Simulasi dengan Beberapa model k-epsilon (k-ε) dan k-omega (k-ω) .................................................................... 55 4.2
Penggambaran Display dan Velocity Vectors didalam coalescer ............................................................................ 68
4.2.1 Display Contours Velocity ................................................ 68 4.2.1 Display Contours Pressure................................................. 70 4.2.2 Velocity Vectors inlet ........................................................ 72 4.2.3 Velocity Vectors Volume Fraction .................................... 73 BAB V Kesimpulan dan Saran ........................................................................... 74 5.1
Kesimpulan ........................................................................ 74
5.2
Saran .................................................................................. 74
DAFTAR PUSTAKA ........................................................................................... 75
xii
DAFTAR GAMBAR
Gambar 1.1
Gambar Instalasi pada Offshore .................................................... 1
Gambar 1.2
Coaleser ......................................................................................... 2
Gambar 1.3
Aliran Coalescer ............................................................................ 3
Gambar 1.4
Flopak media design ..................................................................... 4
Gambar 1.5
Flopak media design ..................................................................... 4
Gambar 2.1
Respon terhadap gaya ................................................................... 7
Gambar 2.2
Klasifikasi aliran ........................................................................... 8
Gambar 2.3
Berbagai daerah aliran lapisan batas diatas plat rata ..................... 9
Gambar 2.4
Variasi kecepatan (satu dimensi) terhadap waktu ....................... 10
Gambar 2.5
Aliran dalam pipa .......................................................................... 11
Gambar 2.6
Massa mengalir kedalam dan keluar elemen fluida ...................... 13
Gambar 2.7
Komponen tegangan pada tiga bidang elemen fluida ................... 15
Gambar 2.8
Komponen tegangan dalam arah x ................................................ 16
Gambar 2.9
Coalescer ....................................................................................... 21
Gambar 2.10 Konfigurasi Coalescer ................................................................... 22 Gambar 2.11 Oleophilic Material ....................................................................... 22 Gambar 2.12 Pelat pemisah. ............................................................................... 24 Gambar 2.13 Pemisah pelat datar bergelombang................................................ 25 Gambar 2.14 Penggabungan tabung pemisah ..................................................... 25 Gambar 2.15 Pemisah jenis ................................................................................. 26 Gambar 3.1
Diagram Alir Prosedur Simulasi FLUENT ................................... 29
Gambar 3.2
Volume
Kendali
Digunakan
Sebagai
Ilustrasi
Diskretisasi
Persamaan Transport Skalar ......................................................... 33 Gambar 3.3
Volume Kendali Digunakan Sebagai Ilustrasi Diskretisasi Persamaan Transport Skalar pada model sel 2D quadrilateral ..... 34
Gambar 3.4
Contoh Hanging Node .................................................................. 44
Gambar 3.5
Hasil Conformal Refinement ......................................................... 45
Gambar 3.6
dengan meenghilangkan titik dan menyetigakan kembali ............ 45 xiii
Gambar 3.7
Hasil Meshing Menggunakan Adapsi
Gambar 3.8
Gambar Coalescer sekat 1.2 inch .................................................. 47
Gambar 3.9
Gambar Coalescer sekat 1 inch ..................................................... 48
.................................. 46
Gambar 3.10 Gambar Coalescer sekat 1.5 inch .................................................. 48 Gambar 3.11 Gambar Mesh Coalescer ............................................................... 50 Gambar 3.12 Model k epsilon standard .............................................................. 51 Gambar 3.13 Model k epsilon RNG ................................................................... 51 Gambar 3.14 Model Mass flow inlet .................................................................. 52 Gambar 3.15 Model Mass flow inlet water......................................................... 53 Gambar 3.16 Model Mass flow inlet oil ............................................................. 53 Gambar 3.17 Model outflow oil .......................................................................... 54 Gambar 3.18 Model outflow water ..................................................................... 54 Gambar 4.1
Perbandingan grafik hasil berapa kandungan minyak yang terjadi pada jarak sekat 1 inch, 1.2 inch dengan model k-epsilon (k-ε) standard ......................................................................................... 56
Gambar 4.2
Perbandingan grafik hasil berapa kandungan minyak yang terjadi pada jarak sekat 1 inch, 1.2 inch dengan model k-epsilon (k-ε) RNG .............................................................................................. 57
Gambar 4.3
Perbandingan grafik hasil berapa kandungan minyak yang terjadi pada jarak sekat 1 inch, 1.2 inch dengan model k-omega (k-ω) standard ......................................................................................... 58
Gambar 4.4
Perbandingan grafik hasil berapa kandungan minyak yang terjadi pada jarak sekat 1 inch, 1.2 inch dengan model k-omega (k-ω) SST ................................................................................................ 59
Gambar 4.5
Perbandingan grafik hasil berapa kandungan minyak yang terjadi pada jarak sekat 1 inch, 1.2 inch dengan model k-epsilon (k-ε) standard ......................................................................................... 60
Gambar 4.6
Perbandingan grafik hasil berapa kandungan minyak yang terjadi pada jarak sekat 1 inch, 1.2 inch dengan model k-epsilon (k-ε) RNG .............................................................................................. 61
xiv
Gambar 4.7
Perbandingan grafik hasil berapa kandungan minyak yang terjadi pada jarak sekat 1 inch, 1.2 inch dengan model k-omega (k-ω) standard ..................................................................................... 62
Gambar 4.8
Gambar grafik perbandingan permodelan pada jarak sekat 1 inch ................................................................................................... 64
Gambar 4.9
Gambar grafik perbandingan permodelan pada jarak sekat 1.2 inch ................................................................................................... 65
Gambar 4.10 Gambar grafik perbandingan permodelan pada jarak sekat 1.5 inch ................................................................................................... 66 Gambar 4.11 Gambar contours velocity sekat 1 inch pada detik ke 180 ........ 68 Gambar 4.12 Gambar contours velocity sekat 1.2 inch pada detik ke 180 ..... 68 Gambar 4.13 Gambar contours velocity sekat 1.5 inch pada detik ke 180 ..... 69 Gambar 4.14 Gambar contours pressure sekat 1 inch pada detik ke 180 ....... 70 Gambar 4.15 Gambar contours pressure sekat 1.2 inch pada detik ke 180 .... 70 Gambar 4.16 Gambar contours pressure sekat 1.5 inch pada detik ke 180 .... 71 Gambar 4.17 Gambar velocity vectors inlet didalam coalecser pada detik ke 180 pada sekat 1 inch ....................................................................... 72 Gambar 4.18 Gambar velocity vectors by volume fraction campuran minyak dan air didalam coalecser pada detik ke 180 pada sekat 1 inch ...... 73
xv
DAFTAR TABEL
Tabel 2.1 Persamaan Pembangun Aliran Fluida Newtonian Kompresibel .........20 Tabel 3.1 Characteristic .......................................................................................49 Tabel 3.2 Jenis Volume Meshing ........................................................................49 Tabel 4.1 Hasil Simulasi Model k-epsilon (k-ε) Standard ...................................55 Tabel 4.2 Hasil Simulasi Model k-epsilon (k-ε) RNG .........................................56 Tabel 4.3 Hasil simulasi model k-omega (k-ω) Standard ....................................57 Tabel 4.4 Hasil simulasi model k-omega (k-ω) SST ...........................................58 Tabel 4.5 Hasil Simulasi Model k-epsilon (k-ε) Standard ...................................60 Tabel 4.6 Hasil Simulasi Model k-epsilon (k-ε) RNG .........................................61 Tabel 4.7 Hasil simulasi model k-omega (k-ω) Standard ....................................62 Tabel 4.8 Hasil Simulasi Model 1 inch ...............................................................64 Tabel 4.9 Hasil Simulasi Model 1.2 inch ............................................................65 Tabel 4.10 Hasil Simulasi Model 1.5 inch ........................................................... 66
xvi
NOMENKLATUR
c
Kecepatan suara
D
Diameter Pipa
m
I
Intensitas turbulen
%
M
Bilangan Mach Number
m massa
kg
P
tekanan
psig
Re
reynold number
t
Waktu
s
Kecepatan V
Kecepatan aliran viskositas Densitas
xvii
BAB 1 PENDAHULUAN I.
LATAR BELAKANG Coalescer digunakan dalam beberapa aplikasi industri. Dalam kilang,
coalescer digunakan dalam pemisahan larutan dari hidrokarbon seperti dalam pemisahan dispersi selama tahap pencucian. Dalam aplikasi lepas pantai coalescer digunakan untuk memisahkan air dari minyak dan juga untuk air dan gas alam glikol dari kondensat. Dalam industri hidrometalurgi, coalescer digunakan dalam pelarut untuk menghilangkan dari zat cair dan juga dalam ekstraksi untuk meningkatkan pemisahan fasa. Coalescer berguna dalam aplikasi proses industri kimia misalnya, dalam menghilangkan air dari pelarut [1].
Gambar 1.1 Gambar process Flow Diagram Produksi Dari sumur reservoir fluida dialirkan ke separator melalui pipa-pipa alir (flow line) untuk dilakukan pemisahan antara air, minyak dan gas bumi. Gas hasil pemisahan, selain dapat langsung dimanfaatkan untuk industri dapat pula digunakan untuk injeksi gas-lift atau pressure maintenance, sedangkan minyak bumi (crude oil) umumnya ditampung terlebih dahulu di pusat pengumpulan minyak sebelum dikirim ke pengilangan atau terminal untuk dikapalkan [2],
1
sedangkan air akan memasuki coalescer karena masi mengandung minyak. Coalescer digunakan untuk memisahkan air dari minyak. Proses pemisahan coalescer oil water separators yaitu air dan campuran minyak masuk separator dan menyebar horizontal, didistribusikan melalui energi dan turbulensi akan menyebar ke perangkatnya. Campuran memasuki media Flopak yang ada pada laminer dan arus sinusoidal dibentuk dan minyak pada permukaan media. Minyak akan menyatu menjadi tetesan yang lebih besar, naik ke atas melalui corrugations, di mana minyak dan air dapat terpisah dan naik ke permukaan air. Pada saat yang bersamaan material padat memasuki media lalu memasuki proses pembentukan (corrugation)
setelah itu
jatuh ke dalam v-
hopper dibawah media flopak [3]. Ada 2 bahan yang digunakan untuk plat coalescer yaitu plat polypropylene dan polietilena densitas tinggi. Polypropylene lebih oleophilic ( karakteristik zat nya lebih menyerupai minyak dari pada air ) dan menjadi bahan yang digunakan dalam produksi minyak secara umum. Beberapa orang industri percaya bahwa polypropylene terlalu oleophilic ( karakteristik zat nya lebih menyerupai minyak dari pada air) dan tidak memungkinkan minyak untuk berpindah ke permukaan air. Bahan Polietilena umumnya digunakan sebagai material aleternatif untuk perbandingan [4].
Gambar 1.2 Coalescer [3] 2
Gambar 1.3 Aliran Coalecer [5]
Media flopak menyediakan pola aliran laminar untuk membuat zona tanpa gerak sehingga memungkinkan tumbukan antara kandungan pada minyak dan permukaan media flopak dengan cara mengurangi kecepatan dan turbulensi aliran limbah. Kontrol terhadap aliran ini menciptakan keadaan ideal untuk terjadinya pemisahan minyak. Berdasarakan desain media flopak, material padat juga akan bertumbukan dengan media flopak dan tetap berada didasar coalescer dalam waktu tertentu, karena minyak lebih ringan dari pada air maka minyak tersebut akan berada diatas coalescer [3] Desain ini di kenal dengan banyak nama yaitu pelat pararel dan pelat bergelombang. Namun desain dan konsepnya umumnya sama. Dan desain ini umumnya kosentrasi pembuangan minyaknya 10 ppm [3].
3
Gambar 1.4 Flopak media design [6]
Gambar 1.5 Flopak media design [3]
II.BATASAN MASALAH Batasan masalah yang diambil dari percobaan ini ada 4 yaitu 1. Coalescer yang digunakan jenis Coalescer 2-fasa minyak dan air. 2. Menggunakan jarak sekat dalam coalescer 1 inch, 1.2 inch dan 1.5 inch. 3. Menggunakan beberapa model k-epsilon (k-ε) dan k-omega (k-ω). 4. Simulasi menggunakan jenis aliran unsteady.
III.TUJUAN Adapun tujuan dari penyusunan Tugas Akhir ini adalah : 1. Untuk memisahkan air dari minyak dan mengetahui berapa kandungan minyak (PPM) yang keluar dari outlet air. 2. Untuk membandingkan beberapa model k-epsilon dan k-omega yang tepat pada saat simulasi.
4
IV.METODE PENELITIAN
Metode penelitian yang digunakan penulis dalam penulisan Tugas Akhir ini adalah: 1.
Studi Pustaka Studi pustaka adalah suatu metode yang dipergunakan dalam penelitian
ilmiah
dengan
memperoleh
informasi
yang
berhubungan dengan hasil-hasil penelitian yang telah dilakukan oleh para peneliti sebelumnya dari literatur yang telah dipublikasikan. 2.
Proses Simulasi Melakukan proses simulasi dengan software FLUENT 6.3.26 dengan input mesh dari GAMBIT 2.3.16 .
3.
Pengolahan dan Analisa Data Data yang diperoleh dari hasil simulasi, penulis melakukan verifikasi
dengan
data
pustaka
dan
melakukan
analisa
perbandingan. 4.
Penyusunan Laporan Setelah tahapan-tahapan diatas selesai dilakukan penyusunan laporan mulai dilakukan, asistensi dilakukan dengan dosen pembimbing
Tugas
Akhir
yang
bersangkutan.
Setelah
mengadakan asistensi dengan dosen dan berdasarkan data-data yang diperoleh, penulis mengambil kesimpulan dan saran mengenai penelitian yang telah dilakukan.
5
V.SISTEMATIKA PENULISAN Dalam pembuatan laporan tugas akhir ini, sistematika penulisan yang dilakukan penulis adalah sebagai berikut :
BAB I
PENDAHULUAN Pada Bab I ini berisi tentang latar belakang, perumusan masalah,
tujuan
penelitian,
pembatasan
masalah,
metodologi, dan sistematika penulisan itu sendiri.
BAB II
LANDASAN TEORI Merupakan landasan teoritis yang digunakan dalam tugas sarjana ini sebagai studi literatur.
BAB III DASAR METODA VOLUME HINGGA DALAM FLUENT 6.3.26 dan Permodelan simulasi Berisi tentang perangkat lunak FLUENT, diskretisasi persamaan-persamaan diferensial dasar serta gambaran umum variabel simulasi yang disediakan oleh FLUENT dan permodelan simulasi. BAB IV HASIL SIMULASI Disini menyajikan hasil simulasi yang telah dilakukan dan memuat tentang verifikasi dari data hasil simulasi dengan literatur yang ada. BAB V
KESIMPULAN DAN SARAN Berisi tentang kesimpulan dan saran yang diambil dari hasil analisa pada bab-bab sebelumnya.
6
BAB II DASAR TEORI 2.1.
Klasifikasi Aliran Fluida adalah zat yang terus menerus mengalami deformasi dibawah
penerapan tegangan geser (tangensial) tidak peduli seberapa kecil tegangan geser. Sehingga fluida terdiri dari cairan dan gas (atau fase uap). Perbedaan antara keadaan fluida dan solid jelas jika anda membandingkan perilaku fluida dan solid. Solid berdeformasi ketika tegangan geser diterapkan, tetapi deformasi yang tidak terus meningkat dengan waktu [7]. Visualisasi konsep fluida dapat dilihat pada Gambar 2.1.
.
a
b
c
Gambar 2.1 Respon terhadap gaya: (a) tidak merespon gaya yang diberikan secara tiba-tiba; (b) beberapa waktu setelah dikenai gaya; (c) setelah lama dikenai gaya [8].
Dari konsep fluida di atas, kemudian dikenal konsep fluida continuum adalah fluida yang mengalir dianggap sebagai suatu kumpulan molekul yang 7
tergabung secara keseluruhan tidak dapat dilihat atau dianalisa secara molekul [8]. Sebagai konsekuensi bahwa fluida continuum adalah setiap sifat fluida (fluid property) diasumsikan sama pada setiap titik pada setiap bagian dan bersifat kontinyu terhadap posisi dan waktu [9]. Salah satu klasifikasi aliran fluida dapat dilihat pada Gambar 2.2.
Aliran Fluida
Inviscid µ=0
Viscous
Laminar
Compressible
Incompressible
Internal
Turbulent
External
Gambar 2.2 Klasifikasi aliran fluida [8].
2.1.1. Aliran Inviscid dan Viscous Aliran fluida berdasarkan viskositasnya dibagi menjadi dua bagian yaitu aliran inviscid dan viscous. Pada aliran inviscid efek dari viskositas (kekentalan) fluida diabaikan (μ = 0). Sebenarnya aliran fluida dengan viskositas sama dengan nol ini tidak ada. Namun untuk menyederhanakan analisa beberapa fenomena aliran mengabaikan viskositas boleh dilakukan. Untuk aliran fluida dimana viskositas sangat penting atau diperhatikan maka aliran itu disebut aliran viscous [9].
8
Gambar 2.3 Berbagai daerah aliran lapisan batas diatas plat rata [10]. Dari Gambar 2.3, terlihat mulai tepi depan plat itu terbentuk suatu daerah dimana pengaruh gaya viskos (Viscous force) makin meningkat. Gaya viskos ini biasanya dihubungkan dengan tegangan geser [10]. τ=μ
du dy
(2.1)
dimana : u = kecepatan dalam arah x
du = perbedaan kecepatan fluida pada arah y dy μ = viskositas fluida 2.1.2. Aliran Laminar dan Aliran Turbulent Berdasarkan struktur alirannya, aliran fluida dibedakan menjadi aliran laminar dan aliran turbulen. Untuk aliran laminar kecepatan pada suatu titik akan tetap
terhadap
waktu.
Sedangkan
aliran
turbulen
kecepatannya
akan
mengindikasikan suatu fluktuasi yang acak. Dalam aliran turbulen, profil kecepatan pada suatu titik dihasilkan dari gerak acak partikel fluida berdasarkan waktu dalam jarak dan arah. Jika kita mengambil kecepatan rata-rata terhadap waktu, maka kecepatan sesaat dapat dihitung dengan menambahkan kecepatan rata-rata dengan kecepatan fluktuasi. Secara grafis dapat ditunjukkan pada Gambar 2.4 [8].
9
(a) Laminar flow
(b) Steady turbulent flow
Gambar 2.4 Variasi kecepatan (satu dimensi) terhadap waktu [8].
2.1.3. Aliran Kompresibel dan Inkompresibel Aliran di mana perbedaan dalam massa jenis dapat diabaikan disebut inkompresibel. Ketika perbedaan massa jenis aliran yang tidak dapat diabaikan, aliran ini disebut kompresibel. Pada kenyataannya tidak ada fluida yang massa jenisnya konstan, tetapi ada beberapa masalah aliran fluida yang dapat disederhanakan dengan menganggap massa jenisnya konstan. Hal ini tidak mengurangi keakuratan solusi yang didapat. Parameter yang menjadi acuan utama untuk menentukan suatu aliran kompresibel atau tidak, dilihat dari nilai Mach Number (M), yang didefinisikan sebagai rasio antara kecepatan aliran lokal terhadap kecepatan suara lokal [8]. M=
V c
Dimana: V = kecepatan aliran (freestream) (m/s) c = kecepatan suara (m/s)
(2.2)
Pada saat M < 0.3, aliran tersebut dianggap sebagai aliran inkompresibel
2.1.4. Aliran Internal dan Eksternal Aliran yang dibatasi oleh suatu permukaan batas seperti pipa atau pembuluh disebut aliran internal. Contohnya seperti pada Gambar 2.5. Aliran mengalir pada benda yang terbenam di dalam fluida yang tak berbatas diistilahkan aliran eksternal.
Aliran internal dan eksternal keduanya dapat berupa aliran
laminer atau turbulen, kompresibel atau inkompresibel [7]. Contoh-contoh aliran 10
eksternal mencakup aliran udara pada pesawat terbang, mobil, gumpalan salju yang turun, atau aliran air disekitar kapal selam dan ikan. Aliran eksternal yang melibatkan udara sering disebut sebagai aerodinamika untuk menunjukkan arti penting dari aliran eksternal yang dihasilkan ketika sebuah objek seperti sebuah pesawat terbang menjelajah atmosfer [10].
Gambar 2.5 Aliran dalam pipa [8]. 2.1.5.
Aliran 2 fasa Aliran 2 fasa adalah aliran yang terjadi dalam suatu sistem yang
mengandung gas dan cairan. Aliran 2 fasa dapat dikelompokan menjadi aliran gas-cair dan cair-cair [9]. Kriteria dari 2 fasa yaitu perbedaan densitas dan viskositas. 2.2
Persamaan Dasar Aliran Fluida Persamaan dasar fluida pada hakikatnya dapat dikembangkan dari ketiga
prinsip dasar fisika yaitu kekekalan massa, kekekalan momentum, dan kekekalan energi. 2.2.1. Persamaan Kekekalan Massa Langkah pertama dalam penurunan persamaan kekekalan massa adalah menuliskan keseimbangan massa untuk elemen fluida. Kelajuan peningkatan massa dalam elemen fluida
=
Neto kelajuan aliran massa ke dalam elemen fluida
Kelajuan peningkatan massa dalam elemen fluida adalah
11
Untuk
mengevaluasi
properties
pada
masing-masing
enam
sisi
permukaan atur, kita menggunakan sebuah pengembangan deret Taylor kira-kira pada titik O. Sebagai contoh, pada sisi kanan
Abaikan suku-suku orde yang lebih tinggi, kita dapat menulis
Dan
Dimana
semua dievaluasi pada titik O. Suku-suku dapat
disamakan pada sisi kiri adalah
Selanjutnya kita perlu menerangkan kelajuan massa aliran melintasi sebuah bidang elemen yang diberikan oleh hasil dari densitas, luas dan komponen kecepatan normal terhadap bidang. Dari Gambar 2.6 dapat dilihat bahwa neto kelajuan aliran massa ke dalam elemen melewati boundarinya diberikan oleh
12
Aliran-aliran yang sejajar ke dalam elemen mengahasilkan sebuah peningkatan massa dalam elemen dan mempunyai sebuah tanda positif dan aliran-aliran yang meninggalkan elemen diberikan sebuah tanda negatif.
z y
x
Gambar 2.6 Massa mengalir kedalam dan keluar elemen fluida [11].
Kelajuan peningkatan massa ke dalam elemen (2.3) sekarang sama dengan neto kelajuan aliran massa ke dalam elemen fluida melintasi bidangnya (2.4). Semua suku-suku hasil akhir keseimbangan massa disusun pada sisi sebelah kiri dengan tanda yang sama dan pernyataan dibagi dengan Ini menghasikan
Karena operator vector, , dalam koordinat rectangular, diberikan oleh
13
maka
dan kekekalan massa boleh ditulis sebagai
Persamaan (2.6) adalah unsteady, kekekalan massa atau persamaan kontinuitas tiga dimensi pada sebuah titik dalam sebuah fluida kompresibel. Suku pertama pada sisi sebelah kiri kelajuan perubahan dalam waktu dari densitas (massa per satuan volume). Suku kedua menjelaskan neto aliran massa keluar dari elemen melintasi boudarinya dan disebut suku konvektif. Untuk fluida inkompresibel (misalnya liquid) densitas adalah konstan dan persamaan (2.6) menjadi [11].
atau dalam penjabarannya
2.2.2. Persamaan Kekekalan Momentum Hukum newton kedua menyatakan bahwa laju perubahan momentum partikel fluida sama dengan jumlah gaya gaya pada partikel.
Tingkat kenaikan momentum partikel fluda = jumlah gaya gaya pada partikel
Kelajuan peningkatan momentum x, y, dan z per satuan volume partikel fluida diberikan oleh
14
.
Kita membedakan dua tipe gaya pada partikel fluida: a) Gaya-gaya permukaan: i.
Gaya tekanan
ii.
Gaya viskos
b) Gaya-gaya badan: i.
Gaya gravitasi
ii.
Gaya sentrifugal
iii.
Gaya coriolis
iv.
Gaya elektromagnetik
Biasa dilakukan menyoroti kontribusi yang disebabkan gaya-gaya permukaan sebagai bagian tersendiri dalam persamaan momentum dan memasukkan gaya-gaya badan sebagai suku source.
Keadaan tegangan dari sebuah elemen fluida didefinisikan dalam suku suku tekanan dan sembilan komponen tegangan viskos ditunjukkan dalam Gambar 2.7. Tekanan, sebuah tekanan normal, di tandai oleh tegangan viskos ditandai oleh . Notasi akhiran yang biasa menandakan arah tegangan viskos. akhiran i dan j dalam
. Tegangan-
digunakan untuk menandakan bahwa
komponen tegangan bekerja dalam arah j pada sebuah permukaan normal kearah i.
15
z y x Gambar 2.7 Komponen tegangan pada tiga bidang elemen fluida [11]. Pertama kita mempertimbangkan komponen-komponen arah x dari gayagaya yang disebabkan komponen tegangan
,
dan
dan tekanan p
ditunjukkan dalam Gambar 2.8. Besarnya hasil akhir gaya dari sebuah tegangan permukaan adalah hasil dari tegangan dan luas. Gaya yang sejajar dengan arah dari sebuah axis koordinat menjadi sebuah tanda positif dan itu dalam arah yang berlawanan dengan sebuah tanda negatif. Neto gaya dalam arah x adalah jumlah komponen – komponen gaya bekerja dalam arah itu pada elemen fluida.
Gambar 2.8. Komponen tegangan dalam arah x [11]. 16
Pada pasangan sisi (timur, barat) kita mempunyai
Gaya total dalam arah x pada pasangan sisi (utara, selatan) adalah
Akhirnya gaya total dalam arah x pada sisi bawah dan atas diberikan oleh
Total gaya per satuan volume pada fluida yang disebabkan tegangan-tegangan permukaan ini adalah sama dengan jumlah dari (2.10a), (2.10b) dan (2.10c) dibagi oleh volume
Tanpa
:
mempertimbangkan
gaya
badan
dalam
lebih
detail
lagi
pengaruhnya secara menyeluruh dapat dimasukkan dengan mendefinisikan sebuah source
dari momentum x per satuan volume per satuan waktu. Komponen x persamaan momentum ditemukan dengan menetapkan
perubahan momentum x dari partikel fluida (2.9) sama dengan total gaya dalam arah x pada elemen yang disebabkan tegangan permukaan (2.11) ditambah kelajuan peningkatan momentum x disebabkan source:
17
Tidak sulit untuk membuktikan bahwa komponen y persamaan momentum diberikan oleh
dan juga komponen z persamaan momentum oleh
Tanda yang dihubungkan dengan tekanan berlawanan terhadap yang dihubungkan dengan tegangan viskos normal, karena konvensi tanda yang biasa digunakan sebuah tegangan tarik yang menjadi tegangan normal positif sehingga tekanannya, yang mana dengan mendefinisikan sebuah tegangan normal tekan, mempunyai tanda negatif [11].
2.2.3. Persamaan Navier–Stokes untuk Sebuah Fluida Newtonian.
Dalam sebuah fluida Newtonian tegangan-tegangan viskos adalah sebanding dengan kelajuan deformasi. Bentuk tiga dimensi hukum Newton viskositas untuk aliran kompresibel melibatkan konstanta kesebandingan: viskositas dinamik (pertama),
, menghubungkan tegangan-tegangan dengan
deformasi linier, viskositas kedua, , menghubungkan tegangan-tegangan dengan deformasi volumetrik. Sembilan komponen tegangan viskos, dari yang enam adalah berdiri sendiri, adalah
18
Tidak banyak diketahui tentang viskositas kedua , karena pengaruhnya kecil dalam praktek. Untuk gas suatu pendekatan yang baik dapat diperoleh dengan
mengambil
nilai
(Schlichting,
1979).
inkompresibel jadi persamaan konservasi massa adalah
Liquid
adalah
dan tegangan
viskos hanya dua kali kelajuan deformasi linier viskositas dinamamiknya. Substitusi tegangan-tegangan geser di atas (2.13) ke dalam (2.12a-c) menghasilkan persamaan yang disebut persamaan Navier-Stokes:
Sering berguna untuk menyusun kembali suku-suku tegangan viskos sebagai berikut:
Tegangan-tegangan viskos dalam komponen-komponen y dan z persamaan dapat dibuat kembali dalam sebuah cara yang sama. Persamaan Navier-Stokes dapat ditulis kembali dalam bentuk yang lebih berguna untuk pengembangan metode volume hingga:
19
Jika kita menggunakan model Newtonian untuk tegangan-tegangan viskos dalam persamaan energi dalam kita peroleh setelah beberapa disusun ulang
Sebuah pengaruh yang disebabkan dalam persamaan energi dalam ini dijelaskan oleh fungsi dissipasi
. Dapat ditunjukkan menjadi sama dengan
Untuk menyimpulkan temuan-temuan hingga kini kita kutip dalam Tabel 2.1 bentuk perbedaaan atau kekekalan system persamaan yang membangun aliran fluida tiga dimensi tergantung waktu dan perpindahan panas sebuah fluida Newtonian kompresibel.
Tabel 2.1 Persamaan Pembangun Aliran Fluida Newtonian Kompresibel [11]. Massa
(2.7)
Momentum x
(2.15a)
Momentum y
(2.15b)
Momentum z
(2.15c)
20
2.3.
Bilangan Reynolds Bilangan Reynolds merupakan parameter tak-berdimmensi yang sangat
terkenal dalam ilmu mekanika fluida. Nama ini diberikan sebagai penghargaan bagi Osborne Reynolds (1842-1912), insinyur dari inggris yang pertama kali mendemontrasikan bahwa kombinasi dari variabel-variabel dapat digunakan sebagai suatu patokan untuk membedakan aliran laminar dengan aliran turbulen. Pada kebanyakan soal aliran fluida, akan kita dapati diameter pipa, D, dan kecepatan,
, demikian juga kerapatan fluida, , dan viskositas, , merupakan
variabel-variabel yang relevan dalam sebuah soal. Dengan variabel tersebut, Bilangan Reynolds
Secara alamiah muncul dari suatu analisis dimensional. Bilangan Reynolds adalah ukuran dari rasio gaya inersia pada suatu elemen fluida terhadap gaya viskositas elemen [8]. 2.4
Coalescer Proses pemisahan coalescer oil water separators yaitu air dan campuran
minyak masuk separator dan menyebar horizontal, didistribusikan melalui energi dan turbulensi akan menyebar ke perangkatnya. Campuran memasuki media Flopak yang ada pada laminer dan arus sinusoidal dibentuk dan minyak pada permukaan media. Minyak akan menyatu menjadi tetesan yang lebih besar, naik ke atas melalui corrugations, di mana minyak dan air dapat terpisah dan naik ke permukaan air.[3]
21
Gambar 2.9 Coalecer [5].
Ada 2 bahan yang digunakan untuk plat coalescer yaitu plat polypropylene dan polietilena densitas tinggi. Polypropylene lebih oleophilic (karakteristik zat nya lebih menyerupai minyak dari pada air ) dan menjadi bahan yang digunakan dalam produksi minyak secara umum. Beberapa orang industri percaya bahwa polypropylene terlalu oleophilic ( karakteristik zat nya lebih menyerupai minyak dari pada air) dan tidak memungkinkan minyak untuk berpindah ke permukaan air. Bahan Polietilena umumnya digunakan sebagai material aleternatif untuk perbandingan [4].
Gambar 2.10 Konfigurasi Coalescer [4]. 22
Gambar 2.11 Oleophilic Material [4]. Banyak
faktor yang harus dipertimbangkan dalam membuat atau
mendesain sistem pemisahan air dan minyak/coalescer yaitu 1. Kondisi dan laju aliran 2. Diperlukan derajat pemisahan 3. Jumlah minyak didalam air inlet 4. Emulsifikasi minyak 5. Metode perbaikan pembuangan minyak. Manfaat dari alat coalescer atau pemisahan air dan minyak dari sudut pandang ekonomi dan lingkungan adalah 1. Ekonomis dan bagus dalam pengoperasiannya 2. Rendah biaya operasi dan pemeliharaan 3. Prinsip kerjanya berdasarkan gravitasi jadi tidak diperlukan alat tambahan 4. Minyak dapat didaur ulang kembali 5. Kemampuan untuk menangani kenaikan debit minyak dan air [12].
2.5
Design coalescer Design yang sudah di lakukan atau yang sudah ada yaitu: 1. Penggabungan plat pemisah. a. Pelat pemisah Pelat pemisah ini telah di gunakan dan berhasil selama bertahun-tahun.
Sistem ini biasanya terbuat dari fiberglass pelat bergelombang yang menjadi satu dan didalamnya frame baja dan stainless steel. Minyak akan naik dan memasuki 23
sistem sampai minyak tersebut mencapai permukaan di atas. Keuntungan dari sistem ini adalah mencakup peningkatan efisiensi minyak. Pemisah jenis ini umumnya dapat di konfigurasikan dalam beberapa orientasi dengan arus tegak lurus terhadap pelat seperi gambar yang di tunjukkan pada gambar 2.12, sejajar dengan pelat secara upflow dan sejajar dengan pelat secara downflow.
Gambar 2.12 Pelat pemisah [13]. b. Pelat pemisah yang datar bergelombang Pelat pemisah datar bergelomang sering menggunakan Horizontal polypropylene oleophilic ( karakteristik zat nya lebih menyerupai minyak dari pada air ) di tumpuk satu di atas yang lain dalam tumpukan vertikal dan diikatkan ke batang dengan kabel. Gambar 2.13 mengilustrasikan gambar sistem pelat datar bergelombang.
24
a. Tampak dalam
b. Tampak samping Gambar 2.13 Pemisah pelat datar bergelombang [13].
c. Sudut pemisah 2.
Penggabungan tabung pemisah Penggabungan tabung pemisah memanfaatkan tabung plastik berlubang
untuk pemisahan di gunakan untuk pemisahan minyak dan air. Pemisah ini biasanya di buat dengan tabung di posisi vertikal tapi ada pula yang di buat
25
dengan posisi horizontal. Keuntungan dari penggunan tipe pemisah ini adalah biaya rendah.
Gambar 2.14 Penggabungan tabung pemisah [13]. 3. Pemisah jenis
Jenis coalescer ini sering digunakan dalam hubungannya dengan gravitasi pemisahan. Coalecser jenis ini terbuat dari bahan stainless steel dan polypropylene. Sistem jenis ini sangat efisien, tapi penggabungan fluida padat bisa mengalami penyumbatan. Jenis ini sering digunakan sebagai tahap kedua setelah tipe plat tahap pertama pemisahan [13].
Gambar 2.15 Pemisah jenis [13]. 26
BAB III METODE VOLUME HINGGA FLUENT 6.3.26 dan PEMODELAN SIMULASI Computational Fluid Dynamics (CFD) adalah ilmu yang mempelajari cara memprediksi aliran fluida, perpindahan panas, reaksi kimia, dan fenomena lainnya dengan menyelesaikan persamaan persamaan matematika ( model matematika ) [14]. CFD memprediksikan aliran berdasarkan [14] : 1. Model matematika ( persamaan diferensial parsial ). Khususnya memecahkan persamaan Navier-Stokes 2. Metode numerik ( teknik solusi dan diskritisasi ) 3. Tools perangkat lunak ( solvers, tools pre- dan postprocessing )
3.1. FLUENT FLUENT adalah program komputer yang memodelkan aliran fluida dan perpindahan panas dalam geometri yang kompleks. FLUENT versi 6.3.26 merupakan salah satu jenis program CFD (Computational Fluid Dynamics) yang menggunakan metode diskritisasi volume hingga. FLUENT memiliki fleksibilitas mesh, sehingga kasus-kasus aliran fluida yang memiliki mesh tidak terstruktur akibat geometri benda yang rumit dapat diselesikan dengan mudah. Selain itu, FLUENT memungkinkan untuk penggenerasian mesh lebih halus atau lebih besar dari mesh yang sudah ada berdasarkan pemilihan solusi aliran [9]. Secara lengkap langkah-langkah FLUENT dalam menyelesaikan suatu simulasi adalah sebagai berikut [14]: 1. Membuat geometri dan mesh pada model. 2. Memilih solver yang tepat untuk model tersebut (2D atau 3D). 3. Mengimpor mesh model (grid). 4. Melakukan pemeriksaan pada mesh model. 5. Memilih formulasi solver. 27
6. Memilih persamaan dasar yang akan dipakai dalam analisa. 7. Menetukan sifat material yang akan dipakai. 8. Menentukan kondisi batas. 9. Mengatur parameter kontrol solusi. 10. Initialize the flow field. 11. Melakukan perhitungan/iterasi. 12. Menyimpan hasil iterasi. 13. Jika diperlukan, memperhalus grid kemudian melakukan iterasi ulang.
28
Mulai
Pemodelan dengan Gambit dan Meshing
Pengecekan Meshing
Tidak Mesh baik (sudut kemiringan mesh < 0.85)?
Data Sifat Fisik
Ya
Penentuan Kondisi batas
Proses Numerik
Tidak
Iterasi Konvergen (kriteria
konvergensi 10-6)
Plot Hasil
Selesai Gambar 3.1 Diagram Alir Prosedur Simulasi FLUENT [14].
29
3.2. SKEMA NUMERIK FLUENT memberikan dua pilihan metode numerik, yaitu metode segregated dan coupled. Kedua metode tersebut dapat digunakan untuk memecahkan persamaan integral kekekalan momentum, massa, dan energy (governing integral equation), serta besaran skalar lainnya seperti turbulensi. Dalam proses pemecahan masalah, baik metode segregated dan coupled memiliki persamaan yaitu menggunakan teknik kontrol volume. Teknik kontrol volume sendiri terdiri dari: Pembagian daerah asal (domain) ke dalam kontrol volume diskrit dengan menggunakan grid komputasi. Integrasi persamaan umum kontrol volume untuk membuat persamaan aljabar dari variabel tak-bebas yang berlainan (discrete dependent variables) seperti kecepatan , tekanan, suhu, dan sebagainya Linearisasi persamaan dan solusi diskritisasi dari resultan sistem persamaan linear untuk menghasilkan nilai taksiran variabel tak-bebas. Pada dasarnya metode segregated dan coupled memiliki persamaan dalam proses diskritisasi yaitu volume berhingga (finite volume), tetapi memiliki perbedaan pada cara pendekatan yang digunakan untuk melinearisasi dan memecahkan suatu permasalahan [9]. 3.2.1. Metode Solusi Segregated Metode ini menyelesaikan persamaan kekekalan massa, momentum, dan energi secara bertahap atau terpisah satu sama lain. Karena persamaan kekekalan massa, momentum, dan energi merupakan persamaan non-linear, beberapa iterasi harus dilakukan secara berulang-ulang sebelum solusi yang konvergen diperoleh. Dalam iterasi terdiri dari beberapa langkah, yaitu [9]: 1. Sifat-sifat fluida diperbarui berdasarkan solusi yang telah dilakukan. Untuk perhitungan awal, sifat-sifat fluida diperbaharui
berdasarkan solusi awal
(initialized solution).
30
2. Persamaan momentum u, v, dan w dipecahkan dengan menggunakan nilainilai tekanan dan fluks massa permukaan, supaya medan kecepatan diperbaharui. 3. Karena kecepatan yang diperoleh dalam tahap yang pertama tidak mungkin memenuhi persamaan kontinuitas secara lokal, persamaan “Poisson-type” untuk koreksi tekanan diturunkan dari persamaan kontinuitas dan persamaan momentum linear. Persamaan koreksi tekanan ini kemudian dipecahkan untuk memperoleh koreksi yang dibutuhkan untuk medan tekanan dan kecepatan serta fluks massa permukaan sampai kontinuitas dipenuhi. 4. Menyelesaikan persamaan-persamaan untuk besaran skalar seperti turbulensi, energi, radiasi dengan menggunakan nilai-nilai variabel lain yang telah diperbaharui. 5. Mengecek konvergensi persamaan. 3.2.2. Metode Solusi Coupled Metode ini menyelesaikan persamaan kekekalan massa, momentum, dan energi secara serempak atau bersamaan (simultaneously). Karena persamaan kekekalan massa, momentum, dan energi merupakan persamaan non-linear, beberapa iterasi harus dilakukan
secara berulang-ulang sebelum solusi yang konvergen
diperoleh. Dalam iterasi terdiri dari beberapa langkah, yaitu [9]: 1. Sifat-sifat fluida diperbaharui berdasarkan solusi yang telah dilakukan. Untuk perhitungan awal, sifat-sifat fluida diperbaharui
berdasarkan solusi awal
(initialized solution). 2. Persamaan kontinuitas, momentum, dan energi jika ada serta besaran-besaran tertentu lainnya dipecahkan secara serempak. 3. Jika ada, persamaan-persamaan skalar seperti turbulensi
dan radiasi
dipecahkan dengan menggunakan nilai yang diperbaharui sebelumnya berdasarkan variable yang lain. 4. Mengecek konvergensi persamaan.
31
3.3. DISKRITISASI (DISCRETIZATION) FLUENT menggunakan suatu teknik berbasis volume kendali untuk mengubah bentuk persamaan umum (governing equation) ke bentuk persamaan aljabar (algebraic equation) agar dapat dipecahkan secara numerik. Teknik kontrol volume ini intinya adalah pengintegralan persamaan diferensial umum untuk setiap olume kendali, sehingga menghasilkan suatu persamaan diskrit yang menetapkan setiap besaran pada suatu basis volume kendali. Diskritisasi persamaan umum dapat diilustrasikan dengan menyatakan persamaan kekekalan kondisi-steady untuk transport suatu besaran skalar. Hal ini ditunjukkan dengan Persamaan 3.1 yang ditulis dalam bentuk integral untuk volume kendali sembarang. Persamaan 3.1 diterapkan untuk tiap volume kendali atau sel dalam daerah asal komputasi (domain). sebagai berikut :
∫ ρ φ v. dA = ∫ Γφ ∇φ . dA + ∫ Sφ dV
(3.1)
V
Suku sumber
Suku konveksi
(source term)
Suku difusi dimana = kerapatan fluida = vektor kecepatan (
dalam 2D)
= vektor area permukaan = koefisien difusi untuk = gradien = sumber
(
dalam 2D)
tiap satuan volume.
Persamaan (3.1) diterapkan untuk tiap volume kendali atau sel dalam daerah asal komputasi (domain). Diskretisasi persamaan (3.1) pada sel tertentu diberikan pada persamaan berikut :
32
(3.2) dimana : = jumlah sisi = nilai
yang dikonveksikan melalui sisi f
= fluks massa yang melalui sisi = luas sisi f, A (= Ax iˆ + Ay ˆj dalam 2D = jumlah ∇φ yang tegak lurus terhadap sisi f = volume sel
Gambar 3.2 mengilustrasikan diskretisasi persamaan perpindahan scalar dengan teknik volume kendali.
Gambar 3.2 Volume Kendali Digunakan Sebagai Ilustrasi Diskretisasi Persamaan Transport Skalar [9].
33
Untuk penggunaan model sel 2D quadrilateral ditunjukkan pada Gambar 3.3 yang merupakan suatu contoh volume kendali.
co A c1 Gambar 3.3 Volume Kendali Digunakan Sebagai Ilustrasi Diskretisasi Persamaan Transport Skalar pada model sel 2D quadrilateral [9]. FLUENT menyimpan nilai-nilai diskrit skalar
pada pusat-pusat sel (c0
dan c1 pada Gambar 3.2 dan Gambar 3.3). Meskipun demikian, nilai-nilai sisi diperlukan untuk suku konveksi dalam Persamaan 3.2 dan harus diinterpolasi dari nilai-nilai
pusat
sel.
Hal
ini
diselesaikan
upwind.Upwinding berarti bahwa nilai sisi
dengan
menggunakan
skema
diturunkan dari besaran-besaran hulu
atau “upwind”, relatif terhadap arah kecepatan tegak lurus ,
dalam Persamaan 3.2
[9]. 3.3.1. First-Order Upwind Ketika menginginkan keakuratan (accuracy) orde-pertama, besaran-besaran sisi sel ditentukan dengan cara mengasumsikan bahwa nilai-nilai pusat-sel pada beberapa variabel medan menggambarkan nilai rata-rata-sel dan berlaku untuk seluruh sel; besaran-besaran sisi identik dengan besaran-besaran sel. Oleh karena itu, ketika first-order upwind dipilih, nilai sisi
diatur sama dengan nilai-pusat
pada
sel upstream[9]. 3.3.2. Second-Order Upwind Scheme Ketika menginginkan keakuratan (accuracy) orde-kedua, besaran-besaran pada sisi sel ditentukan dengan menggunakan suatu pendekatan rekontruksi linear 34
multidimensi. Dalam pendekatan ini, keakuratan orde yang lebih tinggi diperoleh pada sisi-sisi sel melalui ekspansi deret Taylor berdasarkan solusi pusat sel di sekitar sentroid sel. Oleh karena itu, saat second-order upwinding dipilih, nilai sisi dihitung dengan menggunakan persamaan sebagai berikut : (3.3) dimana
dan
merupakan nilai pusat-sel dan gradiennya dalam sel upstream, dan
merupakan vektor perpindahan dari sentroid sel upstream ke sentroid sisi. Formulasi ini membutuhkan penentuan gradient
di setiap sel. Gradien ini dihitung
dengan menggunakan teorima divergensi,dan dalam bentuk diskret ditulis sebagai : (3.4) Oleh karena itu nilai face
dihitung dengan merata-ratakan
dari dua sel
yang berdekatan dengan sisi (face) [9].
3.3.3. Bentuk Linearisasi Persamaan Diskrit Persamaan transport skalar terdiskretisasi (Persamaan 3.2) mengandung variabel skalar
yang tidak diketahui pada pusat sel sebagaimana nilai-nilai yang
tidak diketahui dalam sel-sel tetangga yang ada disekelilingnya. Persamaan ini pada umumnya akan menjadi tidak linear akibat adanya variabel-variabel ini. Bentuk linear Persamaan 3.2 dapat ditulis sebagai berikut : (3.5) dimana subkrip
menunjukkan sel tetangga, dan
koefisien linear untuk
dan
dan
adalah koefisien-
[9].
3.3.4. Under-Relaxation Ketidaklinearan mengakibatkan perubahan
persamaan
yang
dipecahkan
oleh
FLUENT
dapat
yang dihasilkan pada tiap iterasi menjadi tidak teratur.
Tipikal dari adanya under-relaxation adalah mengurangi perubahan
yang
dihasilkan dari setiap iterasi. Dalam bentuk yang sederhana, nilai variable
yang
35
baru dalam sebual sel tergantung kepada nilai
sebelumnya,
, perubahan ∆ φ
yang dihitung, dan faktor under relaxation, α , dirumuskan sebagai berikut :[9] (3.6) 3.4. DISKRITISASI COUPLED SOLVER Formulasi coupled solver pada FLUENT memecahkan persamaan kontinuitas, momentum, energi (jika diperlukan) dan persamaan species transport secara bersamaan sebagai sebuah set, atau vektor, dari persamaan-persamaan tersebut. Persamaan-persamaan untuk besaran-besaran tambahan akan diselesaikan secara bertahap (yaitu segregated dari satu yang lain dan dari set coupled).
3.4.1. Aturan Persamaan-Persamaan dalam Bentuk Vektor Sistem dari persamaan-persamaan untuk fluida single-component, ditulis untuk menguraikan nilai aliran rata-rata, dipindah bentuk integral, bentuk kartesius untuk volume kendali V dengan area permukaan diferensial dA sebagai berikut: (3.7) dimana vektor W, F, dan G didefinisikan sebagai
ρv ρ ρvu ρu W = ρv , F = ρvu ρvu ρw ρvE ρE
+ + + +
0 τ ρiˆ xi ˆ ρj , G = τ yi τ ρkˆ zi ρv τ ij v j
+ q
(3.8)
dan vektor H berisi batas seperti body forces dan energy sources. Dimana ρ, υ, E, dan p adalah densitas, kecepatan, energi total per satuan massa, dan tekanan fluida, secara berturut-turut. τ adalah tekanan yang melekat, dan q adalah fluks panas. Energi total E dihubungkan dengan entalpi total H oleh
dimana (3.9) 36
Persamaan Navier-Stokes dinyatakan dalam persamaan (3.7) (secara numerik) menjadi sangat berat pada Mach number rendah dalam kaitannya dengan perbedaan antara kecepatan fluida v dan acoustic speed c (kecepatan suara). Ini juga dibenarkan untuk aliran inkompresibel, dengan mengabaikan kecepatan fluida, karena gelombang akustik sangat cepat pada fluida inkompresibel (kecepatan suara sangat cepat). Sulitnya perhitungan pada persamaan-persamaan dibawah kondisi tersebut mengakibatkan pada nilai konvergensi yang tidak baik. Kesulitan ini diatasi pada FLUENT’s coupled solver dengan memanfaatkan suatu teknik yang disebut preconditioning (time-derivative)[9]. 3.4.2. Preconditioning Time-derivative preconditioning mengubah bentuk time-derivative dari persamaan (3.7) dengan cara pre-multiplying persamaan tersebut dengan matrik preconditioning. Asal mula matrik preconditioning dimulai dengan mengganti variabel yang tidak bebas pada persamaan (3.7) dari jumlah W untuk variabel primitif Q menggunakan kaidah-rantai sebagai berikut: (3.10) dimana Q adalah vektor {p, u, v, w, T}T dan Jacobian ∂W /∂Q ditunjukkan oleh
ρp ρ u p ∂W = ρ pv ∂Q ρ pw ρ p H
ρT ρ ρT 0 ρ ρT 0 0 ρ ρT ρu ρv ρw ρ T H 0
− δ
0 0
0 0 0
+ ρC p
(3.11)
dimana
ρp =
∂ρ ∂ρ , ρT = ∂p T ∂T
p
dan δ = 1 untuk gas ideal dan δ = 0 untuk fluida inkompresibel.
37
Kita mempersiapkan sistem dengan mengganti matrik Jacobian ∂W /∂Q dengan matrik preconditioning Г jadi sistem preconditioned berubah bentuk menjadi
(3.12) dimana
Θ Θu Γ = Θv Θw ΘH
ρT ρ ρT u 0 ρ ρT u 0 0 ρ ρT u ρu ρv ρw ρ T H 0
− δ
0 0
0 0 0
+ ρC p
(3.13)
Parameter Ө ditunjukkan oleh (3.14)
Hasil eigenvalues dari sistem preconditioned (persamaan 12) adalah
(3.15) Dimana
Untuk gas ideal, β = (γRT)-1 = 1/c2. Jadi, ketika Ur = c (pada kecepatan sonik dan diatasnya), α = 0 dan eigenvalues dari sistem preconditioned memerlukan bentuk 38
dasar, u ± c. Pada kecepatan yang rendah, bagaimanapun juga, Ur → 0, α → 1/2 dan semua eigenvalues menjadi sama seperti u [9].
Flux-Difference Splitting Inviscid flux dari vektor F ditunjukkan pada persamaan (3.12) dievaluasi dengan standar upwind, flux-difference splitting. Dengan pemisahan F menjadi bagianbagian, dimana masing-masing bagian berisi informasi perpindahan arah tertentu (yaitu, informasi karakteristik), dan upwind membedakan fluks yang dipisah pada cara yang tetap dengan eigenvalues yang sesuai, diperoleh sebagai berikut untuk fluks pada masing-masing bagian: F
(3.16)
Dimana δQ adalah perbedaan dari QR - QL. Fluks FR = F(QR) dan FL = F(QL) dihitung menggunakan solusi vektor QR dan QL pada sisi kanan dan kiri dari bagian tersebut. Matrik |Â| didefinisi sebagai (3.17)
dimana Λ adalah matrik diagonal dari eigenvalues dan M adalah matrik yang berhubungan dengan diagonal Г -1A, dimana A adalah fluks Jacobian ∂F/∂Q. Dalam bentuknya sekarang, persamaan (3.16) dapat dipandang sebagai orde kedua perbedaan ditambah matrik pemisah. Batas penambahan matrik pemisah tidak hanya untuk menghasilkan upwinding dari variabel convected, dan tekanan dan kecepatan fluks pada aliran supersonik, tetapi juga gabungan tekanan-kecepatan diperlukan untuk kestabilan dan konvergensi yang efisien pada kecepatan rendah dan aliran inkompresibel [9].
39
3.5 MODEL MULTIFASA Pada aliran multifasa, suatu fasa dapat didefinisikan sebagai suatu tingkat material yang mempunyai respon inersia terhadap aliran dan berinteraksi dengan aliran dan medan potensial yang ada. Misalnya, material partikel padat yang sama tetapi berbeda ukuran dapat diperlukan sebagai fasa yang berbeda karena sekelompok patikel dengan ukuran yang sama akan mempunyai respon dinamik yang sama terhadap medan aliran [14].
FLUENT menyediakan beberapa model untuk aliran multifasa, yaitu: 1. Volume of Fluid ( VOF ) 2. Mixture 3. Eulerian 3.6 MODEL TURBULEN (TURBULENCE MODELING) Aliran turbulen adalah suatu karakteristik yang terjadi karena adanya peningkatan kecepatan aliran. Peningkatan ini mengakibatkan perubahan momentum, energi, dan massa tentunya. Karena terlalu mahalnya untuk melakukan analisa secara langsung dari aliran turbulen yang memiliki skala kecil dengan frekuensi yang tinggi, maka diperlukan suatu manipulasi agar menjadi lebih mudah dan murah. Salah satunya adalah dengan permodelan turbulen (turbulence model). Meskipun demikian, modifikasi persamaan yang meliputi penambahan variabel yang tidak diketahui, dan permodelan turbulen perlu untuk menentukan variabel yang diketahui [9]. FLUENT sendiri menyediakan beberapa permodelan, diantaranya adalah k-ε dan k-ω. 3.6.1. Permodelan k-epsilon (k-ε) 3.6.1.1 Standard Model ini merupakan model turbulensi semi empiris yang lengkap. Walaupun masih sederhana, memungkinkan untuk dua persamaan yaitu kecepatan turbulen (turbulent velocity) dan skala panjang (length scale) ditentukan secara bebas 40
(independent). Model ini dikembangkan oleh Jones dan Launder. Kestabilan, ekonomis (dari segi komputansi), dan akurasi yang cukup memadai membuat model ini sering digunakan dalam simulasi fluida dan perpindahan panas [9]. 3.6.1.2 RNG Model ini diturunkan dengan menggunakan metode statistik yang teliti (teori renormalisasi kelompok). Model ini merupakan perbaikan dari metode kepsilon standard, jadi bentuk persamaan yang digunakan sama. Perbaikan yang dimaksud meliputi [14] : Model RNG memiliki besaran tambahan pada persamaan laju disipasi (epsilon), sehingga mampu meningkatkan akurasi untuk aliran yang terhalang secara tiba-tiba. Efek putaran pada turbulensi juga telah disediakan, sehingga meningkatkan akurasi untuk jenis aliran yang berputar (swirl flow). Menyediakan formulasi analitis untuk bilangan Prandtl turbulen, sementara model k-epsilon standard menggunakan nilai bilangan Prandtl yang ditentukan pengguna (kostan). Model RNG menyediakan formulasi untuk bilangan Reynold rendah, sedang model standard merupakan model untuk Reynold tinggi
3.6.1.3 Realizable Merupakan model pengembangan yang relatif baru dan berbeda dengan model k-epsilon standar dalam dua hal, yaitu: Terdapat formulasi baru untuk memodelkan viskositas turbulen. Sebuah persamaan untuk dissipasi, ε, telah diturunkan dari persamaan yang digunakan untuk menghitung fluktuasi vortisitas rata-rata. Istilah realizable memiliki arti bahwa model tersebut memenuhi beberapa batasan matematis pada bilangan Reynold, konsisten dengan bentuk fisik aliran turbulen. Kelebihannya adalah lebih akurat untuk memprediksi laju penyebaran fluida dari pancaran jet/nosel. Model ini memberikan performa yang bagus untuk aliran 41
yang melibatkan putaran, lapisan batas yang memiliki gradien tekanan yang besar, separasi, dan resirkulasi. Salah satu keterbatasan model realizable k-epsilon adalah terbentuknya viskositas turbulen non-fisik pada kasus dimana domain perhitungan mengandung zona fluida yang diam dan berputar ( multiple reference frame, sliding mesh ). Oleh karena itu, penggunaan model ini pada kasus multiple reference frame dan sliding mesh harus lebih hati-hati [9]. 3.6.2. Permodelan k-omega (k-ω) 3.6.2.1.Standard Model yang terdapat dalam FLUENT merupakan model berdasarkan Wilcox k-omega yang memasukkan beberapa modifikasi untuk menghitung efek aliran pada bilangan Rynold rendah, kompresibilitas, dan penyebaran aliran geser (shear flow). Selain itu, model ini juga mampu diaplikasikan untuk aliran dalam saluran maupun aliran bebas geseran (free shear flow) [14]. 3.6.2.2. SST Model ini dikembangkan oleh Menter untuk memadukan formulasi model komega standard yang stabil dan akurat pada daerah dekat ke dinding dengan model kepsilon yang memiliki kelebihan pada aliran free stream. Model ini mirip dengan komega standard dengan memiliki beberapa perbaikan, yaitu: Model k-omega standard dan k-epsilon yang telah diubah dikalikan dengan suatu fungsi pencampuran dan kedua model digunakan bersama-sama, sehingga lebih akurat untuk daerah dekat dinding maupun untuk aliran yang jauh dari dinding dan free stream flow. Definisi viskositas turbulen dimodifikasi untuk menghitung perubahan tegangan geser turbulen. Konstanta model berbeda dengan model k-omega standar. Melibatkan sebuah besaran dari penurunan damped cross diffusion pada persamaan omega [14]. 42
3.7. ADAPSI GRID FLUENT menyediakan fitur untuk memperbaiki mesh yang telah dibuat, perbaikan mesh ini dapat berupa penghalusan dan atau pengkasaran grid berdasarkan geometri dan solusi data numerik.
Sebagai tambahan, FLUENT memberikan
peralatan untuk membuat dan melihat penambahan dan pengurangan adapsi. Dua keuntungan dari mesh tidak terstruktur dalam FLUENT yaitu: •
Lebih cepat dalam menyusun grid
•
Kemampuan adapsi untuk memperbaiki mesh Dengan menggunakan adapsi, dapat menambah jumlah sel di bagian mesh
yang diinginkan, sehingga menghasilkan hasil yang lebih baik. Penggunaan adapsi dengan benar maka hasil dari mesh akan optimal untuk solusi aliran karena solusi digunakan untuk menentukan dimana grid harus ditambahkan. Dengan kata lain daya komputasi tidak akan terbuang percuma oleh sel-sel yang tidak diperlukan, sama seperti yang terjadi pada pendekatan grid tersusun. 3.7.1 Proses Adapsi Proses adapsi telah dipisahkan dalam dua proses yang terpisah. Pertama sel tunggal ditandai untuk penghalusan atau pengkasaran tergantung dari fungsi adapsi dimana dibuat dari geometri dan atau dari data solusi.
Berikutnya sel tersebut
dihaluskan atau dikasarkan tergantung dari tanda adapsi tersebut. Keuntungan utama yang didapat dari pendekatan ini adalah kemampuan untuk membuat fungsi adapsi yang baik dan bereksperimen dengan fungsi adapsi yang bervariasi tanpa merubah mesh yang telah ada. •
Hanging Node Adaption Grid-grid yang dihasilkan oleh hanging node adaption dikarakteristikkan
oleh titik-titik di atas garis dan dan bidang yang tidak bertemu dipuncak dari semua pembagian sel-sel garis dan bidang seperti digambarkan pada Gambar 3.4.
43
Gambar 3.4 Contoh Hanging Node [9]. Hanging node adaption menyediakan kemampuan untuk beroperasi diatas grid-grid dengan variasi bentuk sel termasuk grid-grid hybrid. Namun meskipun skema hanging node menyediakan fleksibilitas grid dengan signifikan tetapi membutuhakan memori tambahan untuk menjaga hirarki grid yang digunakan untuk rendering dan operasi adapsi grid.
•
Conformal Adaption Proses conformal adaption tidak menciptakan hanging nodes, tetapi semua
pembagian sel-sel garis atau bidang termasuk semua titik di entity tersebut. Proses penghalusan (refinement conformal) menambah titik-titik di garis dan pengkasaran (coarsening conformal) menghilangkan titik-titik dan menyegitigakan kembali ruang yang dihasilkan. Gambar 3.5 dan Gambar 3.6 memperlihatkan hasil dari adapsi conformal.
44
Gambar 3.5 Hasil Conformal Refinement.
Gambar 3.6 Conformal Coarsening dengan Menghilangkan Titik dan Menyegitigakan Kembali [9]. 3.7.2 Adapsi
dan
FLUENT menyediakan tiga pilihan berbeda untuk model turbulensi di dekat dinding (near-wall modeling of turbulance) yaitu standard wall functions, nonequilibrium wall functions, dan the two-layer model. Ada syarat-syarat mesh tertentu untuk tiap pilihan model turbulensi di dekat dinding ini. Tetapi ada kesulitan untuk mengukur syarat-syarat resolusi dekat dinding,
dan
adapsi disediakan 45
untuk memungkinkan ketepatan memperhalus atau mengkasarkan mesh sepanjang dinding selama proses solusi. Hasil dari adapsi menggunakan
dan
dapat dilihat
pada Gambar 3.11. Bila two-layer model digunakan sebagai model turbulensi di dekat dinding maka dapat menggunakan dapat menggunakan
atau
sebagai adapsinya, tetapi bila wall function
.
Gambar 3.7 Hasil Meshing Menggunakan Adapsi Pendekatan untuk menghitung
dan
[9].
untuk batas sel-sel di zona
dinding spesifik, tentukan nilai minimum dan maksimum yang diijinkan ,
atau
dan tandai dan atau adapsi sel-sel yang tepat.
atau
Sel-sel dengan nilai
dibawah batas yang diijinkan akan ditandai untuk dikasarkan dan sel-sel dengan nilai atau
diatas batas maksimum yang diijinkan akan ditandai untuk dihaluskan [9].
46
3.8. Pemodelan Bentuk dan Simulasi Numerik pada Fluent 6.3.26 3.8.1 Pemodelan Bentuk Sebuah coalescer dengan lebar 1.6256 meter, panjang 2.3622 meter, dan tinggi 2.0066 meter dan didalam coalescer terdapat sekat dengan jarak 1 inch, 1.2 inch, dan 1.5 inch. Adapun batas-batas daerah dapat dilihat pada gambar di bawah ini: 3.8.1.1 Bentuk coalescer dengan sekat 1.2 inch outlet oil
inlet
sekat outlet water
1.2 inch sludge Gambar 3.8 Gambar Coalescer sekat 1.2 inch
47
3.8.1.2 Bentuk coalescer dengan sekat 1 inch outlet oil
inlet
sekat outlet water
1 inch sludge Gambar 3.9 Gambar coalescer sekat 1 inch 3.8.1.3 Bentuk coalescer dengan sekat 1.5 inch inlet
sekat
sludge
outlet oil
outlet water
1.5 inch Gambar 3.10 Gambar coalescer sekat 1.5 inch 48
Tabel 3.1 Karakteristik Characteristic
Units
Value
OS – 80
GPM
240
Width
m
1.6256
Hight
m
2.0066
Length
m
2.3622
Inlet
m
0.1524
Nozzle Nominal
Water Outlet
m
0.1524
Diameter
Oil Outlet
m
0.1016
Sludge Outlet
m
0.1016
Operating Pressure
psig
88
Oil
kg/m3
824.3
Water
kg/m3
996.2
Oil
kg/s
0.0032155
Water
kg/s
6.4276845
Model Vessel Dimension
Operating Condition Density Mass flow
3.8.2 Pembentukan Mesh Ukuran mesh yang diterapkan pada model akan mempengaruhi ketelitian analisis CFD. Semakin kecil ukuran mesh pada model, maka hasil yang didapatkan akan semakin teliti, tetapi membutuhkan daya komputasi dan waktu yang lebih lama dibandingkan dengan mesh yang memiliki ukuran lebih besar. Oleh karena itu, besar ukuran mesh harus diatur sedemikian rupa sehingga diperoleh hasil yang teliti dan diusahakan daya komputasi yang dibutuhkan tidak terlalu besar. Jenis volume meshing seperti pada Tabel 3.2 Tabel 3.2 Jenis Volume Meshing Elements
Quad
Type
Pave
Internal size
0.02 49
Gambar 3.11 Gambar Mesh Coalescer
3.8.3 Solver dan Model Turbulensi Simulasi pada tugas akhir ini menggunakan solver unsteady dan menggunakan model k-epsilon, model k-epsilon standard dan k epsilon RNG, dan near wall treatment standard wall function. Gambar 3.12 dan 3.13 menunjukkan pengisian model yang digunakan.
50
Gambar 3.12 Model k epsilon standard
Gambar 3.13 Model k epsilon RNG 51
3.8.4 Pendifinisian Kondisi Batas a. Mass flow inlet
Gambar 3.14 Model Mass flow inlet Dimana turbulen menggunakan rumus:
(3.1)
a. Masukan nilai 3 % pada turbulent intensity Turbuleny intensity dapat dihitung menggunakan rumus:
b. Masukan 0.1524 pada hydraulic diameter
52
b. Mass flow inlet water
Gambar 3.15 Model Mass flow inlet water
c. Mass flow inlet oil
Gambar 3.16 Model Mass flow inlet oil
53
d. Out flow oil
Gambar 3.17 Model outflow oil
Masukkan 0.0005 pada flow rate weightinh, yang didapat dari:
e. Out flow water
Gambar 3.18 Model outflow water Masukkan 0.9995 pada flow rate weightinh, yang didapat dari:
.9995
54
BAB IV HASIL SIMULASI 4.1
Hasil Analisa Simulasi Simulasi coalescer menggunakan model cell antara 1.000 hingga 10.000
dan 10.000 hingga 20.000. Simulasi menggunakan beberapa model k-epsilon (k-ε) dan k-omega (k-ω) yaitu k-epsilon (k-ε) standard, k-epsilon (k-ε) RNG, k-omega (k-ω) standard, dan k-omega (k-ω) SST yang akan di bandingkan dengan hasil oil water separator ( coalescer ) yaitu kandungan minyaknya 10 ppm. 4.1.1
Hasil Simulasi dengan Beberapa Model k-epsilon (k-ε) dan k-omega (k-ω) Simulasi ini dilakukan dengan model k-epsilon (k-ε) dan k-omega (k-ω)
yaitu k-epsilon (k-ε) standard, k-epsilon (k-ε) RNG, k-omega (k-ω) standard dan k-omega (k-ω) SST dengan variasi jarak sekat 1 inch, 1.2 inch, dan 1.5 inch dengan jumlah cell 1.000 hingga 10.000. Pada Gambar 4.1 hingga 4.4 menunjukkan sumbu absis yang menerangkan tentang waktu simulasi sedangkan sumbu ordinat yang menerangkan kandungan minyak. Pada model k-epsilon (k-ε) standard dilakukan waktu simulasi 1320 detik sedangkan untuk model k-epsilon (k-ε) RNG, k-omega (k-ω) standard dan k-omega (k-ω) SST dilakukan waktu simulasi 660 detik. Tabel 4.1 Hasil Simulasi Model k-epsilon (k-ε) Standard Detik 60 120 180 240 300 360 420 480 540 600 660
1 inch (PPM) 0 0 28.9 47.4 55 60.5 64.9 68.8 71.8 74.5 76.7
1.2inch (PPM) 0 0 32.3 49.1 55.8 60.7 64.7 68 71.1 73.7 75.7
1.5inch (PPM) 0 0 34.6 49 55 59 63.1 66.8 69.7 72.1 75.6
Detik 720 780 840 900 960 1020 1080 1140 1200 1260 1320
1 inch (PPM) 78.3 79.5 80.8 81.9 82.9 83.6 84.3 84.9 85.4 85.9 86.6
1.2inch (PPM) 77.3 78.8 80.1 81 81.9 82.6 83.4 84.4 85.2 85.9 86.5
1.5inch (PPM) 75.8 77.5 78.6 79.7 80.8 82 83.5 84.7 85.8 86.6 87.1
55
k-epsilon standard
PPM
100 90 80 70 60 Kandungan 50 minyak 40 30 20 10 0
1 inch 1.2 inch 1.5 inch
0
500
1000
1500 Detik
Waktu Simulasi
Gambar 4.1 Perbandingan grafik hasil berapa PPM yang terjadi pada jarak sekat 1, 1.2, dan 1.5 inch dengan model k-epsilon (k-ε) standard.
Tabel 4.2 Hasil Simulasi Model k-epsilon (k-ε) RNG
Detik 60 120 180 240 300 360 420 480 540 600 660
1 inch (PPM) 0 0 36.3 48.2 53.9 58.8 63 67.1 70.7 74 77.9
1.2inch (PPM) 0 0 39 50.3 55 58.3 62.5 66.2 69.5 72.4 74.6
1.5inch (PPM) 0 0 39.7 49 54.9 56.8 60.9 64.8 67.9 70.5 72.6
56
k-epsilon rng
PPM
90
1 inch
80
1.2 inch
70
1.5 inch
60 Kandungan minyak
50 40 30 20 10 0 0
200
400
Waktu Simulasi
600
800
Detik
Gambar 4.2 Perbandingan grafik hasil berapa PPM yang terjadi pada jarak sekat 1, 1.2, dan 1.5 inch dengan model k-epsilon (k-ε) RNG.
Tabel 4.3 Hasil Simulasi Model k-omega (k-ω) Standard
Detik 60 120 180 240 300 360 420 480 540 600 660
1 inch (PPM) 0 0 27.9 45.5 51.6 57 62.9 67.9 71.7 74.8 77
1.2inch (PPM) 0 0 33 48.4 53.8 57.9 62.7 67 70.6 73.4 75.5
1.5inch (PPM) 0 0 36 48.8 53.9 56.7 61.5 64.9 68.3 70.8 72.5
57
k omega standard
PPM
90 80
1 inch
70
1.2 inch 1.5 inch
60 50 Kandungan minyak
40 30 20 10 0 0
200
400
600
800
Detik
Waktu Simulasi
Gambar 4.3 Perbandingan grafik hasil berapa PPM yang terjadi pada jarak sekat 1, 1.2, dan 1.5 inch dengan model k-omega (k-ω) standard.
\ Tabel 4.4 Hasil Simulasi Model k-omega (k-ω) SST
Detik 60 120 180 240 300 360 420 480 540 600 660
1 inch (PPM) 0 0 30.2 46.8 53.5 58.1 62.6 67.3 70.8 74 76.5
1.2inch (PPM) 0 0 35.2 50 55.8 59.2 62.5 66.2 69.5 72.3 74.5
1.5inch (PPM) 0 0 37.9 50.2 56.5 58.8 61.2 64.2 66.7 69.5 71.4
58
k-omega sst
PPM
90
1 inch 1.2 inch 1.5 inch
80 70 60 Kandungan minyak
50 40 30 20 10 0 0
200
400
600
Waktu Simulasi
800
Detik
Gambar 4.4 Perbandingan grafik hasil berapa PPM yang terjadi pada jarak sekat 1, 1.2, dan 1.5 inch dengan model k-omega (k-ω) SST. Tabel 4.1, 4.2, 4.3, 4.4 dan Gambar 4.1, 4.2, 4.3, 4.4 menujukkan data hasil simulasi dan grafik perbandingan jarak sekat 1 inch, 1.2 inch dan 1.5 inch dengan menggunakan model k-epsilon (k-ε) standard,k-epsilon (k-ε) RNG, komega (k-ω) standard dan k omega (k-ω) SST. Gambar 4.1, 4.2, 4.3 dan 4.4 menunjukkan pada detik ke 60 sampai detik ke 120 belum terlihat kandungan minyak yang didapatkan karena dibutuhkan waktu simulasi untuk mencapai tempat keluarnya air. Pada Gambar 4.1 menunjukkan pada detik ke 180 sampai 1320 sudah mulai terlihat kandugan minyak yang didapatkan tiap jarak sekatnya dan terus mengalami peningkatan kandungan minyak. Pada Gambar 4.2, 4.3, 4.4 pada detik ke 180 sampai 660 sudah mulai terlihat kandungan minyak yang didapatkan tiap jarak sekatnya dan mengalami peningkatan kandungan minyak tiap waktu simulasinya. Dilihat pada detik 180 didapatkan kandungan minyak yang terendah yaitu pada sekat 1 inch.
59
Tabel 4.5 Hasil Simulasi Model k-epsilon (k-ε) Standard
Detik 60 120 180 240 300 360 420 480 540 600 660
1 inch (PPM) 0 0 34.3 49.9 56.2 61.9 67.7 71.6 74.3 76.5 78
1.2inch (PPM) 0 0 38.2 51.3 56.5 61.7 66.8 70.5 73.3 75.5 77
1.5inch (PPM) 0 0 39.5 50.9 56.1 61.2 66.3 69.5 71.9 73.6 75.3
k-epsilon standard
PPM
90
1 inch
80
1.2 inch
70
1.5 inch
60 Kandungan 50 minyak
40 30 20 10 0 0
200
400
600
800
Detik
Waktu Simulasi
Gambar 4.5 Perbandingan grafik hasil berapa PPM yang terjadi pada jarak sekat 1, 1.2, dan 1.5 inch dengan model k-epsilon (k-ε) standard.
60
Tabel 4.6 Hasil Simulasi Model k-epsilon (k-ε) RNG
Detik 60 120 180 240 300 360 420 480 540 600 660
1 inch (PPM) 0 0 44.8 54.4 61.2 64.9 67.3 69.1 70.9 72.8 74.5
1.2inch (PPM) 0 0 45.8 53.3 58.7 61.7 64.2 66.2 68.8 72.2 74.7
1.5inch (PPM) 0 0 44 51.9 56.8 59.9 62.5 66.4 69.9 72.7 74
k-epsilon rng
PPM
80 1 inch
70 60
1.2 inch 1.5 inch
50
Kandungan minyak 40
30 20 10 0 0
200
400
600
800
Detik
Waktu Simulasi
Gambar 4.6 Perbandingan grafik hasil berapa PPM yang terjadi pada jarak sekat 1, 1.2, dan 1.5 inch dengan model k-epsilon (k-ε) RNG.
61
Tabel 4.7 Hasil Simulasi Model k-omega (k-ω) Standard Detik 0 60 120 180 240 300 360 420 480 540 600 660
1 inch (PPM) 0 0 0 34.9 47.9 52.7 59.5 65.8 69.8 72.8 75.2 76.9
1.2inch (PPM) 0 0 0 39.8 49.7 53.4 59.1 64.5 68.1 71.3 73.5 75.2
1.5inch (PPM) 0 0 0 41.9 49.9 53.4 58.6 63.2 66.5 69.1 70.8 72.6
k-omega standard
PPM
90 80
1 inch
70
1.2 inch
60 Kandungan minyak
1.5 inch
50 40 30 20 10 0 0
200
400
600
800
Waktu Simulasi
Detik
Gambar 4.7 Perbandingan grafik hasil berapa PPM yang terjadi pada jarak sekat 1, 1.2, dan 1.5 inch dengan model k-omega (k-ω) standard.
62
Tabel 4.5, 4.6, 4.7 dan Gambar 4.5, 4.6, 4.7 menujukkan hasil data simulasi dan grafik perbandingan jarak sekat 1 inch, 1.2 inch dan 1.5 inch dengan menggunakan model k-epsilon (k-ε) standard, k-epsilon (k-ε) RNG, k-omega (kω) standard dengan jumlah cell 10.000 hingga 20.000. Pada Gambar 4.5 hingga 4.7 menunjukkan
sumbu absis yang menerangkan tentang waktu simulasi
sedangkan sumbu ordinat yang menerangkan kandungan minyak. Gambar 4.5, 4.6, 4.7 menunjukkan pada detik ke 60 sampai 120 belum terlihat kandungan minyak yang didapatkan karena dibutuhkan waktu simulasi untuk mencapai tempat keluarnya air. Pada detik ke 180 sampai 660 mulai terlihat kandungan minyak yang didapatkan dan mengalami peningkatan kandungan minyak tiap waktu simulasinya dan pada detik 180 mendapatkan kandungan minyak yang terendah yaitu jarak sekat 1 inch. Dapat disimpulkan dari berbagai cara yang dilakukan pada simulasi coalescer dan hasil yang didapatkan dengan menggunakan beberapa model kepsilon (k-ε) dan k-omega (k-ω) dengan cell 1.000 hingga 10.000 dan 10.000 hingga 20.000 dapat dilihat bahwa yang mendekati hasil referensi 10 ppm yaitu dengan menggunakan sekat 1 inch pada detik ke 180 yaitu 27 ppm sampai 40 ppm karena sekat tersebut memiliki jarak yang kecil dan juga pengaruh viskositas (kekentalan) karena minyak memiliki viskositas yang lebih besar dari pada air, pada saat minyak melewati sekat tersebut minyak akan menempel di sekat tersebut dan menghambat aliran minyak tersebut sehingga mendapatkan kandungan minyak yang lebih kecil dari pada sekat 1,2 inch dan 1,5 inch.
63
Simulasi ini menggunakan jarak sekat yaitu 1 inch, 1.2 inch dan 1.5 inch dengan menggunakan beberapa model k-epsilon (k-ε) dan k-omega (k-ω). Pada Gambar 4.8 hingga 4.10 menunjukkan sumbu absis yang menerangkan tentang waktu simulasi sedangkan sumbu ordinat yang menerangkan kandungan minyak. Tabel 4.8 Hasil Model Simulasi 1 Inch
Detik 60 120 180 240 300 360 420 480 540 600 660
k-epsilon (k-ε) standard (PPM) 0 0 28.9 47.4 55 60.5 64.9 68.8 71.8 74.5 76.7
k-epsilon (k-ε) RNG (PPM) 0 0 36.3 48.2 53.9 58.8 63 67.1 70.7 74 77.9
k-omega (k-ω) standard (PPM) 0 0 27.9 45.5 51.6 57 62.9 67.9 71.7 74.8 77
k-omega (k-ω) SST (PPM) 0 0 30.2 46.8 53.5 58.1 62.6 67.3 70.8 74 76.5
1 inch
PPM
90 k epsilon standard
80 70
k epsilon rng
60 Kandungan minyak
50
k omega standard
40
k omega sst
30 20 10 0 0
200
400
600
800 Detik
Waktu Simulasi
Gambar 4.8 Gambar grafik perbandingan permodelan pada jarak sekat 1 inch.
64
Tabel 4.9 Hasil Model Simulasi 1.2 Inch
Detik 60 120 180 240 300 360 420 480 540 600 660
k-epsilon(k-ε) standard k-epsilon(k-ε) (PPM) RNG (PPM) 0 0 0 0 32.3 39 49.1 50.3 55.8 55 60.7 58.3 64.7 62.5 68 66.2 71.1 69.5 73.7 72.4 75.7 74.6
k-omega (k-ω) standard (PPM) 0 0 33 48.4 53.8 57.9 62.7 67 70.6 73.4 75.5
k-omega (k-ω) SST (PPM) 0 0 35.2 50 55.8 59.2 62.5 66.2 69.5 72.3 74.5
1.2 inch
PPM
80
k epsilon standard k epsilon rng k omega standard k omega sst
70 60 50 Kandungan minyak
40 30 20 10 0 0
200
400
600
800 Detik
Waktu Simulasi
Gambar 4.9 Gambar grafik perbandingan permodelan pada jarak sekat 1.2 inch.
65
Tabel 4.10 Hasil Model Simulasi 1.5 Inch
Detik 60 120 180 240 300 360 420 480 540 600 660
k-epsilon(k-ε) standard k-epsilon(k-ε) (PPM) RNG (PPM) 0 0 0 0 34.6 39.7 49 49 55 54.9 59 56.8 63.1 60.9 66.8 64.8 69.7 67.9 72.1 70.5 75.6 72.6
PPM
k-omega (k-ω) standard (PPM) 0 0 36 48.8 53.9 56.7 61.5 64.9 68.3 70.8 72.5
k-omega (k-ω) SST (PPM) 0 0 37.9 50.2 56.5 58.8 61.2 64.2 66.7 69.5 71.4
1.5 inch 80
k epsilon standard k epsilon rng k omega standard k omega sst
70 60 50 Kandungan minyak
40 30 20 10 0 0
200
400
600
800
Detik
Waktu Simulasi
Gambar 4.10 Gambar grafik perbandingan permodelan pada jarak sekat 1.5 inch.
66
Tabel 4.8 hingga 4.10 dan Gambar 4.8 hingga 4.10 menunjukkan data hasil simulasi dan grafik perbandingan permodelan tiap jaraknya, maka didapatkan hasilnya yaitu jarak sekat 1 inch yang mendekati hasil referensi kandungan minyak 10 ppm yaitu pada model k-epsilon dengan kandungan minyak 28.9 dan k-omega standard dengan kandungan minyak 27.9 ppm, dengan jarak sekat 1.2 inch pada model k-epsilon standard dengan kandungan minyak 32.3 ppm dan untuk 1.5 inch pada model k-epsilon standard dengan kandungan minyak 34.6 inch. Model k-epsilon lebih mendekati hasil data referensi dibandingkan dengan model yang lain. Hal itu disebabkan oleh kecenderungan model ini memiliki tingkat akurasi yang cukup memadai untuk simulasi fluida.
67
4.2 4.2.1
Penggambaran Display dan velocity Vectors didalam Coalescer Display Contours Velocity
Gambar 4.11 Gambar contours velocity coalescer sekat 1 inch pada detik ke 180.
Gambar 4.12 Gambar contours velocity coalescer sekat 1.2 inch pada detik ke 180. 68
Gambar 4.13 Gambar contours velocity coalescer sekat 1.5 inch pada detik ke 180.
Gambar 4.11, 4.12, 4.13 menunjukkan counturs kecepatan pada tiap jarak sekat 1 inch, 1.2 inch dan 1.5 inch. Pada Gambar 4.11, 4.12 dan 4.13 menunjukkan counturs kecepatan sama akan tetapi pada saat memasuki sekat (media flopak) terjadi perbedaan karena sekat 1 inch memiliki jarak sekat yang lebih kecil sehingga sewaktu memasuki sekat akan tertahan dan kecepatannya lebih kecil dibandingkan dengan jarak sekat 1.2 inch dan 1.5 inch.
69
4.2.2
Display Contours Pressure
Gambar 4.14 Gambar contours pressure coalescer sekat 1 inch pada detik ke 180.
Gambar 4.15 Gambar contours pressure coalescer sekat 1.2 inch pada detik ke 180. 70
Gambar 4.16 Gambar contours pressure coalescer sekat 1.5 inch pada detik ke 180.
Gambar 4.14, 4.15, 4.16 menunjukkan counturs pressure pada tiap jarak sekat 1 inch, 1.2 inch dan 1.5 inch pada detik ke 180. Pada Gambar 4.14, 4.15 dan 4.16 menunjukkan counturs pressure sama pada saat memasuki coalescer tetapi pada saat memasuki media flopak terjadi perbedaan pressure pada sekat 1 inch pressure lebih kecil dari pada sekat 1.2 inch dan 1.5 karena jarak sekatnya lebih kecil dan pada sisi outlet air terjadi perbedaan counturs pressure.
71
4.2.3
Velovity Vectors Inlet
Gambar 4.17 Velocity vector di dalam coalescer pada detik ke 180 pada sekat 1 inch.
Gambar 4.17 menunjukan vector kecepatan fluida dalam tanki coalescer dimana terlihat kecepatan maksimum dapat mencapai 0.0747 m/s, kecepatan tinggi terjadi pada sisi inlet, sisi outlet air. Kecepatan mengalami penurunan pada sisi masuk karena terdapat diffuser yang mengarah kebawah, kemudian fluida akan melewati flopak media agar kotoran-kotoran yang terkandung dalam fluida dapat terpisah dan kotoran tersebut akan tersaring dan tertampung di tangki kotoran kemudian dikeluarkan melalui pipa kotoran, selain memisahkan menyaring kotoran flpak media berfungsi untuk melaminarkan aliran fluida guna untuk mempercepat pemisahan fluida.
72
4.2.4
Velovity Vectors Volume Fraction
Gambar 4.18 Velocity vectors by volume fraction campuran minyak dan air di dalam coalescer pada detik ke 660 pada sekat 1 inch.
Gambar 4.18 menunjukan vector kecepatan fraksi volume dengan model 2D. Pada Gambar 4.18 dilihat fenomena aliran fluida yang terjadi di dalalm tangki coalescer dimana pemisahan antara air dan minyak terjadi.
73
BAB V KESIMPULAN DAN SARAN
5.1
Kesimpulan
Dari simulasi yang telah dilakukan dan dibahas maka dapat disimpulkan bahwa: 1. Variasi jarak sekat 1 inch, 1.2 inch dan 1.5 inch di dalam coalescer didapatkan hasil yaitu jarak sekat 1 inch mendapatkan hasil kandungan minyak yang lebih kecil dibandingkan jarak sekat 1.2 inch dan 1.5 inch. Jarak sekat 1 inch mendapatkan hasil kandungan minyak sebesar 27.9 PPM dan mendekati hasil referensi sebesar 10 PPM. 2. Jarak sekat ( flopak media ) sangat berpengaruh pada hasil kandungan minyak yang didapatkan di dalam coalescer.
5.2
Saran 1. Simulasi ini memerlukan ketelitian pada saat membuat geometri dan meshing sehingga menghasilan hasil yang akurat. 2. Jarak sekat ( flopak media ) harus diperkecil sehingga mendapatkan hasil yang sesuai dengan referensi.
74
DAFTAR PUSTAKA 1. Amoanu, Daniel, Studies of Liquid-liquid Separation in Fibre Coalescer, LAPPEENRANTA UNIVERSITY OF TECHNOLOGY, October 2008. 2. http://www.scribd.com/doc/24278391/Catatan-Kuliah-Teknik-Produksi 3. http://www.wash-water-treatment.com/OS2007.pdf 4. Grubich, Jeffrey, dan De Guzman Angelo, Designing Coalescing Oil/Water Separator for Use at Army Washracks, Engineer Research and Development Center, December 2000. 5. http://en.wikipedia.org/wiki/API_oil-water_separator 6. http://wtcommodities.com/coalescing.php 7. Gerhart, P, M., and Richard J. G., 1985, Fundamentals of Fluid Mechanics, Addison-Wesley Publishing Company, Inc., USA, 509530.Tuakia,
Firman,
Dasar-dasar
CFD
Menggunakan
FLUENT,
Informatika, Bandung, 2008. 8. Fox, R. W., and Alan, T. M. D., 1994, Introduction to Fluid Mechanics, Fourth Edition, SI Version, John Wiley & Sons, Inc., Canada, 2-33. 9. FLUENT, Inc, Manual FLUENT Documentation, FLUENT, Inc, 1998. 10. Munson, B. R., Young, D. F., dan Okiishi, T. H., 2008, Fluid Dynamics, Jilid 2, Erlangga, Jakarta, 111-473. 11. Versteeg, H. K., and Malalasakera, W., 1995, an Introduction to Computational Fluid Dynamic, Longman Scientific & Technical, England, 11-25. 12. Mohr. Kirby S. Introduction to Separation of Oil adn Water, Oklahoma, USA, 2001-2008. 13. Mohr. Kirby S. How Oil Water Seaparator Work and How To Be Use Them, 2001-2010. 14. Tuakia, Firman, Dasar-dasar CFD Menggunakan FLUENT, Informatika, Bandung, 2008.
LAMPIRAN
Gambar A1. Gambar contours velocity coalescer sekat 1 inch pada detik ke 60.
Gambar A2. Gambar contours velocity coalescer sekat 1 inch pada detik ke 120.
Gambar A3. Gambar contours velocity coalescer sekat 1 inch pada detik ke 180.
Gambar A4. Gambar contours velocity coalescer sekat 1.2 inch pada detik ke 60.
Gambar A5. Gambar contours velocity coalescer sekat 1.2 inch pada detik ke 120.
Gambar A6. Gambar contours velocity coalescer sekat 1.2 inch pada detik ke 180.
Gambar A7. Gambar contours velocity coalescer sekat 1.5 inch pada detik ke 60.
Gambar A8. Gambar contours velocity coalescer sekat 1.5 inch pada detik ke 180.
Gambar A9. Gambar contours pressure coalescer sekat 1 inch pada detik ke 60.
Gambar A10. Gambar contours pressure coalescer sekat 1 inch pada detik ke 120.
Gambar A11. Gambar contours pressure coalescer sekat 1 inch pada detik ke 180.
Gambar A12. Gambar contours pressure coalescer sekat 1.2 inch pada detik ke 60.
Gambar A13. Gambar contours pressure coalescer sekat 1.2 inch pada detik ke 120.
Gambar A14. Gambar contours pressure coalescer sekat 1.2 inch pada detik ke 180.
Gambar A15. Gambar contours pressure coalescer sekat 1.5 inch pada detik ke 60.
Gambar A16. Gambar contours pressure coalescer sekat 1.5 inch pada detik ke 180.
Gambar A17. Velocity vector di dalam coalescer pada detik ke 180 pada sekat 1 inch.
Gambar A18. Velocity vector di dalam coalescer pada detik ke 180 pada sekat 1.2 inch.
Gambar A19. Velocity vector di dalam coalescer pada detik ke 180 pada sekat 1.5 inch.
Gambar A20. Velocity vectors by volume fraction campuran minyak dan air di dalam coalescer pada detik ke 660 pada sekat 1 inch.
Gambar A21. Velocity vectors by volume fraction campuran minyak dan air di dalam coalescer pada detik ke 660 pada sekat 1.2 inch.
Gambar A22. Velocity vectors by volume fraction campuran minyak dan air di dalam coalescer pada detik ke 660 pada sekat 1.5 inch.
