A Particle Module for the PLUTO Code: III -- Dust

Abstract

The implementation of a new particle module describing the physics of dust grains coupled to the gas via drag forces is the subject of this work. The proposed particle-gas hybrid scheme has been designed to work in Cartesian as well as in cylindrical and spherical geometries. The numerical method relies on a Godunov-type second-order scheme for the fluid and an exponential midpoint rule for dust particles which overcomes the stiffness introduced by the linear coupling term. Besides being time-reversible and globally second-order accurate in time, the exponential integrator provides energy errors which are always bounded and it remains stable in the limit of arbitrarily small particle stopping times yielding the correct asymptotic solution. Such properties make this method preferable to the more widely used semi-implicit or fully implicit schemes at a very modest increase in computational cost. Coupling between particles and grid quantities is achieved through particle deposition and field-weighting techniques borrowed from Particle-In-Cell simulation methods. In this respect, we derive new weight factors in curvilinear coordinates that are more accurate than traditional volume- or area-weighting. A comprehensive suite of numerical benchmarks is presented to assess the accuracy and robustness of the algorithm in Cartesian, cylindrical and spherical coordinates. Particular attention is devoted to the streaming instability which is analyzed in both local and global disk models. The module is part of the PLUTO code for astrophysical gas-dynamics and it is mainly intended for the numerical modeling of protoplanetary disks in which solid and gas interact via aerodynamic drag.

arXiv comment22 pages, 13 figures

References in corpus (34)

• PLUTO: a Numerical Code for Computational Astrophysics
• Streaming Instabilities in Protoplanetary Disks
• First Results from High Angular Resolution ALMA Observations Toward the HL Tau Region
• The Disk Substructures at High Angular Resolution Project (DSHARP): I. Motivation, Sample, Calibration, and Overview
• The PLUTO Code for Adaptive Mesh Computations in Astrophysical Fluid Dynamics
• Protoplanetary Disk Turbulence Driven by the Streaming Instability: Non-Linear Saturation and Particle Concentration
• Linear and nonlinear evolution of the vertical shear instability in accretion discs
• Dust and Gas in the disc of HL Tauri: Surface density, dust settling, and dust-to-gas ratio
• Disks ARound TTauri Stars with Sphere (DARTTS-S) I: Sphere / IRDIS Polarimetric Imaging of 8 prominent TTauri Disks
• Protoplanetary Disk Turbulence Driven by the Streaming Instability: Linear Evolution and Numerical Methods
• MPI-AMRVAC for Solar and Astrophysics
• Dust Migration and Morphology in Optically Thin Circumstellar Gas Disks
• Dust diffusion in protoplanetary discs by magnetorotational turbulence
• Particle-Gas Dynamics with Athena: Method and Convergence
• Particle Concentration At Planet Induced Gap Edges and Vortices: I. Inviscid 3-D Hydro Disks
• High-order conservative reconstruction schemes for finite volume methods in cylindrical and spherical coordinates
• Asymptotically stable numerical method for multispecies momentum transfer: gas and multifluid-dust test suite and implementation in FARGO3D
• Dusty gas with one fluid
• Dust-trapping Rossby vortices in protoplanetary disks
• The VLA view of the HL Tau Disk - Disk Mass, Grain Evolution, and Early Planet Formation
• A Particle Module for the PLUTO code: II - Hybrid Framework for Modeling Non-thermal emission from Relativistic Magnetized flows
• Particle dynamics in discs with turbulence generated by the vertical shear instability
• Computing the dust distribution in the bowshock of a fast moving, evolved star
• A conservative orbital advection scheme for simulations of magnetized shear flows with the PLUTO code
• Models of the circumstellar medium of evolving, massive runaway stars moving through the Galactic plane
• A Particle Module for the PLUTO Code: I - an implementation of the MHD-PIC equations
• Bow shock nebulae of hot massive stars in a magnetized medium
• Dust Settling in Magnetorotationally-Driven Turbulent Discs I: Numerical Methods and Evidence for a Vigorous Streaming Instability
• Integration of Particle-Gas Systems with Stiff Mutual Drag Interaction
• Streaming instability in the quasi-global protoplanetary discs
• Particle accretion onto planets in discs with hydrodynamic turbulence
• The properties of the inner disk around HL Tau: Multi-wavelength modeling of the dust emission
• A Hybrid Scheme for Gas-Dust Systems Stiffly Coupled via Viscous Drag
• Dynamical friction for supersonic motion in a homogeneous gaseous medium

Cited by in corpus (17)

• Thresholds for Particle Clumping by the Streaming Instability
• MPI-AMRVAC 3.0: updates to an open-source simulation framework
• 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
• A Multi-Fluid Dust Module in Athena++: Algorithms and Numerical Tests
• The coexistence of the streaming instability and the vertical shear instability in protoplanetary disks: Planetesimal formation thresholds explored in two-dimensional global models
• Pebble-driven migration of low-mass planets in the 2D regime of pebble accretion
• Mixing is easy: New insights for cosmochemical evolution from pre-stellar core collapse
• Streaming instability in a global patch simulation of protoplanetary disks
• Inhibited destruction of dust by supernova in a clumpy medium
• Formation of dust filaments in the diffuse envelopes of molecular clouds
• PyPLUTO: a data analysis Python package for the PLUTO code
• Dynamics of small, constant size particles in a protoplanetary disk with an embedded protoplanet
• Tracking Dusty Cloud Crushed by a Hot Flow
• Developing a Non-Newtonian Fluid Model for Dust, for Application to Astrophysical Flows
• Infrared and X-ray emission of a supernova remnant in a clumpy medium
• Gas dynamics around dust asymmetries in turbulent disks