Numerical Relativistic Hydrodynamics Based on the Total Variation Diminishing Scheme

Eunwoo Choi,Dongsu Ryu

DOI: https://doi.org/10.1016/j.newast.2005.06.010

2005-07-05

Abstract:This paper describes a multidimensional hydrodynamic code which can be used for studies of relativistic astrophysical flows. The code solves the special relativistic hydrodynamic equations as a hyperbolic system of conservation laws based on the total variation diminishing (TVD) scheme. It uses a new set of conserved quantities and employs an analytic formula for transformation from the conserved quantities in the reference frame to the physical quantities in the local rest frame. Several standard tests, including relativistic shock tubes, a relativistic wall shock, and a relativistic blast wave, are presented to demonstrate that the code captures discontinuities correctly and sharply in ultrarelativistic regimes. The robustness and flexibility of the code are demonstrated through test simulations of the relativistic Hawley-Zabusky shock and a relativistic extragalactic jet.

Astrophysics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the relativistic hydrodynamics problems involved in high - energy astrophysics. Specifically, the author has developed a multi - dimensional special relativistic hydrodynamics code, aiming to study relativistic astrophysical flows with high Lorentz factors. Such problems are common in many high - energy astrophysical phenomena, for example:

- **Superluminal motion**: The apparently superluminal motion observed in microquasars in the Milky Way and extragalactic radio sources associated with active galactic nuclei requires an intrinsic beam velocity exceeding 0.9c to be explained.
- **Extreme jet velocities**: In some powerful extragalactic radio sources, the true beam velocities of matter ejected from the galactic nuclei exceed 0.98c.
- **Supernova explosions**: In the early stages of supernova explosions, the rapid expansion process also requires a relativistic description.
- **Gamma - ray bursts**: In the fireball model of gamma - ray bursts, the internal energy of the gas is converted into bulk kinetic energy, resulting in a relativistic outflow with a high Lorentz factor (> 100).

Since these relativistic hydrodynamics problems are highly nonlinear and inherently complex, they can usually only be studied through numerical simulations. Therefore, the main goal of this paper is to develop a multi - dimensional special relativistic hydrodynamics code based on the Total Variation Diminishing (TVD) scheme to capture and accurately describe fluid phenomena under these extreme conditions. Through a series of standard tests, such as relativistic shock tubes, relativistic wall shocks, relativistic detonation waves, etc., the validity and accuracy of the code have been verified. In addition, by simulating relativistic HLL (Harten - Lax - van Leer) shocks and relativistic extragalactic jets, the robustness and flexibility of the code in dealing with actual astrophysical problems have been demonstrated.