Laplace Transform Analytic Element Method for Transient Groundwater Flow Simulation

Abstract

The Laplace transform analytic element method (LT-AEM), applies the traditionally steady-state analytic element method (AEM) to the Laplace-transformed diffusion equation (Furman and Neuman, 2003). This strategy preserves the accuracy and elegance of the AEM while extending the method to transient phenomena. The approach taken here utilizes eigenfunction expansion to derive analytic solutions to the modified Helmholtz equation, then back-transforms the LT-AEM results with a numerical inverse Laplace transform algorithm. The two-dimensional elements derived here include the point, circle, line segment, ellipse, and infinite line, corresponding to polar, elliptical and Cartesian coordinates. Each element is derived for the simplest useful case, an impulse response due to a confined, transient, single-aquifer source. The extension of these elements to include effects due to leaky, unconfined, multi-aquifer, wellbore storage, and inertia is shown for a few simple elements (point and line), with ready extension to other elements. General temporal behavior is achieved using convolution between these impulse and general time functions; convolution allows the spatial and temporal components of an element to be handled independently.Comparisons are made between inverse Laplace transform algorithms; the accelerated Fourier series approach of de Hoog et al. (1982) is found to be the most appropriate for LT-AEM applications. An application and synthetic examples are shown for several illustrative forward and parameter estimation simulations to illustrate LT-AEM capabilities. Extension of LT-AEM to three-dimensional flow and non-linear infiltration are discussed.