[論文レビュー] Maximization of the Steklov eigenvalues with a diameter constraint

Abdelkader Al Sayed, Beniamin Bogosel+2 more|arXiv (Cornell University)|Apr 29, 2020

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

ひとこと要約

この論文は、直径制約下でのステクロフ固有値の最大化を検討し、体積または周囲制約とは異なり、球は1番目の非自明固有値に対しても最適化の最大化者でないことを証明している。形状最適化技術と数値アルゴリズムを用いて、最適な凸領域の存在を確立し、高次の固有値に対しては、最適形状がより平らになり、最適解で複数の固有値が現れることを示している。数値結果により、$ k = 1,2,3 $ に対して非凸の最大化者であることが確認された。

ABSTRACT

In this paper, we address the problem of maximizing the Steklov eigenvalues with a diameter constraint. We provide an estimate of the Steklov eigenvalues for a convex domain in terms of its diameter and volume and we show the existence of an optimal convex domain. We establish that balls are never maximizers, even for the first non-trivial eigenvalue that contrasts with the case of volume or perimeter constraints. Under an additional regularity assumption, we are able to prove that the Steklov eigenvalue is multiple for the optimal domain. We illustrate our theoretical results by giving some optimal domains in the plane thanks to a numerical algorithm.

研究の動機と目的

• 固定された直径制約 $D(\Omega) = d_0$ の下でステクロフ固有値 $\sigma_k(\Omega)$ の最大化という形状最適化問題に取り組む。
• 直径制約下でのステクロフ固有値最大化問題に対して、最適な凸領域の存在を確立する。
• 体積または周囲制約下での既知の最大性とは対照的に、球が直径制約下でも $\sigma_k(\Omega)$ の最大化者であるかどうかを検証する。
• 特に固有値の重複度と、最適解で直径制約が正確に2点で達成されることを含め、最適領域の構造を分析する。
• 凸性を仮定しない数値アルゴリズムを開発・検証し、非凸の最大化者を探索可能にする。

提案手法

• ステクロフ固有値の正homogeneity（$\sigma_k(t\Omega) = \sigma_k(\Omega)/t$）を用いて、直径制約問題を $D(\Omega)\sigma_k(\Omega)$ を最大化する等価な非制約問題に再定式化する。
• 形状微分理論を用いて、ステクロフ固有値および直径関数の両方の勾配情報を導出し、数値最適化に用いる。
• 固定区間 $[-D/2, D/2]$ 上の2つの関数 $f_1, f_2$ を用いた平面領域のパラメータ化を行い、離散的な凸性・凹性制約を課して凸性を保証する。
• FreeFEMを用いて有限要素法により、離散化された領域上でのステクロフ固有値問題を数値的に解く。
• $\sigma_k(\Omega)D(\Omega)$ を最大化する勾配型最適化アルゴリズムを実装し、$f_1$ および $f_2$ の摂動が境界の法線変化に対応するようにする。
• 2番目の数値的手法では凸性制約を緩和し、非凸の最大化者を探索。その結果、最適形状が必ずしも凸であるとは限らないことが確認された。

実験結果

リサーチクエスチョン

• RQ1 任意の次元において、直径制約下で1番目の非自明固有値 $\sigma_1(\Omega)$ の最大化者としての球が成立するか？
• RQ2 直径制約下での $\sigma_k(\Omega)$ の最適領域において、最適解で複数の固有値が現れるか？
• RQ3 最適領域において直径制約が正確に2点で達成されるか？この性質は異なる固有値においても保持されるか？
• RQ4 $k$ が増加するにつれて $\sigma_k(\Omega)$ の最適形状はどのように変化するか？より平らになるか？
• RQ5 非凸領域が凸領域よりも高い $\sigma_k(\Omega)D(\Omega)$ 値を達成できるか？これは凸性仮定なしでの最大化者存在の意味にどのような示唆を与えるか？

主な発見

• 任意の $k \geq 1$ に対して、球は直径制約下での $\sigma_k(\Omega)$ の最大化者でない。これは、体積または周囲制約下での最大性とは対照的である。
• 1番目の固有値に対しては、最適領域 $\omega_1^*$ が $\sigma_1(\omega_1^*)D(\omega_1^*) = 2.13623$ を満たし、数値的形状は角付近でわずかな非凸性を示している。
• $k=2$ の場合、最適領域は $\sigma_2(\omega_2^*)D(\omega_2^*) = 4.92925$ を達成し、境界の一部で明確な非凸性の兆候が見られる。
• $k=3$ の場合、最適領域は $\sigma_3(\omega_3^*)D(\omega_3^*) = 7.76108$ を達成し、凸の場合とは顕著に異なる形状を示している。
• 各 $k$ に対して最適領域は、正則性仮定の下で多重ステクロフ固有値を持つことが証明された。
• すべての数値的解において、直径制約が正確に2点で達成されており、理論的期待と整合的である。