A new code for the numerical simulation of relativistic flows on supercomputers by means of a low-dissipation scheme

Abstract

A new code to simulate special relativistic hydrodynamic flows on supercomputer architectures with distributed memory is described. The code is based on a combination of Godunov’s method and a piecewise parabolic method with a local stencil. This approach has good conservation properties, correctly reproduces shock waves, and ensures high accuracy on smooth solutions and low dissipation on discontinuities. Only a local computation stencil is needed for the piecewise parabolic reconstruction of the solution. The code scalability is 94% on a cluster, Intel Xeon X5670 NKS-30T, with 768 cores. The results of code verification using a relativistic jet problem and computational experiments on the evolution of a galactic jet are presented. Program summaryProgram title: RHydroBox3DCPC Library link to program files:http://dx.doi.org/10.17632/xskchgynfy.1Licensing provisions: GNU General Public License 3Programming language: C++Nature of problem: The special relativistic hydrodynamic flows.Solution method: The combination of Godunov’s method and a piecewise parabolic method on a local stencil.Additional comments including restrictions and unusual features: External library – MPI.