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.     

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.     

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.     

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.     

Entropy generation in natural convection in a symmetrically and uniformly heated vertical channel

In this study numerical predictions of local and global entropy generation rates in natural convection in air in a vertical channel symmetrically heated at uniform heat flux are reported. Results of entropy generation analysis are obtained by solving the entropy generation equation based on the velocity and temperature data. The analyzed regime is two-dimensional, laminar and steady state. The numerical procedure expands an existing computer code on natural convection in vertical channels. Results in terms of fields and profiles of local entropy generation, for various Rayleigh number, Ra, and aspect ratio values, L/b, are given. The distributions of local values show different behaviours for the different Ra values. A correlation between global entropy generation rates, Rayleigh number and aspect ratio is proposed in the ranges 10³ ⩽ Ra ⩽ 10⁶ and 5 ⩽ L/b ⩽ 20.     

Steady forced convection heat transfer from a heated circular cylinder to power-law fluids

Forced convection heat transfer to incompressible power-law fluids from a heated circular cylinder in the steady cross-flow regime has been investigated numerically by solving the momentum and thermal energy equations using a finite volume method and the QUICK scheme on a non-uniform Cartesian grid. The dependence of the average Nusselt number on the Reynolds number (5 ⩽ Re ⩽ 40), power-law index (0.6 ⩽ n ⩽ 2) and Prandtl number (1 ⩽ Pr ⩽ 1000) has been studied in detail. The numerical results are used to develop simple correlations as functions of the pertinent dimensionless variables. In addition to the average Nusselt number, the effects of Re, Pr and n on the local Nusselt number distribution have also been studied to provide further physical insights. The role of the two types of thermal boundary conditions, namely, constant temperature and uniform heat flux on the surface of the cylinder has also been presented.     

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.     

Natural convection in a nanofluid-filled eccentric annulus with constant heat flux wall: A lattice Boltzmann study with immersed boundary method

A numerical investigation of natural convection in a Cu–water nanofluid-filled eccentric annulus with constant heat flux wall is presented. The governing equations of the flow and temperature fields are solved by lattice Boltzmann method (LBM), and the Dirichlet and Neumann boundary conditions are treated using the immersed boundary method (IBM). Influences of the Rayleigh number (10³ ≤ Ra ≤ 10⁷), eccentricity (ε = −0.625,0 and 0.625), nanoparticles volume fraction (0 ≤ ϕ ≤ 0.03) and radial ratio (rr = 2.33,2.6 and 3) on the streamlines, isotherms and Nusselt number are studied. It is found that the inclusion of the nanoparticles into pure fluid changes the flow pattern. And the Nusselt number has a positive relationship with nanoparticle volume fraction, Rayleigh number and radial ratio. Also, it can be confirmed that Nusselt number in the case with negative eccentricity (ε = −0.625) is larger than the others.     

Oscillatory double-diffusive mixed convection in a two-dimensional ventilated enclosure

By starting from a steady flow configuration based on the work of Deng et al. [Qi-Hong Deng, Jiemin Zhou, Chi Mei, Yong-Ming Shen, Fluid, heat and contaminant transport structures of laminar double-diffusive mixed convection in a two-dimensional ventilated enclosure, Int. J. Heat Mass Transfer 47 (2004) 5257–5269], a numerical investigation was conducted to analyse the unsteady double-diffusive mixed convection in two-dimensional ventilated room due to heat and contaminant sources. Owing to the large number of parameters, the results are reported only for a constant buoyancy ratio N equal to 1. The flow is found to be oscillatory for a fixed Reynolds number (700 ⩽ Re ⩽ 1000) when the Grashof number is varied in a wide range (10³ ⩽ Gr ⩽ 10⁶). Results of the simulations show that the onset of the oscillatory indoor airflow occurs for couples (Re, Gr) values that can be correlated as Re = aGr^b.     

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.     

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.     

Jet impingement interaction with cross flow in horizontal porous layer under thermal non-equilibrium conditions

In the present article the jet impingement cooling of heated portion of a horizontal surface immersed in a thermally non-equilibrium porous layer is considered for investigation numerically with the presence of a cross flow. The mathematical model is derived for steady, two-dimensional laminar flow based on Darcy model and two-energy equation for fluid and solid phases. A parametric study is carried out by varying the following parameters: cross flow to jet flow velocity ratio parameter (0 ⩽ M ⩽ 1); porosity scaled thermal conductivity ratio parameter (0.1 ⩽ K_r ⩽ 1000); heat transfer coefficient parameter (0.1 ⩽ H ⩽ 1000); Péclet number (1 ⩽ Pe ⩽ 1000) and Rayleigh number (10 ⩽ Ra ⩽ 100). The total average Nusselt number is defined based on the overall thermal conductivity, which is assumed to be the arithmetic mean of the porosity scaled thermal conductivity of the fluid and solid phases. The total average Nusselt number as well as the average Nusselt number for both fluid and solid phases is presented for different governing parameters. It is found that the presence of a weak cross flow in a jet impinging jet may degrade the heat transfer. The results show that the average Nusselt number calculated from the thermal equilibrium model are the maximum possible values and these values can be reproduced by large values of H × K_r.     

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.     

Comparative evaluation of the performance of sinuous symmetrically-heated natural convection channels

Buoyancy-driven flows established in vertical channels are investigated numerically. Two different sinuous morphologies for a vertical channel are proposed as alternatives to a typical (straight) vertical channel; one is a sinus-shaped channel, and the other is a three-segment channel. The objective is to carry out a comparative evaluation of the thermal and dynamic performance of the proposed configurations, comparing their performance with that corresponding to a typical vertical channel. The low-Reynolds k − ω turbulence model is employed to simulate the transitional or fully turbulent flow. The average Nusselt number and the dimensionless mass flow rate are obtained for Rayleigh numbers ranging 10² ≤ Ra_H ≤ 10¹⁵, with an aspect ratio of the channel equal to 0.1. Sizeable increases in the Nusselt number and the mass flow rate are found for the sinuous configurations, under given conditions.     

Natural convection effects on heat and mass transfer in a curvilinear triangular cavity

Finite element method is used in this study to analyze the effects of buoyancy ratio and Lewis number on heat and mass transfer in a triangular cavity with zig-zag shaped bottom wall. Buoyancy ratio is defined as the ratio of Grashof number of solutal and thermal. Inclined walls of the cavity have lower temperature and concentration according to zig-zag shaped bottom wall. Enclosed space consists mostly of an absorber plate and two inclined glass covers that form a cavity. Both high temperature and high concentrations are applied to bottom corrugated wall. Computations were done for different values of buoyancy ratio (?10 ? Br ? 10), Lewis number (0.1 ? Le ? 20) and thermal Rayleigh number (10⁴ ? Ra_T ? 10⁶). Streamlines, isotherms, iso-concentration, average Nusselt and Sherwood numbers are obtained. It is found that average Nusselt and Sherwood numbers increase by 89.18% and 101.91% respectively as Br increases from ?10 to 20 at Ra_T = 10⁶. Also, average Nusselt decreases by 16.22% and Sherwood numbers increases by 144.84% as Le increases from 0.1 to 20 at this Rayleigh number.     

Magnetic field and internal heat generation effects on the free convection in a rectangular cavity filled with a porous medium

A numerical investigation of the steady magnetohydrodynamics free convection in a rectangular cavity filled with a fluid-saturated porous medium and with internal heat generation has been performed. A uniform magnetic field, inclined at an angle γ with respect to the horizontal plane, is externally imposed. The values of the governing parameters are the inclined angle γ = 0, π/6, π/4 and π/2, Hartmann number Ha = 0, 1, 5, 10 and 50, Rayleigh number Ra = 10, 100, 10³ and 10⁵, and the aspect ratio a = 0.01, 0.2, 0.5 and 1 (square cavity). It is shown that the intensity of the core convection is considerably affected by the considered parameters. It is also found that the local Nusselt number Nu_Y decreases on the bottom wall as γ increases (magnetic field changes its direction from the horizontal to the vertical direction) and vice versa for the top wall of the cavity. The reported results are in good agreement with the available published work in the literature.     

Simulation of high Rayleigh number natural convection in a square cavity using the lattice Boltzmann method

A thermal lattice Boltzmann method based on the BGK model has been used to simulate high Rayleigh number natural convection in a square cavity. The model uses the double populations approach to simulate hydrodynamic and thermal fields. The traditional lattice Boltzmann method on a uniform grid has unreasonably high grid requirements at higher Rayleigh numbers which renders the method impractical. In this work, the interpolation supplemented lattice Boltzmann method has been utilized. This is shown to be effective even at high Rayleigh numbers. Numerical results are presented for natural convection in a square cavity with insulated horizontal walls and isothermal vertical walls maintained at different temperatures. Very fine grids (wall y⁺ < 0.3) have been used for the higher Rayleigh number simulations. A universal structure is shown to exist in the mean velocity turbulent boundary layer profile for y⁺ < 10. This agrees extremely well with previously reported experimental data. The numerical results (for Rayleigh numbers up to 10¹⁰) are in very good agreement with the benchmark results available in the literature. The highlight of the calculations is that no turbulence model has been employed.     

Radiative effects on natural convection in vertical convergent channels

Natural convection in air, in a convergent channel, uniformly heated at the principal walls, is experimentally investigated, in order to analyze the effects of the radiative heat transfer. Results in terms of wall temperature profiles as a function of the walls inclination angle, the spacing between the walls, the heat flux, are given for two values of the wall emissivity. Flow visualization is carried out to show the peculiar pattern of the flow between the plates in several configurations. The comparison between two wall emissivity values, 0.10 and 0.90, shows that the effect of thermal radiation is more pronounced for larger convergence angles. For a wall emissivity equal to 0.90 and for small values of the minimum channel spacing, heat transfer in slightly convergent vertical channels is stronger than in a vertical parallel channel. Flow visualization points out a recirculating zone in the upper part of the channel for small values of the minimum channel spacing and for converging angles equal to 5° and 10°. Nusselt numbers and dimensionless maximum temperatures are then evaluated and correlated to the Rayleigh number, in the investigated range from 5 to 5 × 10⁸ and 0° ? θ ? 10°. A very good agreement between experimental data and correlations is observed for the dimensionless parameters based on the maximum channel spacing. Comparisons between experimental and numerical data are also performed and a good relationship is observed.     

Final steady flow near a stagnation point on a vertical surface in a porous medium

This paper investigates the large time (final state flow) solutions for unsteady mixed convection boundary layer flow near a stagnation point on a vertical surface embedded in a Darcian fluid-saturated porous medium. Through numerical computations Nazar et al. [R. Nazar, N. Amin, I. Pop, Unsteady mixed convection boundary layer flow near the stagnation point on a vertical surface in a porous medium, Int. J. Heat Mass Transfer 47 (2004) 2681–2688] concluded that for values of the mixed convection parameter λ > −1, the governing boundary value problem (BVP) had a unique solution. If λ_c ≈ −1.4175 < λ ⩽ −1 two solutions were reported, and if λ < λ_c then no solutions were found. The purpose of this note is to provide further mathematical and numerical analysis of this problem. We prove existence of a solution to the governing BVP for all λ > −1. We also present numerical evidence that a second solution exists for λ > −1, thus giving dual solutions for all λ > λ_c. It is also proven that if λ < −2.9136 no solution to the BVP exists. Finally, a stability analysis is performed to show that solutions on the upper branch are linearly stable while those on the lower branch are linearly unstable.     

Conjugate heat transfer in open cavities with a discrete heater at its optimized position

In this article, we determined optimum position of a discrete heater by maximizing the conductance and then studied heat transfer and volume flow rate with the discrete heater at its optimum position in open cavities. Continuity, Navier–Stokes and energy equations are solved by finite difference-control volume numerical method. The relevant governing parameters were: the Rayleigh numbers from 10⁶ to 10¹², the Prandtl number, Pr = 0.7, the cavity aspect ratio, A = H/L from 0.5 to 2, the wall thickness l/L from 0.05 to 0.15, the heater size h/L from 0.15 to 0.6, and the conductivity ratio k_r from 1 to 50. We found that the global conductance is an increasing function of the Rayleigh number, the conductivity ratio, and a decreasing function of the wall thickness. Best thermal performance is obtained by positioning the discrete heater at off center and slightly closer to the bottom. The Nusselt number and the volume flow rate in and out the open cavity are an increasing function of the Rayleigh number and the wall thickness, and a decreasing function of the conductivity ratio. The Nusselt number is a decreasing function of the cavity aspect ratio and the volume flow rate is an increasing function of it.