Three Iterative Finite Element Methods for the Stationary Smagorinsky Model