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
arXiv:0812.3752 · doi:10.1016/j.jcp.2009.12.016
Abstract
Motivated by the problem of solving the Einstein equations, we discuss high order finite difference discretizations of first order in time, second order in space hyperbolic systems.Particular attention is paid to the case when first order derivatives that can be identified with advection terms are approximated with non-centered finite difference operators.We first derive general properties of these discrete operators, then we extend a known result on numerical stability for such systems to general order of accuracy.As an application we analyze the shifted wave equation, including the behavior of the numerical phase and group speeds at different orders of approximations. Special attention is paid to when the use of off-centered schemes improves the accuracy over the centered schemes.
36 pages, 10 figures
References in corpus (10)
- Evolution of Binary Black Hole Spacetimes
- Accurate Evolutions of Orbiting Black-Hole Binaries Without Excision
- Gravitational wave extraction from an inspiraling configuration of merging black holes
- Accurate black hole evolutions by fourth-order numerical relativity
- Reducing phase error in long numerical binary black hole evolutions with sixth order finite differencing
- Well-posedness of formulations of the Einstein equations with dynamical lapse and shift conditions
- Status of black-hole-binary simulations for gravitational-wave detection
- Testing numerical relativity with the shifted gauge wave
- Discrete boundary treatment for the shifted wave equation
- Finite difference schemes for second order systems describing black holes
Cited by in corpus (10)
- Compact binary evolutions with the Z4c formulation
- SENR/NRPy+: Numerical Relativity in Singular Curvilinear Coordinate Systems
- GRAthena++: puncture evolutions on vertex-centered oct-tree AMR
- 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
- Stability of the puncture method with a generalized BSSN formulation
- Summation by Parts and Truncation Error Matching on Hyperboloidal Slices
- Spherically symmetric black hole spacetimes on hyperboloidal slices
- Free evolution of the hyperboloidal initial value problem in spherical symmetry