[論文レビュー] A Multilevel Correction Type of Adaptive Finite Element Method for Eigenvalue Problems

Hehu Xie|arXiv (Cornell University)|Jan 11, 2012

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

ひとこと要約

本稿では、固有値問題を解くための多層補正に基づく適応有限要素法（AFEM）を提案する。計算コストを著しく削減するため、粗いメッシュ上で小スケールの固有値問題を解き、細かいメッシュ上で関連する境界値問題を解く。この手法は最適収束率と計算量を達成し、境界値問題に対する標準的なAFEMと同等の効率性を示す。

ABSTRACT

A type of adaptive finite element method for the eigenvalue problems is proposed based on the multilevel correction scheme. In this method, adaptive finite element method to solve eigenvalue problems involves solving associated boundary value problems on the adaptive partitions and small scale eigenvalue problems on the coarsest partitions. Hence the efficiency of solving eigenvalue problems can be improved to be similar to the adaptive finite element method for the associated boundary value problems. The convergence and optimal complexity is theoretically verified and numerically demonstrated.

研究の動機と目的

• 2階楕円型固有値問題を、計算コストを低減させつつ効率的に解くための適応有限要素法を開発すること。
• 高価な固有値問題の解法を適応的メッシュ細分化プロセスから分離するために、多層補正フレームワークを活用すること。
• 提案されたAFEMの理論的収束性と最適計算量を、固有値問題の文脈で確立すること。
• 数値実験を通じて、特異点を有する領域に対しても、固有値および固有関数の両方で最適収束率が達成されることを示すこと。
• 非凸領域および不規則領域における数値実験を通じて、標準的直接AFEMと比較して、本手法の効率性を検証すること。

提案手法

• 本手法は、解を2つの成分に分離する多層補正スキームを採用する。1つは最も粗いメッシュ上での小スケール固有値問題、もう1つは細かいメッシュ上での関連境界値問題の適応的有限要素解法である。
• 適応的グリッド細分化は、残差項および振動項に基づく後験的誤差推定子によって駆動され、固有関数近似における高誤差領域に注目して細分化が行われる。
• 誤差推定子に基づき、パラメータ $\theta \in (0,1)$ を用いたマーク戦略により、細分化対象の要素が選択され、最適収束性が保証される。
• 理論的分析は、先行研究（例：[20, 27]）で確立された、固有値問題の有限要素近似とその関連境界値問題との関係に依存する。
• 収束性および最適性は、誤差推定子の準最適性および信頼性を含む、適応的有限要素理論の技術を用いて証明される。
• 数値実験では、$L^2$-ノルム正規化された固有関数が用いられ、収束をモニタリングするために後験的誤差推定子 $\eta_h(u_h, \Omega)$ が計算される。

実験結果

リサーチクエスチョン

• RQ1 多層補正戦略を適応的有限要素法による固有値問題の解法に効果的に適用することで、計算コストを低減できるか？
• RQ2 提案手法は、境界値問題に対する標準的AFEMと同等の最適収束率を、固有値および固有関数両方で達成できるか？
• RQ3 固有関数が特異性を示す再入角を持つ領域では、標準的線形要素が最適収束性を達成できないが、本手法はそのような領域でも有効に機能するか？
• RQ4 最小固有値だけでなく、複数の固有値ペアを解く場合でも、最適計算量を維持できるか？
• RQ5 収束速度および誤差低減の観点から、直接AFEMと比較して、本アルゴリズムの性能はどの程度か？

主な発見

• 提案された多層補正AFEMは、固有値および固有関数の両方で最適収束率を達成し、境界値問題に適用された標準的AFEMと同等の誤差減少率を示す。
• $L$-形状領域における数値実験では、再入角に特異性を有するが、最初の固有値近似が理論的予測と整合する収束率を示す。
• 例2では、後験的誤差推定子 $\eta_h(u_h, \Omega)$ と固有値誤差の両方が、20回の適応的反復において最適収束性を示し、理論的分析の妥当性が確認された。
• 例3では、変数係数の拡散行列と非定数ポテンシャルを有する場合でも、最初の固有値ペアおよび5つの最小固有値ペアに対して、最適収束性が達成された。
• 適応的アルゴリズムCと直接AFEMの比較により、提案手法は収束速度が直接AFEMと同等またはそれを上回り、特に固有値および固有関数誤差の低減において効率的であることが示された。
• 合計自由度の総数が望ましい精度に対して最適な割合で増加することから、本手法は最適計算量を維持しており、理論的計算量境界が確認された。