[論文レビュー] Two-scale extensions for non-periodic coefficients

Vsevolod Laptev|ArXiv.org|Dec 6, 2005

Advanced Mathematical Modeling in Engineering参考文献 13被引用数 6

ひとこと要約

本稿は、楕円型PDEにおける非周期的で急激に振動する係数に対して、局所的周期的近似を用いた均質化を可能にする二尺度拡張フレームワークを導入する。$ a(x,y) $ を $ a(x, x/\bar{u}) = a_M(x) $ となるように構築することにより、二尺度収束を保証し、$ a_M(x) $ が周期的でない場合でも、有効な均質化問題の厳密な導出を可能にする。

ABSTRACT

We consider non-homogeneous media with properties which can be characterized by rapidly oscillated coefficients. For such coefficients we define a notion of two-scale extension, present several ways to construct two-scale extensions, discuss their properties and relation to homogenization

研究の動機と目的

• 実材料に一般的に見られるが、古典的周期的均質化には適用できない非周期的で急激に振動する係数を有するPDEの均質化の課題に対処すること。
• 与えられた非周期的係数 $ a_M(x) $ に対して、特定のスケール $ \bar{\varepsilon} $ で一致する二尺度拡張 $ a(x,y) $ を体系的に構築する手法を開発すること。
• 構築された二尺度拡張が二尺度収束の文脈における適切なテスト関数として適格であることを保証し、解の収束証明を可能にすること。
• このような拡張を用いて導かれた均質化極限問題が、周期的ケースと同一の構造的性質および収束性を保持することを示すこと。
• 特に緩やかに変化する領域や既知の統計的均質性を有する領域において、効率的なセル問題の解法を実現する実用的な計算戦略を提供すること。

提案手法

• $ a(x,y) $ を $ y $ に関して $ Y $-周期的とし、$ a(x, x/\bar{\varepsilon}) = a_M(x) $ となるように二尺度拡張を定義する。ここで $ \bar{\varepsilon} $ は微細構造に対応する固定スケールである。
• 連続的拡張（$ \Omega $ 上での補間を用いたもの）と離散的拡張（$ \Omega $ を部分領域に分割し、局所的セル問題を解くもの）の2種類の拡張を構築する。
• 二尺度収束理論を用いて、係数 $ a^\varepsilon(x) = a(x, x/\varepsilon) $ を有する問題列 $ \mathcal{P}^\varepsilon $ の解 $ u_\varepsilon $ の極限を分析する。
• 二尺度収束におけるテスト関数としての適格性を確認することで、二尺度拡張の適格性を確立し、$ u_\varepsilon \to u_0 $ の収束を保証する。
• 有効係数 $ A(x) $ を周期セル $ Y $ 上での $ a(x,y) $ の平均として導出し、連続的拡張では $ A(x) $ が連続的であり、離散的拡張では領域 $ \Omega_k $ の中心で独立にセル問題を解くことで、$ A(x) $ が区分定数となるようにする。
• 2階楕円型方程式にこのフレームワークを適用し、適切な条件下で $ H^1 $ 収束が成り立つこと、かつ均質化解 $ u_0 $ が元の解 $ u $ を近似できることを示す。

実験結果

リサーチクエスチョン

• RQ1 非周期的係数 $ a_M(x) $ が、ある $ Y $-周期的関数 $ a(x,y) $ を用いて $ a(x, x/\bar{\varepsilon}) $ の形に表現可能であるか。これにより、二尺度法を用いた均質化が可能になるか。
• RQ2 構築された二尺度拡張 $ a(x,y) $ が、二尺度収束の意味で適格であるための条件は何か。
• RQ3 連続的拡張と離散的拡張の選択が、有効係数 $ A(x) $ の正則性および計算コストにどのように影響するか。
• RQ4 非周期的ケースにおいて、均質化解 $ u_0 $ が元の解 $ u $ をどの程度近似できるか。収束が成り立つ条件は何か。
• RQ5 $ a_M(x) $ の緩やかな変化領域、周期的領域、または既知の平均値を有する領域を活用することで、計算コストの削減が達成可能か。

主な発見

• 任意の可測かつ有界な係数 $ a_M(x) $ に対して、$ \bar{\varepsilon} $ を $ a(x, x/\bar{\varepsilon}) = a_M(x) $ を満たすように選べる限り、二尺度拡張 $ a(x,y) $ が存在し、非周期的設定下での均質化を可能にする。
• 連続的拡張により、平均化係数 $ A(x) $ が連続的となる。これは、セル問題の積分における振動項の小ささを用いて証明される。
• 離散的拡張により、$ A(x) $ は区分定数となる。$ A(x) $ は、部分領域 $ \Omega_k $ の中心で独立に $ N $ 個のセル問題を解くことで計算され、計算コストを顕著に削減する。
• 小さな $ \bar{\varepsilon} $ に対しては、有効係数 $ A(x) $ は $ a_M(x) $ の振動的性質を継承するが、$ \bar{\varepsilon} $ が大きくなるにつれて滑らかになり、空間的に一様な領域では定数に近づく可能性がある。
• 拡張の適格性が保証される限り、$ a(x,y) $ に滑らかさの仮定がなくても、$ u_\varepsilon $ の二尺度収束 $ u_\varepsilon \to u_0 $ が成立する。これにより、周期的ケースを超えた均質化結果の有効性が拡張される。
• 1次元および2次元における数値的検証により、均質化解 $ u_0 $ と微細スケール解 $ u $ の間に良好な一致が得られ、$ H^1 $-補正近似により精度が向上することが示された。