An efficient and accurate numerical scheme is presented for the three-dimensional Navier-Stokes equations in primitive variables in a cylinder. The scheme is based on a spectral-Galerkin approximation for the space variables and a second-order projection scheme for time. The new spectral-projection scheme is implemented to simulate unsteady incompressible flows in a cylinder.