An efficient laplace transform solution for multiaquifer systems

Publisher Website
Google Scholar
Abstract
An efficient numerical algorithm is developed for groundwater flow in multiaquifer systems. The method is based on Herrera's formulation which eliminates the need of explicitly modeling aquitards. Laplace transform is applied to the governing integrodifferential equations, which removes not only the time dependence in the solution system, but also the integral expression. The final system is a set of nontransient, linear partial differential equations in two spatial dimensions. A regular finite difference method is used for its solution. The solution in the Laplace transform domain is then inverted to the time domain using either a slow or a quick scheme, depending on the requirement of speed and accuracy. Examples of a two‐aquifer one‐aquitard system involving radial and regional geometries are examimed.