In this paper we present a discontinuous least-squares spectral element method for Stokes equations with primitive variable formulation on both smooth and non-smooth domains. We propose an exponentially accurate numerical scheme based on the stability estimates and implement it on different polygonal domains. Preconditioned conjugate gradient method is used to obtain the numerical solution.