[论文解读] Minimal geodesics along volume preserving maps, through semi-discrete optimal transport

Quentin Mérigot, Jean‐Marie Mirebeau|arXiv (Cornell University)|May 13, 2015

Geometric Analysis and Curvature Flows参考文献 2被引用 4

一句话总结

本文提出一种数值方法，利用半离散最优传输计算体积保持映射群中的最小测地线，实现了不可压缩流体欧拉方程非经典、多值解的首次数值重建。该方法利用布里内尔的广义极分解，在广义流的松弛设定下收敛至解，揭示了类似分形的支撑结构和超越经典流场范围的非确定性粒子轨迹。

ABSTRACT

We introduce a numerical method for extracting minimal geodesics along the group of volume preserving maps, equipped with the L2 metric, which as observed by Arnold solve Euler's equations of inviscid incompressible fluids. The method relies on the generalized polar decomposition of Brenier, numerically implemented through semi-discrete optimal transport. It is robust enough to extract non-classical, multi-valued solutions of Euler's equations, for which the flow dimension is higher than the domain dimension, a striking and unavoidable consequence of this model. Our convergence results encompass this generalized model, and our numerical experiments illustrate it for the first time in two space dimensions.

研究动机与目标

• 开发一种稳健的数值方法，用于在 $L^2$ 度量下计算体积保持微分同胚群中的最小测地线，其对应于不可压缩流体欧拉方程的解。
• 通过在广义流的松弛设定下工作，解决在 $d \geq 3$ 维中经典极小化器失效的问题，其中粒子轨迹可能发生分裂和交叉。
• 通过半离散最优传输数值实现布里内尔的广义极分解，即使在经典解不存在时也能提取最小测地线。
• 提供二维空间中非经典、多值流体流的首次数值图示，其支撑结构具有类似分形的特性。

提出的方法

• 该方法将测地线问题表述为在连续路径空间 $\Omega = C^0([0,1], X)$ 上的概率测度上的凸优化问题，通过最小化作用量 $\int_0^1 |\dot{\omega}(t)|^2 dt$，并满足不可压缩性和边界约束条件。
• 采用布里内尔的广义极分解，将体积保持映射分解为梯度映射与测度保持映射的复合，从而通过半离散最优传输实现数值计算。
• 半离散最优传输框架对目标测度（最终流体构型）进行离散化，同时保持源测度（初始构型）为连续形式，从而高效计算传输计划。
• 通过求解涉及压力作为拉格朗日乘子的对偶问题来计算广义流，其在松弛的广义流设定下建立了收敛性。
• 通过从最优传输计划中采样来重构粒子轨迹，并通过聚类和盒维数估计方法可视化非确定性行为。
• 该方法已实现于开源软件中，网址为 https://github.com/mrgt/EulerSemidiscrete，支持可复现性与进一步实验。

实验结果

研究问题

• RQ1 当高维中由于振荡行为导致经典解不存在时，能否在体积保持映射群中数值计算最小测地线？
• RQ2 在二维空间中，非经典、多值流体流（粒子分裂且轨迹交叉）的行为如何？能否通过数值方法捕捉此类行为？
• RQ3 此类广义流的支撑结构具有何种几何特性？其是否如理论分析所建议的那样表现出类似分形的性质？
• RQ4 压力场的正则性如何与重构流体轨迹的非确定性特征相关联？
• RQ5 当内部区域出现全局非确定性时，数值解在边界附近在多大程度上仍保持确定性？

主要发现

• 当 $t_{\text{max}} < 1$ 时，计算得到的解恢复了单位正方形中的经典贝尔特拉米流，具有确定性的粒子轨迹和光滑的压力场。
• 当 $t_{\text{max}} \geq 1$ 时，解变为非经典且非确定性，粒子轨迹发生交叉和分裂，表明广义流行为的出现。
• 解的支撑结构的盒维数从 $D=2$（当 $t_{\text{max}}=0.9$ 时）增加至 $t_{\text{max}} \in \{1.1, 1.3, 1.5\}$ 时严格介于 2 和 3 之间的值，表明其具有分形结构。
• 在区域边界附近，粒子运动仍近似保持确定性且呈顺时针方向，而内部则表现出高度非确定性的逆时针混合行为。
• 压力梯度场在初始和最终时刻附近显示出正则性降低的迹象，与理论结果 $L^2_{\text{loc}}(]0,T[, \text{BV}(X))$ 一致。
• 该算法成功重构了广义流，涵盖圆盘反转和贝尔特拉米流两种情形，后者在 $t_{\text{max}} = 1.5$ 时表现出更强的非确定性行为。