[论文解读] Approximation of length minimization problems among compact connected sets

Matthieu Bonnivard, Antoine Lemenant+1 more|arXiv (Cornell University)|Mar 9, 2014

Computational Geometry and Mesh Generation参考文献 38被引用 44

一句话总结

本文提出了一种用于二维紧致连通集上长度最小化问题的新型Ambrosio-Tortorelli型逼近方法，通过引入加权测地距离项以在极限情况下强制实现连通性。该方法通过Γ-收敛性成功逼近了Steiner问题、平均距离问题以及p-柔度能量问题，数值验证表明相场的零水平集能涌现出连通的极小化子。

ABSTRACT

In this paper we provide an approximation \\`a la Ambrosio-Tortorelli of some classical minimization problems involving the length of an unknown one-dimensional set, with an additional connectedness constraint, in dimension two. We introduce a term of new type relying on a weighted geodesic distance that forces the minimizers to be connected at the limit. We apply this approach to approximate the so-called Steiner Problem, but also the average distance problem, and finally a problem relying on the p-compliance energy. The proof of convergence of the approximating functional, which is stated in terms of Gamma-convergence relies on technical tools from geometric measure theory, as for instance a uniform lower bound for a sort of average directional Minkowski content of a family of compact connected sets.

研究动机与目标

• 解决二维连通性约束下长度最小化问题缺乏数值逼近方法的问题。
• 开发一种相场方法，通过新颖的加权测地距离项在极限情况下强制实现连通性。
• 将Ambrosio-Tortorelli框架扩展至Steiner问题、平均距离问题及p-柔度能量最小化等新问题。
• 利用几何测度论工具，证明所提泛函在ε → 0时关于原始几何问题的Γ-收敛性。
• 提供一种可数值实现的算法，通过迭代优化与次梯度推进法生成连通极小化子。

提出的方法

• 提出一种新泛函，结合标准Modica-Mortola项与加权测地距离项：$ \frac{1}{4\varepsilon}\int_{\Omega}(1-\varphi)^2dx + \varepsilon\int_{\Omega}|\nabla\varphi|^2dx + \frac{1}{c_\varepsilon}\sum_{i=1}^N d_\varphi(x_i,x_1) $，其中$ d_\varphi $为由$ \varphi $加权的测地距离。
• 将加权测地距离定义为$ d_\varphi(x,y) = \inf\left\{ \int_\gamma \varphi \, d\mathcal{H}^1 \,;\, \gamma \text{ 连接 } x \text{ 与 } y \right\} $，可通过快速推进法实现数值计算。
• 利用性质：若$ \sum_{i=1}^N d_\varphi(x_i,x_1) = 0 $，则$ \{\varphi = 0\} $为路径连通且包含所有$ x_i $，从而在极限下实现连通性约束。
• 应用Γ-收敛性理论，证明当$ \varepsilon \to 0 $时，逼近泛函收敛于原始几何问题，依赖于方向性Minkowski上密度的统一下界。
• 实现一种数值算法，利用次梯度推进法同时计算$ d_\varphi $及其上梯度，使凹项的优化成为可能。
• 采用递减$ \varepsilon $的迭代优化策略，每一步均以前一解为初始值，以避免局部极小值并促进收敛至连通极小化子。

实验结果

研究问题

• RQ1 能否构建一种相场方法，以逼近紧致连通集约束下的长度最小化问题？
• RQ2 如何在无显式拓扑约束的扩散界面模型中编码连通性约束？
• RQ3 加权测地距离项的引入是否能确保相场零水平集在$ \varepsilon \to 0 $极限下保持连通？
• RQ4 该框架能否扩展至其他几何问题，如平均距离问题与p-柔度能量最小化？
• RQ5 所提数值方法在不同问题配置与初始条件下是否具有鲁棒性，能稳定生成连通极小化子？

主要发现

• 所提泛函$ G^{h}_{\Lambda,\varepsilon} $结合Modica-Mortola项与加权测地距离项，在$ \varepsilon \to 0 $时Γ-收敛于原始几何问题，确保极小化子的收敛性。
• 数值模拟中$ \varepsilon = 0.05 $时，$ \phi $的零水平集保持连通并包含指定点$ y_0 $，即使$ y_0 $偏离中心位置亦成立。
• 随着长度惩罚参数$ \Lambda $增大，连通集$ \{\phi = 0\} $的长度减小，与预期物理行为一致。
• 采用递减$ \varepsilon $的迭代优化策略，成功避免局部极小值，并在多次运行中生成一致且连通的极小化子。
• 次梯度推进法可高效计算$ d_\varphi $及其上梯度，使凹项的优化成为可行。
• 该方法成功逼近了三种不同几何问题：Steiner问题、平均距离问题及p-柔度能量最小化问题，均在连通性约束下完成。