[論文レビュー] Reliability of astrophysical jet simulations in 2D: On inter-code reliability and numerical convergence

Martin Krause, M. Camenzind|ArXiv.org|Oct 22, 2001

Computational Fluid Dynamics and Aerodynamics参考文献 12被引用数 19

ひとこと要約

本研究は、NIRVANA MHDコードを用いた2次元銀河的ジェットシミュレーションの数値的収束性とコード間信頼性を評価する。Lind et al. 1989およびKössl & Müller 1988の2つの既存コードと比較した結果、流体的状況では約100 ppbでグローバルジェット特性が収束するが、磁化ジェットシミュレーションでは移動するマッハディスクのため収束せず、高解像度の計算で小スケールのケルビン＝ヘルムホルツ不安定性およびレイノルズ＝ターレイユ不安定性といった新たな不安定性が顕在化し、質量の取り込みが強化された。

ABSTRACT

In the present paper, we examine the convergence behavior and inter-code reliability of astrophysical jet simulations in axial symmetry. We consider both, pure hydrodynamic jets and jets with a dynamically significant magnetic field. The setups were chosen to match the setups of two other publications, and recomputed with the MHD code NIRVANA. We show that NIRVANA and the two other codes give comparable, but not identical results. We find that some global properties of a hydrodynamical jet simulation, like e.g. the bow shock velocity, converge at 100 points per beam radius (ppb) with NIRVANA. The situation is quite different after switching on the toroidal magnetic field: In this case, global properties converge even at 10 ppb. In both cases, details of the inner jet structure and especially the terminal shock region are still insufficiently resolved, even at our highest resolution of 70 ppb in the magnetized case and 400 ppb for the pure hydrodynamic jet. In the case of our highest resolution simulation, we can report two new features: First, small scale Kelvin-Helmholtz instabilities are excited at the contact discontinuity next to the jet head. This slows down the development of the long wavelength Kelvin-Helmholtz instability and its turbulent cascade to smaller wavelengths. Second, the jet head develops Rayleigh-Taylor instabilities which manage to entrain an increasing amount of mass from the ambient medium with resolution. This region extends in our highest resolution simulation over 2 jet radii in the axial direction.

研究の動機と目的

• NIRVANA MHDコードを用いた2次元銀河的ジェットシミュレーションの数値的収束性を評価すること。
• NIRVANA、Lind et al. 1989、Kössl & Müller 1988の複数コード間での結果比較を通じて、コード間信頼性を評価すること。
• 人工粘性および数値的解像度がジェット構造および衝撃波ダイナミクスのシミュレーションに与える影響を調査すること。
• グローバルジェットパラメータおよび内部構造（終端衝撃波および接触不連続面など）の収束に必要な解像度要件を特定すること。
• 流体的および磁化ジェット状況下で、解像度を変化させた際の不安定性（ケルビン＝ヘルムホルツおよびレイノルズ＝ターレイユ）の発生を検討すること。

提案手法

• NIRVANA MHDコードを用いて、空間解像度を10～400ポイント/ビーム半径の範囲で変化させた、ベンチマーク用ジェットシミュレーション（流体的および磁化）を再計算した。
• Lind et al. (1989)およびKössl & Müller (1988)の先行研究と同一の条件を再現するため、軸対称性（2次元）を適用した。
• 有限体積法を用いたMHDソルバーに、Rusanovフラックスと人工粘性を適用し、衝撃波および不連続面を扱った。
• 解像度レベルごとに、バウシュック速度やマッハディスク位置といったグローバル特性の収束を追跡した。
• 時系列データおよび等高線プロットの分析を通じて、ジェット先端およびコッコウ内の不安定性の発展および構造的進化を特定した。
• 高解像度シミュレーション（流体的では最大400 ppb、磁化では70 ppb）を実施し、小スケール乱流や質量の取り込みといった解像度依存の特徴を検出した。

実験結果

リサーチクエスチョン

• RQ1 異なるMHDコードが同一の銀河的ジェット設定に対して収束する結果をどれほど得られるか？
• RQ2 数値的解像度が、バウシュック速度やマッハディスク位置といったグローバルジェットパラメータの収束にどのように影響するか？
• RQ3 人工粘性が、特に衝撃波および接触不連続面の解像度において、コード出力の差異を生じる役割を果たすか？
• RQ4 高解像度シミュレーションは、低解像度では見逃される小スケールの不安定性や強化された質量の取り込みといった新たな物理的特徴を明らかにするか？
• RQ5 磁化ジェットシミュレーションにおいて、NIRVANAではマッハディスクが移動しており収束しないにもかかわらず、グローバルパラメータには表面的な収束が見られるのはなぜか？

主な発見

• NIRVANAコードにおける流体的ジェットシミュレーションのグローバル特性は、約100ポイント/ビーム半径（ppb）で収束する。
• 磁化ジェットシミュレーションでは、グローバル特徴（例：バウシュックの形状や速度）が10 ppb程度でも収束するが、マッハディスクが流入境界に向かって移動しており、収束していないことが示唆される。
• 最高解像度（流体的では400 ppb、磁化では70 ppb）に達した場合、長波長モードの成長を抑制し、乱流スケーリングの遅延を引き起こす小スケールケルビン＝ヘルムホルツ不安定性がジェット先端に顕在化する。
• 最高解像度の磁化ジェットシミュレーションでは、レイノルズ＝ターレイユ不安定性が発達し、周囲媒体からの質量が増加して取り込まれ、軸方向に2ジェット半径にわたって広がる。
• 高解像度であっても、内部ジェット構造および終端衝撃波領域は依然として十分に解像されておらず、流体的ケースで400 ppbでも乱流的進化が継続している。
• 磁化ケースにおけるFLOWとNIRVANAの差異（特にNIRVANAにおける収束しないマッハディスク）は、従来のシミュレーション結果に疑問を呈し、磁場の影響下ではコード依存の収束行動が顕著である可能性を示唆している。