Numerical simulations of solute transport in highly heterogeneous formations: A comparison of alternative numerical schemes