Provably Positive Discontinuous Galerkin Methods for Multidimensional Ideal Magnetohydrodynamics

Abstract

The density and pressure are positive physical quantities in magnetohydrodynamics (MHD). Design of provably positivity-preserving (PP) numerical schemes for ideal compressible MHD is highly desirable, but remains a challenge especially in the multidimensional cases. In this paper, we first develop uniformly high-order discontinuous Galerkin (DG) schemes which provably preserve the positivity of density and pressure for multidimensional ideal MHD. The schemes are constructed by using the locally divergence-free DG schemes for the symmetrizable ideal MHD equations as the base schemes, a PP limiter to enforce the positivity of the DG solutions, and the strong stability preserving methods for time discretization. The significant innovation is that we discover and rigorously prove the PP property of the proposed DG schemes by using a novel equivalent form of the admissible state set and very technical estimates. Several two-dimensional numerical examples further confirm the PP property, and demonstrate the accuracy, effectiveness and robustness of the proposed PP methods.

Submitted for publication on Jan. 31, 2018

References in corpus (4)

• A robust numerical scheme for highly compressible magnetohydrodynamics: Nonlinear stability, implementation and tests
• Positivity-Preserving Analysis of Numerical Schemes for Ideal Magnetohydrodynamics
• Design of Provably Physical-Constraint-Preserving Methods for General Relativistic Hydrodynamics
• Runge-Kutta discontinuous Galerkin methods for the special relativistic magnetohydrodynamics

Cited by in corpus (12)

• Provably Positive High-Order Schemes for Ideal Magnetohydrodynamics: Analysis on General Meshes
• Positivity-Preserving Entropy-Based Adaptive Filtering for Discontinuous Spectral Element Methods
• Entropy Symmetrization and High-Order Accurate Entropy Stable Numerical Schemes for Relativistic MHD Equations
• Physical-constraints-preserving Lagrangian finite volume schemes for one- and two-dimensional special relativistic hydrodynamics
• A Physical-Constraint-Preserving Finite Volume WENO Method for Special Relativistic Hydrodynamics on Unstructured Meshes
• Is the Classic Convex Decomposition Optimal for Bound-Preserving Schemes in Multiple Dimensions?
• On physical-constraints-preserving schemes for special relativistic magnetohydrodynamics with a general equation of state
• High order asymptotic preserving finite difference WENO schemes with constrained transport for MHD equations in all sonic Mach numbers
• Positivity-preserving entropy filtering for the ideal magnetohydrodynamics equations
• Monolithic parabolic regularization of the MHD equations and entropy principles
• Physical Constraint Preserving Higher Order Finite Volume Schemes for Divergence-Free Astrophysical MHD and RMHD
• An Alternative Finite Difference WENO-like Scheme with Physical Constraint Preservation for Divergence-Preserving Hyperbolic Systems