Thermoelastic fracture solutions using distributions of singular influence functions—II: Numerical modelling of thermally self-driven cracks