This paper investigates turbulent flows with or without polymer additives in open channels and pipes. Equations of mean velocity, root mean square of velocity fluctuations, and energy spectrum are derived, in which the shear stress deficit model is used and the non-Newtonian properties are represented by the viscoelasticity ␣ * . The obtained results show that, with ␣ * increment, ͑1͒ the streamwise velocity fluctuations is increased, ͑2͒ the wall-normal velocity fluctuation is attenuated, ͑3͒ the Reynolds stress is reduced, and ͑4͒ there is a redistribution of energy from high frequencies to the low frequencies for the streamwise component, but dimensionless distribution over all frequencies almost remains the same as that in Newtonian fluid flows. Good agreement between the derived equations and experimental data in small drag-reduction regime is achieved, which indicates that the present model is workable for Newtonian/non-Newtonian fluid turbulent flows.