Gaseous Dynamical Friction: a Numerical Study of Extended Perturbers

Ben Morton,Sadegh Khochfar,Jose Oñorbe
2024-09-12
Abstract:The process of momentum and energy transfer from a massive body moving through a background medium, known as dynamical friction (DF), is key to our understanding of many astrophysical systems. We present a series of high-resolution simulations of gaseous DF using Lagrangian meshless finite mass hydrodynamics solvers in the state-of-the-art multiphysics code Gizmo, the moving-mesh MUSCL scheme in Arepo, and the PPM solver in the Enzo adaptive mesh refinement code. We setup simulations of extended perturbers moving with Mach $0.2 \leq \mathcal{M} \leq 3$ and investigate at which radial distances from the perturber they recover the linear point-mass perturber solutions for the DF drag force, the radial structure of the wake, and their time evolution. The various different solvers agree with each other and we find that the wake structure is recovered, at the $5\%$ level, beyond $r_\mathrm{min}=4r_\mathrm{s}$, with $r_\mathrm{s}$ being the softening scale of the Plummer perturber. The drag force from the gravitationally induced wake is recovered, at this level, beyond this radius. Numerical convergence is reached where the effective initial spatial resolution is $0.2r_s$. This is roughly equivalent to 5 cells per gravitional softening scale in a grid-based code. The value of $r_\mathrm{min}$ provides a natural way to compare different codes and allows us to estimate the difference in DF force that will be experienced. Our results further predict that semianalytic models using the point-mass linear DF estimate from Ostriker (1999) will overestimate the DF forces on the extended perturber by as much as $25\%$, for Mach numbers close to 1. Finally, we show that DF is typically in the linear regime for most subhaloes in hosting haloes $<10^{11}$ M$_{\odot}$ but non-linear in more massive host haloes.
Astrophysics of Galaxies,Cosmology and Nongalactic Astrophysics,Instrumentation and Methods for Astrophysics
What problem does this paper attempt to address?