A pseudo-spectral solver with multigrid acceleration for the numerical prediction of incompressible non-isothermal flows is presented. The spatial discretization is based on a Chebyshev collocation method on Gauss-Lobatto points and for the discretization in time the second-order backward differencing scheme (BDF2) is employed. The multigrid method is invoked at the level of algebraic system solving within a pressure-correction method. The approach combines the high accuracy of spectral methods with efficient solver properties of multigrid methods. The capabilities of the proposed scheme are illustrated by a buoyancy driven cavity flow as a standard benchmark case. To cite this article: K.