[论文解读] Sphenix: Smoothed Particle Hydrodynamics for the next generation of galaxy formation simulations

Josh Borrow, Matthieu Schaller+2 more|Durham Research Online (Durham University)|Dec 7, 2020

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

一句话总结

Sphenix 是一种用于星系形成模拟的下一代光滑粒子流体动力学（SPH）方案，采用密度-能量运动方程，具备可变人工黏性和导热系数，并引入一种新型限制器以抑制激波处的导热。该方案在关键流体动力学测试中实现了稳定且收敛的结果，已集成至 Swift 代码中并可投入生产使用，能够在保持低计算成本的同时实现精确的反馈建模。

ABSTRACT

Smoothed Particle Hydrodynamics (SPH) is a ubiquitous numerical method for solving the fluid equations, and is prized for its conservation properties, natural adaptivity, and simplicity. We introduce the Sphenix SPH scheme, which was designed with three key goals in mind: to work well with sub-grid physics modules that inject energy, be highly computationally efficient (both in terms of compute and memory), and to be Lagrangian. Sphenix uses a Density-Energy equation of motion, along with variable artificial viscosity and conduction, including limiters designed to work with common sub-grid models of galaxy formation. In particular, we present and test a novel limiter that prevents conduction across shocks, preventing spurious radiative losses in feedback events. Sphenix is shown to solve many difficult test problems for traditional SPH, including fluid mixing and vorticity conservation, and it is shown to produce convergent behaviour in all tests where this is appropriate. Crucially, we use the same parameters within Sphenix for the various switches throughout, to demonstrate the performance of the scheme as it would be used in production simulations.

研究动机与目标

• 开发一种SPH方案，使其能与注入能量至星际介质的亚网格反馈模型稳健耦合。
• 确保在大规模宇宙学模拟中具有高计算效率，兼顾内存与计算资源使用。
• 保持完全的拉格朗日框架，以实现与N体模拟中引力求解器的无缝耦合。
• 解决长期存在的SPH问题，如人工表面张力、涡度守恒性差以及虚假热传导。
• 开发一种稳定、收敛且经过最少调优即可在多种天体物理测试问题中投入生产的方案。

提出的方法

• 该方案采用SPH的密度-能量形式，以粒子密度和内能作为主要求解变量。
• 使用Balsara（1989）开关实现可变人工黏性，并引入限制器以防止在平稳流动中过度激活。
• 引入人工导热以减少接触间断面处的误差并促进相混合，同时采用一种新型限制器以抑制激波处的导热。
• 人工黏性和导热系数通过Sod激波管测试标定，并在所有其他测试中保持固定。
• 该方案已集成至Swift宇宙学模拟代码中，采用FMM引力求解器，并以开源方式发布。
• 采用改进的平滑核函数，扩展了‘翼部’以提升梯度精度并减少数值噪声。

实验结果

研究问题

• RQ1 单一固定SPH参数集是否能在多种流体动力学测试问题中实现稳定且收敛的结果？
• RQ2 在反馈事件中，引入具有激波阻断限制器的人工导热是否能防止虚假能量损失？
• RQ3 与标准SPH方案相比，Sphenix在捕捉激波、涡旋和接触间断面方面的表现如何？
• RQ4 Sphenix在复杂多相流（如nIFTy星系团和Evrard坍缩）中是否能保持精度与收敛性？
• RQ5 人工导热在星系团模拟中对相混合与熵核心形成有多大改善作用？

主要发现

• Sphenix在全部七个流体动力学测试中均表现出收敛行为，包括Sod激波管、Gresho涡旋和Kelvin-Helmholtz不稳定性测试。
• 新型导热限制器成功防止了激波处的虚假能量损失，实现了稳定反馈注入而无非物理解冷却现象。
• 该方案在固定参数下实现了比Cullen & Dehnen（2010）方案更快的收敛速度，且能稳定演化涡旋。
• 在nIFTy星系团测试中，Sphenix由于人工导热有效促进相混合，成功生成了具有核心结构的熵剖面，符合物理预期。
• Evrard坍缩测试表明，Sphenix能准确保持激波结构与能量守恒，无虚假加热或冷却现象。
• Sphenix现已成为Swift中的默认SPH方案，具备开源可用性与极低的每粒子计算成本，可支持可扩展的星系形成模拟。