Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userStratified and vertically-shearing streaming instabilities in protoplanetary disks
104
0
0.0
(
0
)
Authors
Min-Kai Lin
Ask ChatGPT about the research
No Arabic abstract
Under the right conditions, the streaming instability between imperfectly coupled dust and gas is a powerful mechanism for planetesimal formation as it can concentrate dust grains to the point of gravitational collapse. In its simplest form, the streaming instability can be captured by analyzing the linear stability of unstratified disk models, which represent the midplane of protoplanetary disks. We extend such studies by carrying out vertically-global linear stability analyses of dust layers in protoplanetary disks. We find the dominant form of instability in stratified dust layers is one driven by the vertical gradient in the rotation velocity of the dust-gas mixture, but also requires partial dust-gas coupling. These vertically-shearing streaming instabilities grow on orbital timescales and occur on radial length scales $sim10^{-3}H_mathrm{g}$, where $H_mathrm{g}$ is the local pressure scale height. The classic streaming instability, associated with the relative radial drift between dust and gas, occur on radial length scales $sim10^{-2}H_mathrm{g}$, but have much smaller growth rates than vertically-shearing streaming instabilities. Including gas viscosity is strongly stabilizing and leads to vertically-elongated disturbances. We briefly discuss the potential effects of vertically-shearing streaming instabilities on planetesimal formation.
rate research
Read More
We use the Fokker-Planck equation and model the dispersive dynamics of solid particles in annular protoplanetary disks whose gas component is more massive than the particle phase. We model particle--gas interactions as hard sphere collisions, determine the functional form of diffusion coefficients, and show the existence of two global unstable modes in the particle phase. These modes have spiral patterns with the azimuthal wavenumber $m=1$ and rotate slowly. We show that in ring-shaped disks, the phase space density of solid particles increases linearly in time towards an accumulation point near the location of pressure maximum, while instabilities grow exponentially. Therefore, planetesimals and planetary cores can be efficiently produced near the peaks of unstable density waves. In this mechanism, particles migrating towards the accumulation point will not participate in the formation of planets, and should eventually form a debris ring like the main asteroid belt or classical Kuiper belt objects. We present the implications of global instabilities to the formation of ice giants and terrestrial planets in the solar system.
Large scale vortices could play a key role in the evolution of protoplanetary disks, particularly in the dead-zone where no turbulence associated with magnetic field is expected. Their possible formation by the subcritical baroclinic instability is a complex issue due to the vertical structure of the disk and to the elliptical instability.} {In two-dimensional disks the baroclinic instability is studied as a function of the thermal transfer efficiency. In three-dimensional disks we explore the importance of radial and vertical stratification on the processes of vortex formation and amplification.} {Numerical simulations are performed using a fully compressible hydrodynamical code based on a second order finite volume method. We assume a perfect gas law in inviscid disk models in which heat transfer is due to either relaxation or diffusion.} {In 2D, the baroclinic instability with thermal relaxation leads to the formation of large-scale vortices, which are unstable with respect to the elliptic instability. In the presence of heat diffusion, hollow vortices are formed which evolve into vortical structures with a turbulent core. In 3D, the disk stratification is found to be unstable in a finite layer which can include the mid-plane or not. When the unstable layer contains the mid-plane, the 3D baroclinic instability with thermal relaxation is found to develop first in the unstable layer as in 2D, producing large-scale vortices. These vortices are then stretched out in the stable layer, creating long-lived columnar vortical structures extending through the height of the disk. They are also found to be the source of internal vortex layers that develop across the whole disk along baroclinic critical layer surfaces, and form new vortices in the upper region of the disk.} {In three-dimensional disks, vortices can survive for a very long time if the production of vorticity by...
The streaming instability is a popular candidate for planetesimal formation by concentrating dust particles to trigger gravitational collapse. However, its robustness against physical conditions expected in protoplanetary disks is unclear. In particular, particle stirring by turbulence may impede the instability. To quantify this effect, we develop the linear theory of the streaming instability with external turbulence modelled by gas viscosity and particle diffusion. We find the streaming instability is sensitive to turbulence, with growth rates becoming negligible for alpha-viscosity parameters $alpha gtrsim mathrm{St} ^{1.5}$, where $mathrm{St}$ is the particle Stokes number. We explore the effect of non-linear drag laws, which may be applicable to porous dust particles, and find growth rates are modestly reduced. We also find that gas compressibility increase growth rates by reducing the effect of diffusion. We then apply linear theory to global models of viscous protoplanetary disks. For minimum-mass Solar nebula disk models, we find the streaming instability only grows within disk lifetimes beyond $sim 10$s of AU, even for cm-sized particles and weak turbulence ($alphasim 10^{-4}$). Our results suggest it is rather difficult to trigger the streaming instability in non-laminar protoplanetary disks, especially for small particles.
In the recent years, sub/mm observations of protoplanetary disks have discovered an incredible diversity of substructures in the dust emission. An important result was the finding that dust grains of mm size are embedded in very thin dusty disks. This implies that the dust mass fraction in the midplane becomes comparable to the gas, increasing the importance of the interaction between the two components there. We address this problem by means of numerical 2.5D simulations in order to study the gas and dust interaction in fully global stratified disks. To this purpose, we employ the recently developed dust grain module in the PLUTO code. Our model focuses on a typical T Tauri disk model, simulating a short patch of the disk at 10 au which includes grains of constant Stokes number of $St=0.01$ and $St=0.1$, corresponding to grains with sizes of 0.9 cm and 0.9 mm, respectively, for the given disk model. By injecting a constant pebble flux at the outer domain, the system reaches a quasi steady state of turbulence and dust concentrations driven by the streaming instability. For our given setup and using resolutions up to 2500 cells per scale height we resolve the streaming instability, leading to local dust clumping and concentrations. Our results show dust density values of around 10-100 times the gas density with a steady state pebble flux between $3.5 times 10^{-4}$ and $2.5 times 10^{-3} M_{rm Earth}/mathit{year}$ for the models with $mathit{St}=0.01$ and $mathit{St}=0.1$. The grain size and pebble flux for model $mathit{St}=0.01$ compares well with dust evolution models of the first million years of disk evolution. For those grains the scatter opacity dominates the extinction coefficient at mm wavelengths. These types of global dust and gas simulations are a promising tool for studies of the gas and dust evolution at pressure bumps in protoplanetary disks.
The streaming instability is a leading candidate mechanism to explain the formation of planetesimals. Yet, the role of this instability in the driving of turbulence in protoplanetary disks, given its fundamental nature as a linear hydrodynamical instability, has so far not been investigated in detail. We study the turbulence that is induced by the streaming instability as well as its interaction with the vertical shear instability. For this purpose, we employ the FLASH Code to conduct two-dimensional axisymmetric global disk simulations spanning radii from $1$ au to $100$ au, including the mutual drag between gas and dust as well as the radial and vertical stellar gravity. If the streaming instability and the vertical shear instability start their growth at the same time, we find the turbulence in the dust mid-plane layer to be primarily driven by the streaming instability. It gives rise to vertical gas motions with a Mach number of up to ${sim}10^{-2}$. The dust scale height is set in a self-regulatory manner to about $1%$ of the gas scale height. In contrast, if the vertical shear instability is allowed to saturate before the dust is introduced into our simulations, then it continues to be the main source of the turbulence in the dust layer. The vertical shear instability induces turbulence with a Mach number of ${sim}10^{-1}$ and thus impedes dust sedimentation. Nonetheless, we find the vertical shear instability and the streaming instability in combination to lead to radial dust concentration in long-lived accumulations which are significantly denser than those formed by the streaming instability alone. Thus, the vertical shear instability may promote planetesimal formation by creating weak overdensities that act as seeds for the streaming instability.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
