We present a new methodology for the numerical resolution of the hydrodynamicsof incompressible viscid newtonian fluids. It is based on the Navier-Stokesequations and we refer to it as the vorticity projection method. The method is robust enough to handle complex and convoluted configurationstypical to the motion of biological structures in viscous fluids.Although the method is applicable to three dimensions, we address herein detail only the two dimensional case. We provide numerical data forsome test cases to which we apply the computational scheme.