A Particle Module for the PLUTO Code: III -- Dust
A. Mignone, M. Flock, B. Vaidya

TL;DR
This paper introduces a new dust particle module for the PLUTO code, employing a stable, second-order accurate numerical scheme suitable for various geometries, enhancing modeling of dust-gas interactions in astrophysical disks.
Contribution
It presents a novel particle-gas hybrid scheme with an exponential midpoint integrator, improving stability and accuracy in dust-gas coupling simulations across multiple coordinate systems.
Findings
The scheme is second-order accurate in time and space.
It remains stable for small particle stopping times.
Numerical benchmarks demonstrate robustness and accuracy.
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…
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
