We present recent results using adaptive finite element methods, based on a posteriori error estimates, to compute various output functionals for incompressible flow problems in 3d, for both laminar and turbulent flows. The a posteriori error estimates are based on the solution of an associated dual problem with data connected to the output functional we want to compute.