Unsaturated Soil Seepage Analysis Using a Rational Transformation Method with Under-Relaxation
Publication: International Journal of Geomechanics
Volume 8, Issue 3
Abstract
The finite-element method provides a convenient and effective means for solving problems of seepage in unsaturated soils. However, convergence difficulties exist in numerical simulations of unsaturated flow analyses because of the high nonlinearity of the soil hydraulic properties. This technical note presents a combination approach consisting of a rational function transformation method and a common under-relaxation technique to solve the -based form of Richards equation. Numerical studies show that this combined method can use a larger time step and corresponding oscillation-free mesh size to produce acceptable results and also converge to a stable solution quickly in each time step.
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.
Diersch, H.-J. G., and Perrochet, P. (1999). “On the primary variable switching technique for simulating unsaturated-saturated flows.” Adv. Water Resour., 23, 271–301.
GEO-SLOPE International Ltd. (2004). Seepage modeling with SEEP/W: An engineering methodology, Calgary, Alb., Canada.
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, 285–294.
Ju, S. H., and Kung, K. J. S. (1997). “Mass types, element orders and solution schemes for the Richards equation.” Comput. Geosci., 23(2), 175–187.
Karthikeyan, M., Tan, T. S., and Phoon, K. K. (2001). “Numerical oscillation in seepage analysis of unsaturated soils.” Can. Geotech. J., 38, 639–651.
Kirkland, M. R., Hills, R. G., and Wierenga, P. J. (1992). “Algorithms for solving Richards’ equation for variably saturated soils.” Water Resour. Res., 28(14), 2049–2058.
McBride, D., Cross, M., Croft, N., Bennett, C., and Gebhardt, J. (2006). “Computational modeling of variably saturated flow in porous media with complex three-dimensional geometries.” Int. J. Numer. Methods Fluids, 50, 1085–1117.
Mualem, Y. (1976). “A new model for predicting the hydraulic conductivity of unsaturated porous media.” Water Resour. Res., 12, 513–522.
Murti, V., Valliappan, S., and Khalili-Naghadeh, N. (1981). “Time step constraints in finite-element analysis of the Poisson type equation.” Comput. Struct., 31, 269–273.
Pan, L., and Wierenga, P. J. (1995). “A transformed pressure head-based approach to solve Richards’ equation for variably saturated soils.” Water Resour. Res., 31(4), 925–931.
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.
Rank, E., Katz, C., and Werner, H. (1983). “On the importance of the discrete maximum principle in transient analysis using finite-element methods.” Int. J. Numer. Methods Eng., 19, 1771–1782.
Rathfelder, K., and Abriola, L. M. (1994). “Mass conservative numerical solutions of the head-based Richards equation.” Water Resour. Res., 30(9), 2579–2586.
Ross, P. J. (1990). “Efficient numerical methods for infiltration using Richards’ equation.” Water Resour. Res., 26(2), 279–290.
Tan, T. S., Phoon, K. K., and Chong, P. C. (2004). “Numerical study of finite-element method based solutions for propagation of wetting fronts in unsaturated soil.” J. Geotech. Geoenviron. Eng., 130(3), 254–263.
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. T. (1980). “A close-form equation for predicting the hydraulic conductivity of unsaturated soils.” Soil Sci., 44(5), 892–898.
Vermeer, P. A., and Verruijt, A. (1981). “An accuracy condition for consolidation by finite elements.” Int. J. Numer. Analyt. Meth. Geomech., 5, 1–14.
Warrick, A. W., Lomen, D. O., and Yates, S. R. (1985). “A generalized solution to infiltration.” Soil Sci. Soc. Am. J., 49, 34–38.
Williams, G. A., and Miller, C. T. (1999). “An evaluation of temporally adaptive transformation approaches for solving Richards’ equation.” Adv. Water Resour., 22(14), 831–840.
Williams, G. A., Miller, C. T., and Kelley, C. T. (2000). “Transformation approaches for simulating flow in variably saturated porous media.” Water Resour. Res., 36(4), 923–934.
Wu, Y. S., and Forsyth, P. A. (2001). “On the selection of primary variables in numerical formulation for modeling multiphase flow in porous media.” J. Contam. Hydrol., 48(3–4), 277–304.
Information & Authors
Information
Published In
Copyright
© 2008 ASCE.
History
Received: May 3, 2006
Accepted: Oct 18, 2007
Published online: May 1, 2008
Published in print: May 2008
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.