[论文解读] 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

一句话总结

本文研究在直径约束下最大化斯泰克洛夫特征值的问题，证明了球体永远不是最大值点——即使对于第一个非平凡特征值也是如此，这与体积或周长约束下的情况相反。通过形状优化技术和数值算法，本文建立了最优凸区域的存在性，并表明更高特征值的最优形状趋于扁平化，且在最优解处出现多重特征值。数值结果证实，对于 $ 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)$ 的最大值点，这与它在体积或周长约束下已知的最大性形成对比。
• 分析最优区域的结构，特别是特征值的重数以及直径约束在恰好两点处达到饱和的性质。
• 开发并验证一种无需假设凸性的数值算法，以计算最优区域，从而探索非凸最大值点。

提出的方法

• 利用斯泰克洛夫特征值的正齐次性（$\sigma_k(t\Omega) = \sigma_k(\Omega)/t$），将直径约束问题转化为等价的无约束问题，涉及 $D(\Omega)\sigma_k(\Omega)$。
• 应用形状导数理论，对斯泰克洛夫特征值和直径泛函分别推导梯度信息，用于数值优化。
• 采用两个函数 $f_1, f_2$ 在固定区间 $[-D/2, D/2]$ 上对平面区域进行参数化，并施加离散的凸性/凹性约束以保证凸性。
• 在 FreeFEM 中使用有限元方法，对离散化区域上的斯泰克洛夫特征值问题进行数值求解。
• 实施一种梯度型优化算法，以最大化 $\sigma_k(\Omega)D(\Omega)$，其中 $f_1$ 和 $f_2$ 的扰动对应于法向边界变化。
• 在第二种数值方法中放松凸性约束，以探索非凸最大值点，结果证实最优形状不一定是凸的。

实验结果

研究问题

• RQ1 在任意维度下，球体是否是直径约束下第一个非平凡斯泰克洛夫特征值 $\sigma_1(\Omega)$ 的最大值点？
• RQ2 在直径约束下，$\sigma_k(\Omega)$ 的最优区域在最优解处是否会出现多重特征值？
• RQ3 直径约束是否总能在最优区域的恰好两点处达到饱和？这一性质是否在不同特征值间保持一致？
• RQ4 随着 $k$ 的增加，$\sigma_k(\Omega)$ 的最优形状如何变化？它们是否会变得越来越扁平？
• RQ5 非凸区域是否可能使 $\sigma_k(\Omega)D(\Omega)$ 的值超过凸区域？这在无凸性假设下对最大值点的存在性意味着什么？

主要发现

• 对于任意 $k \geq 1$，球体在直径约束下永远不是 $\sigma_k(\Omega)$ 的最大值点，这与它在体积或周长约束下的已知最大性相矛盾。
• 对于第一个特征值，最优区域满足 $\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$ 的最优区域均具有多重斯泰克洛夫特征值。
• 所有数值解中，直径约束均在恰好两点处达到饱和，与理论预期一致。