[論文レビュー] Inverse estimates for elliptic boundary integral operators and their application to the adaptive coupling of FEM and BEM

Markus Aurada, Michael Feischl+4 more|arXiv (Cornell University)|Nov 19, 2012

Advanced Numerical Methods in Computational Mathematics参考文献 34被引用数 5

ひとこと要約

本稿は、楕円境界要素法（BEM）における4つの古典的境界積分作用素について逆推定を確立し、h-適応有限要素法-BEM結合アルゴリズムの収束解析を可能にする。主な貢献は、重み付き残差誤差推定子を用いた適応的FEM-BEMスキームの収束を証明し、特異的摂動問題において最適収束率で数値的に検証することにある。

ABSTRACT

We prove inverse-type estimates for the four classical boundary integral operators associated with the Laplace operator. These estimates are used to show convergence of an h-adaptive algorithm for the coupling of a finite element method with a boundary element method which is driven by a weighted residual error estimator.

研究の動機と目的

• ラプラス作用素に関連する4つの古典的境界積分作用素の逆型推定を確立すること。
• 特異的解を有する線形伝達問題に対するh-適応的FEM-BEM結合の収束解析を可能にすること。
• FEM-BEM結合フレームワークにデータ近似項を含めた信頼性のある後験的誤差推定子を拡張すること。
• 推定子低減技術を用いて重み付き残差誤差推定子を駆動とする適応的アルゴリズムの収束を証明すること。
• 理論的結果を数値実験により検証し、適応的グリッド・リファインメント下での最適収束率を示すこと。

提案手法

• 局所的楕円的正則性理論を用いて、単層作用素、二重層作用素、特異的作用素、結合ポテンシャル作用素の新しい逆推定を導出する。
• これらの推定を用いて、FEM-BEM結合の文脈において重み付き残差誤差推定子の効率性と信頼性を証明する。
• AFEMにおける推定子低減フレームワークをFEM-BEM設定に適応し、誤差推定子の収縮性を確立する。
• FEM-BEM結合にコスタベルの対称的結合形式を用いることで、対称性と安定性を保証する。
• マークパラメータ θ = 0.25 を用いて適応的アルゴリズムを実装し、均一リファインメントと比較する。
• 離散的境界積分作用素の構築および線形方程式系の直接解法にMATLABのHILBERT BEMライブラリを用いる。

実験結果

リサーチクエスチョン

• RQ1 一般の多面体表面において、4つの古典的境界積分作用素の逆推定を導出可能か？
• RQ2 コーナー特異性を有する問題に対して、適応的FEM-BEM結合アルゴリズムは収束するか？
• RQ3 FEM-BEM結合において、重み付き残差誤差推定子をデータ近似項を含める形に拡張可能か？
• RQ4 推定子低減を用いた適応的アルゴリズムは、最適収束率に到達可能か？
• RQ5 特異的解を有する問題において、適応的リファインメントと均一リファインメントの収束特性はどのように比較できるか？

主な発見

• 局所的楕円的正則性を用いて、多面体表面においても有効な4つの古典的境界積分作用素の逆推定を確立した。
• θ = 0.25 を用いた適応的FEM-BEM結合アルゴリズムは、ℓ → ∞ のとき誤差推定子 ϱℓ → 0 に収束する。
• 数値実験により最適収束率が確認された：FEMでは α = p/2、BEMでは β = 3/2 + q（p, q は多項式次数）。
• 均一リファインメントでは、再入角における特異性のため、多項式次数にかかわらず劣化した収束率 α = 2/7 を示す。
• 適応的リファインメントにより最適収束率が回復され、p = 1, q = 1 の場合、全体の誤差はFEM部分が支配的となり、α = 1/2 を達成する。
• 誤差推定子 ϱℓ はエネルギー誤差の信頼性のある上界を提供し、体積項および境界項の両方が最適収束率で収束する。