In this paper, we present a finite difference method for the implementation of the rotation of a circular cylinder in the incompressible flow field by solving the two-dimensional unsteady Navier-Stokes equations. The approach is to use staggered grid method so that the accuracy and order of convergence of the associated algorithms can be maintained. The proposed method is easy to be implemented and is effective. A set of simulations for the flow dynamics is provided to show the computational results.