The paper deals with shape optimization of systems governed by the Stokes flow with threshold slip boundary conditions.Stability of solutions to the state problem with respect to a class of domains is established. For computational purposes the slip and impermeability conditions are regularized.
Numerical results are presented.