In this article we present efficient numerical methods for the Navier Stoke
s equations with slip boundary conditions. A rst method is based on a saddl
e-point formulation of the slip boundary condition. A congruent gradient (C
G) method is applied to the Schur complement operator in order to solve the
problem. We present two preconditioners for the CG-method which result in
convergence rates independent of the mesh size. For a second method the sli
p boundary condition is enforced pointwise for nodal values of the velocity
at boundary nodes.