We present a (physically based) Lagrangian approach to simulate diffusive mixing processes on the pore scale beyond perfectly mixed conditions. Results show the feasibility of the approach for reproducing measured mixing times and concentrations of isotopes over pore sizes and that typical shapes of breakthrough curves (normally associated with non-uniform transport in heterogeneous soils) may also occur as a result of imperfect subscale mixing in a macroscopically homogeneous soil matrix.