[論文レビュー] Low Regularity Primal-Dual Weak Galerkin Finite Element Methods for Convection-Diffusion Equations

Chunmei Wang, Ludmil Zikatanov|arXiv (Cornell University)|Jan 20, 2019

Advanced Numerical Methods in Computational Mathematics参考文献 29被引用数 5

ひとこと要約

本稿は、低正則性仮定の下で対流拡散方程式に対して、離散弱勾配を介して原始変数と随伴方程式を結合する混合変分形式を用いた、原始的双対弱ガラーキン（PDWG）有限要素法を提案する。この手法は、$\epsilon \in [0, \frac{1}{2})$ の $H^\epsilon$-ノルムにおいて最適収束を達成し、不連続解や対流優勢領域でも安定性と精度を示す。

ABSTRACT

We propose a numerical method for convection-diffusion problems under low regularity assumptions. We derive the method and analyze it using the primal-dual weak Galerkin (PDWG) finite element framework. The Euler-Lagrange formulation resulting from the PDWG scheme yields a system of equations involving not only the equation for the primal variable but also its adjoint for the dual variable. We show that the proposed PDWG method is stable and convergent. We also derive a priori error estimates for the primal variable in the $H^ε$-norm for $ε\in [0,\frac12)$. A series of numerical tests that validate the theory and are presented as well.

研究の動機と目的

• 対流優勢な状況下で、標準ガラーキン法が示す数値的不安定性や振動を解消すること。
• 正確な解が低正則性（例えば、不連続性や制限されたソボレフ滑らかさ）を示す場合でも、安定かつ正確に保てる強固な有限要素法の開発。
• 原始変数の $H^\epsilon$-ノルムにおける事前誤差推定を、$\epsilon \in [0, \frac{1}{2})$ の範囲で提示し、滑らかさが低い解への適用可能性を拡張すること。
• 幾何的柔軟性を高めるために、一般の多面体メッシュ（非単体要素を含む）を用いることを可能にすること。
• 特定の状況下では安定化パラメータ $\gamma$ を省略可能であり、実装の堅牢性と簡潔さが向上することを示唆すること。

提案手法

• 原始変数と随伴変数を含む混合変分形式に基づく、オイラー＝ラグランジュ方程式を用いて対流拡散問題を定式化する。
• 低正則性解を扱えるように、多面体メッシュ上で作用する離散弱勾配を導入する原始的双対弱ガラーキン（PDWG）フレームワークを採用する。
• $C^0$-連続有限要素空間を用いる：原始変数には $P_k(T)$、そのトレースには $P_{k-1}(\partial T)$、随伴変数には $P_s(T)$（$s = k-1$ または $s = k-2$）を設定する。
• 原始方程式と随伴方程式を結合することで安定性と可解性を向上させる安定化項を組み込み、対称的かつ適切に定義された離散系を構築する。
• 非単体要素を含む一般の多面体分割に適用可能であり、メッシュ設計の柔軟性を高める。
• 関数解析および近似理論を用いて誤差推定を導出する。正確な解の正則性が最小限である条件下で、$\epsilon \in [0, \frac{1}{2})$ の $H^\epsilon$-ノルムに注目する。

実験結果

リサーチクエスチョン

• RQ1 正確な解が $H^\epsilon$-正則性（$\epsilon < \frac{1}{2}$）を示すような場合でも、有限要素法が対流拡散問題に対して安定かつ収束性を保つことができるか？
• RQ2 PDWGフレームワークは、このような低正則性条件下でも $H^\epsilon$-ノルムにおける最適誤差推定を可能にするか？
• RQ3 $s = k-1$ の場合、安定化パラメータ $\gamma$ を必要とせずに、収束性と精度を維持できるか？
• RQ4 特異点や境界層を有する複雑な領域（例：亀裂あり領域や L 字型領域）において、数値的にどのように性能を発揮するか？
• RQ5 PDWG法は、不連続または弱い正則性を示す解を有する対流優勢問題に対して、最適収束率を維持しながらどれほど効果的に処理できるか？

主な発見

• 提案されたPDWG法は、$\epsilon \in [0, \frac{1}{2})$ の $H^\epsilon$-正則性を有する解を含む、低正則性仮定下でも安定かつ収束性を示す。
• 正確な解の正則性が最小限である条件下で、原始変数の $H^\epsilon$-ノルムにおける事前誤差推定が導出され、最適収束次数が確認された。
• 数値実験により、理論的予測以上または同等の性能を示し、対流優勢領域でも一貫した結果が得られた。
• 安定化パラメータ $\gamma = 0$ の場合でも収束性と精度が維持されることから、特定の条件下では安定化項が不要である可能性が示唆された。
• 亀裂や L 字型領域などの複雑な領域において、境界層や不連続性を正確に捉え、正確な等高線図および表面プロットを生成した。
• 多面体メッシュと $C^0$-連続有限要素の使用により、幾何的柔軟性が向上した一方で、収束性や安定性は損なわれなかった。