Physical-constraints-preserving central discontinuous Galerkin methods for special relativistic hydrodynamics with a general equation of state

Abstract

The ideal gas equation of state (EOS) with a constant adiabatic index is a poor approximation for most relativistic astrophysical flows, although it is commonly used in relativistic hydrodynamics. The paper develops high-order accurate physical-constraints-preserving (PCP) central discontinuous Galerkin (DG) methods for the one- and two-dimensional special relativistic hydrodynamic (RHD) equations with a general EOS. It is built on the theoretical analysis of the admissible states for the RHD and the PCP limiting procedure enforcing the admissibility of central DG solutions. The convexity, scaling and orthogonal invariance, and Lax-Friedrichs splitting property of the admissible state set are first proved with the aid of its equivalent form, and then the high-order central DG methods with the PCP limiting procedure and strong stability preserving time discretization are proved to preserve the positivity of the density, pressure, and specific internal energy, and the bound of the fluid velocity, maintain the high-order accuracy, and be $L^1$-stable. The accuracy, robustness, and effectiveness of the proposed methods are demonstrated by several 1D and 2D numerical examples involving large Lorentz factor, strong discontinuities, or low density or pressure etc.

32 pages, 13 figures

References in corpus (1)

• Equation of State in Numerical Relativistic Hydrodynamics

Cited by in corpus (14)

• Provably Positive High-Order Schemes for Ideal Magnetohydrodynamics: Analysis on General Meshes
• Positivity-Preserving Analysis of Numerical Schemes for Ideal Magnetohydrodynamics
• Design of Provably Physical-Constraint-Preserving Methods for General Relativistic Hydrodynamics
• Entropy stable adaptive moving mesh schemes for 2D and 3D special relativistic hydrodynamics
• Runge-Kutta discontinuous Galerkin methods for the special relativistic magnetohydrodynamics
• 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
• On physical-constraints-preserving schemes for special relativistic magnetohydrodynamics with a general equation of state
• On the relativistic compressibility conditions
• Second-order accurate genuine BGK schemes for the ultra-relativistic flow simulations
• High-order discontinuous Galerkin schemes with subcell finite volume limiter and adaptive mesh refinement for a monolithic first-order BSSNOK formulation of the Einstein-Euler equations
• Constraints Preserving Lax-Wendroff Flux Reconstruction for Relativistic Hydrodynamics with General Equations of State