Evaluation of Hopscotch Method for Transient Ground-Water Flow
Publication: Journal of Hydraulic Engineering
Volume 126, Issue 8
Abstract
The hopscotch finite-difference technique is shown to be a fast and accurate way to simulate transient, saturated, ground-water flow in relatively typical but heterogeneous 2D and 3D domains. The odd-even hopscotch (OEH) and line hopscotch methods are reviewed, and their implementation for saturated ground-water flow is presented. The OEH scheme, which is a second-order accurate explicit process, is efficient, requiring only six floating point operations per mesh node and time step, and is unconditionally stable (for saturated ground-water flow). Numerical experiments on typical 2D meshes (2,500 nodes) with synthetic, randomly heterogeneous hydraulic conductivity, suggest that the OEH process is approximately 1.5 times faster than the alternating direction implicit method and 3–4 times faster than the Crank-Nicolson implicit method using preconditioned conjugate gradient iteration. Similar experiments on medium-sized 3D meshes (87,500 nodes) suggest that the OEH process is between 7 and 10 times faster than the Crank-Nicolson preconditioned conjugate gradient method. Although the numerical results presented illustrate only typical test problem performance, they nevertheless clearly indicate promise for using OEH to simulate transient ground-water flow in 2D and, especially, 3D heterogeneous domains requiring fine spatial meshes.
Get full access to this article
View all available purchase options and get full access to this article.
References
1.
Ababou, R., McLaughlin, D., Gelhar, L., and Tompson, A. F. B. (1989). “Numerical simulation of three-dimensional saturated flow in randomly heterogeneous porous media.” Transport in Porous Media, 4, 549–565.
2.
Anderson, E., et al. ( 1992). LAPACK user's guide, Society for Industrial and Applied Mathematics, Philadelphia.
3.
Bear, J. (1979). Hydraulics of groundwater, McGraw-Hill, New York.
4.
Carrera, J., and Neuman, S. P. (1986). “Estimation of aquifer parameters under transient and steady-state conditions: 2. Uniqueness, stability, and solution algorithms.” Water Resour. Res., 22(2).
5.
Davidson, B., Vichnevetsky, R., and Wang, H. T. (1978). “Numerical techniques for estimating best-distributed Manning's roughness coefficients for open estuarial river systems.” Water Resour. Res., 14(3), 777–789.
6.
DeGoede, E. F., and Ten Thije Boonkkamp, J. H. M. (1990). “Vectorization of the odd-even hopscotch scheme and the alternating direction implicit scheme for the two-dimensional Burgers equations.” J. Scientific and Statistical Computing, 11(2), 354–367.
7.
Dongarra, J., et al. (1991). Solving linear systems on vector and shared memory computers, Society for Industrial and Applied Mathematics, Philadelphia.
8.
Doughty, C., Long, J. C. S., Hestir, K., and Benson, S. M. (1994). “Hydrologic characterization of heterogeneous geological media with an inverse method based on iterated function systems.” Water Resour. Res., 30(6), 1721–1745.
9.
Gane, C. R., and Gourlay, A. R. (1977). “Block hopscotch procedures for second order parabolic differential equations.” J. Inst. for Mathematics and its Applications, 19, 205–216.
10.
Gelhar, L. (1986). “Stochastic subsurface hydrology from theory to application.” Water Resour. Res., 22(9).
11.
Golub, G., and Van Loan, C. (1983). Matrix computations, Johns Hopkins University Press, Baltimore.
12.
Gordon, P. (1965). “Nonsymmetric difference equations.” J. Soc. for Industrial and Appl. Mathematics, 13.
13.
Gourlay, A. R. (1970). “Hopscotch: A fast second-order partial differential equation solver.” J. Inst. for Mathematics and its Applications, 6, 375–390.
14.
Gourlay, A. R., and McGuire, G. R. (1971). “General hopscotch algorithm for the numerical solution of partial differential equations.” J. Inst. for Mathematics and its Applications, 7, 216–227.
15.
Gourlay, A. R., and McKee, S. (1977). “The construction of hopscotch methods for parabolic and elliptic equations in two space dimensions with a mixed derivative.” J. Computational and Appl. Mathematics, 3(3), 201–206.
16.
Gourlay, A. R., and McKee, S. (1979). “Hopscotch methods for elliptic partial differential equations.” J. Computational and Appl. Mathematics, 5(2), 103–110.
17.
Henrici, P. (1962). Discrete variable methods in ordinary differential equations, Wiley, New York.
18.
Hill, M. (1990). “Solving groundwater flow problems by conjugate-gradient methods and the strongly implicit procedure.” Water Resour. Res., 26(9), 1961–1969.
19.
Huyakorn, P. S., and Pinder, G. F. (1983). Computational methods in subsurface flow, Academic, San Diego.
20.
Kershaw, D. S. (1978). “The incomplete Cholesky-conjugate gradient method for the iterative solution of systems of linear equations.” J. Computational Phys., 26, 43–65.
21.
Lapidus, L., and Pinder, G. F. (1982). Numerical solution of partial differential equations in science and engineering, Wiley, New York.
22.
McCormick, S. F., ed. (1987). Multigrid methods, Society for Industrial and Applied Mathematics, Philadelphia.
23.
McDonald, M. G., and Harbaugh, A. W. (1984). A modular three-dimensional finite-difference groundwater flow model, Scientific Publications, Washington, D.C.
24.
Meyer, P. D., Valocchi, A. J., Ashby, S. F., and Saylor, P. E. (1989). “A numerical investigation of the conjugate gradient method as applied to three-dimensional groundwater flow problems in randomly heterogeneous porous media.” Water Resour. Res., 25(6), 1440–1446.
25.
Oppe, T. ( 1992). STENCIL user's guide: A package for solving large sparse linear systems by various iterative methods, Supercomputer Computations Research Institute, Florida State University, Tallahassee, Fla.
26.
Peaceman, D. W., and Rachford, H. H., Jr. (1955). “The numerical solution of parabolic and elliptic differential equations.” J. Soc. for Industrial and Appl. Mathematics, 3(1), 28–41.
27.
Ségol, G. (1992). “The hopscotch algorithm for three-dimensional simulation.”J. Hydr. Engrg., ASCE, 118(3), 385–406.
28.
Ten Thije Boonkkamp, J. H. M. (1988). “The odd-even hopscotch pressure-correction scheme for the incompressible Navier-Stokes equations.” J. Scientific and Statistical Computing, 9(2).
29.
Ten Thije Boonkkamp, J. H. M., and Verwer, J. G. (1987). “On the odd-even hopscotch scheme for the numerical integration of time-dependent partial differential equations.” Appl. Numer. Mathematics, 3, 183–193.
30.
Theis, C. V. (1935). “The relation between the lowering of the piezometric surface and the rate and duration of discharge of a well using groundwater storage.” Trans. Am. Geophys. Union, 2, 519–524.
31.
Tompson, A. F. B., Ababou, R., and Gelhar, L. W. (1989). “Implementation of the three-dimensional turning bands random field generator.” Water Resour. Res., 25(10).
32.
Trescott, P. C., and Larson, S. P. (1977). “Comparison of iterative methods of solving two-dimensional groundwater flow equations.” Water Resour. Res., 13(1), 125–136.
Information & Authors
Information
Published In
Journal of Hydraulic Engineering
Volume 126 • Issue 8 • August 2000
Pages: 615 - 626
History
Received: Feb 27, 1995
Published online: Aug 1, 2000
Published in print: Aug 2000
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.