Solute transport in a semi-infinite homogeneous aquifer with a fixed point source concentration

This study derives an analytical solution of the advection–dispersion (AD) equation commonly used to describe the transport of pollutants in a semi-infinite homogeneous aquifer. When an extra constant source term is added to the AD equation, it changes the solution of the equation. The AD equation is solved analytically using Laplace transform. Also, the equation is solved numerically using an explicit finite difference method and its stability condition is presented with the aid of matrix method. For the solution of the AD equation the following considerations are made: (1) The dispersion and velocity are considered as time-dependent; (2) dispersion is expressed as directly proportional to the square of velocity; (3) there is also diffusion; (4) there is some initial concentration and the aquifer domain is, therefore, not pollutant-free; (5) There is a time-dependent exponentially decreasing input source; and (6) the concentration gradient is assumed to be zero at the exit boundary. It is found that the contaminant concentration decreases with time contrary to what happens when the extra term is not included.