ExaGRyPE: Numerical General Relativity Solvers Based upon the Hyperbolic PDEs Solver Engine ExaHyPE
arXiv:2406.11626 · doi:10.1016/j.cpc.2024.109435
Abstract
ExaGRyPE describes a suite of solvers for numerical relativity, based upon ExaHyPE 2, the second generation of our Exascale Hyperbolic PDE Engine. The presented generation of ExaGRyPE solves the Einstein equations in the CCZ4 formulation under a 3+1 foliation and focuses on black hole spacetimes. It employs a block-structured Cartesian grid carrying a higher-order order Finite Difference scheme with adaptive mesh refinement, it facilitates massive parallelism combining message passing, domain decomposition and task parallelism, and it supports the injection of particles into the grid as data probes or tracers. We introduce the ExaGRyPE-specific building blocks within ExaHyPE 2, and discuss its software architecture and compute-n-feel: For this, we formalize the creation of any specific simulation with ExaGRyPE as a sequence of lowering operations, where abstract logical tasks are successively broken into finer tasks until we obtain an abstraction level that can be mapped onto a C++ executable. The overall program logic is fully specified via a domain-specific Python interface, we map this logic onto a more detailed set of numerical tasks, subsequently lower this representation onto technical tasks that the underlying ExaHyPE engine uses to parallelize the application, before eventually the technical tasks are mapped onto task graphs including the actual PDE term evaluations, initial conditions, boundary conditions, and so forth. We end up with a rigorous separation of concerns which shields ExaGRyPE users from technical details and hence simplifies the development of novel physical models. We present the simulations and data for the gauge wave, static single black holes and rotating binary black hole systems, demonstrating that the code base is mature and usable. However, we also uncover domain-specific numerical challenges that need further study by the community in future work.
References in corpus (33)
- The Dynamics of General Relativity
- The Einstein Toolkit: A Community Computational Infrastructure for Relativistic Astrophysics
- Gauge conditions for long-term numerical black hole evolutions without excision
- Calibration of Moving Puncture Simulations
- A single-domain spectral method for black hole puncture data
- Constraint damping in the Z4 formulation and harmonic gauge
- Conformal and covariant formulation of the Z4 system with constraint-violation damping
- Numerical simulation of orbiting black holes
- Simulating coalescing compact binaries by a new code SACRA
- Compact binary evolutions with the Z4c formulation
- Accurate black hole evolutions by fourth-order numerical relativity
- General-covariant evolution formalism for Numerical Relativity
- A Survey of High Level Frameworks in Block-Structured Adaptive Mesh Refinement Packages
- Geometry and Regularity of Moving Punctures
- GRChombo : Numerical Relativity with Adaptive Mesh Refinement
- Toward standard testbeds for numerical relativity
- Distributing many points on spheres: minimal energy and designs
- A symmetry-breaking mechanism for the Z4 general-covariant evolution system
- An explicit harmonic code for black-hole evolution using excision
- ExaHyPE: An Engine for Parallel Dynamically Adaptive Simulations of Wave Problems
- Quasi-equilibrium binary black hole sequences for puncture data derived from helical Killing vector conditions
- Implementation of standard testbeds for numerical relativity
- Conformal and covariant Z4 formulation of the Einstein equations: strongly hyperbolic first-order reduction and solution with discontinuous Galerkin schemes
- Massively Parallel Simulations of Binary Black Hole Intermediate-Mass-Ratio Inspirals
- GRAthena++: puncture evolutions on vertex-centered oct-tree AMR
- On GLM curl cleaning for a first order reduction of the CCZ4 formulation of the Einstein field equations
- Lessons for adaptive mesh refinement in numerical relativity
- The Peano software - parallel, automaton-based, dynamically adaptive grid traversals
- Reducing reflections from mesh refinement interfaces in numerical relativity
- Enclave Tasking for Discontinuous Galerkin Methods on Dynamically Adaptive Meshes
- Adaptive mesh refinement in binary black holes simulations
- Spherical accretion of collisional gas in modified gravity I: self-similar solutions and a new cosmological hydrodynamical code
- SYCL compute kernels for ExaHyPE