Analysis of Numerical Algorithms for Computing Rapid Momentum Transfers between the Gas and Dust in Simulations of Circumstellar Disks Full article
Journal |
Astronomy Reports
ISSN: 1063-7729 , E-ISSN: 1562-6881 |
||||||||
---|---|---|---|---|---|---|---|---|---|
Output data | Year: 2018, Volume: 62, Number: 7, Pages: 455-468 Pages count : 14 DOI: 10.1134/s1063772918060069 | ||||||||
Tags | SMOOTHED PARTICLE HYDRODYNAMICS; PROTOPLANETARY DISKS; GRAVITATIONAL-INSTABILITY; SEMIIMPLICIT APPROACH; SOLID PARTICLES; 2-FLUID DUST; MIXTURES; VORTICES; DYNAMICS; EVOLUTION | ||||||||
Authors |
|
||||||||
Affiliations |
|
Funding (1)
1 | Russian Science Foundation | 17-12-01168 |
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 the methods used to simulate the dust dynamics. A grid-based, 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.
Cite:
Stoyanovskaya O.P.
, Vorobyov E.I.
, Snytnikov V.N.
Analysis of Numerical Algorithms for Computing Rapid Momentum Transfers between the Gas and Dust in Simulations of Circumstellar Disks
Astronomy Reports. 2018. V.62. N7. P.455-468. DOI: 10.1134/s1063772918060069 WOS Scopus РИНЦ OpenAlex
Analysis of Numerical Algorithms for Computing Rapid Momentum Transfers between the Gas and Dust in Simulations of Circumstellar Disks
Astronomy Reports. 2018. V.62. N7. P.455-468. DOI: 10.1134/s1063772918060069 WOS Scopus РИНЦ OpenAlex
Original:
Стояновская О.П.
, Воробьев Э.И.
, Снытников В.Н.
Анализ численных алгоритмов расчета быстрого обмена импульсом между газом и пылью при моделировании околозвездных дисков
Астрономический журнал. 2018. Т.95. №7. С.487-500. DOI: 10.1134/S0004629918060063 РИНЦ OpenAlex
Анализ численных алгоритмов расчета быстрого обмена импульсом между газом и пылью при моделировании околозвездных дисков
Астрономический журнал. 2018. Т.95. №7. С.487-500. DOI: 10.1134/S0004629918060063 РИНЦ OpenAlex
Files:
Full text from publisher
Dates:
Submitted: | Dec 29, 2017 |
Accepted: | Jan 26, 2018 |
Published print: | Jul 1, 2018 |
Published online: | Jul 5, 2018 |
Identifiers:
Web of science: | WOS:000437442300005 |
Scopus: | 2-s2.0-85049586134 |
Elibrary: | 35774756 |
OpenAlex: | W3104756557 |