[论文解读] Improving Newton's method performance by parametrization: the case of Richards equation

Konstantin Brenner, Clément Cancès|arXiv (Cornell University)|Jul 6, 2016

Advanced Numerical Methods in Computational Mathematics参考文献 42被引用 54

一句话总结

本论文通过基于参数化的重构方法，对 Richards 方程进行重新表述，以稳定求解非饱和多孔介质流中出现的退化非线性系统时 Newton 法的收敛行为。通过引入一种新的主要变量（基于通量的参数化），该方法确保了二次收敛和鲁棒的质量守恒，相较于标准的压力-饱和度公式，在迭代次数和精度方面表现显著更优，尤其在网格细化和材料属性变化的情况下表现突出。

ABSTRACT

The nonlinear systems obtained by discretizing degenerate parabolic equations may be hard to solve, especially with Newton's method. In this paper, we apply to Richards equation a strategy that consists in defining a new primary unknown for the continuous equation in order to stabilize Newton's method by parametrizing the graph linking the pressure and the saturation. The resulting form of Richards equation is then discretized thanks to a monotone Finite Volume scheme. We prove the well-posedness of the numerical scheme. Then we show under appropriate non-degeneracy conditions on the parametrization that Newton\\^as method converges locally and quadratically. Finally, we provide numerical evidences of the efficiency of our approach.

研究动机与目标

• 为解决在求解离散化 Richards 方程产生的非线性系统时，Newton 法收敛性能差的问题。
• 在渗透率和饱和度函数出现退化的情况下，提升 Newton 法的鲁棒性与效率。
• 通过压力-饱和度关系的新型参数化，确保质量守恒与二次收敛。
• 通过数值实验展示基于通量的（τ）公式相较于标准压力基（u）公式的优越性。

提出的方法

• 通过引入新主变量 τ（即饱和度函数的倒数）来重构 Richards 方程，以稳定 Newton 迭代。
• 采用单调的有限体积格式对参数化后的 Richards 方程进行离散化，以确保稳定性和一致性。
• 在参数化满足适当的非退化条件下，证明离散格式的适定性。
• 在参数化满足光滑性与非退化性假设的前提下，建立参数化系统下 Newton 法的局部二次收敛性。
• 实现一种带线搜索的不精确 Newton 方法，利用参数化在每次迭代中保持质量守恒。
• 基于参数化残差的雅可比矩阵构建线性化系统，以求解 Newton 更新子问题。

实验结果

研究问题

• RQ1 对压力-饱和度关系进行参数化，是否能改善 Newton 法在 Richards 方程求解中的收敛行为？
• RQ2 所提出的参数化是否相比标准公式具有更好的质量守恒性能？
• RQ3 在新公式下，Newton 法的性能如何受网格细化和材料参数 β 的影响？
• RQ4 参数化系统是否能在保持适定性与稳定性的同时维持二次收敛性？
• RQ5 在迭代次数与误差传播方面，τ 公式相较于 u 公式的相对效率如何？

主要发现

• 与 u 公式相比，τ 公式将达到收敛所需的 Newton 迭代次数减少了高达一个数量级，尤其在网格细化时表现显著。
• 当 ε = 10⁻⁸ 时，u 公式的相对饱和度误差主要受质量守恒误差主导（最高达 10⁻⁴），而 τ 公式的质量守恒精度达到 10⁻¹⁵ 量级。
• τ 公式在不同 β 值下均保持鲁棒性能，而 u 公式的迭代次数随 β 增大而显著上升。
• 理论分析证实，在参数化满足非退化条件时，参数化系统的 Newton 法具有局部二次收敛性。
• 数值结果表明，由于在选定参数区间内 s(τ) 的线性性，τ 公式可实现机器精度的质量守恒。
• 采用 τ 公式并取 ε = 10⁻¹⁶ 计算的参考解表现出可忽略的质量误差，验证了其作为高精度基准的适用性。