Quadrature error estimates on non-matching grids in a fictitious domain framework for fluid-structure interaction problems

技术摘要：非匹配网格下虚构域框架内流固耦合问题的求积误差估计

问题陈述
本文研究了基于分布式拉格朗日乘子的虚构域框架下，流固耦合（FSI）问题的数值离散化。在该框架中，流体和固体分别使用两个完全独立的、非匹配的网格进行离散。固体在参考域 $B$ 上通过拉格朗日框架描述，并映射到固定欧拉流体域 $\Omega$ 内的当前位置。流体速度与固体变形之间的运动学耦合通过定义在固体参考域上的拉格朗日乘子 $\lambda$ 来实现。

一个关键的计算挑战在于组装耦合（界面）矩阵。该矩阵需要对流体基函数（经固体映射后）与固体基函数的乘积在固体单元上进行积分。由于这些函数定义在非匹配网格上，精确积分需要计算流体网格与映射后的固体网格之间的几何交集，以构建复合求积规则。这一过程计算成本极高，特别是由于耦合项必须在每个时间步都进行组装。或者，另一种近似方法是直接在固体单元上使用标准求积规则进行积分，而不进行网格交集处理，这会引入求积误差。本文研究了这种近似积分策略的理论有效性及其误差估计。

方法论
作者分析了一个稳态 Stokes 型 FSI 模型，但指出其结果对更一般的非线性情况也具有代表性。变分形式涉及寻找流体速度 $u$、压力 $p$、固体变形映射 $X$ 和拉格朗日乘子 $\lambda$，使之满足一个鞍点系统。

研究重点在于用于定义耦合项 $c(\cdot, \cdot)$ 的两种特定选择：

1. $c_0$（对偶配对）： 基于 $L^2(B)$ 内积，其中乘子空间是 $H^1(B)$ 的对偶空间。
2. $c_1$（标量积）： 基于 $H^1(B)$ 内积，其中乘子空间本身为 $H^1(B)$。

研究方法通过以下步骤进行：

• 离散化： 为所有变量引入有限元空间。流体采用 Bercovier–Pironneau 元（$P1$-iso-$P2/P1$），而固体映射和乘子采用连续分片线性元。
• 近似积分： 作者定义了离散双线性形式，其中耦合项使用固体网格上的标准求积规则进行计算，忽略了复合后的流体基函数相对于固体网格具有分片多项式性质且存在不连续性的事实。
• 抽象误差分析： 遵循 inf-sup 条件的框架，作者推导出了一个一般误差估计（定理 2），该估计界定了精确离散解与使用不精确积分得到的解之间的差异。该界限取决于求积误差泛函。
• 求积误差估计： 核心技术贡献在于证明了耦合项的特定估计。作者利用了复合函数 $v_h(X)$ 在固体单元上是连续且分片线性的这一事实，从而获得了特定的分数阶 Sobolev 空间（$H^{1+s}$）正则性属性。他们利用插值理论和反向不等式来界定误差。
• 稳定性证明： 作者证明了即使在耦合项通过近似构造的情况下，离散问题仍然是适定的（满足离散 inf-sup 条件），前提是满足特定的网格尺寸条件。

主要贡献与结果

1. 具有近似耦合的适定性： 本文证明了即使在耦合矩阵使用近似求积进行组装时，离散问题也是稳定且适定的，无需进行网格交集处理。
2. 求积误差估计：
   • 对于 基于 $L^2$ 的耦合 ($c_0$)，求积误差由 $O(h_B^{3/2} |\log h_B|)$ 界定。作者表明，只要固体网格尺寸 $h_B \to 0$，该近似方法即可产生最优收敛速率。
   • 对于 基于 $H^1$ 的耦合 ($c_1$)，由于梯度项的存在，误差分析更为复杂。误差界限包含依赖于网格尺寸比例的项。具体而言，实现最优收敛不仅要求 $h_B \to 0$，还要求比率 $h_B/h_\Omega \to 0$（其中 $h_\Omega$ 是流体网格尺寸）。如果保持 $h_B/h_\Omega$ 为常数，该方法可能会失去最优性。
3. 反向不等式： 作者推导并利用了特定的反向不等式，用以处理误差估计中拉格朗日乘数的范数，特别是在对偶范数情形下。
4. 数值验证： 文中提出了四个数值测试以验证理论估计：
   • 测试 1 & 2： 比较了 $c_0$ 和 $c_1$ 在不同网格细化策略（常数比例与递减比例）下的收敛速率。结果证实，$c_0$ 在近似积分下具有鲁棒性，而 $c_1$ 则需要 $h_B/h_\Omega \to 0$ 才能维持精度。
   • 测试 3： 使用了在界面处存在速度梯度跳跃的解（典型的 FSI 特征），证实了由于解的奇异性导致的预期次优收敛速率 $O(h^{1/2})$。
   • 测试 4： 研究了高阶元（$P2/P1$），结果显示虽然 $c_0$ 保持鲁棒，但 $c_1$ 仍需精确积分或特定的网格比例才能达到最优速率。

意义与主张
本文声称为在虚构域 FSI 方法中使用近似积分提供了严谨的理论依据，这种做法在计算效率方面经常被采用。作者证明了：

• “不精确”方法不会损害数值方案的稳定性。
• 基于 $L^2$ 的耦合 ($c_0$) 特别适合近似积分，因为它不需要流体网格和固体网格以特定的比例进行细化即可保持精度。
• 基于 $H^1$ 的耦合 ($c_1$) 虽然可能提供不同的稳定性特性，但它对网格提出了更严格的要求（$h_B/h_\Omega \to 0$），以避免求积误差主导解的误差。

作者总结道，其理论结果证实了其先前工作 [12, 13] 中的实验观察，为根据特定的耦合形式和网格配置来选择精确（网格交集）或近似（非交集）组装技术提供了坚实的数学基础。