Natural convection in a square cavity filled with a porous medium: Effects of various thermal boundary conditions

Natural convection flows in a square cavity filled with a porous matrix has been studied numerically using penalty finite element method for uniformly and non-uniformly heated bottom wall, and adiabatic top wall maintaining constant temperature of cold vertical walls. Darcy–Forchheimer model is used to simulate the momentum transfer in the porous medium. The numerical procedure is adopted in the present study yields consistent performance over a wide range of parameters (Rayleigh number Ra, 10³ ⩽ Ra ⩽ 10⁶, Darcy number Da, 10⁻⁵ ⩽ Da ⩽ 10⁻³, and Prandtl number Pr, 0.71 ⩽ Pr ⩽ 10) with respect to continuous and discontinuous thermal boundary conditions. Numerical results are presented in terms of stream functions, temperature profiles and Nusselt numbers. Non-uniform heating of the bottom wall produces greater heat transfer rate at the center of the bottom wall than uniform heating case for all Rayleigh numbers but average Nusselt number shows overall lower heat transfer rate for non-uniform heating case. It has been found that the heat transfer is primarily due to conduction for Da ⩽ 10⁻⁵ irrespective of Ra and Pr. The conductive heat transfer regime as a function of Ra has also been reported for Da ⩾ 10⁻⁴. Critical Rayleigh numbers for conduction dominant heat transfer cases have been obtained and for convection dominated regimes the power law correlations between average Nusselt number and Rayleigh numbers are presented.     

Effects of thermal boundary conditions on natural convection flows within a square cavity

A numerical study to investigate the steady laminar natural convection flow in a square cavity with uniformly and non-uniformly heated bottom wall, and adiabatic top wall maintaining constant temperature of cold vertical walls has been performed. A penalty finite element method with bi-quadratic rectangular elements has been used to solve the governing mass, momentum and energy equations. The numerical procedure adopted in the present study yields consistent performance over a wide range of parameters (Rayleigh number Ra, 10³ ⩽ Ra ⩽ 10⁵ and Prandtl number Pr, 0.7 ⩽ Pr ⩽ 10) with respect to continuous and discontinuous Dirichlet boundary conditions. Non-uniform heating of the bottom wall produces greater heat transfer rates at the center of the bottom wall than the uniform heating case for all Rayleigh numbers; however, average Nusselt numbers show overall lower heat transfer rates for the non-uniform heating case. Critical Rayleigh numbers for conduction dominant heat transfer cases have been obtained and for convection dominated regimes, power law correlations between average Nusselt number and Rayleigh numbers are presented.     

Thermal lattice-BGK model based on large-eddy simulation of turbulent natural convection due to internal heat generation

To simulate turbulent convection at high Rayleigh number (Ra), we propose a new thermal lattice-BGK (LBGK) model based on large eddy simulation (LES). Two-dimensional numerical simulations of natural convection with internal heat generation in a square cavity were performed at Ra from 10⁶ to 10¹³ with Prandtl numbers (Pr) at 0.25 and 0.60. Simulation results indicate that our model is fit to simulate high Ra flow for its better numerical stability. At Ra = 10¹³, a global turbulent has occurred. With a further increase in Ra, the flow will arrive in a fully turbulence regime. The Nusselt–Rayleigh relationship is also discussed.     

Steady natural convection flows in a square cavity with linearly heated side wall(s)

The present numerical study deals with natural convection flow in a closed square cavity when the bottom wall is uniformly heated and vertical wall(s) are linearly heated whereas the top wall is well insulated. Non-linear coupled PDEs governing the flow have been solved by penalty finite element method with bi-quadratic rectangular elements. Numerical results are obtained for various values of Rayleigh number (Ra) (10³ ⩽ Ra ⩽ 10⁵) and Prandtl number (Pr) (0.7 ⩽ Pr ⩽ 10). Results are presented in the form of streamlines, isotherm contours, local Nusselt number and the average Nusselt as a function of Rayleigh number.     

Heat flow visualization for natural convection in rhombic enclosures due to isothermal and non-isothermal heating at the bottom wall

Analysis has been carried out for the energy distribution and thermal mixing in steady laminar natural convective flow through the rhombic enclosures with various inclination angles, φ for various industrial applications. Simulations are carried out for various regimes of Prandtl (Pr) and Rayleigh (Ra) numbers. Dimensionless streamfunctions and heatfunctions are used to visualize the flow and energy distribution, respectively. Multiple flow circulations are observed at Pr = 0.015 and 0.7 for all φs at Ra = 10⁵. On the other hand, two asymmetric flow circulation cells are found to occupy the entire cavity for φ = 75° at higher Pr (Pr = 7.2 and 1000) and Ra (Ra = 10⁵). Heatlines are found to be parallel circular arcs connecting the cold and hot walls for the conduction dominant heat transfer at Ra = 10³. The enhanced convective heat transfer is explained with dense heatlines and convective loop of heatlines at Ra = 10⁵. Heatlines clearly demonstrate that the left wall receives heat from the bottom wall as heatlines directly connect both the walls whereas the convective heat circulation cells play lead role to distribute the heat along the right wall, especially for smaller φs. On the other hand, the heat flow is evenly distributed to both side walls at higher φs via convection as well as direct conductive transport. Significant convective heat transfer from the bottom hot wall to the left cold wall occurs for φ = 30° cavity whereas the heat transfer to the right cold wall is maximum for φ = 75° irrespective of Pr. Average Nusselt number studies also show that φ = 30° cavity gives maximum heat transfer rate from the bottom to left wall irrespective of Pr in isothermal heating case. On the other hand, enhanced thermal mixing occurs at φ = 75° for both isothermal and non-isothermal heating strategies except at Pr = 0.015 in isothermal heating case.     

Critical condition for Rayleigh-Bénard convection of Bingham fluids in rectangular enclosures

Two-dimensional steady-state numerical simulations have been conducted for laminar Rayleigh-Bénard convection of Bingham fluids in rectangular enclosures to analyse the critical Rayleigh number Ra_crit for which convection ceases to influence the thermal transport and thermal conduction becomes the principal heat transfer mechanism. The influences of Bingham number Bn on the critical Rayleigh number Ra_crit have been investigated for different values of aspect ratio (height: length) AR (ranging from 1/4 to 4) and nominal Prandtl number Pr (ranging from 10 to 500) for both constant wall temperature (CWT) and constant wall heat flux (CWHF) boundary conditions for the horizontal walls. It has been found that Ra_crit increases with increasing values of Bn and AR, regardless of the boundary condition. The values of Ra_crit have been found to be greater in the case of CWT boundary condition than in the CWHF configuration for AR ≤ 1, whereas an opposite trend is obtained for AR > 1 for Bingham fluids. Additionally, Ra_crit has been found be insensitive to the change of Pr for Newtonian fluids (i.e. Bn = 0), whereas Ra_crit increases with increasing Pr for Bingham fluids irrespective of the boundary condition. A detailed scaling analysis has also been performed to elucidate the effects of Bn ,Pr , AR on Ra_crit for Bingham fluids. The results of scaling analysis and numerical findings have been utilised to propose a new correlation for Ra_crit for both Newtonian and Bingham fluids in the case of both CWT and CWHF boundary conditions.     

Numerical investigation of several physical and geometric parameters in the natural convection into trapezoidal cavities

Natural convection in trapezoidal cavities, especially those with two internal baffles in conjunction with an insulated floor, inclined top surface, and isothermal left-heated and isothermal right-cooled vertical walls, has been investigated numerically using the Element based Finite Volume Method (EbFVM). In numerical simulations, the effect of three inclination angles of the upper surface as well as the effect of the Rayleigh number (Ra), the Prandtl number (Pr), and the baffle’s height (H_b) on the stream functions, temperature profiles, and local and average Nusselt numbers has been investigated. A parametric study was performed for a wide range of Ra numbers (10³ ? Ra ? 10⁶) H_b heights (H_b = H¹/3, 2H¹/3, and H¹), Pr numbers (Pr = 0.7, 10 and 130), and top angle (θ) ranges from 10 to 20. A correlation for the average Nusselt number in terms of Pr and Ra numbers, and the inclination of the upper surface of the cavity is proposed for each baffle height investigated.     

Investigation of natural circulation in cavities with uniform heat generation for different Prandtl number fluids

Natural convection in enclosures with uniform heat generation and isothermal side walls is studied here. For the rectangular enclosure, two-dimensional conservation equations are solved using SIMPLE algorithm. Parametric studies are conducted to examine the effects of orientation of the cavity, fluid properties (Pr number), and aspect ratio for Rayleigh numbers up to 10⁶. For a horizontal square cavity, the flow becomes periodically oscillating at Ra = 5 × 10⁴ and chaotic at Ra = 8 × 10⁵. With a slight increase in the inclination angle, the oscillations die and for inclination angles greater than 15⁰, the flow attain a steady state over a range of Ra. It is found that for tall cavities (aspect ratio > 1), the steady-state solution is obtained for all values of Ra considered here. However, for wide cavities (aspect ratio < 1), an oscillatory flow regime is observed. The maximum temperature within the cavity is calculated for the range of Ra, aspect ratio and Pr number. Correlations for the maximum cavity temperature is presented here. The values of critical Rayleigh number at which the convection sets in the rectangular cavity are also studied and two distinct criteria are determined to evaluate the critical Rayleigh number. Further, a three-dimensional simulation is performed for a cubic cavity. It is found that the steady state solutions are obtained for all Rayleigh number, except at Ra = 10⁶. This is in contrast to the predictions for a two-dimensional square cavity, which has an oscillatory zone from Ra = 5 × 10⁴ onwards.     

Natural convection heat transfer enhancement using adiabatic block: Optimal block size and Prandtl number effect

Numerical methods are used to solve the finite volume formulation of the two-dimensional mass, momentum and energy equations for steady-state natural convection inside a square enclosure. The enclosure consists of adiabatic horizontal walls and differentially heated vertical walls, but it also contains an adiabatic centrally-placed solid block. The aim of the study is to delineate the effect of such a block on the flow and temperature fields. The parametric study covers the range 10³ ⩽ Ra ⩽ 10⁶ and is done at three Pr namely, 0.071, 0.71 and 7.1. In addition the effect of increasing the size (characterized by the solidity Φ) of the adiabatic block is ascertained. It is found that the wall heat transfer increases, with increase in the Φ, until it reaches a critical value Φ = Φ_OPT, where the wall heat transfer attains its maximum. Further increases in the block size beyond Φ_OPT, reduces the wall heat transfer, for as the block size becomes larger than the conduction dominant core size it reduces the thermal mass of the convecting fluid. A steady-state heat transfer enhancement of 10% is observed for certain Ra and Pr values. Useful correlations predicting this optimum block size and the corresponding maximum heat transfer as a function of Ra and Pr are proposed; these predict within ±3%, the numerical results.     

Numerical research of nature convective heat transfer enhancement filled with nanofluids in rectangular enclosures

Heat transfer enhancement utilizing nanofluids in a two-dimensional enclosure is investigated for various pertinent parameters. The Khanafer's model is used to analyze heat transfer performance of nanofluids inside an enclosure taking into account the solid particle dispersion. Transport equations are model by a stream function-vorticity formulation and are solved numerically by finite-difference approach. Based upon the numerical predictions, the effects of Rayleigh number (Ra) and aspect ratio (AR) on the flow pattern and energy transport within the thermal boundary layer are presented. The diameter of the nanoparticle d_p is taken as 10 nm in nanofluids. The buoyancy parameter is 10³ ≤ Ra ≤ 10⁶ and aspect ratios (AR) of two-dimensional enclosure are 1/2, 1, 2. Results show that increasing the buoyancy parameter and volume fraction of nanofluids cause an increase in the average heat transfer coefficient. Finally, the empirical equation was built between average Nusselt number and volume fraction.     

Effect of optical properties on oscillatory hydromagnetic double-diffusive convection within semitransparent fluid

The effect of radiative heat transfer on the hydromagnetic double-diffusive convection in two-dimensional rectangular enclosure is studied numerically for fixed Prandtl, Rayleigh, and Lewis numbers, Pr = 13.6, Ra = 10⁵, Le = 2. Uniform temperatures and concentrations are imposed along the vertical walls while the horizontal walls are assumed to be adiabatic and impermeable to mass transfer. The influences of the optical thickness and scattering albedo of the semitransparent fluid on heat and mass transfer with and without magnetic damping are depicted. When progressively varying the optical thickness, multiple solutions are obtained which are steady or oscillatory accordingly to the initial conditions. the mechanisms of the transitions between steady compositionally dominated flow and unsteady thermally dominated flow are analyzed.     

Turbulent Rayleigh–Bénard convection of water in cubical cavities: A numerical and experimental study

Experimental measurements and numerical simulations of natural convection in a cubical cavity heated from below and cooled from above are reported at turbulent Rayleigh numbers using water as a convective fluid (Pr = 6.0). Direct numerical simulations were carried out considering the Boussinesq approximation with a second-order finite volume code (10⁷ ⩽ Ra ⩽ 10⁸). The particle image velocimetry technique was used to measure the velocity field at Ra = 10⁷, Ra = 7 × 10⁷ and Ra = 10⁸ and there was general agreement between the predicted time averaged local velocities and those experimentally measured if the heat conduction through the sidewalls was considered in the simulations.     

Bénard convection from a circular cylinder in a packed bed

Bénard convection around a circular heated cylinder embedded in a packed bed of spheres is studied numerically. The Forchheimer–Brinkman–extended Darcy momentum model with the Local Thermal Non-Equilibrium energy model is used in the mathematical formulation for the porous layer. The governing parameters considered are the Rayleigh number (10³ ≤ Ra ≤ 5 × 10⁷) and the thermal conductivity ratio (0.1 ≤ k_r ≤ 10,000). The structural properties of the packed bed are kept constant as: cylinder-to-particle diameter ratio D/d = 20 and porosity ε = 0.5, while the Prandtl number is fixed at Pr = 0.71. It is found that the presence of the porous medium suppresses significantly the strong free convection produced in the empty enclosure, and reduces considerably the high intensity of the pair of vortices generated behind the cylinder. Also, the results show that the porous medium can play the role of insulator or enhancer of heat transfer from the heat source, depending mainly on their thermal conductivities regardless of the Rayleigh number.     

Simulation of turbulent natural convection in a porous cylindrical annulus using a macroscopic two-equation model

This work presents numerical computations for laminar and turbulent natural convection within a horizontal cylindrical annulus filled with a fluid saturated porous medium. Computations covered the range 25 < Ra_m < 500 and 3.2 × 10⁻⁴ > Da > 3.2 × 10⁻⁶ and made use of the finite volume method. The inner and outer walls are maintained at constant but different temperatures. The macroscopic k–ε turbulence model with wall function is used to handle turbulent flows in porous media. First, the turbulence model is switched off and the laminar branch of the solution is found when increasing the Rayleigh number, Ra_m. Subsequently, the turbulence model is included and calculations start at high Ra_m, merging to the laminar branch for a reducing Ra_m. This convergence of results as Ra_m decreases can be seen as an estimate of the so-called laminarization phenomenon. Here, a critical Rayleigh number was not identified and results indicated that when the porosity, Prandtl number, conductivity ratio between the fluid and the solid matrix and Ra_m are kept fixed, the lower the Darcy number, the higher is the difference of the average Nusselt number given by the laminar and turbulent models.     

Steady natural convection flow in a square cavity filled with a porous medium for linearly heated side wall(s)

In this paper natural convection flows in a square cavity filled with a porous matrix has been investigated numerically when the bottom wall is uniformly heated and vertical wall(s) are linearly heated whereas the top wall is well insulated. Darcy–Forchheimer model without the inertia term is used to simulate the momentum transfer in the porous medium. Penalty finite element method with bi-quadratic rectangular elements is used to solve the non-dimensional governing equations. Numerical results are presented for a range of parameters (Rayleigh number Ra, 10³ ⩽ Ra ⩽ 10⁶, Darcy number Da, 10⁻⁵ ⩽ Da ⩽ 10⁻³, and Prandtl number Pr, 0.2 ⩽ Pr ⩽ 100) in terms of stream functions and isotherm contours, and local and average Nusselt numbers.     

Asymptotic analysis of heat transfer in turbulent Rayleigh–Bénard convection

Based on asymptotic considerations a heat transfer law for turbulent Rayleigh–Bénard convection is found that differs from existing correlations which often are of a power law type with respect to their Rayleigh number dependence. From the asymptotic temperature profile, derived by matching temperature gradients in the overlap region of the wall layer and the core layer, a Nusselt number follows which includes a logarithmic term. This correlation is in good agreement with data from highly accurate Rayleigh–Bénard experiments for Rayleigh numbers between 10⁵ and 10¹⁵ and Prandtl numbers larger than 0.5. It is an alternative to existing power laws or more complicated correlations for Nu = Nu(Ra,Pr).     

Natural convection in an annulus between coaxial horizontal cylinders with internal heat generation

Natural convection of gas (Pr = 0.7) between two horizontal coaxial cylinders with uniform internal heat generation is numerically investigated. Such a problem plays an important role in the analysis of lasers in which a similar circuit of a heat-conducting path is used. It has been found that the behavior of the system critically depends on three parameters including the inverse relative gap width σ(=diameter of the inner cylinder/gap width), the Rayleigh number Ra and the modified Rayleigh number Ra_T which describes heat generation. In the case Ra_T = 0 our results coincide with the known ones. It has been established that in such a system there exist two types of fluid flow for low Rayleigh numbers with different vortex structure. Optimization of the corresponding coaxial laser system has been analyzed.     

Laminar natural convection in a square cavity: Low Prandtl numbers and large density differences

Steady natural convection at low Prandtl numbers caused by large density differences in a square cavity heated through the side walls is investigated numerically and theoretically. An appropriate dimensionless parameter characterizing the density differences of the working fluid is identified by the Gay-Lussac number. The Boussinesq assumption is achieved when the Gay-Lussac number tends to zero. The Nusselt number is derived for the ranges in Rayleigh number 10 ? Ra ? 10⁸, in Prandtl number 0.0071 ? Pr ? 7.1 and in Gay-Lussac number 0 ? Ga < 2. The effects of the Rayleigh, Prandtl and Gay-Lussac numbers on the Nusselt number are discussed on physical grounds by means of a scale analysis. Finally, based on physical arguments, a heat transfer correlation is proposed, valid for all Prandtl and Gay-Lussac number ranges addressed.     

Natural convection in a rectangular enclosure containing an oval-shaped heat source and filled with Fe3O4/water nanofluid

This work concerns with the study of natural convection heat transfer in rectangular cavities with an inside oval-shaped heat source filled with Fe₃O₄/water nanofluid. The finite element method is employed to solve the governing equations for this problem. Average Nusselt numbers are presented for a wide range of Rayleigh number (10³ ≤ Ra ≤ 10⁵), volume fraction of nanoparticles (0 ≤ ϕ ≤ 14%), and four different size and shapes of the heat source. Depending on concentration of the nanoparticle, geometry of the heat source, and the value of Rayleigh number different behaviors are monitored for average Nusselt numbers. Configuration of the heat source dictates a significant change on the behavior of the average Nusselt number, while addition of the nanoparticles has a negative effect on the magnitude of Nusselt number for this problem.     

Natural convection and fluid flow in inclined enclosure with a corner heater

The phenomena of natural convection in an inclined square enclosure heated via corner heater have been studied numerically. Finite difference method is used for solving momentum and energy equations in the form of stream function–vorticity. One wall of the enclosure is isothermal but its temperature is colder than that of heaters while the remaining walls are adiabatic. The numerical procedure adopted in this analysis yields consistent performance over a wide range of parameters; Rayleigh number, Ra (10³ ? Ra ? 10⁶); Prandtl number, Pr (0.07 ? Pr ? 70); dimensionless lengths of heater in x and y directions (0.25 ? hx ? 0.75, 0.25 ? hy ? 0.75); and inclination angle, ? (0° ? ? ? 270°). It is observed that heat transfer is maximum or minimum depending on the inclination angle and depending on the length of the corner heaters. The effect of Prandtl number on mean Nusselt number is more significant for Pr < 1.