[論文レビュー] Convergence and Optimal Complexity of Adaptive Finite Element Methods

Lianhua He, Aihui Zhou|arXiv (Cornell University)|Feb 4, 2010

Advanced Numerical Methods in Computational Mathematics参考文献 27被引用数 7

ひとこと要約

本稿は、非対称、非線形、および有界でない係数固有値問題を含む広範な楕円型PDEのクラスに適用される適応型有限要素法（AFEM）の収束性と最適な複雑度を、それらを対称楕円型問題に結びつける摂動フレームワークを通じて確立する。主な貢献は、3つのモデル問題における数値実験によって検証された一般化された収束性および複雑度理論である。

ABSTRACT

In this paper, we study adaptive finite element approximations in a perturbation framework, which makes use of the existing adaptive finite element analysis of a linear symmetric elliptic problem. We prove the convergence and complexity of adaptive finite element methods for a class of elliptic partial differential equations. For illustration, we apply the general approach to obtain the convergence and complexity of adaptive finite element methods for a nonsymmetric problem, a nonlinear problem as well as an unbounded coefficient eigenvalue problem.

研究の動機と目的

• 対称問題に限らない広範な楕円型PDEクラスに適用可能な、適応型有限要素法の一般化された収束性および複雑度理論の構築を目的とする。
• 既存のAFEM文献において非対称および非線形楕円型問題の複雑度解析が欠落している問題に対処すること。
• 特異性と非コンパクト性のため困難である有界でない係数固有値問題に対し、最適な複雑度結果を拡張すること。
• 適応メッシュの改良戦略における後騒ぎ誤差推定器の効率性と信頼性を、代表的な問題における数値実験を通じて検証すること。
• 後騒ぎ誤差推定器が、適応メッシュ生成戦略において効果的かつ信頼性を持って機能することを示すこと。

提案手法

• 一般楕円型問題を基準となる対称線形楕円型問題に結びつける摂動的議論を構築し、既知の収束性および複雑度結果の移転を可能にする。
• 摂動フレームワークを用いて、非対称および非線形問題における準直交性および後騒ぎ誤差推定器の信頼性・有効性を導出する。
• 各適応ループにおける誤差および推定器の低減を証明することで、適応型有限要素法の収束性を確立する。
• 必要な精度に対して自由度の増加が最適な割合で行われることを示すことにより、最適な複雑度を導出する。
• 一般理論を3つの具体的な問題に適用する：非対称PDE、非線形PDE（Poisson-Boltzmann型）、有界でない係数固有値問題（Kohn-Sham方程式）。
• 線形および二次有限要素を用いた適応型有限要素スキームを実装し、誤差指標がメッシュの細分化を制御する。

実験結果

リサーチクエスチョン

• RQ1 非対称楕円型問題において、先行研究で複雑度解析が欠落していた状況でも、AFEMの収束性および複雑度を確立できるか？
• RQ2 一般化された摂動フレームワークを用いて、最適なレートで非線形楕円型問題の収束性および複雑度結果を拡張できるか？
• RQ3 量子化学で生じるような有界でない係数を有する固有値問題に対しても、理論が適用可能か？
• RQ4 これらの多様な問題クラスにおいて、後騒ぎ誤差推定器は実際の応用でどのように機能するか？
• RQ5 適応メッシュは、自由度と誤差の観点から最適な収束レートを達成するか？

主な発見

• 提示された摂動に基づくフレームワークは、非対称、非線形、および有界でない係数固有値問題に適用された適応型有限要素法の収束性と最適な複雑度を成功裏に証明した。
• 数値結果から、誤差 $ parallel u - u_h \nparallel_{1}$ が後騒ぎ誤差推定器に比例するレートで収束することが示され、推定器の有効性が確認された。
• 線形有限要素では、誤差の収束レートが最適な傾き $-1/3$ に一致し、二次要素では $-2/3$ に一致しており、最適な近似性質が裏付けられた。
• 適応的二次有限要素を用いたヘリウム原子の基底状態エネルギーは、100,000未満の自由度で $-2.8346$ a.u.（参考値）に達した。
• 適応メッシュは、断面図で示されるように特異性および境界層付近に集中して細分化され、誤差の局所化が有効に実現された。
• 3つの数値例すべてにおいて誤差と推定器の比例関係が観察されたことから、後騒ぎ誤差推定器は有効かつ信頼性があることが裏付けられた。