  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.
