The computation of natural convection in the annulus between horizontal concentric cylinders using spectral element method in polar coordinates is presented in this paper.