A Robust TVD-WENO Scheme for Conservation Laws

Authors: A. Abdalla, A. Kaltayev

Abstract:

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.

Keywords: WENO scheme, TVD schemes, smoothness indicators, multi-resolution.

Digital Object Identifier (DOI): doi.org/10.5281/zenodo.1328098

References:

[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.