Particle velocity interpolation in block‐centered finite difference groundwater flow models

A block-centered, finite difference model of two-dimensional groundwater flow yields velocity values at the midpoints of interfaces between adjacent blocks. Method of characteristics, random walk and particle-tracking models of solute transport require velocities at arbitrary particle locations within the finite difference grid. Particle path lines and travel times are sensitive to the spatial interpolation scheme employed, particularly in heterogeneous aquifers. This paper briefly reviews linear and bilinear interpolation of velocity and introduces a new interpolation scheme. Linear interpolation of velocity is consistent with the numerical solution of the flow equation and preserves discontinuities in velocity caused by abrupt (blocky) changes in transmissivity or hydraulic conductivity. However, linear interpolation yields discontinuous and somewhat unrealistic velocities in homogeneous aquifers. Bilinear interpolation of velocity yields continuous and realistic velocities in homogeneous and smoothly heterogeneous aquifers but does not preserve discontinuities in velocity at abrupt transmissivity boundaries. The new scheme uses potentiometric head gradients and offers improved accuracy for nonuniform flow in heterogeneous aquifers with abrupt changes in transmissivity. The new scheme is equivalent to bilinear interpolation in homogeneous media and is equivalent to linear interpolation where gradients are uniform. Selecting the best interpolation scheme depends, in part, on the conceptualization of aquifer heterogeneity, that is, whether changes in transmissivity occur abruptly or smoothly.

[1]  D. W. Pollock.  Semianalytical Computation of Path Lines for Finite‐Difference Models , 1988 .

[2]  George F. Pinder,et al.  Finite difference model for aquifer simulation in two dimensions with results of numerical experiments , 1976 .

[3]  T. A. Prickett,et al.  A "random-walk" solute transport model for selected groundwater quality evaluations , 1981 .

[4]  C. Axness,et al.  Three‐dimensional stochastic analysis of macrodispersion in aquifers , 1983 .

[5]  David B. Grove,et al.  The use of Galerkin finite-element methods to solve mass-transport equations , 1977 .

[6]  A. Davis.  Deterministic Modeling of Dispersion in Heterogeneous Permeable Media , 1986 .

[7]  Randel Haverkamp,et al.  A Comparison of Numerical Simulation Models For One-Dimensional Infiltration1 , 1977 .

[8]  Arlen W. Harbaugh,et al.  A modular three-dimensional finite-difference ground-water flow model , 1984 .

[9]  A. El-Kadi.  Remedial actions under variability of hydraulic conductivity , 1988 .

[10]  J. R. Macmillan,et al.  Stochastic analysis of spatial variability in subsurface flows: 2. Evaluation and application , 1978 .

[11]  J. Bredehoeft,et al.  Computer model of two-dimensional solute transport and dispersion in ground water , 1978 .

[12]  Franklin W. Schwartz,et al.  Mass transport: 1. A stochastic analysis of macroscopic dispersion , 1980 .

[13]  R. Ababou,et al.  Three-dimensional flow in random porous media , 1988 .

[14]  S. P. Neuman,et al.  Eulerian-Lagrangian Methods for Advection-Dispersion , 1982 .

[15]  D. W. Peaceman,et al.  Numerical Calculation of Multidimensional Miscible Displacement by the Method of Characteristics , 1964 .

[16]  W. G. Gray,et al.  Finite Element Simulation in Surface and Subsurface Hydrology , 1977 .

[17]  Charles A. Appel,et al.  A note on computing finite difference interblock transmissivities , 1976 .