A robust and efficient two-dimensional numerical model for the simulation of river pollutant and thermal dynamics