A Multi-Fluid Dust Module in Athena++: Algorithms and Numerical Tests

Abstract

We describe the algorithm, implementation and numerical tests of a multifluid dust module in the Athena++ magnetohydrodynamic (MHD) code. The module can accommodate an arbitrary number of dust species interacting with the gas via aerodynamic drag (characterized by the stopping time), with a number of numerical solvers. In particular, we describe two second-order accurate, two-stage, fully-implicit solvers that are stable in stiff regimes including short stopping time and high dust mass loading, and they are paired with the second-order explicit van-Leer and Runge-Kutta gas dynamics solvers in Athena++, respectively. Moreover, we formulate a consistent treatment of dust concentration diffusion with dust back-reaction, which incorporates momentum diffusion and ensures Galilean invariance. The new formulation and stiff drag solvers are implemented to be compatible with most existing features of Athena++, including different coordinate systems, mesh refinement, shearing-box and orbital advection. We present a large suite of test problems, including the streaming instability in linear and nonlinear regimes, as well as local and global setting, which demonstrate that the code achieves the desired performance. This module will be particularly useful for studies of dust dynamics and planet formation in protoplanetary disks.

arXiv comment32 pages, 14 figures, ApJS resubmitted after addressing referee's comments

References in corpus (83)

• Astropy: A Community Python Package for Astronomy
• Rapid planetesimal formation in turbulent circumstellar discs
• Streaming Instabilities in Protoplanetary Disks
• Rapid growth of gas-giant cores by pebble accretion
• Gas- and dust evolution in protoplanetary disks
• Particle Stirring in Turbulent Gas Disks: Including Orbital Oscillations
• The effect of gas drag on the growth of protoplanets -- Analytical expressions for the accretion of small bodies in laminar disks
• A Major Asymmetric Dust Trap in a Transition Disk
• FARGO: a fast eulerian transport algorithm for differentially rotating disks
• The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers
• Protoplanetary Disk Turbulence Driven by the Streaming Instability: Non-Linear Saturation and Particle Concentration
• Closed-form expressions for particle relative velocities induced by turbulence
• A comparative study of disc-planet interaction
• Phantom: A smoothed particle hydrodynamics and magnetohydrodynamics code for astrophysics
• Observations of Protoplanetary Disk Structures
• Gaps, Rings, and Non-Axisymmetric Structures in Protoplanetary Disks - From Simulations to ALMA Observations
• Dynamics of Solids in the Midplane of Protoplanetary Disks: Implications for Planetesimal Formation
• Evidence of fast pebble growth near condensation fronts in the HL Tau protoplanetary disk
• Protoplanetary Disk Turbulence Driven by the Streaming Instability: Linear Evolution and Numerical Methods
• Concentrating small particles in protoplanetary disks through the streaming instability
• Sintering-induced Dust Ring Formation in Protoplanetary Disks: Application to the HL Tau Disk
• FARGO3D: A new GPU-oriented MHD code
• How to form planetesimals from mm-sized chondrules and chondrule aggregates
• The great dichotomy of the Solar System: small terrestrial embryos and massive giant planet cores
• MPI-AMRVAC for Solar and Astrophysics
• MPI-AMRVAC 2.0 for Solar and Astrophysical Applications
• Dust Growth and Settling in Protoplanetary Disks and Disk Spectral Energy Distributions. I. Laminar Disks
• Multiple Disk Gaps and Rings Generated by a Single Super-Earth
• The gas temperature in the surface layers of protoplanetary disks
• Two-Component Secular Gravitational Instability in a Protoplanetary Disk: A Possible Mechanism for Creating Ring-Like Structures
• Vortex generation in protoplanetary disks with an embedded giant planet
• The Effect of the Radial Pressure Gradient in Protoplanetary Disks on Planetesimal Formation
• Particle-Gas Dynamics with Athena: Method and Convergence
• Thresholds for Particle Clumping by the Streaming Instability
• Particle Concentration At Planet Induced Gap Edges and Vortices: I. Inviscid 3-D Hydro Disks
• Implementation of the Shearing Box Approximation in Athena
• The co-orbital corotation torque in a viscous disk: numerical simulations
• High-resolution simulations of planetesimal formation in turbulent protoplanetary discs
• On the formation of multiple concentric rings and gaps in protoplanetary disks
• Secular Instability and Planetesimal Formation in the Dust Layer
• The formation of rings and gaps in magnetically coupled disk-wind systems: ambipolar diffusion and reconnection
• Asymptotically stable numerical method for multispecies momentum transfer: gas and multifluid-dust test suite and implementation in FARGO3D
• Potential Vorticity Evolution of a Protoplanetary Disk with An Embedded Protoplanet
• Dusty gas with SPH - I. Algorithm and test suite
• Type I Planet Migration in Nearly Laminar Disks
• Evidence for universality in the initial planetesimal mass function
• Effects of dust feedback on vortices in protoplanetary disks
• Dusty gas with one fluid
• Multiple Disk Gaps and Rings Generated by a Single Super-Earth: II. Spacings, Depths, and Number of Gaps, with Application to Real Systems
• Density Waves Excited by Low-Mass Planets in Protoplanetary Disks II: High-Resolution Simulations of the Nonlinear Regime
• Planet-driven spiral arms in protoplanetary disks: I. Formation mechanism
• Including Dust Coagulation in Hydrodynamic Models of Protoplanetary Disks: Dust Evolution in the Vicinity of a Jupiter-mass Planet
• Gas and dust hydrodynamical simulations of massive lopsided transition discs - II. Dust concentration
• An Origin of Multiple Ring Structure and Hidden Planets in HL Tau: A Unified Picture by Secular Gravitational Instability
• Density Waves Excited by Low-Mass Planets in Protoplanetary Disks I: Linear Regime
• Ring formation and dust dynamics in wind-driven protoplanetary discs: global simulations
• Dust Transport in MRI Turbulent Disks: Ideal and Non-ideal MHD with Ambipolar Diffusion
• An explicit scheme for multifluid magnetohydrodynamics
• Global Three-Dimensional Simulations of Outer Protoplanetary Disks with Ambipolar Diffusion
• Planet-driven spiral arms in protoplanetary disks: II. Implications
• Individual and collective behavior of dust particles in a protoplanetary nebula
• The Athena Astrophysical MHD Code in Cylindrical Geometry
• On the coexistence of the streaming instability and the vertical shear instability in protoplanetary disks
• Revised description of dust diffusion and a new instability creating multiple rings in protoplanetary disks
• The Gravitational instability of solids assisted by gas drag: slowing by turbulent mass diffusivity
• DUSTYBOX and DUSTYWAVE: Two test problems for numerical simulations of two fluid astrophysical dust-gas mixtures
• Ubiquitous Instabilities of Dust Moving in Magnetized Gas
• Integration of Particle-Gas Systems with Stiff Mutual Drag Interaction
• Streaming instability in the quasi-global protoplanetary discs
• Dust Settling and Clumping in MRI Turbulent Outer Protoplanetary Disks
• Snow-lines can be thermally unstable
• Pebble Accretion in Turbulent Protoplanetary Disks
• Planet-induced Vortices with Dust Coagulation in Protoplanetary Disks
• A Particle Module for the PLUTO Code: III -- Dust
• Meso-scale Instability Triggered by Dust Feedback in Dusty Rings: Origin and Observational Implications
• PIERNIK mhd code - a multi-fluid, non-ideal extension of the relaxing-TVD scheme (I)
• Effects of Ringed Structures and Dust Size Growth on Millimeter Observations of Protoplanetary Disks
• Vortex instabilities triggered by low-mass planets in pebble-rich, inviscid protoplanetary discs
• PIERNIK mhd code - a multi-fluid, non-ideal extension of the relaxing-TVD scheme (II)
• Morphological signatures induced by dust back reaction in discs with an embedded planet
• Migrating Low-Mass Planets in Inviscid Dusty Protoplanetary Discs
• Dust Transport in Protoplanetary Disks with Wind-driven Accretion
• A Staggered Semi-Analytic Method for Simulating Dust Grains Subject to Gas Drag

Cited by in corpus (20)

• MPI-AMRVAC 3.0: updates to an open-source simulation framework
• Three-dimensional Global Simulations of Type-II Planet-disk Interaction with a Magnetized Disk Wind: I. Magnetic Flux Concentration and Gap Properties
• Atmospheric Recyling of Volatiles by Pebble-Accreting Planets
• Three-temperature radiation hydrodynamics with PLUTO: Tests and applications to protoplanetary disks
• Chemical footprints of giant planet formation. Role of planet accretion in shaping the C/O ratio of protoplanetary disks
• Solving for the 2D Water Snowline with Hydrodynamic Simulations. Emergence of gas outflow, water cycle and temperature plateau
• TriPoD: Tri-Population size distributions for Dust evolution. Coagulation in vertically integrated hydrodynamic simulations of protoplanetary disks
• Three-Dimensional Dust Stirring by a Giant Planet Embedded in a Protoplanetary Disk
• Dust-gas dynamics driven by the streaming instability with various pressure gradients
• On Secular Gravitational Instability in Vertically Stratified Disks
• Particle clustering in turbulence: Prediction of spatial and statistical properties with deep learning
• Dynamics of small, constant size particles in a protoplanetary disk with an embedded protoplanet
• Filament formation due to diffusive instabilities in dusty protoplanetary disks
• A new multifluid method for dusty astrophysical flows. Application to turbulent protostellar collapses
• Evading the dust fragmentation barrier with the streaming instability in protoplanetary disks
• Hydrodynamical simulations of the vertical shear instability with dynamic dust and cooling rates in protoplanetary disks
• Nonlinear Evolution of the Unstratified Polydisperse Dust Settling Instability
• A multi-fluid approach for polydisperse pebble accretion: From particles to fluids, establishing the multifluid framework
• Planetesimal gravitational collapse in a gaseous environment: Thermal and dynamic evolution
• Modeling Multiphase Galactic Outflows: A Multifluid Moving Mesh Approach