Mathematical Forms and Numerical Schemes for the Solution of Unsaturated Flow Equations
Publication: Journal of Irrigation and Drainage Engineering
Volume 138, Issue 1
Abstract
Prediction of water infiltration into the soil, fluid movement in the unsaturated soils, and groundwater recharge are important problems in different fields of science and engineering. Moreover, the transfer of the different pollutants (e.g., pesticides) from ground surface to groundwater occurs through the unsaturated zone. In this paper, a comprehensive evaluation of different finite difference schemes (e.g., fully implicit, Crank-Nicolson, and Runge-Kutta) is presented for the solution of head-based and mixed forms of the Richard’s equation. Two examples of water infiltration in very dry and relatively wet unsaturated soils are used for the evaluation of schemes. In addition, the effects of various approximations of moisture capacity function, convergence criteria, and time stepping methods on the performance of the schemes are investigated, and the results showed their significant influences on mass balance, number of iterations, and convergence condition of the numerical schemes. The results of numerical simulations showed that the generally mixed form has better performance than the head-based form. In addition, the Crank-Nicolson scheme showed better results than the Runge-Kutta scheme, but both have convergence problems. In general, the comparison of numerical methods showed that the fully implicit scheme has the best performance among various finite difference schemes and can be selected as a reliable scheme with acceptable solutions for different infiltration problems in unsaturated soils.
Get full access to this article
View all available purchase options and get full access to this article.
Acknowledgments
This work was supported by a grant of Shahrood University of Technology, Shahrood, Iran (Project number: UNSPECIFIED21014).
References
Allen, M. B., and Murphy, C. L. (1985). “A finite element collocation method for variably saturated flows in porous media.” Numer. Methods Partial Differential Equations, 1(3), 229–239.
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., Bouloutas, B. T., and Zarba, R. L. (1990). “A general mass conservative solution for the unsaturated flow equation.” Water Resour. Res., 26(7), 1483–1496.
Cooley, R. L. (1983). “Some new procedures for numerical solution of variably saturated flow problems.” Water Resour. Res., 19(5), 1271–1285.
El-Kadi, A. I., and Ling, G. (1993). “The Courant and Peclet number criteria for the numerical solution of Richards equation.” Water Resour. Res., 29(10), 3485–3494.
Haverkamp, R., Vauclin, M., Touma, J., Wierenga, P. J., and Vachaud, G. (1977). “A comparison of numerical simulation models for one-dimensional infiltration.” Soil Sci. Soc. Am. J., 41(2), 285–294.
Hills, R. G., Porro, I., Hudson, D. B., and Wierenga, P. J. (1989). “Modeling one-dimensional infiltration into very dry soils I. Model development and evaluation.” Water Resour. Res., 25(6), 1259–1269.
Huang, K., Mohanty, B. P., and van Genuchten, M. Th. (1996). “A new convergence criterion for the modified Picard iteration method to solve the variably saturated flow equation.” J. Hydrol., 178(1-4), 69–91.
Huyakorn, P. S., and Wadsworth, T. D. (1985). FLAMINCO: A three dimensional finite element code for analyzing water flow and solute transport in saturated-unsaturated porous media, Geo. Trans. Inc., Herndon, VA.
Kaluarachchi, J. J., and Parker, J. C. (1989). “An efficient finite element method for modeling multiphase flow.” Water Resour. Res., 25(1), 43–54.
Magnuson, S. O., Baca, R. G., and Sondrup, A. J. (1990). “Independent verification and benchmark testing of the PORFLO-3 computer code.” Version 1.0, Rep. EGG-BG-9175, Idaho National Laboratory, Idaho Falls, ID.
Milly, P. C. D. (1985). “A mass-conservative procedure for time stepping in models of unsaturated flow.” Adv. Water Resour., 8(1), 32–36.
Rathfelder, K., and Abriola, L. M. (1994). “Mass conservative numerical solutions of the head based Richards equation.” Water Resour. Res., 30(9), 2579–2586.
Ray, C., and Mohanty, B. P. (1992). “Some numerical investigations of the Richards equation.” ASAE Paper 92-2586, ASAE, St. Joseph, MI.
Richards, L. A. (1931). “Capillary conduction of liquids through porous mediums.” J. Appl. Physics, 1(5), 318–333.
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.
van Genuchten, M. Th. (1980). “A closed form equation for predicting the hydraulic conductivity of saturated soils.” Soil Sci. Soc. Am. J., 44(5), 892–898.
van Genuchten, M. Th. (1982). “A comparison of numerical solutions of one-dimensional unsaturated-saturated flow and mass transport equations.” Adv. Water Resour., 5(1), 47–55.
Zarba, R. L. (1988). “A numerical investigation of unsaturated flow.” MS thesis, Dept. of Civil Eng., Mass. Inst. of Technol., Cambridge, MA.
Information & Authors
Information
Published In
Journal of Irrigation and Drainage Engineering
Volume 138 • Issue 1 • January 2012
Pages: 63 - 72
Copyright
© 2012 American Society of Civil Engineers.
History
Received: Jun 13, 2010
Accepted: May 5, 2011
Published online: May 6, 2011
Published in print: Jan 1, 2012
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.