A very simple and efficient finite element method is introduced for two and three dimensional viscous incompressible flows using the vorticity formulation. This method relies on recasting the traditional finite element method in the spirit of the high order accurate finite difference methods introduced by the authors in another work. Optimal accuracy of arbitrary order can be achieved using standard finite element or spectral elements. The method is convectively stable and is particularly suited for moderate to high Reynolds number flows.