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.

arXiv commentChapter 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 (22)

• 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
• HARM: A Numerical Scheme for General Relativistic Magnetohydrodynamics
• First M87 Event Horizon Telescope Results. VIII. Magnetic Field Structure near The Event Horizon
• A Measurement of the Electromagnetic Luminosity of a Kerr Black Hole
• First Sagittarius A* Event Horizon Telescope Results. V. Testing Astrophysical Models of the Galactic Center Black Hole
• Universal Interferometric Signatures of a Black Hole's Photon Ring
• The Event Horizon General Relativistic Magnetohydrodynamic Code Comparison Project
• Primitive Variable Solvers for Conservative General Relativistic Magnetohydrodynamics
• The Disc-Jet Symbiosis Emerges: Modeling the Emission of Sagittarius A* with Electron Thermodynamics
• A Second Order Godunov Method for Multidimensional Relativistic Magnetohydrodynamics
• Discriminating Accretion States via Rotational Symmetry in Simulated Polarimetric Images of M87
• An Extended Magnetohydrodynamics Model for Relativistic Weakly Collisional Plasmas
• Evolution of Accretion Discs around a Kerr Black Hole using Extended Magnetohydrodynamics
• Robust Recovery of Primitive Variables in Relativistic Ideal Magnetohydrodynamics
• Iharm3D: Vectorized General Relativistic Magnetohydrodynamics
• Parthenon -- a performance portable block-structured adaptive mesh refinement framework
• Black Hole Magnetic Fields and Their Imprint on Circular Polarization Images
• An Extension of the Athena++ Code Framework for Radiation-Magnetohydrodynamics in General Relativity Using a Finite-Solid-Angle Discretization
• Unraveling Twisty Linear Polarization Morphologies in Black Hole Images
• grim: A Flexible, Conservative Scheme for Relativistic Fluid Theories