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
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.
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
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.
Referee Report
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)
- [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.
- [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)
- Clarify the precise differences in implementation cost and accuracy between the new nodal scheme and the two referenced prior DG versions.
- Ensure consistent notation for the nodal basis functions, projection operator, and signal-speed estimates throughout the manuscript.
Simulated Author's Rebuttal
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
-
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
-
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
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
axioms (1)
- standard math Standard properties of discontinuous Galerkin discretizations and HLL-type numerical fluxes for hyperbolic systems
Reference graph
Works this paper leans on
-
[1]
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
work page 2018
-
[2]
Dinshaw S. Balsara,Self-adjusting, positivity preserving high order schemes for hydrody- namics and magnetohydrodynamics, Journal of Computational Physics231(2012), no. 22, 7504–7517
work page 2012
-
[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
work page 2021
-
[4]
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
work page 1999
-
[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
work page 1998
-
[6]
,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,...
work page 2006
-
[7]
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
work page 2020
-
[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
work page 2007
-
[9]
,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
work page 2010
-
[10]
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
work page 1980
-
[11]
Jesse Chan,On discretely entropy conservative and entropy stable discontinuous Galerkin methods, Journal of Computational Physics362(2018), 346–374
work page 2018
-
[12]
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
work page 2019
-
[13]
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
work page 2016
-
[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
work page 2020
-
[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
work page 2017
-
[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
work page 2020
-
[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
work page 2013
-
[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
work page 2000
- [19]
-
[20]
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
work page 1998
-
[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
work page 2018
-
[22]
Constantine Michael Dafermos,Quasilinear hyperbolic systems with involutions, Archive for Rational Mechanics and Analysis94(1986), no. 4, 373–389
work page 1986
-
[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
work page 2026
-
[24]
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
work page 1989
-
[25]
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
work page 1998
-
[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
work page 2002
-
[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
work page 2014
-
[28]
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
work page 2018
-
[29]
Alexandre Ern and Jean-Luc Guermond,Finite elements I, Texts in Applied Mathematics, vol. 72, Springer, Cham, 02 2021
work page 2021
-
[30]
Charles R. Evans and John F. Hawley,Simulation of magnetohydrodynamic flows: A con- strained transport method, The Astrophysical journal332(1988), 659–677
work page 1988
-
[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
work page 2012
-
[32]
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
work page 2018
-
[33]
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
work page 2005
- [34]
-
[35]
Sigal Gottlieb and Chi-Wang Shu,Total vriation diminishing Runge-Kutta schemes, Math- ematics of Computation67(1998), no. 221, 73–85
work page 1998
- [36]
-
[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
work page 2016
-
[38]
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
work page 2011
-
[39]
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
work page 2019
-
[40]
Rossmanith James Alexander,An unstaggered, high-resolution constrained transport method for magnetohydrodynamic flows, SIAM Journal on Scientific Computing28(2006), no. 5, 1766–1797
work page 2006
-
[41]
Antony Jameson,Eigenvalues, eigenvectors and symmetrization of the magneto- hydrodynamic (MHD) equations, AFOSR Grantees and Contractors Meeting, 06 2006
work page 2006
-
[42]
Pekka Janhunen,A positive conservative method for magnetohydrodynamics based on HLL and Roe methods, Journal of Computational Physics160(2000), no. 2, 649–661
work page 2000
-
[43]
Guang-Shan Jiang and Chi-Wang Shu,Efficient implementation of weighted ENO schemes, Journal of Computational Physics126(1996), no. 1, 202–228
work page 1996
-
[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
work page 2013
-
[45]
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
work page 2017
-
[46]
Fengyan Li and Chi-Wang Shu,Locally divergence-free discontinuous galerkin methods for MHD equations, Journal of Scientific Computing22(2005), no. 1, 413–442
work page 2005
-
[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
work page 2025
-
[48]
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
work page 2022
-
[49]
,An entropy stable essentially oscillation-free discontinuous Galerkin method for solving ideal magnetohydrodynamic equations, Journal of Computational Physics530 (2025), 113911
work page 2025
-
[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
work page 2018
-
[51]
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
work page 2025
-
[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
work page 1994
-
[53]
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
work page 1999
-
[54]
Jianxian Qiu and Chi-Wang Shu,Runge–Kutta discontinuous Galerkin method using WENO limiters, SIAM Journal on Scientific Computing26(2005), no. 3, 907–929
work page 2005
-
[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
work page 2024
-
[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
work page 2023
-
[57]
Chi-Wang Shu,Essentially non-oscillatory and weighted essentially non-oscillatory schemes, Acta Numerica29(2020), 701–762
work page 2020
-
[58]
Eitan Tadmor,Entropy stability theory for difference approximations of nonlinear conser- vation laws and related time-dependent problems, Acta Numerica12(2003), 451–512
work page 2003
- [59]
-
[60]
G´ abor T´ oth,The∇ ·B= 0constraint in shock-capturing magnetohydrodynamics codes, Journal of Computational Physics161(2000), no. 2, 605–652
work page 2000
-
[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
work page 2025
-
[62]
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
work page 2016
-
[63]
Kailiang Wu,Positivity-preserving analysis of numerical schemes for ideal magnetohydro- dynamics, SIAM Journal on Numerical Analysis56(2018), no. 4, 2124–2147
work page 2018
-
[64]
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
work page 2018
-
[65]
,Provably positive high-order schemes for ideal magnetohydrodynamics: analysis on general meshes, Numerische Mathematik142(2019), no. 4, 995–1047
work page 2019
-
[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
work page 2005
-
[67]
,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
work page 2006
-
[68]
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
work page 2010
-
[69]
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
work page 2013
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.