in this paper, numerical simulations of incompressible flows around rotating circular cylinders have been performed. the two-dimensional navier-stokes equations are solved by using a cartesian non-uniform grid. the immersed boundary method (ibm) with the virtual physical model (vpm) was used in order to model the presence of the circular cylinder in the flow. the fractional time step method was used to coupling the pressure and velocity fields. the simulations were carried out for reynolds numbers equals to 60, 100 and 200 for different specific rotations. the effects of rotation on flow characteristics and fluctuating forces were investigated. the strouhal number, obtained by performing the fast fourier transform (fft) of the temporal distribution of the lift coefficient, and the pressure coefficients, were also been calculated. vorticity contours are presented considering different values of the reynolds number and specific rotation. the numerical results obtained are compared to those obtained by other authors and the usefulness of the numerical methodology composed by the combination of the ibm with the vpm to simulate flows in the presence of mobile bodies is highlighted.