SUMMARYThe numerical solution of the time-dependent Navier-Stokes equations in terms of the vorticity and a stream function is a well tested process to describe two-dimensional incompressible ows, both for uid mixing applications and for studies in theoretical uid mechanics. In this paper, we consider the interaction between the unsteady advection-di usion equation for the vorticity, the Poisson equation linking vorticity and stream function and the approximation of the boundary vorticity, examining from a practical viewpoint, global iteration stability and error. Our results show that most schemes have very similar global stability constraints although there may be small stability gains from the choice of method to determine boundary vorticity. Concerning accuracy, for one model problem we observe that there were cases where the boundary vorticity discretization did not propagate to the interior, but for the usual cavity ow all the schemes tested had error close to second order.