Abstract:We present entropy-limited hydrodynamics (ELH): a new approach for the computation of numerical fluxes arising in the discretization of hyperbolic equations in conservation form. ELH is based on the hybridisation of an unfiltered high-order scheme with the first-order Lax-Friedrichs method. The activation of the low-order part of the scheme is driven by a measure of the locally generated entropy inspired by the artificial-viscosity method proposed by Guermond et al. Here, we present ELH in the context of high-order finite-differencing methods and of the equations of general-relativistic hydrodynamics. We study the performance of ELH in a series of classical astrophysical tests in general relativity involving isolated, rotating and nonrotating neutron stars, and including a case of gravitational collapse to black hole. We present a detailed comparison of ELH with the fifth-order monotonicity preserving method MP5, one of the most common high-order schemes currently employed in numerical-relativity simulations. We find that ELH achieves comparable and, in many of the cases studied here, better accuracy than more traditional methods at a fraction of the computational cost (up to 50 % speedup). Given its accuracy and its simplicity of implementation, ELH is a promising framework for the development of new special- and general-relativistic hydrodynamics codes well adapted for massively parallel supercomputers.

What problem does this paper attempt to address?

### Problems the Paper Aims to Solve This paper aims to propose a new numerical method—Entropy-Limited Hydrodynamics (ELH)—to address the numerical simulation issues in relativistic hydrodynamics. Specifically, the ELH method calculates numerical fluxes by combining a high-order unfiltered scheme with a low-order Lax-Friedrichs method, and activates the low-order part through locally generated entropy. This method aims to ensure high accuracy while improving numerical stability and computational efficiency. ### Main Problem Background 1. **Choice of Numerical Methods**: - In the numerical simulation of relativistic hydrodynamics, choosing an appropriate numerical method is crucial. Existing High-Resolution Shock-Capturing (HRSC) techniques, although effective, have some drawbacks such as complexity, dependence on optimized parameters, and potential load imbalance in parallel implementations. 2. **Numerical Stability and Accuracy**: - High-order methods can provide higher accuracy but tend to produce oscillations when dealing with discontinuities like shocks. Low-order methods, while stable, have lower accuracy. Therefore, ensuring numerical stability while maintaining high accuracy is a challenge. 3. **Computational Efficiency**: - Numerical simulations of relativistic hydrodynamics typically require substantial computational resources. Thus, developing a method that ensures accuracy while improving computational efficiency is very important. ### Solution 1. **Entropy-Limited Hydrodynamics (ELH)**: - The ELH method introduces a "shock detector" based on local entropy generation to determine which grid points need low-order contributions. This hybrid approach can restore high-order accuracy in smooth regions while maintaining numerical stability in shock regions. 2. **Numerical Implementation**: - The authors implemented the ELH method within a finite difference framework and applied it to the equations of general relativistic hydrodynamics. The performance of the ELH method was validated through a series of classical astrophysical tests. 3. **Performance Comparison**: - The authors compared the ELH method with the commonly used fifth-order Monotonicity-Preserving method (MP5) and found that ELH not only has comparable or even better accuracy in many cases but also lower computational costs (up to 50% speedup). ### Conclusion As a new numerical method, the ELH method can improve numerical stability and computational efficiency while maintaining high accuracy. This makes ELH a promising framework for developing new special and general relativistic hydrodynamics codes suitable for large-scale parallel supercomputers.

Entropy-limited hydrodynamics: a novel approach to relativistic hydrodynamics

Entropy Stable Adaptive Moving Mesh Schemes for 2D and 3D Special Relativistic Hydrodynamics

Entropy based flux limiting scheme for conservation laws

High-order Accurate Entropy Stable Finite Difference Schemes for One- and Two-Dimensional Special Relativistic Hydrodynamics

High-order Accurate Entropy Stable Schemes for Relativistic Hydrodynamics with General Synge-type Equation of State

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.

General-relativistic hydrodynamics of non-perfect fluids: 3+1 conservative formulation and application to viscous black-hole accretion

Entropy Stable Flux Correction for Hydrostatic Reconstruction Scheme for Shallow Water Flows

Entropy Symmetrization and High-Order Accurate Entropy Stable Numerical Schemes for Relativistic MHD Equations

A molecular dynamics approach to dissipative relativistic hydrodynamics: propagation of fluctuations

High order conservative Lagrangian schemes for two-dimensional radiation hydrodynamics equations in the equilibrium-diffusion limit

Bound Preserving Lax-Wendroff Flux Reconstruction Method for Special Relativistic Hydrodynamics

A dissipative extension to ideal hydrodynamics

New theories of relativistic hydrodynamics in the LHC era

Numerical hydrodynamics in general relativity

Universality classes of relativistic fluid dynamics: Applications

General relativistic hydrodynamics in curvilinear coordinates

Hyperhydrodynamics: Relativistic Viscous Fluids Emerging from Hypermomentum