Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userProtoplanetary Disk Turbulence Driven by the Streaming Instability: Non-Linear Saturation and Particle Concentration
280
0
0.0
(
0
)
Authors
Anders Johansen MPIA Heidelberg
Ask ChatGPT about the research
No Arabic abstract
We present simulations of the non-linear evolution of streaming instabilities in protoplanetary disks. The two components of the disk, gas treated with grid hydrodynamics and solids treated as superparticles, are mutually coupled by drag forces. We find that the initially laminar equilibrium flow spontaneously develops into turbulence in our unstratified local model. Marginally coupled solids (that couple to the gas on a Keplerian time-scale) trigger an upward cascade to large particle clumps with peak overdensities above 100. The clumps evolve dynamically by losing material downstream to the radial drift flow while receiving recycled material from upstream. Smaller, more tightly coupled solids produce weaker turbulence with more transient overdensities on smaller length scales. The net inward radial drift is decreased for marginally coupled particles, whereas the tightly coupled particles migrate faster in the saturated turbulent state. The turbulent diffusion of solid particles, measured by their random walk, depends strongly on their stopping time and on the solids-to-gas ratio of the background state, but diffusion is generally modest, particularly for tightly coupled solids. Angular momentum transport is too weak and of the wrong sign to influence stellar accretion. Self-gravity and collisions will be needed to determine the relevance of particle overdensities for planetesimal formation.
rate research
Read More
We present local simulations that verify the linear streaming instability that arises from aerodynamic coupling between solids and gas in protoplanetary disks. This robust instability creates enhancements in the particle density in order to tap the free energy of the relative drift between solids and gas, generated by the radial pressure gradient of the disk. We confirm the analytic growth rates found by Youdin & Goodman (2005) using grid hydrodynamics to simulate the gas and, alternatively, particle and grid representations of the solids. Since the analytic derivation approximates particles as a fluid, this work corroborates the streaming instability when solids are treated as particles. The idealized physical conditions -- axisymmetry, uniform particle size, and the neglect of vertical stratification and collisions -- provide a rigorous, well-defined test of any numerical algorithm for coupled particle-gas dynamics in protoplanetary disks. We describe a numerical particle-mesh implementation of the drag force, which is crucial for resolving the coupled oscillations. Finally we comment on the balance of energy and angular momentum in two-component disks with frictional coupling. A companion paper details the non-linear evolution of the streaming instability into saturated turbulence with dense particle clumps.
The Streaming Instability (SI) is a mechanism to concentrate solids in protoplanetary disks. Nonlinear particle clumping from the SI can trigger gravitational collapse into planetesimals. To better understand the numerical robustness of the SI, we perform a suite of vertically-stratified 3D simulations with fixed physical parameters known to produce strong clumping. We vary the numerical implementation, namely the computational domain size and the vertical boundary conditions (vBCs), comparing newly-implemented outflow vBCs to the previously-used periodic and reflecting vBCs. We find strong particle clumping by the SI is mostly independent of the vBCs. However, peak particle densities are higher in larger simulation domains due to a larger particle mass reservoir. We report SI-triggered zonal flows, i.e., azimuthally-banded radial variations of gas pressure. These structures have low amplitudes, insufficient to halt particle radial drift, confirming that particle trapping in gas pressure maxima is not the mechanism of the SI. We find that outflow vBCs produce artificially large gas outflow rates at vertical boundaries. However, the outflow vBCs reduce artificial reflections at vertical boundaries, allowing more particle sedimentation, and showing less temporal variation and better convergence with box size. The radial spacing of dense particle filaments is $sim0.15$ gas scale heights ($H$) for all vBCs, which sets the feeding zone for planetesimal growth in self-gravitating simulations. Our results validate the use of the outflow vBCs in SI simulations, even with vertical boundaries close ($leq 0.4H$) to the disk midplane. Overall, our study demonstrates the numerical robustness of nonlinear particle clumping by the SI.
The formation of the first planetesimals and the final growth of planetary cores relies on the abundance of small pebbles. The efficiencies of both the streaming instability (SI) process, suggested to catalyze the early growth of planetesimals, and the pebble-accretion process, suggested to accelerate the growth of planetary cores, depend on the sizes of solids residing in the disk. In particular, these processes were found to be sensitive to size distribution of solids, and efficient planetesimal formation and growth through these channels require a limited pebble size distribution. Here we show that aeolian erosion, a process that efficiently grinds down boulders into a mono-sized distribution of pebbles, provides a natural upper limit for the maximal pebble sizes (in terms of their Stokes number). We find the dependence of this upper limit on the radial separation, disk age, turbulence strength, and the grain-size composition of the boulders in the disk. SI is favorable in areas with a Stokes number less than 0.1, which is found in the inner sub-astronomical-unit regions of the disk. This upper limit shapes the size distribution of small pebbles and thereby catalyzes the early onset of planetesimal formation due to SI, and the later core accretion growth through pebble accretion.
Recent years have seen growing interest in the streaming instability as a candidate mechanism to produce planetesimals. However, these investigations have been limited to small-scale simulations. We now present the results of a global protoplanetary disk evolution model that incorporates planetesimal formation by the streaming instability, along with viscous accretion, photoevaporation by EUV, FUV, and X-ray photons, dust evolution, the water ice line, and stratified turbulence. Our simulations produce massive (60-130 $M_oplus$) planetesimal belts beyond 100 au and up to $sim 20 M_oplus$ of planetesimals in the middle regions (3-100 au). Our most comprehensive model forms 8 $M_oplus$ of planetesimals inside 3 au, where they can give rise to terrestrial planets. The planetesimal mass formed in the inner disk depends critically on the timing of the formation of an inner cavity in the disk by high-energy photons. Our results show that the combination of photoevaporation and the streaming instability are efficient at converting the solid component of protoplanetary disks into planetesimals. Our model, however, does not form enough early planetesimals in the inner and middle regions of the disk to give rise to giant planets and super-Earths with gaseous envelopes. Additional processes such as particle pileups and mass loss driven by MHD winds may be needed to drive the formation of early planetesimal generations in the planet forming regions of protoplanetary disks.
In continuation of previous work, numerical results are presented, concerning relativistically counter-streaming plasmas. Here, the relativistic mixed mode instability evolves through, and beyond, the linear saturation -- well into the nonlinear regime. Besides confirming earlier findings, that wave power initially peaks on the mixed mode branch, it is observed that, during late time evolution wave power is transferred to other wave numbers. It is argued that the isotropization of power in wavenumber space may be a consequence of weak turbulence. Further, some modifications to the ideal weak turbulence limit is observed. Development of almost isotropic predominantly electrostatic -- partially electromagnetic -- turbulent spectra holds relevance when considering the spectral emission signatures of the plasma, namely bremsstrahlung, respectively magneto-bremsstrahlung (synchrotron radiation and jitter radiation) from relativistic shocks in astrophysical jets and shocks from gamma-ray bursts and active galactic nuclei.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
