[論文レビュー] $L^2$-error estimates for finite element approximations of boundary fluxes

Mats G. Larson, André Massing|arXiv (Cornell University)|Jan 27, 2014

Advanced Numerical Methods in Computational Mathematics参考文献 19被引用数 9

ひとこと要約

本稿では、$L^2$-Dirichletデータを備えた離散双対問題に基づく新しい誤差表現式を用いて、境界法線フラックスの有限要素近似に対する準最適な $L^2$-誤差推定値を提示する。本手法は、境界近傍における重み付き離散安定性と異方的補間推定を組み合わせることで、Nitsche法や安定化ラグランジュ乗数法を含むさまざまな弱い強制スキームにおいて、分片的線形要素で $|\text{ln}\hspace{1pt} h|h$ 収束を達成する。

ABSTRACT

We prove quasi-optimal a priori error estimates for finite element approximations of boundary normal fluxes in the $L^2$-norm. Our results are valid for a variety of different schemes for weakly enforcing Dirichlet boundary conditions including Nitsche's method, and Lagrange multiplier methods. The proof is based on an error representation formula that is derived by using a discrete dual problem with $L^2$-Dirichlet boundary data and combines a weighted discrete stability estimate for the dual problem with anisotropic interpolation estimates in the boundary zone.

研究の動機と目的

• 楕円問題における境界法線フラックスの有限要素近似の事前 $L^2$-誤差推定値を導出すること。
• モルタル法を超えて、Nitsche法や安定化ラグランジュ乗数スキームを含む、既存の収束結果を拡張すること。
• 鋭い誤差表現と安定性解析を可能にする、$L^2$-Dirichletデータを備えた離散双対問題フレームワークを確立すること。
• Besov空間フレームワークを避ける代わりに、離散双対問題における安定性推定に解析的手法を用いること。
• マニュフェイクチャドソリューションを用いた数値実験により、理論的収束率を検証すること。

提案手法

• 離散双対問題の解を用いた誤差表現式を用いて、境界フラックス誤差を定式化する。ただし、$L^2$-境界データを備える。
• エネルギーノルムにおける双対問題の重み付き離散安定性推定を確立し、境界近傍における勾配と $L^2$-ノルムを制御する。
• 境界層における異方的補間推定を用いて、主問題の解の補間誤差を抑え込む。
• 変分式において安定化項を導入し、境界項の影響を制御し、弱い強制下での安定性を保証する。
• 離散安定性推定と補間誤差バウンドを組み合わせることで、境界フラックスの最終的な $L^2$-誤差推定値を導出する。
• 均一な三角形メッシュ上でのマニュフェイクチャドソリューションを用いた数値実験により、理論的結果を検証する。$CG(1)$ 要素と $DG(0)$ 乗数を用いる。

実験結果

リサーチクエスチョン

• RQ1 Nitsche法や安定化ラグランジュ乗数スキームを含む、広範な弱いDirichlet条件のクラスに対して、境界フラックスの準最適な $L^2$-誤差推定値を導出できるか。
• RQ2 分片的線形有限要素を用いた場合、境界上の $L^2(\text{boundary})$-ノルムにおける離散境界フラックス近似の収束率は何か。
• RQ3 Besov空間フレームワークに依存せずに、$L^2$-Dirichletデータを備えた離散双対問題を用いて、フラックス誤差を表現・バウンドする方法は何か。
• RQ4 境界層における異方的補間推定と重み付き安定性が、最適収束率を達成するために果たす役割は何か。
• RQ5 数値実験は、さまざまなフラックス近似スキームにおいて、理論的 $|\text{ln}\hspace{1pt} h|h$ 収束率が漸近的領域で確認されるか。

主な発見

• 本稿では、分片的線形有限要素を用いた境界法線フラックス近似に対して、$|\text{ln}\hspace{1pt} h|h$ 階の準最適 $L^2$-誤差推定値を確立した。
• 誤差推定値は、Nitsche法や安定化ラグランジュ乗数スキームを含む、複数の弱いDirichlet条件強制法に対して成立する。
• 離散双対問題がエネルギーノルムにおいて重み付き安定性推定を満たすことが示され、これはフラックス誤差のバウンドに不可欠である。
• 境界層における異方的補間推定は、双対問題における境界層の誤差寄与を制御するために不可欠である。
• 数値結果は、漸近的領域における理論的収束率を確認しており、Nitsche法および安定化ラグランジュ乗数法の両方が1次収束を達成している。
• 前漸近的挙動では最適でない収束率が観察されるが、$h \to 0$ に伴いその影響は消滅し、理論的漸近的推定値が妥当であることが検証された。