Analysis of Methods for Computing the Trajectories of Dust Particles in a Gas-Dust Circumstellar Disk

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 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.

arXiv comment21 pages, 8 figures

References in corpus (18)

• Coagulation, fragmentation and radial motion of solid particles in protoplanetary disks
• Radial Flow of Dust Particles in Accretion Disks
• Dust Filtration by Planet-Induced Gap Edges: Implications for Transitional Disks
• Accelerated planetesimal growth in self-gravitating protoplanetary discs
• Dust diffusion in protoplanetary discs by magnetorotational turbulence
• The minimum mass of detectable planets in protoplanetary discs and the derivation of planetary masses from high resolution observations
• Particle-Gas Dynamics with Athena: Method and Convergence
• A numerical simulation of a "super-Earth" core delivery from ~ 100 AU to ~ 8 AU
• Dust distribution in protoplanetary disks - Vertical settling and radial migration
• Turbulence-Induced Relative Velocity of Dust Particles I: Identical Particles
• Dust and gas mixtures with multiple grain species - a one-fluid approach
• DUSTYBOX and DUSTYWAVE: Two test problems for numerical simulations of two fluid astrophysical dust-gas mixtures
• Embedded protostellar disks around (sub-)solar protostars. I. Disk structure and evolution
• Grand challenges in protoplanetary disc modelling
• Two fluid dust and gas mixtures in SPH: A semi-implicit approach
• Two-fluid dust and gas mixtures in smoothed particle hydrodynamics II: an improved semi-implicit approach
• Effects of photophoresis on the dust distribution in a 3D protoplanetary disc
• Clump formation due to the gravitational instability of a multiphase medium in a massive protoplanetary disc

Cited by in corpus (7)

• The early evolution of viscous and self-gravitating circumstellar disks with a dust component
• Small dust grain dynamics on adaptive mesh-refinement grids. I. Methods
• Two-Fluid Dusty Gas in Smoothed Particle Hydrodynamics: Fast and Implicit Algorithm for Stiff Linear Drag
• Torques felt by solid accreting planets
• Formation of pebbles in (gravito-)viscous protoplanetary disks with various turbulent strengths
• Development and application of fast methods for computing momentum transfer between gas and dust in supercomputer simulation of planet formation
• Gravitational Instability of Gas-Dust Circumnuclear Disks in Galaxies