Abstract:The relativistic hydrodynamics (RHD) equations have three crucial intrinsic physical constraints on the primitive variables: positivity of pressure and density, and subluminal fluid velocity. However, numerical simulations can violate these constraints, leading to nonphysical results or even simulation failure. Designing genuinely physical-constraint-preserving (PCP) schemes is very difficult, as the primitive variables cannot be explicitly reformulated using conservative variables due to relativistic effects. In this paper, we propose three efficient Newton--Raphson (NR) methods for robustly recovering primitive variables from conservative variables. Importantly, we rigorously prove that these NR methods are always convergent and PCP, meaning they preserve the physical constraints throughout the NR iterations. The discovery of these robust NR methods and their PCP convergence analyses are highly nontrivial and technical. As an application, we apply the proposed NR methods to design PCP finite volume Hermite weighted essentially non-oscillatory (HWENO) schemes for solving the RHD equations. Our PCP HWENO schemes incorporate high-order HWENO reconstruction, a PCP limiter, and strong-stability-preserving time discretization. We rigorously prove the PCP property of the fully discrete schemes using convex decomposition techniques. Moreover, we suggest the characteristic decomposition with rescaled eigenvectors and scale-invariant nonlinear weights to enhance the performance of the HWENO schemes in simulating large-scale RHD problems. Several demanding numerical tests are conducted to demonstrate the robustness, accuracy, and high resolution of the proposed PCP HWENO schemes and to validate the efficiency of our NR methods.

A Direct Eulerian GRP Scheme for Relativistic Hydrodynamics: One-dimensional Case

A Direct Eulerian Grp Scheme For Relativistic Hydrodynamics: Two-Dimensional Case

A Direct Eulerian GRP Scheme for Spherically Symmetric General Relativistic Hydrodynamics

A direct Eulerian GRP scheme for radiation hydrodynamical equations in diffusion limit

A Third-Order Accurate Direct Eulerian GRP Scheme for the Euler Equations in Gas Dynamics.

Second-order Direct Eulerian GRP Schemes for Radiation Hydrodynamical Equations

A second-order direct Eulerian GRP scheme for ten-moment Gaussian closure equations with source terms

An adaptive GRP scheme for compressible fluid flows

A Three-Dimensional Finite Volume Method for Conservation Laws in Conjunction with Modified Solution for Nonlinear Coupled Constitutive Relations

The generalized Riemann problem scheme for a laminar two-phase flow model with two-velocities

Genuinely multidimensional physical-constraints-preserving finite volume schemes for the special relativistic hydrodynamics

Efficient solver of relativistic hydrodynamics with implicit Runge-Kutta method

Runge-Kutta discontinuous Galerkin methods with WENO limiter for the special relativistic hydrodynamics

High-order accurate physical-constraints-preserving finite difference WENO schemes for special relativistic hydrodynamics

Accuracy of the Adaptive GRP Scheme and the Simulation of 2-D Riemann Problems for Compressible Euler Equations

Generalized Riemann Problem Method for the Kapila Model of Compressible Multiphase Flows I: Temporal-Spatial Coupling Finite Volume Scheme

Two-stage Fourth-Order Accurate Time Discretizations for 1D and 2D Special Relativistic Hydrodynamics

The Generalized Riemann Problem Method For The Shallow Water Equations With Bottom Topography

Robust Discontinuous Galerkin Methods Maintaining Physical Constraints for General Relativistic Hydrodynamics

Provably convergent Newton-Raphson methods for recovering primitive variables with applications to physical-constraint-preserving Hermite WENO schemes for relativistic hydrodynamics