The paper presents a computational method used to obtain solutions of the buoyancy-driven laminar and turbulent flow and heat transfer in a square cavity with differentially heated side walls. A series of Rayleigh numbers, rangingfrom 10 3 to 10 16 was studied. Donor-cell differencing is used, and mesh-refinement studies have been performed for all Rayleigh numbers considered. The turbulence model used for Rayleigh numbers greater than 10 6 is a ( k ~ ε) two-equation model of turbulence, that includes gravity ~ density gradient interactions. The results are presented in tabular and graphical form, and as correlations of the Nusselt and Rayleigh numbers. Furthermore, the results for Rayleigh numbers up to 10 6 are compared with the benchmark numerical solution of de Vahl Davis.