In this paper a physical model for the anomalous temperature time evolution (decay) observed in complex thermodynamical system in presence of uniform heat source is provided. Measures involving temperatures T with power-law variation in time shows a different evolution of the temperature time rate T˙ (t) with respect to the temperature time-dependence T (t). Indeed the temperature evolution is a power-law increasing function whereas the temperature time rate is a power-law decreasing function of time. Such a behavior may be captured by a physical model that allows for a fast thermal energy diffusion close to the insulated location but must offer more resistance to the thermal energy flux as soon as the distance increases. In this paper this idea has been exploited showing that such thermodynamical system is represented by an heterogeneous onedimensional distributed mass one with power-law spatial scaling of its physical properties. The model yields, exactly a power-law evolution (decay) of the temperature field in terms of a real exponent as that is related to the power-law spatial scaling of the thermodynamical property of the system. The obtained relation yields a physical ground to the formulation of fractional-order generalization of the Fourier diffusion equation. R1) The discrete thermodynamic system described in eqs.(3,4) yields an approximation of powerlaw temporal evolution of the temperature defined by beta as the thermal conductivity c(V) and the specific heat coefficient C(V) varies along masses mj with the relation defined by alpha. The eqs. (3) and (4) do not explicitely contain the parameters alpha and beta. Is that an assumption? or is it a statement ?A1) The balance equations, in eq.(3) for the generic mass and in eq. (4) for the boundary mass m1, have been written in a generic form involving the non-homogeneity of the conductor in the thermal conductivity χ j and in the specific heat C j . Their specific variation, as power-laws of alpha, along the thermodynamic system provided in eqs.(5 a,b) corresponds to a power-law decay of the temperature with beta. It is a challenging problem whether other functional class of parameter variation may correspond to other measured time-varying evolution of the temperature fields. Therefore I decided not to specify the functional class of parameter variation in the generic balance equations. R2) Eq. (5a) and (5b) Why different denominators are different, it should be explained. A2) Thanks for the observation, I kept the same denominator and re-performed derivations in the revised manuscript. R3) Text after Eq.(5a) and (5b) "the real exponent alpha belongs to the interval -1<=alpha<=1". This interval should be -10 and b2<2 should be consider...