Routes to chaos of natural convection flows in vertical channels

The aim of the present study is the analysis of the transition to turbulence of natural convection flows between two infinite vertical plates. For the study of the problem, a number of Direct Numerical Simulations (DNSs) have been performed. The continuity, momentum and energy equations, cast under the Boussinesq assumption, are tackled numerically by means of a pseudospectral method, through which the three-dimensional domain is decomposed with Chebychev polynomials in the wall-normal direction and with Fourier modes in the wall-parallel directions. For low Rayleigh number values, the predictions of the flow regimes are consistent with the classical analytical results and linear stability analyses. In particular, the first bifurcation (Ra ≈ 5800) from the so-called laminar conduction regime to steady convection is correctly captured. By increasing the Rayleigh number beyond a second critical value (Ra ≈ 10200), the flow regime becomes chaotic. This transition to chaos is found to be related with the amplification of spanwise instabilities occurring at scales larger than the channel gap, H. The study of the return of the system from the chaotic regime to the laminar base flow reveals a phenomenon of hysteresis, i.e. the chaotic regime persists even at Ra-values lower than the second critical value. From a numerical point of view, the predicted flow regimes appear to be extremely sensitive to the domain size, grid resolution and perturbation amplitude. These aspects are shown to be of crucial importance for the prediction of the heat transfer performance, and, hence, should be taken into consideration when numerical methods are used for the simulation of real-world problems.