An Adaptive Grid, Implicit Code for Spherically Symmetric, General Relativistic Hydrodynamics in Comoving Coordinates

Matthias Liebendoerfer,Stephan Rosswog,Friedrich-Karl Thielemann

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

2002-04-23

Abstract:We describe an implicit general relativistic hydrodynamics code. The evolution equations are formulated in comoving coordinates. A conservative finite differencing of the Einstein equations is outlined, and artificial viscosity and numerical diffusion are discussed. The time integration is performed with AGILE, an implicit solver for stiff algebrodifferential equations on a dynamical adaptive grid. We extend the adaptive grid technique, known from nonrelativistic hydrodynamics, to the general relativistic application and identify it with the concept of shift vectors in a 3+1 decomposition. The adaptive grid minimizes the number of required computational zones without compromising the resolution in physically important regions. Thus, the computational effort is greatly reduced when the zones are subject to computationally expensive additional processes, such as Boltzmann radiation transport or a nuclear reaction network. We present accurate results in the standard tests for supernova simulations: Sedov's point-blast explosion, the nonrelativistic and relativistic shock tube, the Oppenheimer-Snyder dust collapse, and homologous collapse.

Astrophysics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to handle complex physical processes efficiently and accurately in general - relativistic hydrodynamics simulations under spherical symmetry conditions. Specifically, the authors developed an implicit solver for the numerical simulation of general - relativistic hydrodynamics and introduced an adaptive mesh technique to improve computational efficiency and accuracy. The following is a detailed explanation of the main problems addressed in the paper:

1. **Development of the implicit solver**:
- **Time integration**: The paper proposes an implicit solver (AGILE) for solving stiff algebraic - differential equations. The implicit method can handle physical equilibrium states on long - time scales while ensuring a smooth transition between the dynamic and equilibrium phases.
- **Stability and accuracy**: The implicit method can handle stiff systems, avoiding the instabilities and inaccuracies that explicit methods may exhibit at large time steps.

2. **Adaptive mesh technique**:
- **Reducing computational load**: Through the adaptive mesh technique, more grid points can be concentrated in physically important regions, thereby reducing the total number of grid points required for computation and improving computational efficiency.
- **Dynamic adjustment**: The adaptive mesh can be adjusted dynamically as the physical state changes, ensuring high resolution in key regions without wasting computational resources in regions of known physical equilibrium.

3. **Application in general relativity**:
- **Coordinate selection**: The paper combines the adaptive mesh technique with the concept of shift vectors in general relativity to implement the adaptive mesh technique in the 3 + 1 decomposition.
- **Equations in conservative form**: The paper derives the general - relativistic hydrodynamics equations in conservative form in the comoving coordinate system, ensuring the conservation properties of physical quantities.

4. **Standard tests**:
- **Verification and comparison**: The paper verifies the accuracy and effectiveness of the developed method through standard tests (such as Sedov point explosion, non - relativistic and relativistic shock tubes, Oppenheimer - Snyder dust collapse, homogeneous collapse, etc.) and compares it with analytical solutions and other numerical methods (such as PPM).

In summary, the main objective of this paper is to develop an efficient implicit solver and an adaptive mesh technique for handling complex physical processes in general - relativistic hydrodynamics, especially in astrophysical problems such as supernova simulations. These techniques not only improve computational efficiency but also ensure the stability and accuracy of the numerical solutions.