[論文レビュー] Localization of Elliptic Multiscale Problems

Axel Målqvist, Daniel Peterseim|arXiv (Cornell University)|Oct 4, 2011

Advanced Mathematical Modeling in Engineering被引用数 11

ひとこと要約

本稿は、頂点パッチ上の局所的解から基底関数を構築することで、顕著に異種な係数を有する楕円型問題に対する局所化された一般化された有限要素法を提案する。基底関数の指数的減衰を証明し、$ H \log(1/H) $ のパッチサイズで最適収束を達成する。また、スケール分離や解の正則性を仮定しない、$ \|A^{1/2}\nabla(u - u_H^{\text{ms}})\|_{L^2} \leq C_g H $ のロバストな誤差境界を確立する。

ABSTRACT

This note constructs a local generalized finite element basis for elliptic problems with heterogeneous and highly varying coefficients. The basis functions are solutions of local problems on vertex patches. The error of the corresponding generalized finite element method decays exponentially with respect to the number of element layers in the patches. Hence, on a uniform mesh of size H, patches of diameter H\log(1/H) are sufficient to preserve a linear rate of convergence in H without any pre-asymptotic or resonance effects. The analysis does not rely on regularity of the solution or scale separation in the coefficient. This result motivates new and justifies old classes of variational multiscale methods.

研究の動機と目的

• 顕著に変動する非周期的係数を有する楕円型問題において、細粒度の効果を正確に捉える一般化された有限要素法の開発を目的とする。
• 古典的有限要素法がスケール分解能の欠如によりこのようなマルチスケール設定で失敗することの是正を目的とする。
• 多項式でない、係数に依存する関数からなる基底を頂点パッチに局在化して構築することを目的とする。
• 結果として得られる手法が、局所的計算のみを用いてエネルギーノルムにおいて最適収束率 $ \mathcal{O}(H) $ を達成することを証明することを目的とする。
• 変分的マルチスケール法およびマルチスケール有限要素法に依存する局所化された基底関数に基づく理論的基盤を提供することを目的とする。

提案手法

• 係数依存性を組み込むために、準補間作用素を用いて修正されたノード基底を構築する。
• 局所的基底関数を、同次ディリクレ境界条件を満たす頂点パッチ上の局所的PDEの解として定義する。
• 関連するノードから離れるに従い、局所的基底関数が指数的に減衰することを証明し、有限パッチへの切断が精度を損なわないと正当化する。
• $ H $ に関して線形収束を維持するために、粗メッシュ要素の $ k \approx \log(1/H) $ 層のパッチサイズを用いる。
• これらの局所化された基底関数を組み合わせて適合ガラーキン近似を構築することで、一般化された有限要素法を確立する。
• エネルギーノルムにおける誤差を分析し、局所問題の離散化に関して安定性を証明する。

実験結果

リサーチクエスチョン

• RQ1 パッチ上の局所的解のみを用いて、マルチスケール楕円型問題に対して最適収束を達成する一般化された有限要素法を構築できるか？
• RQ2 頂点パッチ上の局所的PDEの解は指数的減衰を示すか？ これにより、精度を損なわずに有限パッチに切り詰め可能か？
• RQ3 スケール分離や解の正則性の仮定なしに、エネルギーノルムにおいて $ \mathcal{O}(H) $ の収束を達成できるか？
• RQ4 局所問題の離散化に対して誤差境界は安定か？
• RQ5 本分析は、変分的マルチスケール法およびマルチスケール有限要素法における対称的かつ局所化された基底関数の使用を正当化するか？

主な発見

• 一般化された有限要素法の誤差はメッシュサイズ $ H $ に対して線形に減少し、$ \|A^{1/2}\nabla(u - u_H^{\text{ms}})\|_{L^2} \leq C_g H $ を満たす。係数の変動に依存しない。
• 局所的基底関数の指数的減衰が証明され、精度を損なわずにパッチを $ \mathcal{O}(H \log(1/H)) $ の直径に切り詰め可能である。
• スケール分離、周期性、解の正則性を仮定しないで最適収束を達成する。
• パッチサイズ $ k \approx \log(1/H) $ の粗メッシュ要素の層は、線形収束率を維持するのに十分である。
• 本分析は、変分的マルチスケール法における対称的かつ局所化された基底関数の使用を正当化し、適応的パッチの細分化の理論的基盤を提供する。
• 誤差境界は局所問題の離散化に対して安定であり、実装におけるロバスト性を確認する。