Gmunu: toward multigrid based Einstein field equations solver for general-relativistic hydrodynamics simulations

We present a new open-source axisymmetric general relativistic hydrodynamics code Gmunu (general-relativistic multigrid numerical solver) which uses a multigrid method to solve the elliptic metric equations in the conformally flat condition (CFC) approximation on a spherical grid. Most of the existi...

Full description

Saved in:
Bibliographic Details
Published in:Classical and quantum gravity 2020-07, Vol.37 (14), p.145015
Main Authors: Cheong, Patrick Chi-Kit, Lin, Lap-Ming, Li, Tjonnie Guang Feng
Format: Article
Language:English
Subjects:
conformally flat condition
general-relativistic hydrodynamics
multigrid
numerical relativity
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:We present a new open-source axisymmetric general relativistic hydrodynamics code Gmunu (general-relativistic multigrid numerical solver) which uses a multigrid method to solve the elliptic metric equations in the conformally flat condition (CFC) approximation on a spherical grid. Most of the existing relativistic hydrodynamics codes are based on formulations which rely on a free-evolution approach of numerical relativity, where the metric variables are determined by hyperbolic equations without enforcing the constraint equations in the evolution. On the other hand, although a fully constrained-evolution formulation is theoretical more appealing and should lead to more stable and accurate simulations, such an approach is not widely used because solving the elliptic-type constraint equations during the evolution is in general more computationally expensive than hyperbolic free-evolution schemes. Multigrid methods solve differential equations with a hierarchy of discretizations and its computational cost is generally lower than other methods such as direct methods, relaxation methods, successive over-relaxation. With multigrid acceleration, one can solve the metric equations on a comparable time scale as solving the hydrodynamics equations. This would potentially make a fully constrained-evolution formulation more affordable in numerical relativity simulations. As a first step to assess the performance and robustness of multigrid methods in relativistic simulations, we develop a hydrodynamics code that makes use of standard finite-volume methods coupled with a multigrid metric solver to solve the Einstein equations in the CFC approximation. In this paper, we present the methodology and implementation of our code Gmunu and its properties and performance in some benchmarking relativistic hydrodynamics problems.
ISSN:0264-9381
1361-6382
DOI:10.1088/1361-6382/ab8e9c