Adaptive-grid simulation of groundwater flow in heterogeneous aquifers

The prediction of contaminant transport in porous media requires the computation of the flow velocity. This work presents a methodology for high-accuracy computation of flow in a heterogeneous isotropic formation, employing a dual-flow formulation and adaptive gridding. The dual equations, describing the hydraulic head and the streamfunction, are numerically solved through finite element approximations. The application of classic finite-element methods requires a rather large number of nodes to represent suitably the flow in high-contrast formations. We present a mesh-adaptive approach that enhances the accuracy of the numerical flow solution for a given computational effort. We rely on an a posteriori error estimator to identify areas where refinements of the finite element mesh are needed or unrefinements are acceptable. We also demonstrate through numerical experiments that the developed methodology efficiently enhances accuracy through successive mesh adaptation.