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.

CheasePy

FLIPEC, an Ideal MHD free-boundary axisymmetric equilibrium solver in the presence of macroscopic flows.

FEQIS: A free-boundary equilibrium solver for integrated modeling of tokamak plasmas

FreeGSNKE: A Python-based dynamic free-boundary toroidal plasma equilibrium solver

A Multidomain Chebyshev Pseudo-Spectral Method for Fluid Flow and Heat Transfer from Square Cylinders

pyEQUIB Python Package, an addendum to proEQUIB: IDL Library for Plasma Diagnostics and Abundance Analysis

SfePy - Write Your Own FE Application

A new hybrid code (CHIEF) implementing the inertial electron fluid equation without approximation

ExaGRyPE: Numerical General Relativity Solvers Based upon the Hyperbolic PDEs Solver Engine ExaHyPE

celmech: A Python Package for Celestial Mechanics

The muphyII Code: Multiphysics Plasma Simulation on Large HPC Systems

Self-consistent Modeling of Inductively Coupled Plasma Discharges

JefiPIC: A 3-D Full Electromagnetic Particle-in-Cell Simulator Based on Jefimenko's Equations on GPU

AHKASH : a new Hybrid particle-in-cell code for simulations of astrophysical collisionless plasma

Geometric Particle-In-Cell discretizations of a plasma hybrid model with kinetic ions and mass-less fluid electrons

One-Dimensional, One-Phase and Two-Phase Eulerian Explicit Shock Tube Simulation Code

Smilei : A collaborative, open-source, multi-purpose particle-in-cell code for plasma simulation

ExaHyPE: An Engine for Parallel Dynamically Adaptive Simulations of Wave Problems

AHKASH: a new Hybrid particle-in-cell code for simulations of astrophysical collisionless plasma

Pegasus: A New Hybrid-Kinetic Particle-in-cell Code for Astrophysical Plasma Dynamics.

The magnetohydrodynamic-particle-in-cell module in athena plus plus : implementation and code tests