Premium
Finite element method for subsurface hydrology using a mixed explicit‐implicit scheme
Author(s) -
Narasimhan T. N.,
Neuman S. P.,
Witherspoon P. A.
Publication year - 1978
Publication title -
water resources research
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.863
H-Index - 217
eISSN - 1944-7973
pISSN - 0043-1397
DOI - 10.1029/wr014i005p00863
Subject(s) - finite element method , galerkin method , weighting , boundary value problem , computation , aquifer , groundwater flow , mathematics , flow (mathematics) , subsurface flow , rotational symmetry , computer science , groundwater , geotechnical engineering , mathematical optimization , calculus (dental) , geology , algorithm , geometry , mathematical analysis , engineering , structural engineering , medicine , dentistry , radiology
The mixed explicit‐implicit Galerkin finite element method developed previously by the authors is shown to be ideally suited for a wide class of problems arising in subsurface hydrology. These problems include confined saturated flow, unconfined flow under free surface conditions subject to the Dupuit assumption, flow in aquifers which are partly confined and partly unconfined, axisymmetric flow to a well with storage, and flow in saturated‐unsaturated soils. A single computer program, entitled Flump, can now handle all of these problems. The mixed explicit‐implicit solution strategy employed in the program insures a high level of accuracy and computation efficiency in most cases. It eliminates many of the difficulties that groundwater hydrologists have been encountering in trying to simulate extensive aquifer systems by finite elements. Some of the outstanding features of this solution strategy include an automatic control of time step size, reclassification of nodes from explicit to implicit during execution, automatic adjustment of the implicit time‐weighting factor, and the treatment of boundary conditions and source terms as arbitrary functions of time of the state of the system. Five examples are presented to demonstrate the versatility and power of this new approach. A purely physical derivation of the finite element equations which does not rely on the Galerkin formalism is also included in one of the appendices.