A Localized Current-sheet Magnetic-diffusion and Heating Prescription for Ideal-GRMHD Simulations of M87*-like Accretion Flows
- 1Zhejiang University
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; q Ohm = ηJ ² 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 ⁻¹⁵ 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 ⁻³ 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
@article{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},
journal = {The Astrophysical Journal},
year = {2026},
month = nov,
doi = {10.3847/1538-4357/aea6c7},
url = {https://iopscience.iop.org/article/10.3847/1538-4357/aea6c7},
note = {CC BY 4.0},
}
RIS
TY - JOUR
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/11/02
JO - The Astrophysical Journal
DO - 10.3847/1538-4357/aea6c7
UR - https://iopscience.iop.org/article/10.3847/1538-4357/aea6c7
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 ² 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 ⁻¹⁵ 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 ⁻³ 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.
SN - 1538-4357
SN - 0004-637X
ER -