Constraint and gauge shocks in one-dimensional numerical relativity

Abstract

We study how different types of blow-ups can occur in systems of hyperbolic evolution equations of the type found in general relativity. In particular, we discuss two independent criteria that can be used to determine when such blow-ups can be expected. One criteria is related with the so-called geometric blow-up leading to gradient catastrophes, while the other is based upon the ODE-mechanism leading to blow-ups within finite time. We show how both mechanisms work in the case of a simple one-dimensional wave equation with a dynamic wave speed and sources, and later explore how those blow-ups can appear in one-dimensional numerical relativity. In the latter case we recover the well known ``gauge shocks'' associated with Bona-Masso type slicing conditions. However, a crucial result of this study has been the identification of a second family of blow-ups associated with the way in which the constraints have been used to construct a hyperbolic formulation. We call these blow-ups ``constraint shocks'' and show that they are formulation specific, and that choices can be made to eliminate them or at least make them less severe.

19 pages, 8 figures and 1 table, revised version including several amendments suggested by the referee

References in corpus (13)

• On the Numerical Integration of Einstein's Field Equations
• New Formalism for Numerical Relativity
• Extending the lifetime of 3D black hole computations with a new hyperbolic system of evolution equations
• Toward standard testbeds for numerical relativity
• Towards an understanding of the stability properties of the 3+1 evolution equations in general relativity
• Hyperbolicity of the BSSN system of Einstein evolution equations
• Hyperbolic slicings of spacetime: singularity avoidance and gauge shocks
• The appearance of coordinate shocks in hyperbolic formalisms of General Realtivity
• On the well posedness of the Baumgarte-Shapiro-Shibata-Nakamura formulation of Einstein's field equations
• Ill-posedness in the Einstein equations
• Exploiting gauge and constraint freedom in hyperbolic formulations of Einstein's equations
• Energy Norms and the Stability of the Einstein Evolution Equations
• Adjusted ADM systems and their expected stability properties: constraint propagation analysis in Schwarzschild spacetime

Cited by in corpus (7)

• Continuum and Discrete Initial-Boundary-Value Problems and Einstein's Field Equations
• Well-posedness of formulations of the Einstein equations with dynamical lapse and shift conditions
• Numerical relativity for Horndeski gravity
• Spherical symmetry as a test case for unconstrained hyperboloidal evolution II: gauge conditions
• Are gauge shocks really shocks?
• Generalized harmonic spatial coordinates and hyperbolic shift conditions
• Matters of Gravity, The Newsletter of The Topical Group in Gravitation of the American Physical Society, Volume 27, Spring 2006