Improved numerical stability of stationary black hole evolution calculations
arXiv:gr-qc/0209066 · doi:10.1103/PhysRevD.66.084026
Abstract
We experiment with modifications of the BSSN form of the Einstein field equations (a reformulation of the ADM equations) and demonstrate how these modifications affect the stability of numerical black hole evolution calculations. We use excision to evolve both non-rotating and rotating Kerr-Schild black holes in octant and equatorial symmetry, and without any symmetry assumptions, and obtain accurate and stable simulations for specific angular momenta J/M of up to about 0.9M.
13 pages, 11 figures, 1 typo in Eq. (20) corrected