[論文レビュー] Divergence-Free WENO Reconstruction-Based Finite Volume Scheme for Solving Ideal MHD Equations on Triangular Meshes

Zhiliang Xu, Dinshaw S. Balsara|arXiv (Cornell University)|Oct 4, 2011

Computational Fluid Dynamics and Aerodynamics参考文献 31被引用数 8

ひとこと要約

本稿では、非構造三角メッシュ上で理想磁気流体力学（MHD）方程式を解くための高次精度で発散なしの有限体積スキームを提案する。制約伝送（CT）フレームワークと、新規の重み付きフラックス補間法および発散なしのWENO再構成を組み合わせることで、正確に発散なしの磁界を維持し、3次精度を達成するとともに、衝撃波や乱流を含む複雑なMHDテストケースにおいても頑健であることを示している。

ABSTRACT

In this paper, we introduce a high-order accurate constrained transport type finite volume method to solve ideal magnetohydrodynamic equations on two-dimensional triangular meshes. A new divergence-free WENO-based reconstruction method is developed to maintain exactly divergence-free evolution of the numerical magnetic field. A new weighted flux interpolation approach is also developed to compute the z-component of the electric field at vertices of grid cells. We also present numerical examples to demonstrate the accuracy and robustness of the proposed scheme.

研究の動機と目的

• 非構造三角メッシュ上で高次精度で発散なしの有限体積法を構築すること。これまでの高次手法は構造格子に限定されていた。
• 非構造メッシュ上での離散的発散なし磁界の維持という課題に取り組む。これは、算術平均化の欠如により、セル中心の磁界再構成が非自明であるためである。
• 強い衝撃波や高磁界圧力を含むシミュレーションにおいて、負の圧力を防ぐポジティブ性補正戦略を用いて、物理的一致性を保つ。
• WENO再構成と重み付きフラックス補間を用いて、非構造三角メッシュ上に高次精度で一般化された制約伝送（CT）法を実現する。
• ローター、爆風波、Orszag-Tang渦の問題を含むベンチマークMHD問題において、本手法の精度、安定性、非振動的挙動を実証する。

提案手法

• 制約伝送（CT）フレームワークを採用し、磁界をセル辺に、電界をセル頂点に配置することで、磁界の離散的発散なし進化を強制する。
• セル頂点における電界のz成分を正確にかつ一貫して計算するための新しい重み付きフラックス補間法を導入する。
• 辺上の磁界値からセル中心の磁界を再構成する新規の発散なしWENO再構成法を開発し、共有辺をまたがる法線成分の連続性を保つ。
• 密度、運動量、エネルギーといったセル中心の保存変数に対して、標準のWENO再構成を用いて高次精度を達成する。
• 数値誤差によって負の圧力が生じる領域において、局所的なポジティブ性補正を適用し、エントロピーに基づいて全エネルギー密度を更新することで物理的一致性を維持する。
• 非構造三角メッシュ上に高次空間再構成とルンゲ・クッタ時間積分を組み合わせた有限体積スキームを実装し、保存則と精度を保証する。

実験結果

リサーチクエスチョン

• RQ1 非構造三角メッシュ上での理想MHDシミュレーションにおいて、高次（3次）精度の有限体積スキームが正確に発散なしの磁界を維持できるか？
• RQ2 非構造メッシュ上の制約伝送フレームワークにおいて、セル頂点における電界のz成分を正確かつ一貫して計算する方法は何か？
• RQ3 非構造三角グリッド上での共有辺をまたがる磁界法線成分の一貫性を保証する発散なしWENO再構成法を構築できるか？
• RQ4 提案手法は、衝撃波や乱流のような複雑なMHD流れにおいて、物理的圧力の正の性質を維持し、不適切な振動を避けることができるか？
• RQ5 本手法は3次元や一般の非構造メッシュへどの程度拡張可能か？

主な発見

• 提案された2次および3次精度スキームは滑らかな流れのテストで期待される収束次数を達成しており、3次スキームは不連続性の解像度が優れていることが示された。
• 発散なしWENO再構成により、共有辺をまたがる磁界法線成分が連続的であることが保証され、機械精度に近い離散的発散なし条件が維持された。
• 重み付きフラックス補間法により、CTフレームワーク内での磁界発散なし性を維持するために不可欠な頂点における電界の正確かつ一貫した計算が可能になった。
• 高磁界圧力・低ガス圧力の爆風波問題において、局所的エントロピーに基づくポジティブ性補正により負の圧力を効果的に防止し、物理的一致性を維持した。
• Orszag-Tang渦のシミュレーションでは、非振動的かつ複雑な乱流的流れ構造を正しく捉え、領域全体で正の密度および圧力を維持した。
• ローター問題のシミュレーションでは、密度に劣化が生じず、ポジティブ性補正の適用も不要であり、強い回転流においても頑健であることが示された。