Simulation of propagating fronts in geothermal reservoirs with the implicit Leonard total variation diminishing scheme