[论文解读] A variational proof of partial regularity for optimal transportation maps

Michaël Goldman, Félix Otto|arXiv (Cornell University)|Apr 18, 2017

Nonlinear Partial Differential Equations参考文献 12被引用 9

一句话总结

该论文提出了一种新颖的变分法证明，针对 Hölder 连续密度之间的最优传输映射的局部正则性，避免使用经典的最大值原理技术。通过在 Benamou-Brenier 的欧拉形式中构造一个竞争解，并结合基于调和逼近的 Campanato 迭代，证明了传输映射在零测度集外为 $C^{1,\beta}$，正则性通过直接在 $C^{1,\alpha}$ 层次上的 $\varepsilon$-正则性结果实现。

ABSTRACT

We provide a new proof of the known partial regularity result for the optimal transportation map (Brenier map) between two sets. Contrary to the existing regularity theory for the Monge-Amp{è}re equation, which is based on the maximum principle, our approach is purely variational. By constructing a competitor on the level of the Eulerian (Benamou-Brenier) formulation, we show that locally, the velocity is close to the gradient of a harmonic function provided the transportation cost is small. We then translate back to the Lagrangian description and perform a Campanato iteration to obtain an $ε$-regularity result.

研究动机与目标

• 提供一种全新的、纯粹基于变分法的局部正则性证明，避免依赖最大值原理。
• 直接在 $C^{1,\alpha}$ 层次上建立 $\varepsilon$-正则性结果，而非通过中间的 $C^{0,\alpha}$ 估计。
• 证明当源密度和目标密度均为 $C^{0,\alpha}$ 且远离零时，传输映射在零测度集外为 $C^{1,\alpha}$。
• 通过直接在传输映射上工作，利用位移凸性和边界层构造，统一分析框架。

提出的方法

• 在 Benamou-Brenier 欧拉形式中，通过求解具有受控通量边界条件的泊松方程，构造一个竞争解。
• 基于对偶的边界层构造导出迹估计，将问题简化为控制边界数据。
• 证明当传输能量较小时，速度场在更小球体内接近某个调和函数的梯度，通过调和逼近实现。
• 通过调和函数的二阶泰勒展开，实现“通过倾斜改善平坦度”，从而实现迭代精化。
• 利用仿射变换下的不变性进行 Campanato 迭代，实现能量和 $C^{1,\alpha}$ 正则性的几何衰减。
• 将欧拉近似结果回传至拉格朗日框架，得出传输映射在 Lebesgue 点邻域内为 $C^{1,\alpha}$ 的结论。

实验结果

研究问题

• RQ1 能否通过完全基于变分法的方法建立最优传输映射的局部正则性，而无需依赖最大值原理或 Alexandrov 解理论？
• RQ2 在能量较小且具有 Hölder 连续性假设下，最优传输映射是否在 $C^{1,\alpha}$ 层次上存在 $\varepsilon$-正则性结果？
• RQ3 能否在 Benamou-Brenier 形式中利用调和逼近，推导出传输映射的迭代正则性改进？
• RQ4 位移凸性与边界层构造之间的相互作用，如何促成欧拉框架中有效竞争解的构造？
• RQ5 Campanato 迭代在多大程度上可被调整以直接获得 $C^{1,\alpha}$ 正则性，从而跳过中间的 $C^{0,\alpha}$ 步骤？

主要发现

• 传输映射 $T$ 在全测度开集 $E' \subseteq \text{supp}\, \rho_0$ 及其像集 $F' \subseteq \text{supp}\, \rho_1$ 上为 $C^{1,\alpha}$，从而确立了部分 $C^{1,\alpha}$ 正则性。
• 直接在 $C^{1,\alpha}$ 层次上证明了 $\varepsilon$-正则性定理：若在某球体内能量与 Hölder 半范数足够小，则 $T$ 在更小球体内为 $C^{1,\alpha}$。
• 关键步骤是调和逼近结果：在欧拉形式中，当能量较小时，速度场接近某个调和函数的梯度。
• 证明依赖于通过泊松方程与边界层构造构建的变分竞争解，其中迹估计 (3.10) 起到核心作用。
• 由于问题的自然标度性，Campanato 迭代直接作用于 $C^{1,\alpha}$ 层次，避免了中间的 $C^{0,\alpha}$ 估计。
• 结果在仿射变换下保持不变，这对平坦度的迭代改进至关重要。