无参数预条件如何处理近不可压缩多孔弹性与弹性问题？

锁定区域为何使线性系统接近奇异

论文研究线性多孔弹性的 Biot 模型，其中固体位移与流体压力相互耦合，空间上采用最低阶弱 Galerkin 方法离散，时间上采用隐式 Euler 格式。每个时间步都会产生一个需要用迭代法求解的鞍点线性系统。

困难来自锁定区域：随着 Lamé 参数 λ 增大，固体趋于不可压缩，位移的散度趋于零。此时散度项在第一个方程中占主导，在纯 Dirichlet 边界条件下，所得的代数系统接近奇异。论文还提到相关的 Poisson 现象：无法保持正确渐近行为的格式，通常会随着网格加密而出现收敛性恶化。

这一点重要，是因为无论离散格式多么有效地避免了锁定，接近奇异的系统对迭代求解器来说都很难处理。因此论文把问题分成两个：得到无锁定的离散格式，以及在 λ 很大时高效求解由此产生的线性系统。

正则化来自恒等式而非新增模型

作者引入一个数值压力变量，把原来的两场系统改写为三场格式。随后在第三个块方程中加入形如 −ρ w wᵀ z_h = 0 的项，其中 w 在论文中给出定义，ρ 是正则化常数。

关键在于，条件 wᵀ z_h = 0 是原格式本身固有的恒等式。因此，加入一个在精确解上取零的项，可以在不扰动待求解的情况下正则化代数系统。论文指出，正则化后的系统保持原解不变，并且所得 Schur 补的特征值被正数一致地上下界住，远离零。

正则化参数取为 ρ = C λ min(M°_p)，其中 C 是不依赖网格尺寸 h 和物理参数的常数。数值实验中作者使用 ρ = 0.1 λ min(M°_p)，也测试了 ρ = 0.01 λ min(M°_p) 和 ρ = λ min(M°_p)。由于该常数与问题本身绑定，而不是按网格逐个调节，这种方法被称为无参数预条件。

块 Schur 补预条件重新变得有效

针对正则化后的三场系统，论文采用标准的非精确块 Schur 补预条件。作者构造了两种变体：块对角预条件 P_d 和块三角预条件 P_t，二者都由前导块、中间块和对角压力块的近似构成。

分析表明，用 P_d 预条件时 MINRES 的收敛、用 P_t 预条件时 GMRES 的收敛，在收敛因子和渐近误差常数两方面都基本不依赖网格尺寸 h、锁定参数 λ 以及组合 κΔt。论文针对正则化系统给出这一结论，并说明三场预条件对应的证明因篇幅原因省略。

实现时需要求逆对角块。前导块和中间块用 CG 求解，并采用带阈值丢弃、丢弃容差为 10⁻³ 的不完全 Cholesky 分解；构造该分解时省略 ρ w wᵀ 项。压力块是对角矩阵，求逆很简单。作者报告 PCG 通常只需几次迭代即可收敛。

λ 与网格变化时迭代次数基本持平

二维和三维实验取 c₀ = 1、μ = 1、κ ∈ {10⁻⁶, 10⁻³, 1}、λ ∈ {1, 10², 10⁴, 10⁶}、Δt ∈ {10⁻⁶, 10⁻³, 10⁻²}，相对容差为 10⁻⁸，最多 1000 次迭代，初始猜测为零向量，GMRES 重启参数为 30。

在二维正则化结果中，无论 λ、Δt、κ 和 h 如何变化，GMRES 在所有测试情形下都约 25 次迭代收敛。MINRES 的迭代次数更高，但同样稳定：例如在 κ = 1、Δt = 10⁻⁶ 时，网格从 918 个单元增至 58,608 个单元，迭代次数从 44 变为 48；λ 从 1 升到 10⁶ 时，迭代次数从 44 变为 58。三维表格在更大的迭代次数上呈现同样规律，GMRES 约 24 至 45 次，MINRES 约 44 至 88 次。

论文指出，λ 从 1 增至 10² 时迭代次数略有变化，而在 λ = 10²、10⁴、10⁶ 时基本保持不变。作者还报告，求解时间随网格尺寸按 O(N^{3/2}) 增长，这与外层迭代次数基本不变、内层 PCG 迭代次数按 O(h⁻¹) 增长且每次迭代代价为 O(N) 的理论一致。

正则化主要稳定了 λ 增大时的 GMRES

论文在二维算例中比较了正则化系统与未正则化系统。当 λ = 1 时，未正则化格式所需的迭代次数与正则化格式几乎相同。λ 更大时，未正则化格式需要略多的迭代，且次数随 λ 增大而增加。

最明显的对比在 GMRES：正则化后它在所有测试情形下约 25 次迭代收敛，而未正则化时 GMRES 的迭代次数随 λ 增大而增加，λ 从 1 升到 10⁶ 时几乎翻倍。对 MINRES 而言，正则化带来的改善不如 GMRES 明显，但迭代次数比未正则化时更稳定。论文还报告，正则化后 GMRES 的收敛时间略短于未正则化。

这一比较是论文中对判断最有价值的结果：正则化不只是理论工具，它确实阻止了外层求解器的迭代次数随锁定参数增长，在三角预条件的 GMRES 路径上尤其明显。

渗透率跳变有测试，但分析假设 κ 为常数

第二个实验让渗透率发生跳变：在 0.25 ≤ x ≤ 0.75 上取 κ = κ₁，其余区域取 κ = 1，其中 κ₁ ∈ {10⁻¹⁰, 10⁻⁶, 10⁻²}，Δt = 10⁻⁶，λ ∈ {1, 10², 10⁴, 10⁶}，α ∈ {1, 10⁻², 10⁻⁴}。论文指出，由于使用三角网格，κ 的跳变可能出现在单元内部。

MINRES 和 GMRES 的收敛分析只适用于 κ 为常数的情形。块预条件在算法上仍适用于可变甚至不连续的 κ，该实验检验它们是否仍表现良好。在其他参数固定时，迭代次数随 κ₁ 变化基本不变，并且在网格加密下保持稳定。固定 κ₁ 时，λ 和 α 的变化只带来适度变化：在 58,608 个单元的网格上，α = 10⁻²、λ = 10⁶ 时 MINRES 需要 52 次迭代收敛，而 α = 10⁻⁴、λ = 1 时为 26 次。

作者还考察了内层 PCG 精度对外层求解器的影响。当内层 PCG 相对容差不大于 10⁻⁶ 时，MINRES 迭代次数保持稳定，这比外层 10⁻⁸ 的容差更宽松；为保证稳健性，他们把内层容差设为外层的 0.01 倍。论文把推广到高阶弱 Galerkin 离散和空间变化参数列为未来研究方向。

重点是保持解的修正而非新离散格式

论文的核心主张比通用求解方案更窄：对于用最低阶弱 Galerkin 方法和隐式 Euler 离散的线性多孔弹性，锁定区域的近奇异性可以通过取自原系统恒等式的正则化消除，且精确解保持不变。其收益是，标准非精确块 Schur 补预条件随后能让 MINRES 和 GMRES 的收敛基本不依赖网格尺寸和锁定参数。

指导读者判断的证据是正则化与未正则化运行的对比。正则化对 GMRES 帮助最大，在二维中 λ 从 1 到 10⁶ 时把它保持在约 25 次迭代，而未正则化的次数在该范围内几乎翻倍。MINRES 获益较少，在 λ = 1 时两种格式几乎无法区分。

适用范围是明确的：收敛分析覆盖 κ 为常数的情形，跳变实验只是数值验证；测试使用最低阶弱 Galerkin 方法与隐式 Euler，并对位移和压力都施加 Dirichlet 边界条件；正则化常数取为 λ min(M°_p) 的倍数，测试的几个倍数给出可比的迭代次数。在这一范围内，论文的优先事项很清楚：在不改变解的前提下修正代数上的近奇异性，再让标准块预条件承担求解工作。

本文引用的文献