“…An unconditionally stable Crank-Nicolson finite difference scheme was used by Ravi et al (2014) to analyze the constant and longitudinal dispersion profile of contaminants. In his studies, Ravi et al (2014) did not take into account the time-dependent dispersions coefficients, the decay rate constant and the zero-order production rate coefficient of solute in the liquid phase as done by Das et al (2017) , Guleria et al (2020) with a linear dispersion advection equation model. Our objective is to exploit the more stable fourth-order Runge-Kutta method (RK4) to evaluate the profile of C (X, T) contaminants through a porous medium and to conduct a comparative study between the profiles of the contaminants obtained from an asymptotic and linear dispersion coefficient, and on the other hand, from the initial and boundary conditions used in the work of Das et al (2017) .…”