Numerical simulation of unstable preferential flow during water infiltration into heterogeneous dry soil
Name
water-12-00909-v2.pdf
Size
79.94 MB
Format
Adobe PDF
Checksum (MD5)
f7af5742638856995dea67d33ba743b0
Author(s)
Cueto-Felgueroso, Luis
Suarez-Navarro, María José
Fu, Xiaojing
Juanes, Ruben
Date Issued
March 23, 2020
Journal
Water
Publisher
Multidisciplinary Digital Publishing Institute
Citation
Cueto-Felgueroso, Luis, María José Suarez-Navarro, Xiaojing Fu, and Ruben Juanes, "Numerical simulation of unstable preferential flow during water infiltration into heterogeneous dry soil." Water 12, 3 (Mar. 2020): no. 909 doi 10.3390/w12030909 ©2020 Author(s)
Version
Final published version
Abstract
Water infiltration and unsaturated flow through heterogeneous soil control the distribution of soil moisture in the vadose zone and the dynamics of groundwater recharge, providing the link between climate, biogeochemical soil processes and vegetation dynamics. Infiltration into dry soil is hydrodynamically unstable, leading to preferential flow through narrow wet regions (fingers). In this paper we use numerical simulation to study the interplay between fingering instabilities and soil heterogeneity during water infiltration. We consider soil with heterogeneous intrinsic permeability. Permeabilities are random, with point Gaussian statistics, and vary smoothly in space due to spatial correlation. The key research question is whether the presence of moderate or strong heterogeneity overwhelms the fingering instability, recovering the simple stable displacement patterns predicted by most simplified model of infiltration currently used in hydrological models from the Darcy to the basin scales. We perform detailed simulations of constant-rate infiltration into soils with isotropic and anisotropic intrinsic permeability fields. Our results demonstrate that soil heterogeneity does not suppress fingering instabilities, but it rather enhances its effect of preferential flow and channeling. Fingering patterns strongly depend on soil structure, in particular the correlation length and anisotropy of the permeability field. While the finger size and flow dynamics are only slightly controlled by correlation length in isotropic fields, layering leads to significant finger meandering and bulging, changing arrival times and wetting efficiencies. Fingering and soil heterogeneity need to be considered when upscaling the constitutive relationships of multiphase flow in porous media (relative permeability and water retention curve) from the finger to field and basin scales. While relative permeabilities remain unchanged upon upscaling for stable displacements, the inefficient wetting due to fingering leads to relative permeabilities at the field scale that are significantly different from those at the Darcy scale. These effective relative permeability functions also depend, although less strongly, on heterogeneity and soil structure. ©2020 Keywords: unsaturated flow; infiltration; preferential flow; multiphase flow modeling; soil heterogeneity
MIT Department
Massachusetts Institute of Technology. Department of Civil and Environmental Engineering
Terms of Use
Creative Commons Attribution
Persistent DSpace Link
DOI of Published Version
https://doi.org/10.3390/w12030909
