Relativistic Hydrodynamic Flows Using Spatial and Temporal Adaptive Structured Mesh Refinement

Peng Wang,Tom Abel,Weiqun Zhang

DOI: https://doi.org/10.1086/529434

2008-01-01

Abstract:Astrophysical relativistic flow problems require high resolution three-dimensional numerical simulations. In this paper, we describe a new parallel three-dimensional code for simulations of special relativistic hydrodynamics (SRHD) using both spatially and temporally structured adaptive mesh refinement (AMR). We used the method of lines to discretize the SRHD equations spatially and a total variation diminishing (TVD) Runge-Kutta scheme for time integration. For spatial reconstruction, we have implemented piecewise linear method (PLM), piecewise parabolic method (PPM), third order convex essentially non-oscillatory (CENO) and third and fifth order weighted essentially non-oscillatory (WENO) schemes. Flux is computed using either direct flux reconstruction or approximate Riemann solvers including HLL, modified Marquina flux, local Lax-Friedrichs flux formulas and HLLC. The AMR part of the code is built on top of the cosmological Eulerian AMR code {\sl enzo}. We discuss the coupling of the AMR framework with the relativistic solvers. Via various test problems, we emphasize the importance of resolution studies in relativistic flow simulations because extremely high resolution is required especially when shear flows are present in the problem. We also present the results of two 3d simulations of astrophysical jets: AGN jets and GRB jets. Resolution study of those two cases further highlights the need of high resolutions to calculate accurately relativistic flow problems.

Astrophysics

What problem does this paper attempt to address?

This paper aims to solve the problem of high - resolution three - dimensional numerical simulation in relativistic hydrodynamics. Specifically, the paper describes a new parallel three - dimensional code for simulating special relativistic hydrodynamics (SRHD) using the Adaptive Mesh Refinement (AMR) technique in space and time. The focus of the paper is as follows:

1. **High - resolution Requirement**: Relativistic hydrodynamics problems require extremely high resolution, especially in the presence of shear flow. The paper emphasizes the importance of this through various test problems.

2. **Methods and Techniques**:
- Spatial discretization of SRHD equations is carried out using the Method of Lines, and time integration is performed using the Total Variation Diminishing (TVD) Runge - Kutta scheme.
- Multiple spatial reconstruction methods are implemented, including the Piecewise Linear Method (PLM), the Piecewise Parabolic Method (PPM), the third - order Convex Essentially Non - oscillatory (CENO) and the third - and fifth - order Weighted Essentially Non - oscillatory (WENO) schemes.
- Flux calculations use either direct flux reconstruction or approximate Riemann solvers, such as HLL, the modified Marquina flux, the local Lax - Friedrichs flux formula and HLLC.

3. **Adaptive Mesh Refinement (AMR)**: The AMR part is built based on the cosmological Euler AMR code Enzo. The paper discusses how to couple the AMR framework with the relativistic solver.

4. **Application Examples**: The paper presents the simulation results of two three - dimensional astrophysical jets, namely the Active Galactic Nucleus (AGN) jet and the Gamma - Ray Burst (GRB) jet. These simulations further highlight the importance of high resolution in accurately calculating relativistic hydrodynamics problems.

In conclusion, the main goal of this paper is to develop an efficient numerical tool that can handle complex relativistic hydrodynamics problems, especially those astrophysical phenomena that require extremely high resolution.