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.     

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.     

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.     

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.     

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 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.     

Unsteady conjugate natural convection in a square enclosure filled with a porous medium

Mathematical simulation of unsteady natural convection modes in a square cavity filled with a porous medium having finite thickness heat-conducting walls with local heat source in conditions of heterogeneous heat exchange with an environment at one of the external boundaries has been carried out. Numerical analysis was based on Darcy–Forchheimer model in dimensionless variables such as a stream function, a vorticity vector and a temperature. The special attention was given to analysis of Rayleigh number effect Ra = 10⁴, 10⁵, 10⁶, of Darcy number effect Da = 10^?5, 10^?4, 10^?3, ∞, of the transient factor effect 0 < τ < 1000 and of the heat conductivity ratio k_2,1 = 3.7 × 10^?2, 5.7 × 10^?4, 6.8 × 10^?5 on the velocity and temperature fields. The influence scales of the defining parameters on the average Nusselt number have been detected.     

Jet impingement cooling of a constant heat flux horizontal surface in a confined porous medium: Mixed convection regime

In the present study, numerical investigation of jet impingement cooling of a constant heat flux horizontal surface immersed in a confined porous channel is performed under mixed convection conditions with the limitation of the Darcy model. The results are presented in the mixed convection regime with wide ranges of the governing parameters: Péclet number (1 ? Pe ? 1000), Rayleigh number (10 ? Ra ? 100), half jet width (0.1 ? D ? 1.0), and the distance between the jet and the heated portion (0.1 ? H ? 1.0). It is found that the average Nusselt number increases with increase in either Rayleigh number or jet width for high values of Péclet number. The average Nusselt number also increases with decrease in the distance between the jet and the heated portion. The correlation for Nu_avg in the forced convection regime is suggested. It is shown that mixed convection mode can cause minimum average Nusselt number unfavorably due to counteraction of jet flow against buoyancy driven flow. Hence, careful consideration must be given while designing a system of jet impingement cooling through porous medium.     

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.     

Numerical simulation of natural convection heat transfer in a porous cavity heated from below using a non-Darcian and thermal non-equilibrium model

The present paper investigates the numerical simulation of steady laminar incompressible natural convection heat transfer in an enclosed cavity that is filled with a fluid-saturated porous medium. The bottom wall is subjected to a relatively higher temperature than the top wall while the vertical walls are considered to be insulated. The flow field is modeled upon incorporating different non-Darcian effects, such as the convective term, Brinkman effect and Forchhiemer quadratic inertial effect. Moreover the two-equation model is used to separately account for the local fluid and solid temperatures. The numerical solution is obtained through the application of the finite volume method. The appraisals of the sought objectives are performed upon identifying key dimensionless groups of parameters. These dimensionless groups along with their operating domains are: Rayleigh number 1 ⩽ Ra ⩽ 400, Darcy number 10⁻⁴ ⩽ Da ⩽ 10⁻³, effective fluid-to-solid thermal conductivity ratio 0.1 ⩽ κ ⩽ 1.0, and the modified Biot number 1 ⩽ χ ⩽ 100. The non-Darcian effects are first examined over a broad range of Rayleigh number. Next, the implications of the group of parameters on the flow circulation intensity, local thermal non-equilibrium (LTNE) and average Nusselt number are highlighted and pertinent observations are documented.     

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.     

Laminar free convection inside an inclined L-shaped enclosure

We have numerically reported the buoyancy induced flow and heat transfer characteristics inside an inclined L-shaped enclosure. A control volume based Finite-Volume method is applied to discretize the governing equations with collocated variable arrangement. SIMPLE algorithm is used and the system of equations is solved by Stone's SIP solver with full multigrid acceleration. Results are presented in the form of the average Nusselt number for a range of inclination angle, θ = 0°–360°; Rayleigh number, Ra = 1–10⁵; and aspect ratio, A = 0.1–0.5.     

Jet impingement cooling of a horizontal surface in an unconfined porous medium: Mixed convection regime

Two-dimensional slot jet impingement cooling of an isothermal horizontal surface immersed in an unconfined porous medium is simulated numerically to gain insight into thermal characteristics under mixed convection conditions with the limitation of the Darcy model. The jet direction is considered to be perpendicular from the top to the horizontal heated element; therefore, the jet flow and the buoyancy driven flow are in opposite directions. The results are presented in the mixed convection regime with wide ranges of the governing parameters: Péclet number (1 ? Pe ? 1000), Rayleigh number (10 ? Ra ? 100), half jet width (0.1 ? D ? 0.5), and the distance between the jet and the heated portion (0.1 ? H ? 1.0). It is found that the average Nusselt number increases with increase in either Rayleigh number or jet width for high values of Péclet number. The average Nusselt number also increases with decrease in the distance between the jet and the heated portion. It is shown that mixed convection mode can cause minimum average Nusselt number at two values of Péclet number and a maximum average Nusselt number occurs in between theses two Péclet numbers at higher Rayleigh number due to counteraction of jet flow against buoyancy driven flow. Hence careful consideration must be given while designing a system of jet impingement cooling through porous medium.     

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 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.     

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.     

Critical Rayleigh number of natural convection in high porosity anisotropic horizontal porous layers

The critical Rayleigh number Ra_C at the onset of natural convection was studied by linear stability analysis for high porosity anisotropic horizontal porous layers of uniformly arraying vertical thin circular wires stretched across a hot and a cold surface. Navier–Stokes equations including flow resistance by the wires were solved since Darcy’s law cannot be applied due to high porosity. Ratio of permeability in the horizontal direction to the vertical direction is constant to be 0.5 but ratio of effective thermal diffusivity ξ is changed dependent on wire materials. The critical Rayleigh number Ra_C for the case of ξ = 1 obtained by the analysis agreed well with the experiment.     

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.     

Laminar natural convection of non-Newtonian power-law fluids between concentric circular cylinders

The two-dimensional steady-state natural convection of power-law fluids is studied numerically between two concentric horizontal cylinders with different constant temperatures. The governing equations are discretized using finite volume technique based on second order upwind and are solved using the SIMPLE algorithm. The effects of Rayleigh number (10³ ≤ Ra ≤ 10⁵) and Prandtl number (10 ≤ Pr ≤ 10³) on the dimensionless velocity and temperature are investigated for both pseudoplastic and dilatant fluids. Also the mean Nusselt number for various values of governing parameters is obtained and discussed. The results indicate that with increasing the power-law index from 0.6 to 1.4, the mean Nusselt number decreases. In the best case among the range of parameters considered here the heat transfer rate for pseudo-plastic fluid (n = 0.6) is 170% higher than the Newtonian one and for dilatant fluid (n = 1.4) the heat transfer rate is 43% lower than the Newtonian fluid. So the pseudoplastic and dilatant fluids are more efficient than Newtonian fluids for cooling and insulating purposes, respectively. It is shown that as the Rayleigh number increases the cooling effect of pseudoplastic fluid and the insulating effect of dilatant fluid become more pronounced.