A modified form of Carey'ls method for predicting the heat transfer to a turbulent falling film is presented. The new model is successfully tested against experimental data as well as other models. This method utilizes an optimization technique along with the Rombcrg integration to determine two universal constants for the modified turbulent diffusivisity model. The results show that in spite of its relative simplicity, the new model is a fairly accurate one, having an overall average error of less than 9%.