Influences of aspect ratio and wall boundary condition on laminar Rayleigh–Bénard convection of Bingham fluids in rectangular enclosures

• Published in: International journal of numerical methods for heat & fluid flow 2017-02, Vol.27 (2), p.310-333
• Main Authors: Yigit, Sahin, Chakraborty, Nilanjan
• Format: Article
• Language: English
• Subjects: Aspect ratio; Boundary conditions; Computational fluid dynamics; Computer simulation; Conservation equations; Convection; Cooling; Dependence; Enclosures; Energy conservation; Finite volume method; Flow resistance; Fluid flow; Fluids; Heat flux; Heat transfer; Height; Momentum; Prandtl number; Rayleigh number; Rheology; Scaling; Simulation; Temperature; Thermal conductivity; Transport; Viscosity; Wall temperature; Yield stress
• ISSN: 0961-5539; 1758-6585
• DOI: 10.1108/HFF-09-2015-0366

Purpose This paper aims to investigate the aspect ratio (AR; ratio of enclosure height:length) dependence of steady-state Rayleigh–Bénard convection of Bingham fluids within rectangular enclosures for both constant wall temperature and constant wall heat flux boundary conditions. A nominal Rayleigh number range 103 ≤ Ra ≤ 105 (Ra defined based on the height) for a single representative value of nominal Prandtl number (i.e. Pr = 500) has been considered for 1/4 ≤ AR ≤ 4. Design/methodology/approach The bi-viscosity Bingham model is used to mimic Bingham fluids for Rayleigh–Bénard convection of Bingham fluids in rectangular enclosures. The conservation equations of mass, momentum and energy have been solved in a coupled manner using the finite volume method where a second-order central differencing scheme is used for the diffusive terms and a second-order up-wind scheme is used for the convective terms. The well-known semi-implicit method for pressure-linked equations algorithm is used for the coupling of the pressure and velocity. Findings It has been found that buoyancy-driven flow strengthens with increasing nominal Rayleigh number Ra, but the convective transport weakens with increasing Bingham number Bn, because of additional flow resistance arising from yield stress in Bingham fluids. The relative contribution of thermal conduction (advection) to the total thermal transport strengthens (diminishes) with increasing AR for a given set of values of Ra and Pr for both Newtonian and Bingham fluids for both boundary conditions, and the thermal transport takes place purely because of conduction for tall enclosures. Originality/value Correlations for the mean Nusselt number Nu ¯ have been proposed for both boundary conditions for both Newtonian and Bingham fluids using scaling arguments, and the correlations have been demonstrated to predict Nu ¯ obtained from simulation data for 1/4 ≤ AR ≤ 4, 103 ≤ Ra ≤ 105 and Pr = 500.