Free and constrained symplectic integrators for numerical general relativity

DOI10.1088/0264-9381/25/22/225018MaRDI QIDQ3544867zbMATH OpenOpenAlexWikidataFDO

Authors Ronny Richter, Christian Lubich

Publication date 8 December 2008

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

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

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)

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.

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

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)