Numerical Study of Finite Element Method Based Solutions for Propagation of Wetting Fronts in Unsaturated Soil
Publication: Journal of Geotechnical and Geoenvironmental Engineering
Volume 130, Issue 3
Abstract
The accurate prediction of the propagation of a wetting front in an unsaturated soil subjected to surficial infiltration is of practical importance to many geotechnical and geoenvironmental problems. The finite element method is the most common solution technique as the hydraulic soil properties are highly nonlinear. Two important issues are often found to create difficulties in such analyses. First, numerical oscillations are usually observed in the calculated pore pressures at the wetting front. Second, when a reasonable mesh size and time step are used, the elevation of the wetting front may be seriously overpredicted. This paper is focused on the second issue. The under-relaxation (UR) technique used in the iterative process within each time step is found to have a serious impact on rate of convergence with refinement in mesh size and time step. Two different techniques are typically used; the first evaluates the hydraulic conductivity using an average of heads calculated from the preceding time node and the most recent iteration of the current time node (UR1), and the second evaluates the hydraulic conductivity using the average of heads calculated from the two most recent iterations of the current time nodes (UR2). The study shows that UR1, which is adopted in programs such as SEEP/W, ensures that the solution converges rapidly to a stable solution within a time step, but may converge to the wrong wetting front at a given elapsed time unless a sufficiently refined mesh is used. UR2 converges much more slowly within a time step, but the error in the wetting front is smaller than that generated by UR1.
Get full access to this article
View all available purchase options and get full access to this article.
References
Celia, M. A., Bouloutas, E., and Zarba, R. L.(1990). “A general mass-conservative numerical solution for the unsaturated flow equation.” Water Resour. Res., 26(7), 1483–1496.
Chong, P. C. (2000). “Convergence study of transient seepage in unsaturated soils.” M. Eng. thesis, The National Univ. of Singapore.
Cooley, R. L.(1983). “Some new procedures for numerical solution of variably saturated flow problems.” Water Resour. Res., 19, 1271–1285.
Fredlund, D. G., and Rahardjo, H. (1993). Soil mechanics for unsaturated soils, Wiley, New York.
Geo-Slope. (2000). User’s guide for SEEP/W and SLOPE/W, Version 4.23, Geo-Slope International, Calgary, Canada.
Karthikeyan, M., Tan, T. S., and Phoon, K. K.(2001). “Numerical oscillation in seepage analysis of unsaturated soils and its effect on slope stability.” Can. Geotech. J. 38(3), 639–651.
Leij, F. J., Russell, W. B., and Lesch, S. M.(1996). “Closed-form expressions for water retention and conductivity data.” Ground Water, 35(5), 848–858.
Mualem, Y.(1976). “A new model for predicting the hydraulic conductivity of unsaturated porous media.” Water Resour. Res., 12, 513–522.
Neuman, S. P.(1973). “Saturated-unsaturated seepage by finite elements.” J. Hydraul. Div., Am. Soc. Civ. Eng., 99(12), 2233–2250.
Paniconi, C., and Putti, M.(1994). “A comparison of Picard and Newton iteration in the numerical solution of multidimensional variably saturated flow problems.” Water Resour. Res., 30(12), 3357–3374.
Press, W. H., Flannery, B. P., Teukolsky, S. A., and Vetterling, W. T. (1986). Numerical recipes, Cambridge Univ. Press, Cambridge, U.K.
Rahardjo, H., and Leong, E. C. (1997). “Soil water characteristic curves and flux boundary problems.” Unsaturated Soil Engineering (GSP 68), ASCE, Reston, Va., 88–112.
Shouse, P. J., Russell, W. B., Burden, D. S., Selim, H. M., Sisson, J. B., and van Genuchten, M. Th.(1995). “Spatial variability of soil water retention functions in a silt loam soil.” Soil Sci., 159(1), 1–12.
Sweeney, D. J., and Robertson, P. K.(1979). “A fundamental approach to slope stability problems in Hong Kong.” Hong Kong Inst. Eng. J., 7(10), 35–44.
Thomas, H. R., and Zhou, Z.(1997). “Minimum time-step size for diffusion problem in FEM analysis.” Int. J. Numer. Methods Eng., 40, 3865–3880.
van Genuchten, M. Th.(1980). “A closed-form equation for predicting hydraulic conductivity of unsaturated soils.” Soil Sci. Soc. Am. J., 44(5), 892–898.
Warrick, A. W., Lomen, D. O., and Yates, S. R.(1985). “A generalized solution to infiltration.” Soil Sci. Soc. Am. J., 49, 34–38.
Information & Authors
Information
Published In
Copyright
Copyright © 2004 American Society of Civil Engineers.
History
Received: Oct 23, 2000
Accepted: May 16, 2002
Published online: Feb 19, 2004
Published in print: Mar 2004
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.