We propose a method to enclose solutions for the stationary Navier-Stokes equation in nonconvex polygonal domains. Our method is based on an infinite dimensional Newton-type formulation by using the finite element method with constructive error estimates andfixed point theorems. Numerical examples related to the step flow problems in L-shape domain are presented. © 2007 by the University of Notre Dame. All rights reserved.