On physical-constraints-preserving schemes for special relativistic magnetohydrodynamics with a general equation of state
This work provides theoretical guarantees for numerical schemes in relativistic MHD with general equations of state, addressing a known bottleneck in extending previous methods.
The paper extends physical-constraints-preserving (PCP) schemes for special relativistic magnetohydrodynamics to general equations of state and general meshes, proving that the PCP property holds under a milder discrete divergence-free condition.
The paper studies the physical-constraints-preserving (PCP) schemes for multi-dimensional special relativistic magnetohydrodynamics with a general equation of state (EOS) on more general meshes. It is an extension of the work [Math. Models Methods Appl. Sci., 27:1871-1928, 2017] which focuses on the ideal EOS and uniform Cartesian meshes. The general EOS without a special expression poses some additional difficulties in discussing the mathematical properties of admissible state set with the physical constraints on the fluid velocity, density and pressure. Rigorous analyses are provided for the PCP property of finite volume or discontinuous Galerkin schemes with the Lax-Friedrichs (LxF) type flux on a general mesh with non-self-intersecting polytopes. Those are built on a more general form of generalized LxF splitting property and a different convex decomposition technique. It is shown in theory that the PCP property is closely connected with a discrete divergence-free condition, which is proposed on the general mesh and milder than that in [Math. Models Methods Appl. Sci., 27:1871-1928, 2017].