Variance methods for thermo load of elements of IC engines