A Localized Current-Sheet Magnetic-Diffusion and Heating Prescription for Ideal-GRMHD Simulations of M87*-like Accretion Flows
Abstract
Ideal-GRMHD calculations provide the bulk accretion flows used in black-hole emission models. A localized dissipation prescription can add a resolved map of current-sheet activity while retaining that dynamical reference. We implement such a prescription in Athena++, using a smooth threshold on the dimensionless current proxy αJ to select the magnetic-diffusion coefficient. The constrained-transport update is coupled to a conservative energy-flux correction. Gas heating follows from primitive recovery; qOhm=ηJ2 is recorded as a diagnostic, without a second energy source. Static Fourier and traveling Alfvén tests measure the evolved thermal gain, and an open-boundary Schwarzschild test closes the Killing-energy ledger. The relative total-energy residual is below 2×10−15 in these tests. A two-resolution Harris matrix gives localized rate exponents 0.510±0.027 and 0.488±0.028, consistent with 1/2 over the specified early interval. The η=10−3 rate exceeds the same-resolution ideal control by a factor of 5.96. In matched axisymmetric M87*-like runs, localized diffusion changes the early mean magnetic flux, accretion rate, and normalized flux by less than 0.02%, with a radial-density distance of 0.005%. During the late active phase, a common-restart on/off comparison keeps these mean shifts below 1.6%, within the intrinsic temporal variability, and gives a density distance of 0.29%. Uniform diffusion at the tested coefficient reduces the raw flux and accretion rate by 85% and 99.8%. These comparisons establish a spatially selective diffusion prescription with explicit energy accounting and a quantified, small bulk-flow response.
Related
- 3D simulations of magnetospheric accretion in T Tauri stars I. Disk truncation, stellar torques, and application to observations
- Three-dimensional simulations of accretion disks in pre-CE systems -- II. Accretion efficiency and angular momentum transport
- The long-term evolution of a black hole accretion torus: a comparison between GRMHD simulations and relativistic thin disk theory
- Scaled three-dimensional simulations of giant-planet accretion and disk clearing
- 3D simulations of magnetospheric accretion in T Tauri stars: I. Disk truncation, stellar torques, and application to observations
BibTeX
@misc{dai2026localized,
title = {{A Localized Current-Sheet Magnetic-Diffusion and Heating Prescription for Ideal-GRMHD Simulations of M87*-like Accretion Flows}},
author = {Songyan Dai and Renyue Cen},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.12377},
eprint = {2609.12377},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.12377},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - A Localized Current-Sheet Magnetic-Diffusion and Heating Prescription for Ideal-GRMHD Simulations of M87*-like Accretion Flows
AU - Songyan Dai
AU - Renyue Cen
PY - 2026
DA - 2026/09/11
PB - arXiv
DO - 10.48550/arxiv.2609.12377
UR - https://arxiv.org/abs/2609.12377
AB - Ideal-GRMHD calculations provide the bulk accretion flows used in black-hole emission models. A localized dissipation prescription can add a resolved map of current-sheet activity while retaining that dynamical reference. We implement such a prescription in Athena++, using a smooth threshold on the dimensionless current proxy α_J to select the magnetic-diffusion coefficient. The constrained-transport update is coupled to a conservative energy-flux correction. Gas heating follows from primitive recovery; q_Ohm=ηJ^2 is recorded as a diagnostic, without a second energy source. Static Fourier and traveling Alfvén tests measure the evolved thermal gain, and an open-boundary Schwarzschild test closes the Killing-energy ledger. The relative total-energy residual is below 2×10^−15 in these tests. A two-resolution Harris matrix gives localized rate exponents 0.510±0.027 and 0.488±0.028, consistent with 1/2 over the specified early interval. The η=10^−3 rate exceeds the same-resolution ideal control by a factor of 5.96. In matched axisymmetric M87*-like runs, localized diffusion changes the early mean magnetic flux, accretion rate, and normalized flux by less than 0.02%, with a radial-density distance of 0.005%. During the late active phase, a common-restart on/off comparison keeps these mean shifts below 1.6%, within the intrinsic temporal variability, and gives a density distance of 0.29%. Uniform diffusion at the tested coefficient reduces the raw flux and accretion rate by 85% and 99.8%. These comparisons establish a spatially selective diffusion prescription with explicit energy accounting and a quantified, small bulk-flow response.
ER -