[論文レビュー] Localized bases for finite dimensional homogenization approximations with non-separated scales and high-contrast

Houman Owhadi, Lei Zhang|arXiv (Cornell University)|Nov 3, 2010

Advanced Mathematical Modeling in Engineering参考文献 35被引用数 25

ひとこと要約

本稿では、分離スケールやエルゴード性を仮定せず、非分離スケールおよび高対比係数を有する楕円型・放物型・放物型PDEの解に対する有限次元近似空間を構築する局所的でメッシュに依存しない手法を提案する。ソース項が $ L^2 $ に属する場合、解空間が $ H^1 $ にコンパクトに埋め込まれることを活用し、局所PDEを解くことで局所的基底関数を構築する。これにより、基底関数がサブドメイン上にサポートされ、サイズ $ O(h^{\beta} \text{ln}(1/h)) $ の領域で、$ O(h^{-d}) $ 次元の空間において $ \tilde{O}(h^{2-2\beta}) $ の精度が達成される。$ \beta \ne 1 $ である。バッファゾーンを含めることで、高対比媒体においても精度が保たれる。

ABSTRACT

We construct finite-dimensional approximations of solution spaces of divergence form operators with $L^\infty$-coefficients. Our method does not rely on concepts of ergodicity or scale-separation, but on the property that the solution space of these operators is compactly embedded in $H^1$ if source terms are in the unit ball of $L^2$ instead of the unit ball of $H^{-1}$. Approximation spaces are generated by solving elliptic PDEs on localized sub-domains with source terms corresponding to approximation bases for $H^2$. The $H^1$-error estimates show that $\mathcal{O}(h^{-d})$-dimensional spaces with basis elements localized to sub-domains of diameter $\mathcal{O}(h^α\ln \frac{1}{h})$ (with $α\in [1/2,1)$) result in an $\mathcal{O}(h^{2-2α})$ accuracy for elliptic, parabolic and hyperbolic problems. For high-contrast media, the accuracy of the method is preserved provided that localized sub-domains contain buffer zones of width $\mathcal{O}(h^α\ln \frac{1}{h})$ where the contrast of the medium remains bounded. The proposed method can naturally be generalized to vectorial equations (such as elasto-dynamics).

研究の動機と目的

• スケール分離やエルゴード性に依存しない $ L^∞ $ 係数を有する発散型PDEの有限次元近似法の開発。
• 古典的均質化法が失敗する高対比およびマルチスケール媒体における解空間の近似に挑戦すること。
• 正確な $ H^1 $ ノルム近似を達成する、局所的で計算効率の良い基底関数の構築。
• 同じ局所的基底フレームワークを用いて、放物型および双曲型問題へと拡張すること。
• 係数行列の最小固有値に依存しない、厳密な $ H^1 $-誤差推定を確立すること。

提案手法

• ソース項 $ g ∈ L^2(Ω) $ の場合、解空間が $ H^1_0(Ω) $ にコンパクトに埋め込まれることを活用し、$ H^{-1}(Ω) $ ではなく、有限次元近似を可能にする。
• サブドメイン上で $ g $-に似たソース項を有する局所的楕円型PDEを解くことで、近似基底を構築し、局所的サポートを保証する。
• エネルギーノルムと同値であり、係数 $ a $ に依存しないフラックスノルム $ \|\cdot\|_{a\text{-flux}} $ を用いて誤差推定を導出する。
• [10]のトランスファープロパティを適用し、局所解のサポート外での減衰を制御することで、指数的減衰を伴う局所化を可能にする。
• 高対比媒体において、精度を保つために、幅 $ O(h^\beta \ln(1/h)) $ のバッファゾーンを導入する。
• 同じ局所的基底構築法を用いて、ひずみ動力学などのベクトル型方程式へと一般化する。

実験結果

リサーチクエスチョン

• RQ1 スケール分離やエルゴード性を仮定せず、マルチスケールPDEのための有限次元近似空間を構築できるか？
• RQ2 高対比媒体においても、メッシュサイズ $ h $ に従って $ H^1 $-誤差が減少するような局所的基底関数をどのように構築できるか？
• RQ3 与えられた $ H^1 $-ノルムの精度を達成するための最適な局所的サブドメインのサイズは何か？
• RQ4 同じ局所的基底フレームワークを、放物型および双曲型方程式などの時間依存問題へと拡張できるか？
• RQ5 係数 $ a $ に高対比が存在する場合、局所的近似の精度にどのように影響し、それをどのように緩和できるか？

主な発見

• ソース項が $ L^2(Ω) $ に属する場合、PDEの解空間は $ H^1_0(Ω) $ にコンパクトに埋め込まれ、$ H^1 $-ノルムにおいて任意の精度で有限次元近似が可能である。
• $ O(h^{-d}) $ 次元の近似空間で、基底関数が直径 $ O(h^\beta \ln(1/h)) $ のサブドメインに局所化され、$ \beta \in [1/2, 1) $ である場合、$ H^1 $-ノルムにおいて $ O(h^{2-2\beta}) $ の精度が達成される。
• 高対比媒体では、局所的サブドメインに、対比が有界である幅 $ O(h^\beta \ln(1/h)) $ のバッファゾーンを含めることで、$ O(h^{2-2\beta}) $ の精度を維持できる。
• フラックスノルム誤差推定は $ \lambda_{\min}(a) $ および $ \lambda_{\max}(a) $ に依存しないため、広範な係数範囲においてロバストである。
• 本手法は、弾性動力学などのベクトル型方程式へ自然に一般化され、同じ収束率を保つ。
• ソース項 $ g \in H^{-\nu} $、$ \nu < 1 $ の下で解空間の強いコンパクト性が保たれることにより、古典的均質化、数値的均質化、低次元モデル化の統一的基盤が提供される。