On the Effect of Constraint Enforcement on the Quality of Numerical Solutions in General Relativity

Abstract

In Brodbeck et al 1999 it has been shown that the linearised time evolution equations of general relativity can be extended to a system whose solutions asymptotically approach solutions of the constraints. In this paper we extend the non-linear equations in similar ways and investigate the effect of various possibilities by numerical means. Although we were not able to make the constraint submanifold an attractor for all solutions of the extended system, we were able to significantly reduce the growth of the numerical violation of the constraints. Contrary to our expectations this improvement did not imply a numerical solution closer to the exact solution, and therefore did not improve the quality of the numerical solution.

14 pages, 9 figures, accepted for publication in Phys. Rev. D

Cited by in corpus (7)

• Constraint damping in the Z4 formulation and harmonic gauge
• Advantages of modified ADM formulation: constraint propagation analysis of Baumgarte-Shapiro-Shibata-Nakamura system
• Controlling the Growth of Constraints in Hyperbolic Evolution Systems
• Adjusted ADM systems and their expected stability properties: constraint propagation analysis in Schwarzschild spacetime
• On the Hyperbolicity and Stability of $3+1$ Formulations of Metric $f(R)$ Gravity
• Tips for implementing multigrid methods on domains containing holes
• Constraint Damping in First-Order Evolution Systems for Numerical Relativity