English

A Staggered Semi-Analytic Method for Simulating Dust Grains Subject to Gas Drag

Earth and Planetary Astrophysics 2020-01-08 v2 Astrophysics of Galaxies Computational Physics

Abstract

Numerical simulations of dust-gas dynamics are one of the fundamental tools in astrophysical research, such as the study of star and planet formation. It is common to find tightly coupled dust and gas in astrophysical systems, which demands that any practical integration method be able to take time steps Δt\Delta t much longer than the stopping time tst_{\rm s} due to drag. A number of methods have been developed to ensure stability in this stiff (Δtts\Delta t\gg t_{\rm s}) regime, but there remains large room for improvement in terms of accuracy. In this paper, we describe an easy-to-implement method, the "staggered semi-analytic method" (SSA), and conduct numerical tests to compare it to other implicit and semi-analytic methods, including the 2nd2^{\rm nd} order implicit method and the Verlet method. SSA makes use of a staggered step to better approximate the terminal velocity in the stiff regime. In applications to protoplanetary disks, this not only leads to orders-of-magnitude higher accuracy than the other methods, but also provides greater stability, making it possible to take time steps 100 times larger in some situations. SSA is also 2nd2^{\rm nd} order accurate and symplectic when Δtts\Delta t \ll t_{\rm s}. More generally, the robustness of SSA makes it applicable to linear dust-gas drag in virtually any context.

Keywords

Cite

@article{arxiv.1909.02006,
  title  = {A Staggered Semi-Analytic Method for Simulating Dust Grains Subject to Gas Drag},
  author = {Jeffrey Fung and Dhruv Muley},
  journal= {arXiv preprint arXiv:1909.02006},
  year   = {2020}
}

Comments

Accepted to ApJ Supplement

R2 v1 2026-06-23T11:05:47.663Z