Modeling groundwater flow to elliptical lakes and through multi-aquifer elliptical inhomogeneities

Two new analytic element solutions are presented for steady flow problems with elliptical boundaries. The first solution concerns groundwater flow to shallow elliptical lakes with leaky lake beds in a single-aquifer. The second solution concerns groundwater flow through elliptical cylinder inhomogeneities in a multi-aquifer system. Both the transmissivity of each aquifer and the resistance of each leaky layer may differ between the inside and the outside of an inhomogeneity. The elliptical inhomogeneity may be bounded on top by a shallow elliptical lake with a leaky lake bed. Analytic element solutions are obtained for both problems through separation of variables of the Laplace and modified-Helmholtz differential equations in elliptical coordinates. The resulting equations for the discharge potential consist of infinite sums of products of exponentials, trigonometric functions, and modified-Mathieu functions. The series are truncated but still fulfill the differential equation exactly; boundary conditions are met approximately, but up to machine accuracy provided enough terms are used. The head and flow may be computed analytically at any point in the aquifer. Examples are given of uniform flow through an elliptical lake, a well pumping near two elliptical lakes, and uniform flow through three elliptical inhomogeneities in a multi-aquifer system. Mathieu functions may be applied in a similar fashion to solve other groundwater flow problems in semi-confined aquifers and leaky aquifer systems with elliptical internal or external boundaries.

[1]  Michael G. Trefry,et al.  Surface water‐groundwater interaction near shallow circular lakes: Flow geometry in three dimensions , 2000 .

[2]  Aldo Fiori,et al.  Effective Conductivity of an Isotropic Heterogeneous Medium of Lognormal Conductivity Distribution , 2003, Multiscale Model. Simul..

[3]  Mark Bakker.  Steady groundwater flow through many cylindrical inhomogeneities in a multi-aquifer system , 2003 .

[4]  E. I. Anderson.  Conformal mapping of groundwater flow fields with internal boundaries , 2002 .

[5]  Mark Bakker,et al.  Analytic elements for multiaquifer flow , 2003 .

[6]  Charles R. Fitts.  Modeling three‐dimensional flow about ellipsoidal inhomogeneities with application to flow to a gravel‐packed well and flow through lens‐shaped inhomogeneities , 1991 .

[7]  J. Maxwell.  A Treatise on Electricity and Magnetism , 1873, Nature.

[8]  L. Townley,et al.  Influence of regional setting on the interaction between shallow lakes and aquifers , 2002 .

[9]  Two exact solutions for a cylindrical inhomogeneity in a multi-aquifer system , 2002 .

[10]  Karl Bandilla,et al.  Steady two‐dimensional groundwater flow through many elliptical inhomogeneities , 2004 .

[11]  E. I. Anderson.  An analytical solution representing groundwater–surface water interaction , 2003 .

[12]  Irene A. Stegun,et al.  Handbook of Mathematical Functions. , 1966 .

[13]  Fayez A. Alhargan,et al.  Algorithm 804: subroutines for the computation of Mathieu functions of integer orders , 2000, TOMS.

[14]  Aldo Fiori,et al.  Flow and transport in highly heterogeneous formations: 3. Numerical simulations and comparison with theoretical results , 2003 .

[15]  C. J. Hemker,et al.  Steady groundwater flow in leaky multiple-aquifer systems , 1984 .

[16]  M. Anderson,et al.  Seasonal reversals of groundwater flow around lakes and the relevance to stagnation points and lake budgets , 1981 .

[17]  A. Kacimov.  Three-dimensional groundwater flow to a lake: an explicit analytical solution , 2000 .

[18]  S. P. Neuman,et al.  Laplace-transform analytic element solution of transient flow in porous media , 2003 .

[19]  Randal J. Barnes,et al.  Two-dimensional flow through large numbers of circular inhomogeneities , 1999 .

[20]  L. Townley,et al.  A framework for quantitative analysis of surface water-groundwater interaction: Flow geometry in a vertical section , 1994 .

[21]  G. Bruggeman.  Analytical solutions of geohydrological problems , 1999 .

[22]  E. Veling,et al.  Elliptical inhomogeneities in groundwater flow-an analytical description , 1987 .

[23]  R. Hunt,et al.  Simulating Ground Water‐Lake Interactions: Approaches and Insights , 2003, Ground water.

[24]  Otto D Strack,et al.  Groundwater mechanics , 1989 .

[25]  H. M. Haitjema,et al.  Modeling double aquifer flow using a comprehensive potential and distributed singularities: 2. Solution for inhomogeneous permeabilities , 1981 .

[26]  Fayez A. Alhargan,et al.  Algorithms for the computation of all Mathieu functions of integer orders , 2000, TOMS.

[27]  D. Spencer,et al.  Field theory handbook: Including coordinate systems, differential equations, and their solutions , 1988 .

[28]  Henk M. Haitjema,et al.  Modeling Three‐Dimensional Flow in Confined Aquifers by Superposition of Both Two‐ and Three‐Dimensional Analytic Functions , 1985 .

[29]  C. Maas,et al.  The use of matrix differential calculus in problems of multiple-aquifer flow , 1986 .

[30]  Randal J. Barnes,et al.  High-order line elements in modeling two-dimensional groundwater flow , 1999 .

[31]  T. C. Winter.  Numerical simulation analysis of the interaction of lakes and ground water , 1976 .

[32]  C. R. Fitts.  Modeling three-dimensional flow about ellipsoidal inhomogeneities with application to flow to a gravel-packed well and flow through lens-shaped inhomogeneities : Water Resour Res V27, N5, May 1991, P815–824 , 1991 .

[33]  Definition of a capture zone for shallow water table lakes , 1988 .

[34]  Victor A Kelson,et al.  Improving a Regional Model Using Reduced Complexity and Parameter Estimation , 2002, Ground water.

[35]  Randal J. Barnes,et al.  Three-dimensional flow through large numbers of spheroidal inhomogeneities , 1999 .

[36]  T. N. Olsthoorn,et al.  Regional groundwater modeling of the Yucca Mountain site using analytic elements , 1999 .