The hysteresis effect on the vortex induced vibration (VIV) on a circular cylinder is investigated by the numerical solution of the two-dimensional Reynolds averaged Navier-Stokes equations. An upwind and total variation diminishing (TVD) conservative scheme is used to solve the governing equations written in curvilinear coordinates and the k-ɛ turbulence model is used to simulate the turbulent flow in the wake of the body. The cylinder is supported by a spring and a damper and free to vibrate in the transverse direction. In previous work, numerical results for the amplitude of oscillation and vortex shedding frequency were compared to experimental data obtained from the literature to validate the code for VIV simulations. In the present work, results of practical interest are presented for the power absorbed by the system, phase angle, amplitude, frequency, and lift coefficient. The numerical results indicate that the hysteresis effect is observed only when the frequency of vortex shedding gets closer to the natural frequency of the structure in air.