Numerical analysis of the Biot equations coupled to frictional contact mechanics

技术摘要：耦合摩擦接触力学的 Biot 方程数值分析

问题表述
本文研究了受流体饱和的孔隙粘弹性介质在受到刚性障碍物摩擦接触时的数学建模与数值近似问题。该物理系统由准静态 Biot 方程控制，该方程将线性弹性（具体为 Kelvin-Voigt 粘弹性模型）与达西流（Darby flow）耦合在一起。接触力学通过边界上特定子集 ($\Gamma_c$) 的两个特定约束进行建模：

1. 法向顺应性（Normal Compliance）： 一种本构关系，其中法向接触应力与穿透深度成正比，从而避免了非光滑的 Signorini 条件。
2. 库伦摩擦（Coulomb Friction）： 一种控制切向运动的法则，当切向应力低于摩擦界限（与法向应力成正比）时，物体保持粘着状态；否则发生滑动。

由此产生的数学模型是一个耦合了线性偏微分方程（质量平衡方程）与非线性变分不等式（带有接触条件的动量平衡方程）的系统。主要的挑战在于流体与力学的耦合、问题的时变性以及摩擦接触条件引入的非线性。

方法论
作者针对该问题提出并分析了一种全离散数值方案：

• 空间离散： 使用保形有限元。位移空间使用 $k_1$ 阶拉格朗日单元进行近似，压力空间使用 $k_2$ 阶单元。
• 时间离散： 采用隐式欧拉法进行时间步进。
• 存在性证明策略： 为了建立离散解的存在性与唯一性，作者采用了一种新颖的理论方法。他们首先考虑一个简化的问题，其中接触项被“滞后”处理（实际上将摩擦模型简化为 Tresca 模型）。对于这个简化问题，他们利用固定应力分裂方案（fixed-stress splitting scheme）——这通常是用于解耦流体与力学的实用算法——作为一种理论工具，通过收缩映射论证来证明解的存在性。随后，应用来自接触力学文献的固定点论证来恢复完整的库伦摩擦模型。
• 误差分析： 误差估计的推导结合了现有关于 Biot 方程文献（Girault 等人，Phillips 和 Wheeler）以及随时间变化的变分不等式（Sofonea 和 Han）的技术。一个关键的分析步骤是证明流体与力学之间的耦合项在稳定性与误差分析中可以相互抵消，从而防止其降低估计精度。

主要结果

1. 适定性（Well-posedness）： 本文确立了离散解的存在性与唯一性，前提是时间步长满足与双线性形式的强制性常数及接触泛函的 Lipschitz 常数相关的特定条件。
2. 稳定性： 推导出了一个稳定性估计，表明只要时间步长足够小以处理接触项的非线性，离散解在不依赖于网格尺寸 ($h$) 和时间步长 ($\Delta t$) 的情况下保持有界。
3. 先验误差估计： 作者推导了位移在 $H^1$ 范数下和压力在 $L^2$ 范数下的误差估计。理论收敛阶数为：
   • 时间： 一阶收敛 ($O(\Delta t)$)。
   • 空间： 压力的收敛阶为最优阶 ($O(h^s)$)。位移的收敛阶相对于标准弹性问题降低了 $\sqrt{h}$ 倍，对于线性单元，其阶数约为 $O(h^{r-0.5})$。这种降低归因于接触条件降低了解的正则性。
4. 数值实验： 使用 FEniCSx 软件进行了数值实验。结果证实了理论预测：
   • 压力和位移均观察到一阶时间收敛。
   • 压力的空间收敛为最优。
   • 位移的空间收敛略有降低（线性单元的观测阶数约为 0.8），这与接触导致正则性降低的理论预期一致。

意义与主张
本文声称填补了孔隙弹性领域的一个特定研究空白。虽然关于带接触的弹性体以及不带接触的孔隙弹性已有大量文献，但将随时间变化的 Biot 方程与摩擦接触约束（特别是法向顺应性和库联盟摩擦）相结合的研究，此前尚未进行过先验误差估计的分析。

作者强调，其工作将此前已建立的（针对非多孔介质的）粘弹性接触问题的可处理分析，扩展到了耦合的孔隙粘弹性设置中。这项工作的意义在于：

• 为耦合流体与摩擦接触在多孔介质中的数值模拟提供了严谨的数学基础。
• 证明了在结合孔隙弹性与接触力学技术时，耦合项不会阻碍误差估计的推导。
• 提供了一个具有证明的稳定性与收敛率的全离散方案，并通过数值实验进行了验证。

本文在应用方面保持谦逊，严格专注于所提数值方案的数学分析与验证，而非提出新的物理应用或超出所推导估计范围的未来扩展。