An Alternative Finite Difference WENO-Like Scheme with Physical Constraint Preservation for Divergence-Preserving Hyperbolic Systems

Abstract

Alternative finite difference Weighted Essentially Non-Oscillatory (AFD-WENO) schemes allow us to very efficiently update hyperbolic systems even in complex geometries. Recent innovations in AFD-WENO methods allow us to treat hyperbolic systems with non-conservative products almost as efficiently as conservation laws. However, some PDE systems, like computational electrodynamics (CED), magnetohydrodynamics (MHD), and relativistic magnetohydrodynamics (RMHD), have involution constraints that require divergence-free or divergence-preserving evolution of vector fields. In such situations, a Yee-style collocation of variables proves indispensable, and that collocation is retained in this work. In previous works, only higher order finite volume (FV) discretization of such involution-constrained systems was possible. In this work, we show that substantially more efficient AFD-WENO methods have been extended to encompass divergence-preserving hyperbolic PDEs. Our method retains the Yee-style collocation of normal components of the divergence-free/preserving vector field. However, the variables that require zone-centered evolution are evolved with AFD-WENO methods. Since those variables make up the bulk of the primal variables for the PDE of interest, this results in substantial savings in computational complexity. Even the volumetric reconstruction of the divergence-free/preserving vector field is bypassed. Instead, we realize that any divergence-preserving update of a vector field must have a general form. This general form looks closely like the familiar induction equation that is well known in CED or MHD. We exploit the generality of that form to extract the edge-centered variables that are needed in the update of the facially averaged vector field components. The two-dimensional (2D) Riemann solver is used to provide us with multidimensionally stabilized versions of these update terms. The generality of our approach is demonstrated by the fact that problems in CED, MHD, and RMHD can all be solved by the same general AFD-WENO algorithm that is presented here. Spatial accuracies up to the ninth order of accuracy are demonstrated. Several stringent test problems from CED, MHD, and RMHD are shown. We also show that the algorithm takes well to the physical constraint-preserving (PCP) formulation of AFD-WENO schemes that was presented by the authors. The efficient and time-explicit PCP strategy for divergence-preserving PDEs that we have presented here extends the applicability of our method to very stringent MHD and RMHD problems.

Keywords

Hyperbolic PDEs / Divergence free numerical schemes / Physical contraint preservation (PCP) / Computational electrodynamics (CED) / Finite difference Weighted Essentially Non-Oscillatory (FD-WENO) methods / 35Q85 / 35Q75 / 35Q35 / 35Q61 / 35L75 / 35L40 / 85A30 / 65M06 / 65M22

References

[1] Abgrall R. J. Comput. Phys., 2018, 372: 640-666.

[2] Abgrall R. J. Comput. Phys., 2022, 453110955

[3] Alfvén, H.: On the existence of electromagnetic-hydrodynamic waves. Arkiv. Mat. Astron. Fys. 29B(2), 1–7 (1942)

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

[5] Anile AM. Relativistic Fluids and Magneto-Fluids, 1989, Cambridge. Cambridge University Press

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

[7] Balsara DS. Total variation diminishing scheme for adiabatic and isothermal magnetohydrodynamics. Astrophys. J. Suppl. Ser., 1998, 116(1): 133-153.

[8] Balsara DS. Divergence-free adaptive mesh refinement for magnetohydrodynamics. J. Comput. Phys., 2001, 174: 614-648.

[9] Balsara DS. Total variation diminishing scheme for relativistic magnetohydrodynamics. Astrophys. J. Suppl., 2001, 132: 83-101.

[10] Balsara DS. Second-order-accurate schemes for magnetohydrodynamics with divergence-free reconstruction. Astrophys. J. Supplement, 2004, 151: 149-184.

[11] Balsara DS. Divergence-free reconstruction of magnetic fields and WENO schemes for magnetohydrodynamics. J. Comput. Phys., 2009, 228: 5040-5056.

[12] Balsara DS. Multidimensional HLLE Riemann solver: application to Euler and magnetohydrodynamic flows. J. Comput. Phys., 2010, 229: 1970-1993.

[13] Balsara DS. A two-dimensional HLLC Riemann solver for conservation laws: application to Euler and magnetohydrodynamic flows. J. Comput. Phys., 2012, 231: 7476-7503.

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

[15] Balsara DS. Multidimensional Riemann problem with self-similar internal structure—part I—application to hyperbolic conservation laws on structured meshes. J. Comput. Phys., 2014, 277: 163-200.

[16] Balsara DS. Three dimensional HLL Riemann solver for structured meshes: application to Euler and MHD flow. J. Comput. Phys., 2015, 295: 1-23.

[17] Balsara DS, et al. . High order accuracy divergence-free scheme for the electrodynamics of relativistic plasmas with multidimensional Riemann solvers. J. Comput. Phys., 2016, 318: 169-200.

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

[19] Balsara DS, Bhoriya D, Shu C-W, Kumar H. Efficient alternative finite difference WENO schemes for hyperbolic conservation laws. Commun.Appl. Math. Comput., 2024

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

[21] Balsara, D.S., Bhoriya, D., Singh, C., Kumar, H., Kappeli, R., Gatti, F.: Physical constraint preserving higher order finite volume schemes for divergence-free astrophysical MHD and RMHD, accepted. Astrophys. J. (2025)

[22] Balsara DS, Bhoriya D, Zanotti O, Dumbser M. Well-balanced high order finite difference WENO schemes for a first-order Z4 formulation of the Einstein field equations, submitted. Astrophys. J. Suppl. Ser., 2024

[23] Balsara DS, Dumbser M. Multidimensional Riemann problem with self-similar internal structure—part II—application to hyperbolic conservation laws on unstructured meshes. J. Comput. Phys., 2015, 287: 269-292.

[24] Balsara DS, Dumbser M, Abgrall R. Multidimensional HLL and HLLC Riemann solvers for unstructured meshes—with application to Euler and MHD flows. J. Comput. Phys., 2014, 261: 172-208.

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

[26] Balsara DS, Käppeli R. Von Neumann stability analysis of globally divergence-free RKDG and PNPM schemes for the induction equation using multidimensional Riemann solvers. J. Comput. Phys., 2017, 336: 104-127.

[27] Balsara DS, Käppeli R. Von Neumann stability analysis of globally divergence-free RKDG and PNPM schemes for the Maxwell Equations. J. Comput. Phys., 2019, 376: 1108-1137.

[28] Balsara DS, Kim J. A subluminal relativistic magnetohydrodynamics scheme with ADER-WENO predictor and multidimensional Riemann solver-based corrector. J. Comput. Phys., 2016, 312: 357-384.

[29] Balsara DS, Kumar R, Chandrashekar P. Globally divergence-free DG schemes for ideal compressible MHD at all orders. Commun. Appl. Math. Comput. Sci., 2021, 16(1): 59-98.

[30] Balsara DS, Meyer C, Dumbser M, Du H, Xu Z. Efficient implementation of ADER schemes for Euler and magnetohydrodynamical flows on structured meshes—comparison with Runge-Kutta methods. J. Comput. Phys., 2013, 235: 934-969.

[31] Balsara DS, Nkonga B. Formulating multidimensional Riemann solvers in similarity variables—part III : a multidimensional analogue of the HLLI Riemann solver for conservative hyperbolic systems. J. Comput. Phys., 2017, 346: 25-48.

[32] 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-2516.

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

[34] Balsara DS, Sarris CR. A systematic approach to adaptive mesh refinement for computational electrodynamics. IEEE J. Multiscale Multiphys. Comput.Tech., 2023, 8: 82-96.

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

[36] Balsara DS, Spicer DS. A staggered mesh algorithm using high order Godunov fluxes to ensure solenoidal magnetic fields in magnetohydrodynamic simulations. J. Comput. Phys., 1999, 149: 270-292.

[37] Balsara DS, Taflove A, Garain S, Montecinos G. Computational electrodynamics in material media with constraint-preservation, multidimensional Riemann solvers and sub-cell resolution – part I, second-order FVTD schemes. J. Comput. Phys., 2017, 349: 604-635.

[38] Balsara DS, Taflove A, Garain S, Montecinos G. Computational electrodynamics in material media with constraint-preservation, multidimensional Riemann solvers and sub-cell resolution – part II, higher-order FVTD schemes. J. Comput. Phys., 2018, 354: 613-645.

[39] Bhoriya D, Balsara DS, Chandrashekar P, Shu CW. Physical constraint preserving method for alternative finite difference WENO with non-conservative products. J. Sci. Comput., 2024

[40] Bhoriya D, Balsara DS, Florinski VA, Kumar H. Going beyond the MHD approximation: physics-based numerical solution of the CGL equations, to appear. Astrophys. J., 2024

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

[42] Buchmüller P, Dreher J, Helzel C. Improved accuracy of high-order WENO finite volume methods on Cartesian grids. J. Sci. Comput., 2014, 61: 343-368.

[43] Cai, C., Qiu, J., Wu, K.: Provably convergent and robust Newton-Raphson method: a new dawn in primitive variable recovery for relativistic MHD. arXiv:2404.05531 (2024)

[44] Castro M, Gallardo JE, Pares 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.

[45] Chew GF, Goldberger ML, Low FE. The Boltzmann equation and the one-fluid hydromagnetic equations in the absence of particle collisions. Proc. R. Soc. London Ser. A Math. Phys. Sci., 1956, 236: 112-118. DOI:

[46] Cockburn B, Hou S, Shu C-W. TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws IV: the multidimensional case. J. Comput. Phys., 1990, 54: 545-581

[47] Cockburn B, Lin SY, Shu C-W. TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws III: one-dimensional systems. J. Comput. Phys., 1989, 84: 90-113.

[48] Cockburn B, Shu CW. TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws II: general framework. Math. Comput., 1989, 52: 411-435. DOI:

[49] Cockburn B, Shu C-W. The Runge-Kutta discontinuous Galerkin method for conservation laws V: multidimensional systems. J. Comput. Phys., 1998, 141: 199-224.

[50] Cravero I, Semplice M. On the accuracy of WENO and CWENO reconstructions of third order on nonuniform meshes. J. Sci. Comput., 2016, 67(3): 1219-1246.

[51] Dai W, Woodward PR. On the divergence-free condition and conservation laws in numerical simulations for supersonic magnetohydrodynamic flows. Astrophys. J., 1998, 494: 317-335.

[52] Donnert JMF, Jang H, Mendygral P. WENO-WOMBAT: scalable fifth-order constrained-transport magnetohydrodynamics for astrophysical applications. Astrophys. J. Suppl. Ser., 2019, 241: 23.

[53] Dumbser M, Balsara DS, Toro EF, Munz C-D. A unified framework for the construction of quadrature-free one-step finite-volume and discontinuous galerkin schemes. J. Comput. Phys., 2008, 227: 8209-8253.

[54] Dumbser M, Zanotti O, Hidalgo A, Balsara DS. ADER-WENO finite volume schemes with space-time adaptive mesh refinement. J. Comput. Phys., 2013, 248: 257-286.

[55] Gurski, K.F.: An HLLC-type approximate Riemann solver for ideal magnetohydrodynamics. SIAM J. Sci. Comput. 25(6), 2165–2187 (2004)

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

[57] Hazra A, Chandrashekar P, Balsara DS. Globally constraint-preserving FR/DG scheme for Maxwell's equations at all orders. J. Comput. Phys., 2019, 394: 298-328.

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

[59] van der Holst B, Keppens R, Meliani Z. A multidimensional grid-adaptive relativistic magnetofluid code. Comput. Phys. Commun., 2008, 179: 617-627.

[60] Hu XY, Adams NA, Shu CW. Positivity-preserving method for high-order conservative schemes solving compressible Euler equations. J. Comput. Phys., 2013, 242: 169-180.

[61] Jeffrey A, Taniuti A. Nonlinear Wave Propagation, 1964, New York. Academic Press

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

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

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

[65] Komissarov SS. A Godunov-type scheme for relativistic magnetohydrodynamics. Monthly Not. R. Astron. Soc., 1999, 303: 343-366.

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

[67] Levy D, Puppo G, Russo G. Compact central WENO schemes for multidimensional conservation laws. SIAM J. Sci. Comput., 2000, 22: 656-672.

[68] McCorquodale P, Colella P. A high order finite volume method for conservation laws onlogically refined grids. Commun. App. Math. Comput. Sci., 2011, 6: 1.

[69] Mignone A, Bodo G. An HLLC Riemann solver for relativistic flows–II. Magnetohydrodynamics. Monthly Not. R. Astron. Soc., 2006, 368: 1040-1054.

[70] Orszag SA, Tang C-M. Small-scale structure of two-dimensional magnetohydrodynamic turbulence. J. Fluid Mech., 1979, 90(1): 129-143.

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

[72] Roe PL, Balsara DS. Notes on the eigensystem of magnetohydrodynamics. J. SIAM Appl. Math., 1996, 56: 57-67.

[73] Ryu D, Miniati F, Jones TW, Frank A. A divergence-free upwind code for multidimensional magnetohydrodynamic flows. Astrophys. J., 1998, 509: 244-255.

[74] Semplice M, Coco A, Russo G. Adaptive mesh refinement for hyperbolic systems based on third-order compact WENO reconstruction. J. Sci. Comput., 2016, 66(2): 692-724.

[75] Seo J, Ryu D. HOW-MHD: a high-order WENO-based magnetohydrodynamic code with a high-order constrained transport algorithm for astrophysical applications. Astrophys. J., 2023, 953: 39.

[76] Shi J, Hu C, Shu C-W. A technique for treating negative weights in WENO schemes. J. Comput. Phys., 2002, 175(108): 127

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

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

[79] Singh C, Yadav A, Bhoriya D, Kumar H, Balsara DS. Entropy stable finite difference schemes for Chew, Goldberger and Low anisotropic plasma flow equations. Commun. Appl. Math. Comput., 2024

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

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

[82] van der Vegt JJW, van der Ven H. Space–time discontinuous Galerkin finite element method with dynamic grid motion for inviscid compressible flows I. General formulation. J. Comput. Phys., 2002, 182: 546-585.

[83] van der Ven H, van der Vegt JJW. Space–time discontinuous Galerkin finite element method with dynamic grid motion for inviscid compressible flows II. Efficient flux quadrature. Comput. Methods Appl. Mech. Eng., 2002, 191: 4747-4780.

[84] K., Wu, & C.-W. Shu, . A provably positive discontinuous Galerkin method for multidimensional ideal magnetohydrodynamics. SIAM J. Sci. Comput., 2018, 40(5): 1302-1329.

[85] K., Wu, C.-W. Shu, . Entropy symmetrization and high-order accurate entropy stable numerical schemes for relativistic MHD equations. SIAM J. Sci. Comput., 2020, 42(4): 2230-2261.

[86] Xu, Z., Shu, C.-W.: A well-balanced conservative high-order alternative finite difference WENO (A-WENO) method for the shallow water equations. Adv. Water Resour. 196, 104898 (2025)

[87] Yee KS. Numerical solution of initial boundary value problems involving Maxwell equation in an isotropic media. IEEE Trans. Antenna Propagation, 1966, 14: 302.

[88] del Zanna L, Zanotti O, Bucciantini N, Londrillo P. ECHO: a Eulerian conservative high-order scheme for general relativistic magnetohydrodynamics and magnetodynamics. Astron. Astrophys., 2007, 473: 11.

[89] Zhang R, Zhang M, Shu C-W. On the order of accuracy and numerical performance of two classes of finite volume WENO schemes. Commun. Comput. Phys., 2011, 9(807): 827

[90] Zhu, J., Qiu, J.: A new fifth order finite difference WENO scheme for solving hyperbolic conservation laws. J. Comput. Phys. 38, 110–121 (2016)

[91] 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)

RIGHTS & PERMISSIONS

The Author(s)

Just Accepted

This article has successfully passed peer review and final editorial review, and will soon enter typesetting, proofreading and other publishing processes. The currently displayed version is the accepted final manuscript. The officially published version will be updated with format, DOI and citation information upon launch. We recommend that you pay attention to subsequent journal notifications and preferentially cite the officially published version. Thank you for your support and cooperation.