[論文レビュー] A rigorous Multiscale Method for semi-linear elliptic problems

Patrick Henning, Axel Målqvist+1 more|arXiv (Cornell University)|Nov 27, 2012

Advanced Mathematical Modeling in Engineering参考文献 28被引用数 3

ひとこと要約

本稿では、高振動的で不均一な係数を有する半線形楕円問題を解くためのきめ細やかなスケール手法を提示する。局所的な微小スケール計算をサイズ H log(H⁻¹) のパッチで行い、一般化された有限要素基底を構築することで、粗いメッシュサイズ H に対して線形 H¹収束を達成する。また、得られるマルチスケール系を解くために減衰付きニュートン法を採用する。

ABSTRACT

In this paper we propose and analyze a new Multiscale Method for solving semi-linear elliptic problemswith heterogeneous andhighly variable coefficient functions. For this purposewe construct a generalized finiteelement basisthat spansalow dimensionalmultiscale space. The basis is assembled by performing localized linear fine-scale computations in small patches that have a diameter of order Hlog(H −1) where H is the coarse mesh size. Without any assumptions on the type of the oscillations in the coefficients, we give a rigorous proof for a linear convergence of the H 1-error with respect to the coarse mesh size. To solve the arising equations, we propose an algorithm that is based on a damped Newton scheme in the multiscale space.

研究の動機と目的

• 高変動的かつ不均一な係数を有する半線形楕円問題のための頑健なマルチスケール手法を開発すること。
• 係数の振動に関する制限的な仮定を課さずに、細かいスケールの特徴を正確に捉える低次元のマルチスケール空間を構築すること。
• 係数の正則性に関する最小限の仮定の下で、この手法の収束を厳密に解析すること。
• 得られるマルチスケール系を解くための効率的な解法アルゴリズムを、減衰付きニュートン法を用いて設計すること。

提案手法

• 粗いメッシュサイズ H に対し、直径 H log(H⁻¹) のパッチで局所的な微小スケール問題を解くことにより、一般化された有限要素基底を構築する。
• 不均一な係数の微小スケール挙動を網羅する局所解を組み合わせることで、マルチスケール空間を構成する。
• 元の問題の本質的特徴を保ちつつ、次元が削減されたマルチスケール空間におけるグローバルなマルチスケール問題を定式化する。
• ガレルキン離散化によって得られるマルチスケール空間内での非線形系を、減衰付きニュートン法を用いて解く。
• ステップ長を制御することで、反復解法の収束を保証し、グローバル収束性を維持する。
• 粗いメッシュサイズ H に対する H¹誤差の線形収束を確立するための厳密な誤差解析を実施する。

実験結果

リサーチクエスチョン

• RQ1 周期性やスケール分離の仮定を課さずに、高振動的係数を有する半線形楕円問題のためのマルチスケール手法を構築可能か？
• RQ2 マルチスケール基底構築における精度と計算効率のバランスを考慮した場合、最適なパッチサイズは何か？
• RQ3 提案手法は、粗いメッシュサイズ H に対して H¹ノルムで線形収束を達成するか？
• RQ4 得られる非線形マルチスケール系を解くための頑健かつ効率的な反復解法をどのように設計できるか？
• RQ5 係数の変動に関する制限的な仮定を課さずに、収束解析を厳密に確立できるか？

主な発見

• H¹誤差の収束は、係数の振動の種別に依存せず、粗いメッシュサイズ H に対して線形収束を達成する。
• マルチスケール基底は、サイズ H log(H⁻¹) のパッチにおける局所的計算によって構築され、計算効率とスケーラビリティを確保する。
• 係数の振動の性質（周期性やスケール分離など）に関する仮定を必要としない。
• 減衰付きニュートン法により、マルチスケール定式化から生じる非線形系に対する頑健な収束性が保証される。
• 得られるマルチスケール空間は低次元であり、正確な解の近似に必要な本質的な細かいスケールの特徴を捉えている。
• 収束解析は厳密であり、係数関数に対する仮定を最小限にとどめる条件下でも成立する。