[论文解读] Discretization of the 3D Monge-Ampere operator, between Wide Stencils and Power Diagrams

Jean‐Marie Mirebeau|arXiv (Cornell University)|Mar 3, 2015

Geometry and complex manifolds参考文献 1被引用 9

一句话总结

该论文通过结合宽模板有限差分与幂图的几何洞察，提出了一种在结构化网格上对三维Monge-Ampère算子进行单调且一致的离散化方法，实现了高精度与鲁棒的收敛性。该方法保证了阻尼牛顿求解器的全局收敛，并在处理各向异性和奇异解方面，尤其在三维情况下，优于现有方案。

ABSTRACT

We introduce a monotone (degenerate elliptic) discretization of the Monge-Ampere operator, on domains discretized on cartesian grids. The scheme is consistent provided the solution hessian condition number is uniformly bounded. Our approach enjoys the simplicity of the Wide Stencil method, but significantly improves its accuracy using ideas from discretizations of optimal transport based on power diagrams. We establish the global convergence of a damped Newton solver for the discrete system of equations. Numerical experiments, in three dimensions, illustrate the scheme efficiency.

研究动机与目标

• 开发一种在结构化网格上对3D Monge-Ampère算子进行单调且一致的离散化方法，克服现有非单调或低精度方案的局限性。
• 通过结合最优传输与幂图的几何原理，提升宽模板方法的精度，同时保持局部、基于网格的实现方式。
• 确保对离散系统 $\mathcal{D}_V u = \rho$ 的阻尼牛顿求解器实现全局收敛，从而鲁棒地求解具有狄利克雷边界条件的Monge-Ampère PDE。
• 在具有挑战性的测试用例（包括奇异和高度各向异性的解）中展示该方法的有效性，其中标准方案会失效或收敛缓慢。
• 证明在Hessian矩阵条件数有界时，一致性得以保持，从而在3D中可靠逼近粘性解。

提出的方法

• 该方案使用用户定义的 $V$ 模板（$\mathbb{Z}^3$ 中的向量集合），包含对称与带符号的排列，通过可变步长的二阶有限差分构造离散Hessian近似。
• 提出一种对称化、一致的Hessian二次型 $\langle e, (\nabla^2 u) e \rangle$ 近似，采用非均匀模板以适应域边界。
• 离散Monge-Ampère算子 $\mathcal{D}_V u$ 定义为由这些二阶差分构造的对称正定矩阵的行列式，确保退化椭圆性（单调性）。
• 该方法利用幂图的几何洞察提升角分辨率与一致性，同时将相互作用局部化至邻近网格点，避免复杂的全局几何谓词。
• 在网格上对离散系统 $\mathcal{D}_V u = \rho$ 应用阻尼牛顿求解器，其全局收敛性在所提出的单调性与一致性条件下得到保证。
• 该方案在维度上具有独立性，但特别针对3D域设计并验证，通过半定规划验证了在Hessian矩阵条件数有界时的一致性。

实验结果

研究问题

• RQ1 能否构造一种单调且一致的3D Monge-Ampère算子离散化方法，避免依赖全局几何结构（如幂图）的同时保持精度？
• RQ2 模板大小与形状的选择如何影响各向异性或退化Hessian矩阵下离散格式的一致性与收敛性？
• RQ3 在奇异或高度各向异性的解上，该方法在精度与收敛速率方面相较于宽模板与有限差分方案的优越程度如何？
• RQ4 Hessian条件数与所提方案一致性的关系是什么？有界条件数是否足以保证收敛至粘性解？
• RQ5 能否通过局部、基于网格的方法实现与幂图等全局几何方法相当的精度，而无需任意精度算术或符号扰动？

主要发现

• 所提出的方案 ${\cal D}_V$ 在小模板下对所有满足 $\operatorname{Tr}(M)/\det(M)^{1/3} \leq 7.8$ 的正定Hessian矩阵 $M$ 实现一致性，大模板下可扩展至 $11.9$。
• 在二次测试用例中，非单调有限差分方案 ${\cal D}^{\text{FD}}$ 在高分辨率下无法收敛，而所提方案保持分辨率无关的一致性。
• 在平滑圆锥测试用例中，尽管计算成本更高，所提方案使用大模板时的 $L^\infty$ 误差显著低于宽模板方案，尤其在高分辨率下表现更优。
• 在奇异测试用例中，所提方案在 $50^3$ 网格上实现的 $L^\infty$ 误差约为宽模板方案的1/20，且小模板与大模板的误差相同，表明对模板大小具有鲁棒性。
• 在 $50^3$ 网格上的计算时间显示，所提方案耗时4.7–9.7分钟（小至大模板），优于宽模板方案（6.4–65分钟），并与有限差分方案（6.5分钟）相当。
• 阻尼牛顿求解器对所提方案实现了全局收敛，且该方法成功恢复了Hessian条件数有界的弱Alexandrov解，而基于幂图的方法则需全局几何计算，无法实现。