[論文レビュー] Exact solution of the 1D Riemann problem in Newtonian and relativistic hydrodynamics

F. D. Lora-Clavijo, J. P. Cruz-Pérez+2 more|arXiv (Cornell University)|Mar 16, 2013

Cosmology and Gravitation Theories参考文献 2被引用数 11

ひとこと要約

本稿は、ニュートン的および相対論的流体力学における1次元リーマン問題の正確な解の詳細な段階的導出および実装を提示する。リーマン不変量、衝撃波および接触不連続性の条件、および根の探索技術を用いた包括的な数値的レシピを提供し、天体物理学的シミュレーションにおける数値的流体力学コードの正確な検証を可能にする。

ABSTRACT

Some of the most interesting scenarios that can be studied in astrophysics, contain fluids and plasma moving under the influence of strong gravitational fields. To study these problems it is required to implement numerical algorithms robust enough to deal with the equations describing such scenarios, which usually involve hydrodynamical shocks. It is traditional that the first problem a student willing to develop research in this area is to numerically solve the one dimensional Riemann problem, both Newtonian and relativistic. Even a more basic requirement is the construction of the exact solution to this problem in order to verify that the numerical implementations are correct. We describe in this paper the construction of the exact solution and a detailed procedure of its implementation.

研究の動機と目的

• ニュートン流体力学における1次元リーマン問題の正確解を、明確で明示的かつ実装可能な手順として提供すること。
• 既存の文献でより体系的に提示されていない相対論的リーマン問題に対しても、同じ詳細なアプローチを拡張すること。
• 数値的流体力学コードの実装を検証するための実用的参考資料として提供すること。
• リーマン不変量および特徴的速度を用いて、両状態における衝撃波および接触不連続性の条件の構築を明確にすること。
• 高分解能衝撃捕捉法およびスムーズ・パーティクル流体力学法の強固なテストを、正確なベンチマーク解を用いて可能にすること。

提案手法

• 密度、速度、圧力の項で表されるフラックスを用いて、原始変数および保存変数を用いた保存形でのオイラー方程式を定式化する。
• フラックスヤコビアンの固有構造を導出し、特徴的速度（v−a、v、v+a）および対応する固有ベクトルを特定する。
• 衝撃波を貫いてランキン＝ユゴニオ条件を適用し、不連続性をはさんで状態を関連付ける。衝撃速度およびジャンプ条件を用いる。
• 接触不連続性の条件（接触面をはさんで速度が等しい）を用いて、スター状態の圧力に関する超越方程式を導出する。
• 根の探索アルゴリズムを用いて、スター状態の圧力に関する超越方程式を解き、解構造の完全な再構成を可能にする。
• 左状態、衝撃波、接触不連続性、第二の衝撃波、右状態の6つの領域にわたる完全な解を構築し、各領域における圧力、速度、密度の明示的表現を与える。

実験結果

リサーチクエスチョン

• RQ1 ニュートン流体力学における1次元リーマン問題の正確解を、コード検証の目的で体系的かつ実装可能な形で導出・実装するにはどうすればよいか？
• RQ2 リーマン問題解における相対論的補正の主要な点は何か。また、ニュートン的場合とどのように異なるか？
• RQ3 複数の波を伴う状況において、リーマン不変量および衝撃ジャンプ条件を用いてスター状態の変数をどのように決定できるか？
• RQ4 相対論的状況において、接触不連続性の条件から生じる超越方程式を解くために必要な数値的手法は何か？
• RQ5 ニュートン的および相対論的両状態において、衝撃波および接触不連続性を含むすべての領域にわたる解構造をどのように再構成できるか？

主な発見

• ニュートン的1次元リーマン問題の正確解は、特徴的速度およびリーマン不変量を用い、衝撃波および接触不連続性の条件を適用してスター状態を決定することで構築される。
• 相対論的状況では、接触不連続性をはさんで速度が一致する条件から導かれる超越方程式を根探索法で解き、スター状態の圧力を求める。
• 衝撃波の速度は相対論的ジャンプ条件を用いて計算され、ローレンツ因子および圧力・速度のジャンプに依存して表される。
• 解構造は6つの明確な領域から成る：左状態、衝撃波、接触不連続性、第二の衝撃波、右状態。各領域では原始変数が定数となる。
• 任意の時刻における解の完全な再構成が可能であり、各領域における圧力、速度、密度の明示的表現が与えられる。
• このアプローチは直接実装可能であり、ニュートン的および相対論的天体物理学における数値的流体力学コードの強固なベンチマークとして機能する。