In this chapter we develop a numerical model of the transport of chemical mass within flowing groundwater. Starting with a derivation of the advection–dispersion equation, we consider how the equation may be cast in discrete form as a finite difference equation. We further discuss the numerical stability of the calculation procedure and the numerical dispersion inherent in the discrete solution, concluding with a worked example.