An official website of the United States government
Abstract Next-generation gravitational wave detectors such as Cosmic Explorer, the Einstein Telescope, and LISA, demand highly accurate and extensive gravitational wave (GW) catalogs to faithfully extract physical parameters from observed signals. However, numerical relativity (NR) faces significant challenges in generating these catalogs at the required scale and accuracy on modern computers, as NR codes do not fully exploit modern GPU capabilities. In response, we extend NRPy, a Python-based NR code-generation framework, to develop NRPyEllipticGPU—a CUDA-optimized elliptic solver tailored for the binary black hole (BBH) initial data problem. NRPyEllipticGPU is the first GPU-enabled elliptic solver in the NR community, supporting a variety of coordinate systems and demonstrating substantial performance improvements on both consumer-grade and HPC-grade GPUs. We show that, when compared to a high-end CPU, NRPyEllipticGPU achieves on a high- end GPU up to a sixteenfold speedup in single precision while increasing double- precision performance by a factor of 2–4. This performance boost leverages the GPU’s superior parallelism and memory bandwidth to achieve a compute-bound application and enhancing the overall simulation efficiency. As NRPyEllipticGPU shares the core infrastructure common to NR codes, this work serves as a practical guide for developing full, CUDA-optimized NR codes.
more »
« less
NRPyElliptic: A Fast Hyperbolic Relaxation Elliptic Solver for Numerical Relativity, I: Conformally Flat, Binary Puncture Initial Data
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 ∼12x slower than TwoPunctures at setting up full-3D BBH ID, NRPyElliptic requires only ≈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.
more »
« less
- Award ID(s):
- 1726534
- PAR ID:
- 10315154
- Date Published:
- Journal Name:
- ArXivorg
- ISSN:
- 2331-8422
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract There are few observed high-mass X-ray binaries (HMXBs) that harbor massive black holes (BHs), and none are likely to result in a binary black hole (BBH) that merges within a Hubble time; however, we know that massive merging BBHs exist from gravitational-wave (GW) observations. We investigate the role that X-ray and GW observational selection effects play in determining the properties of their respective detected binary populations. We find that, as a result of selection effects, detectable HMXBs and detectable BBHs form at different redshifts and metallicities, with detectable HMXBs forming at much lower redshifts and higher metallicities than detectable BBHs. We also find disparities in the mass distributions of these populations, with detectable merging BBH progenitors pulling to higher component masses relative to the full detectable HMXB population. Fewer than 3% of detectable HMXBs host BHs >35M⊙in our simulated populations. Furthermore, we find the probability that a detectable HMXB will merge as a BBH system within a Hubble time is ≃0.6%. Thus, it is unsurprising that no currently observed HMXB is predicted to form a merging BBH with high probability.more » « less
-
Abstract Fixed-point fast sweeping WENO methods are a class of efficient high-order numerical methods to solve steady-state solutions of hyperbolic partial differential equations (PDEs). The Gauss-Seidel iterations and alternating sweeping strategy are used to cover characteristics of hyperbolic PDEs in each sweeping order to achieve fast convergence rate to steady-state solutions. A nice property of fixed-point fast sweeping WENO methods which distinguishes them from other fast sweeping methods is that they are explicit and do not require inverse operation of nonlinear local systems. Hence, they are easy to be applied to a general hyperbolic system. To deal with the difficulties associated with numerical boundary treatment when high-order finite difference methods on a Cartesian mesh are used to solve hyperbolic PDEs on complex domains, inverse Lax-Wendroff (ILW) procedures were developed as a very effective approach in the literature. In this paper, we combine a fifth-order fixed-point fast sweeping WENO method with an ILW procedure to solve steady-state solution of hyperbolic conservation laws on complex computing regions. Numerical experiments are performed to test the method in solving various problems including the cases with the physical boundary not aligned with the grids. Numerical results show high-order accuracy and good performance of the method. Furthermore, the method is compared with the popular third-order total variation diminishing Runge-Kutta (TVD-RK3) time-marching method for steady-state computations. Numerical examples show that for most of examples, the fixed-point fast sweeping method saves more than half CPU time costs than TVD-RK3 to converge to steady-state solutions.more » « less
-
null (Ed.)We propose and study a one-dimensional \begin{document}$$ 2\times 2 $$\end{document} hyperbolic Eulerian system with local relaxation from critical threshold phenomena perspective. The system features dynamic transition between strictly and weakly hyperbolic. For different classes of relaxation we identify intrinsic critical thresholds for initial data that distinguish global regularity and finite time blowup. For relaxation independent of density, we estimate bounds on density in terms of velocity where the system is strictly hyperbolic.more » « less
-
Abstract It has been proposed that some black holes (BHs) in binary black hole (BBH) systems are born from “hierarchical mergers” (HMs), i.e., earlier mergers of smaller BHs. These HM products have spin magnitudes χ ∼ 0.7, and, if they are dynamically assembled into BBH systems, their spin orientations will sometimes be antialigned with the binary orbital angular momentum. In fact, as Baibhav et al. showed, ∼16% of BBH systems that include HM products will have an effective inspiral spin parameter, χ eff < −0.3. Nevertheless, the LIGO–Virgo–KAGRA (LVK) gravitational-wave (GW) detectors have yet to observe a BBH system with χ eff ≲ −0.2, leading to upper limits on the fraction of HM products in the population. We fit the astrophysical mass and spin distribution of BBH systems and measure the fraction of BBH systems with χ eff < −0.3, which implies an upper limit on the HM fraction. We find that fewer than 26% of systems in the underlying BBH population include HM products (90% credibility). Even among BBH systems with primary masses m 1 = 60 M ⊙ , the HM fraction is less than 69%, which may constrain the location of the pair-instability mass gap. With 300 GW events (to be expected in the LVK’s next observing run), if we fail to observe a BBH with χ eff < −0.3, we can conclude that the HM fraction is smaller than 2.5 − 2.2 + 9.1 % .more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- The DOI is not currently available.
