Free and constrained symplectic integrators for numerical general relativity

DOI10.1088/0264-9381/25/22/225018zbMATH Open1152.83309arXiv0807.0734OpenAlexW3098051418WikidataQ125950460 ScholiaQ125950460MaRDI QIDQ3544867FDOQ3544867

Authors: Ronny Richter, Christian Lubich 

Publication date: 8 December 2008

Published in: Classical and Quantum Gravity (Search for Journal in Brave)

Abstract: We consider symplectic time integrators in numerical General Relativity and discuss both free and constrained evolution schemes. For free evolution of ADM-like equations we propose the use of the Stoermer-Verlet method, a standard symplectic integrator which here is explicit in the computationally expensive curvature terms. For the constrained evolution we give a formulation of the evolution equations that enforces the momentum constraints in a holonomically constrained Hamiltonian system and turns the Hamilton constraint function from a weak to a strong invariant of the system. This formulation permits the use of the constraint-preserving symplectic RATTLE integrator, a constrained version of the Stoermer-Verlet method. The behavior of the methods is illustrated on two effectively 1+1-dimensional versions of Einstein's equations, that allow to investigate a perturbed Minkowski problem and the Schwarzschild space-time. We compare symplectic and non-symplectic integrators for free evolution, showing very different numerical behavior for nearly-conserved quantities in the perturbed Minkowski problem. Further we compare free and constrained evolution, demonstrating in our examples that enforcing the momentum constraints can turn an unstable free evolution into a stable constrained evolution. This is demonstrated in the stabilization of a perturbed Minkowski problem with Dirac gauge, and in the suppression of the propagation of boundary instabilities into the interior of the domain in Schwarzschild space-time.

Full work available at URL: https://arxiv.org/abs/0807.0734

Recommendations

• The applicability of constrained symplectic integrators in general relativity
• Strongly hyperbolic Hamiltonian systems in numerical relativity: formulation and symplectic integration
• Forward symplectic integrators for solving gravitational few-body problems
• Symbolic and numerical analysis in general relativity with open source computer algebra systems
• NUMERICAL RELATIVITY
• Numerical relativity
• Numerical relativity
• scientific article; zbMATH DE number 3888440
• Numerical relativity using a generalized harmonic decomposition
• A connection approach to numerical relativity

Mathematics Subject Classification ID

Einstein's equations (general structure, canonical formalism, Cauchy problems) (83C05) Software, source code, etc. for problems pertaining to relativity and gravitational theory (83-04)

Cited In (5)

• Strongly hyperbolic Hamiltonian systems in numerical relativity: formulation and symplectic integration
• The applicability of constrained symplectic integrators in general relativity
• Continuum and discrete initial-boundary value problems and Einstein's field equations
• Constructing of constraint preserving scheme for Einstein equations
• On the linearization of Regge calculus

This page was built for publication: Free and constrained symplectic integrators for numerical general relativity

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q3544867)