Accurate numerical schemes are proposed for solving incompressible Navier-Stokes equa-tions for 2D or 3D fluid flow problems. These are based on low-storage Runge-Kutta schemes for temporal discretization and fourth and sixth order compact schemes for spatial discretization. The incompressibility requirement is satisfied by solving a Poisson equation for pressure, with the same compact scheme used for discretization to ensure consistent global accuracy. The accuracy of the present procedure is demonstrated by application to several pertinent benchmark problems.