An adaptive remeshing procedure based on interpolation error has been developed for solving laminar natural convection problems. A simple relation between the Rayleigh number and the equilibration constant is introduced to increase the efficiency of the method. Unstructured meshes have been regenerated based on an initial solution and used in the calculations. This avoids the expensive mesh sensitivity study. Two typical natural convection problems are solved to demonstrate the present technique. Excellent performance has been observed at moderate and high Rayleigh numbers although present method gives good results over the whole range of Rayleigh numbers considered.