Admissible state and physical constraints preserving schemes for relativistic magnetohydrodynamic equations

Abstract

This paper first studies the admissible state set $\mathcal{G}$ of relativistic magnetohydrodynamics (RMHD). It paves a way for developing physical-constraints-preserving (PCP) schemes for RMHD equations with the solutions in $\mathcal{G}$. To overcome the difficulties arising from the extremely strong nonlinearities and no explicit formulas of the primitive variables and flux vectors with respect to the conservative vector, two equivalent forms of $\mathcal{G}$ with explicit constraints on the conservative vector are skillfully discovered. The first is derived by analyzing roots of several polynomials and transferring successively them, and further used to prove the convexity of $\mathcal{G}$ with the aid of semi-positive definiteness of the second fundamental form of a hypersurface. While the second is derived based on the convexity and then used to show the orthogonal invariance of $\mathcal{G}$. The Lax-Friedrichs (LxF) splitting property does not hold generally for nonzero magnetic field, but by a constructive inequality and pivotal techniques, we discover the generalized LxF splitting properties, combining the convex combination of some LxF splitting terms with a discrete divergence-free condition of magnetic field. Based on above analyses, several one- and two-dimensional PCP schemes are then studied. In 1D case, a first-order accurate LxF type scheme is first proved to be PCP under the CFL condition, and then high-order PCP schemes are proposed via a PCP limiter. In 2D case, the discrete divergence-free condition and PCP property are analyzed for a first-order LxF scheme, and two sufficient conditions are derived for high-order accurate PCP schemes. Our analysis reveals in theory for the first time that the discrete divergence-free condition is closely connected with the PCP property. Several numerical examples demonstrate theoretical findings and performance of numerical schemes.

50 pages, 11 figure

References in corpus (7)

• ECHO: an Eulerian Conservative High Order scheme for general relativistic magnetohydrodynamics and magnetodynamics
• An HLLC Solver for Relativistic Flows -- II. Magnetohydrodynamics
• Relativistic MHD with Adaptive Mesh Refinement
• A multidimensional grid-adaptive relativistic magnetofluid code
• Design of Provably Physical-Constraint-Preserving Methods for General Relativistic Hydrodynamics
• Physical-constraints-preserving central discontinuous Galerkin methods for special relativistic hydrodynamics with a general equation of state
• Runge-Kutta discontinuous Galerkin methods for the special relativistic magnetohydrodynamics

Cited by in corpus (13)

• Provably Positive High-Order Schemes for Ideal Magnetohydrodynamics: Analysis on General Meshes
• Positivity-Preserving Analysis of Numerical Schemes for Ideal Magnetohydrodynamics
• Provably Positive Discontinuous Galerkin Methods for Multidimensional Ideal Magnetohydrodynamics
• Entropy stable adaptive moving mesh schemes for 2D and 3D special relativistic hydrodynamics
• Entropy Symmetrization and High-Order Accurate Entropy Stable Numerical Schemes for Relativistic MHD Equations
• High-order accurate entropy stable nodal discontinuous Galerkin schemes for the ideal special relativistic magnetohydrodynamics
• Physical-constraints-preserving Lagrangian finite volume schemes for one- and two-dimensional special relativistic hydrodynamics
• High-order accurate entropy stable adaptive moving mesh finite difference schemes for special relativistic (magneto)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
• Recent developments in mathematical aspects of relativistic fluids
• Second-order accurate BGK schemes for the special relativistic hydrodynamics with the Synge equation of state