In this paper, three-dimensional numerical simulations are carried out to investigate and analyze the gradual effects of thermal buoyancy strength on laminar flow of an incompressible viscous fluid and heat transfer rate inside a 180° curved channel of square cross-section. The governing equations of continuity, momentum and energy balance are obtained and solved numerically using finite volume method. The effect of Dean number, De, and Richardson number, Ri, on dimensionless velocity profiles and Nusselt number are examined for the conditions: De = 125 to 150, Ri = 0 to 2 at Pr = 1. The mean results are illustrated in terms of streamline and isotherm contours to interpret the flow behaviors and its effect on heat transfer rate. Dimensionless velocity profiles and the local Nusselt number at the angle 0° and 90° are presented and discussed. Also, the average Nusselt number on surfaces of curved duct is computed. The obtained results showed that by adding thermal buoyancy to computed domain, some early Dean vortices are observed at the angle 0° and new sort are observed at 90°. Furthermore, increase in Dean number increases the heat transfer rate. In other hand, increase in Richardson number decreases the average Nusselt number of 180° curved duct.