[論文レビュー] Shape optimization of a weighted two-phase Dirichlet eigenvalue

Idriss Mazari, Grégoire Nadin+1 more|arXiv (Cornell University)|Jan 9, 2020

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

ひとこと要約

本稿は、ドリフト項を含む重み付き2相ディリクレ固有値作用素の形状最適化を検討し、領域が球でない限り、正則な最小化子は存在しないことを証明する。球の場合、新規のホモジエナイゼーションに基づく手法と2階形状微分の単調性原理を用いて、中心に配置された球対称なバン・バイ・バン分布（中心資源）が最適であることを確立する。この手法により安定性解析が大幅に簡素化され、球対称性のもとで中心配置が最小化子であることが確認される。

ABSTRACT

Let $m$ be a bounded function and $\alpha$ a nonnegative parameter. This article is concerned with the first eigenvalue $\lambda\_\alpha(m)$ of the drifted Laplacian type operator $\mathcal L\_m$ given by $\mathcal L\_m(u)= -\operatorname{div} \left((1+\alpha m) abla u ight)-mu$ on a smooth bounded domain, with Dirichlet boundary conditions. Assuming uniform pointwise and integral bounds on $m$, we investigate the issue of minimizing $\lambda\_\alpha(m)$ with respect to $m$. Such a problem is related to the so-called "two phase extremal eigenvalue problem" and arises naturally, for instance in population dynamics where it is related to the survival ability of a species in a domain. We prove that unless the domain is a ball, this problem has no "regular" solution. We then provide a careful analysis in the case of a ball by: (1) characterizing the solution among all radially symmetric resources distributions, with the help of a new method involving a homogenized version of the problem; (2) proving in a more general setting, a stability result for the centered distribution of resources with the help of a monotonicity principle for second order shape derivatives which significantly simplifies the analysis.

研究の動機と目的

• ドリフト項を含む重み付き2相ディリクレ作用素の第一固有値の最小化子の存在と構造を調査すること。
• スペクトル最適化問題における最適リソース分布が正則かバン・バイ・バンかを特定すること、特に径対称性下での場合を対象とすること。
• 球領域における中心の径対称リソース分布が最小化子としての安定性を確立すること。
• ホモジエナイゼーションと2階形状微分解析を組み合わせた新規手法を構築・適用し、安定性証明を簡素化すること。
• 解析を高次元に拡張し、強制性結果の一般化可能性を議論すること。

提案手法

• 非球形領域における正則最小化子の欠如を分析するため、Murat-TartarおよびCox-Liptonの手法を適応する。
• 球内における最適径対称分布を特徴付けるために、固有値問題のホモジエナイズド版を導入する。
• 2階形状微分の単調性原理を適用し、複雑な固有関数展開を必要とせずに中心分布の安定性を証明する。
• 体積制約下でのラグランジアンの第一および第二階形状微分を計算し、最適性条件を分析する。
• 高次元における球面調和関数展開を用いて2次元を超える一般化を図り、第二微分形式の対角化を実施する。
• 陰関数定理を用いて、L²およびW¹,²ノルムにおける密度パラメータに関する固有対のC∞-滑らか性を確立する。

実験結果

リサーチクエスチョン

• RQ1 重み付き2相ディリクレ作用素の第一固有値が正則最小化子をもつのはどのような条件下か？
• RQ2 球領域において、径対称でバン・バイ・バンのリソース分布が最適か？
• RQ3 明示的な固有関数計算を伴わずに、中心配置の安定性を証明できるか？
• RQ4 径対称性下での第二階形状微分の振る舞いはどのようにか？最適性の証明に利用可能か？
• RQ5 2次元の手法を3次元を含む高次元に拡張可能か？

主な発見

• 領域が球でない限り、問題は正則最小化子をもたない。これは、最適分布が一般にバン・バイ・バンであることを示唆する。
• 球領域の場合、中心に配置された径対称バン・バイ・バン分布 m∗₀ は、すべての径対称リソース分布の中で最適である。
• 2階形状微分の単調性原理を用いた安定性証明により、従来の手法と比較して解析が大幅に簡素化される。
• 第二階形状微分に対する強制性推定が確立され、二次形式に一様な下界 C > 0 が存在し、局所的最適性が保証される。
• この手法は高次元（例：3次元）にまで拡張可能であり、球面調和関数展開により第二微分の対角化が可能となり、最大値原理の議論により強制性の下界が保たれる。
• 固有対 (uα,m, λα(m)) が密度パラメータ m に関してL²の意味でC∞-滑らかであることが示され、最適化フレームワークの解析的安定性が保証される。