Heat diffusion processes are generally modeled based on Fourier’s law to estimate how the temperature propagates inside a body. This type of modeling leads to a parabolic partial differential equation, which predicts an infinite thermal wave speed of propagation. However, experimental evidence shows that diffusive processes occur with a finite velocity of thermal propagation in many applications. In this paper, we develop a mathematical formulation to predict the finite speed of heat propagation in multidimensional phase change problems. The model generalizes the enthalpy formulation by adding a hyperbolic term. The governing equations are simulated by the finite element method. The proposed model is first verified by comparing numerical and experimental results illustrating the difference between the infinite and finite propagation velocity for heat inside biological tissues. Then, the results of the two and three-dimensional numerical solution of the continuous steel casting process are presented. We will illustrate that the effects of the initial conditions vanish faster when using the parabolic equation, while they persist in the hyperbolic modeling approach. The results demonstrate significant differences in the initial thermal dynamics and at the solid-liquid interface position when adding the hyperbolic term. The changes are more noticeable in the regions of the steel beam where rapid heat loss and, consequently, faster phase change occur.