Action principle for Numerical Relativity evolution systems
Abstract
A Lagrangian density is provided, that allows to recover the Z4 evolution system from an action principle. The resulting system is then strongly hyperbolic when supplemented by gauge conditions like '1+log' or 'freezing shift', suitable for numerical evolution. The physical constraint can be imposed just on the initial data. The corresponding canonical equations are also provided. This opens the door to analogous results for other numerical-relativity formalisms, like BSSN, that can be derived from Z4 by a symmetry-breaking procedure. The harmonic formulation can be easily recovered by a slight modification of the procedure. This provides a mechanism for deriving both the field evolution equations and the gauge conditions from the action principle, with a view on using symplectic integrators for a constraint-preserving numerical evolution. The gauge sources corresponding to the 'puncture gauge' conditions are identified in this context.
Cite
@article{arxiv.1008.0747,
title = {Action principle for Numerical Relativity evolution systems},
author = {C. Bona and C. Bona-Casas and C. Palenzuela},
journal= {arXiv preprint arXiv:1008.0747},
year = {2010}
}
Comments
Revised version, includes explicit expresions for gauge sources corresponding to 1+log and gamma-driver gauge conditions ('punctures' gauge)