Second order explicit splitting scheme for fluid-poroelastic structure interaction problems

技术摘要：用于流体–孔隙弹性结构相互作用问题的二阶显式分裂格式

问题陈述
本文研究了流体–孔隙弹性结构相互作用（FPSI）问题的数值模拟，具体针对固定区域内的时变 Stokes–Biot 系统。这些问题出现在诸如生物组织灌注和通过变形多孔材料的流动等背景中。该领域的一个核心挑战在于如何高效且稳定地处理界面耦合条件，即强制执行法向通量连续性和法向及切向应力平衡。虽然单块（monolithic）方法具有鲁棒性，但它们在每个时间步都需要求解大型耦合系统。分块（partitioned）方法在计算上更具吸引力，因为它们允许独立使用专门的流体（Stokes）和孔隙弹性（Biot）子问题求解器。然而，显式或松耦合的分块格式往往难以维持稳定性，特别是在界面耦合采用显式处理时。

方法论
作者提出了一种用于 Stokes–Biot 问题的全离散、二阶显式分裂格式。该方法论依赖于以下组成部分：

1. Robin 重构： 界面条件使用惩罚参数 $\gamma > 0$（切向）和 $L > 0$（法向）被重写为 Robin 形式。这种重构实现了子问题的解耦，同时保持了适用于分块算法的一致耦合机制。
2. 时间离散： 该方案采用二阶后向差分公式（BDF2）进行子问题方程的时间演化。
3. 显式外推： 为了实现一种全显式的、可并行化的格式，在 $t_{n+1}$ 时刻所需的 Robin 边界条件数据是通过基于先前时间步（$t_n$ 和 $t_{n-1}$）的二阶 Adams–Bashforth (AB2) 外推法进行近似。
4. 算法结构： 在每个时间步，Stokes 和 Biot 子问题被独立且并行地求解。结构位移与速度之间的运动学关系通过 Biot 子问题中的 BDF2 关系进行处理。

主要贡献
本研究的主要贡献在于理论与分析层面：

• 离散稳定性分析： 作者为所提出的格式提供了严格的稳定性证明。通过利用 BDF2 能量恒等式和外推界面项的精细分解，他们建立了一个封闭的稳定性界限。该界限在抛物型 Courant–Friedbachs–Lewy (CFL) 条件下成立，该条件将时间步长 $\Delta t$ 与网格尺寸 $h$ 联系起来（具体为 $\Delta t \leq c^* h^2$）。
• 先验误差估计： 本文利用基于投影的框架推导了先验误差估计。这包括：
  • 使用 Fortin 投影处理流体速度的散度自由约束。
  • 使用 Ritz 型投影处理孔隙弹性位移和孔隙压力。
  • 对由 BDF2 时间离散、AB2 界面外推以及投影运动学关系产生的相容性缺陷进行了仔细分析。
• 收敛率： 分析表明，流体速度、结构速度、孔隙压力和弹性位移的总误差在整体能量范数下被限制在 $C(h^k + \Delta t^2)$ 内，其中 $1 \leq k \leq 3$。这证实了二阶时间精度和最优阶空间收敛性。

数值结果
理论发现通过两个数值实验进行了验证：

1. 固定区域基准测试（制造解）： 使用已知精确解的测试用例来验证收敛率。结果证实了所有变量（流体速度/压力、结构位移/速度以及孔隙压力）均具有预测的二阶时间收敛性。在精细网格上的空间收敛性测试显示出与所选有限元空间（用于 Stokes 的 Taylor–Hood $P_2/P_1$ 和用于 Biot 变量的 $P_2$ 元）一致的最优阶率。
2. 移动区域示例（Navier–Stokes–Biot）： 为了展示该分块策略在固定域 Stokes–Biot 设置之外的适用性，作者将该方案应用于由耦合 Biot 系统的 Navier–Stokes 方程控制的 2D 血液流动问题。通过使用任意拉格朗日–欧拉（ALE）公式，模拟成功捕捉了压力波传播和结构变形，其结果与现有文献一致。

意义与主张
本文声称提供了一种数学上经过证明的、具有二阶精度且可并行的流体–孔隙弹性相互作用显式分块方法。其意义在于克服了 FPSI 问题中通常与显式耦合相关的稳定性挑战。通过证明在抛物型 CFL 条件下的稳定性并推导出严格的误差界限，这项工作确立了显式分裂可以在不牺牲数学严谨性的情况下兼顾效率与精度。

作者指出，尽管稳定性与误差分析严格针对固定域 Stokes–Biot 系统，但数值示例表明该策略足以扩展到更一般的移动域流体问题，例如涉及 Navier–Stokes 方程的问题。他们总结道，该方案在模块化、并行效率和可证明的精度之间提供了一个实际的折衷方案，尽管他们也谦虚地承认，在放宽 CFL 限制或将该框架扩展到几何非线性移动界面问题方面仍需进一步研究。