Extended Lifetime in Computational Evolution of Isolated Black Holes
arXiv:gr-qc/0307055 · doi:10.1007/s10701-005-6477-x
Abstract
Solving the 4-d Einstein equations as evolution in time requires solving equations of two types: the four elliptic initial data (constraint) equations, followed by the six second order evolution equations. Analytically the constraint equations remain solved under the action of the evolution, and one approach is to simply monitor them ({\it unconstrained} evolution). The problem of the 3-d computational simulation of even a single isolated vacuum black hole has proven to be remarkably difficult. Recently, we have become aware of two publications that describe very long term evolution, at least for single isolated black holes. An essential feature in each of these results is {\it constraint subtraction}. Additionally, each of these approaches is based on what we call "modern," hyperbolic formulations of the Einstein equations. It is generally assumed, based on computational experience, that the use of such modern formulations is essential for long-term black hole stability. We report here on comparable lifetime results based on the much simpler ("traditional") - formulation. We have also carried out a series of {\it constrained} 3-d evolutions of single isolated black holes. We find that constraint solution can produce substantially stabilized long-term single hole evolutions. However, we have found that for large domains, neither constraint-subtracted nor constrained - evolutions carried out in Cartesian coordinates admit arbitrarily long-lived simulations. The failure appears to arise from features at the inner excision boundary; the behavior does generally improve with resolution.
20 pages, 6 figures
References in corpus (7)
- A constrained scheme for Einstein equations based on Dirac gauge and spherical coordinates
- Strongly hyperbolic second order Einstein's evolution equations
- Improved numerical stability of stationary black hole evolution calculations
- An Axisymmetric Gravitational Collapse Code
- Optimal Constraint Projection for Hyperbolic Evolution Systems
- Novel finite-differencing techniques for numerical relativity: application to black hole excision
- Toward stable 3D numerical evolutions of black-hole spacetimes
Cited by in corpus (11)
- Numerical Relativity Using a Generalized Harmonic Decomposition
- From Geometry to Numerics: interdisciplinary aspects in mathematical and numerical relativity
- Adaptive Mesh Refinement for Coupled Elliptic-Hyperbolic Systems
- Mixed Hyperbolic - Second-Order Parabolic Formulations of General Relativity
- Numerical performance of the parabolized ADM (PADM) formulation of General Relativity
- Boundary conditions for the Einstein-Christoffel formulation of Einstein's equations
- Well-posed constrained evolution of 3+1 formulations of General Relativity
- Modifying the Einstein Equations off the Constraint Hypersuface
- Constraint Damping in First-Order Evolution Systems for Numerical Relativity
- Hyperbolicity and Constrained Evolution in Linearized Gravity
- Momentum constraint relaxation