HERACLES++: A multidimensional Eulerian code for exascale computing

Numerical simulations of multidimensional astrophysical fluids present considerable challenges. However, the development of exascale computing has significantly enhanced computational capabilities, motivating the development of new codes that can take full advantage of these resources. In this article, we introduce HERACLES++ , a new hydrodynamics code with high portability, optimized for exascale machines with different architectures and running efficiently both on CPUs and GPUs. The code is Eulerian and employs a Godunov finite-volume method to solve the hydrodynamics equations, which ensures accuracy in capturing shocks and discontinuities. It includes different Riemann solvers, equations of state, and gravity solvers. It works in Cartesian and spherical coordinates, either in 1D, 2D, or 3D, and uses passive scalars to handle gases with several species. The code accepts a user-supplied heating or cooling term to treat a variety of astrophysical contexts. In addition to the usual series of benchmarking tests, we use HERACLES++ to simulate the propagation of a supernova shock in a red-supergiant star envelope, from minutes after core collapse until shock emergence. In 1D, the results from HERACLES++ are in agreement with those of V1D for the same configuration. In 3D, the Rayleigh-Taylor instability develops and modifies the 1D picture by introducing density and composition fluctuations as well as turbulence. The focus on a wedge, rather than the full solid angle, and the ability to run HERACLES++ with a large number of GPUs allow for long-term simulations of 3D supernova ejecta with a sub-degree resolution. Future developments will extend HERACLES++ to a radiation-hydrodynamics code.