A solution to the overdamping problem when simulating dustgas mixtures with smoothed particle hydrodynamics
Abstract
We present a fix to the overdamping problem found by Laibe & Price when simulating strongly coupled dustgas mixtures using two different sets of particles using smoothed particle hydrodynamics. Our solution is to compute the drag at the barycentre between gas and dust particle pairs when computing the drag force by reconstructing the velocity field, similar to the procedure in Godunovtype solvers. This fixes the overdamping problem at negligible computational cost, but with additional memory required to store velocity derivatives. We employ slope limiters to avoid spurious oscillations at shocks, finding the van Leer Monotonized Central limiter most effective.
 Publication:

Monthly Notices of the Royal Astronomical Society
 Pub Date:
 July 2020
 DOI:
 10.1093/mnras/staa1366
 arXiv:
 arXiv:2005.06562
 Bibcode:
 2020MNRAS.495.3929P
 Keywords:

 hydrodynamics;
 methods: numerical;
 protoplanetary discs;
 dust;
 extinction;
 Astrophysics  Instrumentation and Methods for Astrophysics;
 Astrophysics  Earth and Planetary Astrophysics;
 Physics  Computational Physics
 EPrint:
 6 pages, 5 figures, accepted to MNRAS