Modeling solute advection coupled with sorption kinetics in heterogeneous formations