Higher-Order Finite Difference Scheme With TVD Filter Algorithm. Part 1: Simulations of Scalar Advective Dominant Problems.

Abstract

Engquist et al. (1989) proposed non-linear TVD filters. When combined with a traditional higher-order finite difference scheme, these filters can simulate shock or high concentration gradient problems with no spurious oscillations. Among the filter proposed is the filter algorithm 2.2. However, this algorithm flattens extrema that are not the results of overshooting and consequency the scheme reduces to a low order of accuracy locally around smooth extrema. Modification of the TVD filter algorithm 2.2 has been proposed in this paper to overcome this problem. Several conservative finite difference schemes are considered for testing the TVD filter. Non-conservative schemes consisting of 4th and 6th -order Runge-Kutta method are also evaluated. The modified filter has been tested to simulate seven test cases, including a pure advection of scalar profiles, a pure advection with variable velocity, two inviscid burger equations, an advection-diffusion equation with variable velocity and dispersion, advection of three solid bodies in rotating fluid around a square of a side length of 2, and a two-dimensional advection-diffusion equation. The numerical experiments showed that applying the modified TVD filter, combined with the higher-order non-TVD finite difference schemes for solving the advection equation, can produce accurate solutions with no oscillations and no clipping effect extrema. Keywords : Non-Oscillatory schemes, TVD filter, advection, higher-order accurate Abstrak Engquist dkk. (1989) mengusulkan filter TVD non-linear. Ketika dikombinasikan dengan skema beda hingga orde tinggi tradisional, filter ini dapat mensimulasikan kejutan atau masalah gradien konsentrasi tinggi tanpa osilasi palsu. Di antara filter yang diusulkan adalah algoritma filter 2.2. Namun, algoritma ini meratakan ekstrem yang bukan merupakan hasil dari “overshooting” dan konsekuensinya skema tersebut direduksi ke tingkat akurasi yang rendah secara lokal di sekitar ekstrem halus. Modifikasi algoritma filter TVD 2.2 telah diusulkan dalam makalah ini untuk mengatasi masalah ini. Beberapa skema konservatif dipertimbangkan untuk menguji filter TVD. Skema non-konservatif metode Runge-Kutta orde ke-4 dan ke-6 juga dievaluasi. Filter yang dimodifikasi diuji untuk mensimulasikan tujuh kasus, termasuk adveksi murni profil skalar, adveksi murni dengan kecepatan bervariasi, dua persamaan burger inviscid, persamaan adveksi-difusi dengan kecepatan bervariasi dan dispersi, adveksi tiga benda padat berputar oleh aliran di bidang persegi dengan panjang sisi 2, dan persamaan adveksi-difusi dua dimensi. Eksperimen numerik menunjukkan bahwa penerapan filter TVD yang dimodifikasi, dikombinasikan dengan skema beda hingga non-TVD orde tinggi untuk menyelesaikan persamaan adveksi, dapat menghasilkan solusi yang akurat tanpa osilasi dan tanpa efek kliping Kata Kunci : skema tanpa osilasi, Filter TVD, advection, orde-akurasi tinggi