Linear poroelasticity with solid incompressibility: consistent formulation and scalable numerical solution

技术摘要：具有固体不可压缩性的线性孔隙弹性理论

问题陈述
本文研究了固体相不可压缩的孔隙弹性介质的数值建模问题，这一条件在软组织和水凝胶的研究中十分重要。虽然全非线性孔隙弹性方程具有普适性，但其原初公式（primal formulations）往往存在不一致性问题，特别是在强制执行固体不可压缩性时，压力会随变形而变化。传统的处理方法通常依赖于总不可压缩性（$J=1$）或准不可压缩性，而作者旨在严谨地推导一种仅针对固体相不可压缩的模型，以避免以往非线性孔隙弹性原初公式中出现的不一致性。

方法论
作者提出了一种在固体不可压缩假设下对全非线性孔隙弹性方程进行一致性线性化的方法。该方法通过以下步骤进行：

1. 模型推导： 从具有拉格朗日孔隙率和亥姆霍兹势的非线性模型出发，作者通过拉格朗日乘子强制执行固体不可压缩性。为了避免由三阶项引起的不一致性以及计算 $L^2$ 乘子拉普拉斯算子的需求，作者采用了混合压力公式。这涉及定义一个“总压力”变量（$\mu$）并在平衡点附近对系统进行线性化。
2. 线性化系统： 生成的模型是一个四场系统，涉及：
   • 位移 ($d$)
   • 总孔隙压力 ($\mu$)
   • 孔隙率 ($\phi$)
   • 不可压缩性拉格朗日乘子 ($\lambda$)
     该系统通过隐式时间离散化被表述为一个半离散弱问题，从而产生一个双鞍点结构。
3. 数学分析： 作者确立了半离散和全离散方案的适定性（well-posedness）。他们验证了算子的有界性，并证明了无穷小（inf-sup）条件（Ladyzhenskaya–Babuška–Brezzi 条件）。一个关键发现是，离散的 inf-sup 条件仅要求拉格朗日乘子的测试空间 ($V_{\lambda, h}$) 是孔隙率空间的子空间 ($V_{\phi, h}$)。这使得使用最低阶有限元（FE）对（例如，位移采用 $P_1$ 连续，孔隙率采用 $P_0$ 不连续）成为可能，而无需复杂的稳定化或高阶单元。
4. 与 Biot 模型的等效性： 论文证明了所提出的四场系统在完全不可压缩极限下，在修改后的体积拉梅参数下，与标准的二场 Biot 模型在数学上是等效的。
5. 数值求解策略： 为了高效求解线性化系统，作者提出了一种基于 Schur 补形式化的可扩展块预条件子。他们利用固定应力迭代法（被解释为 Uzawa 迭代）来近似 Schur 补。该预条件子采用分块划分策略：
   • 位移块通过代数多网格（AMG）处理。
   • 剩余的块（压力、孔隙率、乘子）利用基于 AMG 的椭圆压力算子和用于质量矩阵的 Jacobi 迭代进行块 Jacobi 处理。

核心贡献

• 一致性公式： 一个严谨执行固体不可压缩性的线性化模型，解决了标准非线性孔隙弹性原初公式中的不一致问题。
• 最低阶稳定性： 该公式允许使用稳定的最低阶有限元近似（$P_1-P_0$ 型对），这在计算上非常高效，且无需界面拉格朗日乘子或复杂的稳定化项。
• 可扩展求解器： 开发了一种基于固定应力近似的鲁棒、可扩展预条件子，该预条件器能保持最优的迭代次数，且不受网格尺寸或问题维度的影响（已在 2D 和 3D 中得到测试）。
• 理论验证： 证明了所提离散方案的适定性和收敛阶数（$O(h^{k+1})$）。

结果
通过数值测试验证了理论主张：

• 等效性测试： 在稳态 Mandel 问题对比中，所提出的四场模型产生的计算结果与标准不可压缩 Biot 模型完全一致（差异在机器精度数量级内）。
• 收敛性测试： 在 2D 和 3D 中的制造解（manufactured solution）测试确认了位移、压力、孔隙率和拉格朗日乘子的理论收敛速率，涵盖了一阶（$k=0$）和二阶（$k=1$）方案。
• 可扩展性测试： 在 3D 溶胀基准测试中评估了预条件器的性能。随着自由度（DoFs）从约 1,000 增加到超过 3,300 万，平均 GMRES 迭代次数保持恒定（约为 60–70 次），展示了求解器的可扩展性。

意义
本文声称，这项工作为模拟不可压缩固体相的孔隙弹性提供了一个一致且高效的计算框架。通过直接从非线性方程推导模型并利用总压力变量，作者弥合了复杂非线性物理与实用、稳定的线性求解器之间的鸿沟。使用最低阶单元的能力降低了计算成本，而基于固定应力的预条件器确保了求解策略在处理大规模 3D 问题时的可扩展性。作者指出，虽然该模型目前是线性化的，但其框架旨在与底层非线性物理保持一致，未来的工作将探索将此策略应用于全非线性及耦合系统（例如热-孔隙弹性）。