Efficient Alternative Finite Difference WENO Schemes for Hyperbolic Systems with Non-conservative Products

Abstract

Higher order finite difference Weighted Essentially Non-oscillatory (WENO) schemes for conservation laws represent a technology that has been reasonably consolidated. They are extremely popular because, when applied to multidimensional problems, they offer high order accuracy at a fraction of the cost of finite volume WENO or DG schemes. They come in two flavors. There is the classical finite difference WENO (FD-WENO) method (Shu and Osher in J. Comput. Phys. 83: 32–78, 1989). However, in recent years there is also an alternative finite difference WENO (AFD-WENO) method which has recently been formalized into a very useful general-purpose algorithm for conservation laws (Balsara et al. in Efficient alternative finite difference WENO schemes for hyperbolic conservation laws, submitted to CAMC, 2023). However, the FD-WENO algorithm has only very recently been formulated for hyperbolic systems with non-conservative products (Balsara et al. in Efficient finite difference WENO scheme for hyperbolic systems with non-conservative products, to appear CAMC, 2023). In this paper, we show that there are substantial advantages in obtaining an AFD-WENO algorithm for hyperbolic systems with non-conservative products. Such an algorithm is documented in this paper. We present an AFD-WENO formulation in a fluctuation form that is carefully engineered to retrieve the flux form when that is warranted and nevertheless extends to non-conservative products. The method is flexible because it allows any Riemann solver to be used. The formulation we arrive at is such that when non-conservative products are absent it reverts exactly to the formulation in the second citation above which is in the exact flux conservation form. The ability to transition to a precise conservation form when non-conservative products are absent ensures, via the Lax-Wendroff theorem, that shock locations will be exactly captured by the method. We present two formulations of AFD-WENO that can be used with hyperbolic systems with non-conservative products and stiff source terms with slightly differing computational complexities. The speeds of our new AFD-WENO schemes are compared to the speed of the classical FD-WENO algorithm from the first of the above-cited papers. At all orders, AFD-WENO outperforms FD-WENO. We also show a very desirable result that higher order variants of AFD-WENO schemes do not cost that much more than their lower order variants. This is because the larger number of floating point operations associated with larger stencils is almost very efficiently amortized by the CPU when the AFD-WENO code is designed to be cache friendly. This should have great, and very beneficial, implications for the role of our AFD-WENO schemes in the Peta- and Exascale computing. We apply the method to several stringent test problems drawn from the Baer-Nunziato system, two-layer shallow water equations, and the multicomponent debris flow. The method meets its design accuracy for the smooth flow and can handle stringent problems in one and multiple dimensions. Because of the pointwise nature of its update, AFD-WENO for hyperbolic systems with non-conservative products is also shown to be a very efficient performer on problems with stiff source terms.

Keywords

Hyperbolic PDEs / Numerical schemes / Conservation laws / Stiff source terms / Finite difference WENO / 65M06 / 65M22 / 65N06 / 35Q35 / 76T10 / 76M20

References

[1] Andrianov N, Warnecke G. The Riemann problem for the Baer-Nunziato two-phase flow model. J. Comput. Phys., 2004, 212: 434-464

[2] Baer MR, Nunziato JW. A two-phase mixture theory for the deflagration-to-detonation transition (DDT) in reactive granular materials. Int. J. Multiph. Flow, 1986, 12: 861-889

[3] Balsara DS. Total variation diminishing algorithm for adiabatic and isothermal magnetohydrodynamics. Ap. J. Supp., 1998, 116: 133-153

[4] Balsara DS. Self-adjusting, positivity preserving high order schemes for hydrodynamics and magnetohydrodynamics. J. Comput. Phys., 2012, 231: 7504-7517

[5] Balsara DS, Bhoriya D, Shu C-W, Kumar H. Efficient finite difference WENO scheme for hyperbolic systems with non-conservative products. Commun. Appl. Math. Comput., 2023

[6] Balsara, D.S., Bhoriya, D., Shu, C.-W., Kumar, H.: Efficient alternative finite difference WENO schemes for hyperbolic conservation laws. (2024). arxiv:2403.01264v1

[7] Balsara DS, Garain S, Florinski V, Boscheri W. An efficient class of WENO schemes with adaptive order for unstructured meshes. J. Comput. Phys., 2020, 404: 109062

[8] Balsara DS, Garain S, Shu C-W. An efficient class of WENO schemes with adaptive order. J. Comput. Phys., 2016, 326: 780-804

[9] Balsara DS, Rumpf T, Dumbser M, Munz C-D. Efficient, high-accuracy ADER-WENO schemes for hydrodynamics and divergence-free magnetohydrodynamics. J. Comput. Phys., 2009, 228: 2480

[10] Balsara DS, Samantaray S, Subramanian S. Efficient WENO-based prolongation strategies for divergence-preserving vector fields. Commun. Appl. Math. Comput., 2023, 5(1): 428-484

[11] Balsara DS, Shu C-W. Monotonicity preserving weighted non-oscillatory schemes with increasingly high order of accuracy. J. Comput. Phys., 2000, 160: 405-452

[12] Borges R, Carmona M, Costa B, Don WS. An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws. J. Comput. Phys., 2008, 227(63101-3211

[13] Castro M, Costa B, Don WS. High order weighted essentially non-oscillatory WENO-Z schemes for hyperbolic conservation laws. J. Comput. Phys., 2011, 230: 1766-1792

[14] Castro M, Gallardo JM, Parés C. High order finite volume schemes based on reconstruction of states for solving hyperbolic systems with nonconservative products. Applications to shallow-water systems. Math. Comput., 2006, 75(255): 1103-1134

[15] Castro MJ, Pardo A, Parés C, Toro EF. On some fast well-balanced first order solvers for nonconservative systems. Math. Comput., 2010, 79: 1427-1472

[16] Chiochetti, S., Müller, C.: A solver for stiff finite-rate relaxation in Baer-Nunziato two-phase flow models. (2020). arxiv: 2002.09381

[17] Colella P, Woodward PR. The piecewise parabolic method (PPM) for gas-dynamical simulations. J. Comput. Phys., 1984, 54(1): 174-201

[18] Coquel F, Marmignon C, Rai P, Renac F. An entropy stable high-order discontinuous Galerkin spectral element method for the Baer-Nunziato two-phase flow model. J. Comput. Phys., 2021, 431: 110135

[19] Dumbser M, Balsara DS. A new efficient formulation of the HLLEM Riemann solver for general conservative and non-conservative hyperbolic systems. J. Comput. Phys., 2016, 304: 275-319

[20] Dumbser M, Boscheri W. High-order unstructured Lagrangian one-step WENO finite volume schemes for non-conservative hyperbolic systems: applications to compressible multi-phase flows. Comput. Fluids, 2013, 86: 405-432

[21] Dumbser M, Castro M, Parés C, Toro EF. ADER schemes on unstructured meshes for non-conservative hyperbolic systems: applications to geophysical flows. Comput. Fluids, 2009, 38: 1731-1748

[22] Dumbser M, Hidalgo A, Castro M, Parés C, Toro EF. FORCE schemes on unstructured meshes II: non-conservative hyperbolic systems. Comput. Methods Appl. Mech. Engrg., 2010, 199: 625-647

[23] Dumbser M, Hidalgo A, Zanotti O. High-order space-time adaptive ADER-WENO finite volume schemes for non-conservative hyperbolic systems. Comput. Methods Appl. Mech. Engrg., 2014, 268: 359-387

[24] Dumbser M, Iben U, Ioriatti M. An efficient semi-implicit finite volume method for axially symmetric compressible flows in compliant tubes. Appl. Numer. Math., 2015, 89: 24-44

[25] Gao Z, Fang L-L, Wang B-S, Wang Y, Don WS. Seventh and ninth orders characteristic-wise alternative WENO finite difference schemes for hyperbolic conservation laws. Comput. Fluids, 2020, 202: 104519

[26] Gerolymos GA, Sénéchal D, Vallet I. Very high order WENO schemes. J. Comput. Phys., 2009, 228: 8481-8524

[27] Harten A, Engquist B, Osher S, Chakravarthy S. Uniformly high order essentially non-oscillatory schemes III. J. Comput. Phys., 1987, 71: 231-303

[28] Henrick AK, Aslam TD, Powers JM. Mapped weighted essentially non-oscillatoriy schemes: achieving optimal order near critical points. J. Comput. Phys., 2006, 207: 542-567

[29] Jiang G-S, Shu C-W. Efficient implementation of weighted ENO schemes. J. Comput. Phys., 1996, 126: 202-228

[30] 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., 2013, 35(2): A1137-A1160

[31] Jiang Y, Shu C-W, Zhang M. Free-stream preserving finite-difference schemes on curvilinear meshes. Methods Appl. Anal., 2014, 21(1): 001-030

[32] 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., 2012, 231: 3561-3586

[33] Liu X-D, Osher S, Chan T. Weighted essentially non-oscillatory schemes. J. Comput. Phys., 1994, 115: 200-212

[34] Merriman B. Understanding the Shu-Osher conservative finite difference form. J. Sci. Comput., 2003, 19(1/2/3309

[35] Pareschi L, Russo G. Implicit-explicit Runge-Kutta schemes and applications to hyperbolic systems with relaxation. J. Sci. Comput., 2005, 25: 129

[36] Pelanti M, Bouchut F, Mangeney A. A Roe-type scheme for two-phase shallow granular flows over variable topography. Math. Model. Numer. Anal., 2008, 42: 851-885

[37] Pitman EB, Le L. A two-fluid model for avalanche and debris flows. Philos. Trans. R. Soc. A, 2005, 363: 1573-1601

[38] Rhebergen S, Bokhove O, van der Vegt JJW. Discontinuous Galerkin finite element methods for hyperbolic nonconservative partial differential equations. J. Comput. Phys., 2008, 227: 1887-1922

[39] Saurel R, Abgrall R. A multiphase Godunov method for compressible multifluid and multiphase flows. J. Comput. Phys., 1999, 150: 425-467

[40] Schwendeman DW, Wahle CW, Kapila AK. The Riemann problem and a high-resolution Godunov method for a model of compressible two-phase flow. J. Comput. Phys., 2006, 212: 490-526

[41] Shu C-W. High order weighted essentially non-oscillatory schemes for convection dominated problems. SIAM Rev., 2009, 51: 82-126

[42] Shu C-W. Essentially non-oscillatory and weighted essentially non-oscillatory schemes. Acta Numer., 2020, 29: 701-762

[43] Shu C-W, Osher SJ. Efficient implementation of essentially non-oscillatory shock capturing schemes. J. Comput. Phys., 1988, 77: 439-471

[44] Shu C-W, Osher SJ. Efficient implementation of essentially non-oscillatory shock capturing schemes II. J. Comput. Phys., 1989, 83: 32-78

[45] Spiteri RJ, Ruuth SJ. A new class of optimal high-order strong-stability-preserving time-stepping schemes. SIAM J. Numer. Anal., 2002, 40: 469-491

[46] Spiteri RJ, Ruuth SJ. Non-linear evolution using optimal fourth-order strong-stability-preserving Runge-Kutta methods. Math. Comput. Simul, 2003, 62: 125-135

[47] Tokareva SA, Toro EF. HLLC-type Riemann solver for the Baer-Nunziato equations of compressible two-phase flow. J. Comput. Phys., 2010, 229: 3573-3604

[48] Woodward P, Colella P. The numerical simulation of two-dimensional fluid flow with strong shocks. J. Comput. Phys., 1984, 54: 115-173

[49] Zheng F, Shu C-W, Qiu J. A high order conservative finite difference scheme for compressible two-medium flows. J. Comput. Phys., 2021, 445: 110597

[50] Zhu J, Qiu J. A new fifth order finite difference WENO scheme for solving hyperbolic conservation laws. J. Comput. Phys., 2016, 318: 110-121

[51] Zhu J, Shu C-W. A new type of multi-resolution WENO schemes with increasingly higher order of accuracy. J. Comput. Phys., 2018, 375: 659-683

Funding

National Science Foundation(NSF-AST-2009776)

National Aeronautics and Space Administration(NASA-2020-1241)

DST(VJR/2018/00129)

NDI(Grant)

RIGHTS & PERMISSIONS

Shanghai University