[论文解读] A Cut Finite Element Method with Boundary Value Correction

Erik Burman, Peter Hansbo+1 more|arXiv (Cornell University)|Jul 11, 2015

Advanced Numerical Methods in Computational Mathematics参考文献 8被引用 4

一句话总结

该论文提出了一种带有边界值校正的切割有限元方法，仅使用分段线性边界近似即可实现最优收敛性，克服了虚构域方法中的几何复杂性。通过在法线方向引入泰勒展开以校正边界条件，该方法在能量范数和$L^2$范数下均实现了最优误差估计，而无需高精度的边界表示。

ABSTRACT

In this contribution we develop a cut finite element method with boundary value correction of the type originally proposed by Bramble, Dupont, and Thomee. The cut finite element method is a fictitious domain method with Nitsche type enforcement of Dirichlet conditions together with stabilization of the elements at the boundary which is stable and enjoy optimal order approximation properties. A computational difficulty is, however, the geometric computations related to quadrature on the cut elements which must be accurate enough to achieve higher order approximation. With boundary value correction we may use only a piecewise linear approximation of the boundary, which is very convenient in a cut finite element method, and still obtain optimal order convergence. The boundary value correction is a modified Nitsche formulation involving a Taylor expansion in the normal direction compensating for the approximation of the boundary. Key to the analysis is a consistent stabilization term which enables us to prove stability of the method and a priori error estimates with explicit dependence on the meshsize and distance between the exact and approximate boundary.

研究动机与目标

• 为解决切割有限元方法（CutFEM）中在切割单元上进行精确积分的计算挑战。
• 在仅使用边界分段线性近似的情况下，实现最优收敛速率，从而简化几何计算。
• 提出一种带有边界值校正的稳定Nitsche格式，即使在边界近似误差存在的情况下，仍能保持稳定性和最优收敛性。
• 证明具有显式依赖于网格尺寸和边界近似误差的先验误差估计。
• 证明即使在低阶边界表示下，最优收敛性依然可实现，避免了标准CutFEM中所需的$O(h^{p+1})$阶边界精度要求。

提出的方法

• 通过法线方向的泰勒展开，利用近似边界上的函数值和导数值来近似精确边界上的解，实现边界值校正。
• 在近似边界上采用Nitsche型弱形式施加狄利克雷条件，并通过泰勒展开对几何误差进行校正。
• 引入一致的稳定项，控制边界附近的函数变化，确保乘子的强制性以及最优条件数$O(h^{-2})$。
• 使用最近点映射定义精确边界，实现对近似区域的一致积分。
• 采用分段线性距离函数定义离散区域，简化了几何计算。
• 推导并分析了一种结合Nitsche方法、边界校正与稳定性的对称变分格式。

实验结果

研究问题

• RQ1 是否可以在仅使用边界分段线性近似的情况下，在CutFEM中实现最优收敛性，从而避免高阶几何重构？
• RQ2 在切割有限元方法中，法线方向的泰勒展开如何提升狄利克雷边界条件的精度？
• RQ3 当边界任意切割网格时，何种稳定机制可确保稳定性和最优条件数？
• RQ4 在边界近似误差存在的情况下，所提方法是否仍能保持能量范数和$L^2$范数下的最优误差估计？
• RQ5 对于更高阶多项式近似（如$P^2$、$P^3$），该方法是否仅需少量泰勒展开项即可实现最优收敛性？

主要发现

• 即使边界仅通过分段线性表示近似，该方法在能量范数和$L^2$范数下均实现了最优收敛阶。
• 对于$P^2$和$P^3$单元，仅需泰勒展开的前两项即可实现最优收敛速率。
• 稳定项确保线性系统的条件数保持在$O(h^{-2})$，且与背景网格中边界的位置无关。
• 数值结果证实了$P^2$和$P^3$单元的最优收敛率分别为2和4，与理论预测一致。
• 即使在部分切割的单元中，该方法也无需将离散解外推至活动区域之外，仍能保持稳定性和精度。
• 误差估计显式依赖于网格尺寸$h$和精确边界与近似边界之间的距离，从而支持精确的收敛性分析。