Numerical Issues Arising in Determination of Interlayer Conductivities in Layered Unsaturated Soils
Publication: International Journal of Geomechanics
Volume 17, Issue 3
Abstract
The one-dimensional (1D) Richards equation is used to model infiltration flow in an unsaturated two-layer soil. Imposing continuities of both flux and pressure head at the interface yields a nonlinear equation for determining interface conductivities. The authors show that multiple solutions of this nonlinear interface equation may exist if the spatial discretization is not fine enough around the interface, in particular as sharp wetting fronts pass through the interface, or for flow across highly dissimilar materials. Three hydraulic models, the Gardner model (G), the Mualem-van Genuchten model (MvG), and the Fredlund-Xing-Leong-Rahardjo model (FXLR), are investigated to demonstrate the nonuniqueness of solutions to the interface problem. For the simplest G model, a full mathematical analysis of the interface problem in terms of two parameters depending on local hydraulic conditions and mesh size is possible. For more advanced models, the interface equation can only be analyzed numerically. In all cases, the interface equation exhibits multiple solutions under suitable conditions, and the impact on the overall numerical solutions is studied. In particular, both a fixed-point iteration and a Newton iteration are used to solve the interface equation numerically.
Get full access to this article
View all available purchase options and get full access to this article.
Acknowledgments
This study was supported by the National Science Foundation (NSF) under Grant CMMI-0825089. The opinions, conclusions, and interpretations expressed in this paper are those of the authors, and not necessarily of the NSF. There are no data-sharing issues, because all of the numerical simulations are provided in the figures produced by solving the 1D Richards equation. The initial setups of the numerical experiments are all stated in the paper and are also available from the authors on request.
References
Brown, P. N., Byrne, G. D., and Hindmarsh, A. C. (1989). “VODE: A variable-coefficient ODE solver.” SIAM J. Sci. Stat. Comput., 10(5), 1038–1051.
Brunone, B., Ferrante, M., Romano, N., and Santini, A. (2003). “Numerical simulations of one-dimensional infiltration into layered soils with the Richards equation using different estimates of the interlayer conductivity.” Vadose Zone J., 2(2), 193–200.
DVODE [Computer software]. Center for Applied Scientific Computer, Lawrence Livermore National Laboratory, Livermore, CA.
Dye, H. B. (2008). “Moisture movement through expansive soil and impact on performance of residential structures.” Ph.D. thesis, Arizona State Univ., Tempe, AZ.
Fredlund, D. G., and Xing, A. (1994). “Equations for the soil-water characteristic curve.” Can. Geotech. J., 31(4), 521–532.
Gardner, W. (1958). “Some steady-state solutions of the unsaturated moisture flow equation with application to evaporation from a water table.” Soil Sci., 85(4), 228–232.
Hills, R., Porro, I., Hudson, D., and Wierenga, P. (1989). “Modeling one-dimensional infiltration into very dry soils: 1. Model development and evaluation.” Water Resour. Res., 25(6), 1259–1269.
HYDRUS [Computer software]. PC-Progress s.r.o., Prague, Czech Republic.
Kirkland, M. R., Hills, R., and Wierenga, P. (1992). “Algorithms for solving Richards’ equation for variably saturated soils.” Water Resour. Res., 28(8), 2049–2058.
Leong, E., and Rahardjo, H. (1997). “Permeability functions for unsaturated soils.” J. Geotech. Geoenviron. Eng., 1118–1126.
Matthews, C. J., Braddock, R. D., and Sander, G. (2004). “Modeling flow through a one-dimensional multi-layered soil profile using the method of lines.” Environ. Model. Assess., 9(2), 103–113.
Romano, N., Brunone, B., and Santini, A. (1998). “Numerical analysis of one-dimensional unsaturated flow in layered soils.” Adv. Water Resour., 21(4), 315–324.
Schaudt, K., and Morrill, J. (2002). “On the treatment of heterogeneous interfaces in soil.” Geophys. Res. Lett., 29(10), 1–34.
Simunek, J., Huang, K., and van Genuchten, M. T. (1998). “The HYDRUS code for simulating the one-dimensional movement of water, heat, and multiple solutes in variably-saturated media.” Res. Rep. 144, U.S. Salinity Lab, Riverside, CA.
SVFlux [Computer software]. SoilVision Systems Ltd., Saskatoon, SK, Canada.
Srivastava, R., and Yeh, T. C. J. (1991). “Analytical solutions for one-dimensional, transient infiltration toward the water table in homogeneous and layered soils.” Water Resour. Res, 27(5), 753–762.
Vadose/W [Computer software]. GEO-SLOPE, Calgary, AB, Canada.
van Genuchten, M. T. (1980). “A closed-form equation for predicting the hydraulic conductivity of unsaturated soils.” Soil Sci. Soc. Am. J., 44(5), 892–898.
Zha, Y., Yang, J., Shi, L., and Song, X. (2013). “Simulating one-dimensional unsaturated flow in heterogeneous soils with water content-based Richards equation.” Vadose Zone J., 12(2).
Information & Authors
Information
Published In
International Journal of Geomechanics
Volume 17 • Issue 3 • March 2017
Copyright
© 2016 American Society of Civil Engineers.
History
Received: Sep 25, 2015
Accepted: May 27, 2016
Published online: Aug 19, 2016
Discussion open until: Jan 19, 2017
Published in print: Mar 1, 2017
Authors
Metrics & Citations
Metrics
Citations
Download citation
If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.