A solution to the overdamping problem when simulating dust-gas mixtures with smoothed particle hydrodynamics

Abstract

We present a fix to the overdamping problem found by Laibe & Price (2012) when simulating strongly coupled dust-gas 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 Godunov-type 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.

6 pages, 5 figures, accepted to MNRAS

References in corpus (25)

• SPLASH: An interactive visualisation tool for Smoothed Particle Hydrodynamics simulations
• Smoothed Particle Hydrodynamics and Magnetohydrodynamics
• Phantom: A smoothed particle hydrodynamics and magnetohydrodynamics code for astrophysics
• Inviscid SPH
• Particle-Gas Dynamics with Athena: Method and Convergence
• Asymptotically stable numerical method for multispecies momentum transfer: gas and multifluid-dust test suite and implementation in FARGO3D
• Dusty gas with SPH - I. Algorithm and test suite
• Dusty gas with one fluid
• A comparison between grid and particle methods on the statistics of driven, supersonic, isothermal turbulence
• A fast and explicit algorithm for simulating the dynamics of small dust grains with smoothed particle hydrodynamics
• Dusty gas with SPH - II. Implicit timestepping and astrophysical drag regimes
• DUSTYBOX and DUSTYWAVE: Two test problems for numerical simulations of two fluid astrophysical dust-gas mixtures
• Dusty gas with one fluid in smoothed particle hydrodynamics
• Two fluid dust and gas mixtures in SPH: A semi-implicit approach
• MULTIGRAIN: A smoothed particle hydrodynamics algorithm for multiple small dust grains and gas
• Integration of Particle-Gas Systems with Stiff Mutual Drag Interaction
• Small dust grain dynamics on adaptive mesh-refinement grids. I. Methods
• Enforcing dust mass conservation in 3D simulations of tightly-coupled grains with the Phantom SPH code
• A Hybrid Scheme for Gas-Dust Systems Stiffly Coupled via Viscous Drag
• Two-fluid dust and gas mixtures in smoothed particle hydrodynamics II: an improved semi-implicit approach
• Accurate tracer particles of baryon dynamics in the adaptive mesh refinement code Ramses
• Two-Fluid Dusty Gas in Smoothed Particle Hydrodynamics: Fast and Implicit Algorithm for Stiff Linear Drag
• Super-Earths in the TW Hya disc
• The impact of numerical oversteepening on the fragmentation boundary in self-gravitating disks
• A Staggered Semi-Analytic Method for Simulating Dust Grains Subject to Gas Drag

Cited by in corpus (7)

• Dust growth, fragmentation and self-induced dust traps in PHANTOM
• The role of the drag force in the gravitational stability of dusty planet-forming disc -- II. Numerical simulations
• Dust Traffic Jams in Inclined Circumbinary Protoplanetary Discs I. Morphology and Formation Theory
• Fast method to simulate dynamics of two-phase medium with intense interaction between phases by smoothed particle hydrodynamics: gas-dust mixture with polydisperse particles, linear drag, one-dimensional tests
• The protoplanetary disc around HD 169142: circumstellar or circumbinary?
• A smoothed particle hydrodynamics algorithm for multigrain dust with separate sets of particles
• Dust trapping in protoplanetary discs after stellar flybys