English

Finding Apparent Horizons in Numerical Relativity

General Relativity and Quantum Cosmology 2009-07-10 v4 Astrophysics

Abstract

This paper presents a detailed discussion of the ``Newton's method'' algorithm for finding apparent horizons in 3+1 numerical relativity. We describe a method for computing the Jacobian matrix of the finite differenced H(h)H(h) function by symbolically differentiating the finite difference equations, giving the Jacobian elements directly in terms of the finite difference molecule coefficients used in computing H(h)H(h). Assuming the finite differencing scheme commutes with linearization, we show how the Jacobian elements may be computed by first linearizing the continuum H(h)H(h) equations, then finite differencing the linearized (continuum) equations. We find this symbolic differentiation method of computing the H(h)H(h) Jacobian to be {\em much} more efficient than the usual numerical perturbation method, and also much easier to implement than is commonly thought. When solving the discrete H(h)=0H(h) = 0 equations, we find that Newton's method generally converges very rapidly. However, if the initial guess for the horizon position contains significant high-spatial-frequency error components, Newton's method has a small (poor) radius of convergence. This is {\em not} an artifact of insufficient resolution in the finite difference grid; rather, it appears to be caused by a strong nonlinearity in the continuum H(h)H(h) function for high-spatial-frequency error components in hh. Robust variants of Newton's method can boost the radius of convergence by O(1) factors, but the underlying nonlinearity remains, and appears to worsen rapidly with increasing initial-guess-error spatial frequency. Using 4th~order finite differencing, we find typical accuracies for computed horizon positions in the 10510^{-5} range for Δθ=π/250\Delta\theta = \frac{\pi/2}{50}.

Keywords

Cite

@article{arxiv.gr-qc/9508014,
  title  = {Finding Apparent Horizons in Numerical Relativity},
  author = {Jonathan Thornburg},
  journal= {arXiv preprint arXiv:gr-qc/9508014},
  year   = {2009}
}

Comments

38 pages text + postscript rotated-table (1 page) + postscript figures (10 pages), REVTeX macros, slightly revised from previous versions to clarify presentation, files also available from ftp://black_hole.physics.ubc.ca/outgoing/thornburg/horizon.tex.gz ftp://black_hole.physics.ubc.ca/outgoing/thornburg/horizon.ps.gz ftp://black_hole.physics.ubc.ca/outgoing/thornburg/horizon_tabsr.ps.gz ftp://black_hole.physics.ubc.ca/outgoing/thornburg/horizon_figs.ps.gz