SUMMARYWe propose the level set method of topology optimization as a viable, robust and efficient alternative to density-based approaches in the setting of fluid flow. The proposed algorithm maintains the discrete nature of the optimization problem throughout the optimization process, leading to significant advantages over density-based topology optimization algorithms. Specifically, the no-slip boundary condition is implemented directly-this is accurate, removes the need for interpolation schemes and continuation methods, and gives significant computational savings by only requiring flow to be modeled in fluid regions. Topological sensitivity information is utilized to give a robust algorithm in two dimensions and familiar two-dimensional power dissipation minimization problems are solved successfully. Computational efficiency of the algorithm is also clearly demonstrated on large-scale three-dimensional problems.