Two-Fluid Dusty Gas in Smoothed Particle Hydrodynamics: Fast and Implicit Algorithm for Stiff Linear Drag Full article
| Journal |
Astronomy and Computing
ISSN: 2213-1337 |
||||||||
|---|---|---|---|---|---|---|---|---|---|
| Output data | Year: 2018, Volume: 25, Pages: 25-37 Pages count : 13 DOI: 10.1016/j.ascom.2018.08.004 | ||||||||
| Tags | Protoplanetary discs; Hydrodynamics; Epstein drag simulation | ||||||||
| Authors |
|
||||||||
| Affiliations |
|
Funding (1)
| 1 | Ministry of Science and Higher Education of the Russian Federation | ГЗ-2017-2020 |
Abstract:
Simulation of the dynamics of dust–gas circumstellar discs is crucial in understanding the mechanisms ofplanet formation. The dynamics of small grains in the disc is stiffly coupled to the gas, while the dynamicsof 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 placehigh requirements on the numerical methods for disc simulations. In particular, when gas and dust aresimulated with two different fluids, explicit methods require very small timestep (must be less than duststopping time tstopduring which the velocity of a solid particle is equalized with respect to the gas velocity)to obtain solution, while some implicit methods require high temporal resolution to obtain acceptableaccuracy. Moreover, recent studies underlined that for Smoothed particle hydrodynamics (SPH) whenthe gas and the dust are simulated with different sets of particles only high spatial resolution h < cststopguarantees 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 integrationof linear (Epstein) drag and (2) exact conservation of local linear momentum. We derived formulas formonodisperse dust–gas in two-fluid SPH and tested the new method on problems with known analyticalsolutions. We found that our method is a promising alternative for the previously developed two-fluidSPH scheme in case of stiff linear drag thanks to the fact that spatial resolution condition h < cststopis notrequired anymore for accurate results.
Cite:
Stoyanovskaya O.P.
, Glushko T.A.
, Snytnikov N.V.
, Snytnikov V.N.
Two-Fluid Dusty Gas in Smoothed Particle Hydrodynamics: Fast and Implicit Algorithm for Stiff Linear Drag
Astronomy and Computing. 2018.
V.25. P.25-37. DOI: 10.1016/j.ascom.2018.08.004
WOS
Scopus
РИНЦ
OpenAlex
Two-Fluid Dusty Gas in Smoothed Particle Hydrodynamics: Fast and Implicit Algorithm for Stiff Linear Drag
Astronomy and Computing. 2018.
V.25. P.25-37. DOI: 10.1016/j.ascom.2018.08.004
WOS
Scopus
РИНЦ
OpenAlex
Files:
Full text from publisher
Dates:
| Submitted: | May 31, 2018 |
| Accepted: | Aug 13, 2018 |
| Published online: | Aug 29, 2018 |
| Published print: | Oct 1, 2018 |
Identifiers:
| Web of science: | WOS:000451390700003 |
| Scopus: | 2-s2.0-85053066025 |
| Elibrary: | 35726002 |
| OpenAlex: | W3106448031 |