玻尔百科几十年来，经典的有限元法（FEM）一直是计算工程的基石，但它在一个重要的约束下运行：计算网格必须与物体的几何形状完美对齐。这一要求，通常被称为“网格的束缚”，对于涉及移动或演化特征的问题，如扩展的裂纹或变化的流体界面，成为了一个主要的瓶颈，因为在这些问题中，持续且昂贵的网格重划分是不可避免的。本文探讨了一种革命性的替代方案：非贴体有限元法（Unfitted Finite Element Methods, UFEMs），这是一类旨在摆脱这一限制的技术。
通过允许物理几何体简单地“切割”固定的、非贴体的背景网格，UFEMs 为一类新的复杂问题提供了前所未有的灵活性和效率。本文深入探讨了使这种自由成为可能的优雅数学和计算思想。在第一部分“原理与机制”中，我们将探讨非贴体边界带来的核心挑战，并发现那些确保准确性和鲁棒性的巧妙解决方案——例如水平集函数、Nitsche 方法和鬼点罚函数稳定化。我们还将对比两种领先的哲学思想：CutFEM 和 XFEM。随后的“应用与跨学科联系”部分将展示这些方法的变革性影响，演示它们如何用于解决从断裂力学和岩土力学到拓扑优化和多相流的真实世界问题。
要真正领会非贴体有限元法的精妙之处，我们必须首先理解它们所解决的问题。想象一下，你是一名工程师，任务是模拟一个裂纹在复杂机械部件中的扩展过程。几十年来，解决此类问题的主力一直是有限元法（FEM）。在其经典形式中，FEM 是应用数学的杰作，但它遵循一个严格的规则：计算网格——将物体细分为如三角形或四面体等小型、简单形状的镶嵌——必须贴合所有的几何特征。单元的边缘必须与部件的边界、任何内部界面，以及最麻烦的裂纹本身完美对齐。
这个约束就是我们可能称之为“网格的束缚”的东西。随着裂纹在材料中生长和蜿蜒，你不得不在模拟的每一步都扔掉旧网格，生成一个全新的网格。这个过程，即网格重划分（remeshing），不仅计算成本高昂，而且众所周知地困难且容易出错。为复杂、演化的几何体生成高质量的网格是一门玄学，而将解数据从旧网格投影到新网格会引入不准确性。其成本可能高得惊人。对于一个在固定距离上扩展的裂纹，模拟步数随着我们加密网格而增加（例如，与 成正比，其中 是未知数的数量）。如果每一步都涉及昂贵的网格重划分和求解大型方程组，总计算量可能会急剧膨胀，使得高保真度模拟变得异常缓慢。如果我们能够摆脱这种束缚呢？如果我们能让物理过程在一个简单的、固定的背景网格上展开，而这个网格对其包含的复杂几何形状一无所知呢？
这就是所有非贴体方法核心的那个激进且解放性的想法。我们从一个简单的、结构化的背景网格开始——可以想象一个均匀的正方形网格或一个规则的三角形晶格——它填充了我们物体周围的一个包围盒。然后，我们简单地让我们物体的真实、复杂的几何形状“切割”这个网格。网格本身不改变。裂纹扩展，界面移动，但底层的网格保持固定和平静。
然而，这个美好的想法立即给我们带来了两个深远的挑战。首先，如果我们的物体边界不再与单元边缘对齐，我们如何描述它的位置，更不用说在其上施加如指定温度或位移等物理定律了？其次，当边界切掉一个单元极小的碎片时会发生什么？这些“小切割”是否会产生数值上的小恶魔，破坏我们模拟的稳定性？非贴体方法的发展史就是为这两个问题发现巧妙答案的历史。
让我们来解决第一个难题。为了处理一个在网格中看不到的边界，我们需要一张地图。这张地图通常由一个水平集函数（level-set function）提供，这是一个光滑的函数 ，其零等值线 定义了我们的边界或界面。例如， 可以是到边界的有符号距离，物体内部的点有 ，外部的点有 。通过简单地在我们背景网格的节点上评估这个函数并进行插值，我们创建了一个离散的、分片多项式的几何近似 。这个近似效果非常好：对于一个尺寸为 的网格上的线性插值，真实边界 和近似边界 之间的距离通常是 阶，而边界法向量的误差是 阶。
现在我们有了地图，但我们如何在这个浮动的边界 上施加像 这样的条件呢？我们不能简单地抓住节点并设置它们的值，因为 上很可能没有节点。答案在于一种非常灵活的技术，称为 Nitsche 方法。我们不是“强行”（通过直接代入）施加条件，而是“弱行”地施加，即修改我们方程的变分形式或弱形式。Nitsche 公式在边界 上增加了两种类型的积分项：
Nitsche 方法就像一组软弹簧，将解拉向边界上期望的状态。一个关键特性是它的对称性；得到的方程组保持对称，这在计算上是理想的。这种对称性是公式结构的代数属性，并不会因为使用近似的法向量 而被破坏，这是一个常见的误解。该方法提供了一种鲁棒且数学上可靠的方式来处理任意位置界面上的边界条件。
现在来看第二个，更隐蔽的难题。当我们的边界 从一个背景单元上切下一个微小的、薄片状的部分时会发生什么？这就是臭名昭著的“小切割”问题。一个有效体积相对于其总尺寸极小的单元在数值上会变得病态。定义在其上的基函数变得几乎线性相关，导致最终方程组出现灾难性的病态。方法的稳定性常数可能会退化，意味着解可能会被大的、非物理的振荡所污染。这个问题是非贴体方法早期发展的一个主要障碍。
解决方案是一种被称为鬼点罚函数稳定化（ghost penalty stabilization）的数值魔法 [@problem_id:2609375, @problem_id:2609389]。这个名字非常形象。这个稳定化项并不作用于物理边界 。相反，它作用于背景网格在切割单元附近的内部面上。它惩罚解的导数在这些“鬼”面上的跳跃。通过这样做，它将小切割部分上行为不佳、几乎不稳定的自由度与其在单元主体中稳定、行为良好的邻居弱耦合起来。这就像对一个摇摇晃晃的组件说：“你必须像你的邻居一样行事！”这强制了一定程度的光滑性，并防止解产生剧烈的振荡。
至关重要的是，无论边界如何切割网格，这种鬼点罚函数都能统一地恢复方法的稳定性。由此产生的方法是鲁棒和可靠的，并且系统[矩阵的条件数](@article_id:305575)不再受切割位置的影响。当设计得当时，这种稳定化也是一致的，这意味着它对于精确解会消失，并且不会破坏方法的精度阶数。这种从内部、远离问题边界本身进行稳定化的聪明想法是一个分水岭，它使得鲁棒、高阶的非贴体方法成为现实。
有了这些核心工具——用于几何的水平集、用于边界条件的 Nitsche 方法和用于稳定性的鬼点罚函数——构建完整的数值方法的舞台已经搭建好。两种主流哲学应运而生，为实现将模拟从网格中解放出来的同一目标提供了不同的路径。
切割有限元法（CutFEM）从字面上拥抱了切割的思想。它使用传统 FEM 中标准的、简单的多项式基函数。然而，当一个界面 将域划分为 和 时，CutFEM 认为函数空间是“断开的”。它在 部分的网格上定义一组多项式函数，并在 部分的网格上定义一组完全独立的函数。在被界面切割的单元上，这意味着你有重复的自由度，从而允许数值解在界面上有一个自然的跳跃。
在这个世界里，界面条件（如解或其通量的跳跃）并不是内建于基函数中的。相反，它们都是使用 Nitsche 型项在界面 上弱施加的。小切割问题始终存在，因此鬼点罚函数稳定化是确保鲁棒性的一个必要组成部分。CutFEM 是一个通用而强大的框架，它依赖于一个简单的底层函数空间，将所有复杂性都放在变分问题的公式化中（即 Nitsche 和鬼点罚函数项）。
扩展有限元法（XFEM）遵循一条不同的、也许更具艺术雄心的路径。它不是断开函数空间，而是“富集”它。XFEM 的核心引擎是单位分解法（Partition of Unity Method, PUM）。标准的有限元形函数 有一个特殊的性质：在任何点 ，它们的和为一（）。XFEM 利用这个性质，将关于解的专门知识直接构建到基函数中。
想象一下，你知道你的解有一个特定的奇异特征，由一个富集函数 描述。这可能是一个跳跃、一个扭折或一个奇异点。在 XFEM 中，你通过简单地将标准形函数 乘以你的特殊函数 来创建新的基函数。得到的近似解看起来像： 乘积 将特殊行为 局限在节点 具有影响的区域。你只为靠近特征的节点 添加这种富集。这就像给你标准的乐高积木一套特殊用途的附件，你只在需要的地方使用它们。
对于断裂力学，XFEM 以其强大的功能而闻名。使用了两种类型的富集：
所以我们有两种优雅的方法：CutFEM 使用简单的函数但复杂的公式来连接断开的部分，而 XFEM 使用更复杂的、富集的函数空间，这个空间已经了解了解的特殊特征。
这种新获得的自由并非没有代价。非贴体方法的威力依赖于面对和解决一系列在传统 FEM 中不存在的新挑战。
首先，是积分问题。例如，在 XFEM 中出现的被积函数远非标准数值求积法则所设计的光滑多项式。在裂纹附近，被积函数可能是不连续的，甚至是奇异的。对一个行为像 的函数应用标准的高斯求积法则简直是灾难；误差将是巨大的，并且不会随着网格加密而减小。解决方案需要更复杂的积分技术。要么必须将被切割的单元划分为被积函数光滑的子三角形，要么设计专门的坐标变换或定制的求积法则，这些法则是明确为处理奇异性而量身定做的。
其次，一个被称为混合误差（blending error）的微妙问题可能会出现在 XFEM 中。在富集区域的外边缘，存在一些节点被富集而另一些没有的单元。这些被称为混合单元（blending elements）。在这些单元中，允许富集空间完美表示富集函数 的单位分解“技巧”失效了。局部基函数不再能形成一个完整的单位分解，导致逼近能力的丧失和收敛速度的下降。这需要进一步的修正，例如使用“斜坡函数”（ramp functions），使富集在富集区域的边界处平滑地衰减到零，以恢复最佳精度。
我们的旅程始于渴望逃离贴体网格的束缚。通过拥抱在固定背景网格上切割的简单思想，我们踏上了一条引导我们发现一系列美丽而强大的数学思想的道路。我们学会了用水平集来描述隐藏的几何，用 Nitsche 方法在稀薄的空气中施加边界条件，以及用鬼点罚函数驱除小切割的恶魔。我们看到了这些工具如何被组装成不同但相关的哲学，如 CutFEM 和 XFEM，它们要么使用巧妙的公式，要么使用富集的函数来捕捉复杂的物理现象。
这些方法代表了计算科学与工程领域的一次范式转变。虽然详细分析可能表明，完整模拟的总渐近成本有时可能与网格重划分方法相似，但其实际效益是变革性的。通过消除复杂、脆弱且易错的网格重划分过程，非贴体方法为模拟以前无法想象的复杂问题打开了大门：从流固耦合的复杂舞蹈到材料的灾难性破坏，再到未来主义设计的拓扑优化。它们为我们提供了一种新的、更灵活、更强大的语言来描述物理世界。
现在我们已经拆解了非贴体有限元法的内部构造，看到了富集和水平集的齿轮如何转动，真正的乐趣开始了。科学中的一个强大思想不仅仅是放在架子上供人欣赏的美物；它是一把能打开新大门的钥匙。其力量的真正衡量标准是它允许我们进入的房间的数量和种类。从贴体网格的束缚中解放出来后，让我们开始一场冒险，探索现在我们力所能及的广阔的科学和工程问题新领域。
这些方法的最初动机是棘手的断裂问题。裂纹天生具有拓扑上的挑战性。它们出现、生长、转向、停止，有时还会分叉。对于要求网格贴合几何的传统有限元法来说，移动的裂纹是一个反复出现的噩梦，需要持续且昂贵的网格重划分。非贴体方法彻底改变了游戏规则。
那么，我们如何知道我们这个花哨的新工具是否好用呢？我们做任何优秀科学家都会做的事：在老牌大师的主场上，用它来挑战他们。在断裂力学中，一个告诉我们驱动裂纹所需能量的关键量是 -积分。几十年来，工程师们一直使用巧妙的、专门的技术，比如带有奇特的“四分之一点”单元的网格，来高精度地计算这个值。当我们将扩展有限元法（XFEM）应用于同一问题时，我们发现了一个美妙的结果：用 XFEM 计算的 -积分值与理论预测以及那些经典、可信方法的结果非常吻合。这不仅适用于简单的弹性材料，也适用于裂纹尖端附近的材料开始屈服和塑性变形的更复杂情况。这让我们相信，我们的方法不仅仅是一个聪明的技巧；它是一种合法、强大且准确的科学工具。有了这个工具，我们就可以轻松地提取其他重要的断裂参数，如能量释放率 以及单个的应力强度因子 和 ，这都得益于像域相互作用积分这样的优雅技术。
但世界比材料中一个简单的空洞要复杂得多。在微观层面，当裂纹试图张开时，通常有未断裂的韧带和分子力抵抗分离。这种行为由“内聚区模型”捕捉，这些模型将裂纹面描述为被一种取决于它们分开距离的牵引力法则维系在一起。UFEMs 以非凡的优雅处理了这个问题。我们只需富集近似解以允许位移跳跃，然后在我们的方程中增加一项，以解释这些内聚力所做的功。这使我们能够以更高的保真度模拟断裂过程，捕捉到简单模型无法看到的现象。
然而，真正的魔力出现在物体开始移动时。裂纹的静态图片是一回事；其生命过程的模拟是另一回事。想象一个平板中的裂纹突然受到拉伸。它会生长吗？多快？朝哪个方向？在动态 XFEM 模拟中，计算机在每个微小的时间步计算裂纹尖端的能量释放率 。然后，它参考一个材料定律，该定律告诉它材料的断裂韧性 ，这个值甚至可能取决于裂纹的速度 。模拟接着强制执行动态断裂的基本法则：只有当可用能量等于所需能量时，裂纹才会生长，即 。通过在每一步求解这个简单的方程，模拟确定了裂纹的速度，并根据应力确定其方向。然后更新水平集函数，裂纹延伸，整个过程重复。我们不再仅仅是分析一个结构；我们正在目睹其破坏过程随时间展开。
这个过程的最终结局是什么？对于一个正在扩展的裂纹来说，通常是分叉——一个单一裂纹分裂成两个或更多，使材料破碎的戏剧性时刻。这是终极的拓扑挑战，对于基于网格的方法来说绝对是恐怖。对于非贴体方法，其简单性几乎令人扫兴。当裂纹尖端的物理学规定应该形成一个分支时，模拟只是简单地引入一个新的水平集函数来表示新的裂纹分支。富集函数在局部被复制，每个分支一个，模拟在同一个固定的背景网格上继续进行。我们可以观察到一个裂纹一次又一次地分叉，形成一个复杂的断裂网络，而所有这些都无需经历网格重划分的计算痛苦。这就是 UFEMs 提供的自由。
不连续性的概念是普适的，一旦我们有了处理它的工具，我们就会开始在各处看到它。“裂纹”只是分离不同物质状态或场的界面的一种表现形式。
让我们把目光投向脚下的大地。在岩土力学中，当土壤或岩石在剪切作用下破坏时，通常会形成一个“剪切带”——一个强烈的、局部化变形的薄区域。这不是一个开放的空隙，但它是位移场中的一个不连续性。我们可以精确地将其建模为一个内聚界面，就像我们复杂的断裂模型一样，并使用 XFEM 来捕捉其物理过程。在一些简单的情况下，该方法是如此强大，以至于一个富集单元就可以完美地复制整个地质层的行为，这表明物理现象真正被富集所捕捉，而不是通过精细的网格。
现在让我们转向流体的世界。气泡与周围水之间的边界是一个移动的界面。这个界面由表面张力的微妙物理学所支配，根据杨-拉普拉斯定律，它会在边界上产生压力（或扩散通量）的跳跃。一个非贴体 FEM 可以被公式化来处理界面上流体性质（如粘度）的跳跃，以及这种特殊的通量跳跃条件。这为模拟复杂的多相流打开了大门，从反应器中上升的气泡到生物系统中细胞和囊泡的动力学。
或者考虑航天器重返大气层时的炽热下降。强烈的热量导致隔热罩的材料烧蚀或燃烧掉。航天器的表面是一个移动的边界，但它不是分离两种材料——它正在退缩并消失。这是一个经典的“斯特凡问题”，它也可以用 UFEM 框架优雅地建模。水平集跟踪后退的表面，而能量平衡的物理学，包括烧蚀的潜热，被施加在这个移动的边界上。航天器结构的底层网格保持固定，而其表面则优雅地消失在计算的虚空中。
在所有这些情况中——剪切带、气泡和烧蚀表面——故事都是一样的：一个强大的、通用的处理固定网格上界面的思想，使我们能够跨越学科，解决那些表面上看起来彼此关系不大的问题。
到目前为止，我们已经用我们新获得的自由来分析世界本来的样子。但是，一个新工具最深远的应用往往是创造前所未有的东西。
想象一下，你得到一块坚固的材料，被要求用有限的材料量，从中雕刻出最坚固的结构来支撑某个载荷。这就是“拓扑优化”问题。几十年来，这个问题的解决方案都是“模糊的”，将结构表示为一种灰度密度场。但有了 UFEMs，我们可以将材料和空隙之间的边界表示为一个由水平集定义的清晰、明确的界面。然后，我们可以让计算机在数学优化算法的引导下移动这个边界，演化形状。边界可以合并、分裂、形成孔洞。从网格的束缚中解放出来后，优化可以发现复杂的、有机的、通常是非直观的设计，这些设计拥有非凡的效率和优雅。我们正在使用 UFEMs 不是为了分析一个设计，而是为了演化一个设计。
我们可以将这个创造过程再向前推进一步，深入到微观层面。现代复合材料（如碳纤维）的性能取决于其组成部分——纤维和基体——的复杂排列。要设计一种新材料，我们首先需要从其微观结构预测其宏观性能（如刚度或导电性）。这就是“均匀化理论”的目标。挑战在于微观结构在几何上可能极其复杂。使用 UFEMs，我们可以建立一个虚拟实验室。我们创建一个微观结构的代表性“单元”，由于我们的方法是非贴体的，我们可以轻松地表示即使是最复杂的纤维或夹杂物排列。然后我们在这个单元上进行虚拟测试，以计算宏观材料的有效性能。这使我们能够在计算机上快速设计和测试新材料，甚至在实验室制造它们之前。我们不再仅仅是设计一个零件；我们正在设计构成它的物质本身。
从桥梁的灾难性倒塌到生物反应器中无声的气泡，从飞机机翼的设计到新型复合材料的发明，连接这些不同世界的线索是一个单一、统一的数学思想的力量。非贴体有限元法的美妙之处不仅在于其公式的巧妙，更在于它们赋予我们的深远自由。通过将对物理世界的描述与计算网格的人为构造分离开来，我们可以自由地跟随物理学走向任何地方，无论旅程在拓扑上变得多么复杂。这是一个有力的提醒：在科学和工程中，最深刻的洞见和最通用的工具往往来自于看到隐藏在复杂世界表面之下的简单、统一的原则。