REVIEW 24 cited by
AthenaK: A Performance-Portable Version of the Athena++ AMR Framework
Not yet reviewed by Pith; the record is open.
This paper has not been read by Pith yet. Machine review is queued; the pith claim, tier, and objections will appear here once it completes.
SPECIMEN: schema-true, not a live event
T0 review · schema-true
One-sentence machine reading of the paper's core claim.
pith:XXXXXXXX · record.json · timestamp
Signed reviews
read the original abstract
We describe AthenaK: a new implementation of the Athena++ block-based adaptive mesh refinement (AMR) framework using the Kokkos programming model. Finite volume methods for Newtonian, special relativistic (SR), and general relativistic (GR) hydrodynamics and magnetohydrodynamics (MHD), and GR-radiation hydrodynamics and MHD, as well as a module for evolving Lagrangian tracer or charged test particles (e.g., cosmic rays) are implemented using the framework. In two companion papers we describe (1) a new solver for the Einstein equations based on the Z4c formalism and (2) a GRMHD solver in dynamical spacetimes also implemented using the framework, enabling new applications in numerical relativity. By adopting Kokkos, the code can be run on virtually any hardware, including CPUs, GPUs from multiple vendors, and emerging ARM processors. AthenaK shows excellent performance and weak scaling, achieving over one billion cell updates per second for hydrodynamics in three-dimensions on a single NVIDIA Grace Hopper processor and with a typical parallel efficiency of 80% on 65536 AMD GPUs on the OLCF Frontier system. Such performance portability enables AthenaK to leverage modern exascale computing systems for challenging applications in astrophysical fluid dynamics, numerical relativity, and multimessenger astrophysics.
Forward citations
Cited by 24 Pith papers
-
Galactic Dust Polarization in Turbulent Multiphase ISM: On the Origin of the $EE/BB$ Asymmetry
Dust polarization from the unstable neutral medium best matches Planck's EE/BB asymmetry, suggesting it dominates the high-latitude polarized dust foreground.
-
The Interplay of Parametric and Magnetorotational Instabilities in Oscillatory Shear Flows
In forced oscillatory shear-flow simulations, MRI turbulence quenches the parametric instability for weak forcing, while for stronger forcing the parametric instability emerges and dominates vertical angular-momentum ...
-
Radiation GRMHD Models of Accretion onto Stellar-Mass Black Holes: I. Survey of Eddington Ratios
Full-transport radiation GRMHD simulations show super-Eddington black hole accretion is geometrically thick, drives strong outflows, and radiates with very low efficiency (below about 0.5% at 150 times Eddington).
-
Relativistic gas accretion onto supermassive black Hole binaries from inspiral through merger
A new code hand-off enables a 3D GRMHD simulation of an equal-mass nonspinning supermassive black hole binary from 20M separation through merger and postmerger, showing minidisk dissolution, declining accretion, and a...
-
Plasmoid-Mediated 2D Magnetic Reconnection in Partially Ionized Plasmas
In weakly ionized plasma, neutral-ion decoupling suppresses large-scale coalescence, producing a chain of small plasmoids with extreme ion density pile-ups while the reconnection rate stays slow.
-
GRACE: An Open-Source Framework for GPU-Accelerated Numerical Relativity
GRACE is a validated, open-source, Kokkos+p4est GPU-portable framework that evolves ideal GRMHD with constrained transport self-consistently coupled to Z4c Einstein equations on fixed or adaptive meshes.
-
Multiphase gas in Circumgalactic cloud complexes: Insights from kiloparsec-scale Magnetohydrodynamic Turbulence Simulations
Cold-gas survival in the circumgalactic medium is set by t_cool/t_mix; small-scale cold gas requires dense, quiescent regions, not turbulent cascade patches.
-
$\textit{BMAD}$-Circumbinary Magnetically Arrested Disks around Stellar or Black Hole Binaries: Hot Accretion Flows, Disk Properties, and Angular Momentum Transfer
Circumbinary accretion disks can enter a magnetically arrested state, and in weakly cooled or adiabatic regimes the resulting magnetic flux eruptions may drive the binary orbit to shrink.
-
Self-Similar Cosmic-Ray Transport in High-Resolution Magnetohydrodynamic Turbulence
In a 10,240^3 MHD turbulence simulation, cosmic-ray diffusion is nearly energy-independent over two decades in energy because rare, large excursions set the transport rate.
-
Turbulence in Magnetised Neutron Stars
Magnetic-field instabilities in a neutron star drive a Kolmogorov-like turbulent cascade and a 145 Hz oscillation in the highest-resolution runs to date, while net magnetic helicity remains near zero.
-
A "Neural" Riemann solver for Relativistic Hydrodynamics
Small neural networks trained on exact solutions can replace the iterative root-finding in a relativistic Riemann solver, giving exact-like accuracy about 14 times faster in 1D tests.
-
The plunging region of a thin accretion disc around a Schwarzschild black hole
A 3D plasma simulation of a thin disc confirms the analytic geodesic model for the plunging region inside the ISCO, finding a 5.3% angular momentum loss and rejecting constant alpha-disc models there.
-
Idealized Global Models of Accretion Disks with Strong Toroidal Magnetic Fields
A zero-net-vertical-flux accretion disk with a strong toroidal magnetic field settles into a state with midplane magnetic and gas pressures comparable, and accretion stress alpha about 0.1, when the gas scale height i...
-
Cyclic Zoom: Multiscale GRMHD Modeling of Black Hole Accretion and Feedback
Cyclic zoom simulations bridge horizon-to-galaxy scales in GRMHD, yielding a Bondi accretion suppression near (10 rg/rB)^1/2 and feedback efficiencies of 1% (a=0) and 10% (a~0.9).
-
The Type Ia Supernova and Asymptotic Giant Branch Stellar Ejecta-regulated Interstellar Medium of Massive Galaxies
Resolved supernova heating and stellar mass loss in 3D simulations naturally produce cool-core-like and non-cool-core-like interstellar media whose properties are set by the confining circumgalactic pressure.
-
Evolution of Supernova Remnants in a Cloudy Multiphase Interstellar Medium
A supernova remnant in a cloudy interstellar medium expands at the same rate as in a uniform medium but loses much more energy to cooling at cloud interfaces, reducing its hot gas mass and momentum.
-
Effects of Differential Rotation on the Maximum Mass of Neutron Stars
Differential rotation substantially boosts hyperonic neutron-star maximum mass to the GW190814 range, yet a soft SU(6) EoS still fails the PSR J0740+6620 constraint at 346 Hz.
-
Magnetic field dynamics in isolated neutron stars with an external dipole field
Long-term numerical relativity simulations find that neutron star magnetic fields relax to stable mixed configurations with toroidal energy fraction ≲10% within one Alfvén time after Tayler instability saturation.
-
SYCL for Energy-Efficient Numerical Astrophysics: the case of DPEcho
Using the DPEcho SYCL GR-MHD benchmark, the authors measure energy efficiency across 12 CPU and GPU systems and report that GPUs deliver far more cell updates per joule than CPUs.
-
AthenaK simulations of the binary black hole merger GW150914
A new open-source GPU code, AthenaK, reproduces the GW150914 merger: remnant mass within 0.01%, spin within 0.02%, and waveform phase within about 0.35 radians of established simulations.
-
Origin of the Multi-Phase Interstellar Medium: the Effects of Turbulence and Magnetic Field
Turbulence and magnetic fields together sustain a substantial fraction of gas at intermediate temperatures in the simulated multiphase interstellar medium.
-
Characterizing 3D Magnetic Fields and Turbulence in H I Clouds
A neural network maps H I spectral cubes to 3D magnetic field orientation, strength, sonic and Alfven Mach numbers, and it is applied to two FAST-survey clouds in Monoceros.
-
Not-Quite-Transcendental Functions For Logarithmic Interpolation of Tabulated Data
Second-order NQT functions, a C^1 piecewise quadratic approximation to log2, preserve second-order convergence of linear interpolation on logarithmic grids while running faster than true logarithms.
-
The Kratos Framework for Heterogeneous Astrophysical Simulations: Fundamental Infrastructures and Hydrodynamics
Kratos is a new GPU-optimized astrophysical hydrodynamics framework using mixed precision and stream-aware communication, verified on standard shock and instability benchmarks.
Discussion (0). Continue with ORCID to comment.