Vertical shear instability in two-moment radiation-hydrodynamical simulations of irradiated protoplanetary disks I. Angular momentum transport and turbulent heating (2312.06882v1)

Abstract: We studied the linear and nonlinear evolution of the Vertical Shear Instability (VSI) in axisymmetric models of protoplanetary disks, focusing on the transport of angular momentum, the produced temperature perturbations, and the applicability of local stability conditions. We modeled the gas-dust mixture via high-resolution two-moment (M1) radiation-hydrodynamical simulations including stellar irradiation with frequency-dependent opacities. We found that, given sufficient depletion of small grains (with a dust-to-gas mass ratio of $10\%$ of our nominal value of $10^{-3}$ for $<0.25$ $\mu$m grains), the VSI can operate in surface disk layers while being inactive close to the midplane, resulting in a suppression of the VSI body modes. The VSI reduces the initial vertical shear in bands of approximately uniform specific angular momentum, whose formation is likely favored by the enforced axisymmetry. Similarities with Reynolds stresses and angular momentum distributions in 3D simulations suggest that the VSI-induced angular momentum mixing in the radial direction may be predominantly axisymmetric. The stability regions in our models are well explained by local stability criteria, while the employment of global criteria is still justifiable up to a few scale heights above the midplane, at least as long as VSI modes are radially optically thin. Turbulent heating produces only marginal temperature increases of at most $0.1\%$ and $0.01\%$ in the nominal and dust-depleted models, respectively, peaking at a few (approximately three) scale heights above the midplane. We conclude that it is unlikely that the VSI can, in general, lead to any significant temperature increase since that would either require it to efficiently operate in largely optically thick disk regions or to produce larger levels of turbulence than predicted by models of passive irradiated disks.