We introduce a new method for computing a posteriori bounds on engineering outputs from finite element discretizations of the incompressible Stokes equations. The method results from recasting the output problem as a minimization statement without resorting to an error formulation. The minimization statement engenders a duality relationship which we solve approximately by Lagrangian relaxation. We demonstrate the method for a stabilized equal-order approximation of Stokes flow, a problem to which previous output bounding methods do not apply. The conceptual framework for the method is quite general and shows promise for application to stabilized nonlinear problems, such as Burger's equation and the incompressible Navier-Stokes equations, as well as potential for compressible flow problems.Singapore-MIT Alliance (SMA