Efficient Alternative Finite Difference WENO Schemes for Hyperbolic Conservation Laws

Abstract

Higher order finite difference Weighted Essentially Non-Oscillatory (FD-WENO) schemes for conservation laws are extremely popular because, for multidimensional problems, they offer high order accuracy at a fraction of the cost of finite volume WENO or DG schemes. Such schemes come in two formulations. The very popular classical FD-WENO method (Shu and Osher J Comput Phys 83: 32–78, 1989) relies on two reconstruction steps applied to two split fluxes. However, the method cannot accommodate different types of Riemann solvers and cannot preserve free stream boundary conditions on curvilinear meshes. This limits its utility. The alternative FD-WENO (AFD-WENO) method can overcome these deficiencies, however, much less work has been done on this method. The reasons are three-fold. First, it is difficult for the casual reader to understand the intricate logic that requires higher order derivatives of the fluxes to be evaluated at zone boundaries. The analytical methods for deriving the update equation for AFD-WENO schemes are somewhat recondite. To overcome that difficulty, we provide an easily accessible script that is based on a computer algebra system in Appendix A of this paper. Second, the method relies on interpolation rather than reconstruction, and WENO interpolation formulae have not been documented in the literature as thoroughly as WENO reconstruction formulae. In this paper, we explicitly provide all necessary WENO interpolation formulae that are needed for implementing the AFD-WENO up to the ninth order. The third reason is that the AFD-WENO requires higher order derivatives of the fluxes to be available at zone boundaries. Since those derivatives are usually obtained by finite differencing the zone-centered fluxes, they become susceptible to a Gibbs phenomenon when the solution is non-smooth. The inclusion of those fluxes is also crucially important for preserving the order property when the solution is smooth. This has limited the utility of the AFD-WENO in the past even though the method per se has many desirable features. Some efforts to mitigate the effect of finite differencing of the fluxes have been tried, but so far they have been done on a case by case basis for the PDE being considered. In this paper we find a general-purpose strategy that is based on a different type of the WENO interpolation. This new WENO interpolation takes the first derivatives of the fluxes at zone centers as its inputs and returns the requisite non-linearly hybridized higher order derivatives of flux-like terms at the zone boundaries as its output. With these three advances, we find that the AFD-WENO becomes a robust and general-purpose solution strategy for large classes of conservation laws. It allows any Riemann solver to be used. The AFD-WENO has a computational complexity that is entirely comparable to the classical FD-WENO, because it relies on two interpolation steps which cost the same as the two reconstruction steps in the classical FD-WENO. We apply the method to several stringent test problems drawn from Euler flow, relativistic hydrodynamics (RHD), and ten-moment equations. The method meets its design accuracy for smooth flow and can handle stringent problems in one and multiple dimensions.

Funding

National Science Foundation(DMS-2309249); National Aeronautics and Space Administration(80NSSC22K0628); Department of Science and Technology, Government of Rajasthan(VJR/2018/00129); University of Notre Dame(Travel Grant)