2025-11-13T00:13:15.031611

Ghost stabilisation for cut finite element exterior calculus

Di Pietro, Droniou, Nilsson

We introduce the cut finite element method in the language of finite element exterior calculus, by formulating a stabilisation -- for any form degree -- that makes the method robust with respect to the position of the interface relative to the mesh. We prove that the $L^2$-norm on the physical domain augmented with this stabilisation is uniformly equivalent to the $L^2$-norm on the ``active'' mesh that contains all the degrees of freedom of the finite element space (including those external to the physical domain). We show how this CutFEEC method can be applied to discretize the Hodge Laplace equations on an unfitted mesh, in any dimension and any topology. A numerical illustration is provided involving a conforming finite element space of $H^{\text{curl}}$ posed on a filled torus, with convergence and condition number scaling independent of the position of the boundary with respect to the background mesh.

academic

Ghost stabilisation for cut finite element exterior calculus

基本信息

论文ID: 2510.14772
标题: Ghost stabilisation for cut finite element exterior calculus
作者: Daniele A. Di Pietro (IMAG, CNRS, Montpellier), Jérôme Droniou (IMAG, CNRS & Monash University), Erik Nilsson (IMAG, CNRS, Montpellier)
分类: math.NA cs.NA
发表时间: October 17, 2025
论文链接: https://arxiv.org/abs/2510.14772

摘要

本文在有限元外微积分(FEEC)的框架下引入了切割有限元方法，通过为任意形式度数制定稳定化项，使方法对界面相对于网格的位置具有鲁棒性。作者证明了物理域上增强稳定化的 $L^2$ 范数与包含所有有限元空间自由度的"活跃"网格上的 $L^2$ 范数一致等价。该CutFEEC方法可应用于在不拟合网格上离散化任意维数和拓扑的Hodge Laplace方程。数值验证在填充环面上使用 $H^{\text{curl}}$ 一致有限元空间，收敛性和条件数的尺度独立于边界相对于背景网格的位置。

研究背景与动机

问题背景

传统有限元方法的局限性: 传统有限元方法要求网格与域的几何边界一致，这对复杂几何形状或演化界面的网格生成带来挑战。
切割有限元方法的需求: CutFEM允许界面切割背景网格的单元，简化了网格生成和自适应过程，但需要适当的稳定化技术确保稳定性和精度。
混合问题的挑战: 当几何以特殊方式切割网格时（如 $|T \cap \Omega| \ll |T|$ ），会观察到相关线性系统条件数的显著恶化。

研究动机

理论统一: 将CutFEM技术与有限元外微积分(FEEC)框架相结合，为微分形式的数值方法提供统一理论。
稳定性保证: 开发适用于任意形式度数的稳定化项，确保方法对界面位置的鲁棒性。
拓扑保持: 在离散层面保持de Rham复形的同调结构，这对数值格式的稳定性至关重要。

核心贡献

提出了CutFEEC框架: 将切割有限元方法扩展到有限元外微积分，适用于任意维数和形式度数。
设计了通用稳定化项: 为任意 $k$ -形式开发了ghost penalty稳定化项，这是对现有混合ghost penalty项的推广。
建立了范数等价性: 证明了物理域上增强稳定化的 $L^2$ 范数与活跃网格上 $L^2$ 范数的一致等价性。
应用于Hodge Laplace方程: 展示了该方法在不拟合网格上离散化Hodge Laplace方程的有效性。
数值验证: 在填充环面上验证了方法的收敛性和条件数独立性。

方法详解

任务定义

本文研究在不拟合网格上求解Hodge Laplace方程的混合变分形式：

\begin{cases} (d^{k-1}\delta^k + \delta^{k+1}d^k)\eta = f - \pi f & \text{in } \Omega \\ \pi\eta = 0 & \text{in } \Omega \end{cases}

其中 $d^k$ 是外导数， $\delta^k$ 是余微分， $\pi$ 是调和形式空间上的正交投影。

核心技术框架

1. 不拟合网格设置

活跃网格: $T_h = \{T \in T_{0,h} : T \cap \Omega \neq \emptyset\}$
活跃域: $\Omega_h = (\bigcup_{T \in T_h} T)^\circ$
切割单元: $T_h^{\text{cut}} = \{T \in T_h : T \not\subset \Omega\}$
稳定化面: $F_h^\partial = \{F \in \bigcup_{T \in T_h^{\text{cut}}} F_T : F \not\subset \partial\Omega_h\}$

2. 微分形式的切向和法向部分

对于 $k$ -形式 $\omega \in \Lambda^k T$ 和面 $F \subset \partial T$ ：

切向部分: $\omega_\parallel = (\pi_F^* \circ \gamma)\omega$ 法向部分: $\omega_\perp = \omega|_F - \omega_\parallel$

其中 $\gamma = \iota^*$ 是标准迹算子， $\gamma_n$ 是法向迹算子。

3. Ghost Penalty稳定化项

对于 $\omega \in V_h^{k,r}$ ，定义跳跃： $[\omega] = \omega_1|_F - \omega_2|_F$

关键性质: $[\omega] = [\omega_\perp]$ （只有法向部分参与跳跃）

稳定化项定义为： $s(\omega,\zeta) = \sum_{F \in F_h^\partial} \sum_{\ell=0}^r \eta h_F^{2\ell+1} \int_F \left(\gamma_n[\nabla_n^{(\ell)}\omega] \wedge \star \gamma_n[\nabla_n^{(\ell)}\zeta] + \gamma[\nabla_n^{(\ell)}\omega] \wedge \star \gamma[\nabla_n^{(\ell)}\zeta]\right)$

其中 $\nabla_n^{(\ell)}$ 是 $\ell$ 阶法向方向导数， $\eta > 0$ 是罚参数。

4. Ghost内积和范数

$(\omega,\zeta)_s = (\omega,\zeta)_\Omega + s(\omega,\zeta)$ $\|\omega\|_s = \sqrt{(\omega,\omega)_s}$

理论分析

主要定理（范数等价性）

定理1: 范数 $\|\cdot\|_{\Omega_h}$ 和 $\|\cdot\|_s$ 在 $V_h^{k,r}$ 上一致等价。

证明思路:

局部控制（引理3）: 通过Taylor展开和跳跃项控制切割单元上的范数
全局控制（引理4）: 利用cut-to-uncut路径将控制扩展到整个活跃网格
逆向估计: 使用标准逆不等式和迹不等式

几何: 填充环面 $\Omega = \{(x,y,z) \in \mathbb{R}^3 : [(\sqrt{x^2+y^2}-0.5)^2 + z^2]^{1/2} \leq 0.25\}$
有限元空间:
- $\sigma_h$ : $P_1^-\Lambda^0\Omega_h$ （Lagrange元）
- $\eta_h$ : $P_1^-\Lambda^1\Omega_h$ （第一类Nédélec元）
网格尺寸: $h = \frac{1}{13}, \frac{1}{26}, \frac{1}{52}$
稳定化参数: $\eta = 1$ ，宏参数 $\delta = 0.25$

最优收敛率: 梯度和旋度的 $L^2$ 误差达到预期的 $O(h)$ 收敛率
场本身收敛: 场 $\sigma_h$ 和 $\eta_h$ 显示出接近 $O(h^2)$ 的更好收敛率
稳定性: 使用ghost稳定化的格式条件数保持在合理水平