[論文レビュー] The Basics of Weak Galerkin Finite Element Methods

Junping Wang, Xiu Ye|arXiv (Cornell University)|Jan 28, 2019

Advanced Numerical Methods in Computational Mathematics参考文献 33被引用数 3

ひとこと要約

本稿は、弱ガレルキン有限要素法（WG-FEM）とハイブリッド可能不連続ガレルキン（HDG）法の間の概念的・手法的差異を明確にし、両者が一部の数学的根拠を共有しているものの、根本的な哲学的基盤が異なることを示している。著者らは理論的分析と数値例を通じて、同じ有限要素空間を用いても、プライマル WG-FEM が既存の HDG 法と等価でないことを証明した。これは、両者の変分定式化に根本的な不一致が存在することを示している。

ABSTRACT

The goal of this article is to clarify some misunderstandings and inappropriate claims made in [6] regarding the relation between the weak Galerkin (WG) finite element method and the hybridizable discontinuous Galerkin (HDG). In this paper, the authors offered their understandings and interpretations on the weak Galerkin finite element method by describing the basics of the WG method and how WG can be applied to a model PDE problem in various variational forms. In the authors' view, WG-FEM and HDG methods are based on different philosophies and therefore represent different methodologies in numerical PDEs, though they share something in common in their roots. A theory and an example are given to show that the primal WG-FEM is not equivalent to the existing HDG [9].

研究の動機と目的

• WG-FEM と HDG 法の等価性について [6] における誤解を是正すること。
• WG-FEM と HDG の間の哲学的・手法的差異を、共通の数学的基盤があるにもかかわらず明確にすること。
• 理論的および具体例を用いて、プライマル WG-FEM が既存の HDG スキームと等価でないことを示すこと。
• 共通の基底関数フレームワークを用いて、ハイブリダイズド混合 WG-FEM と HDG 定式化を体系的に比較すること。
• WG-FEM が HDG の再表現に過ぎないという主張に反論し、各手法の独自の設計原理を強調すること。

提案手法

• 著者らはモデル偏微分方程式に対して3つの変分形式（プライマル、プライマル-混合、混合＝双対混合）を構築し、これにより WG 法と HDG 法の基礎を築いた。
• 各変分形式に対して、離散的弱勾配と安定化項を用いて、対応する弱ガレルキン FEM（プライマル WG-FEM、プライマル-混合 WG-FEM、混合 WG-FEM）を構築した。
• 混合 WG-FEM のハイブリダイズド定式化を、Fraeijs de Veubeke 法を適用することで得た。これは、ハイブリダイズド混合 FEM に類似している。
• 同一の有限要素空間と基底関数を用いて、ハイブリダイズド混合 WG-FEM を HDG 法と直接比較可能な形に再定式化した。
• 両手法の未知数と支配方程式を、同一の離散化フレームワーク下で一致させることで比較を実施した。
• 係数関数 a(x) = 1 + x を用いた反例を構築し、両手法の双線形形式が一致しないことを示した。これにより、等価性の否定が証明された。

実験結果

リサーチクエスチョン

• RQ1 同一の有限要素空間と基底関数を用いる場合、プライマル弱ガレルキン有限要素法は、既存のハイブリッド可能不連続ガレルキン法と等価であるか？
• RQ2 共通の数学的根拠があるにもかかわらず、弱ガレルキン法と HDG 法の根本的な哲学的・手法的差異は何か？
• RQ3 ハイブリダイズド混合 WG-FEM を、未知数と線形方程式系の観点から HDG 法と直接的・公平に比較可能な形に再定式化できるか？
• RQ4 [6] で提示された HDG の第二の定式化は、プライマル WG-FEM の真の再定式化であるか？ もしそうならば、どのような条件下で成立するか？
• RQ5 『2013 年の WG 法は混合法である』という主張は正確か？ 特に変数係数を有する問題においては？

主な発見

• 同一の有限要素空間と基底関数を用いても、プライマル WG-FEM は既存の HDG 法と等価でない。理論的反例によってこれを示した。
• 弱ガレルキン法と HDG 法の双線形形式は根本的に異なる。係数関数 a(x) = 1 + x の特定のテストケースにおいて、恒等式 (9.11) が成り立たない（3/2 (ln 2)² ≠ ln 2 であるため）。
• [6] で提示された HDG の第二の定式化は、特定の有限要素空間および安定化項パラメータの選択下で、プライマル WG-FEM の再定式化であることが示された。
• 著者らは、WG-FEM が HDG の再表現に過ぎないという主張を否定し、WG-FEM が独自の基礎的原理に基づく明確に異なる数値的アプローチであると主張した。
• [6] における『2013 年の WG 法は混合法である』という主張は誤りであり、特に変数拡散係数 a(x) を有する問題において顕著である。
• ハイブリダイズド混合 WG-FEM と HDG は構造的類似性を有するが、その背後にある定式化および解法メカニズムに差異があり、両者が同一の手法ではないことを確認した。