Finite Analytic Numerical Solutions for Unsaturated Flow with Irregular Boundaries
Publication: Journal of Hydraulic Engineering
Volume 119, Issue 11
Abstract
This study illustrates the development of the finite‐analytic (FA) numerical solutions associated with an optimal (variable) time‐weighting factor for water flow in unsaturated porous media with irregular boundaries. This method entails the nine‐point and five‐point elements such that it can simulate quite well the irregular flow domain without using boundary‐fitted transformation. For an element near the irregular boundaries, a five‐point FA algebraic solution is derived; for a regular element, the nine‐point FA algebraic equation is used. In execution of numerical computation, the irregular five‐point boundary element is recasted into the nine‐point element, allowing a regular structured grid system to be used to solve for the water flow in an irregular domain. Two problems are solved in this study. The first is a transient, two‐dimensional flow of water under a strip source infiltration of constant flux with regular boundaries. Uniform clay loam, uniform coarse sand, and layered system of soil are used in the numerical simulation. It is shown that the finite‐analytic numerical solutions conform to accuracy and numerical stability very well. Mass conservation for the FA method also is illustrated. The second problem simulates unsaturated flow of leachate from landfill with irregular boundaries. The water pressure head and content in the unsaturated zone are predicted as a function of time and space.
Get full access to this article
View all available purchase options and get full access to this article.
References
1.
Bear, J. (1979). Hydraulics of groundwater water. McGraw‐Hill, New York, N.Y., 213–224.
2.
Bravo, R. H. (1991). “Development of the three‐dimensional finite analytic method for simulation of fluid flow and conjugate heat transfer,” PhD thesis, Dept. of Mechanical Engineering, the University of Iowa, Iowa City, Iowa.
3.
Celia, M. A., Bouloutas, E. T., and Zarba, R. L. (1990). “A general mass‐conservative numerical solution for the unsaturated flow equation.” Water Resour. Res., 26(7), 1483–1496.
4.
Celia, M. A., Ahuja, L. R., and Pinder, G. F. (1987). “Orthogonal collocation and alternating‐direction procedures for unsaturated flow problems.” Adv. in Water Resour., 10, 178–187.
5.
Chen, C. J., and Chen, H. C. (1984a). “Finite analytic numerical method for unsteady two‐dimensional Navier‐Stokes equation.” J. Comput. Phys., 53(2), 209–226.
6.
Chen, H. C., and Chen, C. J. (1984b). “Development of finite analytic numerical method for unsteady three‐dimensional Navier‐Stokes equations.” ASME Symp. on Computation of Internal Flow, ASME FED, 14, 159–165.
7.
Chen, C. J., and Choi, S. K. (1990). “The finite analytic method and its application—Laminar and turbulent flows past two dimensional and axisymmetric bodies.” IIHR Rep. No. 344‐I, Iowa Institute of Hydraulic Research, University of Iowa, Iowa City, Iowa.
8.
Gambolati, G., Sartoretto, F., and Uliana, F. (1986). “A conjugate gradient finite element model of flow for large multiaquifer system.” Water Resour. Res., 22(7), 1003–1015.
9.
Hwang, J. C., Chen, C. J., Sheikhoslami, M., and Panigrahi, B. K. (1985). “Finite analytic numerical solution for two‐dimensional groundwater solute transport.” Water Resour. Res., 21(90), 1354–1360.
10.
Li, S. G., Ruan, F., and McLaughlin, D. (1992). “A space‐time accurate method for solving solute transport problems.” Water Resour. Res., 28(9), 2297–2306.
11.
Milly, P. C. D. (1985). “A mass‐conservative procedure for time‐stepping in models of unsaturated flow.” Adv. in Water Resour., 8, 32–36.
12.
Obasih, K. M., and Chen, C. J. (1987). “Prediction of laminar and turbulent flows past single and twin hydrofoils.” IIHR Rep. No. 317, Iowa Institute of Hydraulic Research, Univ. of Iowa, Iowa City, Iowa.
13.
Sudicky, E. A., and McLaren, R. G. (1992). “Laplace transform Galerkin technique for large‐scale simulation of mass transport in discretely fractured porous formations.” Water Resour. Res., 28(2), 499–514.
14.
Tien, H. C. (1993). “Finite analytic method for computational domain with irregular boundaries,” PhD thesis, Dept. of Applied Mathematics, the University of Iowa, Iowa City, Iowa.
15.
Tsai, W. F., and Chen, C. J. (1993). “Finite analytic numerical solutions for contaminant transport in two‐dimensional groundwater flow.” XXV IAHR Cong., Tokyo, Japan.
16.
Tsakiris, G. P., Soulis, J. V., and Bellos, C. V. (1991). “Two‐dimensional unsaturated flow in irregular shaped regions using a finite volume method.” Transport in Porous Media, 6, 1–12.
17.
Warrick, A. W., and Lomen, D. O. (1976). “Time‐dependent linearized infiltration: III. strip and disc sources.” Soil Sci. Soc. Am. J., 40, 639–643.
18.
Yeh, G. T., and Ward, D. S. (1980). “FEMWATER: a finite‐element model of water flow through saturated‐unsaturated porous media.” Rep. ORNL‐5567, Oak Ridge Nat. Lab., Oak Ridge, Tenn.
19.
Yeh, G. T. (1987). “FEMWATER: a finite‐element model of water flow through saturated‐unsaturated porous media.” Rep. ORNL‐5567/R1, 1st Rev., Oak Ridge Nat. Lab., Oak Ridge, Tenn.
20.
Zeng, X. J., and Li, W. (1987). “The stability and convergence of finite analytic method for unsteady two‐dimensional convective transport equations.” Turbulence Measurements and Flow Modeling, C. J. Chen, L. D. Chen, and F. M. Holly, eds., Hemisphere Publishing Corp., 427–433.
Information & Authors
Information
Published In
Journal of Hydraulic Engineering
Volume 119 • Issue 11 • November 1993
Pages: 1274 - 1298
Copyright
Copyright © 1993 American Society of Civil Engineers.
History
Received: Jul 23, 1992
Published online: Nov 1, 1993
Published in print: Nov 1993
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.