NUMERICAL EVALUATION OF UNSATURATED FLOW IN POROUS MEDIA USING ITERATIVE TECHNIQUES FOR VARIOUS SOILS
DOI:
https://doi.org/10.11113/mjce.v27.15920Keywords:
Saturation, numerical solution, finite element, iterative techniques, Picard and Newton schemesAbstract
Unrelenting difficulties that arise in forming numerical solutions of water flow in variably saturated problems include controlling the stability of the nonlinear equation solver. Besides, the accuracy of solutions, also of concern is the required computational effort, especially when nonlinear soil hydraulic properties and dry initial conditions are involved. A general numerical algorithm in the context of finite element formulation and Picard and Newton iteration schemes are described to solve Richards’ equation, in which a mass-conservative, head based formulation is proposed to approximate the governing equation, and mass-lumping techniques is employed to improve solution convergence and stability behavior. A difficult one-dimensional test case presented emphasizes different aspects of the performance of the iterative methods and the different issues that can influence their convergence and efficiency. The result is validated with extensive average values of soil hydraulic parameters for 12 soil textural groups.References
Carsel, R.F., Parrish, R.S. (1988). Developing joint probability distributions of soil water retention characteristic. Water Resour. Res. 24(5): 755-76.
Celia, M. A. and Binning, P. (1992). A mass conservative numerical solution for two-phase flow in porous media with application to unsaturated flow. Water Resour. Res., 28(10):
–28.
Cooley, R.L. (1983). Some new procedures for numerical solution of variably saturated Flow problems. Water Resour. Res., 19: 1271-1285.
Huyakorn, P. S.; Thomas, S. D. and Thompson, B. M. (1984). Techniques for making finite elements competitive in modeling flow in variably saturated media. Water Resources Res., 20: 1099-1115.
Huyakorn, P. S.; Springer, E. P., Guvanasen, V. and Wadsworth, T. D. (1986). A three dimensional finite element model for simulating water flow in variably saturated porous media. Water Resources Res., 22: 1790-1808.
Neuman, S. P. (1973). Saturated-unsaturated seepage by finite elements. J. Hydraul. Div. ASCE, 99(HY12): 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 Resources Res., 30: 3357–3374.
Paniconi, C., Aldama, A. A. and Wood, E.F. (1991). Numerical evaluation of iterative and noniterative methods for the solution of the nonlinear Richards equation. Water Resources Res., 27: 1147-1163.
Pan, L., Wierenga, P.J. (1997). Improving numerical modeling of two-dimensional water flow in variably saturated, heterogeneous porous media. Soil. Sci. Soc. Am. J., 61: 335-346.
van Genuchten, M. T., and Nielsen, D. R. (1985). On describing and predicting the hydraulic properties of unsaturated soils. Ann. Geophys., 3(5): 615-628.
Vogel, T., van Genuchten, M.T., Cislerova, M. (2011). Effect of the shape of the soil hydraulic functions near saturation on variably-saturated flow predictions. Advances in Water Resources, 24:133-144.
Williams, G.A., Miller, C.T., Kelley, C.T. (2000). Transformation approaches for simulating flow in variably saturated porous media. Water Resour. Res. 22: 831-840.
