KHARMA: Flexible, Portable Performance for GRMHD

Abstract

KHARMA (an acronym for "Kokkos-based High-Accuracy Relativistic Magnetohydrodynamics with Adaptive mesh refinement") is a new open-source code for conducting general-relativistic magnetohydrodynamic simulations in stationary spacetimes, primarily of accretion systems. It implements among other options the High-Accuracy Relativistic Magnetohydrodynamics (HARM) scheme, but is written from scratch in C++ with the Kokkos programming model in order to run efficiently on both CPUs and GPUs. In addition to being fast, KHARMA is written to be readable, modular, and extensible, separating functionality into "packages," representing, e.g., algorithmic components or physics extensions. Components of the core ideal GRMHD scheme can be swapped at runtime, and additional packages are included to simulate electron temperature evolution, viscous hydrodynamics, and for designing chained multi-scale "bridged" simulations. This chapter presents the computational environment and requirements for KHARMA, features and design which meet these requirements, and finally, validation and performance data.

Chapter submitted to "New Frontiers in GRMHD Simulations of Accreting Black Holes," Springer Nature Singapore, ed. Cosimo Bambi, Yosuke Mizuno, Swarnim Shashank and Feng Yuan. Document released under LA-UR-23-23232

References in corpus (8)

• First M87 Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole
• First M87 Event Horizon Telescope Results. VI. The Shadow and Mass of the Central Black Hole
• First M87 Event Horizon Telescope Results. IV. Imaging the Central Supermassive Black Hole
• First Sagittarius A* Event Horizon Telescope Results. V. Testing Astrophysical Models of the Galactic Center Black Hole
• The Disc-Jet Symbiosis Emerges: Modeling the Emission of Sagittarius A* with Electron Thermodynamics
• Discriminating Accretion States via Rotational Symmetry in Simulated Polarimetric Images of M87
• An Extension of the Athena++ Code Framework for Radiation-Magnetohydrodynamics in General Relativity Using a Finite-Solid-Angle Discretization
• grim: A Flexible, Conservative Scheme for Relativistic Fluid Theories