Entropy-limited hydrodynamics: a novel approach to relativistic hydrodynamics
arXiv:1612.06251 · doi:10.1186/s40668-017-0022-0
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.
26 pages, 26 figures. Small revisions, matches the version to be published
References in corpus (16)
- Binary neutron-star mergers: a review of Einstein's richest laboratory
- THC: a new high-order finite-difference high-resolution shock-capturing code for special-relativistic hydrodynamics
- Tidal effects in binary neutron star coalescence
- General relativistic simulations of compact binary mergers as engines of short gamma-ray bursts
- WHAM: A WENO-based general relativistic numerical scheme I: Hydrodynamics
- Improved constrained scheme for the Einstein equations: An approach to the uniqueness issue
- Challenging the paradigm of singularity excision in gravitational collapse
- Covariant formulations of BSSN and the standard gauge
- SpECTRE: A Task-based Discontinuous Galerkin Code for Relativistic Astrophysics
- On the gravitational radiation from the collapse of neutron stars to rotating black holes
- Discontinuous Galerkin methods for general-relativistic hydrodynamics: formulation and application to spherically symmetric spacetimes
- Design of Provably Physical-Constraint-Preserving Methods for General Relativistic Hydrodynamics
- The trumpet solution from spherical gravitational collapse with puncture gauges
- Critical Phenomena in Neutron Stars I: Linearly Unstable Nonrotating Models
- A new three-dimensional general-relativistic hydrodynamics code
- High-Order Numerical-Relativity Simulations of Binary Neutron Stars
Cited by in corpus (13)
- The Science of the Einstein Telescope
- Gravitational waves and mass ejecta from binary neutron star mergers: Effect of large eccentricities
- Rotating neutron stars with non-barotropic thermal profile
- Entropy Symmetrization and High-Order Accurate Entropy Stable Numerical Schemes for Relativistic MHD Equations
- Full 3D Numerical Relativity Simulations of Neutron Star -- Boson Star Collisions with BAM
- The fate of twin stars on the unstable branch: implications for the formation of twin stars
- A simple, entropy-based dissipation trigger for SPH
- Entropy-limited higher-order central scheme for neutron star merger simulations
- A hybrid approach to long-term binary neutron-star simulations
- Disk formation in the collapse of supramassive neutron stars
- Extremely High-Order Convergence in Simulations of Relativistic Stars
- Minimum Principle on Specific Entropy and High-Order Accurate Invariant Region Preserving Numerical Methods for Relativistic Hydrodynamics
- Eccentricity reduction of binary neutron star initial data with the entropy based flux limiting scheme