A new numerical method for solving the solute transport equation