ANALISIS PROSES PERPINDAHAN PANAS DUA DIMENSI PADA DUA DAERAH ALIRAN FLUIDA YANG TERPISAH OLEH SEBUAH PELAT DENGAN PERANGKAT LUNAK FLUENT V5

44

ISSN 0216 - 3128

V.I.S. Wardhani, dkk.

ANALISIS PROSES PERPINDAHAN PANAS DUA DIMENSI PADA DUA DAERAH ALIRAN FLUIDA YANG TERPISAH OLEH SEBUAH PELAT DENGAN PERANGKAT LUNAK FLUENT V5 V.I.S. Wardhani, Henky P.R, Efrizon Umar Pusat Penelitian Dan Pengembangan Teknik Nuklir

ABSTRAK ANALISIS PROSES PERPINDAHAN PANAS DUA DIMENSI PADA DAERAH ALIRAN FLUIDA YANG TERPISAH OLEH SEBUAH PELAT DENGAN PERANGKAT LUNAK FLUENT V5. Untuk menggambarkan dan menganalisis proses perpindahan panas yang sesungguhnya terjadi pada suatu mesin penukar kalor merupakan persoalan yang cukup rumit. Karena itu perlu dilakukan tahap-tahap penyederhanaan dalam mensimulasikannya. Tahap yang paling sederhana adalah dengan analisis dua dimensi pada dua daerah aliran fluida yang terpisah oleh sebuah pelat seperti yang dilakukan di sini. Analisis dilakukan dengan perangkat lunak FLUENT V5. Hasil analisis yang berupa harga-harga karakteristik fluida seperti distribusi kecepatan maupun distribusi temperatur dapat ditampilkan dalam bentuk grafik. Kata kunci : Fluent V5, mesin penukar kalor, karakteristik fluida.

ABSTRACT THE ANALYSING OF THE HEAT TRANSFER PROCESS TWO DIMENSIONAL ON THE TWO REGIONS FLUID FLOW SEPARATED BY A PLATE, BY USING FLUENT V5 SOFTWARE. Describing and analysing the heat transfer process on the heat exchanger is become rather complicated problem. On simulating this problem, it can be solved in step by step. A simplest case, as is written here, this two regions fluid flow separated by a plate problem is analysed by two dimensional. This is analysed by using Fluent V5 software. The result of this analysing is fluid characteristics, such as velocity and temperature distribution that can be described as a graphics. Key words : Fluent V5, heat exchanger, fluid characteristics.

ke dua daerah aliran fluida yang digambarkan dalam bentuk grafik.

PENDAHULUAN

P

erpindahan panas yang terjadi pada sistem penukar kalor melibatkan proses konduksi dan konveksi. Untuk mengetahui dan memahami proses perpindahan panas yang terjadi, perlu dilakukan analisis terhadap sistem tersebut. Analisis yang demikian cukup rumit oleh karena itu perlu dilakukan penyederhanaan, seperti yang dilakukan pada penulisan di sini, yaitu mengassumsikan terlebih dahulu sebagai masalah dengan dua daerah aliran fluida yang terpisah oleh sebuah pelat dalam bentuk dua dimensi. Langkah ini merupakan langkah awal dari target yang lebih jauh lagi yaitu analisis tiga dimensional pada sistem perpindahan panas yang terjadi pada mesin penukar kalor. Analisis dilakukan dengan menggunakan program Fluent V5. Hasil analisis berupa distribusi kecepatan, maupun distribusi temperatur pada fluida yang mengalir maupun pada pelat pemisah

Fluent V5 adalah suatu program komputer yang ditulis dengan bahasa C yang dapat memodelkan aliran fluida dan proses perpindahan panas yang terjadi pada geometri dua dimensi maupun tiga dimensi. Persamaan atur yang diselesaikan adalah persamaan Navier Stokes yang meliputi tiga hukum kekekalan, yaitu kekekalan massa, kekekalan momentum dan kekekalan energi. Masalah-masalah aliran fluida yang mampu diselesaikan dengan Fluent V5 meliputi geometri 2D/3D, aliran inkom-presibel/kompresibel, tunak/tak tunak, aliran laminar /turbulen, aliran fluida Newtonian/non Newtonian, perpindahan panas konveksi baik konveksi alamiah maupun konveksi paksa, gabungan konduksi/ konveksi, radiasi dan jenis aliran pada lainnya. Dengan melihat ke anekaragaman aliran fluida yang mampu di-selesaikan dengan Fluent V5, maka sangat beralasan bahwa program Fluent V5 ini

Prosiding PPI - PDIPTN 2005 Puslitbang Teknologi Maju - BATAN Yogyakarta, 12 Juli 2005

45

ISSN 0216 - 3128

dapat diaplikasikan pada analisis terhadap prosesproses yang terjadi pada banyak instrument, antara lain mesin penukar kalor, kondensor, evaporator, sistem pemipaan, sistem ventilasi dan sebagainya. [2]

TEORI Karakteristik fluida yang mencerminkan keadaan aliran fluida sebenarnya di alam adalah kompresibel (massa jenis tidak konstan), rotasional (partikel fluida berputar), tak tunak (tidak tetap) dan viskos (terdapat geseran dengan permukaan). Per-samaan atur dinamika fluida yang dapat menggam-barkan semua kondisi di atas dikenal dengan nama persamaan Navier – Stokes[4]. Perangkat lunak FLUENT V5 mampu menyelesaikan persamaan Navier – Stokes dengan dasar persamaan kekekalan massa, kekekalan momentum serta kekekalan energi untuk aliran fluida yang melibatkan proses perpindahan panas. Persamaan kekekalan massa atau lebih dikenal dengan persamaan kontinyuitas dapat dituliskan sebagai berikut[3]:


Dengan  adalah viskositas molekuler, sementara suku ke dua pada ruas kanan adalah pengaruh dari dilatasi volume. Pada aliran fluida yang melibatkan proses perpindahan panas, seperti kasus yang akan dianalisis di sini harus ditambahkan pula persamaan kekekalan energi yang dapat dituliskan dalam bentuk enthalpi h sebagai berikut [3] :    T (  h)  (  ui h)  (k  kt ) t xi xi xi 

 xi

j1

1

Dp ui  ( ik )eff  Sh Dt xk

j j1 

(4)

dengan : k : konduktivitas molekuler kt : konduktivitas turbulen Jj1 : fluks diffusi dari spesies ji Sh : sumber-sumber panas yang didefinisikan Enthalpi h didefinisikan sebagai :

   (  ui )  S m t xi

(1)

h   m j1 h j1

(5)

j1

dengan:  : massa jenis t : waktu x : arah aksial u : kecepatan Sm : massa yang ditambahkan

dengan mj1 adalah fraksi massa dari proses j1 dan h j1 

T

 c p, j

1

dT

(6)

Tref

dengan Tref = 298,15 oK

Persamaan kekekalan momentum dalam arah i dalam kerangka acuan inertia dapat dituliskan sebagai berikut(3) :   p  ij (  ui )  (  ui u j )      g i  Fi t x j xi x j

(2)

FLUENT V5 menyelesaikan kekalan energi (4) dengan temperatur sebagai variable bebas. persamaan (4) ditransformasikan temperatur dengan menggunakan dan (6).

persamaan kemenggunakan Ruas kiri dari dalam bentuk persamaan (5)

Pada daerah padatan (solid), FLUENT V5 menyelesaikan persamaan konduksi sederhana yang melibatkan fluks panas karena konduksi dari sumber-sumber panas volumetrik dalam padatan (solid)[3] :

dengan : p : tekanan statis

ij : tensor tegangan gi : gaya gravitasi Fi : gaya-gaya luar dalam arah i

   T k h  t xi  xi

Tensor tegangan ij diberikan dengan :   u u  ij     i  j  xi   x j

 hj

  2 ui      3 xi ij 

dengan : (3)

 h k

: massa jenis : entahalpi, : konduktivitas


   q 

(7)

T

46

ISSN 0216 - 3128


: temperatur

6.

Memilih jenis persamaan dasar yang harus diselesaikan yaitu laminar atau turbulen (atau tidak viskos), reaksi-reaksi kimia (optional), model-model perpindahan panas dan lain-lain.

7. 8.

Menentukan jenis material Menentukan harga-harga kondisi batas

9.

Mengatur parameter- parameter kontrol solusi

q : sumber panas volumetrik Persamaan (7) diselesaikan secara simultan dengan persamaan kekekalan energi (4) untuk memprediksi proses perpindahan panas gabungan konduksi/ konveksi. Di samping solusi terhadap ke tiga persamaan kekekalan di atas, FLUENT V5 juga menyediakan fasilitas pemodelan turbulensi. Turbulensi di-modelkan dengan model turbulen ”dua persamaan”. Dalam model “dua persamaan” ini, efek turbulen direpresentasikan dengan viskositas isotropik “eddy” atau “turbulen” yang dievaluasi dengan mengguna-kan dua buah besaran yaitu: energi kinetik turbulen (k) dan laju dissipasi (). Ke dua besaran turbulen ini, k dan  diperoleh dari solusi persamaan model transport. FLUENT menyediakan dua pilihan model k -  sebagai berikut:  model k -  standard  model k -  berdasarkan Renormalisasi – Group (RNG) Pemilihan bentuk model turbulen yang digunakan tergantung dari pemahaman terhadap kondisi aliran yang dianalisis.

Langkah-Langkah Perhitungan[2]

10. Menentukan kondisi awal aliran fluida 11. Melakukan proses perhitungan 12. Mempelajari perhitungan

dan

mengevaluasi

hasil-hasil

13. Menyimpan hasil-hasil perhitungan Langkah 1 dan 2 memerlukan perangkat lunak untuk membuat model geometri dan jaringan kisi yaitu perangkat lunak gambit version 1.1. Langkah-langkah selanjutnya dilakukan dengan menggunakan menu-menu yang terdapat dalam perangkat lunak FLUENT V5.

Batasan Masalah Kasus yang akan dianalisis di sini adalah kondisi dua aliran fluida dengan arah yang berlawanan, yang dipisahkan oleh sebuah pelat yang akan dimodelkan secara dua dimensi seperti digambarkan pada Gambar.

Dengan memahami dan mengetahui masalah yang akan diselesaikan, langkah berikutnya adalah memodelkan dan menganalisis bentuk model tersebut dengan menggunakan perangkat lunak Gambit versi 1.1 dan Fluent V 5. Prosedur langkah yang harus dilakukan dapat dituliskan sebagai berikut : 1.

Menggambarkan bentuk geometri model

2.

Membuat jaringan kisi pada geometri model

3.

Menentukan jenis perangkat lunak yang akan digunakan, serta menentukan dimensi model yang akan diselesaikan

4.

Menyimpan jaringan kisi yang terbentuk

Setting Kondisi Batas[2]

5.

Mengevaluasi jaringan kisi yang terbentuk

Kondisi batas pemodelan akan dilakukan dengan menggunakan menu-menu yang terdapat pada FLUENT V5, dapat dilihat pada Gambar 2.

Gambar 1. Kondisi aliran yang akan dimodelkan.


47

ISSN 0216 - 3128


Gambar 2. Kondisi batas pada model. Setting pada inflow 1 adalah VELOCITY yang sangat signifikan. Jumlah titik kisi dalam arah INLET dengan kecepatan 4,5 m/s dan temperatur aliran fluida (arah x) sebanyak 250, sementara 300 K , sementara setting pada inflow 2 adalah dalam arah melintang terhadap aliran ( arah y) VELOCITY INLET dengan kecepatan 4,5 m/s dan sebanyak 240 dengan ratio kerapatan masingtemperatur 318 K. Setting kondisi batas hargamasing sebesar 1.3. Bentuk jaringan kisi dapat harga kecepatan dan temperatur tersebut dilihat pada Gambar 3. merupakan harga pendekatan pada kondisi aliran fluida pendingin primer dan sekunder yang HASIL DAN PEMBAHASAN terjadi pada penukar kalor di reactor Triga 2000. Setting pada dinding dibiarkan bebas dengan fluks Kriteria konvergensi dilakukan terhadap 6 panas sama dengan 0 (nol). Jenis fluida adalah air (enam) buah nilai parameter yaitu kontinyuitas, (H20) sementara pelat pemisah terbuat dari baja. kecepatan arah x, kecepatan arah y, turbulensi (k) ,

Pembuatan Jaringan Kisi[1] Titik-titik pembentuk jaringan kisi pada sisisisi inflow 1 maupun inflow 2 dibuat merapat ke permukaan dinding maupun pelat dengan pertimbangan bahwa parameter-parameter aliran fluida di dekat permukaan dinding ataupun pelat akan mengalami perubahan karakteristik fluida

dan epsilon () masing- masing dengan nilai konvergensi 0,001, sementara nilai konvergensi residual energi sebesar 1 e-6 (3). Dengan nilai-nilai residual tersebut di atas dan jumlah jaringan kisi sebesar 60 000 dapat dilihat bahwa solusi akan konvergen setelah 395 iterasi, seperti dilihat pada Gambar 4.

Gambar 3. Bentuk jaringan kisi yang terjadi.



ISSN 0216 - 3128

48

Gambar 4. Residual vs iterasi. Gambar 5 menunjukkan vektor kecepatan , lapisan tipis di dekat permukaan dinding, di mana di mana kecepatan di daerah fluida bagian atas terjadi perubahan harga kecepatan fluida dari mempunyai arah yang berlawanan dengan hampir mendekati 0 (nol) pada permukaan dinding kecepatan di daerah fluida bagian bawah, sampai mencapai harga kecepatan aliran fluida sementara dibagian pelat tidak terdapat vektor bebas (dalam hal ini sebesar 4.5 m/s) pada daerah kecepatan. yang menjauh dari permukaan dinding. Lapisan tipis tersebut dikenal dengan nama lapisan batas. Dengan membuat jaringan kisi yang sangat Gambar 6 menunjukkan vektor kecepatan yang rapat di permukaan dinding, dapat dilihat berubah karena adanya efek viskositas. timbulnya efek viskositas yaitu terbentuknya

Gambar 5. Vektor kecepatan pada daerah aliran fluida.


49

ISSN 0216 - 3128


Gambar 6. Vektor kecepatan di dekat permukaan. Gambar 7 menunjukkan kontur distribusi temperatur yang berlainan, pada permukaan atas temperatur pada daerah fluida bagian atas, pelat dan bawah pelat. dan pada daerah fluida bagian bawah. Dengan jenis Gambar 8 menunjukkan gambar dua dimenmaterial pelat berupa baja dengan harga konduksional temperatur pelat pada x = 0. Dapat dilihat tivitas panas k = 16,27 W/m-K, harga koefisien bahwa temperatur pada dinding yang bersentuhan perpindahan panas pada tekanan konstan Cp = dengan fluida yang lebih panas akan mempunyai 502,48 J/kg-K dan harga massa jenis  = 8030 harga yang lebih tinggi dan selanjutnya akan 3 kg/m , dapat dilihat terjadinya perubahan menurun untuk daerah permukaan pelat yang temperatur pada pelat karena terdapatnya dua buah bersentuhan dengan fluida yang lebih dingin. aliran fluida yang berlawanan arah dengan harga

Gambar 7. Kontur temperatur.



ISSN 0216 - 3128

50

Gambar 8. Kurva temperatur pada pelat arah melintang terhadap aliran fluida. Gambar 9 dan Gambar 10 menunjukkan aliran fluida yang berlawanan, sementara distribusi temperatur dalam arah aliran fluida pada temperatur fluida yang jauh dari permukaan pelat daerah fluida bagian atas dan daerah fluida bagian hampir tidak terpengaruh. Harga temperatur aliran bawah. fluida pada daerah temperatur fluida yang lebih dingin, yang dekat permukaan pelat akan Dari Gambar 9 dan Gambar 10 dapat mengalami kenaikan, sebaliknya harga temperatur dilihat bahwa temperatur dalam arah aliran fluida aliran fluida pada daerah temperatur fluida yang pada daerah aliran fluida yang dekat dengan lebih panas, yang dekat permukaan pelat akan permukaan pelat akan terpengaruh oleh keberadaan mengalami sedikit penurunan.

Gambar 9. Distribusi temperatur arah x pada daerah aliran fluida atas.


51

ISSN 0216 - 3128


Gambar 10. Distribusi temperatur arah x pada daerah aliran fluida bawah.

KESIMPULAN Dari hasil analisis dengan perangkat lunak Fluent V5, telah dapat diketahui distribusi kecepatan dan temperatur pada dua buah aliran fluida yang mempunyai arah berlawanan dan terpisah oleh keberadaan sebuah pelat. Hasil analisis juga dapat menunjukkan distribusi temperatur pada pelat yang menjadi pemisah ke dua aliran fluida tersebut. Hasil analisis yang lain yaitu dapat dilihat efek viskositas pada kecepatan di dekat permukaan dinding.

DAFTAR PUSTAKA

TANYA JAWAB Endiah Puji H.  Bagaimana dasar pemilihan persamaan dasar laminar/turbulen apabila bilangan Re belum diketahui, ataukah Re dihitung sebagai input?  Dapatkah program Fluent digunakan untuk mengetahui kontur suhu pada dinding reaktor aral vertikal dan horizontal (sebagai pengganti program CFD)? V.I.S. Wardhani

1. ANONYMOUS, GAMBIT Tutorial Guide, Fluent Incorporated, Lebanon, May 1998. 2. ANONYMOUS, FLUENT / UNS & RAMPANT 4.2 User’s Guide vol.1, Lebanon, June 1997. 3. ANONYMOUS, FLUENT / UNS & RAMPANT 4.2 User’s Guide vol.2, Lebanon, June 1997. 4. J.P. HOLMAN, Perpindahan Kalor Edisi ke enam, Penerbit Erlangga, Indonesia, 1998.

 Untuk mengetahui aliran laminar/turbulen, dihitung terlebih dahulu bilangan Re kemudian di inputkan ke dalam Boundary Condition/kondisi batas.  Program fluent adalah untuk analisis CFD baik untuk 2D maupun 3D jadi kemungkinan dapat dipergunakan untuk mengetahui kontur T dinding reactor vertical/horizontal, tergantung cara memodelkannya saja.

5. FRANK KREITH, Prinsip-prinsip Perpindahan Panas, Penerbit Erlangga, Indonesia, 1999, 6. V. INDRIATI SRI WARDHANI, Karakteristik Termohidrolik Reaktor Penelitian Pada Daya 2 MW Dengan Fluks Panas Tak Seragam, Thesis Magister, ITB-Bandung, 1997.

Putu Sukmabuana  Beda code Relap dan Fluent.  Apakah bisa untuk aliran fluida 2 phase.  Untuk semua sirip HE atau sebagian.



ISSN 0216 - 3128

V.I.S. Wardhani  Code Relap bisa dipergunakan untuk analisis thermohidrolik dan neutronik tapi lebih efektif untuk analisis keselamatannya. Fluent menampilkan untuk semua problem yang berhubungan dengan CFD baik 2D/3D.

52

 Bisa untuk aliran fluida 1 phase maupun 2 phase.  Dalam analisis disini hanya diwakili oleh 1 buah pelat saja sudsh dspst mengkover hasil yang diharapkan.


ANALISIS PROSES PERPINDAHAN PANAS DUA DIMENSI PADA DUA DAERAH ALIRAN FLUIDA YANG TERPISAH OLEH SEBUAH PELAT DENGAN PERANGKAT LUNAK FLUENT V5

Recommend Documents