An Explicit Finite Difference Model for Unconfined Aquifers

Abstract

ABSTRACTMost of the current simulation models for unconfined aquifers are based on the assumption that the free surface variation is small so that it can be combined with permeability to reduce the nonlinear Boussinesq equation to a linear partial differential equation (PDE). One of the most obvious reasons for using the linearization assumption is for the ease of numerical solution. This work presents a simpler alternative which permits an easy direct solution of the Boussinesq equation. A forward in time, central in space (FTCS) explicit finite difference method is used in the simulation model. The model was first validated by comparing its results with known analytical solution. It was then applied to an actual situation in which the short‐term responses (from pumping) of an unconfined aquifer were simulated. The study shows that the stability of the model can be easily controlled, and because of the simple algorithm used, the code can be expeditiously developed and economically run on smaller machines. Due to the uncertainties in the calibration results, it is recommended here that more data be collected to improve the calibration before the model is used as a real‐time simulation tool.