pith. sign in

arxiv: 2604.23885 · v1 · submitted 2026-04-26 · 🧮 math.NA · cs.NA· physics.comp-ph

A positivity preserving and entropy stable nodal discontinuous Galerkin scheme for ideal MHD equations

Pith reviewed 2026-05-08 05:23 UTC · model grok-4.3

classification 🧮 math.NA cs.NAphysics.comp-ph
keywords discontinuous Galerkinideal MHDpositivity preservingentropy stablelocally divergence-freeHLL fluxshock capturing
0
0 comments X

The pith

A nodal discontinuous Galerkin scheme for ideal MHD equations preserves positivity, entropy stability, and local divergence-freeness simultaneously.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a discontinuous Galerkin method for solving the ideal magnetohydrodynamic equations that addresses three longstanding numerical difficulties at once: keeping density and pressure positive, satisfying an entropy inequality, and controlling divergence errors in the magnetic field. It merges the locally divergence-free and positivity-preserving features of an earlier modal DG scheme with the entropy-stable properties of a nodal DG scheme. The combination is achieved through an HLL numerical flux that uses entropy-stable signal speed estimates together with a locally divergence-free projection operator; an essentially oscillation-free damping term is added to handle strong shocks. Numerical tests confirm that the resulting method maintains accuracy on smooth solutions and remains robust on discontinuous problems.

Core claim

The authors construct a nodal DG scheme for the ideal MHD system that is simultaneously positivity preserving, entropy stable, and locally divergence-free. The scheme uses an HLL flux equipped with entropy-stable signal speed estimates and a locally divergence-free projection; an essentially oscillation-free damping mechanism is included to treat strong shocks without destroying the three properties.

What carries the argument

HLL numerical flux with entropy-stable signal speed estimates combined with a locally divergence-free projection inside a nodal DG framework.

Load-bearing premise

The HLL flux with entropy-stable signal speeds and the locally divergence-free projection can be combined inside a nodal DG discretization while still preserving positivity, entropy stability, and the divergence-free property even after the oscillation-free damping is applied to strong shocks.

What would settle it

A simulation of a standard strong-shock MHD test problem in which the computed density or pressure becomes negative at any node, or in which the discrete entropy inequality is violated, would show the claim does not hold.

Figures

Figures reproduced from arXiv: 2604.23885 by Chi-Wang Shu, Yue Wu.

Figure 4.1
Figure 4.1. Figure 4.1: Contour plots of the Yee–Sj¨ogreen Riemann problem. 30 equally spaced contours view at source ↗
Figure 4.2
Figure 4.2. Figure 4.2: Contour plots of Orszag–Tang vortex at T = 0.5. 30 equally spaced contours for ρ ∈ [0.08, 0.5]. 14 view at source ↗
Figure 4.3
Figure 4.3. Figure 4.3: Evolution of entropy of Orszag–Tang vortex problem. view at source ↗
Figure 4.4
Figure 4.4. Figure 4.4: Contour plots of Mach number for the first rotor problem at || view at source ↗
Figure 4.5
Figure 4.5. Figure 4.5: Evolution of entropy of the first rotor problem. view at source ↗
Figure 4.6
Figure 4.6. Figure 4.6: Contour plots of the blast wave problem. 30 equally spaced contours for view at source ↗
Figure 4.7
Figure 4.7. Figure 4.7: Results of the cloud-shock interaction problem. For ln view at source ↗
Figure 4.8
Figure 4.8. Figure 4.8: Contour plots of the astrophysical jet example. 40 equally spaced contours for view at source ↗
read the original abstract

Numerically solving magnetohydrodynamic (MHD) equations faces many challenges: avoiding divergence error, maintaining positivity, and satisfying entropy conditions. Among discontinuous Galerkin (DG) schemes, there has been a modal version that is locally divergence-free and positivity preserving and a nodal version that is entropy stable. In this work, we develop a DG scheme that combines the advantages of these two and solves all the three challenges. The key ingredients that bring these two schemes together are an HLL numerical flux with entropy stable signal speed estimates and a locally divergence-free projection. To handle problems with strong shocks, the essentially oscillation-free damping is applied. Various numerical experiments verify the accuracy and robustness of our method.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript develops a nodal discontinuous Galerkin (DG) scheme for the ideal MHD equations that is positivity-preserving, entropy-stable, and locally divergence-free. It combines an HLL numerical flux with entropy-stable signal speed estimates and a locally divergence-free projection to merge advantages from prior modal (locally div-free and positivity-preserving) and nodal (entropy-stable) DG methods. An essentially oscillation-free damping term is added for robustness on strong shocks, with various numerical experiments used to verify accuracy and robustness.

Significance. If the simultaneous preservation of positivity, entropy stability, and exact local div-freeness is rigorously shown to hold with the added damping, the scheme would offer a practical nodal-basis alternative for MHD simulations that addresses the three core numerical challenges in a single framework. This could be useful for high-order computations in plasma and astrophysical applications, building directly on existing DG literature.

major comments (2)
  1. [Section describing the damping term and its integration with the nodal scheme] The description of the essentially oscillation-free damping (introduced for strong shocks) does not provide an explicit entropy balance or positivity estimate showing that its contribution remains compatible with the HLL flux and locally divergence-free projection. Without such control, the three properties cannot be guaranteed to hold simultaneously on strong-shock problems, which is central to the claim.
  2. [Numerical experiments section] The numerical experiments section reports verification of accuracy and robustness but does not include quantitative diagnostics (e.g., maximum divergence error norms, minimum density/pressure values, or entropy production rates) on strong-shock test cases that would confirm all three properties remain intact after damping is activated.
minor comments (2)
  1. Clarify the precise differences in implementation cost and accuracy between the new nodal scheme and the two referenced prior DG versions.
  2. Ensure consistent notation for the nodal basis functions, projection operator, and signal-speed estimates throughout the manuscript.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed comments. We address each major comment below, indicating where revisions will be made to strengthen the manuscript.

read point-by-point responses
  1. Referee: [Section describing the damping term and its integration with the nodal scheme] The description of the essentially oscillation-free damping (introduced for strong shocks) does not provide an explicit entropy balance or positivity estimate showing that its contribution remains compatible with the HLL flux and locally divergence-free projection. Without such control, the three properties cannot be guaranteed to hold simultaneously on strong-shock problems, which is central to the claim.

    Authors: We acknowledge that the manuscript does not derive an explicit entropy balance or positivity estimate for the damping term in combination with the entropy-stable HLL flux and locally divergence-free projection. The damping is constructed to be essentially oscillation-free and is applied in a post-processing step that respects the projection operator, but a complete rigorous proof of simultaneous preservation on strong shocks is not provided. In the revised version we will expand the description of the damping term to clarify its design and expected compatibility, and we will add a remark noting the absence of a full theoretical guarantee while emphasizing the supporting numerical evidence. revision: partial

  2. Referee: [Numerical experiments section] The numerical experiments section reports verification of accuracy and robustness but does not include quantitative diagnostics (e.g., maximum divergence error norms, minimum density/pressure values, or entropy production rates) on strong-shock test cases that would confirm all three properties remain intact after damping is activated.

    Authors: We agree that quantitative diagnostics would provide clearer confirmation. In the revised manuscript we will augment the numerical experiments section with tables reporting maximum divergence error norms, minimum density and pressure values, and entropy production rates for the strong-shock test cases. These diagnostics will be shown both with and without the damping term to verify that the three properties remain intact. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation combines independent prior DG components without self-referential reduction.

full rationale

The paper's central construction relies on an HLL flux with entropy-stable signal speeds and a locally divergence-free projection to merge a modal positivity-preserving scheme with a nodal entropy-stable scheme, plus an oscillation-free damping term for shocks. No equations or steps in the provided abstract or description reduce a claimed prediction or property to a fitted input, self-definition, or load-bearing self-citation chain by construction. The combination is presented as a novel synthesis of separately established techniques, with numerical experiments offered as verification rather than tautological confirmation. This satisfies the default expectation of a non-circular derivation.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Review based solely on abstract; no explicit free parameters, axioms, or invented entities are detailed. The scheme appears to rest on standard DG theory and hyperbolic conservation law assumptions.

axioms (1)
  • standard math Standard properties of discontinuous Galerkin discretizations and HLL-type numerical fluxes for hyperbolic systems
    The abstract invokes these as background for the new combination.

pith-pipeline@v0.9.0 · 5413 in / 1184 out tokens · 55751 ms · 2026-05-08T05:23:30.516306+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

69 extracted references · 69 canonical work pages

  1. [1]

    application to entropy conservative and entropy dissipative schemes, Journal of Computational Physics372(2018), 640–666

    R´ emi Abgrall,A general framework to construct schemes satisfying additional conservation relations. application to entropy conservative and entropy dissipative schemes, Journal of Computational Physics372(2018), 640–666

  2. [2]

    Balsara,Self-adjusting, positivity preserving high order schemes for hydrody- namics and magnetohydrodynamics, Journal of Computational Physics231(2012), no

    Dinshaw S. Balsara,Self-adjusting, positivity preserving high order schemes for hydrody- namics and magnetohydrodynamics, Journal of Computational Physics231(2012), no. 22, 7504–7517

  3. [3]

    Dinshaw S. Balsara, Rakesh Kumar, and Praveen Chandrashekar,Globally divergence-free DG scheme for ideal compressible MHD, Communications in Applied Mathematics and Computational Science16(2021), no. 1, 59–98

  4. [4]

    Balsara and Daniel S

    Dinshaw S. Balsara and Daniel S. Spicer,A staggered mesh algorithm using high order Godunov fluxes to ensure solenoidal magnetic fields in magnetohydrodynamic simulations, Journal of Computational Physics149(1999), no. 2, 270–292

  5. [5]

    5, Springer, Berlin, Heidelberg, 12 1998, pp

    Timothy John Barth,Numerical methods for gasdynamic systems on unstructured meshes, An Introduction to Recent Developments in Theory and Numerics for Conservation Laws (Dietmar Kr¨ oner, Mario Ohlberger, and Christian Rohde, eds.), Lecture Notes in Compu- tational Science and Engineering, vol. 5, Springer, Berlin, Heidelberg, 12 1998, pp. 195–285

  6. [6]

    Nicolaides, and Mikhail Jurievich Shashkov, eds.), The IMA Volumes in Mathematics and its Applica- tions, vol

    ,On the role of involutions in the discontinuous Galerkin discretization of Maxwell and magnetohydrodynamic systems, Compatible Spatial Discretizations (Douglas Norman Arnold, Pavel Blagovestov Bochev, Richard Bruno Lehoucq, Roy A. Nicolaides, and Mikhail Jurievich Shashkov, eds.), The IMA Volumes in Mathematics and its Applica- tions, vol. 142, Springer,...

  7. [7]

    Winters, Gregor J

    Marvin Bohm, Andrew R. Winters, Gregor J. Gassner, Dominik Derigs, Florian J. Hin- denlang, and Joachim Saur,An entropy stable nodal discontinuous Galerkin method for the resistive mhd equations. part I: Theory and numerical verification, Journal of Computa- tional Physics422(2020), 108076

  8. [8]

    I: theoretical framework, Numerische Mathematik108(2007), no

    Fran¸ cois Bouchut, Christian Klingenberg, and Knut Waagan,A multiwave approximate Riemann solver for ideal MHD based on relaxation. I: theoretical framework, Numerische Mathematik108(2007), no. 1, 7–42

  9. [9]

    4, 647–679

    ,A multiwave approximate Riemann solver for ideal MHD based on relaxation II: numerical implementation with 3 and 5 waves, Numerische Mathematik115(2010), no. 4, 647–679

  10. [10]

    3, 426–430

    Jeremiah Uhler Brackbill and Daniel Charles Barnes,The effect of nonzero∇ ·Bon the numerical solution of the magnetohydrodynamic equations, Journal of Computational Physics35(1980), no. 3, 426–430

  11. [11]

    Jesse Chan,On discretely entropy conservative and entropy stable discontinuous Galerkin methods, Journal of Computational Physics362(2018), 346–374

  12. [12]

    5, A2938–A2966

    Jesse Chan, David Cesar Del Rey Fern´ andez, and Mark Huitt Carpenter,Efficient entropy stable Gauss collocation methods, SIAM Journal on Scientific Computing41(2019), no. 5, A2938–A2966

  13. [13]

    2, 1313–1340

    Praveen Chandrashekar and Christian Klingenberg,Entropy stable finite volume scheme for ideal compressible MHD on 2-D Cartesian meshes, SIAM Journal on Numerical Analysis 54(2016), no. 2, 1313–1340. 19

  14. [14]

    Praveen Chandrashekar and Rakesh Kumar,Constraint preserving discontinuous Galerkin method for ideal compressible MHD on 2-D Cartesian grids, Journal of Scientific Computing 84(2020), no. 2, 39

  15. [15]

    Tianheng Chen and Chi-Wang Shu,Entropy stable high order discontinuous Galerkin meth- ods with suitable quadrature rules for hyperbolic conservation laws, Journal of Computa- tional Physics345(2017), 427–461

  16. [16]

    ,Review of entropy stable discontinuous Galerkin methods for systems of conserva- tion laws on unstructured simplex meshes, CSIAM Transactions on Applied Mathematics 1(2020), no. 1, 1–52

  17. [17]

    Yue Cheng, Fengyan Li, Jianxian Qiu, and Liwei Xu,Positivity-preserving DG and central DG methods for ideal MHD equations, Journal of Computational Physics238(2013), 255– 280

  18. [18]

    Bernardo Cockburn, George Em Karniadakis, and Chi-Wang Shu,Discontinuous Galerkin methods, 1 ed., Lecture Notes in Computational Science and Engineering, Springer, Berlin, Heidelberg, 2000

  19. [19]

    1, 90–113

    Bernardo Cockburn, San-Yih Lin, and Chi-Wang Shu,TVB Runge–Kutta local projection discontinuous Galerkin finite element method for conservation laws III: One-dimensional systems, Journal of Computational Physics84(1989), no. 1, 90–113

  20. [20]

    2, 199–224

    Bernardo Cockburn and Chi-Wang Shu,The Runge–Kutta discontinuous Galerkin method for conservation laws V: Multidimensional systems, Journal of Computational Physics141 (1998), no. 2, 199–224

  21. [21]

    Jared Crean, Jason Edward Hicken, David Cesar Del Rey Fern´ andez, David Walter Zingg, and Mark Huitt Carpenter,Entropy-stable summation-by-parts discretization of the Euler equations on general curved elements, Journal of Computational Physics356(2018), 410– 438

  22. [22]

    4, 373–389

    Constantine Michael Dafermos,Quasilinear hyperbolic systems with involutions, Archive for Rational Mechanics and Analysis94(1986), no. 4, 373–389

  23. [23]

    325, Springer, Berlin, Heidelberg, 03 2026

    ,Hyperbolic conservation laws in continuum physics, 5 ed., Grundlehren der math- ematischen Wissenschaften, vol. 325, Springer, Berlin, Heidelberg, 03 2026

  24. [24]

    Dahlburg and J

    Russel B. Dahlburg and J. Michael Picone,Evolution of the Orszag–Tang vortex system in a compressible medium. I. initial average subsonic flow, Physics of Fluids B: Plasma Physics1(1989), no. 11, 2153–2171

  25. [25]

    Woodward,A simple finite difference scheme for multidimen- sional magnetohydrodynamical equations, Journal of Computational Physics142(1998), no

    Wenlong Dai and Paul R. Woodward,A simple finite difference scheme for multidimen- sional magnetohydrodynamical equations, Journal of Computational Physics142(1998), no. 2, 331–369

  26. [26]

    Dedner, Friedemann Kemm, Dietmar Kr¨ oner, Claus-Dieter Munz, T

    Andreas S. Dedner, Friedemann Kemm, Dietmar Kr¨ oner, Claus-Dieter Munz, T. Schnitzer, and Matthias Wesenberg,Hyperbolic divergence cleaning for the MHD equations, Journal of Computational Physics175(2002), no. 2, 645–673

  27. [27]

    David Cesar Del Rey Fern´ andez, Jason Edward Hicken, and David Walter Zingg,Review of summation-by-parts operators with simultaneous approximation terms for the numerical solution of partial differential equations, Computers and Fluids95(2014), 171–196. 20

  28. [28]

    Winters, Gregor J

    Dominik Derigs, Andrew R. Winters, Gregor J. Gassner, Stefanie Walch, and Marvin Bohm,Ideal GLM-MHD: About the entropy consistent nine-wave magnetic field divergence diminishing ideal magnetohydrodynamics equations, Journal of Computational Physics364 (2018), 420–467

  29. [29]

    72, Springer, Cham, 02 2021

    Alexandre Ern and Jean-Luc Guermond,Finite elements I, Texts in Applied Mathematics, vol. 72, Springer, Cham, 02 2021

  30. [30]

    Evans and John F

    Charles R. Evans and John F. Hawley,Simulation of magnetohydrodynamic flows: A con- strained transport method, The Astrophysical journal332(1988), 659–677

  31. [31]

    Ulrik S. Fjordholm, Siddhartha Mishra, and Eitan Tadmor,Arbitrarily high-order accurate entropy stable essentially nonoscillatory schemes for systems of conservation laws, SIAM Journal on Numerical Analysis50(2012), no. 2, 544–573

  32. [32]

    3, 1621–1659

    Pei Fu, Fengyan Li, and Yan Xu,Globally divergence-free discontinuous Galerkin methods for ideal magnetohydrodynamic equations, Journal of Scientific Computing77(2018), no. 3, 1621–1659

  33. [33]

    Gardiner and James M

    Thomas A. Gardiner and James M. Stone,An unsplit Godunov method for ideal MHD via constrained transport, Journal of Computational Physics205(2005), no. 2, 509–539

  34. [34]

    1, 26–34

    Sergei Konstantinovich Godunov,Symmetric form of the equations of magnetohydrody- namics, Numerical Methods of Continuum Mechanics3(1972), no. 1, 26–34

  35. [35]

    221, 73–85

    Sigal Gottlieb and Chi-Wang Shu,Total vriation diminishing Runge-Kutta schemes, Math- ematics of Computation67(1998), no. 221, 73–85

  36. [36]

    1, 89–112

    Sigal Gottlieb, Chi-Wang Shu, and Eitan Tadmor,Strong stability-preserving high-order time discretization methods, SIAM Review43(2001), no. 1, 89–112

  37. [37]

    Jean-Luc Guermond and Bojan Popov,Fast estimation from above of the maximum wave speed in the Riemann problem for the Euler equations, Journal of Computational Physics 321(2016), 908–926

  38. [38]

    10, 3803–3829

    Christiane Helzel, Rossmanith James Alexander, and Bertram Taetz,An unstaggered con- strained transport method for the 3D ideal magnetohydrodynamic equations, Journal of Computational Physics230(2011), no. 10, 3803–3829

  39. [39]

    Hindenlang and Gregor J

    Florian J. Hindenlang and Gregor J. Gassner,A new entropy conservative two-point flux for ideal MHD equations derived from first principles, Talk presented at HONOM, 04 2019

  40. [40]

    5, 1766–1797

    Rossmanith James Alexander,An unstaggered, high-resolution constrained transport method for magnetohydrodynamic flows, SIAM Journal on Scientific Computing28(2006), no. 5, 1766–1797

  41. [41]

    Antony Jameson,Eigenvalues, eigenvectors and symmetrization of the magneto- hydrodynamic (MHD) equations, AFOSR Grantees and Contractors Meeting, 06 2006

  42. [42]

    2, 649–661

    Pekka Janhunen,A positive conservative method for magnetohydrodynamics based on HLL and Roe methods, Journal of Computational Physics160(2000), no. 2, 649–661

  43. [43]

    1, 202–228

    Guang-Shan Jiang and Chi-Wang Shu,Efficient implementation of weighted ENO schemes, Journal of Computational Physics126(1996), no. 1, 202–228

  44. [44]

    Friedemann Kemm,On the origin of divergence errors in MHD simulations and conse- quences for numerical schemes, Communications in Applied Mathematics and Computa- tional Science8(2013), no. 1, 1–38. 21

  45. [45]

    2, 423– 442

    Christian Klingenberg, Frank P¨ orner, and Yinhua Xia,An efficient implementation of the divergence free constraint in a discontinuous Galerkin method for magnetohydrodynamics on unstructured meshes, Communications in Computational Physics21(2017), no. 2, 423– 442

  46. [46]

    1, 413–442

    Fengyan Li and Chi-Wang Shu,Locally divergence-free discontinuous galerkin methods for MHD equations, Journal of Scientific Computing22(2005), no. 1, 413–442

  47. [47]

    Mengqing Liu and Kailiang Wu,Structure-preserving oscillation-eliminating discontinuous galerkin schemes for ideal MHD equations: Locally divergence-free and positivity-preserving, Journal of Computational Physics527(2025), 113795

  48. [48]

    1, A230–A259

    Yong Liu, Jianfang Lu, and Chi-Wang Shu,An essentially oscillation-free discontinuous Galerkin method for hyperbolic systems, SIAM Journal on Scientific Computing44(2022), no. 1, A230–A259

  49. [49]

    ,An entropy stable essentially oscillation-free discontinuous Galerkin method for solving ideal magnetohydrodynamic equations, Journal of Computational Physics530 (2025), 113911

  50. [50]

    Yong Liu, Chi-Wang Shu, and Mengping Zhang,Entropy stable high order discontinuous Galerkin methods for ideal compressible MHD on structured meshes, Journal of Computa- tional Physics354(2018), 163–178

  51. [51]

    353, 1147–1198

    Manting Peng, Zheng Sun, and Kailiang Wu,OEDG: oscillation eliminating discontinuous Galerkin method for hyperbolic conservation laws, Mathematics of Computation94(2025), no. 353, 1147–1198

  52. [52]

    Report ICASE-94-24, NASA Langley Research Center, Hampton, VA, 04 1994

    Kenneth Grant Powell,An approximate Riemann solver for magnetohydrodynamics, Tech. Report ICASE-94-24, NASA Langley Research Center, Hampton, VA, 04 1994

  53. [53]

    Roe, Timur J

    Kenneth Grant Powell, Philip L. Roe, Timur J. Linde, Tamas I. Gombosi, and Darren L. De Zeeuw,A solution-adaptive upwind scheme for ideal magnetohydrodynamics, Journal of Computational Physics154(1999), no. 2, 284–309

  54. [54]

    3, 907–929

    Jianxian Qiu and Chi-Wang Shu,Runge–Kutta discontinuous Galerkin method using WENO limiters, SIAM Journal on Scientific Computing26(2005), no. 3, 907–929

  55. [55]

    Andr´ es Mauricio Rueda-Ram´ ırez and Gregor J. Gassner,A flux-differencing formula for split-form summation by parts discretizations of non-conservative systems: Applications to subcell limiting for magneto-hydrodynamics, Journal of Computational Physics496(2024), 112607

  56. [56]

    Hindenlang, Jesse Chan, and Gregor J

    Andr´ es Mauricio Rueda-Ram´ ırez, Florian J. Hindenlang, Jesse Chan, and Gregor J. Gassner,Entropy-stable Gauss collocation methods for ideal magneto-hydrodynamics, Jour- nal of Computational Physics475(2023), 111851

  57. [57]

    Chi-Wang Shu,Essentially non-oscillatory and weighted essentially non-oscillatory schemes, Acta Numerica29(2020), 701–762

  58. [58]

    Eitan Tadmor,Entropy stability theory for difference approximations of nonlinear conser- vation laws and related time-dependent problems, Acta Numerica12(2003), 451–512

  59. [59]

    1, 25–34

    Eleuterio Francisco Toro, Michael Spruce, and William Speares,Restoration of the contact surface in the HLL-Riemann solver, Shock Waves4(1994), no. 1, 25–34. 22

  60. [60]

    2, 605–652

    G´ abor T´ oth,The∇ ·B= 0constraint in shock-capturing magnetohydrodynamics codes, Journal of Computational Physics161(2000), no. 2, 605–652

  61. [61]

    Lei Wei, Lingling Zhou, and Yinhua Xia,The jump filter in the discontinuous Galerkin method for hyperbolic conservation laws, Journal of Computational Physics520(2025), 113498

  62. [62]

    Winters and Gregor J

    Andrew R. Winters and Gregor J. Gassner,Affordable, entropy conserving and entropy stable flux functions for the ideal MHD equations, Journal of Computational Physics304 (2016), 72–108

  63. [63]

    4, 2124–2147

    Kailiang Wu,Positivity-preserving analysis of numerical schemes for ideal magnetohydro- dynamics, SIAM Journal on Numerical Analysis56(2018), no. 4, 2124–2147

  64. [64]

    5, B1302–B1329

    Kailiang Wu and Chi-Wang Shu,A provably positive discontinuous Galerkin method for multidimensional ideal magnetohydrodynamics, SIAM Journal on Scientific Computing40 (2018), no. 5, B1302–B1329

  65. [65]

    4, 995–1047

    ,Provably positive high-order schemes for ideal magnetohydrodynamics: analysis on general meshes, Numerische Mathematik142(2019), no. 4, 995–1047

  66. [66]

    Helen C. Yee and Bj¨ orn Sj¨ ogreen,Divergence free high order filter methods for the com- pressible MHD equations, Modeling, Simulation and Optimization of Complex Processes (Berlin, Heidelberg) (Hans Georg Bock, Hoang Xuan Phu, Ekaterina Kostina, and Rolf Rannacher, eds.), Springer, Berlin, Heidelberg, 2005, pp. 559–575

  67. [67]

    1, 115–164

    ,Efficient low dissipative high order schemes for multiscale MHD flows, II: Min- imization of∇ ·Bnumerical error, Journal of Scientific Computing29(2006), no. 1, 115–164

  68. [68]

    23, 8918–8934

    Xiangxiong Zhang and Chi-Wang Shu,On positivity-preserving high order discontinuous Galerkin schemes for compressible Euler equations on rectangular meshes, Journal of Com- putational Physics229(2010), no. 23, 8918–8934

  69. [69]

    1, 397–415

    Xinghui Zhong and Chi-Wang Shu,A simple weighted essentially nonoscillatory limiter for Runge–Kutta discontinuous Galerkin methods, Journal of Computational Physics232 (2013), no. 1, 397–415. 23