Two-dimensional hyperbolic heat conduction (HHC) problems with temperature-dependent thermal properties are investigated numerically. The present numerical method involves the hybrid application of the Laplace transform and control-volume methods. The Laplace transform technique is used to remove time-dependent terms, and then the transformed equation is discretized in the space domain by the control-volume formulation. Nonlinear terms induced by temperature-dependent thermal properties are linearized by using the Taylor's series approximation. In general, the numerical solution of the HHC problem has the phenomenon of the jump discontinuity in the vicinity of the thermal wave front. This phenomenon easily causes numerical oscillations in this region. In order to suppress these numerical oscillations, the selection of shape functions is an important task in the present study. The bi-hyperbolic shape function is introduced in the present control-volume formulation. Three examples involving a problem with a...