Abstract:Simulating general relativistic hydrodynamics (GRHD) presents challenges such as handling curved spacetime, achieving high-order shock-capturing accuracy, and preserving key physical constraints (positive density, pressure, and subluminal velocity) under nonlinear coupling. This paper introduces high-order, physical-constraint-preserving, oscillation-eliminating discontinuous Galerkin (PCP-OEDG) schemes with Harten-Lax-van Leer flux for GRHD. To suppress spurious oscillations near discontinuities, we incorporate a computationally efficient oscillation-eliminating (OE) procedure based on a linear damping equation, maintaining accuracy and avoiding complex characteristic decomposition. To enhance stability and robustness, we construct PCP schemes using the W-form of GRHD equations with Cholesky decomposition of the spatial metric, addressing the non-equivalence of admissible state sets in curved spacetime. We rigorously prove the PCP property of cell averages via technical estimates and the Geometric Quasi-Linearization (GQL) approach, which transforms nonlinear constraints into linear forms. Additionally, we present provably convergent PCP iterative algorithms for robust recovery of primitive variables, ensuring physical constraints are satisfied throughout. The PCP-OEDG method is validated through extensive tests, demonstrating its robustness, accuracy, and capability to handle extreme GRHD scenarios involving strong shocks, high Lorentz factors, and intense gravitational fields.

Constraint Propagation of $C^2$-adjusted Formulation - Another Recipe for Robust ADM Evolution System

Constructing of constraint preserving scheme for Einstein equations

ContEvol formalism: possibly a new twist on computational physics

Re-formulating the Einstein equations for stable numerical simulations: Formulation Problem in Numerical Relativity

The formally second-order BDF ADI difference/compact difference scheme for the nonlocal evolution problem in three-dimensional space

Einstein boundary conditions in relation to constraint propagation for the initial-boundary value problem of the Einstein equations

Adjoint-Based Design Optimization of Stability Constrained Systems

A Formally Second-Order BDF Compact ADI Difference Scheme for the Two-Dimensional Fractional Evolution Equation

A Two-Level Linearized Compact ADI Scheme for Two-Dimensional Nonlinear Reaction-Diffusion Equations

Stable $C^1$-conforming finite element methods for a class of nonlinear fourth-order evolution equations

Constraint Energy Minimizing Generalized Multiscale Finite Element Method for Convection Diffusion Equations with Inhomogeneous Boundary Conditions

Constraint-preserved numerical schemes with decoupling structure for the Ericksen–Leslie model with variable density

Stability analysis and improvement of the covariant BSSN formulation against the FLRW spacetime background

Adaptive mesh refinement for coupled elliptic-hyperbolic systems

A second-order exponential integration constraint energy minimizing generalized multiscale method for parabolic problems

A Kinetic Scheme Based On Positivity Preservation For Multi-component Euler Equations

A Higher-Order Compact Adi Method with Monotone Iterative Procedure for Systems of Reaction-Diffusion Equations

Constraint preserving boundary conditions in Bondi-Sachs gauge: a numerical study of stability of pure AdS spacetime

Robust Discontinuous Galerkin Methods Maintaining Physical Constraints for General Relativistic Hydrodynamics

Semi-Implicit-type order-adaptive CAT2 schemes for systems of balance laws with relaxed source term