Dynamic Optimal Transport with Optimal Preferential Paths

技术摘要：具有最优优选路径的动态最优传输

问题表述
本研究探讨了一个定义在域 $\Omega \subset \mathbb{R}^d$ 上的动态最优传输问题，该域由一个体区域（bulk region）和一个低维、不相交且规则的曲线 $\Gamma \subset \Omega$ 组成。该模型描述了质量的传输过程，其中人群可以在体域 $\Omega$ 内移动，也可以沿曲线 $\Gamma$ 移动，并可能在两个区域之间以一定的代价进行质量交换。系统的总质量是守恒的。

动力学过程由一对耦合的连续性方程控制。令 $\rho_t$ 和 $\mu_t$ 分别表示体域和曲线上随时间变化的质量分布，并伴随相应的动量 $J_t$ 和 $V_t$。质量交换通过通量项 $f_t$（法向交换）以及边界项 $G^0_t, G^1_t$（曲线端点的切向交换）进行建模。系统受非线性迁移率函数 $m_\Omega$ 和 $m_\Gamma$ 的约束，这些函数分别调节体域和曲线中的动力学过程。

目标是最小化一个作用泛函 $A_\Gamma$，该泛函聚合了体域中的动能、曲线上的动能以及两者之间质量转移的代价。泛函定义为：
$$A_\Gamma = \int_\Omega \frac{|J_t|^2}{m_\Omega(\rho_t)} d\lambda_\Omega + \alpha_1 \int_\Gamma \frac{|V_t|^2}{m_\Gamma(\mu_t)} d\lambda_\Gamma + \alpha_2 \int_\Gamma \frac{|f_t|^2}{m_\Gamma(\mu_t)} d\lambda_\Gamma + \dots$$
其中 $\alpha_1, \alpha_2 > 0$ 是控制曲线传输相对于体域传输以及进入或离开曲线代价的参数。

本文讨论了两种不同的情景：

1. 固定曲线： $\Gamma$ 的几何形状是固定的，目标是寻找最优的演化测度 $(\rho_t, \mu_t, J_t, V_t, f_t)$。
2. 变化曲线： 曲线的形状也是优化问题中的变量。目标是同时寻找最优曲线和最优传输动力学。

方法论
作者采用变分法中的直接法来建立极小值点的存在性。

• 固定曲线的存在性： 作者定义了一个满足耦合连续性方程（在弱分布意义下）的可容许元组空间。他们证明了作用泛函对于测度的窄收敛（narrow convergence）是凸的且下半连续的。通过建立随时间变化的测度的 Hölder 估计，并利用广义 Arzelà-Ascoli 定理，他们证明了极小化序列的紧致性，从而证明了极小值点的存在性。
• 渐近分析： 本文分析了系统在质量交换代价参数 $\alpha_2$ 和 $\alpha_3$ 趋于零或无穷大时的行为。
  • 当 $\alpha \to \infty$ 时，交换代价变得难以承受，问题解耦为 $\Omega$ 和 $\Gamma$ 上两个独立的最佳传输问题（前提是满足质量兼容性条件）。
  • 当 $\alpha \to 0$ 时，交换变得免费，问题收敛为一个具有零交换代价的耦合传输问题。
• 变化曲线的存在性： 为了处理对曲线形状的优化，作者在作用泛函中引入了一个正则化项，以防止自交并确保正则性。他们利用了最初为闭曲线定义的切点能量（Tangent-Point energy, $E_{tp}^p$），并将其适配于此处处理开曲线的情况。该能量惩罚自交行为，并确保具有有限能量的曲线是嵌入的一维子流形。总正则化泛函 $R(\gamma)$ 包括切点能量、一个防止曲线闭合的对数项（以确保边界条件明确）以及长度惩罚。
  • 通过证明切点能量限制了曲线的 $W^{2-1/p, p}$ Sobolev-Slobodeckij 范数，进而实现紧致性，该范数紧嵌入于 $C^{1,\alpha}$ 中。
  • 通过结合测度的紧致性（命题 4.6）与作用泛函及正则化泛函的下半连续性，证明了联合问题（曲线与测度）极小值点的存在性。

数值方案
作者开发了一种基于原对偶方法和增广拉格朗日算法（ALG2）的数值算法。

• 离散化： 使用有限元方法对时空域进行离散化，其中曲线采用分片线性单元，时空域采用棱柱单元。
• 算法： 对于固定曲线，算法通过交替更新对偶变量（势函数）和原变量（密度与动量）来求解鞍点问题。
• 变化曲线： 为了优化曲线，增广拉格朗日法被嵌入到投影梯度下降循环中。针对曲线控制点的目标泛函梯度通过中心差分进行近似。切点能量使用特定的求积规则进行离散化，以处理奇异性。

主要贡献与结果

1. 严谨的表述： 本文为体域与低维曲线之间的质量传输提供了耦合连续性方程的严谨定义，包括质量交换机制和非线性迁移率。
2. 极小值点的存在性： 建立了对于固定及变化优选路径的动态传输问题极小值点的存在性。
3. 渐近极限： 研究了交换代价消失或趋于无穷大时耦合系统的极限行为，将该模型与具有迁移率的经典非耦合最优传输联系起来。
4. 通过切点能量进行正则化： 本工作将切点能量的应用扩展到最优传输框架下的开曲线，证明了该能量能够强制执行必要的正则性（非自交及 $C^1$ 正则性），从而保证了变化曲线设定下解的存在性。
5. 数值实现： 作者展示了一种功能性的数值方案，能够计算固定曲线及优化曲线下的测地线，并演示了成本参数和几何形状对传输动力学的影响。

意义
本文声称通过将 Benamou-Brenier 框架推广到包含域内低维结构（曲线），为动态最优传输理论做出了贡献。这模拟了在特定路径（如高速公路）上传输更快或更便宜的情景，但访问这些路径需要付出代价。通过允许路径本身被优化，这项工作解决了最优规划中的问题，即基础设施（路径）是根据人群的初始和最终分布来确定的。使用非线性迁移率以及通过切点能量对曲线几何形状进行严谨处理，使这项工作区别于以往侧重于凸对偶或仅关注固定边界的研究。数值结果说明了传输代价与路径几何形状之间的相互作用如何决定质量的演化，为建模复杂的传输网络提供了工具。