Efficient Alternative Finite Difference WENO Schemes for Hyperbolic Conservation Laws

doi:10.1007/s42967-023-00360-z

Abstract

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.

Key words: Hyperbolic PDEs, Numerical schemes, Conservation laws, Finite difference Weighted Essentially Non-Oscillatory (FD-WENO)

CLC Number:

Dinshaw S. Balsara, Deepak Bhoriya, Chi-Wang Shu, Harish Kumar. Efficient Alternative Finite Difference WENO Schemes for Hyperbolic Conservation Laws[J]. Communications on Applied Mathematics and Computation, 2025, 7(6): 2189-2242.

TrendMD

References

[1] Anile, A.M.: Relativistic Fluids and Magneto-Fluids: with Applications in Astrophysics and Plasma Physics. Cambridge University Press, Cambridge (1989)
[2] Arbogast, T., Huang, C.-S., Zhao, X.: Accuracy of WENO and adaptive order WENO reconstructions for solving conservation laws. SIAM J. Numer. Anal. 56(3), 1818–1847 (2018)
[3] Balsara, D.S.: Total variation diminishing algorithm for adiabatic and isothermal magnetohydrodynamics. Astrophys. J. Suppl. 116, 133–153 (1998)
[4] Balsara, D.S.: Second order accurate schemes for magnetohydrodynamics with divergence-free reconstruction. Astrophys. J. Suppl. 151(1), 149–184 (2004)
[5] Balsara, D.S.: Self-adjusting, positivity preserving high order schemes for hydrodynamics and magnetohydrodynamics. J. Comput. Phys. 231, 7504–7517 (2012)
[6] Balsara, D.S., Bhoriya, D., Shu, C.-W., Kumar, H.: Efficient finite difference WENO scheme for hyperbolic systems with non-conservative products. Commun. Appl. Math. Comput. (2023). https://doi.org/10.1007/s42967-023-00275-9
[7] Balsara, D.S., Garain, S., Florinski, V., Boscheri, W.: An efficient class of WENO schemes with adaptive order for unstructured meshes. J. Comput. Phys. 404, 109062 (2020)
[8] Balsara, D.S., Garain, S., Shu, C.-W.: An efficient class of WENO schemes with adaptive order. J. Comput. Phys. 326, 780–804 (2016)
[9] Balsara, D.S., Kim, J.: A subluminal relativistic magnetohydrodynamics scheme with ADER-WENO predictor and multidimensional Riemann solver-based corrector. J. Comput. Phys. 312, 357–384 (2016)
[10] Balsara, D.S., Rumpf, T., Dumbser, M., Munz, C.-D.: Efficient, high-accuracy ADER-WENO schemes for hydrodynamics and divergence-free magnetohydrodynamics. J. Comput. Phys. 228, 2480 (2009)
[11] Balsara, D.S., Samantaray, S., Subramanian, S.: Efficient WENO-based prolongation strategies for divergence-preserving vector fields. Commun. Appl. Math. Comput. 5(1), 428–484 (2022)
[12] Balsara, D.S., Shu, C.-W.: Monotonicity preserving weighted non-oscillatory schemes with increasingly high order of accuracy. J. Comput. Phys. 160, 405–452 (2000)
[13] Berberich, J.P., Käppeli, R., Chandrashekar, P., Klingenberg, C.: High order discretely well-balanced methods for arbitrary hydrostatic atmospheres. Commun. Comput. Phys. 30(3), 666–708 (2021)
[14] Bhoriya, D., Kumar, H.: Entropy-stable schemes for relativistic hydrodynamics equations. Z. Angew. Math. Phys. 71, 1–29 (2020)
[15] Biswas, B., Kumar, H., Bhoriya, D.: Entropy stable discontinuous Galerkin schemes for the special relativistic hydrodynamics equations. Comput. Math. Appl. 112, 55–75 (2022)
[16] Biswas, B., Kumar, H., Yadav, A.: Entropy stable discontinuous Galerkin methods for ten-moment Gaussian closure equations. J. Comput. Phys. 431, 110148 (2021)
[17] Borges, R., Carmona, M., Costa, B., Don, W.S.: An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws. J. Comput. Phys. 227(6), 3101–3211 (2008)
[18] Boscheri, W., Balsara, D.S.: High order direct arbitrary-Lagrangian-Eulerian (ALE) $ {P_N}{P_M} $ schemes on unstructured meshes. J. Comput. Phys. 398, 108899 (2019)
[19] Carlini, E., Ferretti, R., Russo, G.: A weighted essentially nonoscillatory, large time-step scheme for Hamilton-Jacobi equations. SIAM J. Sci. Comput. 27, 1071–1091 (2005)
[20] Castro, M., Costa, B., Don, W.S.: High order weighted essentially non-oscillatory WENO-Z schemes for hyperbolic conservation laws. J. Comput. Phys. 230, 1766–1792 (2011)
[21] Colella, P., Woodward, P.R.: The piecewise parabolic method (PPM) for gas-dynamical simulations. J. Comput. Phys. 54(1), 174–201 (1984)
[22] Cravero, I., Semplice, M.: On the accuracy of WENO and CWENO reconstructions of third order on nonuniform meshes. J. Sci. Comput. 67(3), 1219–1246 (2016)
[23] Duan, J.M., Tang, H.Z.: High-order accurate entropy stable finite difference schemes for one- and two-dimensional special relativistic hydrodynamics. Adv. Appl. Math. Mech. 12(1), 1–29 (2019)
[24] Dumbser, M., Balsara, D.S.: A new efficient formulation of the HLLEM Riemann solver for general conservative and non-conservative hyperbolic systems. J. Comput. Phys. 304, 275–319 (2016)
[25] Gao, Z., Fang, L.-L., Wang, B.-S., Wang, Y., Don, W.S.: Seventh and ninth orders characteristic-wise alternative WENO finite difference schemes for hyperbolic conservation laws. Comput. Fluids 202, 104519 (2020)
[26] Gerolymos, G.A., Sénéchal, D., Vallet, I.: Very high order WENO schemes. J. Comput. Phys. 228, 8481–8524 (2009)
[27] Grosheintz-Laval, L., Käppeli, R.: Well-balanced finite volume schemes for nearly steady adiabatic flows. J. Comput. Phys. 423, 109805 (2020)
[28] Harten, A., Engquist, B., Osher, S., Chakravarthy, S.: Uniformly high order essentially non-oscillatory schemes III. J. Comput. Phys. 71, 231–303 (1987)
[29] He, P., Tang, H.: An adaptive moving mesh method for two-dimensional relativistic hydrodynamics. Commun. Comput. Phys. 11(1), 114–146 (2012)
[30] Henrick, A.K., Aslam, T.D., Powers, J.M.: Mapped weighted essentially non-oscillatory schemes: achieving optimal order near critical points. J. Comput. Phys. 207, 542–567 (2006)
[31] Jiang, G.-S., Shu, C.-W.: Efficient implementation of weighted ENO schemes. J. Comput. Phys. 126, 202–228 (1996)
[32] Jiang, Y., Shu, C.-W., Zhang, M.: An alternative formulation of finite difference ENO schemes with Lax-Wendroff time discretization for conservation laws. SIAM J. Sci. Comput. 35(2), A1137–A1160 (2013)
[33] Jiang, Y., Shu, C.-W., Zhang, M.: Free-stream preserving finite-difference schemes on curvilinear meshes. Methods Appl. Anal. 21(1), 001–030 (2014)
[34] Käppeli, R.: Well-balanced methods for computational astrophysics. Living Rev. Comput. Astrophys. 8(2), 1–88 (2022)
[35] Käppeli, R., Mishra, S.: Well-balanced schemes for the Euler equations with gravitation. J. Comput. Phys. 259, 199–219 (2014)
[36] Kumar, R., Chandrashekar, P.: Simple smoothness indicator and multi-level adaptive order WENO scheme for hyperbolic conservation laws. J. Comput. Phys. 375, 1059–1090 (2018)
[37] Kumar, R., Chandrashekar, P.: Efficient seventh order WENO schemes of adaptive order for hyperbolic conservation laws. Comput. Fluids 190, 49–76 (2019)
[38] Kupka, F., Happenhofer, N., Higueras, I., Koch, O.: Total-variation-diminishing implicit-explicit Runge-Kutta methods for the simulation of double-diffusive convection in astrophysics. J. Comput. Phys. 231, 3561–3586 (2012)
[39] Lax, P.D.: Weak solutions of nonlinear hyperbolic equations and their numerical computation. Commun. Pure Appl. Math. 7(1), 159–193 (1954)
[40] Levy, D., Puppo, G., Russo, G.: Compact central WENO schemes for multidimensional conservation laws. SIAM J. Sci. Comput. 22, 656–672 (2000)
[41] Ling, D., Duan, J.M., Tang, H.Z.: Physical-constraints-preserving Lagrangian finite volume schemes for one- and two-dimensional special relativistic hydrodynamics. J. Comput. Phys. 396, 507–543 (2019)
[42] Liska, R., Wendroff, B.: Comparison of several difference schemes on 1D and 2D test problems for the Euler equations. SIAM J. Sci. Comput. 25(3), 995–1017 (2003)
[43] Liu, X.-D., Osher, S., Chan, T.: Weighted essentially non-oscillatory schemes. J. Comput. Phys. 115, 200–212 (1994)
[44] Martí, M.J., Müller, E.: Numerical hydrodynamics in special relativity. Living Rev. Relat. 6(1), 7 (2003)
[45] Meena, A.K., Kumar, H., Chandrashekar, P.: Positivity-preserving high-order discontinuous Galerkin schemes for ten-moment Gaussian closure equations. J. Comput. Phys. 339, 370–395 (2017)
[46] Merriman, B.: Understanding the Shu-Osher conservative finite difference form. J. Sci. Comput. 19(1/2/3), 309 (2003)
[47] Mignone, A., Bodo, G.: An HLLC Riemann solver for relativistic flows I. Hydrodynamics. Mon. Notices R. Astron. Soc. 364(1), 126–136 (2005)
[48] Pao, S.P., Salas, M.D.: A numerical study of two-dimensional shock-vortex interaction. In: 14th Fluid and Plasma Dynamics Conference, AIAA 81-1205 (1981)
[49] Pareschi, L., Russo, G.: Implicit-explicit Runge-Kutta schemes and applications to hyperbolic systems with relaxation. J. Sci. Comput. 25, 129 (2005)
[50] Sebastian, K., Shu, C.-W.: Multidomain WENO finite difference method with interpolation at sub-domain interfaces. J. Sci. Comput. 19, 405–438 (2003)
[51] Semplice, M., Coco, A., Russo, G.: Adaptive mesh refinement for hyperbolic systems based on third-order compact WENO reconstruction. J. Sci. Comput. 66(2), 692–724 (2016)
[52] Sen, C., Kumar, H.: Entropy stable schemes for ten-moment Gaussian closure equations. J. Sci. Comput. 75(2), 1128–1155 (2018)
[53] Shu, C.-W.: High order weighted essentially non-oscillatory schemes for convection dominated problems. SIAM Rev. 51, 82–126 (2009)
[54] Shu, C.-W.: Essentially non-oscillatory and weighted essentially non-oscillatory schemes. Acta Numer. 29, 701–762 (2020)
[55] Shu, C.-W., Osher, S.J.: Efficient implementation of essentially non-oscillatory shock capturing schemes. J. Comput. Phys. 77, 439–471 (1988)
[56] Shu, C.-W., Osher, S.J.: Efficient implementation of essentially non-oscillatory shock capturing schemes II. J. Comput. Phys. 83, 32–78 (1989)
[57] Sod, G.A.: A survey of several finite difference methods for systems of nonlinear hyperbolic conservation laws. J. Comput. Phys. 27(1), 1–31 (1978)
[58] Spiteri, R.J., Ruuth, S.J.: A new class of optimal high-order strong-stability-preserving time-stepping schemes. SIAM J. Numer. Anal. 40, 469–491 (2002)
[59] Spiteri, R.J., Ruuth, S.J.: Non-linear evolution using optimal fourth-order strong-stability-preserving Runge-Kutta methods. Math. Comput. Simul. 62, 125–135 (2003)
[60] Woodward, P., Colella, P.: The numerical simulation of two-dimensional fluid flow with strong shocks. J. Comput. Phys. 54, 115–173 (1984)
[61] Zheng, F., Shu, C.-W., Qiu, J.: A high order conservative finite difference scheme for compressible two-medium flows. J. Comput. Phys. 445, 110597 (2021)
[62] Zhu, J., Qiu, J.: A new fifth order finite difference WENO scheme for solving hyperbolic conservation laws. J. Comput. Phys. 318, 110–121 (2016)
[63] Zhu, J., Shu, C.-W.: A new type of multi-resolution WENO schemes with increasingly higher order of accuracy. J. Comput. Phys. 375, 659–683 (2018)