Physical mechanism of the convective heat flux increasing in case of mixed boundary conditions in Rayleigh-Bénard convection

Author(s):  
A. Sukhanovskii ◽  
A. Vasiliev
Fluids ◽  
2020 ◽  
Vol 5 (4) ◽  
pp. 185
Author(s):  
Patrick Fischer ◽  
Charles-Henri Bruneau ◽  
Hamid Kellay

Numerical simulations of rotating two-dimensional turbulent thermal convection on a hemisphere are presented in this paper. Previous experiments on a half soap bubble located on a heated plate have been used for studying thermal convection as well as the effects of rotation on a curved surface. Here, two different methods have been used to produce the rotation of the hemisphere: the classical rotation term added to the velocity equation, and a non-zero azimuthal velocity boundary condition. This latter method is more adapted to the soap bubble experiments. These two methods of forcing the rotation of the hemisphere induce different fluid dynamics. While the first method is classically used for describing rotating Rayleigh–Bénard convection experiments, the second method seems to be more adapted for describing rotating flows where a shear layer may be dominant. This is particularly the case where the fluid is not contained in a closed container and the rotation is imposed on only one side of it. Four different diagnostics have been used to compare the two methods: the Nusselt number, the effective computation of the convective heat flux, the velocity and temperature fluctuations root mean square (RMS) generation of vertically aligned vortex tubes (to evaluate the boundary layers) and the energy/enstrophy/temperature spectra/fluxes. We observe that the dynamics of the convective heat flux is strongly inhibited by high rotations for the two different forcing methods. Also, and contrary to classical three-dimensional rotating Rayleigh–Bénard convection experiments, almost no significant improvement of the convective heat flux has been observed when adding a rotation term in the velocity equation. However, moderate rotations induced by non-zero velocity boundary conditions induce a significant enhancement of the convective heat flux. This enhancement is closely related to the presence of a shear layer and to the thermal boundary layer just above the equator.


Author(s):  
Thomas B. Gradinger ◽  
T. Laneryd

Natural-convection cooling with oil or other fluids of high Prandtl number plays an important role in many technical applications such as transformers or other electric equipment. For design and optimization, one-dimensional (1D) flow models are of great value. A standard configuration in such models is flow between vertical parallel plates. Accurate modeling of heat transfer, buoyancy, and pressure drop for this configuration is therefore of high importance but gets challenging as the influence of buoyancy rises. For increasing ratio of Grashof to Reynolds number, the accuracy of one-dimensional models based on the locally forced-flow assumption drops. In the present work, buoyancy corrections for use in one-dimensional models are developed and verified. Based on two-dimensional (2D) simulations of buoyant flow using finite-element solver COMSOL Multiphysics, corrections are derived for the local Nusselt number, the local friction coefficient, and a parameter relating velocity-weighted and volumetric mean temperature. The corrections are expressed in terms of the ratio of local Grashof to Reynolds number and a normalized distance from the channel inlet, both readily available in a one-dimensional model. The corrections universally apply to constant wall temperature, constant wall heat flux, and mixed boundary conditions. The developed correlations are tested against two-dimensional simulations for a case of mixed boundary conditions and are found to yield high accuracy in temperature, wall heat flux, and wall shear stress. An application example of a natural-convection loop with two finned heat exchangers shows the influence on mass-flow rate and top-to-bottom temperature difference.


2017 ◽  
Vol 835 ◽  
pp. 491-511 ◽  
Author(s):  
Dennis Bakhuis ◽  
Rodolfo Ostilla-Mónico ◽  
Erwin P. van der Poel ◽  
Roberto Verzicco ◽  
Detlef Lohse

A series of direct numerical simulations of Rayleigh–Bénard convection, the flow in a fluid layer heated from below and cooled from above, were conducted to investigate the effect of mixed insulating and conducting boundary conditions on convective flows. Rayleigh numbers between $Ra=10^{7}$ and $Ra=10^{9}$ were considered, for Prandtl numbers $\mathit{Pr}=1$ and $\mathit{Pr}=10$. The bottom plate was divided into patterns of conducting and insulating stripes. The size ratio between these stripes was fixed to unity and the total number of stripes was varied. Global quantities, such as the heat transport and average bulk temperature, and local quantities, such as the temperature just below the insulating boundary wall, were investigated. For the case with the top boundary divided into two halves, one conducting and one insulating, the heat transfer was found to be approximately two-thirds of that for the fully conducting case. Increasing the pattern frequency increased the heat transfer, which asymptotically approached the fully conducting case, even if only half of the surface is conducting. Fourier analysis of the temperature field revealed that the imprinted pattern of the plates is diffused in the thermal boundary layers, and cannot be detected in the bulk. With conducting–insulating patterns on both plates, the trends previously described were similar; however, the half-and-half division led to a heat transfer of about a half of that for the fully conducting case instead of two-thirds. The effect of the ratio of conducting and insulating areas was also analysed, and it was found that, even for systems with a top plate with only 25 % conducting surface, heat transport of 60 % of the fully conducting case can be seen. Changing the one-dimensional stripe pattern to a two-dimensional chequerboard tessellation does not result in a significantly different response of the system.


2017 ◽  
Vol 822 ◽  
pp. 1-4 ◽  
Author(s):  
R. Ostilla-Mónico

Natural convection is omnipresent on Earth. A basic and well-studied model for it is Rayleigh–Bénard convection, the fluid flow in a layer heated from below and cooled from above. Most explorations of Rayleigh–Bénard convection focus on spatially uniform, perfectly conducting thermal boundary conditions, but many important geophysical phenomena are characterized by boundary conditions which are a mixture of conducting and adiabatic materials. For example, the differences in thermal conductivity between continental and oceanic lithospheres are believed to play an important role in plate tectonics. To study this, Wang et al. (J. Fluid Mech., vol. 817, 2017, R1), measure the effect of mixed adiabatic–conducting boundary conditions on turbulent Rayleigh–Bénard convection, finding experimental proof that even if the total heat transfer is primarily affected by the adiabatic fraction, the arrangement of adiabatic and conducting plates is crucial in determining the large-scale flow dynamics.


Sign in / Sign up

Export Citation Format

Share Document