On the iterated Crank-Nicolson for hyperbolic and parabolic equations in numerical relativity

Abstract

The iterated Crank-Nicolson is a predictor-corrector algorithm commonly used in numerical relativity for the solution of both hyperbolic and parabolic partial differential equations. We here extend the recent work on the stability of this scheme for hyperbolic equations by investigating the properties when the average between the predicted and corrected values is made with unequal weights and when the scheme is applied to a parabolic equation. We also propose a variant of the scheme in which the coefficients in the averages are swapped between two corrections leading to systematically larger amplification factors and to a smaller numerical dispersion.

7 pages, 3 figures

Cited by in corpus (4)

• Accurate evolutions of inspiralling neutron-star binaries: prompt and delayed collapse to black hole
• WhiskyMHD: a new numerical code for general relativistic magnetohydrodynamics
• Gravitational-Wave Extraction from Neutron Star Oscillations: comparing linear and nonlinear techniques
• An improved formulation of the relativistic hydrodynamics equations in 2D Cartesian coordinates