Abstract
We derive a general formulation of the time domain random walk (TDRW) approach to model the hydrodynamic transport of inert solutes in complex geometries and heterogeneous media. We demonstrate its formal equivalence with the discretized advection-dispersion equation and show that the TDRW is equivalent to a continuous time random walk (CTRW) characterized by space-dependent transition times and transition probabilities. The transition times are exponentially distributed. We discuss the implementation of different concentration boundary conditions and initial conditions as well as the occurrence of numerical dispersion. Furthermore, we propose an extension of the TDRW scheme to account for mobile-immobile multirate mass transfer. Finally, the proposed TDRW scheme is validated by comparison to analytical solutions for spatially homogeneous and heterogeneous transport scenarios.