[论文解读] Gaseous Dynamical Friction: a Numerical Study of Extended Perturbers

B. R. Morton, Sadegh Khochfar+1 more|arXiv (Cornell University)|Mar 29, 2021

Fluid Dynamics Simulations and Interactions参考文献 5被引用 5

一句话总结

本研究利用GIZMO代码中的高分辨率拉格朗日流体动力学模拟，探究了气态动力摩擦现象，发现现代流体动力学格式在马赫数0.7至2之间显著低估了阻力——低估幅度达10%至50%，尤其在M=1附近。该偏差源于马赫锥前缘密度剖面的数值弥散，导致尾迹结构失真并减少动量传递，从而对依赖这些求解器的宇宙学模拟的准确性构成挑战。

ABSTRACT

The process of momentum and energy transfer between a massive body and a background medium it is moving through is known as dynamical friction (DF). It 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 solver, the moving-mesh MUSCL scheme, and the piecewise parabolic method (PPM) solver. We use a set of simulations of massive bodies, modelled as Plummer spheres, moving with Mach $0.2 \leq \mathcal{M} \leq 3$. We investigate at which radial distances from the perturber these solvers recover the linear point mass solution for gaseous DF. We analyse the drag force and the structure and time evolution of the wake. The different solvers agree closely. Numerical convergence is reached when the initial spatial resolution is $0.2r_s$, where $r_s$ is the softening scale of the Plummer sphere. We find that the wake structure and drag force are recovered, at the $5\%$ level, when compared beyond $4r_\mathrm{s}$. Our results predict that models using the standard linear point mass DF solution will overestimate the drag force on extended perturbers by as much as 25\%, for Mach$\sim$1. Finally, we consider DF in the context of galaxy clusters, where dark matter subhaloes move through circumgalactic media. 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.

研究动机与目标

• 评估现代拉格朗日流体动力学求解器在不同马赫数下模拟气态动力摩擦的准确性。
• 识别数值模拟与解析预测之间在大质量扰动体穿过气态介质时阻力系数差异的根源。
• 评估阻力系数的数值偏差是否源于分辨率不足、尾迹结构不准确，或流体动力学格式本身的固有限制。
• 确定该阻力系数偏差对宇宙学模拟的影响，特别是对暗物质子结构并合率的影响。
• 倡导将动力摩擦作为流体动力学代码的标准基准测试，因其具有明确的解析解和物理相关性。

提出的方法

• 使用GIZMO代码中的拉格朗日流体动力学求解器（GSPH与G4PSPH）进行高分辨率模拟，以建模大质量扰动体在均匀气态介质中的运动。
• 模拟覆盖马赫数范围从0.2到3，通过改变引力软化长度和初始条件来测试结果的鲁棒性。
• 基于线性扰动理论推导动力摩擦力的解析预测，特别是Ostriker（1999）和Just & Kegel（1990）的解。
• 通过测量密度剖面和过密特征（尤其是马赫锥前缘）来分析尾迹的径向结构。
• 追踪阻力系数和尾迹形态的时间演化，以评估收敛性和长期行为。
• 通过比较不同半径和马赫数下模拟的阻力系数与解析预测，量化阻力系数的偏差。

实验结果

研究问题

• RQ1 现代拉格朗日流体动力学求解器在多大程度上能准确再现大质量扰动体在气态介质中运动时的解析阻力系数？
• RQ2 为何模拟的阻力系数系统性地低于解析预测，特别是在0.7 < M < 2的马赫数范围内？
• RQ3 马赫锥密度剖面的数值弥散如何影响整体动量传递和阻力计算？
• RQ4 总阻力系数中有多大比例来自远场尾迹（r ≥ 4rs）？其随时间如何演化？
• RQ5 该阻力系数偏差在多大程度上影响宇宙学模拟中暗物质子结构的并合率预测？

主要发现

• 在马赫数0.7至2之间，模拟阻力系数系统性地低于解析预测10%至50%，且在M=1附近偏差最大。
• 该偏差主要源于马赫锥前缘尖锐过密区密度剖面的数值弥散，导致峰值密度降低并向前扩散。
• 在最内层4rs（A=0.1时）或8rs（A=1时）区域，结构失配最为显著，尾迹的锐利过渡被平滑化。
• 远场尾迹（r ≥ 4rs）在数值上被良好再现，其对总阻力系数的贡献随时间增加，但仍不足以弥补内区的偏差。
• 该阻力系数偏差在至少150个穿越时间后依然存在，且在不同求解器（GSPH与G4PSPH）和初始条件下均具鲁棒性。
• 结果表明，使用这些求解器的宇宙学模拟可能系统性低估了暗物质子结构所经历的动力摩擦，从而导致并合率被低估和晕质量函数不准确。