Overlapping Control Volume Method for Solute Transport
Publication: Journal of Hydrologic Engineering
Volume 5, Issue 3
Abstract
In this paper, an overlapping control volume method is presented for the numerical solution of transient 2D solute transport problems in ground water. The method is applicable for nonorthogonal grids and uses an isoparametric formulation for computing the dispersion and for second-order upwinding. Time integration is performed using an implicit approach. Three test cases are considered for comparing the numerical with analytical solutions. The scheme is second order in space and, when combined with Crank-Nicholson, is also second order in time. The results using Crank-Nicholson and full implicit (first-order) time integration methods are compared for problems with a variety of boundary conditions. For diffusion-dominated flows (PΔ ≤ 2) Crank-Nicholson works well, but for convection-dominated flows it produces spurious oscillations due to numerical dispersion errors. These oscillations are controlled by a flux limiter but only for Courant numbers below unity. It is shown that for high Courant and Péclet numbers a slight weighting of the time stepping toward fully implicit is effective against spurious oscillations and offers an optimum compromise between numerical dissipation and dispersion errors for a wide range of Courant and Péclet numbers. The scheme is shown to work on mildly nonorthogonal grids.
Get full access to this article
View all available purchase options and get full access to this article.
References
1.
Batu, V. (1989). “A generalized two-dimensional analytical solution for hydrodynamic dispersion in bounded media with the first-type boundary condition at the sources.” Water Resour. Res., 25(6), 1125–1132.
2.
Batu, V. (1993). “A generalized two-dimensional analytical solute transport model in bounded media for flux-type finite multiple sources.” Water Resour. Res., 29(8), 2881–2892.
3.
Cox, R. A., and Nishikawa, J. (1991). “A new total variation diminishing scheme for the solution of advective-dominant solute transport.” Water Resour. Res., 27(10), 2645–2654.
4.
Freeze, R. A., and Cherry, J. A. (1979). Groundwater, Prentice-Hall, Englewood Cliffs, N.J.
5.
Heinrich, J. C., Huyakorn, P. S., Zienkiewicz, O. C., and Mitchell, A. R. (1977). “An upwind finite element method for two-dimensional convective transport equation.” Int. J. Numer. Methods in Engrg., 11, 131–143.
6.
Hirsch, C. (1991). Numerical computation of internal and external flows, Wiley, New York.
7.
Huyakorn, P. S., and Pinder, G. F. (1983). Computational methods in subsurface flow, Academic, San Diego.
8.
Ijiri, Y., and Karasaki, K. (1994). “A Lagrangian-Eulerian finite element method with adaptive gridding for advection-dispersion problems.” Computational Methods in Water Resour., X, 291–298.
9.
Javendel, I., Doughty, C., and Tsang, C. F. ( 1984). Groundwater transport: Handbook of mathematical models, American Geophysical Union, Washington, D.C., 14–19.
10.
Leonard, B. P. (1979). “Stable and accurate convective modeling procedure based on quadratic upstream interpolation.” Comput. Methods Appl. Mech. Engrg., 19, 59–98.
11.
Peyret, R., and Taylor, T. D. (1983). Computational methods for fluid flow, Springer, New York.
12.
Pinder, G. F., and Shapiro, A. M. (1979). “A new collocation method for the solution of the convection dominated transport equation.” Water Resour. Res., 15(5), 1177–1182.
13.
Putti, M., Yeh, W. W. G., and Mulder, W. A. (1990). “A triangular finite volume approach with high resolution upwind terms for the solution of ground water transport equation.” Water Resour. Res., 26(12), 2865–2880.
14.
Sun, N. Z., and Yeh, W. W. G. (1993). “A proposed upstream weight numerical method for simulating pollutant transport in ground water.” Water Resour. Res., 19(6), 1489–1500.
15.
Van Genuchten, M. Th. ( 1977). “On the accuracy and efficiency of several numerical schemes for solving the convective-dispersive equation.” Finite element in water resources, W. G. Gray, G. F. Pinder, and I. A. Brebbia, eds.
16.
Verma, A. K., and Eswaran, V. (1996). “Overlapping control volume approach for convection-diffusion problems.” Int. J. Numer. Methods in Fluids, 23, 865–882.
17.
Verma, A. K., and Eswaran, V. (1997). “A bounded convection scheme for the overlapping control volume approach.” Int. J. Numer Method in Fluids, 25, 1137–1161.
18.
Wang, C., Sun, N., and Yeh, W. W. G. (1986). “An upstream weight multiple cell balance finite-element method for solving three-dimensional convection-dispersion equations.” Water Resour. Res., 22(11), 1575–1589.
19.
Westerink, J. J., and Shea, D. (1986). “Consistent higher degree Petrov-Galerkin methods for the solution of transient convection-diffusion equation.” Int. J. Numer. Methods in Engrg., 28, 1077–1100.
20.
Yeh, G. T. (1986). “An orthogonal upstream finite-element approach to modeling aquifer contaminant transport.” Water Resour. Res., 22(6), 952–964.
21.
Yeh, G. T. (1990). “A Lagrangian-Eulerian method with zoomable hidden fine-mesh approach for solving advection-dispersion equations.” Water Resour. Res., 26(6), 1133–1144.
22.
Yeh, G. T., and Chang, J. R. (1992). “An exact peak capturing and oscillation-free scheme to solve advection-dispersion transport equations.” Water Resour. Res., 28(11), 2937–2951.
23.
Yu, F. X., and Singh, V. P. (1995). “Improved finite-element method for solute transport.”J. Hydr. Engrg., ASCE, 121(2), 145–158.
Information & Authors
Information
Published In
History
Received: Mar 3, 1998
Published online: Jul 1, 2000
Published in print: Jul 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.