[論文レビュー] A Roe-type Riemann solver based on the spectral decomposition of the equations of Relativistic Magnetohydrodynamics

J. Ibáñez, M. Á. Aloy+4 more|arXiv (Cornell University)|Dec 11, 2010

Cosmology and Gravitation Theories被引用数 3

ひとこと要約

本稿では、フラックスヤコビアンの完全なスペクトル分解に基づき、特異状態でさえも明確に定義された右および左固有ベクトルを保証する、相対論的磁気流体力学（RMHD）の新規ロウ型リーマンスolverを提示する。このスolverは、分散性においてHLLおよびHLLCを上回り、HLLDと同等の精度を達成しながら、多次元シミュレーションにおいて計算効率が向上する。

ABSTRACT

In a recent paper (Antón et al. 2010) we have derived sets of right and left eigenvectors of the Jacobians of the relativistic MHD equations, which are regular and span a complete basis in any physical state including degenerate ones. We present a summary of the main steps followed in the above derivation and the numerical experiments carried out with the linearized (Roe-type) Riemann solver we have developed, and some note on the (non-)convex character of the relativistic MHD equations.

研究の動機と目的

• 物理的状態、特に特異状態を含むすべての状態で、RMHDフラックスヤコビアンの正則かつ完全な右および左固有ベクトルのセットを導出すること。
• 保存変数におけるRMHD方程式のスペクトル分解に基づく線形化（ロウ型）リーマンスローバーを開発すること。
• 標準的な1次元および2次元テスト問題を用いて、HLL、HLLC、HLLDなどの既存スローバーと比較して、スローバーの数値的性能を評価すること。
• 特に特異状態における磁気音速波が真に非線形から線形的特異性へ移行する際の、RMHD方程式の非凸性を調査すること。
• 古典的MHDの戦略（例：Brio & Wu, 1988）を相対論的領域に拡張し、極端な天体物理学的条件下での数値的安定性と精度を確保すること。

提案手法

• 保存変数におけるスペクトル分解を用いて、特異状態を含むすべての物理的状態で正則かつ完全な右および左固有ベクトルのセットを導出。
• 再正規化された固有ベクトルセットに基づく完全波分解（FWD）技術を用い、線形化リーマンスローバー枠組みで数値フラックスを計算。
• 標準的なロウ型フラックス式を採用：$\widehat{{\bf f}}_{j\pm{1\over 2}}=\frac{1}{2}\left({\bf f}({\bf u}_{j\pm{1\over 2}}^{L})+{\bf f}({\bf u}_{j\pm{1\over 2}}^{R})-\sum_{\alpha=1}^{p}|\widetilde{\lambda}_{\alpha}|\Delta\widetilde{\omega}_{\alpha}\widetilde{r}_{\alpha}\right)$、ここで$\widetilde{\lambda}_{\alpha}$は固有値、$\widetilde{r}_{\alpha}$は右固有ベクトル。
• スペクトル分解の基盤として保存変数を用い、中間変換を経由せずにフラックスを直接計算可能にした。
• 標準的な1次元および2次元RMHDテスト問題（ショックチューブ、ブラスト波など）を用いてスローバーの妥当性を検証し、精度、安定性、分散性を評価。
• 固有値の勾配$\vec{\nabla}_{\bf u}\lambda_{\alpha}\cdot{\bf r}_{\alpha}$の計算を通じて、RMHD系の凸性を分析し、磁気音速場が真に非線形性を失う条件を同定。

実験結果

リサーチクエスチョン

• RQ1 物理的状態、特に特異状態を含むすべての状態で、RMHDフラックスヤコビアンの完全かつ正則な右および左固有ベクトルセットを導出可能か？
• RQ2 スペクトル分解に基づくロウ型リーマンスローバーの性能は、HLL、HLLC、HLLDスローバーと比較して、精度および分散性の観点でどのように異なるか？
• RQ3 非凸性、特に特異状態における非凸性の影響は何か。また、波の挙動にどのように影響を与えるか？
• RQ4 FWDベースのスローバーの計算コストは1次元では依然として高価であるが、多次元シミュレーションではその効率が向上するか？
• RQ5 状態方程式の熱力学的性質は、RMHD系の凸性および安定性にどの程度影響を与えるか？

主な発見

• 導出された右および左固有ベクトルのセットは、すべての物理的状態、特に特異状態を含むすべての状態で正則であり、完全な基底を形成する。これにより、特徴的構造における従来の病理的問題が解消された。
• FWDベースのロウ型リーマンスローバーは、標準的なテスト問題においてHLLおよびHLLCスローバーと比較して、より優れた安定性と低い数値的分散性を示した。
• HLLDスローバーと同等の精度を達成しており、複雑な波構造を高精度で捉えることが知られている。
• 1次元ではHLL系スローバーより計算コストが高かったが、多次元シミュレーションでは、より優れた波分解能のおかげで相対的な効率が著しく向上した。
• 分析により、RMHDにおける磁気音速波は特異状態で線形的特異性を示すことが確認され、古典的MHDと同様の非凸性挙動を示す。この移行は、状態方程式およびローレンツ因子に依存する。
• 非凸性は超相対論的領域（$a \to 1$ または $b \to 1$）で発現し、磁気音速場の特異性が凸性破壊の主な要因であることが、$\vec{\nabla}_{\bf u}\lambda_{\alpha}\cdot{\bf r}_{\alpha}$ の理論的解析と一致した。