[論文レビュー] Numerical investigation of the smallest eigenvalues of the p-Laplace operator on planar domains

Jiří Horák|arXiv (Cornell University)|Jun 3, 2011

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

ひとこと要約

本稿は、$p \in [1.1, 10]$ の範囲で、平面領域における $p$-ラプラシアン作用素の2番目に小さい固有値を、Sobolev空間 $W^{1,p}_0(\Omega)$ 内の制約付き勾配降下法およびマウンテンパス法を用いて数値的に調査している。その結果、非回転対称領域（長方形や三角形など）において、第二固有関数の対称性が $p$ に応じて変化し、領域に依存する臨界 $p_0$ で遷移することが明らかになった。一方、ドミンで第二固有関数は、すべての $p$ に対してノードラインに関する奇対称性を保っている。

ABSTRACT

The eigenvalue problem for the p-Laplace operator with p>1 on planar domains with the zero Dirichlet boundary condition is considered. The Constrained Descent Method and the Constrained Mountain Pass Algorithm are used in the Sobolev space setting to numerically investigate the dependence of the two smallest eigenvalues on p. Computations are conducted for values of p between 1.1 and 10. Symmetry properties of the second eigenfunction are also examined numerically. While for the disk an odd symmetry about the nodal line dividing the disk in halves is maintained for all the considered values of p, for rectangles and triangles symmetry changes as p varies. Based on the numerical evidence the change of symmetry in this case occurs at a certain value p_0 which depends on the domain.

研究の動機と目的

• ヒルバート空間からバナッハ空間 $W^{1,p}_0(\Omega)$ への制約付きマウンテンパス法の拡張（$p \neq 2$ の場合）。
• 広範囲（$1.1 \leq p \leq 10$）における $p$ による2番目に小さい固有値の依存関係の計算および解析。
• $p$ の変化に伴う第二固有関数の対称性の変化、特に非回転対称領域における挙動の調査。
• $p$ が 2 から大きく離れた場合の数値的手法の妥当性の検証およびメッシュの細分化による収束性・安定性の評価。
• $p \to 1^+$ および $p \to \infty$ の極限における既知の漸近的挙動と数値結果の比較。

提案手法

• Banach空間 $W^{1,p}_0(\Omega)$（$p \in (1, \infty)$）における制約付き勾配降下法および制約付きマウンテンパス法の適応。
• 制約 $\|u\|_{L^p} = 1$ の下で最小化する汎関数として、レイリー商 $\lambda = \frac{\int_\Omega |\nabla u|^p \, dx}{\int_\Omega |u|^p \, dx}$ を用いる。
• 降下方向 $w_u$ を、$-\Delta_p w = \Delta_p u + \alpha |u|^{p-2}u$ の解として定義し、$w_u \in T_u S$ を満たすように正規化する。
• $(-\Delta_p)^{-1}$ を用いて $|u|^{p-2}u$ を $W^{1,p}_0(\Omega)$ 内の関数に写像し、反復的更新を可能にする。
• 適応的メッシュ細分化を用いた有限要素法により問題を離散化し、収束性を確保する。
• 初期条件 $u(0) \in S$ を満たす時間刻み幅スキーム $\frac{d}{dt}u(t) = w_{u(t)}$ を適用し、臨界点への収束を追跡する。

実験結果

リサーチクエスチョン

• RQ1 平面領域において、$p$ が 1.1 から 10 に変化する際、最小固有値 $\lambda_1$ はどのように変化するか？
• RQ2 第二固有関数 $u_2$ は、さまざまな $p$ 値において、対称性（例えばノードラインに関する奇対称性）を保っているか？
• RQ3 長方形や三角形などの非回転対称領域において、第二固有関数の対称性が変化する臨界 $p_0$ はどの値か？
• RQ4 $p$ が 2 から大きく離れた場合、数値的手法の収束性および安定性はいかがなったか？
• RQ5 計算された固有値および固有関数は、$p \to 1^+$ および $p \to \infty$ の極限において、既知の漸近的挙動とどの程度一致しているか？

主な発見

• 最小固有値 $\lambda_1$ は、$p$ が 1.1 から 10 に増加するに従い、理論的予想と一致して単調に減少する。
• ドミンでは、すべての $p \in [1.1, 10]$ において、第二固有関数 $u_2$ が円を二等分する直径に関する奇対称性を保っている。
• 長方形や三角形では、第二固有関数の対称性が領域の幾何的性質に依存する臨界 $p_0$ で変化しており、固有構造における分岐を示している。
• 制約付きマウンテンパス法は、$p$ が 2 から大きく離れた場合でも第二固有ペアに収束することが確認され、バナッハ空間設定下での頑健性が示された。
• $p \to 1^+$ および $p \to \infty$ の極限における数値結果は、既知の漸近的挙動（$p \to 1^+$ 時に $\lambda_1 \to \infty$、$p \to \infty$ 時に $\lambda_1 \to $ 有限の極限値）と定性的に一致している。
• 逆 $p$-ラプラシアンを用いた降下方向の定義により、臨界点への収束が保証され、収束の証明には $W^{1,p}_0(\Omega)$ 内のコンパクト性およびノルム推定が用いられた。