Abstract: | The advection–dispersion equation with spatially variable coefficients does not have an exact analytical solution and is therefore solved numerically. However, solutions obtained with several of the traditional finite difference or finite element techniques typically exhibit spurious oscillation or numerical dispersion when advection is dominant. The mixing cell and semi-analytical solution methods proposed in this study avoid such oscillation or numerical dispersion when advection dominates. Both the mixing cell and semi-analytical solution methods calculate the spatial step size by equating numerical dispersion to physical dispersion. Because of the spatial variability of the coefficients the spatial step size varies in space. When the time step size Δt→0, the mixing cell method reduces to the semi-analytical solution method. The results of application to two cases show that the mixing cell and semi-analytical solution methods are better than a finite difference method used in the study. © 1998 John Wiley & Sons, Ltd. |