The present study numerically investigates the motion of a solid body suspended in the square enclosure with natural convection. A two-dimensional circular cylinder levitated thermally has been simulated by using thermal lattice Boltzmann method(TLBM) with the direct-forcing immersed boundary method. To deal with the ascending, falling or levitation of a circular cylinder in natural convection, the immersed boundary method is expanded and coupled with the TLBM. The circular cylinder is located at the bottom of a square enclosure with no restriction on the motion and freely migrates due to the Boussinesq approximation which is employed for the coupling between the flow and temperature fields. For different density ratio between the cylinder and the fluid, the motion characteristics of the circular cylinder for various Grashof numbers have been carried out. The Prandtl number is fixed as 0.7.