The paper deals with shape optimization problems governed by the Stokes equation with slip threshold boundary conditions. The stability of solutions to this state problem with respect to an appropriate class of domains is established.
On the basis of this result, existence of optimal shapes is proven. The second part of the paper is devoted to an approximation of the whole optimization problem and to convergence analysis.