The effect of buoyancy on the flow and heat transfer that develop between a horizontal cold surface and an infinite two-dimensional array of open cavities heated from below is studied numerically. In earlier investigations the steady-state features of this problem were studied for the case of unbounded flow above the cavities. The resulting flow pattern was found to be symmetrical with respect to the centerlines of the cavities. In the present work it is shown that the symmetry of the flow can be destroyed due to the presence of an upper wall. The evolutionary path to steady-state flow is examined, and sustained oscillatory behavior has been observed in several cases. The solution structure is governed by five parameters, i.e., the geometric parameters A = 1'/H', B = h'/H', and C = L'/H', the Rayleigh number Ra = g-beta-DELTA-T'H'3/alpha-nu, and the Prandtl number Pr = nu/alpha. For a geometry with A = 1/2, B = 1/4, and C = 1, a complicated solution structure is observed upon increasing the Rayleigh number. For Ra less-than-or-equal-to 4 x 10(3), a steady symmetric two-cell pattern is observed. This pattern becomes asymmetric for 4 x 10(3) < Ra less-than-or-equal-to 9 x 10(4), periodic for Ra less-than-or-equal-to 3 x 10(5), and chaotic above that. The transition to periodic convection occurs at lower Rayleigh numbers with decreasing B.