The non-Fourier model for heat transfer leads to a hyperbolic evolution problem describing the temperature solution. A one-dimensional case is considered for such a propagating heat wave reflected at a boundary. A primary goal is investigation of the effectiveness of numerical solution techniques for the case of a propagating heat front and the influence of different boundary conditions. Finite elements are employed in space and alternative time integration schemes are studied, including ordinary differential equation system integrators. The effect of solution “roughness” on the error and oscillations before and after reflection is examined and rates of convergence are numerically determined.