Abstract:Context. Gas in protoplanetary disks mostly cools via thermal accommodation with dust particles. Thermal relaxation is thus highly sensitive to the local dust size distributions and the spatial distribution of the grains. So far, the interplay between thermal relaxation and gas turbulence has not been dynamically modeled in hydrodynamic simulations of protoplanetary disks with dust. Aims. We aim to study the effects of the vertical shear instability (VSI) on the thermal relaxation times, and vice versa. We are particularly interested in the influence of the initial dust grain size on the VSI and whether the emerging turbulence is sustained over long timescales. Methods. We ran three axisymmetric hydrodynamic simulations of a protoplanetary disk including four dust fluids that initially resemble MRN size distributions of different initial grain sizes. From the local dust densities, we calculated the thermal accommodation timescale of dust and gas and used the result as the thermal relaxation time of the gas in our simulation. We included the effect of dust growth by applying the monodisperse dust growth rate and the typical growth limits. Results. We find that the emergence of the VSI is strongly dependent on the initial dust grain size. Coagulation also counteracts the emergence of hydrodynamic turbulence in our simulations, as shown by others before. Starting a simulation with larger grains (100 μm) generally leads to a less turbulent outcome. While the inner disk regions (within ∼70 au) develop turbulence in all three simulations, we find that the simulations with larger particles do not develop VSI in the outer disk. Conclusions. Our simulations with dynamically calculated thermal accommodation times based on the drifting and settling dust distribution show that the VSI, once developed in a disk, can be sustained over long timescales, even if grain growth is occurring. The VSI corrugates the dust layer and even diffuses the smaller grains into the upper atmosphere, where they can cool the gas. Whether the instability can emerge for a specific stratification depends on the initial dust grain sizes and the initial dust scale height. If the grains are initially ≳100 μm and if the level of turbulence is initially assumed to be low, we find no VSI turbulence in the outer disk regions.

On corrugation mode radial wavelengths of the vertical shear instability

On corrugation mode radial wavelengths of the vertical shear instability

Vertical shear instability with partially reflecting boundary conditions

Vertical shear instability in two-moment radiation-hydrodynamical simulations of irradiated protoplanetary disks II. Secondary instabilities and stability regions

A high resolution simulation of protoplanetary disk turbulence driven by the vertical shear instability

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

Thermal Structure Determines Kinematics: Vertical Shear Instability in Stellar Irradiated Protoplanetary Disks

Parametric instability in warped astrophysical discs: growth, saturation and feedback

Two saturated states of the vertical shear instability in protoplanetary disks with vertically varying cooling times

Global Simulations of the Vertical Shear Instability with Non-ideal Magnetohydrodynamical Effects

Global Simulations of the Vertical Shear Instability with Nonideal Magnetohydrodynamic Effects

The vertical shear instability in poorly ionised, magnetized protoplanetary discs

Thermal instabilities in accretion disks II: Numerical Experiments for the Goldreich-Schubert-Fricke Instability and the Convective Overstability in disks around young stars

Vertical shear instability with dust evolution and consistent cooling times

A minimal model for vertical shear instability in protoplanetary accretion disks

Spiral density waves and vertical circulation in protoplanetary discs

Vertical shear instability with dust evolution and consistent cooling times. On the importance of the initial dust distribution

Tearing Mode with Strong Flow Shear in the Viscosity-Dominated Limit

A parameter study of strato-rotational low-frequency modulations: impacts on momentum transfer and energy distribution