TY - JOUR
T1 - Transition region from turbulent to dead zone in protoplanetary disks
T2 - Local shearing box simulations
AU - Pucci, Fulvia
AU - Tomida, Kengo
AU - Stone, James
AU - Takasao, Shinsuke
AU - Ji, Hantao
AU - Okamura, Shoichi
N1 - Publisher Copyright:
© 2021 Institute of Physics Publishing. All rights reserved.
PY - 2021/1/20
Y1 - 2021/1/20
N2 - The dynamical evolution of protoplanetary disks is of key interest for building a comprehensive theory of planet formation and to explain the observational properties of these objects. Using the magnetohydrodynamics code Athena++, with an isothermal shearing box setup, we study the boundary between the active and dead zone, where the accretion rate changes and mass can accumulate. We quantify how the turbulence level is affected by the presence of a non-uniform Ohmic resistivity in the radial x direction that leads to a region of inhibited turbulence (or dead zone). Comparing the turbulent activity to that of ideal simulations, the turbulence-inhibited area shows density fluctuations and magnetic activity at its boundaries, driven by energy injection from the active (ideal) zone boundaries. We find magnetic dissipation to be significantly stronger in the ideal regions, and the turbulence penetration through the boundary of the dead zone is determined by the value of the resistivity itself, through the Ohmic dissipation process, though the thickness of the transition does not play a significant role in changing the dissipation. We investigate the 1D spectra along the shearing direction: magnetic spectra appear flat at large scales both in ideal as well as resistive simulations, though a Kolmogorov scaling over more than one decade persists in the dead zone, suggesting the turbulent cascade is determined by the hydrodynamics of the system: magnetorotational instability dynamo action is inhibited where sufficiently high resistivity is present.
AB - The dynamical evolution of protoplanetary disks is of key interest for building a comprehensive theory of planet formation and to explain the observational properties of these objects. Using the magnetohydrodynamics code Athena++, with an isothermal shearing box setup, we study the boundary between the active and dead zone, where the accretion rate changes and mass can accumulate. We quantify how the turbulence level is affected by the presence of a non-uniform Ohmic resistivity in the radial x direction that leads to a region of inhibited turbulence (or dead zone). Comparing the turbulent activity to that of ideal simulations, the turbulence-inhibited area shows density fluctuations and magnetic activity at its boundaries, driven by energy injection from the active (ideal) zone boundaries. We find magnetic dissipation to be significantly stronger in the ideal regions, and the turbulence penetration through the boundary of the dead zone is determined by the value of the resistivity itself, through the Ohmic dissipation process, though the thickness of the transition does not play a significant role in changing the dissipation. We investigate the 1D spectra along the shearing direction: magnetic spectra appear flat at large scales both in ideal as well as resistive simulations, though a Kolmogorov scaling over more than one decade persists in the dead zone, suggesting the turbulent cascade is determined by the hydrodynamics of the system: magnetorotational instability dynamo action is inhibited where sufficiently high resistivity is present.
UR - http://www.scopus.com/inward/record.url?scp=85100858211&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85100858211&partnerID=8YFLogxK
U2 - 10.3847/1538-4357/abc9c0
DO - 10.3847/1538-4357/abc9c0
M3 - Article
AN - SCOPUS:85100858211
SN - 0004-637X
VL - 907
JO - Astrophysical Journal
JF - Astrophysical Journal
IS - 1
M1 - abc9c0
ER -