A remedy for constraint growth in numerical relativity: the Maxwell case
A remedy for constraint growth in numerical relativity: the Maxwell case
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.
4025-4040
Calabrese, G.
ffefb9b7-fca3-41ac-89ea-1a7a85dbd7d9
2004
Calabrese, G.
ffefb9b7-fca3-41ac-89ea-1a7a85dbd7d9
Calabrese, G.
(2004)
A remedy for constraint growth in numerical relativity: the Maxwell case.
Classical and Quantum Gravity, 21, .
(doi:10.1088/0264-9381/21/17/002).
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.
This record has no associated files available for download.
More information
Published date: 2004
Identifiers
Local EPrints ID: 40852
URI: http://eprints.soton.ac.uk/id/eprint/40852
ISSN: 0264-9381
PURE UUID: 266c13c1-fb6b-4342-83c9-800e6b4964f4
Catalogue record
Date deposited: 11 Jul 2006
Last modified: 15 Mar 2024 08:23
Export record
Altmetrics
Contributors
Author:
G. Calabrese
Download statistics
Downloads from ePrints over the past year. Other digital versions may also be available to download e.g. from the publisher's website.
View more statistics