In this work, we consider the heat equation coupled with Stokes equations under threshold type boundary condition. The conditions for existence and uniqueness of the weak solution are made clear. Next we formulate the finite element problem, recall the conditions of its solvability, and study its convergence by making use of Babuska–Brezzi's conditions for mixed problems. Third we formulate an Uzawa's type iterative algorithm that separates the fluid from heat conduction, study its feasibility, and convergence. Finally the theoretical findings are validated by numerical simulations.