Results of direct numerical simulation (DNS) of incompressible homogeneous and isotropic forced turbulence with grid points up to 1024 3 are presented with emphasis on the time dependence of the turbulence statistics. The Taylor micro scale Reynolds number at a quasi-stationary state is 283 in the DNS with 1024 3 grid points, and the energy spectrum exhibits an universal equilibrium range. The results suggest that for R-lambda > 280 or so, the energy dissipation rate normalized by the integral length scale and root-mean-square velocity is almost independent of R-lambda and is 0.44 - 0.45.