skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Iterated Crank-Nicolson method 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.

Authors:
 [1];  [2];  [3];  [4];  [5]
  1. Department of Physics, Udine University, Udine (Italy)
  2. (Germany)
  3. Max-Planck-Institut fuer Gravitationsphysik, Albert Einstein Institut, Golm (Germany)
  4. (Italy)
  5. (United States)
Publication Date:
OSTI Identifier:
20776739
Resource Type:
Journal Article
Resource Relation:
Journal Name: Physical Review. D, Particles Fields; Journal Volume: 73; Journal Issue: 4; Other Information: DOI: 10.1103/PhysRevD.73.044001; (c) 2006 The American Physical Society; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS; ALGORITHMS; AMPLIFICATION; CORRECTIONS; COSMOLOGY; GENERAL RELATIVITY THEORY; ITERATIVE METHODS; MATHEMATICAL SOLUTIONS; PARTIAL DIFFERENTIAL EQUATIONS; STABILITY

Citation Formats

Leiler, Gregor, Max-Planck-Institut fuer Gravitationsphysik, Albert Einstein Institut, Golm, Rezzolla, Luciano, SISSA, International School for Advanced Studies and INFN, Trieste, and Department of Physics and Astronomy, Louisiana State University, Baton Rouge, Louisiana 70803-4001. Iterated Crank-Nicolson method for hyperbolic and parabolic equations in numerical relativity. United States: N. p., 2006. Web. doi:10.1103/PhysRevD.73.044001.
Leiler, Gregor, Max-Planck-Institut fuer Gravitationsphysik, Albert Einstein Institut, Golm, Rezzolla, Luciano, SISSA, International School for Advanced Studies and INFN, Trieste, & Department of Physics and Astronomy, Louisiana State University, Baton Rouge, Louisiana 70803-4001. Iterated Crank-Nicolson method for hyperbolic and parabolic equations in numerical relativity. United States. doi:10.1103/PhysRevD.73.044001.
Leiler, Gregor, Max-Planck-Institut fuer Gravitationsphysik, Albert Einstein Institut, Golm, Rezzolla, Luciano, SISSA, International School for Advanced Studies and INFN, Trieste, and Department of Physics and Astronomy, Louisiana State University, Baton Rouge, Louisiana 70803-4001. Wed . "Iterated Crank-Nicolson method for hyperbolic and parabolic equations in numerical relativity". United States. doi:10.1103/PhysRevD.73.044001.
@article{osti_20776739,
title = {Iterated Crank-Nicolson method for hyperbolic and parabolic equations in numerical relativity},
author = {Leiler, Gregor and Max-Planck-Institut fuer Gravitationsphysik, Albert Einstein Institut, Golm and Rezzolla, Luciano and SISSA, International School for Advanced Studies and INFN, Trieste and Department of Physics and Astronomy, Louisiana State University, Baton Rouge, Louisiana 70803-4001},
abstractNote = {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.},
doi = {10.1103/PhysRevD.73.044001},
journal = {Physical Review. D, Particles Fields},
number = 4,
volume = 73,
place = {United States},
year = {Wed Feb 15 00:00:00 EST 2006},
month = {Wed Feb 15 00:00:00 EST 2006}
}
  • No abstract prepared.
  • A stability analysis of the classical Crank-Nicolson-Galerkin (CNG) scheme applied to the one-dimensional solute transport equation is proposed on the basis of two fairly different approaches. Using a space-time eigenvalue analysis, it is shown, at least for subsurface hydrology applications, that the CNG scheme is theoretically stable under the condition PeCr [le] 2, where Pe and Cr are the mesh Peclet and Courant numbers. Then, to assess the computational stability of the scheme, the amplification matrix is constructed, and its norm is displayed in the (Pe, Cr) space. The results indicate that the norm of the amplification matrix is nevermore » smaller than unity and exhibits a hyperbolic nature analogous to the above theoretical condition. 21 refs., 2 figs.« less
  • We propose and analyze a linear stabilization of the Crank-Nicolson Leapfrog (CNLF) method that removes all time step/CFL conditions for stability and controls the unstable mode. It also increases the SPD part of the linear system to be solved at each time step while increasing solution accuracy. We give a proof of unconditional stability of the method as well as a proof of unconditional, asymptotic stability of both the stable and unstable modes. As a result, we illustrate two applications of the method: uncoupling groundwater-surface water flows and Stokes flow plus a Coriolis term.
  • Two new formulations of general relativity are introduced. The first one is a parabolization of the Arnowitt-Deser-Misner formulation and is derived by the addition of combinations of the constraints and their derivatives to the right-hand side of the Arnowitt-Deser-Misner evolution equations. The desirable property of this modification is that it turns the surface of constraints into a local attractor because the constraint propagation equations become second-order parabolic independently of the gauge conditions employed. This system may be classified as mixed hyperbolic--second-order parabolic. The second formulation is a parabolization of the Kidder-Scheel-Teukolsky formulation and is a manifestly mixed strongly hyperbolic--second-order-parabolic setmore » of equations, bearing thus resemblance to the compressible Navier-Stokes equations. As a first test, a stability analysis of flat space is carried out and it is shown that the first modification exponentially damps and smoothes all constraint-violating modes. These systems provide a new basis for constructing schemes for long-term and stable numerical integration of the Einstein field equations.« less