Three-dimensional numerical simulation of a uniform incompressible viscous flow around a stationary circular cylinder was conducted. The CFX-4 software was used to calculate the hydrodynamic characteristics of the flow and the finite volume method for incompressible Navier-Stokes equations was employed in the program. The simulation of the flow was performed for Re=10 3 and Re=10 4 respectively within the sub-critical region. In order to overcome numerical instability for the high Reynolds number flows, a quadratic upwind scheme was incorporated for the Navier-Stokes equations. The periodicity boundary condition was used at the ends of the cylinder. It was found that the evolution of the lift and drag coefficients in each plane along the cylinder span is different. Comparison between the predicted results based on the three-dimensional and the two-dimensional analysis was also given. It is concluded that at the high Reynolds number the effect of three-dimensionality of the flow around the circular cylinder is remarkable, and in addition hydrodynamic coefficients with of 3-D simulation are less than those given by 2-D simulation. Three-dimensional numerical simulation of a uniform incompressible viscous flow around a stationary circular cylinder was conducted. The CFX-4 software was used to calculate the hydrodynamic characteristics of the flow and the finite volume method for incompressible Navier-Stokes equations was employed in the program . The simulation of the flow was performed for Re = 10 3 and Re = 10 4 respectively within the sub-critical region. In order to overcome numerical instability for the high Reynolds number flows, a quadratic upwind scheme was incorporated for the Navier-Stokes equations. The periodicity boundary condition was used at the ends of the cylinder. It was found that the evolution of the lift and drag coefficients in each plane along the cylinder span is different. Comparison between the predicted results based on the three-dimensional and the two-dimensional analysis was also given. It is concluded that at the high Reynolds number the effect of three-dimensionality of the flow aroun d the circular cylinder is remarkable, and in addition hydrodynamic coefficients with of 3-D simulation are less than those given by 2-D simulation.