Aims. The long-term evolution of a circumstellar disk starting from its formation and ending in the T Tauri phase was simulated numerically with the purpose of studying the evolution of dust in the disk with distinct values of viscous α-parameter and dust fragmentation velocity v frag . Methods. We solved numerical hydrodynamics equations in the thin-disk limit, which are modified to include a dust component consisting of two parts: sub-micron-sized dust and grown dust with a maximum radius ar. The former is strictly coupled to the gas, while the latter interacts with the gas via friction. The conversion of small to grown dust, dust growth, and dust self-gravity are also considered. Results. We found that the process of dust growth known for the older protoplanetary phase also holds for the embedded phase of disk evolution. The dust growth efficiency depends on the radial distance from the star -ar is largest in the inner disk and gradually declines with radial distance. In the inner disk, ar is limited by the dust fragmentation barrier. The process of small-to-grown dust conversion is very fast once the disk is formed. The total mass of grown dust in the disk (beyond 1 AU) reaches tens or even hundreds of Earth masses already in the embedded phase of star formation and even a greater amount of grown dust drifts in the inner, unresolved 1 AU of the disk. Dust does not usually grow to radii greater than a few cm. A notable exception are models with α ≤ 10 −3 , in which case a zone with reduced mass transport develops in the inner disk and dust can grow to meter-sized boulders in the inner 10 AU. Grown dust drifts inward and accumulates in the inner disk regions. This effect is most pronounced in the α ≤ 10 −3 models where several hundreds of Earth masses can be accumulated in a narrow region of several AU from the star by the end of embedded phase. The efficiency of grown dust accumulation in spiral arms is stronger near corotation where the azimuthal velocity of dust grains is closest to the local velocity of the spiral pattern. In the framework of the adopted dust growth model, the efficiency of small-to-grown dust conversion was found to increase for lower values of α and v frag .

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. snyt@catalysis.ru 4 The commonly used term for this ratio, the "stopping time", can lead to confusion, since, in the approximation considered here, a particle does not acquire the velocity of the gas, but instead a constant velocity relative to the gas velocity.

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 tstop is less than or comparable to the rotation period, are considered. The methods are analyzed from the point of view of their suitability for computing the motions of small bodies, including dust grains less than 1 µm in size, 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 restrictions on the time step size τ due to gas drag in computations of the motion of grains of any size are presented. For the conditions of a circumstellar disk, the error in the velocity computations obtained when using some stable methods becomes unacceptably large when the time step size is τ > tstop. 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 computations 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.