Toward Making the Constraint Hypersurface an Attractor in Free Evolution

Abstract

There is an abundance of empirical evidence in the numerical relativity literature that the form in which the Einstein evolution equations are written plays a significant role in the lifetime of numerical simulations. This paper attempts to present a consistent framework for modifying any system of evolution equations by adding terms that push the evolution toward the constraint hypersurface. The method is, in principle, applicable to any system of partial differential equations which can be divided into evolution equations and constraints, although it is only demonstrated here through an application to the Maxwell equations.

6 pages, 3 figures, 1 table. Uses REVTeX4

References in corpus (2)

• On the Stability of the Iterated Crank-Nicholson Method in Numerical Relativity
• Illustrating Stability Properties of Numerical Relativity in Electrodynamics

Cited by in corpus (9)

• 3D simulations of Einstein's equations: symmetric hyperbolicity, live gauges and dynamic control of the constraints
• Constraint Propagation of $C^2$-adjusted Formulation - Another Recipe for Robust ADM Evolution System
• A model problem for the initial-boundary value formulation of Einstein's field equations
• Excising das All: Evolving Maxwell waves beyond scri
• Constraint propagation equations of the 3+1 decomposition of f(R) gravity
• Mixed Hyperbolic - Second-Order Parabolic Formulations of General Relativity
• A remedy for constraint growth in Numerical Relativity
• Numerical performance of the parabolized ADM (PADM) formulation of General Relativity
• Constraint Propagation of $C^2$-adjusted Formulation II --- Another Recipe for Robust Baumgarte-Shapiro-Shibata-Nakamura Evolution System ---