Finite differencing second order systems describing black hole spacetimes

Abstract

Keeping Einstein's equations in second order form can be appealing for computational efficiency, because of the reduced number of variables and constraints. Stability issues emerge, however, which are not present in first order formulations. We show that a standard discretization of the second order ``shifted'' wave equation leads to an unstable semi-discrete scheme if the shift parameter is too large. This implies that discretizations obtained using integrators such as Runge-Kutta, Crank-Nicholson, leap-frog are unstable for any fixed value of the Courant factor. We argue that this situation arises in numerical relativity, particularly in simulations of spacetimes containing black holes, and discuss several ways of circumventing this problem. We find that the first order reduction in time based on ``ADM'' type variables is very effective.

Cited by in corpus (12)

• Numerical Relativity Using a Generalized Harmonic Decomposition
• Simulation of Binary Black Hole Spacetimes with a Harmonic Evolution Scheme
• An explicit harmonic code for black-hole evolution using excision
• Harmonic Initial-Boundary Evolution in General Relativity
• Hyperbolicity of second-order in space systems of evolution equations
• Testing numerical relativity with the shifted gauge wave
• Numerical stability for finite difference approximations of Einstein's equations
• Modeling the Black Hole Excision Problem
• Discrete boundary treatment for the shifted wave equation
• Finite difference schemes for second order systems describing black holes
• Excising a boosted rotating black hole with overlapping grids
• Exact boundary conditions in numerical relativity using multiple grids: scalar field tests