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

Karim Shariff,Orkan M. Umurhan

2024-11-06

Abstract:A high resolution fourth-order Padé scheme is used to simulate locally isothermal 3D disk turbulence driven by the vertical shear instability (VSI) using 268.4 M points. In the early non-linear period of axisymmetric VSI, angular momentum transport by vertical jets creates correlated N-shaped radial profiles of perturbation vertical and azimuthal velocity. This implies dominance of positive perturbation vertical vorticity layers and a recently discovered angular momentum staircase with respect to radius ($r$). These features are present in 3D in a weaker form. The 3D flow consists of vertically and azimuthally coherent turbulent shear layers with small vortices with all three components are active. Previously observed large persistent vortices in the interior of the domain driven by the Rossby wave instability are absent. We speculate that this is due to a weaker angular momentum staircase in 3D in the present simulations than in a previous simulation. The turbulent viscosity parameter $\alpha(r)$ increases linearly with $r$. At intermediate resolution, the value of $\alpha(r)$ at mid-radius is close to that of a previous simulation. The specific kinetic energy spectrum with respect to radial wavenumber has a power law region with exponent $-1.84$, close to the value $-2$ expected for shear layers. The spectrum with respect to azimuthal wavenumber has a $-5/3$ region and lacks a $-5$ region reported in an earlier study. Finally, it is found axisymmetric VSI has artifacts at late times, including a very strong angular momentum staircase, which in 3D is present weakly in the upper layers of the disk.

Solar and Stellar Astrophysics,Earth and Planetary Astrophysics,Fluid Dynamics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is about the characteristics of turbulence driven by the Vertical Shear Instability (VSI) in protoplanetary disks (PPDs) and its impact on the planet - formation process. Specifically, the research aims to:

1. **Understand the structure and intensity of VSI - driven turbulence**: Through high - resolution numerical simulations, researchers hope to gain in - depth understanding of the structure and intensity of VSI - driven turbulence at different scales. This includes how turbulence affects angular momentum transport, what types of vortex structures are formed, and how these structures evolve over time.

2. **Explore the role of VSI in the aggregation of particles into planetesimals**: The turbulence generated by VSI may play a crucial role in the aggregation of micron - sized particles into kilometer - sized planetesimals. Therefore, understanding the specific characteristics of this turbulence helps to explain the formation mechanism of planetesimals in the early solar system.

3. **Verify previous theories and simulation results**: Previous studies have shown that VSI may be one of the most important turbulence - generating mechanisms in protoplanetary disks, especially in the mid - plane region of the disk. This study hopes to verify these theoretical predictions through high - resolution simulations and reveal the behavioral differences of VSI in a three - dimensional environment.

4. **Analyze the variation law of the turbulent viscosity parameter α**: It has been found that the turbulent viscosity parameter α increases linearly with the radius. Through simulations, this parameter can be measured more accurately, and its potential impact on mass transport and disk evolution can be explored.

5. **Explore the interaction between VSI and other instabilities**: In addition to VSI, there are other types of instabilities in protoplanetary disks (such as Rossby wave instability). Studying the interaction between VSI and these instabilities helps to comprehensively understand the dynamic processes in the disk.

### Key issues in the abstract

- **Turbulence structure**: The paper describes the turbulence structure driven by VSI, especially the N - shaped radial profile formed by vertical jets and the positive - perturbed vertical vorticity layer.
- **Angular momentum ladder**: The angular momentum ladder phenomenon has been discovered, especially more obvious in the axisymmetric case.
- **Three - dimensional effects**: Three - dimensional simulations show that although there are similar features, the intensity is weaker, and no large - scale persistent vortices previously reported have been observed.
- **Turbulent viscosity α parameter**: The study shows that the α parameter increases linearly with the radius and is consistent with previous simulation results at medium resolution.
- **Kinetic energy spectrum**: The kinetic energy spectrum of the radial wavenumber shows a power - law distribution with an exponent of - 1.84, close to the expected shear - layer value; while the kinetic energy spectrum of the azimuthal wavenumber shows a power - law distribution of - 5/3.

In summary, through high - resolution numerical simulations, this paper discusses in detail the characteristics of VSI - driven turbulence and its impact on the formation of planetesimals in protoplanetary disks, fills the gaps in existing research, and provides important reference data for further research.