[论文解读] Parameter-robust convergence analysis of fixed-stress split iterative method for multiple-permeability poroelasticity systems

Qingguo Hong, Johannes Kraus+2 more|arXiv (Cornell University)|Dec 31, 2018

Advanced Numerical Methods in Computational Mathematics参考文献 23被引用 5

一句话总结

本文提出了一种用于多渗透率poroelasticity（MPET）系统的固定应力分裂迭代方法，证明其线性收敛性，且收缩率与所有物理参数无关。该方法通过引入所有流体压力之和实现稳定化，确保参数鲁棒性收敛，数值实验验证了其相对于预处理MinRes求解器的优越效率和鲁棒性，适用于广泛参数变化。

ABSTRACT

We consider flux-based multiple-porosity/multiple-permeability poroelasticity systems describing multiple-network flow and deformation in a poro-elastic medium, sometimes also referred to as MPET models. The focus of the paper is on the convergence analysis of the fixed-stress split iteration, a commonly used coupling technique for the flow and mechanics equations in poromechanics. We formulate the fixed-stress split method in the present context and prove its linear convergence. The contraction rate of this fixed-point iteration does not depend on any of the physical parameters appearing in the model. The theory is confirmed by numerical results which further demonstrate the advantage of the fixed-stress split scheme over a fully implicit method relying on norm-equivalent preconditioning.

研究动机与目标

• 开发一种针对多渗透率poroelasticity（MPET）系统的鲁棒迭代耦合方法，使其在广泛物理参数变化下仍保持有效性。
• 证明在MPET背景下，固定应力分裂方法的线性收敛性，且收缩率与模型参数无关。
• 建立理论和计算证据，表明固定应力方案优于采用范数等价预条件的完全隐式方法。
• 将固定应力分裂框架从Biot模型扩展至具有多个流体网络的更复杂MPET系统。
• 通过在不同渗透率、刚度和储水系数下进行的数值实验，验证理论收敛速率。

提出的方法

• 将固定应力分裂方法表述为MPET系统的定点迭代，解耦流动方程与力学方程。
• 通过引入所有流体压力之和实现稳定化，确保鲁棒性，并支持参数鲁棒收敛性分析。
• 收敛性分析基于先前关于基于通量的MPET系统的稳定性结果，利用能量范数推导收缩估计。
• 采用变分形式推导MPET系统的弱形式，支持有限元离散化与迭代求解。
• 使用标准有限元方法对位移和压力进行离散，通量通过Darcy定律计算。
• 数值实验将固定应力方案与预处理MinRes求解器进行比较，测量残差降低10^8所需的迭代次数。

实验结果

研究问题

• RQ1 固定应力分裂方法能否扩展至具有多个流体网络的MPET系统，并证明其收敛性与物理参数无关？
• RQ2 何种稳定化策略可确保固定应力迭代在MPET系统中的参数鲁棒性？
• RQ3 固定应力方案的收敛率与采用范数等价预条件的完全隐式方法相比如何？
• RQ4 固定应力方法在渗透率、刚度和储水系数极端变化下是否仍保持高效与鲁棒？
• RQ5 理论收缩率估计能否在不同网格细化和参数范围内通过数值实验验证？

主要发现

• 固定应力分裂方法在MPET模型中表现出线性收敛性，收缩率与所有物理参数（包括渗透率、刚度和储水系数）无关。
• 对于四网络MPET模型，固定应力方案在所有参数变化下仅需10至14次迭代即可将残差降低10^8，即使在极端值（如λ·10^8或K₃·10^10）下亦然。
• 预处理MinRes求解器需要显著更多的迭代次数——最粗网格下高达45次——尤其在渗透率或刚度广泛变化时，表明其鲁棒性较差。
• 该方法在所有测试参数范围内均保持稳定与高效，包括λ（最高达基值的10^8倍）、K₃（最高达基值的10^10倍）和K（最高达基值的10^2倍）的极端值。
• 数值结果证实了理论收缩率估计，所有参数变化和网格细化下迭代次数均稳定在低水平（10至14次）。
• 在所有测试场景中，固定应力方案均优于采用范数等价预条件的MinRes方法，展现出更优的鲁棒性与效率。