The effect of boundary conditions on the onset of chaos in Rayleigh–Bénard convection using energy-conserving Lorenz models

Abstract

The influence of 16 boundary conditions on linear and nonlinear stability analyses of Rayleigh–Bénard system is reported. A Stuart–Landau amplitude equation for the Rayleigh–Bénard system between stress-free, isothermal boundary conditions is derived and the procedure used in this derivation serves as guidance for constructing an appropriate Fourier–Galerkin expansion for the other 15 boundary conditions to derive a generalized Lorenz model. The influence of the boundary conditions comes within the coefficients of the generalized Lorenz model. It is shown that the obtained generalized Lorenz model is energy conserving and for certain boundary conditions it retains features of the classical Lorenz model. Further, the principle of exchange of stabilities is shown to be valid for the present problem and hence it is the steady-state, linearized version of the generalized Lorenz model which yields an analytical expression for the Rayleigh number. On minimizing this expression with respect to wave number the critical Rayleigh number at which the onset of regular convective motion occurs in the form of rolls is determined for all 16 boundary conditions. It is found that these values are in good agreement with those of previous investigations leading to the conclusion that the chosen minimal Fourier–Galerkin expansion is a valid one. Exhibition of chaotic motion in the generalized Lorenz system at the Hopf Rayleigh number is studied. The phase-space plots which indicate a clear-cut visualization of the transition from regular convective motion to chaotic motion in the generalized Lorenz system are presented. Further, existence of a developing region for chaos (mildly chaotic motion) and windows of periodicity are captured using the bifurcation diagrams. It is concluded from the phase-space plots and the bifurcation diagrams that the generalized Lorenz model for certain sets of boundary conditions retains all the features of the classical Lorenz model. Such a conclusion cannot be made by a linear stability analysis and the need thus for a nonlinear analysis is highlighted in the paper.