The present study is devoted to propose a hybrid Green's function method to investigate the hyperbolic heat conduction problems. The difficulty of the numerical solutions of hyperbolic heat conduction problems is the numerical oscillation in the vicinity of sharp discontinuities. In the present study, we have developed a hybrid method combined the Laplace transform, Green's function and epsilon-algorithm acceleration method for solving time dependent hyperbolic heat conduction equation. From one- to three-dimensional problems, six different examples have been analyzed by the present method. It is found from these examples that the present method is in agreement with the Tsai-tse Kao's solutions [Tsai-tse Kao, Non-Fourier heat conduction in thin surface layers, J. Heat Transfer 99 (1977) 343-345] and does not exhibit numerical oscillations at the wave front. The propagation of the two- and three-dimensional thermal wave becomes so complicated because it occur jump discontinuities, reflections and interactions in these numerical results of the problem and it is difficult to find the analytical solutions or the result of other study to compare with the solutions of the present method. (C) 2009 Elsevier Ltd. All rights reserved.