Skema Diskritisasi dalam CFD: First Order vs Second Order
Skema Diskritisasi dalam CFD: First Order vs Second Order
Dalam Computational Fluid Dynamics (CFD), persamaan aliran fluida seperti persamaan Navier–Stokes tidak dapat langsung diselesaikan dalam bentuk diferensial kontinu. Oleh karena itu, persamaan tersebut harus diubah menjadi bentuk diskrit yang dapat dihitung oleh komputer. Proses ini dikenal sebagai diskritisasi.
Salah satu aspek penting dalam diskritisasi adalah skema orde (order of accuracy) yang digunakan untuk menghitung gradien, konveksi, maupun difusi pada grid komputasi. Dua skema yang paling umum digunakan adalah First Order dan Second Order.
Memahami perbedaan kedua metode ini sangat penting karena akan mempengaruhi akurasi, stabilitas, dan kebutuhan komputasi dari simulasi CFD.
Apa itu Skema Diskritisasi?
Diskritisasi adalah proses mengubah persamaan diferensial kontinu menjadi persamaan aljabar pada titik-titik grid atau mesh.
Sebagai contoh, turunan pertama suatu variabel dapat dituliskan dalam bentuk pendekatan numerik. Salah satu bentuk sederhana adalah forward difference:
Pendekatan seperti ini memiliki tingkat akurasi tertentu yang disebut order of accuracy. Order ini menunjukkan seberapa besar error yang muncul ketika ukuran grid (Δx) diperkecil.
Semakin tinggi order, semakin kecil error numerik yang dihasilkan.
First Order Scheme
Skema First Order merupakan pendekatan diskritisasi dengan tingkat akurasi orde satu. Error numeriknya berbanding lurus dengan ukuran grid.
Error∼O(Δx)
Dalam CFD, contoh skema first order yang paling sering digunakan adalah:
-
First Order Upwind Scheme
-
First Order Backward Difference
Pada First Order Upwind, nilai variabel pada sebuah face dihitung berdasarkan arah aliran fluida. Jika aliran datang dari kiri, maka nilai diambil dari sel sebelah kiri (upstream cell).
Karakteristik First Order
Beberapa karakteristik utama skema ini antara lain:
1. Stabilitas tinggi
Skema first order sangat stabil secara numerik sehingga sering digunakan pada tahap awal simulasi untuk membantu konvergensi.
2. Numerical diffusion besar
Kelemahan utama metode ini adalah munculnya numerical diffusion yang cukup besar. Hal ini menyebabkan gradien tajam pada aliran menjadi lebih “halus” dari kondisi sebenarnya.
Contohnya pada simulasi:
-
shock wave
-
boundary layer
-
shear layer
gradien fisik yang seharusnya tajam bisa terlihat lebih menyebar.
3. Cocok untuk simulasi awal
Banyak praktisi CFD menggunakan first order terlebih dahulu untuk mendapatkan solusi awal sebelum beralih ke skema yang lebih akurat.
Second Order Scheme
Skema Second Order memiliki akurasi yang lebih tinggi dibandingkan first order. Error numeriknya sebanding dengan kuadrat ukuran grid.
Error∼O(Δx2)
Artinya, ketika ukuran mesh diperkecil, error berkurang jauh lebih cepat dibandingkan metode first order.
Beberapa contoh skema second order dalam CFD:
-
Second Order Upwind
-
Central Difference Scheme
-
QUICK Scheme
Pada metode ini, nilai variabel pada face dihitung menggunakan lebih banyak titik grid sehingga gradien aliran dapat direpresentasikan dengan lebih akurat.
Karakteristik Second Order
1. Akurasi lebih tinggi
Second order mampu menangkap gradien aliran dengan lebih baik, terutama pada kasus:
-
aliran turbulen
-
shear layer
-
vorteks
-
shock
2. Numerical diffusion lebih kecil
Karena akurasinya lebih tinggi, fenomena numerical diffusion jauh lebih kecil dibandingkan first order.
3. Lebih sensitif terhadap mesh
Second order membutuhkan kualitas mesh yang lebih baik. Mesh yang buruk dapat menyebabkan:
-
osilasi numerik
-
ketidakstabilan simulasi
Pendekatan seperti ini memiliki tingkat akurasi tertentu yang disebut order of accuracy. Order ini menunjukkan seberapa besar error yang muncul ketika ukuran grid (Δx) diperkecil.
Semakin tinggi order, semakin kecil error numerik yang dihasilkan.
First Order Scheme
Skema First Order merupakan pendekatan diskritisasi dengan tingkat akurasi orde satu. Error numeriknya berbanding lurus dengan ukuran grid.
Error∼O(Δx)
Dalam CFD, contoh skema first order yang paling sering digunakan adalah:
-
First Order Upwind Scheme
-
First Order Backward Difference
Pada First Order Upwind, nilai variabel pada sebuah face dihitung berdasarkan arah aliran fluida. Jika aliran datang dari kiri, maka nilai diambil dari sel sebelah kiri (upstream cell).
Karakteristik First Order
Beberapa karakteristik utama skema ini antara lain:
1. Stabilitas tinggi
Skema first order sangat stabil secara numerik sehingga sering digunakan pada tahap awal simulasi untuk membantu konvergensi.
2. Numerical diffusion besar
Kelemahan utama metode ini adalah munculnya numerical diffusion yang cukup besar. Hal ini menyebabkan gradien tajam pada aliran menjadi lebih “halus” dari kondisi sebenarnya.
Contohnya pada simulasi:
-
shock wave
-
boundary layer
-
shear layer
gradien fisik yang seharusnya tajam bisa terlihat lebih menyebar.
3. Cocok untuk simulasi awal
Banyak praktisi CFD menggunakan first order terlebih dahulu untuk mendapatkan solusi awal sebelum beralih ke skema yang lebih akurat.
Second Order Scheme
Skema Second Order memiliki akurasi yang lebih tinggi dibandingkan first order. Error numeriknya sebanding dengan kuadrat ukuran grid.
Error∼O(Δx2)
Artinya, ketika ukuran mesh diperkecil, error berkurang jauh lebih cepat dibandingkan metode first order.
Beberapa contoh skema second order dalam CFD:
-
Second Order Upwind
-
Central Difference Scheme
-
QUICK Scheme
Pada metode ini, nilai variabel pada face dihitung menggunakan lebih banyak titik grid sehingga gradien aliran dapat direpresentasikan dengan lebih akurat.
Karakteristik Second Order
1. Akurasi lebih tinggi
Second order mampu menangkap gradien aliran dengan lebih baik, terutama pada kasus:
-
aliran turbulen
-
shear layer
-
vorteks
-
shock
2. Numerical diffusion lebih kecil
Karena akurasinya lebih tinggi, fenomena numerical diffusion jauh lebih kecil dibandingkan first order.
3. Lebih sensitif terhadap mesh
Second order membutuhkan kualitas mesh yang lebih baik. Mesh yang buruk dapat menyebabkan:
-
osilasi numerik
-
ketidakstabilan simulasi