Solving the Einstein Constraints Numerically on Compact Three-Manifolds Using Hyperbolic Relaxation

Abstract

The effectiveness of the hyperbolic relaxation method for solving the Einstein constraint equations numerically is studied here on a variety of compact orientable three-manifolds. Convergent numerical solutions are found using this method on manifolds admitting negative Ricci scalar curvature metrics, i.e. those from the $H^3$ and the $H^2\times S^1$ geometrization classes. The method fails to produce solutions, however, on all the manifolds examined here admitting non-negative Ricci scalar curvatures, i.e. those from the $S^3$, $S^2\times S^1$, and the $E^3$ classes. This study also finds that the accuracy of the convergent solutions produced by hyperbolic relaxation can be increased significantly by performing fairly low-cost standard elliptic solves using the hyperbolic relaxation solutions as initial guesses.

arXiv comment7 pages, 8 figures, to appear in PRD

References in corpus (5)

• Optimal Constraint Projection for Hyperbolic Evolution Systems
• Hyperbolic Relaxation Method for Elliptic Equations
• Solving Partial Differential Equations Numerically on Manifolds with Arbitrary Spatial Topologies
• NRPyElliptic: A Fast Hyperbolic Relaxation Elliptic Solver for Numerical Relativity, I: Conformally Flat, Binary Puncture Initial Data
• Simple Numerical Solutions to the Einstein Constraints on Various Three-Manifolds