Efficient Approach to Solving Transient Unsaturated Flow Problems. II: Numerical Solutions
Publication: International Journal of Geomechanics
Volume 17, Issue 7
Abstract
Although the finite-difference method (FDM) has been commonly used to numerically solve Richard’s equation, numerical difficulties are often encountered, even for comparatively simple problems. To minimize convergence problems, comparatively small discretization and time steps have often been adopted to solve this highly nonlinear equation, resulting in significant computational costs. To overcome these difficulties, this paper presents an efficient approach to solving Richard’s equation that combines two numerical techniques: the FDM and the cubic interpolated pseudoparticle (CIP) method. The FDM is used to solving the diffusive flow component of Richard’s equation, the convergence of which can be controlled by adopting time steps corresponding to Neumann’s number under 0.5. In contrast, the CIP method is used to solve the advective flow component of the equation. The CIP method is found to be particularly suitable for facilitating convergence and eliminating the presence of spurious results when the Courant number is under 1.0. Analytical solutions for transient unsaturated flow problems, developed in a companion paper, allow comparison between the predictions obtained using the proposed numerical approach and the exact solutions. Use of the newly developed algorithm is found to be particularly accurate and stable for solving Richard’s equation, being clearly superior to the use of the traditional FDM. After validating the new numerical approach using the boundary conditions and hydraulic functions for which analytical solutions have been developed, the new numerical scheme was subsequently implemented to address more general unsaturated flow problems. In particular, the new numerical approach was extended to solve unsaturated flow problems involving complex soil hydraulic functions as well as different boundary conditions. Comparisons are presented to illustrate the accuracy of the new numerical approach even when extended to incorporate the use of complex hydraulic functions for which there are no analytical solutions. The efficient, validated numerical schemes presented in this paper are found to be well suited for solving complex unsaturated flow problems.
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 institutions: 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
Bear, J. (1979). Hydraulics of groundwater, McGraw-Hill, New York.
Bouloutas, E. T. (1989). “Improved numerical approximations for flow and transport in the unsaturated zone.” Ph.D. thesis, Dept. of Civil Engineering, Massachusetts Institute of Technology, Cambridge, MA.
Brooks, R. H., and Corey, A. T. (1964). Hydraulic properties of porous media, Colorado State Univ., Fort Collins, CO.
Cavalcante, A. L. B., and Zornberg, J. G. (2016). “Numerical schemes to solve advective contaminant transport problems with linear sorption and first order decay.” Electron. J. Geotech. Eng., 21(5), 2043–2060.
Cavalcante, A. L. B., and Zornberg, J. G. (2017). “Efficient approach to solving transient unsaturated flow problems. I: Analytical solutions.” Int. J. Geomech., 04017013.
Celia, M. A., Ahuja, L. R., and Pinder, G. F. (1987). “Orthogonal collocation and alternating-direction procedures for unsaturated flow problems.” Adv. Water Resour., 10(4), 178–187.
Celia, M. A., and Bouloutas, E. T. (1990). “A general mass-conservative numerical solution for the unsaturated flow equation.” Water Resour. Res., 26(7), 1483–1496.
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.
Hundsdorfer, W. H., and Verwer, D. B. (2003). “Numerical solution of time-dependent advection-diffusion-reaction equations.” Springer series in computational mathematics, Vol. 33, Springer, New York.
Mathematica [Computer software]. Wolfram, Champaign, IL.
Philip, J. R. (1969). “Theory of infiltration.” Advances in hydroscience, V. T. Chow, ed., Vol. 5, Academic Press, New York, 215–296.
Richards, L. A. (1931). “Capillary conduction of liquids through porous mediums.” J. Appl. Phys., 1(5), 318–333.
Smith, G. D. (1985). Numerical solution of partial differential equations: Finite difference methods, 3rd Ed., Oxford University Press, New York, 67–68.
Takewaki, H., Nishiguchi, A., and Yabe, T. (1985). “Cubic interpolated pseudo-particle method (CIP) for solving hyperbolic-type equations.” J. Comput. Phys., 61(2), 261–268.
Takewaki, H., and Yabe, T. (1987). “The cubic-interpolated pseudo particle (CIP) method: Application to nonlinear and multi-dimensional hyperbolic equations.” J. Comput. Phys., 70(2), 355–372.
Van Genuchten, M. T. (1980). “A closed form equation for predicting the hydraulic conductivity of unsaturated soils.” Soil Sci. Am. J., 44(5), 892–898.
Information & Authors
Information
Published In
Copyright
© 2017 American Society of Civil Engineers.
History
Received: Apr 18, 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.