[論文レビュー] Fast estimation of the maximum wave speed in the Riemann problem for the Euler equations

Jean‐Luc Guermond, Bojan Popov|arXiv (Cornell University)|Nov 9, 2015

Geophysics and Gravity Measurements被引用数 3

ひとこと要約

本稿では、体積係数を有する状態方程式を用いたEuler方程式のリーマン問題における最大波速の推定に向け、高速かつ3次収束性を示すアルゴリズムを提示する。このアルゴリズムは、所定の精度内での上界への収束を保証し、γ ∈ (1, 5/3] に属する気体を用いたシミュレーションにおいて顕著な高速化を実現する。

ABSTRACT

This paper is concerned with the construction of a fast algorithm for computing the maximum speed of propagation in the Riemann solution for the Euler system of gas dynamics with the co-volume equation of state. The novelty in the algorithm is that it stops when a guaranteed upper bound for the maximum speed is reached with a prescribed accuracy. The convergence rate of the algorithm is cubic and the bound is guaranteed for gasses with the co-volume equation of state and the heat capacity ratio $\gamma$ in the range $(1,5/3]$

研究の動機と目的

• 体積係数を有する状態方程式を用いたEuler系のリーマン解における最大波速を推定する高速かつ信頼性の高い手法の開発。
• 計算された上界が、ユーザーが指定した許容誤差内に真の最大波速を含むように保証すること。
• 波速推定の高速化を通じて、双曲型保存則シミュレーションの計算効率を向上させること。
• 実ガスの非理想状態を示すγ ∈ (1, 5/3] の範囲にまで、高速波速推定の適用範囲を拡張すること。
• 反復計算中に厳密な誤差境界を維持しつつ、3次収束率を達成すること。

提案手法

• 本アルゴリズムは、体積係数を有する状態方程式を用いたEuler方程式のリーマン解の構造に基づく反復スキームを採用する。
• 最大波速を支配する特徴方程式に対する根の探索戦略を適用し、上界への収束を保証する。
• 上界が所定の精度内に達した段階で反復を停止する収束基準を組み込む。
• 状態方程式の凸性および単調性の性質を活用することで、3次収束率を保証する。
• 熱容量比γをパラメータとして用い、数学的境界を保証するための必要な数学的構造を維持するため、有効範囲をγ ∈ (1, 5/3] に制限する。
• 数値的安定性と効率性に優れ、有限体積法または有限差分法のソルバーへの統合に適している。

実験結果

リサーチクエスチョン

• RQ1 リーマン問題における最大波速を、上界が保証される範囲で効率的に推定する方法は何か？
• RQ2 体積係数を有する状態方程式を用いたEuler方程式における波速推定で達成可能な収束率は何か？
• RQ3 ユーザーが定めた許容誤差内に上界が達成されるように保証できる高速アルゴリズムを構築可能か？
• RQ4 アルゴリズムの数学的保証を維持するためには、γにどのような制限が必要か？
• RQ5 実ガスシミュレーションにおいて、標準的手法と比較して、本アルゴリズムの速度と信頼性はどのように異なるか？

主な発見

• 本アルゴリズムは、最大波速推定において3次収束率を達成し、線形または2次収束法と比較して顕著に優れた性能を示す。
• 計算された上界は、体積係数を有する状態方程式を満たし、γ ∈ (1, 5/3] に属するすべてのガスについて、真の最大波速を数学的に含むことが保証される。
• 上界が所定の精度内に達した段階ですぐに反復を停止するため、計算効率が確保される。
• 非理想状態を示す実ガス（体積係数を有する状態方程式でモデル化）に対しても本アルゴリズムは適用可能である。
• 波速推定の高速化により、双曲型PDEソルバーにおける時間積分が高速化され、精度を損なわずに行える。
• 本手法は、γ ∈ (1, 5/3] のすべての値に対して安定かつ信頼性が保証され、二原子および多原子ガスを含む広範な物理的ガスに適用可能である。