-
Formation of pebbles in (gravito-)viscous protoplanetary disks with various turbulent strengths
Authors:
Eduard I. Vorobyov,
Vardan G. Elbakyan,
Anders Johansen,
Michiel Lambrechts,
Aleksandr M. Skliarevskii,
O. P. Stoyanovskaya
Abstract:
Aims. Dust plays a crucial role in the evolution of protoplanetary disks. We study the dynamics and growth of initially sub-$μm$ dust particles in self-gravitating young protoplanetary disks with various strengths of turbulent viscosity. We aim to understand the physical conditions that determine the formation and spatial distribution of pebbles when both disk self-gravity and turbulent viscosity…
▽ More
Aims. Dust plays a crucial role in the evolution of protoplanetary disks. We study the dynamics and growth of initially sub-$μm$ dust particles in self-gravitating young protoplanetary disks with various strengths of turbulent viscosity. We aim to understand the physical conditions that determine the formation and spatial distribution of pebbles when both disk self-gravity and turbulent viscosity can be concurrently at work. Methods. We perform the thin-disk hydrodynamics simulations of self-gravitating protoplanetary disks over an initial time period of 0.5 Myr using the FEOSAD code. Turbulent viscosity is parameterized in terms of the spatially and temporally constant $α$-parameter, while the effects of gravitational instability on dust growth is accounted for by calculating the effective parameter $α_{\rm GI}$. We consider the evolution of dust component including momentum exchange with gas, dust self-gravity, and also a simplified model of dust growth. Results. We find that the level of turbulent viscosity strongly affects the spatial distribution and total mass of pebbles in the disk. The $α=10^{-2}$ model is viscosity-dominated, pebbles are completely absent, and dust-to-gas mass ratio deviates from the reference 1:100 value no more than by 30\% throughout the disk extent. On the contrary, the $α=10^{-3}$ model and, especially, the $α=10^{-4}$ model are dominated by gravitational instability. The effective parameter $α+α_{\rm GI}$ is now a strongly varying function of radial distance. As a consequence, a bottle neck effect develops in the innermost disk regions, which makes gas and dust accumulate in a ring-like structure. Abridged.
△ Less
Submitted 2 December, 2022;
originally announced December 2022.
-
Simulations of Dynamical Gas-Dust Circumstellar Disks: Going Beyond the Epstein Regime
Authors:
Olga P. Stoyanovskaya,
Fedor A. Okladnikov,
Eduard I. Vorobyov,
Yaroslav N. Pavlyuchenkov,
Vitaliy V. Akimkin
Abstract:
In circumstellar disks, the size of dust particles varies from submicron to several centimeters, while planetesimals have sizes of hundreds of kilometers. Therefore, various regimes for the aerodynamic drag between solid bodies and gas can be realized in these disks, depending on the grain sizes and velocities: Epstein, Stokes, and Newton, as well as transitional regimes between them. For small bo…
▽ More
In circumstellar disks, the size of dust particles varies from submicron to several centimeters, while planetesimals have sizes of hundreds of kilometers. Therefore, various regimes for the aerodynamic drag between solid bodies and gas can be realized in these disks, depending on the grain sizes and velocities: Epstein, Stokes, and Newton, as well as transitional regimes between them. For small bodies moving in the Epstein regime, the time required to establish the constant relative velocity between the gas and bodies can be much less than the dynamical time scale for the problem - the time for the rotation of the disk about the central body. In addition, the dust may be concentrated in individual regions of the disk, making it necessary to take into account the transfer of momentum between the dust and gas. It is shown that, for a system of equations for gas and monodisperse dust, a semi-implicit first-order approximation scheme in time in which the interphase interaction is calculated implicitly, while other forces, such as the pressure gradient and gravity are calculated explicitly, is suitable for stiff problems with intense interphase interactions and for computations of the drag in non-linear regimes. The piece-wise drag coefficient widely used in astrophysical simulations has a discontinuity at some values of the Mach and Knudsen numbers that are realized in a circumstellar disk. A continuous drag coefficient is presented, which corresponds to experimental dependences obtained for various drag regimes.
△ Less
Submitted 17 February, 2021;
originally announced February 2021.
-
Development and application of fast methods for computing momentum transfer between gas and dust in supercomputer simulation of planet formation
Authors:
Olga P. Stoyanovskaya,
Vitaly V. Akimkin,
Eduard I. Vorobyov,
Tatiana A. Glushko,
Yaroslav N. Pavlyuchenkov,
Valeriy N. Snytnikov,
Nikolay V. Snytnikov
Abstract:
Circumstellar discs, from which planetary systems are formed, consist of gas, dust and solids. Simulations of self-consistent dynamics of gas, dust and solids in circumstellar discs is a challenging problem. In the paper we present fast algorithms for computing the drag force (momentum transfer) between solid phase and gas. These algorithms (a) are universal and applicable to dust and solids with…
▽ More
Circumstellar discs, from which planetary systems are formed, consist of gas, dust and solids. Simulations of self-consistent dynamics of gas, dust and solids in circumstellar discs is a challenging problem. In the paper we present fast algorithms for computing the drag force (momentum transfer) between solid phase and gas. These algorithms (a) are universal and applicable to dust and solids with any sizes smaller than the mean free path of gas molecules, (b) can be used to calculate the momentum transfer between dust and gas instead of one-way effect, as it is done in many models, (c) can perform simulations, without a loss in accuracy, with the time step determined by gas-dynamic parameters rather than by drag force, and (d) are compatible with the widely used parallel algorithms for solving 3D equations of gas dynamics, hydrodynamic equations for dust, and the collisionless Boltzmann equation for large bodies. Preliminary results of supercomputer simulation of the gas-dust disc dynamics within the developed approach are reported.
△ Less
Submitted 15 November, 2018;
originally announced November 2018.
-
Two-Fluid Dusty Gas in Smoothed Particle Hydrodynamics: Fast and Implicit Algorithm for Stiff Linear Drag
Authors:
Olga P. Stoyanovskaya,
Tatiana A. Glushko,
Nikolay V. Snytnikov,
Valeriy N. Snytnikov
Abstract:
Simulation of the dynamics of dust-gas circumstellar discs is crucial in understanding the mechanisms of planet formation. The dynamics of small grains in the disc is stiffly coupled to the gas, while the dynamics of grown solids is decoupled. Moreover, in some parts of the disc the concentration of the dust is low (dust to gas mass ratio is about 0.01), while in other parts it can be much higher.…
▽ More
Simulation of the dynamics of dust-gas circumstellar discs is crucial in understanding the mechanisms of planet formation. The dynamics of small grains in the disc is stiffly coupled to the gas, while the dynamics of grown solids is decoupled. Moreover, in some parts of the disc the concentration of the dust is low (dust to gas mass ratio is about 0.01), while in other parts it can be much higher. These factors place high requirements on the numerical methods for disc simulations. In particular, when gas and dust are simulated with two different fluids, explicit methods require very small timestep (must be less than dust stopping time $t_{\rm stop}$ during which the velocity of a solid particle is equalized with respect to the gas velocity) to obtain solution, while some implicit methods requires high temporal resolution to obtain acceptable accuracy. Moreover, recent studies underlined that for Smoothed particle hydrodynamics (SPH) when the gas and the dust are simulated with different sets of particles only high spatial resolution $h<c_{\rm s} t_{\rm stop}$ guaranties suppression of numerical overdissipation due to gas and dust interaction.
To address these problems, we developed a fast algorithm based on the ideas of (1) implicit integration of linear (Epstein) drag and (2) exact conservation of local linear momentum. We derived formulas for monodisperse dust-gas in two-fluid SPH and tested the new method on problems with known analytical solutions. We found that our method is a promising alternative for the previously developed two-fluid SPH scheme in case of stiff linear drag thanks to the fact that spatial resolution condition $h<c_{\rm s} t_{\rm stop}$ is not required anymore for accurate results.
△ Less
Submitted 15 November, 2018;
originally announced November 2018.
-
Modeling Circumstellar Disc Fragmentation and Episodic Protostellar Accretion with Smoothed Particle Hydrodynamics in Cell
Authors:
Olga P. Stoyanovskaya,
Nikolay V. Snytnikov,
Valeriy N. Snytnikov
Abstract:
We discuss the ability of the smoothed particle hydrodynamics (SPH) method combined with a grid-based solver for the Poisson equation to model mass accretion onto protostars in gravitationally unstable protostellar discs. We scrutinize important features of coupling the SPH with grid-based solvers and numerical issues associated with (1) large number of SPH neighbours and (2) relation between grav…
▽ More
We discuss the ability of the smoothed particle hydrodynamics (SPH) method combined with a grid-based solver for the Poisson equation to model mass accretion onto protostars in gravitationally unstable protostellar discs. We scrutinize important features of coupling the SPH with grid-based solvers and numerical issues associated with (1) large number of SPH neighbours and (2) relation between gravitational softening and hydrodynamic smoothing length.
We report results of our simulations of razor-thin disc prone to fragmentation and demonstrate that the algorithm being simple and homogeneous captures the target physical processes - disc gravitational fragmentation and accretion of gas onto the protostar caused by inward migration of dense clumps.
In particular, we obtain two types of accretion bursts: a short-duration one caused by a quick inward migration of the clump, previously reported in the literature, and the prolonged one caused by the clump lingering at radial distances on the order of 15-25~au. The latter is culminated with a sharp accretion surge caused by the clump ultimately falling on the protostar.
△ Less
Submitted 4 September, 2018;
originally announced September 2018.
-
Analysis of Methods for Computing the Trajectories of Dust Particles in a Gas-Dust Circumstellar Disk
Authors:
Olga P. Stoyanovskaya,
Valeriy N. Snytnikov,
Eduard I. Vorobyov
Abstract:
A systematic analysis of methods for computing the trajectories of solid-phase particles applied in modern astrophysics codes designed for modeling gas-dust circumstellar disks has been carried out for the first time. The motion of grains whose velocities are determined mainly by the gas drag, that is, for which the stopping time or relaxation time for the velocity of the dust to the velocity of t…
▽ More
A systematic analysis of methods for computing the trajectories of solid-phase particles applied in modern astrophysics codes designed for modeling gas-dust circumstellar disks has been carried out for the first time. The motion of grains whose velocities are determined mainly by the gas drag, that is, for which the stopping time or relaxation time for the velocity of the dust to the velocity of the gas t_{stop} is less than or comparable to the rotation period, are considered. The methods are analyzed from the point of view of their stability for computing the motions of small bodies, including dust grains less than 1 micron, which are strongly coupled to the gas. Two test problems are with analytical solutions. Fast first order accurate methods that make it possible to avoid additional restriction on the time step due to gas drag in computations of the motion of grains of any size are presented. For the conditions of a circumstellar disks, the error in the velocity computations obtained when using some stable methods become unacceptably large when the time step size τ>t_{stop}. For the radial migration of bodies that exhibit drifts along nearly Keplerian orbits, an asymptotic approximation, sometimes called the short friction time approximation or drift flux model, gives a relative error for the radial velocity equals to St^2 where St is the Stokes number, the ratio of the stopping time of the body to some fraction of the rotation period (dynamical time scale) in the disk.
△ Less
Submitted 4 September, 2018;
originally announced September 2018.
-
Analysis of Numerical Algorithms for Computing Rapid Momentum Transfers between the Gas and Dust in Simulations of Circumstellar Disks
Authors:
O. P. Stoyanovskaya,
E. I. Vorobyov,
V. N. Snytnikov
Abstract:
Approaches used in modern numerical simulations of the dynamics of dust and gas in circumstellar disks are tested. The gas and dust are treated like interpenetrating continuous media that can exchange momentum. A stiff coupling between the gas and dust phases is typical for such disks, with the dust stopping time much less than the characteristic dynamical time scale. This imposes high demands on…
▽ More
Approaches used in modern numerical simulations of the dynamics of dust and gas in circumstellar disks are tested. The gas and dust are treated like interpenetrating continuous media that can exchange momentum. A stiff coupling between the gas and dust phases is typical for such disks, with the dust stopping time much less than the characteristic dynamical time scale. This imposes high demands on the methods used to simulate the dust dynamics. A grid, piecewise-parabolic method is used as the basic algorithm for solving the gas-dynamical equations. Numerical solutions obtained using various methods to compute the momentum exchanges are presented for the case of monodisperse dust. Numerical solutions are obtained for shock tube problem and the propagation of sound waves in a gas-dust medium. The studied methods are compared in terms of their ability to model media with (a) an arbitrary (short or long) dust stopping time, and (b) an arbitrary dust concentration in the gas (varying the dust to gas mass ratio from 0.01 to 1). A method for computing the momentum exchange with infinite-order accuracy in time is identified, which makes it possible to satisfy the conditions (a) and (b) with minimal computational costs. A first-order method that shows similar results in the test computations is also presented. It is shown that the proposed first-order method for monodisperse dust can be extended to a regime when the dust is polydisperse; i.e., a regime represented by several fractions with different stopping times. Formulas for computing the gas and dust velocities for polydisperse dust with each fraction exchanging momentum with the gas are presented.
△ Less
Submitted 8 August, 2018;
originally announced August 2018.
-
Clump formation due to the gravitational instability of a multiphase medium in a massive protoplanetary disc
Authors:
Valeriy N. Snytnikov,
Olga P. Stoyanovskaya
Abstract:
Planetary systems form in gas-dust protoplanetary discs via the growth of solid bodies. In this paper, we show that the most intriguing stage of such growth --- namely, the transformation of 1-10 m boulders into kilometre-sized planetesimals --- can be explained by a mechanism of gravitational instability. The present work focused on the origin of self-gravitating clumps in which planetesimal form…
▽ More
Planetary systems form in gas-dust protoplanetary discs via the growth of solid bodies. In this paper, we show that the most intriguing stage of such growth --- namely, the transformation of 1-10 m boulders into kilometre-sized planetesimals --- can be explained by a mechanism of gravitational instability. The present work focused on the origin of self-gravitating clumps in which planetesimal formation could take place. Our computer simulations demonstrated that such clumps of gas and boulders formed due to the development of a two-phase instability. This instability revealed a so-called 'mutual influence effect' in the protoplanetary disc, where the dynamics of the system were determined by the collisionless collective motion of a low-mass subdisc composed of primary solids. We found that a 0.1 $c_s$ velocity dispersion in the boulder subdisc was sufficient to cause the formation of self-gravitating clumps of gas and boulders. In such regimes, the time needed for the formation of the collapsing objects was less than the boulders' dissipation time in the density waves of the medium.
△ Less
Submitted 3 October, 2012;
originally announced October 2012.