🔢 数学

Hybrid CutFEM-IFEM model for embedded and immersed interfaces with triply-cut cells for the 2D Poisson problem

本文提出了一种用于二维泊松问题的混合 CutFEM-IFEM 方法，该方法将非拟合有限元与浸没界面技术相结合，以在单一网格上高效处理复杂的流体-流体和流体-固体界面，并利用三重切割单元的复合多项式表示和富集基函数，以实现与拟合网格相当的最佳收敛性和高精度。

原作者： Darsh Nathawani, Christopher Kees

发布于 2026-10-30

📖 1 分钟阅读🧠 深度阅读

CC BY 4.0

原作者： Darsh Nathawani, Christopher Kees

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

想象一下，试图模拟水流经复杂岩石结构时的流动情况，同时还要追踪两种不同性质水体（例如淡水和海水，或油和水）之间的边界。在工程学和物理学领域，这些相互作用不断发生，从海浪冲击防波堤到空气和水绕过船只的运动。为了理解这些现象，科学家使用数学模型来预测流体的行为。然而，一个主要的障碍始终是网格，即用于划分计算空间的微小形状组成的网格。传统上，这个网格必须经过完美的设计，以契合每一处岩石的曲线和每一条流体交界线。如果岩石移动或水线偏移，整个网格都必须重建，这是一个缓慢、脆弱且计算成本极高的过程。

为了解决这个问题，研究人员开发了“非匹配”（unfitted）方法，在这种方法中，网格保持固定，物体只需直接穿过网格即可。出现了两种强大的技术来处理这些切割：一种被称为浸没界面法（immersed interface method），它擅长捕捉两种流体之间属性的剧烈跳变；另一种被称为切分有限元法（cut finite element method），它非常适合处理流体停止的固体边界。长期以来，这两种方法各行其道。挑战在于如何将它们合并成一个单一且稳健的系统，能够处理在一个存在两种不同流体相互作用的区域内，还坐落着一个固体对象的情况，且无需重建网格。

在这项工作中，路易斯安那州立大学的研究人员成功地将这两种技术合并为一个混合模型。他们创建了一种模拟三相问题的新方法：两种不同的流体和一个固体对象，所有这些都在同一个不变的背景网格上进行相互作用。这项创新的关键在于如何处理那些同时被流体边界和固体对象切割的细碎、杂乱的网格单元。当一个网格单元同时被两条不同的线切割时，会产生一个复杂的、由三部分组成的形状。团队开发了一种利用“等效多项式”（equivalent polynomials）的数学技巧，用于计算这些棘手的“三重切割”单元内部的物理特性，而无需将其分解为更小的简单形状。这使得计算机可以准确地对方程进行积分，即使在几何形状高度不规则的情况下也是如此。

研究人员在一个包含圆形固体对象和一条分隔两种流体的直线场景中测试了他们的方法。他们使用了线性元素和二次元素进行模拟，这是用于近似解的不同数学复杂度水平。结果表明，他们的混合方法达到了该类数学模型所能达到的最高精度，与需要完美匹配网格的传统方法精度相匹配。至关重要的是，即使在两种流体之间的性质差异极大（材料常数之比高达一千）的情况下，这种高精度依然得以保持。这表明该方法足够稳健，可以处理现实世界中材料差异巨大的场景，例如空气与水或油与水。

研究还探讨了该方法如何处理不同类型的边界。当固体对象具有标准的“无滑移”（no-slip）条件时，该方法的表现非常完美。然而，当他们引入一种模拟热传递的更复杂的边界条件时，最复杂数学元素的精度受到了轻微限制。这种局限性并非由于新方法的缺陷，而是因为固体对象的曲线边缘在网格中是用直线表示的。研究人员指出，若要在未来实现完美的精度，需要对曲线进行更高级的表示。

最终，这项研究为那些需要模拟复杂、移动界面，却又不想承受不断重新划分网格负担的工程师和科学家提供了一个强大的工具。通过结合现有两种方法的优势并解决了三重切割单元的特定问题，该团队证明了在涉及多相的问题中实现最优精度和效率是可能的。该方法能够处理强间断和复杂几何形状，其稳定性足以媲美传统方法，为模拟海岸及环境工程中流体与固体的动态相互作用提供了一条可靠的路径。

技术摘要：用于嵌入式与浸没界面的混合 CutFEM-IFEM 模型

问题陈述
本文研究了涉及三种不同相态（两种流体相 $\Omega_a, \Omega_b$ 和一种固体相 $\Omega_s$ ）的二维 Poisson 问题的数值解。该区域包含两类界面：

浸没界面 ( $\Gamma_f$ )： 分隔两种流体相，其特征是材料属性（ $\mu$ ）、解 ( $u$ ) 以及通量 ( $\mu \partial u / \partial n$ ) 存在不连续性。
嵌入界面 ( $\Gamma_s$ )： 代表固体区域的边界，该区域从计算域中移除，需要施加边界条件（Dirichlet 或 Robin）。

挑战在于如何在单一、固定的背景网格上模拟这些相互作用，其中两个界面可能穿过同一个单元，从而产生“三重切割单元”（triply-cut cells）。传统的网格拟合方法在界面移动或变形时计算成本高昂且脆弱。因此，更倾向于使用非拟合方法（unfitted methods），但现有技术往往难以在同一个切割单元内同时处理强不连续性（通过浸没有限元法，IFEM）和嵌入边界（通过 CutFEM），特别是在积分和一致性方面。

方法论
作者提出了一种混合 CutFEM–IFEM 公式，旨在利用这两种方法在单一非拟合网格上的优势：

基函数与富集：
- 对于被流体-流体界面 ( $\Gamma_f$ ) 切割的单元，该方法采用 IFEM 基函数。这些基函数的构造使其在界面处连续，但在单元边上是不连续的，并精确满足跳跃条件。
- 为了处理非齐次跳跃（给定的跳跃值 $J_s$ 和 $J_f$ ），试函数通过特定的“提升”（lift）函数（用于解跳跃的 $j_h$ 和用于通量跳跃的 $q_h$ ）进行扩充。这些富集函数能够重现精确的跳跃，确保误差不依赖于不连续性的量级。
- 对于嵌入固体边界 ( $\Gamma_s$ )，使用 Nitsche 方法 弱化施加边界条件（Dirichlet 或 Robin），而不增加自由度。
弱形式与稳定性：
- 弱形式包含了用于嵌入边界的标准 Nitsche 项，以及用于针对小切割比例进行稳定的 ghost penalty 项。
- 针对三重切割单元引入了一项关键修改。由于固体边界移除了两个流体切割单元之间共享的部分边，标准的 IFEM 边积分会变得不一致。作者推导了特定的一致性修正项 ( $c(u_h, v_h)$ )，以补偿剩余流体部分边上的通量，从而恢复最优收敛性。
三重切割单元的积分：
- 本文扩展了等效多项式技术，以处理被多个界面切割的单元。
- 该方法并非将单元细分为子三角形，而是用等效多项式替换 Heaviside 函数和 Dirac delta 函数。
- 对于三重切割单元，作者根据切割情况（基于孤立节点以及是否存在三点）对八种可能的拓扑配置进行了分类。他们推导了用于计算复杂几何结构中各相（固体、流体 $a$ 、流体 $b$ ）内矩匹配过程（moment-fitting process）右端项向量的解析公式。

核心贡献

混合公式： 提供了一个统一的框架，在单一网格上结合了用于流体-流体界面的 IFEM 和基于 Nitsche 的用于流体-固体边界的 CutFEM，适用于线性 ( $p=1$ ) 和二次 ( $p=2$ ) 元。
三重切割积分： 扩展了等效多项式积分方法，使其能够处理被两个不同界面相交的单元，涵盖了包括三点在内的所有拓扑配置。
一致性修正： 提出了一种针对三重切割单元的新型修正项，该项考虑了由嵌入边界导致的截断边，从而确保了方法的一致性并实现最优收敛。
精确跳跃表示： 利用富集函数处理非齐次解跳跃和通量跳跃，使误差量级与跳跃大小脱钩。

数值结果
该方法通过具有不同系数对比度（最高 $\mu_b/\mu_a = 1000$ ）和跳跃条件（齐次与非齐次）的人造解进行了验证。

收敛性： 对于线性及二次元，该方法在 $L^2$ 范数下均实现了 $p+1$ 阶的最优收敛速率。
精度： 尽管使用单一背景网格，其精度与显式拟合两个界面的网格相当。
鲁棒性： 在高系数对比度下仍能保持稳定和精确。
局限性： 当在弯曲嵌入边界上应用 Robin 条件时，二次收敛速率 ( $p=2$ ) 会退化为二阶。作者将其归因于水平集函数中对弯曲边界的线性表示，指出需要更高阶的边界表示才能恢复 $p=2$ 的最优速率。

意义与主张
本文声称为复杂的多相问题提供了一种鲁棒且高效的数值工具，这类问题通常涉及界面的移动与相互作用。通过结合 IFEM 和 CutFEM，该方法避免了重新生成网格的需求，同时保持了高精度并能处理强不连续性。作者强调，这是首个在同一切割单元内，将处理浸没界面的 IFEM 处理方式与处理嵌入边界的 Nitsche 处理方式相结合的工作，特别解决了三重切割单元的积分与一致性挑战。结果表明，这种混合方法既具备拟合网格的精度，又具备非拟合方法的灵活性，使其适用于涉及移动流体-流体界面和复杂固体边界的应用场景。

技术摘要：用于嵌入式与浸没界面的混合 CutFEM-IFEM 模型

类似论文