A vorticity-velocity method for the numerical solution of 3D incompressible flows