A Robust TVD-WENO Scheme for Conservation Laws

The ultimate goal of this article is to develop a robust and accurate numerical method for solving hyperbolic conservation laws in one and two dimensions. A hybrid numerical method, coupling a cheap fourth order total variation diminishing (TVD) scheme [1] for smooth region and a Robust seventh-order weighted non-oscillatory (WENO) scheme [2] near discontinuities, is considered. High order multi-resolution analysis is used to detect the high gradients regions of the numerical solution in order to capture the shocks with the WENO scheme, while the smooth regions are computed with fourth order total variation diminishing (TVD). For time integration, we use the third order TVD Runge-Kutta scheme. The accuracy of the resulting hybrid high order scheme is comparable with these of WENO, but with significant decrease of the CPU cost. Numerical demonstrates that the proposed scheme is comparable to the high order WENO scheme and superior to the fourth order TVD scheme. Our scheme has the added advantage of simplicity and computational efficiency. Numerical tests are presented which show the robustness and effectiveness of the proposed scheme.

• A. Abdalla
• A. Kaltayev

• WENO scheme
• TVD schemes
• smoothness indicators
• multi-resolution.

[1] J. Shi and E. F. Toro, "Fully discrete high resolution schemes for
hyperbolic conservation laws," Int. J. Num. Method Fluids., vol. 23,
PP. 241-269, 1996.
[2] Y. Shen and G. Zhay, "A Robust Seventh-order WENO Scheme and Its
Applications" AIAA Paper 2008-0757.
[3] X.D. Liu, S. Osher, and T. Chan, "Weighted essentially non-oscillatory
schemes" J. Comput. Phys., vol. 126, pp. 200-212, 1994.
[4] G.S. Jiang, and C.W. Shu,"Efficient implementation of weighted ENO
schemes" J. Comput. Phys., vol. 126, pp. 202-228, 1996.
[5] D.S. Balsara and C.-W. Shu, "Monotonicity Preserving weighted
essentially non-oscillatory schemes with increasingly high order of
accuracy," J.Comput.Phys., vol. 160, pp. 405-452, 2000.
[6] Z.F. Xu and C.-W. Shu, "Anti-diffusive flux corrections for high order
finite difference WENO schemes," J.Comput.Phys., vol. 205, pp. 458-
485, 2005.
[7] A. Harten, "Adaptive multi-resolution schemes for shock
computations," Comput Phys., vol. 115, PP. 319-338, 1994.
[8] M. A. Alves, P. Cruz, A. Mendes, F. D. Magalhaes, F. T. Pinho, and P.
J. Oliveira, "Adaptive multi-resolution approach for solution of
hyperbolic PDEs," Comput Methods Appl. Mech Eng., vol. 191,PP.
3909-3928, 2002.
[9] S. Gotlieb and C. W. Shu, "Total variation diminishing Runge-Kutta
schemes," Math Comp., vol. 67, PP. 73-85, 1998.
[10] A. Harten, High resolution schemes for hyperbolic conservation laws, J
Comput Phys., vol. 49, PP. 357-393, 1983.
[11] Glaister P. "An approximate linearized Riemann solver for Euler
equations for real gases". Journal of Computational Physics., vol. 74,
PP. 382-408, 1988.