[論文レビュー] Asymptotics of Dirichlet eigenvalues and eigenfunctions of the Laplacian on thin domains in R^d

D. I. Borisov, Pedro Freitas|ArXiv.org|Aug 17, 2009

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

ひとこと要約

本稿では、1つの方向にスケーリングすることで (d−1) 次元の集合に近づく d 次元の細長い領域におけるラプラシアンの第1ディリクレ固有値について、3項の漸近展開を導出する。境界の幾何構造とスペクトル解析に基づく特異摂動法を用い、領域の曲率とスケーリングパラメータの関数として明示的な係数を導出する。特に、d 次元の楕円体に応用することで、細い極限に近い精度の高い固有値近似が得られる。

ABSTRACT

We consider the Laplace operator with Dirichlet boundary conditions on a domain in R^d and study the effect that performing a scaling in one direction has on the eigenvalues and corresponding eigenfunctions as a function of the scaling parameter around zero. This generalizes our previous results in two dimensions and, as in that case, allows us to obtain an approximation for Dirichlet eigenvalues for a large class of domains, under very mild assumptions. As an application, we derive a three--term asymptotic expansion for the first eigenvalue of d-dimensional ellipsoids.

研究の動機と目的

• 2次元領域におけるスケーリング下での固有値漸近挙動に関する先行研究を、任意の次元 d ≥ 2 に拡張すること。
• スケーリングパラメータ ε → 0 の極限において、ℝ^d 内の細い領域におけるラプラシアンの第1ディリクレ固有値の明示的漸近展開を導出すること。
• 展開における最初の3つの非ゼロ係数（ε^{-2}、ε^{-1}、定数項に対応）を、領域の幾何学的およびスペクトル的データを用いて計算すること。
• 一般枠組みを d 次元の楕円体に適用し、その主固有値について定量的な3項展開を導出すること。

提案手法

• d 次元領域の1方向（x_d）に対してスケーリング変換を施し、ε → 0 の極限で細い領域に還元すること。
• 境界層解析と ε のべき級数による漸近展開を用いて、固有値展開を導出すること。
• 領域を、高さ関数 h_±(x') に依存する (d−1) 次元の基本集合の摂動とモデル化すること。
• 無限大ストリップ ℝ^{d−1} × (−H₀/2, H₀/2) 上の極限作用素に対してスペクトル理論を適用し、変数分離とベッセル型固有関数を用いること。
• 基底状態 ψ₀ と曲率項 H₂、H₄ を含む内積を用いて、展開における係数 c₀、c₂、c₄ を計算すること。
• 楕円体の幾何構造を用いて、H₀、H₂、H₄ および基底状態 ψ₀ の明示的表現を計算し、係数の閉形式評価を可能にする。

実験結果

リサーチクエスチョン

• RQ1 d 次元領域が1軸に沿って (d−1) 次元の集合にスケーリングされるとき、ラプラシアンの第1ディリクレ固有値はどのように漸近的に振る舞うか？
• RQ2 このような細い領域における主固有値の漸近展開における最初の3つの非ゼロ係数の明示的形は何か？
• RQ3 一般の漸近枠組みを、d 次元の楕円体のような特定の幾何形状に適用して固有値展開を計算できるか？
• RQ4 曲率と境界の幾何構造は、高次元における固有値展開係数にどのように影響を与えるか？
• RQ5 細い極限における固有関数の構造はどのようなものか？ また、それは固有値漸近挙動とどのように関連しているか？

主な発見

• 第1ディリクレ固有値は3項の漸近展開を満たす：λ₁(ε) = π²/(4a_d²ε²) + π/(2a_d)Σ(1/a_i)ε⁻¹ + c₄ε⁰ + O(ε)，ε → 0 のとき。
• 係数 c₀^{(1,1)} = π²/(4a_d²) は、スケーリング軸上の1次元ラプラシアンに起因する ε^{-2} 項に対応する。
• 係数 c₂^{(1,1)} = π/(2a_d)Σ(1/a_i) は、基本方向における半軸の逆数の和に比例する。
• d 次元の楕円体に対しては、定数項 c₄^{(1,1)} は (x_i x_j / a_i a_j)² と基底状態 ψ₀ を含むガウス積分として明示的に計算される。
• 領域の境界に滑らかさの仮定が成り立つ限り、展開は有効である。奇数次曲率項 H₃ ≡ 0 であることにより、高次係数の計算が簡略化される。
• 本手法は、標準的摂動理論が失敗する場合でさえも、明示的な誤差項を伴う固有値近似を体系的に計算可能にする。