Numerical solutions for unsteady flow in unconfined aquifers

Abstract
Two numerical methods for solving the problem of unsteady flow in unconfined aquifers are studied. They are an explicit finite difference method (FDM), and the finite element method (FEM). The FEM is further subdivided into three schemes: vertical displacement approach, explicit scheme (FEM1), normal velocity approach, explicit scheme (FEM2), and vertical displacement approach, implicit scheme (FEM3). Results from the above methods are compared with experimental results from a sand box model. Various factors affecting the accuracy and numerical stability are investigated. Results indicate that, for a similar degree of accuracy, the FEM requires less computational effort than the explicit FDM. Amongst the three FEM schemes, FEM3 appears to be most attractive as it is the most stable and economical of the three schemes compared.