We investigate the phenomenon of drag reduction in a viscoelastic fluid model of dilute polymer solutions. By means of direct numerical simulations of the three-dimensional turbulent Kolmogorov flow we show that drag reduction takes place above a critical Reynolds number Re(c). An explicit expression for the dependence of Re(c) on polymer elasticity and diffusivity is derived. The values of the...