Testing numerical relativity with the shifted gauge wave
arXiv:gr-qc/0511154 · doi:10.1088/0264-9381/23/16/S03
Abstract
Computational methods are essential to provide waveforms from coalescing black holes, which are expected to produce strong signals for the gravitational wave observatories being developed. Although partial simulations of the coalescence have been reported, scientifically useful waveforms have so far not been delivered. The goal of the AppleswithApples (AwA) Alliance is to design, coordinate and document standardized code tests for comparing numerical relativity codes. The first round of AwA tests have now being completed and the results are being analyzed. These initial tests are based upon periodic boundary conditions designed to isolate performance of the main evolution code. Here we describe and carry out an additional test with periodic boundary conditions which deals with an essential feature of the black hole excision problem, namely a non-vanishing shift. The test is a shifted version of the existing AwA gauge wave test. We show how a shift introduces an exponentially growing instability which violates the constraints of a standard harmonic formulation of Einstein's equations. We analyze the Cauchy problem in a harmonic gauge and discuss particular options for suppressing instabilities in the gauge wave tests. We implement these techniques in a finite difference evolution algorithm and present test results. Although our application here is limited to a model problem, the techniques should benefit the simulation of black holes using harmonic evolution codes.
Submitted to special numerical relativity issue of Classical and Quantum Gravity
References in corpus (16)
- Evolution of Binary Black Hole Spacetimes
- Numerical Relativity Using a Generalized Harmonic Decomposition
- Constraint damping in the Z4 formulation and harmonic gauge
- General-covariant evolution formalism for Numerical Relativity
- Well-Posed Initial-Boundary Evolution in General Relativity
- Toward standard testbeds for numerical relativity
- Gravitational wave extraction and outer boundary conditions by perturbative matching
- Stable characteristic evolution of generic 3-dimensional single-black-hole spacetimes
- Summation by parts and dissipation for domains with excised regions
- Implementation of standard testbeds for numerical relativity
- Harmonic Initial-Boundary Evolution in General Relativity
- 3D simulations of Einstein's equations: symmetric hyperbolicity, live gauges and dynamic control of the constraints
- Modeling the Black Hole Excision Problem
- Discrete boundary treatment for the shifted wave equation
- The discrete energy method in numerical relativity: Towards long-term stability
- Finite differencing second order systems describing black hole spacetimes
Cited by in corpus (16)
- A New Generalized Harmonic Evolution System
- Simulation of Binary Black Hole Spacetimes with a Harmonic Evolution Scheme
- Continuum and Discrete Initial-Boundary-Value Problems and Einstein's Field Equations
- An explicit harmonic code for black-hole evolution using excision
- Key Elements of Robustness in Binary Black Hole Evolutions using Spectral Methods
- Problems which are well-posed in a generalized sense with applications to the Einstein equations
- Implementation of standard testbeds for numerical relativity
- Harmonic Initial-Boundary Evolution in General Relativity
- From Geometry to Numerics: interdisciplinary aspects in mathematical and numerical relativity
- Testing the Accuracy and Stability of Spectral Methods in Numerical Relativity
- Long term black hole evolution with the BSSN system by pseudospectral methods
- Constraint-preserving Sommerfeld conditions for the harmonic Einstein equations
- Constraint-preserving boundary treatment for a harmonic formulation of the Einstein equations
- Finite Difference Methods for Second Order in Space, First Order in Time Hyperbolic Systems and the Linear Shifted Wave Equation as a Model Problem in Numerical Relativity
- Boundary Conditions for the Gravitational Field
- Numerical performance of the parabolized ADM (PADM) formulation of General Relativity