期刊
JOURNAL OF HYDROLOGY
卷 523, 期 -, 页码 119-127出版社
ELSEVIER SCIENCE BV
DOI: 10.1016/j.jhydrol.2015.01.053
关键词
Richards' equation; Finite volume predictor-correcter; Semi-implicit; Mass-conservative
资金
- US National Science Foundation, EPSCoR program [EPS-1135483]
- CI-WATER project
- EPSCoR
- Office Of The Director [1135483] Funding Source: National Science Foundation
Numerical solution of the Richards' equation (RE) in variably saturated soils continues to be a challenge due to its highly non-linear behavior. This is particularly true as soils approach saturation and the behavior of the fundamental partial differential equation changes from elliptic to parabolic. In this paper, a finite volume predictor-corrector method with adaptive time-stepping was developed to solve the 1D vertical RE. The numerical method was mass-conservative and non-iterative. In the predictor step, the pressure head-based form of the RE was solved using the cell-centered finite volume method and the pressure head was updated. In the corrector step, the soil water content was calculated by solving the mixed form RE. Five different schemes to evaluate the inter-cell hydraulic conductivity were investigated. The robustness and accuracy of the numerical model were demonstrated through simulation of experimental tests, including free drainage, field infiltration into wet and dry soils, and laboratory infiltration with falling water table. Numerical results were compared against laboratory measurements, simulation results from the Hydrus-1D program, or analytical solution when available. Results showed that the developed scheme is robust and accurate in simulating variably saturated flows with various boundary conditions. The arithmetic mean and Szymkiewicz's mean of inter-cell hydraulic conductivity performed better than other methods especially in the case of infiltration into very dry soil. (C) 2015 Elsevier B.V. All rights reserved.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据