Numerical simulations are performed for laminar natural convection heat transfer from a centered conducting body enclosed in a square cavity. A high accuracy unsteady numerical method is used, combining the unique condition of the pressure, the convergent solutions and the stream-function value of the inside heat-conducting body are given simultaneously. Two examples are simulated with this numerical method and compared with the experimental results. The results of the numerical solutions are consistent with the experimental results. It shows that the numerical method is valid and feasible.