PHYSICAL-CONSTRAINT-PRESERVING CENTRAL DISCONTINUOUS GALERKIN METHODS FOR SPECIAL RELATIVISTIC HYDRODYNAMICS WITH A GENERAL EQUATION OF STATE
- School of Mathematical Sciences, Peking University, Beijing 100871 (China)
- HEDPS, CAPT and LMAM, School of Mathematical Sciences, Peking University, Beijing 100871 (China)
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 (RHD). This paper develops high-order accurate, physical-constraints-preserving (PCP), central, discontinuous Galerkin (DG) methods for the one- and two-dimensional special RHD equations with a general EOS. It is built on our theoretical analysis of the admissible states for RHD and the PCP limiting procedure that enforce the admissibility of central DG solutions. The convexity, scaling invariance, orthogonal invariance, and Lax–Friedrichs splitting property of the admissible state set are first proved with the aid of its equivalent form. 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, specific internal energy, and the bound of the fluid velocity, maintain high-order accuracy, and be L {sup 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/pressure, etc.
- OSTI ID:
- 22661391
- Journal Information:
- Astrophysical Journal, Supplement Series, Vol. 228, Issue 1; Other Information: Country of input: International Atomic Energy Agency (IAEA); ISSN 0067-0049
- Country of Publication:
- United States
- Language:
- English
Similar Records
Sparse invariant domain preserving discontinuous Galerkin methods with subcell convex limiting
CosmosDG: An hp -adaptive Discontinuous Galerkin Code for Hyper-resolved Relativistic MHD