This investigation concerns numerical calculation of fully developed turbulent forced convective heat transfer and fluid flow in ducts over a wide range of Reynolds numbers. The low Reynolds number version of a non-linear eddy viscosity model is combined with a two-equation heat flux model with the eddy diffusivity concept. The model can theoretically be used for a range of Prandtl numbers or a range of at different fluids. The computed results compare satisfactory with the. available experiment. Based on existing DNS data and calculations in this work the ratio between the time-scales (temperature to velocity) is found to be approximately 0.7. In light of this assumption an algebraic scalar flux model with variable diffusivity is presented.