[論文レビュー] High order well-balanced finite volume methods for multi-dimensional systems of hyperbolic balance laws

Jonas P. Berberich, Praveen Chandrashekar+1 more|arXiv (Cornell University)|Mar 12, 2019

Computational Fluid Dynamics and Aerodynamics被引用数 4

ひとこと要約

本稿では、任意の既知の定常解または時間依存解（解析的・離散的）を正確に維持する一般化された高次 well-balanced 有限体積フレームワークを提示する。この手法は、解からの偏差を進化させることで、目標解の正確な維持を実現する。本手法は高次精度を保ち、コード変更が最小限で、高次スキームですら計算コストがわずかに5–30%増加するにとどまる。

ABSTRACT

We introduce a general framework for the construction of well-balanced finite volume methods for hyperbolic balance laws. We use the phrase well-balancing in a broader sense, since our proposed method can be applied to exactly follow any solution of any system of hyperbolic balance laws in multiple spatial dimensions and not only time independent solutions. The solution has to be known a priori, either as an analytical expression or as discrete data. The proposed framework modifies the standard finite volume approach such that the well-balancing property is obtained and in case the method is high order accurate, this is maintained under our modification. We present numerical tests for the compressible Euler equations with and without gravity source term and with different equations of state, and for the equations of compressible ideal magnetohydrodynamics.

研究の動機と目的

• 多次元空間における任意の双曲型バランス法のシステムに適用可能な、高次 well-balanced 有限体積法を構築する一般フレームワークの開発。
• 特定の状態方程式や静水圧的仮定に制限されない、任意の既知の解（静的、定常的、時間依存的）の正確な保存を可能にする。
• 既存の有限体積ソルバーと互換性を保ちつつ、ソース項離散化および再構成戦略の変更のみで、高次精度を維持する。
• 圧縮性Euler方程式に重力が作用する系や理想MHD系など、複雑な系、非自明な定常状態（非ゼロ速度を伴う）においても、ロバストで効率的な性能を示す。
• 異なる次数の精度および解の種別に対して、well-balancing 変更による計算コスト増加を定量的に評価する。

提案手法

• 解そのものではなく、既知の目標解からの偏差を進化させる形に有限体積スキームを再定式化することで、目標解の正確な維持を保証する。
• 静水圧再構成の原則を任意の目標解に適応し、既知の解の構造に基づいた再構成を実施する。
• 離散レベルでフラックスとソース項が正確にバランスされるように、well-balanced な数値積分戦略を用いてソース項を修正する。
• 任意の整合性のある再構成、数値フラックス、積分則、時間積分スキームと互換性があり、任意の高次精度を実現可能である。
• 構造化直交座標系および極座標系のグリッドに実装可能であり、静的および時間発展する目標解をサポートする。
• 離散的目標解（例：数値シミュレーションから得たもの）に対しては、事前に計算されたデータ配列を用いることで、精度と効率を維持する。

実験結果

リサーチクエスチョン

• RQ1 任意の双曲型バランス法の解（系や解の種別にかかわらず）を正確に維持できる一般化された有限体積フレームワークを設計できるか？
• RQ2 解が単純な静的状態でない場合に、well-balanced スキームで高次精度をどのように維持できるか？
• RQ3 時間依存解や複雑な定常解に対して well-balancing を強制する際の計算コスト増加はどの程度か？
• RQ4 圧縮性状態方程式や非自明な速度場（例：差動回転星）を有する系に対しても、本手法を適用可能か？
• RQ5 well-balancing 変更の計算コストは、次数の増加に伴ってどのように変化するか？

主な発見

• 圧縮性Euler方程式に重力が作用する系や理想MHD系において、複雑な状態方程式や非ゼロ速度場を伴う静的・定常的解を正確に維持する。
• 1次元静的解に対しては、1次および2次スキームでは計算コストが約20%増加したが、7次スキームではわずかに5%増加にとどまった。
• 極座標系のケプラーディスク（定常的・時間不変）に対しては、1次スキームで20%未満、2次スキームで10%未満のコスト増加を示した。
• 時間依存解（摂動を伴うEuler波）に対しては、1次スキームで最大30%のウォールクロックタイム比を示したが、3次スキームでは約15%にとどまった。
• すべての数値実験において、理論的期待と一致する観測収束率を示し、高次精度が維持された。
• 主にソース項と再構成戦略の調整のみで、既存の有限体積コードへの実装が容易であり、最小限のコード変更で実現可能である。