The numerical modeling of the fully developed Poiseuille flow of a Newtonian fluid in a square section with slip yield boundary condition at the wall is presented. The stick regions in outer corners and the slip region in the center of the pipe faces are exhibited. Numerical computations cover the complete range of the dimensionless number describing the slip yield effect, from a full slip to a full stick flow regime. The resolution of variational inequalities describing the flow is based on the augmented Lagrangian method and a finite element method. The localization of the stick-slip transition points is approximated by an anisotropic auto-adaptive mesh procedure. The singular behavior of the solution at the neighborhood of the stick-slip transition point is investigated.