Computational Fluid Dynamics (CFD) adalah sebuah analisis system yang melibatkan aliran fluida, perpindahan panas dan fenomena terkait seperti reaksi kimia dengan cara simulasi berbasis computer. Teknik ini sangat handal dan meliputi cakupan luas dalam area industry dan non industry. Beberapa contohnya yaitu :

· Aerodinamika pesawat dan kendaraan : lift dan drag

· Hidrodinamika kapal

· Pembangkit Tenaga : pembakaran dalam mesin IC dan turbin gas

· Mesin turbo : aliran dalam laluan rotating, diffuser dsb.

· Rekayasa electrical dan electronic : pendingin peralatan termasuk microchip

· Rekayasa proses kimia : mixing dan separation, polymer moulding

· Lingkungan internal dan eksternal gedung : beban angin dan pendinginan/veentilasi

· Rekayasa kelautan : beban struktur off-shore

· Rekayasa lingkungan : distribusi polutan dan anak sungai

· Oceanografi dan hidrologi : aliran sungai, muara, laut

· Meteorologi : prediksi cuaca

· Rekayasa biomedis : aliran darah melalui arteri dan vena

Keunggulan/keuntungan teknik analisis ini dibandingkan dengan pendekatan eksperimen dalam sebuah desain system fluida yaitu :

· Reduksi substansial waktu dan biaya untuk desain baru

· Kemampuan studi system yang tidak mampu dikontrol dengan eksperimen (misal system yang sangat luas)

· Kemampuan studi system dalam kondisi berbahaya pada dan di luar batas kinerja normal (seperti pada studi keselamatan dan scenario kecelakaan)

· Detil hasil yang lebih banyak secara praktis

Sistem analisis CFD tersusun dari beberapa algoritma numeric (disebut Code), yang terbagi oleh 3 elemen sebagai berikut :

Pre-processor

Merupakan bagian input suatu problem fluida ke sebuah program CFD melalui interface dan tranformasi lanjut ke dalam sebuah bentuk yang sesuai untuk solver. Langkah-langkah pengguna dalam tahap pre-processing yaitu :

· Definisi geometri region analisa : domain komputasional

· Pembuatan grid : pemecahan domain menjadi beberapa sub domain yang lebih kecil dan non overlapping : sebuah grid (mesh) atau volume atur/elemen

· Pemilihan fenomena fisik dan kimia yang perlu dimodelkan

· Definisi properties fluida

· Spesikasikan kondisi batas yang sesuai pada sel-sel yang berhimpit dengan batas domain

Solusi sebuah problem fluida (kecepatan, tekanan, temperature dsb) didefinisikan di setiap nodal di dalam masing-masing sel. Akurasi sebuah solusi CFD ditentukan oleh jumlah sel dalam grid. Secara umum, semakin besar jumlah sel semakin baik akurasi solusi. Baik akurasi solusi dan biaya hardaware computer serta lama kalkulasi tergantung kepada halusnya/rapatnya grid. Mesh-mesh optimal sering merupakan non-uniform : lebih rapat pada area di mana variasi-variasi banyak terjadi dari poin ke poin dan lebih jarang pada region dengan perubahan yang sedikit. Kemampuan teknik (self) adaptive meshing telah membantu pengembangan CFD guna otomatikal penghalusan grid untuk area dengan variasi yang padat. Sekitar 50% waktu proyek CFD di industry tercurah pada pendefinisian geometri domain dan penyusunan grid. Guna meningkatkan produktivitas pengguna code-code utama sekarang termasuk interface jenis CAD dan/atau fasilitas import data dari pemodelan surface dan meshing seperti PATRAN dan I-DEAS. Pre-prosesor hingga saat ini juga membantu user mengakses data library properties fluida umum dan fasilitas memasukkan model proses fisikal dan kimikal (model turbulence, perpindahan kalor radiatif, pembakaran) bersama persamaan aliran fluida utama.

Solver

Terdapat 3 macam teknik solusi numeric : beda hingga (finite difference), elemen hingga (finite element) dan metode spectral. Kerangka utama metode numeric untuk dasar sebuah solver terdiri dari langkah :

· Aproksimasi variable-variabel aliran yang tidak diketahui dengan fungsi-fungsi sederhana.

· Diskretisasi dengan substitusi aproksimasi ke dalam persamaan atur aliran dan manipulasi matematis lanjut.

· Solusi persamaan-persamaan aljabar. Perbedaan utama di antara ketiga macam teknik adalah pada cara aproksimasi variable-variabel aliran dan proses diskretisasi.

Metode Beda Hingga

Menggambarkan variable tidak diketahui Φ sebuah problem aliran dengan cara sampel-sampel titik pada titik-titik nodal sebuah grid dari garis koordinat. Ekspansi Deret Taylor terpotong sering dipakai untuk membangun aproksimasi-aproksimasi beda hingga derivative Φ dalam suku-suku sampel-sampel titik Φ di masing-masing titik grid dan tetangga terdekat. Derivatif tersebut muncul dalam persamaan atur digantikan oleh beda hingga menghasilkan persamaan aljabar untuk nilai-nilai Φ di setiap titik grid. Smith (1985) memberikan sebuah perhitungan komperhensif dari seluruh aspek metode beda hingga.

Metode Elemen Hingga

Menggunakan fungsi-fungsi potong (piecewise) sederhana (missal linier atau kuadratik) pada elemen-elemen untuk menggambarkan variasi-variasi local variable aliran yang tidak diketahui Φ. Persamaan atur terpenuhi secara tepat oleh solusi eksak Φ. Jika fungsi-fungsi aproksimasi potong untuk Φ disubstitusikan ke dalam persamaan, terdapat sebuah ketidak pastian hasil (residual) yang didefinisikan untuk mengukur kesalahan. Kemudian residual (kesalahan) diminimalkan melalui sebuah pengalian dengan sebuah set fungsi berbobot dan mengintegrasikannya. Hasilnya diperoleh sekumpulan persamaan aljabar untuk koefisien-koefisien tak diketahui dari fungsi-fungsi aproksimasi. Teori elemen hingga awalnya dikembangkan untuk analisis tegangan struktur. Untuk sebuah standar aplikasi fluida adalah Zienkiewics dan Taylor (1991).

Metode Spektral

Mengaproksimasikan variable Φ dengan deret Fourier terpotong atau deret Polinomial Chebyshev. Aproksimasi tidak secara local namun valid di semua domain komputasional, mengganti tak diketahui dalam persamaan atur dengan deret-deret terpotong. Batasan yang membawa ke persamaan aljabar untuk seluruh koefisien deret Fourier dan Chebyshev diberikan oleh konsep residual berbobot mirip dengan elemen hingga atau membuat fungsi aproksimasi serupa dengan solusi eksak pada sebuah nilai dari titik-titik grid. Informasi selebihnya dapat diperoleh dalam Gottlieb dan Orszag (1977).

Metode Volume Hingga (Finite Volume)

Awalnya dikembangkan untuk special formulasi beda hingga, algoritma numeric terdiri dari langkah :

· Intergrasi persamaan atur aliran fluida di seluruh volume atur (hingga) dari domain solusi

· Diskretisasi dengan substitusi beragam aproksimasi beda hingga untuk suku-suku persamaan terintegrasi proses aliran seperti konveksi, difusi dan sumber. Akan dikonversikan persamaan integral menjadi sebuah istem persamaan aljabar.

· Solusi persamaan-persamaan aljabar dengan metode iterative

Langkah awal, integrasi volume atur, membedakan metode volume hingga dari seluruh teknik CFD. Hasilnya menggambarkan konservasi (eksak) properties relevan di setiap sel ukuran hingga. Relasi yang jelas antara algoritma numeric dan prinsip konservasi fisis dasar memberikan sebuah ketertarikan dan konsep yang lebih mudah bagi para enginer.

Konservasi variable umum aliran Φ contohnya sebuah komponen kecepatan atau entalpi, dalam sebuah volume hingga dapat digambarkan sebagai keseimbangan di antara bermacam proses berkecenderungan menambah atau mengurangi. Dengan kata lain :

Post-processor

Hasil perhitungan modul solver berupa nilai-nilai numerik (angka-angka) variabel-variabel dasar aliran seperti komponen-komponen kecepatan, tekanan, temperatur dan fraksi-fraksi masa. Dalam modul post-processor nilai-nilai numerik ini diolah agar pengguna dapat dengan mudah membaca dan menganalisis hasil-hasil perhitungan CFD. Hasil-hasil ini dapat disajikan dalam bentuk grafis-grafis ataupun kontur-kontur distribusi parameter-parameter aliran fluida. Selain itu juga, modul post-processor menghitung parameter-parameter desain seperti koefisien gesek, Cd, Cl, Fluks panas , Gaya-gaya yang dikembangkan aliran fluida, Torsi, Daya dan lain sebagainya.