The Tempered Finite Element Method

技术摘要：调温有限元法 (The Tempered Finite Element Method, TFEM)

问题陈述
有限元法 (FEM) 是计算物理和工程领域占主导地位的数值技术，因其处理复杂几何形状的能力（通过非结构化网格）而闻名。然而，标准的 FEM 公式在遇到特定的网格配置时会表现出脆性，特别是包含“帽状单元”（角度接近 180° 的三角形）排列成带状的情况。虽然孤立的退化单元（如针状单元或孤立的帽状单元）不会阻碍收敛，但帽状单元带会诱发“锁定”现象。在这种配置下，数值解会被迫沿该带呈线性分布，无论真实解如何，这导致有效网格尺寸由该带的长度而非其厚度决定。这导致标准 FEM 出现收敛失效。历史上，“最大角度条件”（要求所有角度远离 180°）被认为是收敛所必需的，尽管近期的研究表明对于孤立单元而言，这一条件并非严格必要。目前的挑战在于如何处理具有零测度或近乎退化元素的网格，这些元素在界面追踪方法（例如 X-Mesh）中非常相关，其中界面通过网格重传进行追踪，可能会产生帽状单元带。

方法论
作者提出了调温有限元法 (TFEM)，这是一种对经典 FEM 的改进，允许在含有帽状单元带的网格上实现收敛。其核心思想受冶金学中“调温”（tempering）过程的启发：正如调温可以缓解硬化钢的内应力以增加其柔韧性一样，TFEM 在退化带中局部“软化”材料属性，以允许数值解符合连续问题。

对于诸如泊松方程之类的椭圆问题，标准刚度矩阵 $K_{ij}$ 涉及单元映射雅可比行列式 $J$ 的倒数。在退化单元中，$J \to 0$，导致除以零或产生无限大的特征值。TFEM 通过使用参数 $J_{min}$ 将雅可比行列式限制在非零范围内，从而修改局部刚度矩阵：
$$\tilde{K}_{ij} = \int_{\hat{E}} \frac{1}{\max(J, J_{min})} [\dots] d\xi$$
这等效于在弱形式中引入一个可变扩散系数 $d(x)$：
$$d(x) = \begin{cases} 1 & x \in \Omega \setminus B \\ D = \frac{J}{J_{min}} & x \in B \end{cases}$$
其中 $B$ 是退化单元带。选择 $J_{min}$ 时需依赖于网格尺寸 $h$（具体为 $J_{min} \sim h^{d+1}$，其中 $d$ 是空间维度），以平衡降低扩散系数（以避免锁定）与保持近似性质之间的关系。

对于精确的零测度元素（即 $J=0$），作者证明了 TFEM 方案的极限对应于一种带有界面跳跃惩罚项的模态法 (mortaring method)。

核心贡献

1. 收敛性的理论证明： 论文为处理帽状单元带的 TFEM 提供了严谨的误差分析。证明了通过选择最优的调温参数 $D \sim \bar{h}/h^2$（其中 $\bar{h}$ 是带宽度），该方法可以实现 $O(h)$ 的最优 $H^1$ 半范数收敛，达到了标准 FEM 在规则网格上的收敛速率。
2. 与模态法的等价性： 对于精确零宽度的带，作者证明了 TFEM 等价于带有跳跃惩罚项的非协调模态法，为处理退化网格与界面耦合技术之间架起了理论桥梁。
3. 向复杂物理问题的扩展： 该方法成功地从标量拉普拉斯方程扩展到了：
   • 线性弹性力学： 证明了在跨界面应力分析中具有与标准 FEM 相当的精度。
   • 非协调网格模态法： 展示了 TFEM 可以无缝连接不同分辨率的网格，而无需过渡单元。
   • 高阶单元： 验证了只要 $J_{min}$ 随单元阶数进行适当缩放，该方法即可保持最优收敛率（最高达 4 阶多项式）。
   • 平流问题： 将 TFEM 应用于使用迎风 Petrov-Galerkin (SUPG) 稳定化的平流-扩散问题，并指出调温仅需应用于稳定化项。
4. 鲁棒性分析： 数值实验确定了该方法对于关系式 $J_{min} = C h^{d+1}$ 中的常数 $C$ 具有鲁棒性。在广泛的 $C$ 值范围内（跨越近一个数量级），误差仍保持接近最优，且线性系统的条件数保持在可控范围内，尤其是使用标准预条件子时。

结果

• 数值实验： 在具有制造解的 2D 和 3D 泊松问题中，带有帽状单元带的 TFEM 恢复了标准 FEM 在规则网格上的收敛速率。具体而言，在 2D 中选择 $J_{min} \sim h^3$，在 3D 中选择 $J_{min} \sim h^4$，可获得最优的 $L^2$ 和 $H^1$ 误差收敛。
• 误差敏感性： 与 $L^2$ 误差相比，$H^1$ 误差对 $J_{min}$ 的选择不那么敏感，在很宽的参数范围内表现出平台期。
• 条件数： 虽然随着带宽度的减小，刚度矩阵的条件数会增加，但使用不完全 Cholesky 预条件子可确保迭代求解器的收敛性不受带厚度的影响。
• 随机测试： 在 2D 环境下的 100 项随机测试证实，即使采用简单的选择（$C=1$），TFEM 的误差也与标准 FEM 在规则网格上的误差几乎一致。

意义与主张
论文声称，TFEM 通过放宽与退化单元相关的严格网格质量要求，显著增加了有限元分析的灵活性。它允许使用具有零测度或近乎退化元素的网格（例如由动态界面追踪生成的网格），而不会遭受使标准 FEM 失效的锁定现象。

作者强调，该方法极易在现有 FEM 代码中实现，仅需修改雅可比行列式的计算以防止除以零。他们将 TFEM 定位为不是作为高质量网格的标准 FEM 的替代品，而是作为一种鲁棒的扩展，使得使用特定的网格拓扑（帽状单元带）成为可能，而这些拓扑此前被认为无法工作。该方法对于像 X-Mesh 界面追踪方法这样的应用特别有价值，因为在这些应用中，维持共形网格是非常困难或计算成本高昂的。论文结论指出，虽然推导出了参数 $C$ 的理论最优选择，但该方法的实际鲁棒性允许使用广泛的可接受值，使其具有即时适用性。