[論文レビュー] An improved SPH scheme for cosmological simulations

Alexander M. Beck, Giuseppe Murante+10 more|arXiv (Cornell University)|Feb 25, 2015

Fluid Dynamics Simulations and Interactions参考文献 12被引用数 6

ひとこと要約

本稿では、Wendlandカーネル関数の統合、高次勾配とせん断制限を備えた時間依存人工粘性、粒子目覚まし時間刻み制限、および新規の時間依存人工拡散モデルを組み込んだ、宇宙論的シミュレーションに向けた高度な滑らかさ粒子法（SPH）スキームを提示する。この実装により、数値的ノイズが著しく低減され、流体相間のエントロピー混合が向上し、より現実的な銀河およびクラスタ構造（例えば、冷たく、より広がった銀河ディスクやクラスタ内のエントロピー核）が得られる。これは、すべての流体力学的テストにおいて標準的な GADGET-SPH を上回る性能を示す。

ABSTRACT

We present an implementation of smoothed particle hydrodynamics (SPH) with improved accuracy for simulations of galaxies and the large-scale structure. In particular, we combine, implement, modify and test a vast majority of SPH improvement techniques in the latest instalment of the GADGET code. We use the Wendland kernel functions, a particle wake-up time-step limiting mechanism and a time-dependent scheme for artificial viscosity, which includes a high-order gradient computation and shear flow limiter. Additionally, we include a novel prescription for time-dependent artificial conduction, which corrects for gravitationally induced pressure gradients and largely improves the SPH performance in capturing the development of gas-dynamical instabilities. We extensively test our new implementation in a wide range of hydrodynamical standard tests including weak and strong shocks as well as shear flows, turbulent spectra, gas mixing, hydrostatic equilibria and self-gravitating gas clouds. We jointly employ all modifications; however, when necessary we study the performance of individual code modules. We approximate hydrodynamical states more accurately and with significantly less noise than standard SPH. Furthermore, the new implementation promotes the mixing of entropy between different fluid phases, also within cosmological simulations. Finally, we study the performance of the hydrodynamical solver in the context of radiative galaxy formation and non-radiative galaxy cluster formation. We find galactic disks to be colder, thinner and more extended and our results on galaxy clusters show entropy cores instead of steadily declining entropy profiles. In summary, we demonstrate that our improved SPH implementation overcomes most of the undesirable limitations of standard SPH, thus becoming the core of an efficient code for large cosmological simulations.

研究の動機と目的

• 宇宙論的シミュレーションにおける物理的でないエントロピー分布を引き起こす、接触不連続面における混合の悪さや不適切な表面張力といった、長年のSPHの数値的問題に対処すること。
• 数値的ノイズの低減と流体混合の向上により、銀河および大規模構造形成における流体力学的シミュレーションの解像度と精度を向上させること。
• 保存則を維持しながら衝撃波、せん断流れ、乱流不安定性を正確に捉えることが可能な、強固で安定的かつ効率的なSPH実装を開発すること。
• 包括的な流体力学的および宇宙論的テスト問題を用いて、新規スキームを標準SPHおよびEuler形式コードと比較して検証すること。
• 放射冷却、星形成、超新星降着反動を含む物理的整合性のある流体力学ソルバーを提供することで、より現実的な銀河形成シミュレーションを可能にすること。

提案手法

• 粒子サンプリングにおけるノイズ低減とカーネル補間精度の向上のため、Wendlandカーネル関数の実装。
• 滑らかな流れにおける物理的でない散逸を抑える一方で衝撃波を捉えるために、高次勾配再構成とせん断流れ制限を備えた時間依存人工粘性スキームの採用。
• 不活性な粒子が流体サンプリングを歪めることを防ぎ、動的領域における時間分解能を向上させるために、粒子目覚まし時間刻み制限メカニズムの導入。
• 重力に起因する圧力勾配を補正し、流体界面を越えたエントロピー混合を強化するため、時間依存人工拡散モデルの統合。
• 接触不連続面を適切に扱い、多相流における安定性を向上させるために、圧力に基づく運動方程式の定式化の採用。
• 開発者版を用いて完全な互換性と再現性を確保するため、GADGET-3コードベースへの改良SPHスキームの統合。

実験結果

リサーチクエスチョン

• RQ1 新規SPHスキームは、ケルビン＝ヘルムホルツ不安定性やレイノルズ＝ターレイユ不安定性といった流体力学的不安定性の解像度と精度をどのように向上させるか？
• RQ2 時間依存人工拡散モデルは、宇宙論的シミュレーションにおける異なる流体相間のエントロピー混合をどの程度向上させるか？
• RQ3 粒子目覚まし時間刻み制限メカニズムは、強い衝撃波やせん断流れの際の流体サンプリングにおける数値的歪みをどのように低減するか？
• RQ4 新規SPHスキームは、銀河クラスタにおけるガス密度、温度、エントロピーの半径方向プロファイルを、標準GADGET-SPHおよびEulerコードと比較してどの程度正確に再現するか？
• RQ5 新規SPH定式化は、銀河の構造的特性（例えば、ディスク温度やスケール長）にどのような影響を及えるか？

主な発見

• 標準GADGET-SPHと比較して、数値的ノイズが低減され、特に弱いおよび強い衝撃波テストにおいて、より正確な流体力学的状態が得られる。
• このスキームは、標準SPHでしばしば抑制されるケルビン＝ヘルムホルツ不安定性やレイノルズ＝ターレイユ不安定性の発達を的確に捉えている。
• 異なる流体相間のエントロピー混合が顕著に向上し、シミュレートされた銀河クラスタにエントロピー核が形成される一方で、単調に減少するプロファイルが生じない。
• 新規スキームで形成された銀河ディスクは、冷たく、より広がっており、熱的および力学的安定性の向上を示している。
• nIFTy宇宙論比較プロジェクトにおいて、新規スキームはEulerコードおよび最新のSPH実装と同等の性能を示し、密度、温度、エントロピーの半径方向プロファイルにおいて優れた一致を示している。
• カーネル、粘性、時間刻み制御、人工拡散のすべての改良を統合した結果、大規模宇宙論的シミュレーションに適した強固で安定的かつ物理的に整合性のあるソルバーが得られた。