[論文レビュー] A Simulation Study of Ultra-Relativistic Jets -- I. A New Code for Relativistic Hydrodynamics

Jeongbhin Seo, Hyesung Kang+3 more|arXiv (Cornell University)|Jun 8, 2021

Computational Fluid Dynamics and Aerodynamics参考文献 59被引用数 14

ひとこと要約

本稿では、5次精度の有限差分WENOスキームに基づく新しい高次精度特殊相対論的流体力学（RHD）コードを提案する。主な特徴には、SSPRK時間積分、高次元横方向フラックス平均、およびカルブンクル不安定性を抑えるために改変された固有値が含まれる。このコードは、銀河団内媒体（ICM）における超相対論的ジェットのシミュレーションを最適化しており、WENO-Z、SSPRK、およびRC状態方程式が、衝撃波、乱流、ジェットダイナミクスのシミュレーションにおいて安定的かつ高精度であるとされるファイドシャー設定として採用されている。

ABSTRACT

In an attempt to investigate the structures of ultra-relativistic jets injected into the intracluster medium (ICM) and the associated flow dynamics, such as shocks, velocity shear, and turbulence, we have developed a new special relativistic hydrodynamic (RHD) code in the Cartesian coordinates, based on the weighted essentially non-oscillatory (WENO) scheme. It is a finite difference scheme of high spatial accuracy, which has been widely employed for solving hyperbolic systems of conservation equations. The code is equipped with different WENO versions, such as the 5th-order accurate WENO-JS (Jiang & Shu 1996), WENO-Z, and WENO-ZA, and different time integration methods, such as the 4th-order accurate Runge-Kutta (RK4) and strong stability preserving RK (SSPRK), as well as the implementation of the equations of state (EOSs) that closely approximate the EOS of the single-component perfect gas in relativistic regime. In addition, it incorporates a high-order accurate averaging of fluxes along the transverse directions to enhance the accuracy of multi-dimensional problems, and a modification of eigenvalues for the acoustic modes to effectively control the carbuncle instability. Through extensive numerical tests, we assess the accuracy and robustness of the code, and choose WENO-Z, SSPRK, and the EOS suggested in Ryu et al. (2006) as the fiducial setup for simulations of ultra-relativistic jets. The results of our study of ultra-relativistic jets using the code is reported in an accompanying paper (Seo et al. 2021, Paper II).

研究の動機と目的

• 銀河団内媒体（ICM）における超相対論的ジェットのシミュレーションを目的とした高次元で安定性の高いRHDコードの開発。
• 多次元流れにおけるカルブンクル不安定性や横方向拡散といった数値的課題の解決。
• ジェット構造およびダイナミクスの高精度なシミュレーションに最適な数値スキームと状態方程式の組み合わせの特定。

提案手法

• 5次精度の有限差分WENOスキーム（WENO-JS、WENO-Z、WENO-ZA）を用いた高解像度再構成の実装。
• 安定性と衝撃波支配の流れにおける精度を確保するため、4次精度のSSPRKおよびRK4時間積分法の採用。
• 多次元問題における精度向上のため、高次元横方向フラックス平均の統合。
• 音響モードの固有値を改変することで、解像度を損なわずにカルブンクル不安定性を抑制。
• 衝撃領域における断熱指数の変化（ジェットではγ ≈ 1.45、衝撃領域ではγ ≈ 4/3および5/3）を正確にモデル化するため、RC状態方程式（EOS）の採用。
• 1次元、2次元、3次元問題における広範なテストを通じて、精度、安定性、性能の妥当性を検証。

実験結果

リサーチクエスチョン

• RQ1 超相対論的ジェットシミュレーションにおいて、どのWENO再構成スキームが解像度と安定性の最良のバランスを提供するか？
• RQ2 異なる時間積分スキーム（SSPRK対RK4）は、大きな横方向速度を伴う衝撃波のシミュレーションにどのように影響するか？
• RQ3 高次元横方向フラックス平均は、多次元相対論的流体力学における精度を向上させるか？
• RQ4 固有値の修正は、相対論的ジェットのバッシュショック面上でのカルブンクル不安定性をどの程度抑制できるか？
• RQ5 どの状態方程式（ID、TM、RC）が、衝撃加熱および混合過程におけるジェットおよびICM流体の熱力学的進化を最もよく捉えるか？

主な発見

• WENO-Zスキームが解像度と安定性の最良のバランスを提供し、WENO-JS（過剰に拡散的）およびWENO-ZA（極端な条件下で過敏）を上回る性能を示した。
• 大規模な横方向速度を伴う衝撃波支配の流れにおいて、SSPRK時間積分はRK4に比べて優れた結果を示し、安定性を確保するとともに振動を低減した。
• 高次元横方向フラックス平均は、特にココーンやバックフラウインター面上のような複雑な流れ領域で、多次元シミュレーションの精度を向上させた。
• φ ≈ 5–10の固有値修正により、バッシュショック面上のカルブンクル不安定性が効果的に抑制された。高出力ジェット（Qjet ≳ 10⁴⁶ erg s⁻¹）ではφ ≈ 10、低出力ジェットではφ ≈ 5が適切であった。
• RC状態方程式はTMと同等の結果を生み、IDに比べて相対論的完全気体EOSに近い性能を示し、γ ≈ 1.45からγ ≈ 5/3への断熱指数の遷移を正確にモデル化できるため、正確な断熱指数遷移のモデリングのデフォルト選択肢となった。
• WENO-Z、SSPRK、RC EOS、固有値修正を組み合わせたファイドシャー設定を用いたシミュレーションは、KH不安定性、バッシュショック、乱流を示すココーン構造を含む複雑なジェット形状を成功裏に捉えた。