A well balanced finite volume scheme for general relativity

08/06/2021
by   Elena Gaburro, et al.
0

In this work we present a novel second order accurate well balanced (WB) finite volume (FV) scheme for the solution of the general relativistic magnetohydrodynamics (GRMHD) equations and the first order CCZ4 formulation (FO-CCZ4) of the Einstein field equations of general relativity, as well as the fully coupled FO-CCZ4 + GRMHD system. These systems of first order hyperbolic PDEs allow to study the dynamics of the matter and the dynamics of the space-time according to the theory of general relativity. The new well balanced finite volume scheme presented here exploits the knowledge of an equilibrium solution of interest when integrating the conservative fluxes, the nonconservative products and the algebraic source terms, and also when performing the piecewise linear data reconstruction. This results in a rather simple modification of the underlying second order FV scheme, which, however, being able to cancel numerical errors committed with respect to the equilibrium component of the numerical solution, substantially improves the accuracy and long-time stability of the numerical scheme when simulating small perturbations of stationary equilibria. In particular, the need for well balanced techniques appears to be more and more crucial as the applications increase their complexity. We close the paper with a series of numerical tests of increasing difficulty, where we study the evolution of small perturbations of accretion problems and stable TOV neutron stars. Our results show that the well balancing significantly improves the long-time stability of the finite volume scheme compared to a standard one.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
03/12/2018

A second-order well-balanced finite volume scheme for Euler equations with gravity

We present a well-balanced, second order, Godunov-type finite volume sch...
research
07/13/2023

A well-balanced discontinuous Galerkin method for the first–order Z4 formulation of the Einstein–Euler system

In this paper we develop a new well-balanced discontinuous Galerkin (DG)...
research
06/28/2023

A well-balanced semi-implicit IMEX finite volume scheme for ideal Magnetohydrodynamics at all Mach numbers

We propose a second-order accurate semi-implicit and well-balanced finit...
research
02/01/2021

A splitting scheme for the coupled Saint-Venant-Exner model

We present a splitting method for the one-dimensional Saint-Venant-Exner...
research
09/28/2019

Analysis of an Asymptotic Preserving Low Mach Number Accurate IMEX-RK Scheme for the Wave Equation System

In this paper the analysis of an asymptotic preserving (AP) IMEX-RK fini...
research
11/15/2020

Well-balanced algorithms for relativistic fluids on a Schwarzschild background

We design well-balanced numerical algorithms with first-, second-, and t...

Please sign up or login with your details

Forgot password? Click here to reset