[論文レビュー] A Hybridized Weak Galerkin Finite Element Method for the Biharmonic Equation

Chunmei Wang, Junping Wang|arXiv (Cornell University)|Feb 5, 2014

Numerical methods in engineering被引用数 23

ひとこと要約

本稿では、連続性制約を緩和するために要素境界にラグランジュ乗数を組み込むことで、ビアーモンド方程式に対するハイブリッド化された弱ガレルキン有限要素法（HWG）を提案する。この手法は最適な誤差推定を達成し、内部要素自由度を削除することでシュール補行列削減を可能にし、システムサイズを著しく縮小する。これにより、対称正定値の低次元化されたシステムが得られ、効率的な解法アルゴリズムが実現される。

ABSTRACT

This paper presents a hybridized formulation for the weak Galerkin finite element method for the biharmonic equation. The hybridized weak Galerkin scheme is based on the use of a Lagrange multiplier defined on the element boundaries. The Lagrange multiplier is verified to provide a numerical approximation for certain derivatives of the exact solution. An optimal order error estimate is established for the numerical approximations arising from the hybridized weak Galerkin finite element method. The paper also derives a computational algorithm (Schur complement) by eliminating all the unknown variables on each element, yielding a significantly reduced system of linear equations for unknowns on the boundary of each element.

研究の動機と目的

• ビアーモンド方程式に対する弱ガレルキン有限要素法のハイブリッド化定式化を構築し、計算効率を向上させること。
• 要素境界にラグランジュ乗数を導入し、正確な解の微分を近似すること。
• HWGフレームワーク下で数値近似の最適順序誤差推定を確立すること。
• 内部未知数を削除することでシュール補行列系を導出し、グローバルシステムのサイズを低減すること。
• 変数削減と対称正定値低次元化線形システムにより、並列計算が可能な効率的な計算を可能とすること。

提案手法

• 解、そのトレース、要素境界における勾配の多項式近似 $P_k/P_{k-2}/P_{k-2}$ を用いた弱ガレルキン有限要素定式化を採用する。
• 要素界面にラグランジュ乗数を導入し、連続性を弱く強制することでハイブリダイゼーションを実現する。
• シュール補行列技術を用いて各要素の内部未知数をすべて削除し、グローバルシステムを要素境界上の未知数のみに還元する。
• 境界変数のみを含む低次元化された線形システムを導出し、これは対称かつ正定値であるため、効率的な反復解法が可能になる。
• グローバルシステムを2段階のアルゴリズムで解く：まず低次元化されたシュール補行列系を解き、その後各要素で内部値を再構築する。
• 離散的弱ヘッシアンと局所射影作用素を用いて変分定式化を定義し、安定性を確保する。

実験結果

リサーチクエスチョン

• RQ1 ハイブリダイゼーションは、ビアーモンド方程式に対する弱ガレルキン有限要素法に効果的に適用可能であり、計算コストを低減できるか？
• RQ2 要素境界に導入されたラグランジュ乗数は、解およびその微分の安定的かつ収束的数値近似を提供するか？
• RQ3 シュール補行列技術は、ハイブリッド化された弱ガレルキンスキームに成功裏に適用可能であり、大幅に未知数が少ない低次元化システムをもたらすか？
• RQ4 ハイブリッド化された弱ガレルキン法のビアーモンド方程式に対する最適収束速度は何か？
• RQ5 得られた低次元化システムは、特にプリコンディショニングを用いて実際の計算で効率的に解けるか？

主な発見

• ハイブリッド化された弱ガレルキン法は、解およびその微分の両方について最適順序誤差推定を達成しており、理論的収束速度が確認された。
• 要素境界に位置するラグランジュ乗数は、正確な解の法線微分の数値近似を提供し、勾配情報の高精度な回復を可能にする。
• シュール補行列系により、グローバル未知数の数が内部および外部の辺上にのみ制限され、計算コストが著しく低減された。
• 低次元化されたシステムは対称かつ正定値であるため、安定性が保証され、プリコンディショニングを用いた効率的な反復解法が利用可能になる。
• 計算アルゴリズムにより、各要素における局所問題の並列計算が可能となり、スケーラビリティとパフォーマンスが顕著に向上した。
• ハイブリダイゼーションにより、元の弱ガレルキン定式化の整合性と安定性を保持しつつ、計算効率を大幅に向上させた。