[論文レビュー] 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 は、銀河形成シミュレーションを想定した次世代の Smoothed Particle Hydrodynamics (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の渦、ケルビン＝ヘルムホルツ不安定性を含む7つの流体力学的テストすべてで収束性を示した。
• 新規に導入した熱伝導リミッターは、衝撃波を越える人工的エネルギー損失を効果的に防止し、物理的に不自然な冷却を回避する安定したフィードバック注入を可能にした。
• 固定パラメータでも、Cullen & Dehnen (2010) のスキームよりも高速に収束する渦の安定的進化を達成した。
• nIFTyクラスターテストでは、人工的熱伝導による効果的な相の混合のおかげで、物理的期待と一致するコア型エントロピー分布が得られた。
• Evrard崩壊テストでは、Sphenix が正確な衝撃構造とエネルギー保存を維持しており、誤った加熱や冷却が発生しなかった。
• Sphenix は現在、Swift のデフォルトのSPHスキームとなっており、オープンソースで公開されており、粒子あたりの計算コストが低く抑えられており、スケーラブルな銀河形成シミュレーションを可能にしている。