An Efficient Implementation of Flux Formulae in Multidimensional Relativistic Hydrodynamical Codes

Abstract

We derive and analyze a simplified formulation of the numerical viscosity terms appearing in the expression of the numerical fluxes associated to several High-Resolution Shock-Capturing schemes. After some algebraic pre-processing, we give explicit expressions for the numerical viscosity terms of two of the most widely used flux formulae, which implementation saves computational time in multidimensional simulations of relativistic flows. Additionally, such treatment explicitely cancells and factorizes a number of terms helping to amortiguate the growing of round-off errors. We have checked the performance of our formulation running a 3D relativistic hydrodynamical code to solve a standard test-bed problem and found that the improvement in efficiency is of high practical interest in numerical simulations of relativistic flows in Astrophysics.

10 pages, accepted for publication in Computer Physics Communications

Cited by in corpus (10)

• Three-dimensional relativistic simulations of rotating neutron-star collapse to a Kerr black hole
• THC: a new high-order finite-difference high-resolution shock-capturing code for special-relativistic hydrodynamics
• Spectral evolution of superluminal components in parsec-scale jets
• Assessment of a high-resolution central scheme for the solution of the relativistic hydrodynamics equations
• Excision methods for high resolution shock capturing schemes applied to general relativistic hydrodynamics
• Relativistic Magnetohydrodynamics: Renormalized eigenvectors and full wave decomposition Riemann solver
• On the measurements of numerical viscosity and resistivity in Eulerian MHD codes
• Observational Effects of Anomalous Boundary Layers in Relativistic Jets
• Efficient implementation of finite volume methods in Numerical Relativity
• An approximate Riemann solver approach in Physics-Informed Neural Networks for hyperbolic conservation laws