[论文解读] Convergent discrete Laplace-Beltrami operators over surfaces

Jyh-Yang Wu, Mei-Hsiu Chi+1 more|arXiv (Cornell University)|Apr 20, 2010

Advanced Numerical Analysis Techniques参考文献 16被引用 11

一句话总结

本文提出了一种新颖且收敛的算法，用于在三角化曲面上计算离散的拉普拉斯-贝尔特拉米算子，该方法利用局部切多边形和配置方程，将顶点位置与函数值提升至切平面。该方法实现了 O(r) 的收敛速率，其中 r 为网格尺寸，是首个在一般曲面上被证明收敛的此类算法，在精度与速度方面优于现有方法。

ABSTRACT

The convergence problem of the Laplace-Beltrami operators plays an essential role in the convergence analysis of the numerical simulations of some important geometric partial differential equations which involve the operator. In this note we present a new effective and convergent algorithm to compute discrete Laplace-Beltrami operators acting on functions over surfaces. We prove a convergence theorem for our discretization. To our knowledge, this is the first convergent algorithm of discrete Laplace-Beltrami operators over surfaces for functions on general surfaces. Our algorithm is conceptually simple and easy to compute. Indeed, the convergence rate of our new algorithm of discrete Laplace-Beltrami operators over surfaces is $O(r)$ where r represents the size of the mesh of discretization of the surface.

研究动机与目标

• 解决数值几何偏微分方程中，针对一般曲面缺乏可证明收敛的离散拉普拉斯-贝尔特拉米算子的问题。
• 开发一种概念简单且计算高效的算法，确保收敛至连续的拉普拉斯-贝尔特拉米算子。
• 在一般曲面离散化下，建立该离散算子收敛性的数学证明。
• 在收敛速率与计算效率方面，优于现有方法如 Xu 的与 Taubin 的方法。

提出的方法

• 引入局部切多边形的概念，通过顶点处的切平面近似该顶点周围的曲面邻域。
• 通过正交投影将离散网格上的顶点位置与函数值提升至局部切多边形，以保持几何与函数关系的完整性。
• 在二维切多边形上求解配置方程，计算权重以使用类似有限差分的方案近似拉普拉斯算子。
• 通过邻近顶点的加权和定义离散拉普拉斯-贝尔特拉米算子，权重系数由配置方程的解得出。
• 利用泰勒展开与几何提升方法，将切平面上的离散算子与曲面上的连续拉普拉斯算子关联起来。
• 应用格林公式与离散积分，推导出最终的离散算子表达式，确保与连续情形的一致性。

实验结果

研究问题

• RQ1 能否构建一个在一般三角化曲面上收敛于连续算子的离散拉普拉斯-贝尔特拉米算子？
• RQ2 此类离散算子的收敛速率是多少，能否进行数学证明？
• RQ3 与 Xu 的方法和 Taubin 的方法相比，该方法在精度与效率方面表现如何？
• RQ4 通过提升至切平面与配置方程方法，能否实现稳定且一致的离散化？
• RQ5 通过将泰勒展开扩展至二阶以上，是否可能实现更高阶的收敛？

主要发现

• 所提出的离散拉普拉斯-贝尔特拉米算子以 O(r) 的速率收敛至连续算子，其中 r 为网格尺寸，已在主定理中得到证明。
• 数值结果表明，该方法在最大误差方面显著优于 Xu 的方法，尤其在高阶函数下表现更优。
• 对于函数 F1，所提方法的误差为 2.36×10⁻⁵×h⁴，而 Xu 的方法为 2.34×10⁻³×h²，表明收敛速率快得多。
• 该方法比 Xu 的方法更快：在 h₀/2¹⁰ 条件下，域 (c) 的计算时间仅为 0.010–0.015 秒，而 Xu 的方法为 0.024–0.026 秒。
• 该方法在不同类型的网格（包括结构化与非结构化三角化）下均保持 O(r) 的收敛速率。
• 理论分析证实，该提升与配置方法保持了几何精度，从而实现了 O(r) 的收敛性。