A modified Picard iteration scheme for overcoming numerical difficulties of simulating infiltration into dry soil

作者：Zha, Yuanyuan; Yang, Jinzhong; Yin, Lihe; Zhang, Yonggen; Zeng, Wenzhi; Shi, Liangsheng^*

来源：Journal of Hydrology, 2017, 551: 56-69.

DOI：10.1016/j.jhydrol.2017.05.053

• Variably saturated flow
• Richards' equation
• Model robustness
• Picard iteration
• Infiltration

Numerical models based on Richards' equation are often employed to simulate the soil water dynamics. Among them, those Picard iteration models which use the head as primary variable are widely adopted due to their simplicity and capability for handling partially saturated flow conditions. However, it is well-known that those models are prone to convergence failure in some unfavorable flow conditions, especially when simulating infiltration into initially dry soils. Here we analyze the reasons that give rise to the numerical difficulty. Moreover, several modifications to the mass-conservative Picard iteration method are proposed so that numerical difficulty is avoided in these unfavorable flow conditions. Our proposed modifications do not degrade the simulated results, while they lead to more robust convergence performances and cost-effective simulations.

• 出版日期2017-8
• 单位武汉大学; 水资源与水电工程科学国家重点实验室