An adaptive mesh, GPU-accelerated, and error minimized special relativistic hydrodynamics code. (arXiv:2012.11130v2 [astro-ph.HE] UPDATED)
<a href="http://arxiv.org/find/astro-ph/1/au:+Tseng_P/0/1/0/all/0/1">Po-Hsun Tseng</a>, <a href="http://arxiv.org/find/astro-ph/1/au:+Schive_H/0/1/0/all/0/1">Hsi-Yu Schive</a>, <a href="http://arxiv.org/find/astro-ph/1/au:+Chiueh_T/0/1/0/all/0/1">Tzihong Chiueh</a>
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
$7times 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.
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
$7times 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.
http://arxiv.org/icons/sfx.gif
