On Galerkin models for transport in groundwater

The Galerkin finite element model (GFEM) may provide oscillatory results when employed to predict contaminant transport in groundwater unless a very fine mesh is used. Adaptation of a very fine mesh may make the application of the GFEM impractical to field problems. The Petrov—Galerkin finite element models (PGFEMs) can provide oscillation free results for relatively coarser mesh. However, the PGFEM violates the Galerkin principle and introduces large “numerical” dispersion. The objective of this paper has been to develop accurate criteria to improve the applicability of the GFEM to obtain oscillation free accurate results for coarser mesh and compare its performance with that of the PGFEM. It has been shown that the GFEM provides oscillation free accurate results for coarser mesh with Peclet number Pe 20. Further, the GFEM prediction has always been more accurate than the PGFEM for a variety of source configurations and flow fields.