LINEAR STABILITY ANALYSIS AND CFD SIMULATION OF DOUBLE-LAYER RAYLEIGHBÉNARD CONVECTION

Natural convection in superimposed layers of fluids heated from below is commonly observed inmany industrial and natural situations, such as crystal growth, co-extrusion processes and atmospheric flow.The stability analysis of this system reveals a complex dynamic behavior, including the potential multiplicityof stationary states and occurrence of periodic regimes. In this study, a linear stability analysis (LSA) wasperformed to determine the onset of natural convection as a function of imposed boundary conditions,geometrical configuration and specific perturbations. To investigate the effects of the non-linear terms neglectedin LSA, a direct simulation of the full nonlinear problem was performed using computational fluid dynamics(CFD) techniques. The numerical simulation results show an excellent agreement with the LSA results near theonset of convection and an increase in the deviation as the Rayleigh number increases above the critical value.