[論文レビュー] Robust residual-based a posteriori error estimators for mixed finite element methods for fourth order elliptic singularly perturbed problems

Shaohong Du, Runchang Lin+1 more|arXiv (Cornell University)|Sep 15, 2016

Advanced Numerical Methods in Computational Mathematics参考文献 14被引用数 6

ひとこと要約

本稿では、小さな摂動パラメータ ε を持つ4階の特異摂動楕円型問題に適用された混合有限要素法に対する、ロバストな残差ベースの後験的誤差推定器を提案する。著者らは、εに依存する項と依存しない項をバランスさせる新しい誤差測度を提案し、その結果、すべての ε に対して一様に有効である推定器を導出する。数値例により、境界層の効果的な解像度と最適収束率の達成が検証された。

ABSTRACT

We consider mixed finite element approximation of a singularly perturbed fourth-order elliptic problem with two different boundary conditions, and present a new measure of the error, whose components are balanced with respect to the perturbation parameter. Robust residual-based a posteriori estimators for the new measure are obtained, which are achieved via a novel analytical technique based on an approximation result. Numerical examples are presented to validate our theory.

研究の動機と目的

• 小さな ε を持つ4階の特異摂動問題における混合有限要素法のロバストな後験的誤差推定器の設計という課題に取り組む。
• 摂動パラメータ ε と解の正則性の寄与をバランスさせる新しい誤差測度を同定する。
• ε → 0 の極限においてもロバスト（一様に有効）である残差ベースの推定器を開発し、標準的な推定器が ε が小さいと性能を落とすという限界を克服する。
• 理論的推定器を数値的に検証し、境界層を捉え、最適収束率を達成できることを示す。

提案手法

• エネルギーノルム ‖∇ψ‖ と u の H¹-セミノルムを組み合わせた新しい誤差測度 (ε²‖∇ψ‖² + ‖ψ‖² + |u|₁²)¹ᐟ² を導入し、εに依存する成分と依存しない成分をバランスさせる。
• 近似結果に基づく新しい解析的手法を用いて、新しい誤差測度に対する残差ベースの後験的誤差推定器を導出する。
• 要素毎の推定子 ηψ,T と ηu,T に基づく適応的有限要素メッシュの細分化を実施し、ψ における境界層を解像する。
• 効果性指数を用いて推定器の信頼性を検証し、さまざまな ε 値において推定誤差が真の誤差に近いことを保証する。
• 既知または期待される境界層挙動を示す2つのベンチマーク問題に推定器を適用し、マークパラメータ θ を用いた適応的メッシュ細分化を実施する。
• 収束率が ε に依存しないことを示すことで、推定器のロバスト性を検証する。特に、エネルギーノルムにおいて O(h¹ᐟ²) の収束率を示す。

実験結果

リサーチクエスチョン

• RQ1 4階の特異摂動問題において、解の ε に依存する成分と依存しない成分をバランスさせる新しい誤差測度を構築できるか？
• RQ2 この新しい測度から導出された残差ベースの後験的推定器は、摂動パラメータ ε に対してロバストか？
• RQ3 提案された推定器は、ε が非常に小さい場合でも、フラックス ψ = −Δu における鋭い境界層を効果的に捉え、解像できるか？
• RQ4 推定器は、ε ≈ 10⁻⁷ を含むさまざまな ε 値においても最適収束率を維持するか？
• RQ5 推定器の効果性指数は、すべての ε に対して一様に有界か？ これは、誤差推定の信頼性を示唆する。

主な発見

• 提案された誤差測度 (ε²‖∇ψ‖² + ‖ψ‖² + |u|₁²)¹ᐟ² は、ε と解の正則性の寄与を効果的にバランスさせ、層の挙動を捉えられない不均衡なノルムを回避する。
• 残差ベースの後験的誤差推定器は ε に対してロバストであり、テストされたすべての ε 値（ε = 10⁻⁷ を含む）において一様に信頼性を保つ。
• 数値実験により、推定器はエネルギーノルムにおいておおよそ O(h¹ᐟ²) の最適収束率を達成することが示され、ε に依存しない。
• 推定子に基づく適応的メッシュ細分化により、ψ における4つの鋭い境界層が成功裏に解像された。視覚的確認とメッシュ細分化のパターンから裏付けられた。
• ‖ψ − ψh‖ℰ の推定器の効果性指数は、すべての ε 値において 1 に近い（例：0.99 から 1.13 の間）ままであり、推定器の信頼性を確認する。
• ε = 10⁻⁵ の場合、22回の細分化反復後に最小メッシュサイズ h_min(ε) が 6.10×10⁻⁵ に達し、境界層が ε スケールで解像されたことが示された。