Numerical Dispersion dalam CFD
Dalam Computational Fluid Dynamics (CFD), proses penyelesaian persamaan aliran fluida dilakukan dengan metode numerik yang mendiskritisasi ruang dan waktu. Meskipun metode ini memungkinkan komputer menghitung fenomena aliran yang kompleks, proses diskritisasi juga dapat menimbulkan berbagai error numerik.
Selain numerical diffusion, salah satu jenis error yang sering muncul adalah numerical dispersion. Fenomena ini dapat menyebabkan solusi simulasi menunjukkan osilasi atau gelombang yang tidak realistis, terutama pada daerah yang memiliki perubahan variabel yang tajam.
Memahami numerical dispersion sangat penting karena error ini dapat mempengaruhi kualitas hasil simulasi CFD, khususnya pada kasus yang melibatkan transport gelombang, shock wave, atau vorteks.
Pengertian Numerical Dispersion
Numerical dispersion adalah error numerik yang menyebabkan perubahan fase atau posisi gelombang dalam solusi numerik dibandingkan dengan solusi fisik yang sebenarnya.
Dalam simulasi CFD, numerical dispersion biasanya muncul ketika metode diskritisasi menghasilkan kesalahan dalam memprediksi kecepatan propagasi gelombang.
Akibatnya, gelombang atau variasi dalam variabel aliran dapat:
-
bergerak dengan kecepatan yang tidak tepat
-
mengalami distorsi bentuk
-
menampilkan osilasi yang tidak realistis
Fenomena ini paling sering terlihat pada simulasi yang melibatkan transport konveksi atau propagasi gelombang.
Intuisi Sederhana Numerical Dispersion
Bayangkan suatu gelombang tekanan yang bergerak dalam aliran fluida.
Secara fisik, gelombang tersebut akan bergerak dengan kecepatan tertentu dan mempertahankan bentuknya dalam jarak tertentu.
Namun dalam simulasi CFD, metode numerik dapat menyebabkan gelombang tersebut:
-
mengalami pergeseran posisi
-
memiliki fase yang tidak tepat
-
menampilkan osilasi tambahan di sekitar gradien tajam
Osilasi ini sering terlihat sebagai gelombang kecil di sekitar perubahan variabel, meskipun secara fisik fenomena tersebut tidak ada.
Hubungan dengan Persamaan Transport
Numerical dispersion biasanya muncul ketika menyelesaikan persamaan transport konveksi seperti:
Persamaan ini menggambarkan transport suatu variabel ϕ oleh aliran dengan kecepatan uu.
Ketika turunan spasial dalam persamaan tersebut dihitung menggunakan metode diskritisasi tertentu, error dalam aproksimasi turunan dapat menyebabkan perubahan fase gelombang dalam solusi numerik.
Error inilah yang dikenal sebagai numerical dispersion.
Perbedaan Numerical Dispersion dan Numerical Diffusion
Numerical dispersion sering dibandingkan dengan numerical diffusion, karena keduanya merupakan error numerik yang muncul akibat diskritisasi.
Namun kedua fenomena ini memiliki karakteristik yang berbeda.
Numerical diffusion menyebabkan gradien variabel menjadi lebih halus daripada kondisi fisik yang sebenarnya.
Sebaliknya, numerical dispersion menyebabkan osilasi atau gelombang palsu dalam solusi numerik.
Perbedaan utama keduanya dapat diringkas sebagai berikut:
Numerical diffusion
-
menyebabkan solusi menjadi terlalu halus
-
gradien tajam menjadi menyebar
Numerical dispersion
-
menyebabkan osilasi numerik
-
menghasilkan gelombang yang tidak realistis
Dalam banyak kasus simulasi CFD, kedua error ini dapat muncul secara bersamaan tergantung pada metode numerik yang digunakan.
Dampak Numerical Dispersion dalam Simulasi CFD
Numerical dispersion dapat mempengaruhi hasil simulasi dalam berbagai cara.
Beberapa contoh dampaknya antara lain:
Munculnya osilasi di sekitar shock wave
Pada simulasi aliran kompresibel, numerical dispersion dapat menghasilkan osilasi kecil di sekitar shock wave yang tidak ada dalam solusi fisik.
Distorsi gelombang tekanan
Dalam simulasi akustik atau propagasi gelombang tekanan, numerical dispersion dapat menyebabkan gelombang bergerak dengan kecepatan yang tidak tepat.
Ketidakakuratan pada aliran vorteks
Pada simulasi vorteks atau aliran turbulen, error fase dapat menyebabkan struktur vorteks berubah bentuk secara tidak realistis.
Akibatnya, simulasi tidak lagi menggambarkan fenomena fisik secara akurat.
Hubungan dengan Skema Diskritisasi
Besarnya numerical dispersion sangat dipengaruhi oleh skema diskritisasi yang digunakan untuk menghitung suku konveksi.
Beberapa skema memiliki kecenderungan menghasilkan dispersion error yang lebih besar dibandingkan yang lain.
Sebagai contoh:
Central Difference Scheme
Metode ini memiliki numerical diffusion yang rendah, tetapi sering menghasilkan numerical dispersion berupa osilasi pada daerah gradien tajam.
Higher-order schemes
Skema dengan orde tinggi dapat menghasilkan dispersion error jika tidak dilengkapi dengan mekanisme pembatas (limiter).
Upwind schemes
Metode upwind cenderung menghasilkan numerical diffusion yang lebih besar, tetapi biasanya lebih stabil dan memiliki osilasi yang lebih kecil.
Karena itu terdapat trade-off antara diffusion dan dispersion dalam pemilihan metode numerik.
Cara Mengurangi Numerical Dispersion
Dalam praktik simulasi CFD, terdapat beberapa pendekatan yang digunakan untuk mengurangi efek numerical dispersion.
Salah satu cara adalah menggunakan skema konveksi yang memiliki limiter, seperti metode TVD (Total Variation Diminishing). Metode ini dirancang untuk menjaga stabilitas sekaligus menghindari osilasi numerik.
Pendekatan lain adalah meningkatkan kualitas mesh atau menggunakan mesh yang lebih halus sehingga gradien variabel dapat direpresentasikan dengan lebih akurat.
Selain itu, beberapa solver menggunakan bounded schemes yang membatasi nilai variabel agar tidak menghasilkan osilasi yang tidak realistis.
Pentingnya Memahami Numerical Dispersion
Bagi pengguna CFD, memahami numerical dispersion sangat penting karena fenomena ini dapat mempengaruhi interpretasi hasil simulasi.
Tanpa pemahaman yang baik, osilasi numerik dapat disalahartikan sebagai fenomena fisik yang sebenarnya tidak terjadi.
Hal ini terutama penting pada simulasi yang melibatkan:
-
aliran kompresibel
-
propagasi gelombang
-
shock wave
-
akustik
-
transport skalar
Dalam kasus-kasus tersebut, pemilihan skema numerik yang tepat menjadi faktor kunci dalam menghasilkan solusi yang akurat.
Kesimpulan
Numerical dispersion merupakan salah satu jenis error numerik dalam simulasi CFD yang menyebabkan perubahan fase gelombang atau munculnya osilasi yang tidak realistis dalam solusi numerik.
Fenomena ini biasanya muncul akibat keterbatasan metode diskritisasi dalam merepresentasikan transport gelombang secara akurat.
Berbeda dengan numerical diffusion yang membuat gradien menjadi lebih halus, numerical dispersion menghasilkan osilasi atau gelombang palsu dalam hasil simulasi.
Dengan memahami sumber dan karakteristik numerical dispersion, pengguna CFD dapat memilih metode diskritisasi yang tepat sehingga simulasi menjadi lebih stabil, akurat, dan representatif terhadap fenomena fisik yang sebenarnya.