Stable 3-level leapfrog integration in numerical relativity

Details Stable 3-level leapfrog integration in numerical relativity Stable 3-level leapfrog integration in numerical relativity

1998-01-30

arxiv.org/abs/gr-qc/9801110v2

Phys.Rev. D58 (1998) 064022

Astronomy and Astrophysics

Astrophysics

General Relativity and Quantum Cosmology

14 pages and 10 postscript figures; different Hamiltonian constraint normalization used, minor alterations to cost analysis an

Scientific paper

10.1103/PhysRevD.58.064022

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.

Affiliated with

Also associated with

No associations

Rating

Stable 3-level leapfrog integration in numerical relativity does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.

If you have personal experience with Stable 3-level leapfrog integration in numerical relativity, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Stable 3-level leapfrog integration in numerical relativity will most certainly appreciate the feedback.

Rate now

Comments { 0 }

Profile ID: LFWR-SCP-O-554550