[論文レビュー] A high-order Godunov scheme for global 3D MHD accretion disks simulations. I. The linear growth regime of the magneto-rotational instability

Mario Flock, Natalia Dzyurkevich+2 more|arXiv (Cornell University)|Jun 30, 2009

Astrophysics and Star Formation Studies参考文献 56被引用数 26

ひとこと要約

本稿では、プロト惑星系円盤の熱的・化学的進化をモデル化するのに適した曲線座標系における保存的・高次精度なMHDシミュレーションフレームワークを確立する。PLUTOコードを用いた3次元グローバルMHDシミュレーションにおいて、制約伝搬（CT）法における起電力（EMF）の整合的な上流再構成が、数値不安定性を防ぐために不可欠であることを示した。研究では、HLLDおよびRoeのリーマン解法に適切な上流EMF再構成を適用することで、線形成長率が解析的予測と一致し、標準的な算術平均EMF処理を用いる手法を上回ることを確認した。

ABSTRACT

We employ the PLUTO code for computational astrophysics to assess and compare the validity of different numerical algorithms on simulations of the magneto-rotational instability in 3D accretion disks. In particular we stress on the importance of using a consistent upwind reconstruction of the electro-motive force (EMF) when using the constrained transport (CT) method to avoid the onset of numerical instabilities. We show that the electro-motive force (EMF) reconstruction in the classical constrained transport (CT) method for Godunov schemes drives a numerical instability. The well-studied linear growth of magneto-rotational instability (MRI) is used as a benchmark for an inter-code comparison of PLUTO and ZeusMP. We reproduce the analytical results for linear MRI growth in 3D global MHD simulations and present a robust and accurate Godunov code which can be used for 3D accretion disk simulations in curvilinear coordinate systems.

研究の動機と目的

• 3次元グローバル降着円盤における磁気回転不安定性（MRI）の線形成長期をシミュレートする高次精度Godunovスキームの性能を評価すること。
• Godunovスキームにおける古典的制約伝搬（CT）法で生じるEMF再構成の不一致に起因する数値不安定性を特定・解消すること。
• MRI線形成長を基準テストとして、PLUTOコードと有限差分法によるZeusMPコードの数値精度および安定性をベンチマークすること。
• プロト惑星系円盤の熱的・化学的進化をモデル化するのに適した、曲線座標系における強力で保存的かつ高精度なMHDシミュレーションフレームワークを確立すること。
• 高次リーマン解法（HLLD、Roe）に整合的な上流EMF再構成を適用することで、ZeusMPのような有限差分コードと同等の数値拡산レベルを達成できることを示すこと。

提案手法

• 曲線座標系における3次元グローバルシミュレーションに、保存的・多次元・多幾何学的Godunov型MHDコードであるPLUTOコードを用いた。
• 上流制約伝搬（UCT）法を採用し、リーマン解法に整合した上流再構成により、EMFを介して磁場の発散なしを保証した。
• 精度向上と数値拡散の低減を図るため、高次リーマン解法（HLLD、Roe、HLL、Lax-Friedrich）と区分的放物線再構成（PPM）を組み合わせた。
• 結果の妥当性を検証するため、古典的ステガラードCT法と算術平均EMF処理を用いる有限差分コードZeusMPと相互比較を行った。
• グリッドサイズを32×16×16から512×256×256まで変化させ、収束性および数値拡散効果がMRI成長率に与える影響を評価する解像度テストを実施した。
• MRI成長率を測定し、解析的予測と比較することで、数値的精度および安定性を評価した。

実験結果

リサーチクエスチョン

• RQ1 制約伝搬法におけるEMFの古典的算術平均処理は、Godunovスキームを用いたMRIシミュレーションにおいて数値不安定性を引き起こすか？
• RQ2 整合的な上流EMF再構成を施した高次Godunovスキームは、3次元グローバル円盤シミュレーションにおいてMRIの線形成長率を正確に再現できるか？
• RQ3 異なるリーマン解法（HLLD、Roe、HLL、Lax-Friedrich）の数値拡散は、有限差分コードZeusMPと比較してMRI成長シミュレーションにどのように影響するか？
• RQ4 PLUTOコードにおいて上流EMF再構成を用いた場合、MRI成長率の収束性は解像度の向上に伴いどのように変化するか？
• RQ5 リーマン解法におけるアルヴェン特性が、グローバル3次元MHD設定におけるMRIシミュレーションの精度にどの程度影響を与えるか？

主な発見

• 古典的制約伝搬法におけるEMFの算術平均処理は、特にHLLDおよびRoe解法で顕著な数値不安定性を引き起こし、チェスボードパターンとして現れる。
• UCT法による整合的上流EMF再構成により、この不安定性が効果的に抑制され、安定かつ正確なMRIシミュレーションが可能となった。
• 上流EMF再構成を施したHLLDおよびRoeリーマン解法は、解析的MRI成長率γ ≈ 0.7を測定誤差の範囲で再現した。
• 最高解像度（512×256×256）において、すべてのGodunov解法で成長率がわずかに低下したが、これは解像度依存の数値拡散または高波数効果の可能性を示唆した。
• Lax-Friedrich解法は顕著に高い数値拡散を示し、最も成長するMRIモードを解像するのに8セル以上必要であり、低解像度ではZeusMPを下回る性能を示した。
• PPM再構成を用いたHLLD解法は、遅い磁気音速特性を解像しないままでも、Roe解法と同等の性能を示し、収束が最も早く達成された。これは、強健性と効率性を示している。