Abstract
Abstract Implementation of a new particle module describing the physics of dust grains coupled to a 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 that 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.
Published Version (Free)
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.