Numerical stability for finite difference approximations of Einstein's equations
arXiv:gr-qc/0503056 · doi:10.1016/j.jcp.2006.02.027
Abstract
We extend the notion of numerical stability of finite difference approximations to include hyperbolic systems that are first order in time and second order in space, such as those that appear in Numerical Relativity. By analyzing the symbol of the second order system, we obtain necessary and sufficient conditions for stability in a discrete norm containing one-sided difference operators. We prove stability for certain toy models and the linearized Nagy-Ortiz-Reula formulation of Einstein's equations. We also find that, unlike in the fully first order case, standard discretizations of some well-posed problems lead to unstable schemes and that the Courant limits are not always simply related to the characteristic speeds of the continuum problem. Finally, we propose methods for testing stability for second order in space hyperbolic systems.
18 pages, 9 figures
References in corpus (19)
- On the Numerical Integration of Einstein's Field Equations
- General-covariant evolution formalism for Numerical Relativity
- Extending the lifetime of 3D black hole computations with a new hyperbolic system of evolution equations
- On the Stability of the Iterated Crank-Nicholson Method in Numerical Relativity
- Fixing Einstein's equations
- Toward standard testbeds for numerical relativity
- First-order symmetric-hyperbolic Einstein equations with arbitrary fixed gauge
- Strongly hyperbolic second order Einstein's evolution equations
- A symmetry-breaking mechanism for the Z4 general-covariant evolution system
- Hyperbolicity of the BSSN system of Einstein evolution equations
- On the well posedness of the Baumgarte-Shapiro-Shibata-Nakamura formulation of Einstein's field equations
- Symmetric hyperbolicity and consistent boundary conditions for second-order Einstein equations
- Cauchy boundaries in linearized gravitational theory
- Symmetric hyperbolic form of systems of second-order evolution equations subject to constraints
- Ill-posedness in the Einstein equations
- Illustrating Stability Properties of Numerical Relativity in Electrodynamics
- Modeling the Black Hole Excision Problem
- Stability properties of a formulation of Einstein's equations
- Finite differencing second order systems describing black hole spacetimes
Cited by in corpus (27)
- Calibration of Moving Puncture Simulations
- Constraint violation in free evolution schemes: comparing BSSNOK with a conformal decomposition of Z4
- Numerical solution of gravitational dynamics in asymptotically anti-de Sitter spacetimes
- Reducing phase error in long numerical binary black hole evolutions with sixth order finite differencing
- Continuum and Discrete Initial-Boundary-Value Problems and Einstein's Field Equations
- Implementation of standard testbeds for numerical relativity
- Spherical symmetry as a test case for unconstrained hyperboloidal evolution
- Numerical stability of the Z4c formulation of general relativity
- Numerical simulations with a first order BSSN formulation of Einstein's field equations
- Constraint-preserving boundary treatment for a harmonic formulation of the Einstein equations
- Numerical-relativity simulations of the quasi-circular inspiral and merger of non-spinning, charged black holes: methods and comparison with approximate approaches
- Spherical symmetry as a test case for unconstrained hyperboloidal evolution II: gauge conditions
- Fully relativistic three-dimensional Cauchy-characteristic matching for physical degrees of freedom
- Hyperbolicity of General Relativity in Bondi-like gauges
- Nonlinear Dynamics of Quadratic Gravity in Spherical Symmetry
- Stability of the puncture method with a generalized BSSN formulation
- Gauge structure of the Einstein field equations in Bondi-like coordinates
- Numerical modeling of black holes as sources of gravitational waves in a nutshell
- Finite Difference Methods for Second Order in Space, First Order in Time Hyperbolic Systems and the Linear Shifted Wave Equation as a Model Problem in Numerical Relativity
- Numerical convergence of model Cauchy-characteristic extraction and matching
- On the Hyperbolicity and Stability of Formulations of Metric Gravity
- Free and constrained symplectic integrators for numerical general relativity
- Hyperbolicity of High Order Systems of Evolution Equations
- Strongly hyperbolic Hamiltonian systems in numerical relativity: Formulation and symplectic integration
- Well-balanced high order finite difference WENO schemes for a first-order Z4 formulation of the Einstein field equations
- Quantum correlations in a gravitational collapse simulation with SpheriCo.jl
- Extended black hole cosmologies in de Sitter space