[論文レビュー] Error estimates for structure-preserving discretization of the incompressible MHD system

Yicong Ma, Jinchao Xu+1 more|arXiv (Cornell University)|Aug 10, 2016

Advanced Numerical Methods in Computational Mathematics参考文献 25被引用数 4

ひとこと要約

本稿では、エネルギー安定性および物理的保存則を活用して、非圧縮性磁気流体力学（MHD）系の構造保存型有限要素離散化に対する最適な誤差推定値を提示する。適切なCFL条件下では、空間および時間方向で最適な収束次数を達成し、数値実験により理論的予測と一致する収束率が観測された。

ABSTRACT

In this paper, we carry out the error analysis for the structure-preserving discretization of the incompressible MHD system. This system, as a coupled system of Navier-Stokes equations and Maxwell's equations, is nonlinear. We use its energy estimate and the underlying physical structure to facilitate the error analysis. Under certain CFL conditions, we prove the optimal order of convergence. To support the theoretical results, we also present numerical tests.

研究の動機と目的

• 非圧縮性MHD系に適用された構造保存型有限要素法に対する厳密な誤差推定値を確立すること。
• 電場の磁束密度に関するガウスの法則を含む、重要な物理法則を保存するH(curl)およびH(div)適合有限要素を用いた混合有限要素スケームの収束挙動を分析すること。
• 非線形かつ時間に依存するMHD方程式に対して、特定の時間刻み幅およびメッシュサイズの制約下で最適な収束次数を達成すること。
• 製造された解を用いた数値実験と収束率解析を通じて、理論的結果を検証すること。

提案手法

• 本手法は、有限要素外微分計算（FEEC）および離散外微分計算（DEC）に基づく混合有限要素式を採用し、電場にはH(curl)-適合のネデリュク要素、磁場にはH(div)-適合のラヴィアル・トーマス要素を用いる。
• スケームは、磁場の発散がゼロである条件を離散レベルで正確に保持し、ガウスの法則に準拠する。
• 半陽解的時間離散化を用い、非線形対流項は陽的に、線形項は陰的に扱い、CFL条件の下で安定性を維持する。
• MHD系の物理的構造（全エネルギー保存）を基盤としたエネルギー推定値を用いて誤差解析を実施する。
• Sobolevノルムにおける事前推定値および有限要素空間の近似性質に依拠し、速度および磁場勾配の有界性を用いて非線形項を慎重に取り扱う。
• 進化型サドルポイント問題の収束フレームワークを応用し、構造保存性を有する非線形かつ時間に依存する系へ既存理論を拡張する。

実験結果

リサーチクエスチョン

• RQ1 時間および空間離散化の下で、非圧縮性MHD系に対する構造保存型有限要素離散化の最適収束速度は何か？
• RQ2 時間刻み幅 $ k $ およびメッシュサイズ $ h $ が、速度、磁場、電場、圧力変数の収束次数にどのように影響を与えるか？
• RQ3 エネルギー安定性および物理的保存則（例：$ \nabla \cdot \bm{B} = 0 $）を活用して、非線形かつ時間に依存するMHD系に対する厳密な誤差推定値を導出可能か？
• RQ4 提案手法は妥当なCFL条件下で最適な収束次数を達成するか？また、既存手法と比較してどう異なるか？
• RQ5 数値実験は、誤差解析が予測する理論的収束率をどの程度確認できるか？

主な発見

• 時間刻み幅 $ k \leq h^{1/2} $ の条件下で、速度、磁場、電場はそれぞれのエネルギーノルムにおいて最適な収束次数 $ \mathcal{O}(k^2 + h^{2\hat{r}}) $ を達成する。
• 圧力に関しては、$ k \leq h^{3/2} $ の条件下で最適な収束次数 $ \mathcal{O}(k^2 + h^{2r}) $ が達成され、圧力の精度にはより厳しい時間刻み幅制約が課されることが示された。
• 電場の誤差は $ \| \bm{E} - \bm{E}_h \|^2 + k \| \nabla \times (\bm{E} - \bm{E}_h) \|^2 \leq C(k^2 + h^{2\hat{r}}) $ で有界であり、同じCFL条件下で最適収束が確認された。
• 数値実験により、予測された収束率が確認された：観測された収束次数は、メッシュサイズ $ h $ および時間刻み幅 $ k $ に対して概ね1であり、理論的解析と一致した。
• スケームは、離散レベルで磁場の発散がゼロである条件を正確に保持しており、これが解の安定性および精度向上に寄与している。
• 理論的誤差境界は、正確な解の正則性に依存し、$ \| \bm{u} \|_{0,\infty} $、$ \| \nabla \bm{u} \|_{0,3} $、$ \| \bm{B} \|_{0,\infty} $、$ \| \bm{j} \|_{0,\infty} $ の有界性、および領域の幾何構造を含む。