A modified method of characteristics technique and mixed finite elements method for simulation of groundwater solute transport

Abstract

A comprehensive groundwater solute transport simulator is developed based on the modified method of characteristics (MMOC) combined with the Galerkin finite element method for the transport equation and the mixed finite element (MFE) method for the groundwater flow equation. The preconditioned conjugate gradient algorithm is used to solve the two large sparse algebraic system of equations arising from the MMOC and MFE discretizations. The MMOC takes time steps in the direction of flow, along the characteristics of the velocity field of the total fluid. The physical diffusion and dispersion terms are treated by a standard finite element scheme. The crucial aspect of the MMOC technique is that it looks backward in time, along an approximate flow path, instead of forward in time as in many method of characteristics or moving mesh techniques. The MFE procedure involves solving for both the hydraulic head and the specific discharge simultaneously. One order of convergence is gained by the MFE method, as compared with other standard finite element methods, and therefore more accurate velocity fields are simulated. The overall advantages of the MMOC‐MFE method include minimum numerical oscillation or grid orientation problems under steep concentration gradient simulations, and material balance errors are greatly reduced due to a very accurate velocity simulation by the MFE method. In addition, much larger time steps with Courant number well in excess of 1, as compared with the standard Galerkin finite element method, can be taken on a fixed spatial grid system without significant loss of accuracy.