An adaptive mesh, GPU-accelerated, and error minimized special relativistic hydrodynamics code

Abstract

We present a new special relativistic hydrodynamics (SRHD) code capable of handling coexisting ultra-relativistically hot and non-relativistically cold gases. We achieve this by designing a new algorithm for conversion between primitive and conserved variables in the SRHD solver, which incorporates a realistic ideal-gas equation of state covering both the relativistic and non-relativistic regimes. The code can handle problems involving a Lorentz factor as high as $10^6$ and optimally avoid the catastrophic cancellation. In addition, we have integrated this new SRHD solver into the code GAMER (https://github.com/gamer-project/gamer) to support adaptive mesh refinement and hybrid OpenMP/MPI/GPU parallelization. It achieves a peak performance of $7\times 10^7$ cell updates per second on a single Tesla P100 GPU and scales well to 2048 GPUs. We apply this code to two interesting astrophysical applications: (a) an asymmetric explosion source on the relativistic blast wave and (b) the flow acceleration and limb-brightening of relativistic jets.

References in corpus (25)

• A Multi-Code Analysis Toolkit for Astrophysical Simulation Data
• Relativistic Jets in Active Galactic Nuclei
• Superluminal motion of a relativistic jet in the neutron star merger GW170817
• The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers
• Compact radio emission indicates a structured jet was produced by a binary neutron star merger
• A mildly relativistic wide-angle outflow in the neutron star merger GW170817
• Primitive Variable Solvers for Conservative General Relativistic Magnetohydrodynamics
• The Structure of the M87 Jet: A Transition from Parabolic to Conical Streamlines
• An HLLC Riemann Solver for Relativistic Flows: I. Hydrodynamics
• Parsec Scale Properties of Markarian 501
• The Piecewise Parabolic Method for Multidimensional Relativistic Fluid Dynamics
• Equation of State in Relativistic Magnetohydrodynamics: variable versus constant adiabatic index
• The limb-brightened jet of M87 down to 7 Schwarzschild radii scale
• An HLLC Solver for Relativistic Flows -- II. Magnetohydrodynamics
• A wide and collimated radio jet in 3C 84 on the scale of a few hundred gravitational radii
• GAMER: a GPU-Accelerated Adaptive Mesh Refinement Code for Astrophysics
• Limb-Brightened Jet of 3C 84 Revealed by the 43-GHz Very-Long-Baseline-Array Observation
• The Optical Afterglow of GW170817: An Off-axis Structured Jet and Deep Constraints on a Globular Cluster Origin
• The stratified two-sided jet of Cygnus A. Acceleration and collimation
• Equation of State in Numerical Relativistic Hydrodynamics
• GAMER-2: a GPU-accelerated adaptive mesh refinement code -- accuracy, performance, and scalability
• An Improved Exact Riemann Solver for Relativistic Hydrodynamics
• Reconfinement and Loss of Stability in Jets from Active Galactic Nuclei
• The Spatially Uniform Spectrum of the Fermi Bubbles: the Leptonic AGN Jet Scenario
• CAFE: A New Relativistic MHD Code