The paper is devoted to application of the multi-domain algorithms based on the non-overlapping iterative domain decomposition method coupled with the influence matrix method for solution of the incompressible Navier-Stokes equations in vorticity-streamfunction formulation (omega-psi).