The procedure of numerical identification of thermophysical properties of concrete during its hardening is presented. Heat of cement hydration, thermal conductivity and specific heat are determined for purpose of modelling temperature evolution in massive concrete elements. The developed method is based on point temperature measurements in a cylindrical mould and the numerical solution of the inverse heat transfer problem by means of direct search optimization algorithm. The determined thermal characteristics of concrete are not constant and depend on the maturity of concrete. The procedure was verified on set of concrete mixes designed with Portland cement CEM I 42.5R and Portland-slag cement CEM II/B-S 32.5 N. Calcareous fly ash was also used for partial replacement of cement in the mixtures. The obtained results have been compared with experimentally measured temperature in concrete and a fair agreement has been found.