Asymptotically stable numerical method for multispecies momentum transfer: gas and multifluid-dust test suite and implementation in FARGO3D

Abstract

We present an asymptotically and unconditionally stable numerical method to account for the momentum transfer between multiple species. Momentum is conserved to machine precision. This implies that the asymptotic equilibrium corresponds to the velocity of the center of mass. Aimed at studying dust dynamics, we implement this numerical method in the publicly available code FARGO3D. To validate our implementation, we develop a test suite for an arbitrary number of species, based on analytical or exact solutions of problems related to perfect damping, damped sound waves, shocks, local and global gas-dust radial drift in a disk and, linear streaming instability. In particular, we obtain first-order, steady-state solutions for the radial drift of multiple dust species in protoplanetary disks, in which the pressure gradient is not necessarily small. We additionally present non-linear shearing-box simulations of the streaming instability and compare them with previous results obtained with Lagrangian particles. We successfully validate our implementation by recovering the solutions from the test suite to second- and first-order accuracy in space and time, respectively. From this, we conclude that our scheme is suitable, and very robust, to study the self-consistent dynamics of several fluids. In particular, it can be used for solving the collisions between gas and dust in protoplanetary disks, with any degree of coupling.

arXiv comment27 pages, 12 figures; accepted for publication in ApJS

References in corpus (41)

• The NumPy array: a structure for efficient numerical computation
• Streaming Instabilities in Protoplanetary Disks
• Athena: A New Code for Astrophysical MHD
• Protoplanetary Disk Turbulence Driven by the Streaming Instability: Non-Linear Saturation and Particle Concentration
• A comparative study of disc-planet interaction
• Phantom: A smoothed particle hydrodynamics and magnetohydrodynamics code for astrophysics
• Protoplanetary Disk Turbulence Driven by the Streaming Instability: Linear Evolution and Numerical Methods
• FARGO3D: A new GPU-oriented MHD code
• MPI-AMRVAC for Solar and Astrophysics
• Hydromagnetic turbulence in computer simulations
• How Empty are Disk Gaps Opened by Giant Planets?
• Dust flow in gas disks in the presence of embedded planets
• Dust settling in local simulations of turbulent protoplanetary disks
• Particle-Gas Dynamics with Athena: Method and Convergence
• Steady state of dust distributions in disk vortices: Observational predictions and applications to transitional disks
• Simulations of dust-trapping vortices in protoplanetary discs
• Dust Trapping by Vortices in Transitional Disks: Evidence for Non-ideal MHD Effects in Protoplanetary Disks
• Dusty gas with SPH - I. Algorithm and test suite
• On the origin of horseshoes in transitional discs
• Characterizing the variable dust permeability of planet-induced gaps
• Gas and dust hydrodynamical simulations of massive lopsided transition discs - II. Dust concentration
• An opening criterion for dust gaps in protoplanetary discs
• The early evolution of viscous and self-gravitating circumstellar disks with a dust component
• Dust settling and rings in the outer regions of protoplanetary discs subject to ambipolar diffusion
• Dust and gas mixtures with multiple grain species - a one-fluid approach
• Gas and multi-species dust dynamics in viscous protoplanetary discs: the importance of the dust back-reaction
• Torques Induced by Scattered Pebble-flow in Protoplanetary Disks
• DUSTYBOX and DUSTYWAVE: Two test problems for numerical simulations of two fluid astrophysical dust-gas mixtures
• Dust Settling in Magnetorotationally-Driven Turbulent Discs I: Numerical Methods and Evidence for a Vigorous Streaming Instability
• MULTIGRAIN: A smoothed particle hydrodynamics algorithm for multiple small dust grains and gas
• Integration of Particle-Gas Systems with Stiff Mutual Drag Interaction
• Smoothed particle hydrodynamics simulations of gas and dust mixtures
• Streaming instability in the quasi-global protoplanetary discs
• The Wardle Instability in Interstellar Shocks: I. Nonlinear Dynamical Evolution
• Eccentricity excitation and merging of planetary embryos heated by pebble accretion
• A Hybrid Scheme for Gas-Dust Systems Stiffly Coupled via Viscous Drag
• Analysis of Numerical Algorithms for Computing Rapid Momentum Transfers between the Gas and Dust in Simulations of Circumstellar Disks
• Dust segregation in Hall-dominated turbulent protoplanetary disks
• Two-Fluid Dusty Gas in Smoothed Particle Hydrodynamics: Fast and Implicit Algorithm for Stiff Linear Drag
• Long-term and large-scale hydrodynamical simulations of migrating planets
• Dusty disc-planet interaction with dust-free simulations

Cited by in corpus (81)

• Streaming Instability for Particle-Size Distributions
• Streaming Instability in Turbulent Protoplanetary Disks
• Turbulence Regulates the Rate of Planetesimal Formation via Gravitational Collapse
• Probing the Protosolar Disk Using Dust Filtering at Gaps in the Early Solar System
• Streaming Instability with Multiple Dust Species: I. Favourable Conditions for the Linear Growth
• Polydisperse Streaming Instability I. Tightly coupled particles and the terminal velocity approximation
• Predicting the observational signature of migrating Neptune-sized planets in low-viscosity disks
• Protostellar collapse simulations in spherical geometry with dust coagulation and fragmentation
• 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
• Puffed up Edges of Planet-opened Gaps in Protoplanetary Disks. I. hydrodynamic simulations
• Stratified and vertically-shearing streaming instabilities in protoplanetary disks
• A Particle Module for the PLUTO Code: III -- Dust
• Simulations of Dynamical Gas-Dust Circumstellar Disks: Going Beyond the Epstein Regime
• Which planets trigger longer-lived vortices: low-mass or high-mass?
• Impact of Local Pressure Enhancements on Dust Concentration inTurbulent Protoplanetary Discs
• Aiolos- A multi-purpose 1-D hydrodynamics code for planetary atmospheres
• Dust entrainment in magnetically and thermally driven disk winds
• Polydisperse Streaming Instability II. Methods for solving the linear stability problem
• Dust Settling Instability in Protoplanetary Discs
• Three-temperature radiation hydrodynamics with PLUTO: Tests and applications to protoplanetary disks
• Modeling the nonaxisymmetric structure in the HD 163296 disk with planet-disk interaction
• Modeling multi-phase gases in cosmological simulations using compressible multi-fluid hydrodynamics
• The 3D dust and opacity distribution of protoplanets in multi-fluid global simulations
• New growth mechanism of dust grains in protoplanetary disks with magnetically driven disk winds
• Vortex instabilities triggered by low-mass planets in pebble-rich, inviscid protoplanetary discs
• Dusty circumbinary discs: inner cavity structures and stopping locations of migrating planets
• A solution to the overdamping problem when simulating dust-gas mixtures with smoothed particle hydrodynamics
• Dust formation in the outflows of catastrophically evaporating planets
• On the non-axisymmetric fragmentation of rings generated by the Secular Gravitational Instability
• A gap-sharing planet pair shaping the crescent in HD 163296: a disk sculpted by a resonant chain
• Migrating Low-Mass Planets in Inviscid Dusty Protoplanetary Discs
• Streaming instabilities in accreting and magnetized laminar protoplanetary disks
• Pebble-driven migration of low-mass planets in the 2D regime of pebble accretion
• 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
• Gap opening by planets in discs with magnetised winds
• Gap-opening Planets Make Dust Rings Wider
• DBNets: A publicly available deep learning tool to measure the masses of young planets in dusty protoplanetary discs
• TriPoD: Tri-Population size distributions for Dust evolution. Coagulation in vertically integrated hydrodynamic simulations of protoplanetary disks
• Dynamics of dust grains in turbulent molecular clouds. Conditions for decoupling and limits of different numerical implementations
• DPNNet-2.0 Part I: Finding hidden planets from simulated images of protoplanetary disk gaps
• Distinguishing Magnetized Disc Winds from Turbulent Viscosity through Substructure Morphology in Planet-forming Discs
• A Machine Learning model to infer planet masses from gaps observed in protoplanetary disks
• Kinematical Constraint on Eccentricity in the Protoplanetary Disk MWC 758 with ALMA
• Dust-gas dynamics driven by the streaming instability with various pressure gradients
• The steady-state hydrodynamics of a long-lived disc: planetary system architecture and prospects of observing a circumplanetary disc shadow in V4046Sgr
• Self-Sustaining Vortices in Protoplanetary Disks: Setting the Stage for Planetary System Formation
• Polydisperse Formation of Planetesimals: The dust size distribution in clumps
• Using Bayesian Deep Learning to infer Planet Mass from Gaps in Protoplanetary Disks
• How to form compact & other longer-lived planet-induced vortices: VSI, planet migration, or re-triggers, but not feedback
• Nonlinear evolution of streaming instabilities in accreting protoplanetary disks
• Distributions of gas and small and large grains in the LkH$\alpha 330$ disk trace a young planetary system
• Spreading pressure bumps in gas-dust discs can stall planet migration via planet-vortex interactions
• Using planet migration and dust drift to weigh protoplanetary discs
• DBNets2.0: simulation-based inference for planet-induced dust substructures in protoplanetary discs
• A smoothed particle hydrodynamics algorithm for multigrain dust with separate sets of particles
• Substructures in protoplanetary disks imprinted by compact planetary systems
• A Staggered Semi-Analytic Method for Simulating Dust Grains Subject to Gas Drag
• Size-selective accretion of dust onto CPDs: Low CPD masses and filtration of larger grains
• Stability of Dusty Rings in Protoplanetary Discs
• Strong clumping in global streaming instability simulations with a dusty fluid
• Resurgence of CO in a warm bubble around accreting protoplanets and its observability
• Planetary nurseries: vortices formed at smooth viscosity transition
• Dust-void evolution driven by turbulent dust flux can induce runaway migration of Earth-mass planets
• Dynamics of small, constant size particles in a protoplanetary disk with an embedded protoplanet
• A new multifluid method for dusty astrophysical flows. Application to turbulent protostellar collapses
• A multiwavelength analysis of the spiral arms in the protoplanetary disk around WaOph 6
• Dust back-reaction on gas around planets modifies the cold thermal torque
• Enhanced Pebble Drift Across Planet-Opened Gaps in Windy Protoplanetary Disks
• Single Fluid vs. Multifluid: Comparison between single fluid and multifluid dust models for disc planet interactions
• Large dust fractions can prevent the propagation of soundwaves
• Effects of turbulent diffusion and back-reaction on the dust distribution around two resonant planets
• Nonlinear Evolution of the Unstratified Polydisperse Dust Settling Instability
• Evading the dust fragmentation barrier with the streaming instability in protoplanetary disks
• Magnetic clumping of charged dust in the dense interstellar medium
• Streaming instabilities in accreting protoplanetary disks: A parameter study
• A multi-fluid approach for polydisperse pebble accretion: From particles to fluids, establishing the multifluid framework
• Turbulence destroys thermal lobes around Mars-sized planetary embryos
• Not Just Gas: How Solid-Driven Torques Shaped the Migration of the Galilean Moons
• Constraining the properties of the potential embedded planets in the disk around HD 100546
• Modeling Multiphase Galactic Outflows: A Multifluid Moving Mesh Approach