Direct-Formulation Finite Element (DFFE) Method for Groundwater Flow Modeling: Two-Dimensional Case

Abstract

Real-world groundwater modeling deals with heterogeneous and anisotropic aquifer systems. Numerical models, particularly the finite difference and finite element methods, are extensively used for modeling such a complex aquifer system. The finite difference method is an age-old straightforward approach, whereas the finite element method is more idealistic, relatively new, and has emerged as a popular approach. The unique advantages with the finite element method are (i) an irregular-shaped domain can be discretized into elements with a matching boundary, (ii) a node can be connected with the neighboring nodes in any direction. Thus, the finite element analysis allows water/solute to flow more freely toward the surrounding nodes along nodal links. On the other hand, the mathematical formulation of the finite element technique follows a set of rules step by step. Conventionally, this idealistic approach uses shape functions and integration techniques to develop element matrices, which, in turn, are assembled to form the global matrix. This conventional global matrix was found to be inconsistent with the basic partial differential equation representing the mass balance in a nodal domain as a control volume. As a result, the matrix solution, particularly for nonlinear problems (e.g., vadose zone flow, density dependent flow, and so forth) becomes unstable and computationally burdensome. To overcome the limitations of the conventional finite element method for real-world groundwater modeling, a new approach, called the direct-formulation finite element (DFFE) method, has been developed. In this approach, a two-dimensional model domain is discretized into triangular elements. Element matrices are formed directly based on flow principles and material properties of elements without using shape functions and integration process. The DFFE method is robust. The global matrix is consistent, and the solution is stable and efficient. The results obtained from this approach were found in a good agreement with that obtained from analytical and other widely used models.