摘要

A particle-water discrete element based approach to describe water movement in partially saturated granular media is presented and tested. Water potential is governed by both capillary bridges, dominant at low saturations, and the pressure of entrapped air, dominant at high saturations. The approach captures the hysteresis of water retention during wetting and drainage by introducing the local evolution of liquid-solid contact angles at the level of pores and grains. Extensive comparisons against experimental data are presented. While these are made without the involvement of any fitting parameters, the method demonstrates relative high success by achieving a correlation coefficient of at least 82%, and mostly above 90%. For the tested materials with relatively mono-disperse grain size, the hysteresis of water retention during cycles of wetting and drainage has been shown to arise from the dynamics of solid-liquid contact angles as a function of local liquid volume changes.