Mathematical Formulation and Validation of a Mixed Finite Element–Finite Difference Model for Simulating Phreatic Surfaces
Publication: Journal of Hydraulic Engineering
Volume 131, Issue 12
Abstract
The phreatic surface in an unconfined aquifer exists as a movable interface between the saturated and unsaturated zones. The movement of the phreatic surface depends on recharge, hydraulic conductivity, porosity, and horizontal and vertical flows. The location of the phreatic surface helps define the variably saturated flow domain in the subsurface. The variably saturated flow process in the subsurface is described by a parabolic partial differential equation. In this equation, the hydraulic conductivity and soil moisture capacity are used as the subsurface characteristics. The location of the phreatic surface is governed by a first-order partial differential equation. The governing parabolic partial differential equation is solved using a variational finite element formulation. The first order phreatic surface equation is then solved by loosely coupling with the governing parabolic partial differential equation describing the variably saturated flow. In the present study, a two-dimensional space is used to investigate the movement of the phreatic surface in a variably saturated unconfined flow domain. Based on the time-varying solutions of hydraulic heads, the location of the phreatic surface is simulated in a finite two-dimensional space.
Get full access to this article
View all available purchase options and get full access to this article.
References
Ahmed, S. (2001). “Variational formulations for a parabolic partial differential equation describing variably saturated flow.” Thesis, The Univ. of Tennessee, Knoxville, Tenn.
Aral, M. M., and Tang, Y. (1988). “A new boundary element formulation for time-dependent confined and unconfined aquifer problems.” Water Resour. Res., 24(6), 831–842.
Baker, G. A., Bramble, J. H., and Thomée, V. (1976). “Single step Galerkin approximations for parabolic problems.” Math. Comput., 31(140), 818–847.
Bear, J. (1979). Hydraulics of groundwater, McGraw-Hill, New York.
Binning, P., and Celia, M. A. (1996). “A finite volume Eulerian–Lagrangian localized adjoint method for solution of the contaminant transport equations in two-dimensional multiphase flow systems.” Water Resour. Res., 32(1), 103–114.
Du, X., Yao, T., Stone, W. D., and Hendrickx, J. M. H. (2001). “Stability analysis of the unsaturated water flow equation (1). Mathematical derivation.” Water Resour. Res., 37(7), 1869–1874.
Golub, G. H., and Loan, C. F. V. (1996). Matrix computations, The John Hopkins Univ. Press, Baltimore.
Hildebrand, F. (1976). Advanced calculus for applications, Prentice Hall, Englewood Cliffs, N.J.
Huyakorn, P. S., Anderson, P. F., Mercer, J. W., and White, Jr., H. O. (1987). “Saltwater intrusion in aquifers: Development and testing of a three-dimensional finite element model.” Water Resour. Res., 23(2), 293–312.
Konikow, L. F., and Bredehoeft, J. D. (1978). “Computer model of two-dimensional solute transport and dispersion in ground water.” USGS TWRI.
Kuppusamy, T., Sheng, J., Parker, J. C., and Lenhard, R. J. (1987). “Finite-element analysis of multiphase immiscible flow through soils.” Water Resour. Res., 23(4), 625–631.
Lal, A. M. W. (2000). “Numerical errors in groundwater and overland flow models.” Water Resour. Res., 36(5), 1237–1247.
Liggett, J., and Liu, P.L.-F. (1983). “The boundary integral equation method for porous media flow.” George Allen & Unwin.
Luskin, M., and Rannacher, R. (1981). “On the smoothing property of the Galerkin method for parabolic equations.” SIAM (Soc. Ind. Appl. Math.) J. Numer. Anal., 19(1), 93–113.
McDonnald, M. G., and Harbaugh, A. W. (1988). “A modular three-dimensional finite-difference ground-water flow model.” USGS TWRI Chapter 6-A1.
Mosé, R., Siegel, P., Ackerer, P., and Chavent, G. (1994). “Application of the mixed hybrid finite element approximation in a groundwater flow model: Luxury or necessity?” Water Resour. Res., 30(11), 3001–3012.
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.
Pollock, P. W. (1989). “Documentation of computer programs to compute and display pathlines using results from the USGS modular three-dimensional finite-difference ground-water flow model.” USGS Open File Rep. No. 89-391.
Srivastava, R., and Yeh, T.-C.J. (1992). “A three-dimensional numerical model for water flow and transport of chemically reactive solute through porous media under variably saturated conditions.” Adv. Water Resour., 15, 275–287.
Sun, N.-Z. (1999). “A finite cell method for simulating the mass transport process in porous medium.” Water Resour. Res., 35(12), 3649–3662.
Therrien, R. (1992). “Three-dimensional analysis of variably-saturated flow and transport in discretely-fractured porous media.” PhD thesis, Univ. of Waterloo, Ontario, Canada.
Thomée, V. (1980). “Negative norm estimates and superconvergence in Galerkin methods for parabolic problems.” Math. Comput., 34(149), 93–113.
van Genuchten, M. T. (1980). “A closed-form equation for predicting the hydraulic conductivity of unsaturated soils.” Soil Sci. Soc. Am. J., 44, 892–898.
Yeh, G. T. (1987). FEMWATER: A finite element model of water flow through saturated-unsaturated porous media, Oak Ridge National Laboratory, Oak Ridge, Tenn.
Yeh, G. T., and Ward, D. S. (1981). FEMWASTE: A finite-element model of waste transport through porous saturated-unsaturated media, Oak Ridge National Laboratory, Oak Ridge, Tenn.
Yeh, J. T.-C., Srivastava, R., Guzman, A, and Harter, T. (1993). “A numerical model for water flow and chemical transport in variably saturated porous media.” Ground Water, 31(4), 634–644.
Zachmanoglou, E. C., and Thoe, D. W. (1986). Introduction to partial differential equations with applications, Dover, New York.
Zheng, C. (1990). MT3D: A modular three-dimensional transport model for simulation of advection, dispersion, and chemical reactions of contaminants in ground-water systems, U.S. EPA, R. S. Kerr Environmental Research Laboratory, Ada, Okla.
Information & Authors
Information
Published In
Copyright
© 2005 ASCE.
History
Received: Apr 25, 2003
Accepted: Jan 18, 2005
Published online: Dec 1, 2005
Published in print: Dec 2005
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.