Abstract:In this paper we develop a new well-balanced discontinuous Galerkin (DG) finite element scheme with subcell finite volume (FV) limiter for the numerical solution of the Einstein–Euler equations of general relativity based on a first order hyperbolic reformulation of the Z4 formalism. The first order Z4 system, which is composed of 59 equations, is analyzed and proven to be strongly hyperbolic for a general metric. The well-balancing is achieved for arbitrary but a priori known equilibria by subtracting a discrete version of the equilibrium solution from the discretized time-dependent PDE system. Special care has also been taken in the design of the numerical viscosity so that the well-balancing property is achieved. As for the treatment of low density matter, e.g. when simulating massive compact objects like neutron stars surrounded by vacuum, we have introduced a new filter in the conversion from the conserved to the primitive variables, preventing superluminal velocities when the density drops below a certain threshold, and being potentially also very useful for the numerical investigation of highly rarefied relativistic astrophysical flows. Thanks to these improvements, all standard tests of numerical relativity are successfully reproduced, reaching three achievements: (i) we are able to obtain stable long term simulations of stationary black holes, including Kerr black holes with extreme spin, which after an initial perturbation return perfectly back to the equilibrium solution up to machine precision; (ii) a (standard) TOV star under perturbation is evolved in pure vacuum ( ρ=p=0 ) up to t=1000 with no need to introduce any artificial atmosphere around the star; and, (iii) we solve the head on collision of two punctures black holes, that was previously considered un–tractable within the Z4 formalism. Due to the above features, we consider that our new algorithm can be particularly beneficial for the numerical study of quasi normal modes of oscillations, both of black holes and of neutron stars.

Implementation of higher-order absorbing boundary conditions for the Einstein equations

Outer boundary conditions for Einstein's field equations in harmonic coordinates

Boundary conditions for Einstein's field equations: Analytical and numerical analysis

Peptide Conformational Equilibria Computed via a Single-Stage Shifting Protocol

Numerical Relativity Multimodal Waveforms using Absorbing Boundary Conditions

Numerical Tests of Evolution Systems, Gauge Conditions, and Boundary Conditions for 1D Colliding Gravitational Plane Waves

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

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

On boundary conditions for linearised Einstein's equations

Uniform Boundedness for Solutions to the Teukolsky Equation on Schwarzschild from Conservation Laws of Linearised Gravity

Analysis of High-Order Absorbing Boundary Conditions for the Schrödinger Equation

Absorbing Boundary Condition for Nonlinear Euler Equations in Primitive Variables Based on the Perfectly Matched Layer Technique

Formulating the complete initial boundary value problem in numerical relativity to model black hole echoes

Efficient implementation of high order inverse Lax–Wendroff boundary treatment for conservation laws

Adaptive absorbing boundary conditions for Schrödinger-type equations: Application to nonlinear and multi-dimensional problems

A monolithic first--order BSSNOK formulation of the Einstein--Euler equations and its solution with path-conservative finite difference CWENO schemes

Second-order perturbations of the Schwarzschild spacetime: practical, covariant and gauge-invariant formalisms

Second-order perturbations of the Schwarzschild spacetime: Practical, covariant, and gauge-invariant formalisms

Adaptive absorbing boundary conditions for Schrodinger-type equations: application to nonlinear and multi-dimensional problems

Absorbing Boundary Conditions for the Multidimensional Klein-Gordon Equation

Schwarzschild Tests of the Wahlquist-Estabrook-Buchman-Bardeen Tetrad Formulation for Numerical Relativity