[论文解读] On Isogeometric Subdivision Methods for PDEs on Surfaces

Bert Jüttler, Angelos Mantzaflaris+2 more|arXiv (Cornell University)|Mar 12, 2015

Advanced Numerical Analysis Techniques被引用 2

一句话总结

该论文提出了一种基于三角形网格上Loop细分格式的等几何细分方法，用于在曲面上求解椭圆型PDE。研究证明，通过查表实现的中边积分法在高鲁棒性和高效性下实现了最优收敛率，尤其在奇异顶点附近表现优异，相较于高斯积分和重心积分法，在第二和第四阶曲面PDE中均展现出更低的计算成本和更高的稳定性。

ABSTRACT

Subdivision surfaces are proven to be a powerful tool in geometric modeling and computer graphics, due to the great flexibility they offer in capturing irregular topologies. This paper discusses the robust and efficient implementation of an isogeometric discretization approach to partial differential equations on surfaces using subdivision methodology. Elliptic equations with the Laplace-Beltrami and the surface bi-Laplacian operator as well as the associated eigenvalue problems are considered. Thereby, efficiency relies on the proper choice of a numerical quadrature scheme which preserves the expected higher order consistency. A particular emphasis is on the robustness of the approach in the vicinity of extraordinary vertices. In this paper, the focus is on Loop's subdivision scheme on triangular meshes. Based on a series of numerical experiments, different quadrature schemes are compared and a mid-edge quadrature, which is easy-to-implement via lookup tables, turns out to be a preferable choice due to its robustness and efficiency.

研究动机与目标

• 开发一种用于通过细分方法求解曲面上PDE的鲁棒且高效的等几何离散化方法。
• 研究不同数值积分方案对收敛行为和鲁棒性的影响，特别是靠近奇异顶点时的表现。
• 识别一种在计算效率、一致性和准确性之间取得平衡的积分规则，适用于第二和第四阶曲面PDE。
• 通过高效且易于实现的积分方法，实现等几何分析与现有细分建模工具的实际集成。
• 在具有大量奇异顶点的复杂几何体上验证该方法，展示其在实际工程与图形应用中的适用性。

提出的方法

• 使用Loop的细分格式生成除奇异顶点外为C2极限曲面的曲面，从而支持第二和第四阶PDE的配边有限元方法。
• 通过在细分曲面上定义基函数实现等几何分析，确保几何形状的精确表示和更高阶的一致性。
• 应用多种数值积分规则——高斯积分、自适应高斯积分、重心积分和中边积分——对Laplace-Beltrami算子和曲面双调和算子的弱形式进行积分。
• 通过预计算的查表实现中边积分，实现刚度矩阵和质量矩阵的快速且鲁棒的组装。
• 采用基于边迭代器的组装方式以优化性能，尤其针对中边积分，相比非优化方案显著降低计算成本。
• 在基准几何体（如球面、手部模型）上进行收敛性研究，通过逐步细化网格评估收敛阶数和误差范数。

实验结果

研究问题

• RQ1 不同积分规则如何影响等几何细分方法在曲面PDE上的收敛速率和鲁棒性？
• RQ2 中边积分能否在第二和第四阶曲面PDE中均实现最优收敛，尤其是在奇异顶点附近？
• RQ3 在高斯积分、自适应高斯积分、重心积分和中边积分规则之间，计算成本与精度的权衡如何？
• RQ4 该方法在具有大量奇异顶点的复杂曲面上（如具有11,586个顶点的手部模型）表现如何？
• RQ5 所提出的积分方案能否高效集成到现有的细分建模与仿真工作流中？

主要发现

• 中边积分在Laplace-Beltrami问题中实现了L2误差的h^3、H1误差的h^2和H2误差的h^1的最优收敛率，与理论预期一致。
• 在曲面双调和算子问题中，中边积分保持了L2误差的h^2和H1误差的h^1的收敛率，尽管复杂度更高，仍表现出优异鲁棒性。
• 中边积分显著快于高斯积分和自适应高斯积分，例如在Spherical-5-12网格上，Laplace-Beltrami问题的组装时间仅为0.107秒，而高斯积分需1.209秒。
• 尽管实现更简单，重心积分的效率低于中边积分，性能损失达1.5至2倍。
• 自适应高斯积分确保了最高的鲁棒性，但计算开销巨大，因此不适合实时应用。
• 该方法成功计算了复杂曲面上Laplace-Beltrami算子的前24个特征函数，证实其在特征值问题和几何数据分析中的适用性。