[論文レビュー] Two-stage fourth-order accurate time discretizations for 1D and 2D special relativistic hydrodynamics

Yuhuan Yuan, Huazhong Tang|arXiv (Cornell University)|Dec 15, 2017

Advanced Numerical Methods in Computational Mathematics参考文献 33被引用数 7

ひとこと要約

本稿では、直接Euler座標系における一般化リーマン問題（GRP）法と局所的準1次元GRP状態の解析的解法を活用して、1次元および2次元の特殊相対論的流体力学に対する2段階4次精度時間離散化を提案する。このスキームは、衝撃波や複雑な波構造を高精度で捉えることに成功し、標準的なリーマン問題および滑らかな流れの数値実験によって検証されている。

ABSTRACT

This paper studies the two-stage fourth-order accurate time discretization \cite{LI-DU:2016} and applies it to special relativistic hydrodynamical equations. It is shown that new two-stage fourth-order accurate time discretizations can be proposed. With the aid of the direct Eulerian GRP (generalized Riemann problem) methods \cite{Yang-He-Tang:2011,Yang-Tang:2012} and the analytical resolution of the local "quasi 1D" GRP, the two-stage fourth-order accurate time discretizations are successfully implemented for the 1D and 2D special relativistic hydrodynamical equations. Several numerical experiments demonstrate the performance and accuracy as well as robustness of our schemes.

研究の動機と目的

• 衝撃波捕捉においてロバスト性と精度を維持する高次精度数値スキームを、特殊相対論的流体力学（RHD）に対して開発すること。
• 従来のRunge-Kutta法に比べ計算コストが高くなる問題を克服しつつ、2段階4次精度時間離散化フレームワークを1次元および2次元RHD方程式に拡張すること。
• 直接Euler座標系GRP法を用い、局所的準1次元GRP状態の解析的解法を実装することで、計算効率と精度を向上させること。
• 標準的なリーマン問題および滑らかな流れの数値実験を通じて、スキームの性能を検証すること。

提案手法

• GRPフラックス関数から導出された2段階4次精度時間離散化を採用し、全段階におけるリーマンソルバーの呼び出しを回避する。
• 時間微分係数を2次まで求めるために、直接Euler座標系GRP法を用い、局所的リーマン問題の解析的解法によりフラックスの時間微分を計算する。
• 局所的準1次元GRPの解析的解法を用いて、保存変数の高次時間微分を計算し、時間方向に4次精度を達成する。
• 保存形式の1次元および2次元特殊相対論的Euler方程式に、$̳$-法状態方程式を適用してスキームを適用する。
• 特性解析とリーマン不変量を用いて、原始変数の時間微分の発展方程式を導出し、特に音響波および衝撃波のケースに注目する。
• ローレンツ因子と相対論的フラックスを一貫して取り扱うことで、物理的一致性と安定性を確保する。

実験結果

リサーチクエスチョン

• RQ1 2段階4次精度時間離散化は、1次元および2次元特殊相対論的流体力学方程式に効果的に適応可能か？
• RQ2 本スキームは、RHDに対する従来のRunge-Kutta法と比較して、精度と効率の面で優れているか？
• RQ3 直接Euler座標系GRP法に加え、局所的GRPの解析的解法を用いることで、過大な計算コストを伴わずに高次時間精度を達成できるか？
• RQ4 本スキームは、強力な衝撃波や膨張波を含む相対論的流れにおいて、どれほどロバストに動作するか？
• RQ5 本スキームは、標準的な1次元および2次元リーマン問題および滑らかな流れのベンチマークで、どのように性能を発揮するか？

主な発見

• 収束解析により、提案された2段階4次精度時間離散化は、1次元および2次元特殊相対論的流体力学において、時間方向に4次精度を達成していることが確認された。
• 本スキームは、衝撃波や膨張波を含む複雑な波構造を、高分解能かつ最小限の振動で的確に解明できた。
• 数値実験の結果、高ローレンツ因子や強力な衝撃波を含む極端な相対論的状態でも、本手法はロバスト性を維持していることが示された。
• 各段階で全リーマンソルバーの呼び出しを回避するため、標準的な4段階Runge-Kutta法と比較して計算コストが顕著に削減された。
• 滑らかな解において、L1およびL∞ノルムで最適な収束率を示し、理論的な精度順序が裏付けられた。
• 局所的準1次元GRPの解析的解法により、2次までの時間微分を正確かつ効率的に計算でき、高次精度の根幹をなしている。