[論文レビュー] Construction and Analysis of an HDG Method for Incompressible Magnetohydrodynamics

Jeonghun J. Lee, Stephen Shannon+2 more|arXiv (Cornell University)|Feb 16, 2017

Advanced Numerical Methods in Computational Mathematics参考文献 27被引用数 2

ひとこと要約

本稿では、双対 saddle-point 形式から導かれた新規の HDG フラックスを用いて、定常線形化不圧縮性磁気流体力学（MHD）方程式に対するハイブリッド不連続ガレルキン（HDG）法を提案する。この手法は、流体速度および磁場変数に関して最適収束を達成し、他の変数に関しては準最適収束を達成する。これは、事前誤差解析および2次元および3次元の単体メッシュ上での数値実験によって検証されている。

ABSTRACT

We present a hybridized discontinuous Galerkin (HDG) method for stationary linearized incompressible magnetohydrodynamics (MHD) equations. At the heart of the paper is the introduction of an HDG flux of the dual saddle-point form of the MHD equations that facilitates the hybridization of discontinuous Galerkin (DG) method. We carry out the $ extit{a priori}$ error estimates for the proposed HDG method on simplicial meshes in both two- and three-dimensions. The analysis provides optimal convergence for the fluid velocity and the magnetic variables, and quasi-optimal convergence for the remaining quantities. Numerical examples are presented to verify the theoretical findings.

研究の動機と目的

• 定常線形化不圧縮性磁気流体力学（MHD）方程式を安定的かつ高精度に解くための数値法を開発すること。
• 双対 saddle-point 形式における新規 HDG フラックスを導入することで、ハイブリダイズド不連続ガレルキン（HDG）フレームワークを MHD システムへ拡張すること。
• 2次元および3次元の単体メッシュ上における提案された HDG 法の厳密な事前誤差推定を確立すること。
• 2次元および3次元の領域における数値実験を通じて、理論的収束率を検証すること。

提案手法

• 本手法は、線形化 MHD 方程式の双対 saddle-point 形式から導かれたハイブリッド不連続ガレルキン（HDG）式を採用する。
• ハイブリダイゼーションを可能にするために、新規の HDG 数値フラックスを導入し、グローバル自由度を削減しながらも安定性と精度を維持する。
• 手法は2次元および3次元の単体メッシュに適用され、柔軟かつ局所的質量保存性を有する近似を可能にする。
• 要素界面におけるトレース未知数を介してグローバル系を分離することで、効率的な解法戦略が可能になる。
• 関数解析および適合単体三角形分割上の近似理論を用いて、事前誤差推定を導出する。

実験結果

リサーチクエスチョン

• RQ1 線形化不圧縮性 MHD 方程式に対して、ハイブリッド不連続ガレルキン法を効果的に定式化できるか？
• RQ2 2次元および3次元の設定において、HDG 法の流体速度、磁場、および他の変数の収束率はいかほどか？
• RQ3 提案された HDG フラックスは、すべての解成分について最適または準最適収束を保証するか？
• RQ4 理論的誤差推定値は、単体メッシュ上での数値結果とどのように一致するか？

主な発見

• HDG 法は、2次元および3次元の両方において、流体速度および磁場変数に関して最適収束率を達成する。
• 圧力および磁場ポテンシャル変数に関しては、理論的解析で予測された通り、準最適収束が得られる。
• 事前誤差推定は厳密に導出され、単体メッシュ上での数値実験によって検証されている。
• 数値例により、多項式次数やメッシュタイプの異なる条件下でも、理論的収束挙動が確認された。