Abstract: It is quite notorious that the finite element solutions of transient heat conduction equation are very likely oscillating and/or over-bounding ones In order to eliminate the oscillation and/or over-bounding several criteria have been proposed by the author. In this paper a technique for smoothing the oscillation in the numerical results is presented. A high-order discretization scheme in time dimension is also discussed.