NATURAL-CONVECTION FLOWS IN AIR-FILLED, DIFFERENTIALLY HEATED CAVITIES WITH ADIABATIC HORIZONTAL WALLS Shihe Xin LIMSI-CNRS, Orsay Cedex, and De ´ partement de Physique, Universite ´ Paris Sud, Orsay Cedex, France Patrick Le Que ´re ´ LIMSI-CNRS, Orsay Cedex, France Two-dimensional natural convection in air-filled differentially heated cavities with adiabatic horizontal walls has been revisited by using stability analysis algorithms such as Newton’s iteration (steady-state solving), Arnoldi’s method, and the continuation method. We are particularly interested in computing Hopf bifurcation points characterizing the onset of time-dependent flows in these cavities. Aspect ratios of 1–7 have been investigated, and accurate critical points of several unstable modes which were not fully available in the literature have been provided for each cavity. The critical values agree well with the results available and confirm previous observations: for A=1–3 the onset of time-dependent flows is due to the detached flow structure near the exiting corners of vertical boundary layers; for A=4–8 the onset of time-dependent flows results from traveling waves in vertical boundary layers. Investigating mesh dependence showed that the critical values obtained possess four significants figures for A=1–3 and six for A=4–8. Steady-state solution by Newton’s iteration allowed us to discover multiple steady-state solutions at A about 3, which have never been reported before. The first branch, B 1 , of steady-state solutions exists up to a turning point where is born the second branch, B 2 . It ends at another turning point at lower Rayleigh number, and one observes the third branch, B 3 , for increasing Rayleigh number. Solutions on B 1 do not exhibit any detached flow struc- ture, but those on B 2 and B 3 do. The detached flow structure is more pronounced on B 3 than on B 2 . For A=3, B 3 becomes unstable because of the detached flow structure. With increas- ing A, branch B 1 becomes unstable to traveling waves of vertical boundary layers and the stable part of B 3 decreases and disappears completely. Due to different instability mechan- isms observed on branches B 1 and B 3 , there exist multiple time-dependent flows of quite dif- ferent frequencies (for A=3.4 at Ra ¼1.15 10 8 , for example), which is also a novel phenomenon. Multiple steady-state solutions disappear between A=2.8 and 2.9 and are believed to be closely linked to the behavior of the detached flow and the limiting effect of increasing A. Received 29 July 2005; accepted 17 December 2005. The authors thank Dr. L. S. Tuckerman for fruitful discussions and suggestions. Computations have been performed at CINES under research project lim2072 and at IDRIS under research project 50326. Address correspondence to Shihe Xin, LIMSI-CNRS, BP 133, 91403 Orsay Cedex, France. E-mail: xin@limsi.fr 437 Numerical Heat Transfer, Part A, 50: 437–466, 2006 Copyright # Taylor & Francis Group, LLC ISSN: 1040-7782 print=1521-0634 online DOI: 10.1080/10407780600605039