[論文レビュー] An inverse iteration method for obtaining q-eigenpairs of the p-Laplacian in a general bounded domain

Grey Ercole|arXiv (Cornell University)|Oct 14, 2015

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

ひとこと要約

本稿では、p > 1 かつ 1 ≤ q < p である一般の有界領域 Ω ⊂ ℝ^N における p-Laplacian の q-固有対を計算するための逆反復法を提案する。固定点反復とレイリー商型正規化を組み合わせることで、最小の q-固有対への収束が効率的に達成され、一般の領域に対しても収束性と安定性が保証される。

ABSTRACT

Let $Ω$ be a bounded and smooth domain of $\mathbb{R}^{N}$, $N\geq2$, and consider the eigenvalue problem: $-Δ_{p}u=λ\left| u ight| _{L^{q}(Ω)}^{p-q}\left| u ight| ^{q-2}u$ in $Ω,$ $u=0$ on $\partialΩ,$ where $p&gt;1$, $1\leq q

研究の動機と目的

• 一般の有界で滑らかな領域 Ω ⊂ ℝ^N における p-Laplacian の q-固有対を数値的に計算するための手法を開発すること。
• 固有値が固有関数の L^q ノルムに依存するため、標準的なスペクトル法が複雑になる状況に対処すること。
• 領域の幾何構造やパrameter 範囲（p > 1, 1 ≤ q < p）が一般である場合でも、反復スキームの収束性と安定性を保証すること。
• 有限要素法の離散化に適した計算効率が高く、安定したアルゴリズムを提供すること。

提案手法

• 固有関数の L^q ノルムを含む非線形固有値問題として p-Laplacian 固有値問題を定式化する。
• 各ステップでレイリー商型正規化を施した線形化 p-Laplacian 問題を解く逆反復スキームを導入する。
• 固定点反復を用いて、固有関数と固有値の推定値を同時に更新し、L^q ノルムによる正規化を保証する。
• 空間的離散化にガラーキン有限要素法を適用し、一般の領域における数値計算を可能にする。
• 固有値問題の q-斉次構造と整合性を保つ正規化戦略を採用する。
• 領域およびパrameter に関する適切な仮定の下で、反復列が最小の q-固有対に収束することを証明する。

実験結果

リサーチクエスチョン

• RQ1 一般の有界領域における p-Laplacian の最小 q-固有対を計算するための逆反復法を設計できるか？
• RQ2 固有値が固有関数の L^q ノルムに依存する場合、反復固有値ソルバー内でどのように扱えるか？
• RQ3 1 ≤ q < p および p > 1 の場合に、提案された逆反復スキームの収束を保証する条件は何か？
• RQ4 非線形固有値問題の標準的手法と比較して、数値的にどのように性能を発揮するか？
• RQ5 対称性の仮定なしに、一般の滑らかな領域に対しても本手法は頑健で効率的か？

主な発見

• 提案された逆反復法は、任意の有界で滑らかな領域 Ω ⊂ ℝ^N において、p-Laplacian の最小 q-固有対にグローバルに収束する。
• 固有関数の L^q ノルムによる正規化のおかげで、q が 1 に近いか、p が大きい場合でも安定性と精度が維持される。
• 標準的な領域およびパrameter に関する仮定の下で収束が保証され、収束速度はスペクトルギャップに依存する。
• 数値実験により、本手法の頑健性と効率性が、さまざまな領域形状やパrameter 選択において確認された。
• 本アルゴリズムは有限要素法の離散化に適しており、複雑な幾何形状における実用的計算を可能にする。
• 固有値項に含まれる L^q ノルムによる非線形性を効果的に扱い、標準的なパワー法で生じる問題を回避する。