Integration of Particle-Gas Systems with Stiff Mutual Drag Interaction

Abstract

Numerical simulation of numerous mm/cm-sized particles embedded in a gaseous disk has become an important tool in the study of planet formation and in understanding the dust distribution in observed protoplanetary disks. However, the mutual drag force between the gas and the particles can become so stiff, particularly because of small particles and/or strong local solid concentration, that an explicit integration of this system is computationally formidable. In this work, we consider the integration of the mutual drag force in a system of Eulerian gas and Lagrangian solid particles. Despite the entanglement between the gas and the particles under the particle-mesh construct, we are able to devise a numerical algorithm that effectively decomposes the globally coupled system of equations for the mutual drag force and makes it possible to integrate this system on a cell-by-cell basis, which considerably reduces the computational task required. We use an analytical solution for the temporal evolution of each cell to relieve the time-step constraint posed by the mutual drag force as well as to achieve the highest degree of accuracy. To validate our algorithm, we use an extensive suite of benchmarks with known solutions in one, two, and three dimensions, including the linear growth and the nonlinear saturation of the streaming instability. We demonstrate numerical convergence and satisfactory consistency in all cases. Our algorithm can for example be applied to model the evolution of the streaming instability with mm/cm-sized pebbles at high mass loading, which has important consequences for the formation scenarios of planetesimals.

Accepted for publication in the Astrophysical Journal Supplement Series. 21 pages, 15 figures. Fixed cross references for equations

References in corpus (12)

• Asymmetric features in the protoplanetary disk MWC758
• How to form planetesimals from mm-sized chondrules and chondrule aggregates
• Dust flow in gas disks in the presence of embedded planets
• Stellar population synthesis based modelling of the Milky Way using asteroseismology of 13000 Kepler red giants
• Standing on the shoulders of giants: Trojan Earths and vortex trapping in low mass self-gravitating protoplanetary disks of gas and solids
• Dust sedimentation and self-sustained Kelvin-Helmholtz turbulence in protoplanetary disk mid-planes. I. Radially symmetric simulations
• Modeling Dust Emission of HL Tau Disk Based on Planet-Disk Interactions
• A representative particle approach to coagulation and fragmentation of dust aggregates and fluid droplets
• Inferring Planet Mass from Spiral Structures in Protoplanetary Disks
• Turbulent Torques on Protoplanets in a Dead Zone
• Midplane sedimentation of large solid bodies in turbulent protoplanetary discs
• On the Feeding Zone of Planetesimal Formation by the Streaming Instability

Cited by in corpus (35)

• Concentrating small particles in protoplanetary disks through the streaming instability
• Thresholds for Particle Clumping by the Streaming Instability
• Asymptotically stable numerical method for multispecies momentum transfer: gas and multifluid-dust test suite and implementation in FARGO3D
• Diffusion and Concentration of Solids in the Dead Zone of a Protoplanetary Disk
• A thermodynamic view of dusty protoplanetary disks
• Pebble-driven planet formation for TRAPPIST-1 and other compact systems
• Dust settling against hydrodynamic turbulence in protoplanetary discs
• Growth after the streaming instability: from planetesimal accretion to pebble accretion
• Streaming Instability with Multiple Dust Species: I. Favourable Conditions for the Linear Growth
• On the coexistence of the streaming instability and the vertical shear instability in protoplanetary disks
• Grand challenges in protoplanetary disc modelling
• MULTIGRAIN: A smoothed particle hydrodynamics algorithm for multiple small dust grains and gas
• Streaming Instability with Multiple Dust Species: II. Turbulence and Dust-Gas Dynamics at Nonlinear Saturation
• A Multi-Fluid Dust Module in Athena++: Algorithms and Numerical Tests
• The Dynamics of Charged Dust in Magnetized Molecular Clouds
• Analysis of Numerical Algorithms for Computing Rapid Momentum Transfers between the Gas and Dust in Simulations of Circumstellar Disks
• A Particle Module for the PLUTO Code: III -- Dust
• Simulations of Dynamical Gas-Dust Circumstellar Disks: Going Beyond the Epstein Regime
• Formation of Dust Rings and Gaps in Non-ideal MHD Disks Through Meridional Gas Flows
• Two-Fluid Dusty Gas in Smoothed Particle Hydrodynamics: Fast and Implicit Algorithm for Stiff Linear Drag
• Dusty disc-planet interaction with dust-free simulations
• Direct Formation of Planetary Embryos in Self-Gravitating Disks
• Particle Dynamics in 3D Self-gravitating Disks I: Spirals
• A solution to the overdamping problem when simulating dust-gas mixtures with smoothed particle hydrodynamics
• Dust dynamics in RAMSES -- I. Methods and turbulent acceleration
• Morphological signatures induced by dust back reaction in discs with an embedded planet
• Fast method to simulate dynamics of two-phase medium with intense interaction between phases by smoothed particle hydrodynamics: gas-dust mixture with polydisperse particles, linear drag, one-dimensional tests
• Dust-vortex instability in the regime of well-coupled grains
• The Concentration and Growth of Solids in Fragmenting Circumstellar Disks
• Growth after the streaming instability: The radial distance dependence of the planetary growth
• A Staggered Semi-Analytic Method for Simulating Dust Grains Subject to Gas Drag
• Particle clustering in turbulence: Prediction of spatial and statistical properties with deep learning
• Developing a Non-Newtonian Fluid Model for Dust, for Application to Astrophysical Flows
• Large dust fractions can prevent the propagation of soundwaves
• Forming Earth-like and Low-Mass Rocky Exoplanets Through Pebble and Planetesimal Accretion