| Abstract: | More accurate alternatives to the widely used harmonic mean interblock transmissivity are proposed for block-centered finite-difference models of ground-water flow in unconfined aquifers and in aquifers having smoothly varying transmissivity. The harmonic mean is the exact interblock transmissivity for steady-state one-dimensional flow with no recharge if the transmissivity is assumed to be spatially uniform over each finite-difference block, changing abruptly at the block interface. However, the harmonic mean may be inferior to other means if transmissivity varies in a continuous or smooth manner between nodes. Alternative interblock transmissivity functions are analytically derived for the case of steady-state one-dimensional flow with no recharge. The second author has previously derived the exact interblock transmissivity, the logarithmic mean, for one-dimensional flow when transmissivity is a linear function of distance in the direction of flow. We show that the logarithmic mean transmissivity is also exact for uniform flow parallel to the direction of changing transmissivity in a two- or three-dimensional model, regardless of grid orientation relative to the flow vector. For the case of horizontal flow in a homogeneous unconfined or water-table aquifer with a horizontal bottom and with areally distributed recharge, the exact interblock transmissivity is the unweighted arithmetic mean of transmissivity at the nodes. This mean also exhibits no grid-orientation effect for unidirectional flow in a two-dimensional model. For horizontal flow in an unconfined aquifer with no recharge where hydraulic conductivity is a linear function of distance in the direction of flow the exact interblock transmissivity is the product of the arithmetic mean saturated thickness and the logarithmic mean hydraulic conductivity. For several hypothetical two- and three-dimensional cases with smoothly varying transmissivity or hydraulic conductivity, the harmonic mean is shown to yield the least accurate solution to the flow equation of the alternatives considered. Application of the alternative interblock transmissivities to a regional aquifer system model indicates that the changes in computed heads and fluxes are typically small, relative to model calibration error. For this example, the use of alternative interblock transmissivities resulted in an increase in computational effort of less than 3 percent. Numerical algorithms to compute alternative interblock transmissivity functions in a modular three-dimensional flow model are presented and documented. |
| Genre: | USGS Numbered Series |
| ProdID: | 27364 |
| Citation Author: | Goode, D. J.; Appel, C. A. |
| Citation Contributing Office: | |
| Citation Datum: | |
| Citation Day: | |
| Citation Edition: | - |
| Citation Editor: | |
| Citation End Page: | |
| Citation Issue: | |
| Citation Keywords: | |
| Citation Language: | ENGLISH |
| Citation Larger Work Title: | |
| Citation LatN: | |
| Citation LatS: | |
| Citation LonE: | |
| Citation LonW: | |
| Citation Month: | |
| Citation No Pagination: | |
| Citation Number Of Pages: | |
| Citation Online Only Flag: | |
| Citation Phsyical Description: | viii, 79 p. :ill. ;28 cm. |
| Citation Projection: | |
| Citation Public Comments: | |
| Citation Publisher: | U.S. Geological Survey ; |
| Citation Series: | Water-Resources Investigations Report |
| Citation Series Code: | WRI |
| Citation Series Number: | 92-4124 |
| Citation Search Results Text: | Finite-difference interblock transmissivity for unconfined aquifers and for aquifers having smoothly varying transmissivity; 1992; WRI; 92-4124; Goode, D. J.; Appel, C. A. |
| Citation Start Page: | |
| Citation Volume: | |
| Citation Year: | 1992 |
| Type: | citation/reference |
| Text: | Finite-difference interblock transmissivity for unconfined aquifers and for aquifers having smoothly varying transmissivity; 1992; WRI; 92-4124; Goode, D. J.; Appel, C. A. |
| URL (THUMBNAIL): | http://pubs.usgs.gov/wri/1992/4124/report-thumb.jpg |
| URL (DOCUMENT): | http://pubs.usgs.gov/wri/1992/4124/report.pdf |
| Date Other: | Sat, 1 Jan 1994 00:00 -0600 |
| Publisher: | U.S. Geological Survey ; |
