Efficient Approach to Solving Transient Unsaturated Flow Problems. I: Analytical Solutions
Publication: International Journal of Geomechanics
Volume 17, Issue 7
Abstract
Richard’s equation governs the migration of moisture in the soil under unsaturated conditions. Although this differential equation provides a rigorous approach to simulating important infiltration problems, obtaining analytical and numerical solutions to this equation has been a particularly challenging task. This is largely due to the highly nonlinear nature of the soil hydraulic properties, including the moisture retention curve and the hydraulic conductivity function. Whereas analytical solutions of Richard’s equation have been reported for problems involving steady-state conditions and simple hydraulic models, solutions for transient conditions have rarely been obtained. However, such analytical solutions would be particularly valuable, for example, to validate the accuracy of numerical schemes, as well as to facilitate parametric evaluations. A series of analytical solutions of Richard’s equation for unsaturated flow under transient conditions have been developed as part of this study. The solutions involve a variety of initial and boundary conditions. The analytical solutions in this study could be obtained after expressing the governing equation as the addition of advective and diffusive flow components. The solutions consider logarithmic and linear models to represent the soil moisture retention and the hydraulic conductivity functions, respectively. Solutions are also provided for special cases in which either the advective or the diffusive components dominate the flow process, as well as for the steady-state cases. A parametric evaluation was found to provide insight into important characteristics of infiltration problems. In particular, relevant features of an unsaturated flow problem can be explained by evaluating the trends in its advective and diffusive flow components.
Get full access to this article
View all available purchase options and get full access to this article.
Acknowledgments
The authors acknowledge the support of the following agencies: the National Council for Scientific and Technological Development (CNPq) (Project 30449420127), the Coordination for the Improvement of Higher Level Personnel (CAPES) (Project 1431/14-5), the National Science Foundation (CMMI) (Grant 1335456), the University of Brasilia, and the University of Texas, Austin, for funding this research.
References
Basha, H. A. (1999). “Multidimensional linearized nonsteady infiltration with prescribed boundary conditions at the soil surface.” Water Resour. Res., 35(1), 75–83.
Bastian, W. C., and Lapidus, L. (1956). “Longitudinal diffusion in ion exchange and chromatographic columns. Finite column.” J. Phys. Chem., 60(6), 816–817.
Bear, J. (1979). Hydraulics of groundwater, McGraw-Hill, New York.
Biot, M. A. (1941). “General theory of three-dimensional consolidation.” J. Appl. Phys., 12(12), 155–164.
Brenner, H. (1962). “The diffusion model of longitudinal mixing in beds of finite length. Numerical values.” Chem. Eng. Sci., 17(4), 229–243.
Buckingham, E. (1907). “Studies on the movement of soil moisture.” Bull. 38, USDA, Bureau of Soils, Washington, DC.
Cavalcante, A. L. B., and Farias, M. M. (2013). “Alternative solution for advective-dispersive flow of reagent solutes in clay liners.” Int. J. Geomech., 49–56.
Cavalcante, A. L. B., Ozelim, L. C. S. M., Swamee, P. K., and Rathie, P. N. (2013). “Explicit numerical iterative methods applied to the three-parameter infiltration equation.” Soils Rock, V, 36(1), 129–133.
Cavalcante, A. L. B., and Zornberg, J. G. (2017). “Efficient approach to solving transient unsaturated flow problems. II: Numerical solutions.” Int. J. Geomech, 04017014.
Chen, G. J., and Gallipoli, D. (2004). “Steady infiltration from buried point source into heterogeneous cross-anisotropic unsaturated soil.” Int. J. Numer. Anal. Methods Geomech., 28(10), 1033–1055.
Chen, J.-M., Tan, Y.-C., and Chen, C.-H. (2003). “Analytical solutions of one-dimensional infiltration before and after ponding.” Hydrol. Processes, 17(4), 815–822.
Chen J.-M., Tan, Y.-C., Chen, C.-H., and Parlange, J. Y. (2001). “Analytical solutions for linearized Richards eqution with arbitrary time-dependent surface fluxes.” Water Resour. Res., 37(4), 1091–1093.
Cleary, R. W., and Adrian, D. D. (1973). “Analytical solution of the convective-dispersive equation for cation adsorption in soils.” Soil Sci. Soc. Am. J., 37(2), 197–199.
Dell’Avanzi, E., Zornberg, J. G., and Cabral, A. (2004). “Suction profiles and scale factors for unsaturated flow under increased gravitational field.” Soils Found., 44(3), 79–89.
Fityus, S. G., and Smith, D. W. (2001). “Solution of the unsaturated soil moisture equation using repeated transforms.” Int. J. Numer. Anal. Mech. Geomech., 25(15), 1501–1524.
Gardner, W. R. (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.
Gershon, N. D., and Nir, A. (1969). “Effect of boundary conditions of models on tracer distribution in flow through porous mediums.” Water Resour. Res., 5(4), 830–840.
Ghotbi, A. R., Omidvar, M., and Barari, A. (2011). “Infiltration in unsaturated soils—An analytical approach.” Comput. Geotech., 38(6), 77–782.
Hillel, D. (2004). Introduction to environmental soil physics, Elsevier Science, San Diego, 494.
Hogarth, W. L., and Parlange, J. Y. (2000). “Application and improvement of a recent approximate analytical solution of Richards’ equation.” Water Resour. Res., 36(7), 1965–1968.
Hogarth, W. L., Parlange, J. Y., and Braddock, R. D. (1989). “First integrals of the infiltration equation: 2. Nonlinear conductivity.” Soil Sci., 148(3), 165–171.
Hogarth, W. L., Parlange, J. Y., and Norbury, J. (1992). “Addendum to ‘First integrals of the infiltration equation.’” Soil Sci., 154(5), 341–343.
Huang, R. Q., and Wu, L. Z. (2012). “Analytical solutions to 1-D horizontal and vertical water infiltration in saturated/unsaturated soils considering time-varying rainfall.” Comput. Geotech., 39, 66–72.
Jaiswal, D. K., Kumar, A., and Yadav, R. R. (2011). “Analytical solution to the one-dimensional advection-diffusion equation with temporally dependent coefficients.” J. Water Resour. Prot., 3(1), 76–84.
Lapidus, L., and Amundson, N. R. (1952). “Mathematics of adsorption in beds. VI. The effect of longitudinal diffusion in ion exchange and chromatographic columns.” J. Phys. Chem., 56(8), 984–988.
Lindstrom, F. T., Haque, R., Freed, V. H., and Boersma, L. (1967). “The movement of some herbicides in soils. Linear diffusion and convection of chemicals in soils.” Environ. Sci. Technol., 1(7), 561–565.
Mason, M., and Weaver, W. (1924). “The settling of small particles in a fluid.” Phy. Rev., 23(3), 412–426.
Mathematica 10.0 [Computer software]. Wolfram Research, Inc., Champaign, IL.
Narasimhan, T. N. (2004). “Darcy’s law and unsaturated flow.” Vadose Zone J., 3(4), 1059.
Ogata, A., and Banks, R. B. (1961). “A solution of the differential equation of longitudinal dispersion in porous media.” Professional Paper 411-A, A1-A9, USGS, Reston, VA.
Parlance, M. B., Prasad, S. N., Parlange, J. Y., and Romkens, M. J. M. (1992). “Extension of the Heaslet-Alksne technique to arbitrary soil water diffusivities.” Water Resour. Res., 28(10), 2793–2797.
Parlange, J. Y., et al. (1997). “New approximate analytical technique to solve Richards equation for arbitrary surface boundary conditions.” Water Resour. Res., 33(4), 903–906.
Philip, J. R. (1960). “General method of exact solution of the concentration-dependent diffusion equation.” Aust. J. Phys., 13, 1–12.
Philip, J. R. (1969). “Theory of infiltration.” Advances in hydroscience, V. T. Chow, ed., Vol. 5, Academic, New York, 215–296.
Rathie, P. N., Swamee, P. K., Cavalcante, A. L. B., de S. M. Ozelim, L. C. (2012). “Lagrange’s inversion theorem and infiltration.” Int. J. Math. Comput. Phys. Electr. Comput. Eng., 6(7), 386–391.
Ross, P. J., and Parlange, J. Y. (1994). “Comparing exact and numerical solutions of Richards’ equation for one-dimensional infiltration and drainage.” Soil Sci., 157(6), 341–344.
Sander, G. C., Parlange, J. Y., Kühnel, V., Hogarth, W. L., Lockington, D., and O’Kane, J. P. J. (1988). “Exact nonlinear solution for constant flux infiltration.” J. Hydrol., 97(3–4), 341–346.
Swamee, P. K., Rathie, P. N., de S. M. Ozelim, L. C., and Cavalcante, A. L. B. (2014). “Recent advances on solving the three-parameter infiltration equation.” J. Hydrol., 509, 188–192.
Tartakovsky, D. M., Neuman, S. P., and Lu, Z. (1999). “Conditional stochastic averaging of steady state unsaturated flow by means of Kirchoff transformation.” Water Resour. Res., 35(3), 731–745.
Terzaghi, K. (1943). Theoretical soil mechanics, John Wiley and Sons, New York.
Wang, Q. J., and Dooge, J. C. I. (1994). “Limiting cases of water fluxes at the land surface.” J. Hydrol., 155(3–4), 429–440.
Yeh, T. C. J. (1989). “One-dimensional steady state infiltration in heterogeneous soils.” Water Resour. Res., 25(10), 2149–2158.
Zhu, J., and Mohanty, P. (2002). “Analytical solutions for steady state vertical infiltration.” Water Resour., 38(8), 20-1–20-5.
Information & Authors
Information
Published In
Copyright
© 2017 American Society of Civil Engineers.
History
Received: Apr 9, 2016
Accepted: Nov 4, 2016
Published online: Feb 6, 2017
Published in print: Jul 1, 2017
Discussion open until: Jul 6, 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.