A unified approach is proposed for the prediction of heat transfer coefficients in turbulent falling films undergoing heating, evaporation or condensation for both of the cases with or without interfacial shear. A modified van Driest eddy viscosity model, which incorporated a damping factor f and takes into account the effect of variable shear stress, is used to predict the hydrodynamics of turbulent falling films. The calculated film thicknesses are in good agreement with the Nusselt-Brauer correlations for the non-sheared film and the Dukler prediction for highly sheared film. Also, by including a van Driest type turbulent Prandtl number model, the asymptotic heat transfer coefficients are accurately predicted and show better agreement with the extensive literature data and correlations than do most of the existing turbulence models proposed to date.