[論文レビュー] The relativistic Burgers equation on a FLRW background and its finite volume approximation

Tuba Ceylan, Philippe G. LeFloch+1 more|arXiv (Cornell University)|Dec 26, 2015

Computational Fluid Dynamics and Aerodynamics参考文献 6被引用数 3

ひとこと要約

本稿では、圧力のない流体仮定の下で相対論的Euler方程式から導かれる、Friedmann–Lemaître–Robertson–Walker (FLRW)時空背景上での相対論的Burgers方程式を定式化する。空間一様解を保存し、さまざまな曲率パラメータ（k = -1, 0, 1）における衝撃波と希薄化波を正確に捉えるwell-balancedな有限体積スキームを導入し、数値結果からk = -1では特徴速度が増加するため、より速い波動伝播が観測された。

ABSTRACT

A relativistic generalization of the inviscid Burgers equation was proposed by LeFloch, Makhlof, and Okutmustur and then investigated on a Schwarzschild background. Here, we extend their analysis to a Friedmann-Lemaitre-Robertson-Walker (FLRW) background. This problem is more challenging due to the existence of non-trivial spatially homogeneous solutions. First, we derive the relativistic Burgers model of interest and determine its spatially homogeneous solutions. Second, we design a numerical scheme based on the finite volume methodology, which is well-balanced in the sense that spatially homogeneous solutions are preserved at the discrete level of approximation. Numerical experiments demonstrate the efficiency of the proposed method for weak solutions containing shock waves.

研究の動機と目的

• 時空の膨張と空間曲率を考慮したFLRW背景時空上でのBurgers方程式の相対論的一般化を導出すること。
• 曲がった背景の下で非自明となる空間一様解を同定・解析すること。
• 離散的に空間一様解を保存するwell-balancedな有限体積数値スキームを設計すること。
• 弱解（衝撃波と希薄化波を含む）に対して、曲率パラメータ（k = -1, 0, 1）を変化させた数値的妥当性を検証すること。

提案手法

• 圧力がゼロ（p ≡ 0）である圧力のない流体仮定を用いて、曲がったFLRW時空上での相対論的Euler方程式から相対論的Burgers方程式を導出する。これにより幾何学的保存則が得られる。
• スケール因子a(t)と曲率パラメータkを用いてモデルを保存形式で表現し、源項を含む双曲型バランス法則を得る。
• 数値フラックスと再構成技術に基づき、1次および2次精度のwell-balancedな有限体積Godunov型スキームを構築する。
• 数値フラックスと源項が空間一様解と整合することを保証することで、well-balanced性を実装する。
• セル界面における数値フラックスを計算するために、Rusanov型リーマンソルバーを用い、フラックスと源項はセル平均で評価する。
• k = -1, 0, 1の各ケースに対して、衝撃波と希薄化波を含むテストケースにスキームを適用し、収束性とロバスト性を評価する。

実験結果

リサーチクエスチョン

• RQ1 非自明な空間曲率（k = -1, 0, 1）を持つFLRW時空上での相対論的Burgers方程式の挙動はいかなるものか？
• RQ2 FLRW背景上での相対論的Burgersモデルにおける空間一様解の性質は何か？
• RQ3 空間一様解を保存しつつ、衝撃波と希薄化波を正確に捉える有限体積スキームを構築可能か？
• RQ4 曲率パラメータkが、数値解における衝撃波および希薄化波の伝播速度と構造に与える影響は何か？
• RQ5 時間依存かつ曲がった背景上での不連続解に対して、提案されたwell-balancedスキームの収束性とロバスト性はいかがなものか？

主な発見

• 時間依存スケール因子a(t)と曲率kのおかげで、FLRW時空上での相対論的Burgers方程式は非自明な空間一様解を有する。
• 本モデルの特徴速度は(1 - kr²)^{1/2}で与えられ、kが小さくなるほど増加するため、k = -1ではk = 0およびk = 1よりも波動伝播が速くなる。
• 提案された有限体積スキームはwell-balancedであり、離散的に空間一様解を保存する。これは、宇宙論的背景上での長時間シミュレーションにおいて極めて重要である。
• 数値実験により、すべての3つの曲率ケースにおいて衝撃波と希薄化波を正確に捉えるスキームのロバスト性と効率性が確認された。
• Godunovスキームにより2次精度が達成され、希薄化波および衝撃波両方の解に対して効果的に収束する。
• 図3および図4の赤線（k = -1）は、緑線（k = 0）および青線（k = 1）よりも速い収束および伝播を示しており、特徴速度からの解析的予測を裏付けている。