Flow reversals in Rayleigh–Bénard convection with non-Oberbeck–Boussinesq effects

Flow reversals in two-dimensional Rayleigh–Benard convection led by non-Oberbeck–Boussinesq (NOB) effects due to large temperature differences are studied by direct numerical simulation. Perfect gas is chosen as the working fluid and the Prandtl number is 0.71 for the reference state. If NOB effects are included, the flow pattern $P_{11}$ with only one dominant roll often becomes unstable by the growth of the cold corner roll, which sometimes results in cession-led flow reversals. By exploiting the vorticity transport equation, it is found that the asymmetries of buoyancy and viscous forces are responsible for the growth of the cold corner roll because both such asymmetries cause an imbalance between the corner rolls and the large-scale circulation (LSC). The buoyancy force near the cold wall increases and decreases near the hot wall originating from the temperature-dependent isobaric thermal expansion coefficient ${\it\alpha}=1/T$ if NOB effects are included. Moreover, the decreased dissipation due to lower viscosity is favourable for the growth of the cold corner roll, while the increased viscosity further suppresses the growth of the hot corner roll. Finally, it is found that the boundary layer near the cold wall plays an important role in the mass transport from LSC to corner rolls subject to mass conservation.