Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userPolydisperse Streaming Instability III. Dust evolution encourages fast instability
59
0
0.0
(
0
)
Added by
Sijme-Jan Paardekooper
Publication date
2021
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
Planet formation via core accretion requires the production of km-sized planetesimals from cosmic dust. This process must overcome barriers to simple collisional growth, for which the Streaming Instability (SI) is often invoked. Dust evolution is still required to create particles large enough to undergo vigorous instability. The SI has been studied primarily with single size dust, and the role of the full evolved dust distribution is largely unexplored. We survey the Polydispserse Streaming Instability (PSI) with physical parameters corresponding to plausible conditions in protoplanetary discs. We consider a full range of particle stopping times, generalized dust size distributions, and the effect of turbulence. We find that, while the PSI grows in many cases more slowly with a interstellar power-law dust distribution than with a single size, reasonable collisional dust evolution, producing an enhancement of the largest dust sizes, produces instability behaviour similar to the monodisperse case. Considering turbulent diffusion the trend is similar. We conclude that if fast linear growth of PSI is required for planet formation, then dust evolution producing a distribution with peak stopping times on the order of 0.1 orbits and an enhancement of the largest dust significantly above the single power-law distribution produced by a fragmentation cascade is sufficient, along with local enhancement of the dust to gas volume mass density ratio to order unity.
rate research
Read More
Occurring in protoplanetary discs composed of dust and gas, streaming instabilities are a favoured mechanism to drive the formation of planetesimals. The Polydispserse Streaming Instability is a generalisation of the Streaming Instability to a continuum of dust sizes. This second paper in the series provides a more in-depth derivation of the governing equations and presents novel numerical methods for solving the associated linear stability problem. In addition to the direct discretisation of the eigenproblem at second order introduced in the previous paper, a new technique based on numerically reducing the system of integral equations to a complex polynomial combined with root finding is found to yield accurate results at much lower computational cost. A related method for counting roots of the dispersion relation inside a contour without locating those roots is also demonstrated. Applications of these methods show they can reproduce and exceed the accuracy of previous results in the literature, and new benchmark results are provided. Implementations of the methods described are made available in an accompanying Python package psitools.
Damping of the previously discovered resonant drag instability (RDI) of dust streaming in protoplanetary disc is studied using the local approach to dynamics of gas-dust perturbations in the limit of the small dust fraction. Turbulence in a disc is represented by the effective viscosity and diffusivity in equations of motion for gas and dust, respectively. In the standard case of the Schmidt number (ratio of the effective viscosity to diffusivity) Sc = 1, the reduced description of RDI in terms of the inertial wave (IW) and the streaming dust wave (SDW) falling in resonance with each other reveals that damping solution differs from the inviscid solution simply by adding the characteristic damping frequency to its growth rate. RDI is fully suppressed at the threshold viscosity, which is estimated analytically, first, for radial drift, next, for vertical settling of dust, and at last, in the case of settling combined with radial drift of the dust. In the last case, RDI survives up to the highest threshold viscosity, with a greater excess for smaller solids. Once Sc eq 1, a new instability specific for dissipative perturbations on the dust settling background emerges. This instability of the quasi-resonant nature is referred to as settling viscous instability (SVI). The mode akin to SDW (IW) becomes growing in a region of long waves provided that Sc > 1 (Sc < 1). SVI leads to an additional increase of the threshold viscosity.
Recent study suggests that the streaming instability, one of the leading mechanisms for driving the formation of planetesimals, may not be as efficient as previously thought. Under some disc conditions, the growth timescale of the instability can be longer than the disc lifetime when multiple dust species are considered. To further explore this finding, we use both linear analysis and direct numerical simulations with gas fluid and dust particles to mutually validate and study the unstable modes of the instability in more detail. We extend the previously studied parameter space by one order of magnitude in both the range of the dust-size distribution $[T_{s,min}, T_{s,max}]$ and the total solid-to-gas mass ratio $varepsilon$ and introduce a third dimension with the slope $q$ of the size distribution. We find that the fast-growth regime and the slow-growth regime are distinctly separated in the $varepsilon$-$T_{s,max}$ space, while this boundary is not appreciably sensitive to $q$ or $T_{s,min}$. With a wide range of dust sizes present in the disc (e.g. $T_{s,min}lesssim10^{-3}$), the growth rate in the slow-growth regime decreases as more dust species are considered. With a narrow range of dust sizes (e.g. $T_{s,max}/T_{s,min}=5$), on the other hand, the growth rate in most of the $varepsilon$-$T_{s,max}$ space is converged with increasing dust species, but the fast and the slow growth regimes remain clearly separated. Moreover, it is not necessary that the largest dust species dominate the growth of the unstable modes, and the smaller dust species can affect the growth rate in a complicated way. In any case, we find that the fast-growth regime is bounded by $varepsilongtrsim 1$ or $T_{s,max}gtrsim 1$, which may represent the favourable conditions for planetesimal formation.
The recently discovered resonant drag instability (RDI) of dust streaming in protoplanetary disc is considered as the mode coupling of subsonic gas-dust mixture perturbations. This mode coupling is coalescence of two modes with nearly equal phase velocities: inertial wave (IW) having positive energy and a streaming dust wave (SDW) having negative energy as measured in the frame of gas environment being at rest in vertical hydrostatic equilibrium. SDW is a trivial mode produced by the bulk streaming of dust, which transports perturbations of dust density. In this way, settling combined with radial drift of the dust makes possible coupling of SDW with IW and the onset of the instability. In accordance with the concept of the mode coupling, RDI growth rate is proportional to the square root of the coupling term of the dispersion equation, which itself is proportional to mass fraction of dust, $fll 1$. This clarifies why RDI growth rate $propto f^{1/2}$. When SDW has positive energy, its resonance with IW provides an avoided crossing instead of the mode coupling. In the high wavenumber limit RDI with unbounded growth rate $propto f^{1/3}$ is explained by the triple mode coupling, which is coupling of SDW with two IW. It coexists with a new quasi-resonant instability accompanied by bonding of two oppositely propagating low-frequency IW. The mode coupling does not exist for dust streaming only radially in a disc. In this case RDI is provided by the obscured mechanism associated with the inertia of solids.
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
