Direct numerical simulation of instabilities in a two-dimensional near-critical fluid layer heated from below

Abstract

An analysis of the hydrodynamic stability of a fluid near its near critical point – initially at rest and in thermodynamic equilibrium – is considered in the Rayleigh–Bénard configuration, i.e. heated from below. The geometry is a two-dimensional square cavity and the top and bottom walls are maintained at constant temperatures while the sidewalls are insulated. Owing to the homogeneous thermo-acoustic heating (piston effect), the thermal field exhibits a very specific structure in the vertical direction. A very thin hot thermal boundary layer is formed at the bottom, then a homogeneously heated bulk settles in the core at a lower temperature; at the top, a cooler boundary layer forms in order to continuously match the bulk temperature with the colder temperature of the upper wall. We analyse the stability of the two boundary layers by numerically solving the Navier–Stokes equations appropriate for a van der Waals' gas slightly above its critical point. A finite-volume method is used together with an acoustic filtering procedure. The onset of the instabilities in the two different layers is discussed with respect to the results of the theoretical stability analyses available in the literature and stability diagrams are derived. By accounting for the piston effect the present results can be put within the framework of the stability analysis of Gitterman and Steinberg for a single layer subjected to a uniform, steady temperature gradient.