A two-dimensional, boundary-layer program, STAN5, was modified to incorporate a low-Reynolds number version of the K-ε, two-equation turbulence model for the predictions of flow and heat transfer around turbine airfoils. The K-ε, two-equation model with optimized empirical correlations was used to account for the effects of free-stream turbulence and transitional flow. The model was compared with experimental flat plate data and then applied to turbine airfoil heat transfer prediction. A two-zone model was proposed for handling the turbulent kinetic energy and dissipation rate empirically at the airfoil leading edge region. The result showed that the predicted heat transfer coefficient on the airfoil agreed favorably with experimental data, especially for the pressure side. The discrepancy between predictions and experimental data of the suction surface normally occurred at transitional and fully turbulent flow regions.