Axial Groundwater Contaminant Dispersion Modeling for a Finite Heterogeneous Porous Medium

Abstract

In this study, a two-dimensional contaminant transport model with time-varying axial input sources subject to non-linear sorption, decay, and production is numerically solved to find the concentration distribution profile in a heterogeneous, finite soil medium. The axial input sources are assigned on the coordinate axes of the soil medium, with background sources varying sinusoidally with space. The groundwater velocities are considered space-dependent in the longitudinal and transversal directions. Various forms of axial input sources are considered to study their transport patterns in the medium. The alternating direction implicit (ADI) and Crank-Nicolson (CN) methods are applied to approximate the two-dimensional governing equation, and the obtained algebraic system of equations in each case is further solved by MATLAB scripts. Both approximate solutions are illustrated graphically for various hydrological input data. The influence of various hydrogeological input parameters, such as the medium’s porosity, density, sorption conditions, dispersion coefficients, etc., on the contaminant distribution is analyzed. Further, the influence of constant and varying velocity parameters on groundwater contaminant transport is studied. The stability of the proposed model is tested using the Peclet and Courant numbers. Substantial similarity is observed when the approximate solution obtained using the CN method is compared with the finite element method in a special case. The proposed approximate solution is compared with the existing numerical solutions, and an overall agreement of 98–99% is observed between them. Finally, the stability analysis reveals that the model is stable and robust.