In this paper, a numerical study on centrifugal instability with convective heat transfer through a curved square duct is presented by using a spectral method, and covering a wide range of the Dean number ( ) 0 5000 Dn for a tightly coiled square duct of curvature 0.5 The outer and bottom walls of the duct are heated while cooling from the inner and the ceiling. The main objective of this study is to expose combined effects of centrifugal and buoyancy forces on fluid flows through a curved channel. For this purpose, solution structure of the steady solutions is obtained first. As a result, four branches of symmetric/asymmetric steady solutions are obtained. Linear stability of the steady solutions is then investigated. It is found that only the first branch is linearly stable while the other branches are linearly unstable. Unsteady flow behavior, obtained by time evolution calculations, shows that the steady-state flow turns into chaotic flow via periodic and multiperiodic flows, if Dn is increased. Typical contours of secondary flow pattern, stream-wise velocity distribution and temperature profiles are obtained at several values of Dn and it is found that the flow consists of asymmetric two-to four-vortex solutions. The present study shows that convective heat transfer is significantly enhanced by the secondary flow; and the chaotic flow, which occurs at large Dn's, enhances heat transfer more effectively than the steady-state or periodic solutions. Finally, a comparison between the numerical and experimental investigations has been made, and it is found that there is a good agreement between the numerical and experimental investigations.