7.1 GENERAL
对于水或空气作为活动相的严格非饱和带流模拟通常是高度非线性的,需要小的时间步长和多次迭代来解决实际情况。在一些情况下,对于非饱和情况,使用MODFLOW的Picard迭代方案,或者对于非饱和流(空气或水)情况,使用MODFLOW-SURFACT/MODHMS(第2章)的修改Picard方案,或者在不需要进行严格非饱和流模拟的不受限制的含水层中,一步稳态模拟可能无法收敛。提供二次收敛的牛顿-拉夫逊方案可以极大缓解收敛困难。这里讨论的NRB1包将牛顿-拉夫逊方案与用于非线性情况的回溯算法相结合,整合到MODFLOW中,尤其在Picard迭代过多或完全失败时可以获益。
NRB1包包括牛顿-拉夫逊线性化和回溯方案,以稳定牛顿迭代。该公式设计与当前MODFLOW方案兼容,在每次迭代中解决所有节点的水头。
此外,该公式对当前由MODFLOW模拟的任何边界条件均不需要特殊处理。最后,回溯方案限制了任何迭代中残差的增加,同时欠驱动技术有助于解决解的迭代之间的振荡行为。这些问题将在本章的其余部分中讨论。
7.2 FORMULATION OF NEWTON RAPHSON SCHEME
MODFLOW,使用隐式时间离散化,通过矩阵方程表达多孔介质中的流体流动方程:
在这里,上标 k+1 和 k 分别表示当前和先前的时间水平,
Δt 是时间步长,
B 是对角线存储矩阵