A remedy for constraint growth in Numerical Relativity

Abstract

Rapid growth of constraints is often observed in free evolutions of highly gravitating systems. To alleviate this problem we investigate the effect of adding spatial derivatives of the constraints to the right hand side of the evolution equations, and we look at how this affects the character of the system and the treatment of boundaries. We apply this technique to two formulations of Maxwell's equations, the so-called fat Maxwell and the Knapp-Walker-Baumgarte systems, and obtain mixed hyperbolic-parabolic problems in which high frequency constraint violations are damped. Constraint-preserving boundary conditions amount to imposing Dirichlet boundary conditions on constraint variables, which translate into Neumann-like boundary conditions for the main variables. The success of the numerical tests presented in this work suggests that this remedy may bring benefits to fully nonlinear simulations of General Relativity.

10 pages, 5 figures

References in corpus (14)

• On the Numerical Integration of Einstein's Field Equations
• Simulation of merging binary neutron stars in full general relativity: $\Gamma =2$ case
• Fully general relativistic simulation of coalescing binary neutron stars: Preparatory tests
• Simple excision of a black hole in 3+1 numerical relativity
• Well-Posed Initial-Boundary Evolution in General Relativity
• Einstein's Equations with Asymptotically Stable Constraint Propagation
• Well posed constraint-preserving boundary conditions for the linearized Einstein equations
• Constraint-preserving boundary conditions in numerical relativity
• Coordinate Conditions and Their Implementation in 3D Numerical Relativity
• Toward stable 3D numerical evolutions of black-hole spacetimes
• Illustrating Stability Properties of Numerical Relativity in Electrodynamics
• On free evolution of self gravitating, spherically symmetric waves
• Constrained Transport Algorithms for Numerical Relativity. I. Development of a Finite Difference Scheme
• Toward Making the Constraint Hypersurface an Attractor in Free Evolution

Cited by in corpus (9)

• Constraint damping in the Z4 formulation and harmonic gauge
• Constraint-preserving boundary conditions in the Z4 Numerical Relativity formalism
• A model problem for the initial-boundary value formulation of Einstein's field equations
• Discrete boundary treatment for the shifted wave equation
• Constraint propagation equations of the 3+1 decomposition of f(R) gravity
• Excising a boosted rotating black hole with overlapping grids
• Mixed Hyperbolic - Second-Order Parabolic Formulations of General Relativity
• A polyhedral discrete de Rham numerical scheme for the Yang-Mills equations
• On constraint preservation and strong hyperbolicity