Summary: Three different solution methods for the finite-volume discretized incompressible Navier-Stokes equations have been tested: segregated approaches, built around coupling methods such as SIMPLE, SIMPLER and PISO plus a line Gauss-Seidel linear solver; coupled methods, incorporating a penalty formulation to eliminate zero diagonal elements in the coefficient matrices, plus preconditioned GMRES as a linear solver; and a FAS-full multigrid algorithm accelerating a classical segregated method based upon SIMPLE and the line Gauss-Seidel solver. Results demonstrate that the coupled method compares favourably to the segregated technique at small grid sizes but becomes too expensive for large problems. The FAS-full multigrid algorithm outperforms the other two methods when large numbers of nodes are employed in the simulation.