The paper deals with the approximation, convergence analysis and implementation of the Stokes system with threshold slip boundary conditions of Navier type. Based on the fixed-point formulation we prove the existence of a solution for a class of solution-dependent slip functions g satisfying an appropriate growth condition and its uniqueness provided that g is one-sided Lipschitz continuous.
Further we study under which conditions the respective fixed-point mapping is contractive. To discretize the problem we use P1-bubble/P1 elements.
Properties of discrete models in dependence on the discretization parameter are analysed and convergence results are established. In the second part of the paper we briefly describe the duality approach used in computations and present results of a model example.