Automatic Time Stepping with Global Error Control for Groundwater Flow Models
Publication: Journal of Hydrologic Engineering
Volume 13, Issue 9
Abstract
Automatic time stepping with global error control is proposed for the time integration of the diffusion equation to simulate groundwater flow in confined aquifers. The scheme is based on an a posteriori error estimate for the discontinuous Galerkin finite-element methods. A stability factor is involved in the error estimate and it is used to adapt the time step and control the global temporal error for the backward Euler method. The stability factor can be estimated by solving a dual problem. The stability factor is not sensitive to the accuracy of the dual solution and the overhead computational cost can be minimized by solving the dual problem using large time steps. Numerical experiments are conducted to show the application and the performance of the automatic time stepping scheme. Implementation of the scheme can lead to improvements in accuracy and efficiency of groundwater flow models.
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) and by the U.S. Department of Energy’s Office of Science Biological and Environmental Research, Environmental Remediation Sciences Program (ERSP). 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. Oak Ridge National Laboratory is managed by UT-Battelle, LLC, for the U.S. Department of Energy under Contract No. DOEDE-AC05-00OR22725.DOE
References
Bellin, A., and Rubin, Y. (1996). “HYDRO_GEN: A spatially distributed random field generator for correlated properties.” Stochastic Hydrol. Hydraul., 10, 253–278.
Cao, Y., and Petzold, L. (2004). “A posteriori error estimation and global error control for ordinary differential equations by the adjoint method.” SIAM J. Sci. Comput. (USA), 26(2), 359–374.
Diersch, H.-J. G. (2005). FEFLOW reference manual, WASY GmbH, Berlin.
Eriksson, K., Estep, D., Hansbo, P., and Johnson, C. (1996). Computational differential equations, Cambridge University Press, Cambridge, Mass.
Eriksson, K., Johnson, C., and Logg, A. (2004). “Adaptive computational methods for parabolic problems.” Encyclopedia of computational mechanics, E. Stein, R. de Borst, and T. J. R. Hughes, eds., Wiley, New York.
Givoli, D., and Henigsberg, I. (1993). “A simple time-step control scheme.” Commun. Numer. Methods Eng., 9, 873–881.
Harbaugh, A. W., Banta, E. R., Hill, M. C., and McDonald, M. G. (2000). “MODFLOW-2000, the U.S.G.S. modular ground-water model–user guide to modularization concepts and the ground-water flow process.” Open-File Report No. 00-92, U.S.G.S.
Jannelli, A., and Fazio, R. (2006). “Adaptive stiff solvers at low accuracy and complexity.” J. Comput. Appl. Math., 191, 246–258.
Kavetski, D., Binning, P., and Sloan, S. W. (2001). “Adaptive timestepping and error control in a mass conservative numerical solution of the mixed form of Richards equation.” Adv. Water Resour., 24, 595–605.
Krahn, J. (2004). Seepage modeling with SEEP/W, 1st Ed., Geo-Slope International Ltd., Canada.
Lin, H. J., et al. (1997). “FEMWATER: A three-dimensional finite-element computer model for simulating density-dependent flow and transport in variably saturated media.” Technical Rep. No. CHL-97-12, U.S. Army Engineer Research and Development Center.
Logg, A. (2003). “Multi-adaptive Galerkin methods for ODEs. II: Implementation and applications.” SIAM J. Sci. Comput. (USA), 25(4), 1119–1141.
Miller, C. T., Abhishek, C., and Farthing, M. W. (2006). “A spatially and temporally adaptive solution of Richards’ equation.” Adv. Water Resour., 29, 525–545.
Soderlind, G., and Wang, L. (2006). “Evaluating numerical ODE/DAE methods, algorithms and software.” J. Comput. Appl. Math., 185, 244–260.
Tang, G. (2006). “Time integration for groundwater flow and solute transport modeling.” Ph.D. thesis, Northeastern Univ., Boston.
Tang, G., and Alshawabkeh, A. N. (2007). “A semi-analytical time integration for numerical solution of groundwater flow in confined aquifers.” J. Hydrol. Eng., 12(1), 73–82.
Yeh, G. T., et al. (2004).“HYDROGEOCHEM 5.0: A three-dimensional model of coupled fluid flow, thermal transport, and hydrogeochemical transport through variable saturated conditions.” Rep. No. ORNL/TM-2004/107, Oak Ridge National Laboratory, Oak Ridge, Tenn.
Information & Authors
Information
Published In
Copyright
© 2008 ASCE.
History
Received: Oct 23, 2006
Accepted: Jun 4, 2007
Published online: Sep 1, 2008
Published in print: Sep 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.