Turbulence in particle laden midplane layers of planet forming disks

Abstract

We examine the settled particle layers of planet forming disks in which the streaming instability (SI) is thought to be either weak or inactive. A suite of low-to-moderate resolution three-dimensional simulations in a $0.2H$ sized box, where $H$ is the pressure scale height, are performed using PENCIL for two Stokes numbers, \St$=0.04$ and $0.2$, at 1\% disk metallicity. We find a complex of Ekman-layer jet-flows emerge subject to three co-acting linearly growing processes: (1) the Kelvin-Helmholtz instability (KHI), (2) the planet-forming disk analog of the baroclinic Symmetric Instability (SymI), and (3) a later-time weakly acting secondary transition process, possibly a manifestation of the SI, producing a radially propagating pattern state. For \St$=0.2$, KHI is dominant and manifests as off-midplane axisymmetric rolls, while for \St$=0.04$ the axisymmetric SymI mainly drives turbulence. SymI is analytically developed in a model disk flow, predicting that it becomes strongly active when the Richardson number (Ri) of the particle-gas midplane layer transitions below 1, exhibiting growth rates $\le\sqrt{2/\Ri - 2}\cdotΩ$, where $\Omega$ is local disk rotation rate. For fairly general situations absent external sources of turbulence it is conjectured that the SI, when and if initiated, emerges out of a turbulent state primarily driven and shaped by at least SymI and/or KHI. We also find that turbulence produced in $256^3$ resolution simulations are not statistically converged and that corresponding $512^3$ simulations may be converged for \St$=0.2$. Furthermore, we report that our numerical simulations significantly dissipate turbulent kinetic energy on scales less than 6-8 grid points.

55 pages, 27 figures, accepted for publication in ApJ

References in corpus (15)

• Particle Stirring in Turbulent Gas Disks: Including Orbital Oscillations
• How to form planetesimals from mm-sized chondrules and chondrule aggregates
• Dust sedimentation and self-sustained Kelvin-Helmholtz turbulence in protoplanetary disk mid-planes. I. Radially symmetric simulations
• Thresholds for Particle Clumping by the Streaming Instability
• Inertial range scaling in numerical turbulence with hyperviscosity
• Vertical Shearing Instabilities in Radially Shearing Disks: The Dustiest Layers of the Protoplanetary Nebula
• A thermodynamic view of dusty protoplanetary disks
• Local models of astrophysical discs
• A weakly nonlinear analysis of the magnetorotational instability in a model channel flow
• Three-Dimensional Simulations of Kelvin-Helmholtz Instability in Settled Dust Layers in Protoplanetary Disks
• Pebble trapping in vortices: three-dimensional simulations
• Stratified and vertically-shearing streaming instabilities in protoplanetary disks
• Orbital advection with magnetohydrodynamics and vector potential
• On the settling of small grains in dusty discs: analysis and formulas
• A minimal model for vertical shear instability in protoplanetary accretion disks

Cited by in corpus (4)

• Vertical shear instability in two-moment radiation-hydrodynamical simulations of irradiated protoplanetary disks II. Secondary instabilities and stability regions
• Dust-gas dynamics driven by the streaming instability with various pressure gradients
• Formation of the First Planetesimals via the Streaming Instability in Globally Turbulent Protoplanetary Disks?
• Properties of outer solar system pebbles during planetesimal formation from meteor observations