The paper is concerned with fully guaranteed and computable bounds of errors generated by Uzawa type methods for variational problems in the theory of visco-plastic fluids. The respective estimates have two forms. The first form contains global constants (such as the constant in the Friedrichs inequality for the respective domain), and the second one is based upon decomposition of the domain into a collection of subdomains and uses local constants associated with subdomains.