Development of unsteady natural convection in a square cavity under large temperature difference

Abstract

To investigate how the nonuniform fluid density distribution caused by large temperature variations affects the development of unsteady natural convection, we perform a series of direct numerical simulations of two-dimensional compressible natural convection in an air-filled square cavity. The cavity has a hot wall on the left and a cold wall on the right, and two horizontal walls are adiabatic. The simulations are done using a kinetic approach based on a modeled Boltzmann equation, from which the fully compressible Navier–Stokes–Fourier equations are recovered. No Boussinesq approximation or low Mach number approximation is made. An extra source term is introduced to adjust the fluid Prandtl number. Simulations are performed for a range of Rayleigh numbers (107−109) with a fixed dimensionless temperature difference of ε=0.6 to determine the critical Rayleigh number and study the development of unsteady flow. To illustrate the instability mechanism, instantaneous fluctuation field, time trace of temperature, and velocity at selected monitoring points, the spectrum and other statistics are presented and discussed. As expected, significant differences are observed between the instability of compressible natural convection and the Boussinesq-type natural convection. With a large temperature difference, the transition to unsteady flow is asymmetric for the flows near the hot wall and cold wall. For the Rayleigh number range we studied, the cold wall region is dominated by low-frequency impact instability of the boundary thermal jet at the bottom corner. For the hot wall region, besides the upper corner impact instability, a boundary layer instability featuring high-frequency oscillations is observed.