A family of mixed-mixed strongly conservative finite element methods for Biot's model of consolidation

技术摘要：一类用于 Biot 模型的高阶强保守混合-混合有限元方法

问题表述
本研究旨在解决 Biot 固结准静态模型的数值近似问题，该模型描述了线性多孔介质中流体流动与固体变形之间的相互作用。作者关注一种四场表述形式，其主要未知量为有效（弹性）应力张量 ($\sigma$)、固体位移 ($u$)、流体通量（达西速度，$w$）以及流体压力 ($p$)。

控制方程组由动量平衡方程、质量平衡方程、达西定律以及由胡克定律（以柔度形式表示）导出的本构关系组成。模型在有界 Lipschitz 区域 $\Omega \subset \mathbb{R}^d$ ($d=1,2,3$) 内进行分析。在隐式时间离散（例如后向欧拉法）之后，问题简化为一个静态的、两重扰动的鞍点问题。作者指出，在典型的软生物组织物理机制下——其特征为低渗透率和趋于零的受限特定储水常数——标准的数值方案往往会因为违反质量守恒原则而导致体积自锁和伪压力振荡。

方法论与变分框架
为了应对这些挑战，本文提出了一种新型的混合-混合有限元方法。其核心方法论包括：

1. 新型希尔伯特空间： 作者为位移-通量对 $(u, w)$ 引入了一个新的希尔伯特空间，记作 $\bar{V}$。该空间定义为 $\bar{V} := \{(u, w)^T : u, w \in L^2(\Omega; \mathbb{V}), \text{ 且 } \text{div}(u+w) \in L^2(\Omega)\}$。这种选择使得弱形式能够直接强制执行质量和动量守恒，而无需引入诸如总压力或总应力之类的辅助变量。
2. 强保守离散化： 该方法利用针对对称应力张量的 $H(\text{div div}) \cap H(\text{div})$-保连续有限元。这些单元基于 Hu, Ma, Zhang, Chen, 和 Huang 最近的构造，确保了 $\text{div } \sigma$ 的法向连续性得到强制执行。因此，离散应力近似在逐点意义上保持了角动量守恒。
3. 混合-混合结构： 离散化在力学子问题（使用应力作为主变量）和流体子问题（使用通量作为主变量）中均采用了混合表述。所构建的离散空间使得应力空间的散度精确映射到位移空间（具体而言，$\text{div } \Sigma_h = \text{BDM}_{k-1}$），且通量空间的散度精确映射到压力空间。这种结构确保了流体质量平衡在逐点意义上（强形式）得到满足。

主要贡献
本文提出了以下具体贡献：

• 新变分表述： 直接从基本物理原理（动量平衡、质量平衡、达西定律和胡克定律）推导出四场变分表述，无需重写方程或引入非物理变量。
• 适定性分析： 作者证明了该问题在新型希尔伯特空间 $\bar{V}$ 中的连续性以及离散问题的适定性，并利用了参数鲁棒的 inf-sup 条件。稳定性常数被证明与 Lamé 参数（特别是第一 Lamé 参数 $\lambda$）无关，这对于避免在近不可压缩极限下的自锁至关重要。
• 最优误差估计： 为该全保守混合-混合方法建立了最优的先验误差估计。这些估计在应力的 $H(\text{div})$ 范数、与 $\bar{V}$ 及 $L^2$ 范数相关的能量范数，以及压力的 $L^2$ 范数下均达到了最优。
• 压力后处理： 引入了一种高效的后处理程序，将压力近似提升至高一阶的分片多项式空间。该技术提高了压力近似的收敛阶，使其与应力和通量变量的收敛阶相匹配。

结果与主张
理论结果表明，该系列方法具有以下特性：

• 强保守性： 方法在逐点意义上保持角动量和流体质量平衡。
• 无自锁性： 误差界与 Lamé 参数 $\lambda$ 无关，确保了在受限储水常数趋于零和低渗透率情况下的稳健性能。
• 最优性： 推导出的误差估计确认了所有四个物理场的收敛速率均为最优。

意义
本文声称，这项工作为模拟 Biot 固结模型提供了一个稳健的理论和数值框架，特别是在涉及生物组织的实际应用中，这些应用的参数机制往往会导致标准方法的数值不稳定。通过在四场混合-混合框架内使用 $H(\text{div div})$-保连续单元，作者提供了一种既严格遵循物理守恒定律（质量和角动量），又能在广泛的物理参数范围内保持稳定性和最优精度的解决方案。这项工作将此前在 Stokes 和弹性问题中探索的强保守离散化方法扩展到了耦合的多孔弹性语境中。