Design of Provably Physical-Constraint-Preserving Methods for General Relativistic Hydrodynamics

Abstract

The paper develops high-order physical-constraint-preserving (PCP) methods for general relativistic hydrodynamic (GRHD) equations, equipped with a general equation of state. Here the physical constraints, describing the admissible states of GRHD, are referred to the subluminal constraint on the fluid velocity and the positivity of the density, pressure and specific internal energy. Preserving these constraints is very important for robust computations, otherwise violating one of them will lead to the ill-posed problem and numerical instability. To overcome the difficulties arising from the inherent strong nonlinearity contained in the constraints, we derive an equivalent definition of the admissible states. Using this definition, we prove the convexity, scaling invariance and Lax-Friedrichs (LxF) splitting property of the admissible state set $\mathcal{G}$, and discover the dependence of $\mathcal{G}$ on the spacetime metric. Unfortunately, such dependence yields the non-equivalence of $\mathcal{G}$ at different points in curved spacetime, and invalidates the convexity of $\mathcal{G}$ in analyzing PCP schemes. This obstacle is effectively overcame by introducing a new formulation of the GRHD equations. Based on this formulation and the above theories, a first-order LxF scheme is designed on general unstructured mesh and rigorously proved to be PCP under a CFL condition. With two types of PCP limiting procedures, we design high-order, {\em provably} (not probably) PCP methods under discretization on the proposed new formulation. These high-order methods include the PCP finite difference, finite volume and discontinuous Galerkin methods.

17 pages

References in corpus (3)

• Equation of State in Numerical Relativistic Hydrodynamics
• Physical-constraints-preserving central discontinuous Galerkin methods for special relativistic hydrodynamics with a general equation of state
• Bound-Preserving Discontinuous Galerkin Methods for Conservative Phase Space Advection in Curvilinear Coordinates

Cited by in corpus (21)

• 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
• Numerical relativity in spherical coordinates: A new dynamical spacetime and general relativistic MHD evolution framework for the Einstein Toolkit
• Uniformly High-Order Structure-Preserving Discontinuous Galerkin Methods for Euler Equations with Gravitation: Positivity and Well-Balancedness
• Physical-constraints-preserving Lagrangian finite volume schemes for one- and two-dimensional special relativistic hydrodynamics
• High-order accurate entropy stable nodal discontinuous Galerkin schemes for the ideal special relativistic magnetohydrodynamics
• High-order accurate entropy stable adaptive moving mesh finite difference schemes for special relativistic (magneto)hydrodynamics
• Entropy-limited hydrodynamics: a novel approach to 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
• On the relativistic compressibility conditions
• A new first-order formulation of the Einstein equations exploiting analogies with electrodynamics
• Two-stage fourth-order accurate time discretizations for 1D and 2D special relativistic hydrodynamics
• Second-order accurate BGK schemes for the special relativistic hydrodynamics with the Synge equation of state
• Geometric Quasilinearization Framework for Analysis and Design of Bound-Preserving Schemes
• Provably Physical-Constraint-Preserving Discontinuous Galerkin Methods for Multidimensional Relativistic MHD Equations
• A physical-constraints-preserving genuinely multidimensional HLL scheme for the special relativistic hydrodynamics
• Mixmaster Fluids Near the Big Bang
• Minimum Principle on Specific Entropy and High-Order Accurate Invariant Region Preserving Numerical Methods for Relativistic Hydrodynamics