English

Stable 3-level leapfrog integration in numerical relativity

General Relativity and Quantum Cosmology 2009-10-31 v2 Astrophysics

Abstract

The 3-level leapfrog time integration algorithm is an attractive choice for numerical relativity simulations since it is time-symmetric and avoids non-physical damping. In Newtonian problems without velocity dependent forces, this method enjoys the advantage of long term stability. However, for more general differential equations, whether ordinary or partial, delayed onset numerical instabilities can arise and destroy the solution. A known cure for such instabilities appears to have been overlooked in many application areas. We give an improved cure ("deloused leapfrog") that both reduces memory demands (important for 3+1 dimensional wave equations) and allows for the use of adaptive timesteps without a loss in accuracy. We show both that the instability arises and that the cure we propose works in highly relativistic problems such as tightly bound geodesics, spatially homogeneous spacetimes, and strong gravitational waves. In the gravitational wave test case (polarized waves in a Gowdy spacetime) the deloused leapfrog method was five to eight times less CPU costly at various accuracies than the implicit Crank-Nicholson method, which is not subject to this instability.

Keywords

Cite

@article{arxiv.gr-qc/9801110,
  title  = {Stable 3-level leapfrog integration in numerical relativity},
  author = {Kimberly C. B. New and Keith Watt and Charles W. Misner and Joan M. Centrella},
  journal= {arXiv preprint arXiv:gr-qc/9801110},
  year   = {2009}
}

Comments

14 pages and 10 postscript figures; different Hamiltonian constraint normalization used, minor alterations to cost analysis and Conclusions section; accepted for publication in Physical Review D