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

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 $\Delta t$ much longer than the stopping time $t_{\mathrm{s}}$ due to drag. A number of methods have been developed to ensure stability in this stiff ($\Delta t\gg t_{\mathrm{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 $2^{\mathrm{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 $2^{\mathrm{nd}}$ order accurate and symplectic when $\Delta t\ll t_{\mathrm{s}}$. More generally, the robustness of SSA makes it applicable to linear dust-gas drag in virtually any context.

arXiv commentAccepted to ApJ Supplement

References in corpus (29)

• Dust Filtration by Planet-Induced Gap Edges: Implications for Transitional Disks
• Dust flow in gas disks in the presence of embedded planets
• Multiple Disk Gaps and Rings Generated by a Single Super-Earth
• 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
• Can dead zones create structures like a transition disk?
• Dust Density Distribution and Imaging Analysis of Different Ice Lines in Protoplanetary Disks
• Asymptotically stable numerical method for multispecies momentum transfer: gas and multifluid-dust test suite and implementation in FARGO3D
• Dusty gas with SPH - I. Algorithm and test suite
• Effects of dust feedback on vortices in protoplanetary disks
• Multiple Disk Gaps and Rings Generated by a Single Super-Earth: II. Spacings, Depths, and Number of Gaps, with Application to Real Systems
• Gas and dust hydrodynamical simulations of massive lopsided transition discs - II. Dust concentration
• Streaming Instability for Particle-Size Distributions
• Dust distribution in protoplanetary disks - Vertical settling and radial migration
• A thermodynamic view of dusty protoplanetary disks
• Dusty gas with SPH - II. Implicit timestepping and astrophysical drag regimes
• Dust and gas mixtures with multiple grain species - a one-fluid approach
• On the Numerical Robustness of the Streaming Instability: Particle Concentration and Gas Dynamics in Protoplanetary Disks
• Gas and dust hydrodynamical simulations of massive lopsided transition discs - I. Gas distribution
• Dust capture and long-lived density enhancements triggered by vortices in 2D protoplanetary disks
• Pebble dynamics and accretion onto rocky planets. I. Adiabatic and convective models
• Linear growth of streaming instability in pressure bumps
• An analytical model of radial dust trapping in protoplanetary disks
• Integration of Particle-Gas Systems with Stiff Mutual Drag Interaction
• A Hybrid Scheme for Gas-Dust Systems Stiffly Coupled via Viscous Drag
• Two-fluid dust and gas mixtures in smoothed particle hydrodynamics II: an improved semi-implicit approach
• Pebble Accretion in Turbulent Protoplanetary Disks
• Effects of photophoresis on the dust distribution in a 3D protoplanetary disc
• The effect of radiation pressure on spatial distribution of dust inside HII regions

Cited by in corpus (3)

• A Multi-Fluid Dust Module in Athena++: Algorithms and Numerical Tests
• A solution to the overdamping problem when simulating dust-gas mixtures with smoothed particle hydrodynamics
• Dust Dynamics in Transitional Disks: Clumping and Disk Recession