Semianalytical Time Integration for Transient Groundwater Flow in Confined Aquifers
Publication: Journal of Hydrologic Engineering
Volume 12, Issue 1
Abstract
A semianalytical time integration method is proposed for the ordinary differential equations produced by the spatial discretization of the transient groundwater flow equation. Instead of approximating the time derivative by finite difference, the proposed method approximates the exact solution of the ordinary differential equations. The method is unconditionally stable; the accuracy depends only on the approximation accuracy of the stress; for piecewise constant or linear stress (e.g., pumping) in time, the solution can be exact; and the time step size can be as long as a stress period. The tradeoff is the computational cost, which can be reduced by using larger and less variable time step sizes. Two examples are given to show the performance of the semianalytical time integration method.
Get full access to this article
View all available purchase options and get full access to this article.
Acknowledgments
This material is based upon work supported by the National Science Foundation under Grant No. NSFCMS-0093752 (CAREER program). This research is also partially supported by CenSSIS, the Center for Subsurface Sensing and Imaging systems, under the Engineering Research Centers Program of the National Science Foundation (Grant No. NSFEEC-9986821). Any opinions, findings, and conclusions or recommendations expressed in this material are those of the writers and do not necessarily reflect the views of the National Science Foundation.
References
Anderson, M. P., and Woessner, W. W. (1992). Applied groundwater modeling-simulation of flow and advective transport, Academic, San Diego.
Andreu, J., and Sahuquillo, A. (1987). “Efficient aquifer simulation in complex-systems.” J. Water Resour. Plann. Manage., 113(1), 110–129.
Bear, J., and Verruijt, A. (1987). Modeling groundwater flow and pollution, Reidel, Dordrecht, Holland.
Brogan, W. L. (1991). Modern control theory, 3rd Ed., Prentice–Hall, Englewood Cliffs, N.J.
Cook, R., Malkus, D., and Plesha, M. (1989). Concepts and applications of finite element analysis, Wiley, New York.
Crawford, C. R. (1973). “Reduction of a band-symmetric generalized eigenvalue problem.” Commun. ACM, 16, 41–44.
Dunbar, W. S., and Woodbury, A. D. (1989). “Application of the Lanczos algorithm to the solution of the groundwater flow equation.” Water Resour. Res., 25(3), 551–558.
Gambolati, G. (1993). “On time integration of groundwater flow equations by spectral methods.” Water Resour. Res., 29(4), 1257–1267.
George, A., and Liu, J. W. (1981). Computer solution of large sparse positive definite systems, Prentice–Hall, Englewood Cliffs, N.J.
Ginn, T. R., Cushman, J. H., and Houck, M. H. (1990). “A continuous—time inverse operator for groundwater and contaminant transport modeling I, deterministic case.” Water Resour. Res., 26, 241–252.
Golub, G. H., and Van Loan, C. F. (1996). Matrix computations, 3rd Ed., Johns Hopkins University Press, Baltimore.
Harbaugh, A. W., Banta, E. R., Hill, M. C., and McDonald, M. G. (2000). “MODFLOW-2000, The U.S. Geological Survey modular ground-water model User guide to modularization concepts and the ground-water flow process.” Open-File Rep. No. 00-92, U.S. Geological Survey, Center, Miss.
Hill, M. C. (1990). “Preconditioned conjugate-gradient 2 (PCG2), a computer program for solving ground-water flow equations.” Water-Resources Investigations Rep. No. 90-4048, U.S. Geological Survey, Center, Miss.
Kuiper, L. K. (1973). “Analytic solution of spatially discretized groundwater flow equations.” Water Resour. Res., 9(4), 1094–1097.
Moler, C., and Van Loan, C. (2003). “Nineteen dubious ways to compute the exponential of a matrix, twenty-five years later.” SIAM Rev., 45(1), 3–49.
Moulton, D. J., Knapek, S., and Dendy, J. E. (1999). “Multilevel upscaling in heterogeneous porous media.” Rep. No. LA-UR 99-4754 ⟨http://cnls.lanl.gov/Highlights/1999-01⟩ (May 6, 2005).
Pinder, G. F., and Frind, E. O. (1972). “Application of Galerkin’s procedure to aquifer analysis.” Water Resour. Res., 8(1), 108–120.
Reddy, J. N. (1993). An introduction to the finite element method, 2nd Ed., McGraw–Hill, New York.
Sahuquillo, A. (1983). “An eigenvalue numerical technique for solving unsteady linear groundwater models continuous in time.” Water Resour. Res., 19(1), 87–93.
Sidje, R. B. (1998). “Expokit: A software package for computing matrix exponentials.” ACM Trans. Math. Softw., 24(1), 130–156.
Strang, G. (1986). Introduction to applied mathematics, Wellesley-Cambridge Press, Wellesley, Mass.
Sudicky, E. A. (1989). “The Laplace transform Galerkin technique: A time-continuous finite element theory and application to mass transport in groundwater.” Water Resour. Res., 25(8), 1833–1846.
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.
Tsai, W.-F., Shen, C.-Y., Fu, H.-H., and Kou, C.-C. (1999). “Study of parallel computation for ground water solute transport.” J. Hydraul. Eng., 4(1), 49–56.
Umari, A. M., and Gorelick, S. M. (1986). “Evaluation of the matrix exponential for use in ground-water-flow and solute-transport simulations: Theoretical framework.” USGS Water-Resources Investigations Rep. No. 86-4096, U.S. Geological Survey, Center, Miss.
Wang, H. F., and Anderson, M. P. (1982). Introduction to Groundwater modeling-finite difference and finite element methods, W. H. Freeman and Company, San Francisco.
Zhong, W. X. (2004). “On precise integration method.” J. Comput. Appl. Math., 163, 59–78.
Information & Authors
Information
Published In
Copyright
© 2007 ASCE.
History
Received: May 20, 2005
Accepted: Jun 7, 2006
Published online: Jan 1, 2007
Published in print: Jan 2007
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.