Removing geometric details from the computational domain can significantly reduce the complexity of downstream task of meshing and simulation computation, and increase their stability. Proper estimation of the sensitivity analysis error induced by removing such domain details, called defeaturing errors, can ensure that the sensitivity analysis fidelity can still be met after simplification. In this paper, estimation of impacts of removing arbitrarily constrained domain details to the analysis of incompressible fluid flows is studied with applications to fast analysis of incompressible fluid flows in complex environments. The derived error estimator is applicable to geometric details constrained by either Dirichlet or Neumann boundary conditions, and has no special requirements on the outer boundary conditions. Extensive numerical examples were presented to demonstrate the effectiveness and efficiency of the proposed error estimator.