[論文レビュー] Positivity-Preserving Finite Difference WENO Schemes with Constrained Transport for Ideal Magnetohydrodynamic Equations

Andrew Christlieb, Yuan Liu+2 more|arXiv (Cornell University)|Jun 19, 2014

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

ひとこと要約

本稿では、離散的発散自由磁場を維持するための制約伝搬（CT）を用いた、理想磁流体力学（MHD）の高次陽性保存有限差分WENOスキームを提示する。高次WENO-HCLのフラックスに最大原理保存技術に基づくパラメータ化されたフラックスリミッターを適用することにより、CFL安定性限界を低下させることなく、密度および圧力の陽性を保証する。1次元、2次元、3次元のシミュレーションにおいて、高次精度、発散自由磁場、物理的解の陽性を同時に達成する。

ABSTRACT

In this paper, we utilize the maximum-principle-preserving flux limiting technique, originally designed for high order weighted essentially non-oscillatory (WENO) methods for scalar hyperbolic conservation laws, to develop a class of high order positivity-preserving finite difference WENO methods for the ideal magnetohydrodynamic (MHD) equations. Our schemes, under the constrained transport (CT) framework, can achieve high order accuracy, a discrete divergence-free condition and positivity of the numerical solution simultaneously. Numerical examples in 1D, 2D and 3D are provided to demonstrate the performance of the proposed method.

研究の動機と目的

• 低βプラズマ状態下でも解の陽性（密度および圧力）を保つ高次有限差分WENOスキームを、理想MHDに開発すること。
• 時間発展中に離散的発散自由磁場を維持し、非ゼロの発散による数値不安定化を回避すること。
• ベーススキームのCFL安定性限界を損なわせることなく、陽性保存を制約伝搬（CT）と統合する unstaggered フレームワークを構築すること。
• 元々スカラー双曲型保存則に用いられた最大原理保存フラックスリミッティング技術を、多成分かつ非線形な理想MHD方程式系に拡張すること。
• 強い衝撃波、爆発波、渦構造を含む、1次元、2次元、3次元のベンチマーク問題において、提案スキームの頑健性および精度を実証すること。

提案手法

• 数値フラックスの高次WENO再構成に加え、陽性保存のための第一階Lax-Friedrichsフラックスを基本構成要素として用いる。
• Runge-Kutta時間積分の最終段階でパラメータ化されたフラックスリミッターを適用し、最終フラックスのみを変更することで、陽性を強制しつつ、ベーススキームの精度および安定性を維持する。
• 制約伝搬（CT）フレームワークにより、高次有限差分を用いて磁場ポテンシャルを発展させ、離散的に磁場が発散自由であることを保証する。
• 磁場ポテンシャルの発展による発散クリーニング手法を採用し、空間微分をWENOを用いて再構成することで、高次精度を維持する。
• 陽性保存リミッターは、ZhangとShuの圧縮性Euler方程式に対するフレームワークを基盤とし、局所解状態に基づくフラックス制限をMHDに適応したものである。
• 構造化Cartesianメッシュを用いて1次元、2次元、3次元に実装し、流出境界条件を適用し、フラックスには標準的なWENO再構成を用いる。

実験結果

リサーチクエスチョン

• RQ1 理想MHDの高次有限差分WENOスキームは、CFL安定性限界を低下させることなく、密度および圧力の陽性を維持できるか？
• RQ2 制約伝搬（CT）フレームワークを、多次元MHDシミュレーションにおける高次陽性保存フラックスリミッティングと効果的に統合できるか？
• RQ3 提案されたフラックスリミッターは、基礎となるWENOスキームの形式的精度を損なわず、物理的解の境界を保証するか？
• RQ4 爆発波、衝撃管、アルベール波パルスといった複雑なMHD現象を、最小限の発散誤差で解像できるか？
• RQ5 特に3次元において、低βプラズマ状態下で負圧力や負密度が生じる可能性があるが、そのような数値的崩壊を回避できるか？

主な発見

• 提案スキームは、1次元、2次元、3次元の全テストケースにおいて、密度および圧力の陽性を成功裏に維持した。特に球状圧力パルスを伴う3次元爆発波では、リミッターなしでは負圧力が観測された。
• 発散誤差は全シミュレーションで切捨て誤差レベル（O(h^p)）に保たれ、制約伝搬フレームワークが∇·B = 0を効果的に維持していることを確認した。
• 2次元および3次元の爆発問題では、相対的全エネルギー誤差が時間経過に伴い1%未満に保たれ、良好な保存性を示した。
• 滑らかな領域では高次精度を達成し、WENO再構成の形式的次数（例：3次または5次）と整合する収束率を示した。
• CFL制限は、ベースの第一階Lax-Friedrichsスキームと同一であり、陽性リミッターが安定性制限を追加で導入していないことを示した。
• 150³メッシュを用いた3次元爆発問題では、密度および圧力に球対称性が観測され、文献に登録された他の高次スキームと一致する結果を得た。