[論文レビュー] Convergence rate for quasilinear eigenvalue homogenization

Fern Andez Bonder, Juan Pablo Pinasco+1 more|arXiv (Cornell University)|Nov 1, 2012

Advanced Mathematical Modeling in Engineering参考文献 19被引用数 4

ひとこと要約

本稿は、周期的ミクロ構造を有する擬線形楕円型作用素の変分固有値について、G収束およびPr"ufer変換技術を用いて明示的な収束速度を確立する。固有値はパrameter $ \varepsilon $ に対して $ O(\varepsilon) $ の速度で収束することを証明し、その速度は固有値インデックス $ k $ に依存する。線形ケースの結果を非線形 $ p $-ラプラシアン枠組みへ拡張する。

ABSTRACT

In this work we study the homogenization problem for nonlinear eigenvalues of quasilinear elliptic operators. We obtain an explicit order of convergence in $k$ and in $\varepsilon$ for the (variational) eigenvalues.

研究の動機と目的

• 振動係数を有する擬線形楕円型作用素の変分固有値の収束速度を特定すること。
• 既知の線形固有値収束結果を $ 1 < p < \infty $ の非線形 $ p $-ラプラシアン設定へ拡張すること。
• 収束速度が固有値インデックス $ k $ およびスケールパrameter $ \varepsilon $ にどのように依存するかを定量化すること。
• 非線形で周期的な複合材料における固有値ホモジゼーションのための厳密な解析的枠組みを提供すること。

提案手法

• 単調作用素に対するG収束理論を用いて、擬線形固有値問題のホモジゼーション極限を正当化する。
• p-ラプラシアン固有値問題を第一階微分方程式系に変換することで、数値的および解析的取り扱いを可能にするため、Pr"ufer変換を適用する。
• 誤差許容範囲を制御できるように、二分法およびニュートン・ラプソン修正を用いたシューティング法を用いて、数値的に固有値を計算する。
• 周期的重みおよび係数構造下でのPr"ufer位相関数および振幅関数の漸近的解析を通じて、誤差バインドを導出する。
• Weyl型漸近公式 $ \lambda_k \sim k^p $ を用いて、固有値インデックス $ k $ を固有値の大きさに関連づけ、収束推定値を精緻化する。
• 周期的重み $ \rho_\varepsilon(x) = \rho(x/\varepsilon) $ のための固有値および固有関数の数値シミュレーションを通じて、理論的収束速度を検証する。

実験結果

リサーチクエスチョン

• RQ1 $ \varepsilon \to 0 $ のとき、擬線形固有値ホモジゼーション問題における $ k $ 番目の変分固有値の明示的収束次数は何か？
• RQ2 非線形 $ p $-ラプラシアン設定において、収束速度は固有値インデックス $ k $ およびスケールパrameter $ \varepsilon $ にどのように依存するか？
• RQ3 周期的ミクロ構造において、スペクトル全域にわたって一様に収束速度を定量化・境界づけられるか？
• RQ4 固有値列 $ \lambda_k^\varepsilon $ の振動は $ \varepsilon \to 0 $ のときにも持続するか？その挙動はどのように説明できるか？

主な発見

• 本稿は、擬線形固有値ホモジゼーション問題における $ k $ 番目の変分固有値について、収束速度 $ O(\varepsilon) $ を確立する。
• 収束速度は固有値インデックス $ k $ に明示的に依存し、$ k $ が大きい場合に $ Ck^p \varepsilon $ のスケーリングに従う。これはWeylの漸近法則と整合的である。
• 数値結果により、$ \lambda_k^\varepsilon $ が $ \varepsilon \to 0 $ のとき極限固有値 $ \lambda_k $ に収束し、その際の振動は小さくなることが確認された。
• 固有関数 $ u_k^\varepsilon $ は $ \varepsilon $ が小さくなるにつれて $ u_k $ とのノード交差回数が増加するが、差の振幅は減少する。
• Pr"ufer変換により、二分法およびニュートン・ラプソン修正を用いた誤差制御による固有値の高精度な数値計算が可能になった。
• 固有値における周期的振動は、点ごとの不等式では完全に説明できないため、積分重みを用いたSturm型比較定理の必要性が示唆された。