Numerical initial boundary value problem for the generalized conformal field equations

Abstract

In this paper we study a numerical implementation for the initial boundary value formulation for the generalized conformal field equations. We propose a formulation which is well suited for the study of the long-time behaviour of perturbed exact solutions such as a Schwarzschild or even a Kerr black hole. We describe the derivation of the implemented equations which we give in terms of the space-spinor formalism. We discuss the conformal Gauss gauge, and a slight generalization thereof which seems to be particularly useful in the presence of boundaries. We discuss the structure of the equations at the boundary and propose a method for imposing boundary conditions which allow the correct number of degrees of freedom to be freely specified while still preserving the constraints. We show that this implementation yields a numerically well-posed system by testing it on a simple case of gravitational perturbations of Minkowski space-time and subsequently with gravitational perturbations of Schwarzschild space-time. * fbeyer@maths.otago.ac.nz

A new look at the Bondi-Sachs energy-momentum

“…This situation occurs in numerical relativity, when one solves equations which allow full access to null-infinity but which deny the possibility to introduce gauges adapted to I because of some numerically relevant considerations. An example which was the motivation for this work is described in [7].…”

Can Gravitational Waves Halt the Expansion of the Universe?

“…Boundary conditions are imposed using the Simultaneous Approximation Term (SAT) method [27] with τ = 1. This particular selection of numerical methods within COFFEE has proven to be numerically sound for a variety of different systems (see for example [19,28]). In the subsequent situations, all constraints are verified to converge at the expected order everywhere.…”

COFFEE—An MPI-parallelized Python package for the numerical evolution of differential equations

“…COFFEE was specifically developed to compute solutions to a system of hyperbolic partial differential equations (PDEs) that represent Friedrich's conformal field equations [2]. It has been used in eight research projects to numerically study the conformal properties of general relativity, [3,4,5,6,7,8,9,10]. As an illustration of the capabilities of COFFEE, in [10] it was used to solve a system of PDEs in the form of an Initial Boundary Value Problem (IBVP) containing 46 variables and 45 constraints on two different high performance clusters using up to 200 processes.…”

“…For example; as Python is an interpreted language syntax can have a large impact on speed of execution (compare for loops to list comprehensions), there are structural issues with interpreted languages (in the case of Python this is the reason for the GIL), additional overhead in the translation of Python script to Python byte code and then to machine instructions, and a lack of the compile time checks that are found in strongly typed languages. Nevertheless the papers [3,4,5,6,7,8,9,10] demonstrate that COFFEE is capable of solving technically challenging and computationally intensive systems of PDEs. Implementation in Python also has advantages.…”