A computer code for enclosed natural convection in the boundary layer regime

A computer code is developed for the numerical prediction of natural convection in rectangular two-dimensional cavities at high Rayleigh numbers. The governing equations are retained in the primitive variable form. The numerical method is based on finite differences and an ADI scheme. Convective terms may be approximated with either central or hybrid differencing for greater stability. A non-uniform grid distribution is possible for greater efficiency. The pressure is dealt with via a SIMPLE type algorithm and the use of a fast elliptic solver for the solenoidal velocity correction field significantly reduces computing times. Preliminary results indicate that the code is reasonably accurate, robust and fast compared with existing benchmarks and finite difference based codes, particularly at high Rayleigh numbers. Extension to three-dimensional problems and turbulence studies in similar geometries is readily possible and indicated.

The laminar boundary layer regime for natural convection of air in a square cavity

Numerical predictions are obtained for laminar natural convection of air in a square two dimensional cavity at high Rayleigh numbers. Proper resolution of the core reveals weak multi-cellular structure which varies in a complex manner as the effects of convection are increased. The end of the steady laminar regime is numerically estimated to occur at Ra=2.2x10^8.