The onset of zonal modes in two-dimensional Rayleigh–Bénard convection

Abstract

We study the stability of steady convection rolls in two-dimensional Rayleigh–Bénard convection with free-slip boundaries and horizontal periodicity over 12 orders of magnitude in the Prandtl number $(10^{-6} \leq Pr \leq 10^6)$ and 6 orders of magnitude in the Rayleigh number $(8{\rm \pi} ^4 < Ra \leq 10^8)$ . The analysis is facilitated by partitioning our modal expansion into so-called even and odd modes. With aspect ratio $\varGamma = 2$ , we observe that zonal modes (with horizontal wavenumber equal to zero) can emerge only once the steady convection roll state consisting of even modes only becomes unstable to odd perturbations. We determine the stability boundary in the $(Pr,Ra)$ plane and observe remarkably intricate features corresponding to qualitative changes in the solution, as well as three regions where the steady convection rolls lose and subsequently regain stability as the Rayleigh number is increased. We study the asymptotic limit $Pr \to 0$ and find that the steady convection rolls become unstable almost instantaneously, eventually leading to nonlinear relaxation osculations and bursts, which we can explain with a weakly nonlinear analysis. In the complementary large- $Pr$ limit, we observe that the zonal modes at the instability switch off abruptly at a large, but finite, Prandtl number.