A new scheme for the numerical simulation of incompressible flows is presented. The modeling equations are written in conservation-law form. The solution is advanced in time using the implicit Euler method, space discretization is accomplished through finite central difference expressions and the grid is colocated. Artificial dissipation terms are explicitly added by the user. These terms control numerical instabilities and avoid eventual odd-even decoupling. The scheme proposed solves a Poisson equation for the pressure itself, not for a pressure correction. In the actual format, the scheme is composed of three blocks: the first advances in time, simultaneously, the collection of physical equations, that is, momentum and energy, the second corresponds to the Poisson equation, which is iterated accordingly, and the third advances in time the turbulence equations (a k-ε high-Reynolds traditional model). The algorithm is very robust without resorting to any degree of relaxation, and the well-known cumbersome property related to staggered grids is totally removed. The scheme is second-order accurate in space and first order in time, but this low order in time is of no concern because the ultimate objective is the solution of steady state problems. Another very important novelty is the separation of the pressure terms in the physical equations in a pressure flux, what facilitates the tackling of different fluids. The new algorithm has been applied to the laminar and turbulent flows along a twodimensional duct, as well as to more demanding cases, namely, the lid-driven turbulent cavity, the two-dimensional turbulent internal flows about a backwardfacing step and along a duct with a high aspect ratio cavity located at the floor. The results are truly encouraging.