In this paper we present a model of a polluted groundwater site. The model consists of a coupled system of (advection)-diffusion-reaction equations for the groundwater level and the concentration of the pollutant. We use the complete flux scheme for the space discretization in combination with the #-method for time integration and we prove a new stability result for the scheme. Numerical results are computed for the Guarani Aquifer in South America and they show good agreement with results in literature.
Key words: advection-diffusion-reaction equation, finite volume method, complete flux scheme, stability analysis and groundwater flow.