A Moving Least Squares Based Meshless Element-Free Galerkin Method for the Coupled Simulation of Groundwater Flow and Contaminant Transport in an Aquifer

Abstract

In this study, a meshless simulation model based on the element-free Galerkin method (EFGM) is proposed for the coupled simulation of groundwater flow and contaminant transport (GFCT) in the two-dimensions. EFGM discretizes the governing GFCT equations using the Galerkin integral approach and shape functions generated through the moving least squares (MLS) approximation method. It results into the stable system of equations for estimating the groundwater head and contaminant concentration values. EFGM computes the several global matrices using the spatial information of the scattered and unconnected field nodes and hence, terminates the meshing procedure unlike the finite difference (FDM) and finite element (FEM) methods. The Kronecker delta function property is missing in EFGM shape functions and hence, the Dirichlet boundary conditions are enforced with some special methods. We have used the penalty method to include them into the presented EFGM model. This method is straightforward and computationally useful for the real field simulations. The proposed model is first validated using both the one-and two-dimensional analytical solutions. The developed coupled model is then applied to solve the hypothetical and real field problems. EFGM results for both the problems are compared with the MOD-FLOW and MT3DMS results for the head and concentration values respectively. EFGM results are found to be satisfactory for both the cases, and it shows the potential of the proposed approach for the coupled GFCT simulations.

Axial Groundwater Contaminant Dispersion Modeling for a Finite Heterogeneous Porous Medium

“…Further, the initial background source concentration is assumed to be zero (i.e., c 0 = 0). Data from [45] (a 1 = 0.1, a 2 = 0.1, u 0 = 1.05 km/year, v 0 = 0.15 km/year, D x0 = 1.1 km 2 /year, D y0 = 0.1 km 2 /year) and [46] (D x0 = 5 × 0.24 m 2 /day, D y0 = 0.25 m 2 /day, u 0 = 0.24 m/day, v 0 = 0.25 m/day, a 1 = 0, a 2 = 0) were fed into the proposed model, and its output was compared to the output obtained by the data considered by the authors. As evident, a similar overall trend of the contaminant profile-a spatially diminishing concentration trend towards the edge of the domain, is observed between the various existing models and the proposed model for the respective datasets considered.…”