期刊
JOURNAL OF HYDROLOGY
卷 551, 期 -, 页码 56-69出版社
ELSEVIER SCIENCE BV
DOI: 10.1016/j.jhydrol.2017.05.053
关键词
Variably saturated flow; Richards' equation; Model robustness; Picard iteration; Infiltration
资金
- National Natural Science Foundation of China [51609173, 51479144, 51522904]
- Key Laboratory for Groundwater and Ecology in Arid and Semi-Arid Areas, CGS [KLGEAS201601]
- Ministry of Land and Resources of China [201511047]
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. (C) 2017 Elsevier B.V. All rights reserved.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据