Numerical simulation of cellular convection in air

Three-dimensional convection in a Boussinesq fluid confined between horizontal rigid boundaries is studied in a series of numerical experiments. Convection in air, whose Prandtl number Pr = 0·71, is systematically investigated, together with another model for Pr = 1. Convection with a steadily changing mean temperature is also considered. Two-dimensional rolls over the Rayleigh number range 4500 [les ] Ra [les ] 24000 and three-dimensional flow patterns over the range 26000 [les ] Ra ≤ 32000 are shown to be stable in air when the mean temperature of the layer is constant (). Discrete changes in the slope of the heat-flux curve are shown to exist in the ranges \[ 7000\leqslant Ra\leqslant 8000,\quad 12000\leqslant Ra\leqslant 14000\quad{\rm and}\quad 24000\leqslant Ra\leqslant 26000 \] in air. Only the last discrete transition in the heat flux is asSociated with a significant transition in the flow pattern. Two-dimensional rolls with a horizontally asymmetric distribution of upward and downward motions over the range 4500 [les ] Ra [les ] 8000, and three-dimensional flow patterns over the range 10 000 [les ] Ra [les ] 20 000 are shown to be stable when the mean temperature varies with time. The circulation in a three-dimensional cell depends on the sign of the mean temperature change: downward motions occupy the centre of the cell when . Motions start to be time dependent for Ra > 20000. Transitions in the planform are asSociated with discrete changes in the slope of the heat-flux curve. Transitions in both the heat flux and flow pattern depend quantitatively on the Prandtl number.