“…Using the velocity field, each particle moves through the porous media by advection and hydrodynamic dispersion process. The particle's position in an unsteady nonuniform flow can be expressed by the Ito‐Taylor integration scheme (Fayazi & Ghazanfari, ; Hassan & Mohamed, ; Park et al, ; Sun et al, ): where ( x t + Δ t , y t + Δ t ) is the position of a particle at the time ' t + Δ t ' (L); ( x t , y t ) is the position of a particle at time ' t ' (L); r i ( i = 1, 2, 3, 4) is the Gaussian distribution of random variables with zero mean and unit variance.…”