Benard-Rayleigh convection in thin ring layer

A. Fedotov1, Y. Tsitavets1, Y. Levy2
1Belarusian State University, Minsk, Belarus
2HiLASE Centre, Prague, Czech Republic
Veröffentlicht in 2020

Convective instability in a thin liquid layer is observed in many technical systems [1] and is of interest to the self-organization theory. To date, however, only simple problems with the convection between flat plates are well studied [2]. In real systems and devices, convective flows often happen in curved or cylindrical layers. The study of these cases will contribute to practical applications as well as basic theory. In this work, the Benard-Rayleigh convection is addressed in 2 dimensions and in a thin ring layer with gravity perpendicular to the revolution axis. Rayleigh's numbers were kept in the range (1,5 - 4,4)∙10^5. For air, these values are high enough to lead to a quick development of periodic convective cells and low enough to neglect the turbulent subgrid transport. The inner wall of the ring was overheated by ∆T in comparison to the outer wall. COMSOL Multiphysics® provides a convenient interface for the finite-element computations suitable for solving fluid dynamics equations with heat transfer in Boussinesq approximation. For this purpose, preset equations in COMSOL Multiphysics® CFD Module and Heat Transfer Module can be used. Moreover, there is a possibility to track the movement of control fluid volumes with the idealized particles which velocity coincides with the velocity of fluid computed from Navier-Stokes equations. This allows the construction of phase portraits that provide valuable information on dynamic systems. Such functionality is implemented in the Particle Tracing Module. The time integration of partial differential equations was performed with the 3rd order backward differentiation method and the solving of linear equations was carried out by PARDISO algorithm. Computational meshes included 10^4 - 10^5 elements condensing in the vicinity of walls. It was established that convection in a ring layer can be separated by the character of the flow into four spatial regions: upper, transitional, lateral, and bottom. In the upper region classical Benard-Rayleigh cells are observed, in the transitional region the neighboring convective vortexes in cells become asymmetrical, in the lateral area there is a flow without the formation of stable structures. In the bottom region, there is only conductive heat transfer as the heated boundary is on top. In the upper part of the ring layer, the phase portrait is similar to that observed in classic planar Benard-Rayleigh convection. In general, the upper, transitional and bottom regions of the ring layer are characterized by stable phase portraits with periodic envelopes. As one moves along the layer, the phase portrait becomes more chaotic with the most unstable behavior in the lateral region.

References 1. Lappa, M. Rotating Thermal Flows in Natural and Industrial Processes / M. Lappa. – 2012. 2. Borońska, K. and Tuckerman, Laurette S. Extreme multiplicity in cylindrical Rayleigh-Benard convection. II. Bifurcation diagram and symmetry classification, Phys. Rev. E 81, 036321 (2010)

Herunterladen