[論文レビュー] Cosmological AMR MHD with Enzo

David C. Collins, Hao Xu+3 more|arXiv (Cornell University)|Feb 16, 2009

Cosmology and Gravitation Theories参考文献 47被引用数 49

ひとこと要約

本稿では、宇宙論的自己適合メッシュリファインメント（AMR）コードEnzoを、宇宙論的環境下での理想磁気流体力学（MHD）をシミュレートするために拡張したEnzoMHDを提案する。流体フラックスには高次精度のGodunovリーマンソルバーを採用し、磁場の発散をゼロに保つために3次精度の発散自由再構成を用いた制約伝搬法を適用し、重力と宇宙膨張には作用素分割を適用することで、機械精度に近い発散制御と、宇宙論的および非宇宙論的テスト問題においても安定した性能を達成している。

ABSTRACT

In this work, we present MHDEnzo, the extension of the cosmological code Enzo to include the effects magnetic fields through the ideal MHD approximation. We use a higher order Godunov Riemann solver for the computation of interface fluxes. We use two constrained transport methods to compute the electric field from those interface fluxes, which simultaneously advances the induction equation and maintains the divergence of the magnetic field. A third order divergence free reconstruction technique is used to interpolate the magnetic fields in the block structured AMR framework already extant in Enzo. This reconstruction also preserves the divergence of the magnetic field to machine precision. We use operator splitting to include gravity and cosmological expansion. We then present a series of cosmological and non cosmological tests problems to demonstrate the quality of solution resulting from this combination of solvers.

研究の動機と目的

• 宇宙論的AMRコードEnzoに、自己一貫的かつ発散自由な磁場の時間発展を実装するための理想MHDを統合すること。
• ブロック構造AMR環境下で、磁場の発散を機械精度まで維持する数値フレームワークを構築すること。
• 銀河団、星形成、および星間媒体の力学を含む、宇宙構造形成における磁場の正確なシミュレーションを可能にすること。
• 宇宙論的および非宇宙論的テスト問題のスイートを用いて、実装の数値的精度と安定性を検証すること。

提案手法

• 宇宙膨張、重力、および流体力学を含む共動座標系における理想MHD方程式を採用する。
• セル界面におけるフラックス計算に高次精度のGodunovリーマンソルバーを用い、高分解能の流体解を得る。
• 磁場の発散をゼロに保つために、2種類の制約伝搬法を用いて電場を計算し、誘導方程式を進める。
• AMRレベル間での磁場の補間に3次精度の発散自由再構成技術を適用し、発散を機械精度まで維持する。
• 重力と宇宙膨張の処理を別々に行う作用素分割を適用し、時間中心化された更新により整合性を保つ。
• フラックス補正と表面平均化を施した、時間先進型でステガードグリッドのアルゴリズムを実装し、AMRレベル間での保存量の正確な移行を保証する。

実験結果

リサーチクエスチョン

• RQ1 Enzoのようなブロック構造AMRコードは、磁場の発散をゼロに保ちながら、宇宙論的環境下での理想MHDを正確にシミュレートできるか？
• RQ2 高次精度のGodunovスキームと制約伝搬の組み合わせは、適応グリッド上での複雑なMHD流れをどれほど正確に解明できるか？
• RQ3 宇宙論的および非宇宙論的テスト問題において、磁場の時間発展の精度と安定性はいかほどか？
• RQ4 3次精度の発散自由再構成は、AMRのリファイン境界で∇·B = 0を高精度に維持できるか？
• RQ5 作用素分割アプローチは、時間依存シミュレーションにおける重力、膨張、およびMHDの結合をどれほどうまく処理できるか？

主な発見

• EnzoMHDの実装は、3次精度の再構成と制約伝搬を用いることで、すべてのAMRレベルで磁場の発散を機械精度まで維持している。
• コードは標準的なMHDテスト問題（Orszag-Tang渦、ローター、Brio-Wuショック）を正確に解き、滑らかな流れでは2次精度を示し、ショック波の捕捉も安定に実行している。
• 宇宙論的シミュレーションにおいて、EnzoMHDは乱流ダイナモによる磁場強化を正しく捉え、収縮構造における磁場トポロジーの維持も行っている。
• 時間中心化された更新とフラックス補正の適用により、粗いAMRレベルと細かいAMRレベル間での保存量および電場の保存的移行が保証されている。
• 長時間にわたる統合においても、発散∇·Bの発散や物理的に不自然な増大は観測されておらず、安定で正確な性能を示している。
• 本手法により、銀河団、星形成領域、および星間媒体のような天体系における磁場の高精度なシミュレーションが可能となり、数値収束性の検証も完了している。