TY - JOUR
T1 - Finite element method for subsurface hydrology using a mixed explicit‐implicit scheme
AU - Narasimhan, T. N.
AU - Neuman, S. P.
AU - Witherspoon, P. A.
N1 - Copyright:
Copyright 2015 Elsevier B.V., All rights reserved.
PY - 1978/10
Y1 - 1978/10
N2 - 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. This paper is not subject to U.S. copyright. Published in 1978 by the American Geophysical Union.
AB - 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. This paper is not subject to U.S. copyright. Published in 1978 by the American Geophysical Union.
UR - http://www.scopus.com/inward/record.url?scp=0018021457&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=0018021457&partnerID=8YFLogxK
U2 - 10.1029/WR014i005p00863
DO - 10.1029/WR014i005p00863
M3 - Article
AN - SCOPUS:0018021457
VL - 14
SP - 863
EP - 877
JO - Water Resources Research
JF - Water Resources Research
SN - 0043-1397
IS - 5
ER -