The present study investigates numerically the three-dimensional mixed convection flow behavior of a Bingham plastic fluid within a rectangular horizontal duct. The latter has a square cross section, which is symmetrically heated at uniform wall temperature. The effect of the viscous dissipation on the hydrodynamic and thermal fields is taking into account through the energy equation. The finite volume method is used to discretize the governing equations, and the resulting algebraic system is solved iteratively by TDMA algorithm. The simulations are conducted for different values of the Grashof number (10 4^| Gr |^5 9 10 5), the Brinkman number (0^| Br |^10), and the Bingham number (10^| Bn |^20), for Re = 100 and Pr = 50. The results are analyzed through the distribution of the Nusselt number, the temperature and velocity profiles as well as unyielded plugs. For Br = 0, the Grashof number variation affects greatly the heat transfer and the fluid flow, in particular at the duct intermediate region. However, for Br = 0, the Grashof number variation effect is noticeable in both the central and the fully developed regions. Moreover, the combination of Grashof, Brinkman, and Bingham numbers constitutes pertinent parameters control of the unyielded plugs.