Modeling dust dynamics in OpenGadget3 -- I. SPH implementation of the One-Fluid model

TL;DR

This paper implements and tests an improved One-Fluid SPH model for dust-gas dynamics in astrophysics within OpenGadget3, enhancing stability, accuracy, and applicability across various complex flow scenarios.

Contribution

The study introduces specific improvements to the One-Fluid model in OpenGadget3, including time-dependent viscosity, conductivity, and a new dust diffusion treatment, validated through comprehensive benchmarks.

Findings

Model passes all benchmark tests successfully.

Enhanced diffusion term improves handling of high dust-to-gas ratios.

Model remains stable and accurate in complex astrophysical flows.

Abstract

Dust dynamics plays a critical role in astrophysical processes and has been modeled in hydrodynamical simulations using various approaches. Among particle-based methods like Smoothed Particle Hydrodynamics (SPH), the One-Fluid model has proven to be highly effective for simulating gas-dust mixtures. This study presents the implementation of the One-Fluid model in OpenGadget3, introducing improvements to the original formulation. These enhancements include time-dependent artificial viscosity and conductivity, as well as a novel treatment of dust diffusion using a pressure-like term. The improved model is tested using a suite of dust dynamics benchmark problems: DUSTYBOX, DUSTYWAVE, and DUSTYSHOCK, with the latter extended to multidimensional scenarios, as well as a dusty Sedov-Taylor blast wave. Additional tests include simulations of Cold Keplerian Disks, dusty protoplanetary disks, and Kelvin-Helmholtz instabilities to evaluate the model's robustness in more complex flows. The implementation successfully passes all standard benchmark tests. It demonstrates stability and accuracy in both simple and complex simulations. The new diffusion term improves the handling of flows with large dust-to-gas ratios and low drag coefficients, although limitations of the One-Fluid model in these regimes remain. The enhanced One-Fluid model is a reliable and robust tool for simulating dust dynamics in OpenGadget3. While it retains some limitations inherent to the original formulation, the introduced improvements expand its applicability and address some challenges in gas-dust dynamics.