Numerical investigation of laminar Rayleigh-Benard convection of power-law fluids in square cross-sectional cylindrical annular enclosures

Steady-state laminar Rayleigh-Benard convection (i.e. differentially heated horizontal walls heated from below) of power-law fluids in square cross-sectional cylindrical annular enclosures has been numerically investigated under the assumption of axisymmetry. The numerical simulations have been conducted for a range of different values of nominal Rayleigh number Ra, nominal Prandtl number Pr, power-law index n, and internal radius to enclosure height ratio r(i)/L (i.e. 10(3)<= Ra <= 10(5); 10 <= Pr <= 10(4); 0.6 <= n <= 1.8; 0 <= r(i)/L <= 24) for both constant wall temperature (CWT) and constant wall heat flux (CWHF) boundary conditions for differentially heated horizontal walls. It has been found that convective transport is stronger for CWT boundary condition than for CWHF boundary condition for large (small) values of Ra (n) for a given set of values of n (Ra), Pr, and r(i)/L, but an opposite trend is observed for small (large) values of Ra (n). The mean Nusselt number (Nu) over bar (cy) does not show a monotonic increase with increasing (decreasing) Ra (n) especially for small values of ri/L for a given value of Pr due to changes in flow patterns (i.e. number of convection cells). However, the mean Nusselt number (Nu) over bar (cy) and flow patterns for large values of r(i)/L approach those for square enclosures (r(i)/L -> infinity) for both CWT and CWHF boundary conditions. Additionally, the critical Rayleigh number Rant for the onset of convection has been found to be mostly insensitive to the value of r(i)/L for both CWT and CWHF boundary conditions. (C) 2016 Elsevier Ltd. All rights reserved.