Steady-state two-dimensional natural convective heat transfer in horizontal cylindrical annuli was studied by solving the governing equations based on the primitive variables. Emphasis was put on the occurrence of the multiple solutions at a given set of parameter values, and on the determination of the bifurcation points at which those multiple solutions begin to branch out. The multicellular flow pattern from the results of melting process in an isothermally heated horizontal cylinder for high Rayleigh numbers, was used as initial guesses for the field variables. This was succeeded in new bifurcation point to tetracellular solutions for an identical set of parameter variables of previous works. The close examination of flow pattern transition around bifurcation point was also conducted. It was found that the mechanisms of flow transition are different depending on the critical Rayleigh number of bifurcation point.