[论文解读] Cut Finite Element Methods for Partial Differential Equations on Embedded Manifolds of Arbitrary Codimensions

Erik Burman, Peter Hansbo+2 more|arXiv (Cornell University)|Oct 5, 2016

Numerical methods in engineering被引用 4

一句话总结

本文提出了一套统一的理论框架，用于稳定剪切有限元方法求解嵌入在 R^k 中的任意余维数流形上的拉普拉斯-贝尔特拉米方程。该框架引入了三种新颖的稳定化形式，包括一种新型的表面法向梯度惩罚项，实现了与网格-几何对齐无关的最优收敛速率和鲁棒的条件数界，且在 R^3 中的曲线和曲面上通过数值实验得到验证。

ABSTRACT

We develop a theoretical framework for the analysis of stabilized cut finite element methods for the Laplace-Beltrami operator on a manifold embedded in $\\mathbb{R}^d$ of arbitrary codimension. The method is based on using continuous piecewise polynomials on a background mesh in the embedding space for approximation together with a stabilizing form that ensures that the resulting problem is stable. The discrete manifold is represented using a triangulation which does not match the background mesh and does not need to be shape-regular, which includes level set descriptions of codimension one manifolds and the non-matching embedding of independently triangulated manifolds as special cases. We identify abstract key assumptions on the stabilizing form which allow us to prove a bound on the condition number of the stiffness matrix and optimal order a priori estimates. The key assumptions are verified for three different realizations of the stabilizing form including a novel stabilization approach based on penalizing the surface normal gradient on the background mesh. Finally, we present numerical results illustrating our results for a curve and a surface embedded in $\\mathbb{R}^3$.

研究动机与目标

• 开发一个适用于嵌入在 R^k 中的任意余维数流形上稳定剪切有限元方法的一般理论框架。
• 在任意网格-几何相交情况下，建立离散问题的最优先验误差估计和条件数界。
• 验证该框架在三种不同稳定化形式中的适用性，包括一种新型的表面法向梯度稳定化。
• 将几何逼近估计和插值误差界扩展至标准最近点映射不可用的高余维情形。
• 通过 R^3 中曲线和曲面的数值实验，展示方法的鲁棒性和收敛性，涵盖不同网格对齐配置。

提出的方法

• 使用在 R^k 中的背景网格上定义的连续分片多项式有限元，限制在任意切割背景网格的离散流形上。
• 引入一种稳定化形式，确保稳定性与条件数控制，且无需匹配或形状规则的三角剖分。
• 采用三种稳定化机制：(1) 活跃网格面上法向梯度的跃迁，(2) 活跃单元上的完整梯度，以及 (3) 在背景网格上引入一种新型的表面法向梯度惩罚。
• 推导出保证条件数有界和最优收敛速率的稳定化形式的抽象条件。
• 基于仅依赖逼近和正则性假设，建立离散流形的几何逼近估计，避免显式使用最近点投影。
• 在一般余维情形下应用插值误差估计，将先前结果扩展至余维大于一的情形。

实验结果

研究问题

• RQ1 能否为任意余维数流形上的剪切有限元方法建立统一的理论框架？
• RQ2 为实现条件数界和最优误差估计而对稳定化形式提出的关键假设，是否适用于多种稳定化策略？
• RQ3 基于背景网格上表面法向梯度的新型稳定化，是否能比完整梯度稳定化实现更好的一致性和稳定性？
• RQ4 在最近点映射无法显式获得的高余维情形下，几何逼近误差如何表现？
• RQ5 条件数对网格-几何对齐的依赖关系如何？稳定化能否确保所有配置下的鲁棒性？

主要发现

• 所提出的框架对三种稳定化形式（包括新型表面法向梯度稳定化）均保证最优阶先验误差估计和有界的条件数。
• 刚度矩阵的条件数随网格尺寸 h 的变化呈 O(h^{-2}) 的量级，数值结果证实了该界在不同网格对齐配置下均成立。
• 新型表面法向梯度稳定化相比完整梯度稳定化显著降低了的一致性误差，从而支持更灵活的参数选择和更高阶逼近。
• 在 R^3 中的曲面和曲线上的数值实验表明，当施加稳定化时，条件数表现鲁棒，尤其在惩罚参数 τ 足够大时。
• 未稳定化方法在条件数随网格-几何对齐变化时表现出强烈的、不规则的峰值，而稳定化变体在所有配置下均保持良好条件性。
• 平均缩放条件数 h²κ(𝒜) 在 h→0 时保持有界并收敛于常数，证实了所有测试情形和稳定化形式下理论预测的正确性。