Algoritma SIMPLE dalam CFD Simulation
Dalam Computational Fluid Dynamics (CFD), salah satu tantangan utama dalam menyelesaikan persamaan aliran fluida adalah bagaimana menangani hubungan antara kecepatan (velocity) dan tekanan (pressure). Pada aliran fluida inkompresibel, tekanan tidak memiliki persamaan transport langsung seperti kecepatan atau temperatur. Namun tekanan tetap harus dihitung dengan benar agar kondisi kontinuitas massa terpenuhi.
Untuk mengatasi masalah ini, digunakan suatu metode yang dikenal sebagai pressure–velocity coupling algorithm. Salah satu algoritma paling terkenal dan banyak digunakan dalam CFD adalah SIMPLE algorithm.
SIMPLE merupakan singkatan dari Semi-Implicit Method for Pressure Linked Equations, yang diperkenalkan oleh Suhas V. Patankar dan Brian Spalding pada tahun 1972. Algoritma ini menjadi dasar dari banyak solver CFD modern, termasuk pada software seperti OpenFOAM, ANSYS Fluent, dan berbagai kode CFD akademik.
Masalah Pressure–Velocity Coupling
Pada aliran inkompresibel, dua persamaan utama yang harus dipenuhi adalah:
-
Persamaan momentum
-
Persamaan kontinuitas
Persamaan momentum digunakan untuk menghitung kecepatan fluida, sedangkan persamaan kontinuitas memastikan bahwa massa fluida tetap terjaga.
Persamaan kontinuitas dapat dituliskan secara sederhana sebagai:
Persamaan ini menyatakan bahwa tidak ada akumulasi massa dalam suatu elemen kontrol.
Masalahnya adalah tekanan tidak muncul secara eksplisit dalam persamaan kontinuitas, tetapi tekanan sangat mempengaruhi kecepatan dalam persamaan momentum. Oleh karena itu diperlukan suatu metode untuk menyesuaikan tekanan dan kecepatan agar persamaan kontinuitas tetap terpenuhi.
Di sinilah algoritma SIMPLE digunakan.
Ide Dasar Algoritma SIMPLE
Konsep utama dari SIMPLE adalah melakukan iterasi antara tekanan dan kecepatan sampai solusi memenuhi persamaan momentum dan kontinuitas secara bersamaan.
Alih-alih langsung mencari solusi yang sempurna, SIMPLE menggunakan pendekatan bertahap melalui proses koreksi tekanan dan kecepatan.
Proses ini dimulai dengan tebakan awal untuk tekanan. Berdasarkan tekanan tersebut, persamaan momentum digunakan untuk menghitung medan kecepatan sementara.
Namun kecepatan yang diperoleh biasanya belum memenuhi persamaan kontinuitas. Oleh karena itu perlu dilakukan koreksi pada tekanan dan kecepatan.
Proses koreksi ini dilakukan berulang kali sampai kesalahan kontinuitas menjadi sangat kecil.
Langkah-Langkah Algoritma SIMPLE
Secara umum, algoritma SIMPLE berjalan melalui beberapa langkah iteratif.
1. Menebak tekanan awal
Simulasi dimulai dengan tebakan awal distribusi tekanan di seluruh domain.
2. Menyelesaikan persamaan momentum
Dengan tekanan yang diketahui, persamaan momentum diselesaikan untuk memperoleh kecepatan sementara.
Namun kecepatan ini biasanya belum memenuhi persamaan kontinuitas.
3. Menghitung koreksi tekanan
Dari kesalahan dalam persamaan kontinuitas, dibentuk suatu persamaan koreksi tekanan yang menunjukkan bagaimana tekanan harus diperbaiki.
4. Mengoreksi tekanan
Tekanan diperbarui menggunakan nilai koreksi yang diperoleh.
5. Mengoreksi kecepatan
Kecepatan juga diperbarui berdasarkan koreksi tekanan sehingga lebih mendekati kondisi yang memenuhi kontinuitas.
6. Mengulangi proses iterasi
Langkah-langkah tersebut diulang sampai solusi konvergen dan persamaan kontinuitas terpenuhi dengan baik.
Peran Under-Relaxation
Dalam implementasi praktis, algoritma SIMPLE sering menggunakan under-relaxation untuk menjaga stabilitas perhitungan.
Under-relaxation berarti perubahan tekanan atau kecepatan tidak langsung diterapkan sepenuhnya, tetapi hanya sebagian saja pada setiap iterasi.
Tujuannya adalah untuk menghindari perubahan solusi yang terlalu besar yang dapat menyebabkan simulasi menjadi tidak stabil.
Teknik ini sangat umum digunakan dalam solver CFD untuk memperbaiki konvergensi iterasi.
Kelebihan Algoritma SIMPLE
Algoritma SIMPLE memiliki beberapa keunggulan yang membuatnya sangat populer dalam CFD.
Salah satu keunggulan utama adalah stabilitas numeriknya yang baik, terutama untuk simulasi aliran steady-state.
Selain itu, metode ini relatif mudah diimplementasikan dalam metode numerik seperti Finite Volume Method, yang merupakan metode utama dalam banyak kode CFD.
SIMPLE juga cukup fleksibel dan dapat digunakan pada berbagai jenis aliran, mulai dari aliran laminar hingga turbulen.
Karena kelebihan tersebut, algoritma SIMPLE menjadi salah satu algoritma standar dalam banyak solver CFD.
Variasi dari Algoritma SIMPLE
Seiring perkembangan metode numerik, beberapa variasi dari algoritma SIMPLE telah dikembangkan untuk meningkatkan kecepatan konvergensi.
Beberapa variasi yang paling umum antara lain:
SIMPLEC (SIMPLE-Consistent)
Versi ini memodifikasi persamaan koreksi tekanan untuk mempercepat konvergensi.
SIMPLER
Metode ini menggunakan pendekatan yang lebih kuat dalam menghitung tekanan awal sehingga iterasi dapat menjadi lebih efisien.
PISO (Pressure Implicit with Splitting of Operators)
Algoritma ini sering digunakan untuk simulasi transien dan memiliki beberapa langkah koreksi tekanan dalam satu iterasi waktu.
Variasi-variasi ini membantu meningkatkan performa solver CFD pada berbagai jenis simulasi.
Implementasi dalam Software CFD
Algoritma SIMPLE digunakan secara luas dalam berbagai software CFD modern.
Pada OpenFOAM, algoritma SIMPLE digunakan dalam beberapa solver steady-state seperti:
-
simpleFoam
-
buoyantSimpleFoam
-
porousSimpleFoam
Solver ini menggunakan pendekatan SIMPLE untuk menyelesaikan coupling antara tekanan dan kecepatan pada aliran inkompresibel.
Pada software komersial seperti ANSYS Fluent, pengguna biasanya dapat memilih antara metode SIMPLE, SIMPLEC, atau PISO tergantung jenis simulasi yang dilakukan.
Kesimpulan
Algoritma SIMPLE merupakan salah satu metode paling penting dalam simulasi CFD untuk menangani hubungan antara tekanan dan kecepatan pada aliran fluida inkompresibel.
Metode ini bekerja dengan melakukan proses iteratif yang melibatkan perhitungan kecepatan sementara, pembentukan persamaan koreksi tekanan, serta pembaruan tekanan dan kecepatan hingga persamaan kontinuitas terpenuhi.
Karena stabilitas dan kemudahannya dalam implementasi, SIMPLE telah menjadi salah satu algoritma standar dalam banyak solver CFD modern.
Memahami cara kerja algoritma SIMPLE sangat penting bagi engineer dan peneliti yang ingin memahami lebih dalam bagaimana solver CFD menghitung aliran fluida secara numerik.