NRPyElliptic: A Fast Hyperbolic Relaxation Elliptic Solver for Numerical Relativity, I: Conformally Flat, Binary Puncture Initial Data

Abstract

We introduce NRPyElliptic, an elliptic solver for numerical relativity (NR) built within the NRPy+ framework. As its first application, NRPyElliptic sets up conformally flat, binary black hole (BBH) puncture initial data (ID) on a single numerical domain, similar to the widely used TwoPunctures code. Unlike TwoPunctures, NRPyElliptic employs a hyperbolic relaxation scheme, whereby arbitrary elliptic PDEs are trivially transformed into a hyperbolic system of PDEs. As consumers of NR ID generally already possess expertise in solving hyperbolic PDEs, they will generally find NRPyElliptic easier to tweak and extend than other NR elliptic solvers. When evolved forward in (pseudo)time, the hyperbolic system exponentially reaches a steady state that solves the elliptic PDEs. Notably NRPyElliptic accelerates the relaxation waves, which makes it many orders of magnitude faster than the usual constant-wavespeed approach. While it is still $\sim 12$x slower than TwoPunctures at setting up full-3D BBH ID, NRPyElliptic requires only $\approx 0.3\%$ of the runtime for a full BBH simulation in the Einstein Toolkit. Future work will focus on improving performance and generating other types of ID, such as binary neutron star.

15 pages, 10 figures. Matches published version

References in corpus (11)

• Initial data for black hole-neutron star binaries: a flexible, high-accuracy spectral method
• Time Step Size Limitation Introduced by the BSSN Gamma Driver
• Conformal Thin-Sandwich Solver for Generic Initial Data
• Constructing Binary Neutron Star Initial Data with High Spins, High Compactness, and High Mass-Ratios
• Initial data for binary neutron stars with adjustable eccentricity
• Construction of initial data for 3+1 numerical relativity
• The initial value problem as it relates to numerical relativity
• Fully covariant and conformal formulation of the Z4 system in a reference-metric approach: comparison with the BSSN formulation in spherical symmetry
• Numerical method for binary black hole/neutron star initial data: Code test
• A scalable elliptic solver with task-based parallelism for the SpECTRE numerical relativity code
• Initial data for black hole-neutron star binaries, with rotating stars

Cited by in corpus (7)

• The Science of the Einstein Telescope
• CTTK: A new method to solve the initial data constraints in numerical relativity
• Improved Moving-Puncture Techniques for Compact Binary Simulations
• BHaHAHA: A Fast, Robust Apparent Horizon Finder Library for Numerical Relativity
• superB/NRPy: Scalable, Task-Based Numerical Relativity for 3G Gravitational Wave Science
• Solving the Einstein Constraints Numerically on Compact Three-Manifolds Using Hyperbolic Relaxation
• The Parallel Compact Object CALculator: An Efficient General Relativistic Initial Data Solver for Compact Objects