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.     

Effect of nanofluid variable properties on natural convection in enclosures filled with a CuO–EG–Water nanofluid

This work focuses on the study of natural convection heat transfer characteristics in a differentially-heated enclosure filled with a CuO–EG–Water nanofluid for different published variable thermal conductivity and variable viscosity models. The problem is given in terms of the vorticity–stream function formulation and the resulting governing equations are solved numerically using an efficient finite-volume method. Comparisons with previously published work are performed and the results are found to be in good agreement. Various results for the streamline and isotherm contours as well as the local and average Nusselt numbers are presented for a wide range of Rayleigh numbers (Ra = 10³–10⁵), volume fractions of nanoparticles (0 ≤ φ ≤ 6%), and enclosure aspect ratios (½ ≤ A ≤ 2). Different behaviors (enhancement or deterioration) are predicted in the average Nusselt number as the volume fraction of nanoparticles increases depending on the combination of CuO–EG–Water variable thermal conductivity and viscosity models employed. In general, the effects the viscosity models are predicted to be more predominant on the behavior of the average Nusselt number than the influence of the thermal conductivity models. The enclosure aspect ratio is predicted to have significant effects on the behavior of the average Nusselt number which decreases as the enclosure aspect ratio increases.     

Experimental investigation of mixed convection heat transfer in a rectangular channel with discrete heat sources at the top and at the bottom

Mixed convection heat transfer in a top and bottom heated rectangular channel with discrete heat sources has been investigated experimentally for air. The lower and upper surfaces of the channel were equipped with 8 × 4 flush-mounted heat sources subjected to uniform heat flux. Sidewalls, the lower and upper walls were insulated and adiabatic. The experimental study was made for an aspect ratio of AR = 6, Reynolds numbers 955 ≤ Re_Dh ≤ 2220 and modified Grashof numbers Gr* = 1.7 × 10⁷ to 6.7 × 10⁷. From experimental measurements, surface temperature and Nusselt number distributions of the discrete heat sources were obtained for different Grashof numbers. Furthermore, Nusselt number distributions were calculated for different Reynolds numbers. Results show that surface temperatures increase with increasing Grashof number. The row-averaged Nusselt numbers first decrease with the row number and then, due to the increase in the buoyancy affected secondary flow and the onset of instability, they show an increase towards the exit as a result of heat transfer enhancement.     

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.     

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.     

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.     

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.     

Lattice Boltzmann simulation of natural convection heat transfer in an elliptical-triangular annulus

A numerical study for steady-state, laminar natural convection in a horizontal annulus between a heated triangular inner cylinder and cold elliptical outer cylinder was investigated using lattice Boltzmann method. Both inner and outer surfaces are maintained at the constant temperature and air is the working fluid. Study is carried out for Rayleigh numbers ranging from 1.0 × 10³ to 5.0 × 10⁵. The effects of different aspect ratios and elliptical cylinder orientation were studied at different Rayleigh numbers. The local and average Nusselt numbers and percentage of increment heat transfer rate were presented. The average Nusselt number was correlated. The results show that by decreasing the value of aspect ratio and/or increasing the Rayleigh number, the Nusselt number increases. Also the heat transfer rate increases when the ellipse positioned vertically.     

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.     

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 flow of Cu–Water nanofluid in horizontal cylindrical annuli with inner triangular cylinder using lattice Boltzmann method

In this paper the lattice Boltzmann method is used to investigate the effect of nanoparticles on natural convection heat transfer in two-dimensional horizontal annulus. The study consists of an annular-shape enclosure, which is created between a heated triangular inner cylinder and a circular outer cylinder. The inner and outer surface temperatures were set as hot (T_h) and cold temperatures (T_c), respectively and assumed to be isotherms. The effect of nanoparticle volume fraction to the enhancement of heat transfer was examined at different Rayleigh numbers. Furthermore, the effect of vertical, horizontal, and diagonal eccentricities at various locations is examined at Ra = 10⁴. The result is presented in the form of streamlines, isotherms, and local and average Nusselt number. Results show that the Nusselt number and the maximum stream functions increase by augmentation of solid volume fraction. Average Nusselt number increases when the inner cylinder moves downward, but it decreases, when the location of inner cylinder changes horizontally.     

Effect of a magnetic field on natural convection in an open cavity subjugated to water/alumina nanofluid using Lattice Boltzmann method

In this paper, the effect of a magnetic field on natural convection in an open enclosure which subjugated to water/alumina nanofluid using Lattice Boltzmann method has been investigated. The cavity is filled with water and nanoparticles of Al₂O₃ at the presence of a magnetic field. Calculations were performed for Rayleigh numbers (Ra = 10⁴–10⁶), volume fractions of nanoparticles (φ = 0,0.02,0.04 and .0.06) and Hartmann number (0 ≤ Ha ≤ 90) with interval 30 while the magnetic field is considered horizontally. Results show that the heat transfer decreases by the increment of Hartmann number for various Rayleigh numbers and volume fractions. The magnetic field augments the effect of nanoparticles at Rayleigh number of Ra = 10⁶ regularly. Just as the most effect of nanoparticles for Ra = 10⁴ is observed at Ha = 30, so the most influence of nanoparticles occurs at Ha = 60 for Ra = 10⁵.     

Effect of nanofluid variable properties on natural convection in enclosures

In this work, the heat transfer enhancement in a differentially heated enclosure using variable thermal conductivity and variable viscosity of Al₂O₃–water and CuO–water nanofluids is investigated. The results are presented over a wide range of Rayleigh numbers (Ra = 10³–10⁵), volume fractions of nanoparticles (0 ≤ φ ≤ 9%), and aspect ratios (½ ≤ A ≤ 2). For an enclosure with unity aspect ratio, the average Nusselt number of a Al₂O₃–water nanofluid at high Rayleigh numbers was reduced by increasing the volume fraction of nanoparticles above 5%. However, at low Rayleigh numbers, the average Nusselt number was slightly enhanced by increasing the volume fraction of nanoparticles. At high Rayleigh numbers, CuO–water nanofluids manifest a continuous decrease in Nusselt number as the volume fraction of nanoparticles is increased. However, the Nusselt number was not sensitive to the volume fraction at low Rayleigh numbers. The Nusselt number demonstrates to be sensitive to the aspect ratio. It was observed that enclosures, having high aspect ratios, experience more deterioration in the average Nusselt number when compared to enclosures having low aspect ratios. The variable thermal conductivity and variable viscosity models were compared to both the Maxwell-Garnett model and the Brinkman model. It was found that at high Rayleigh numbers the average Nusselt number was more sensitive to the viscosity models than to the thermal conductivity models.     

On thermal control of devices contained in inclined hemispherical cavities with dome oriented downwards and subjected to transient natural convection

In this study, relationships of Nusselt–Rayleigh–Fourier type are proposed for the case of air-filled hemispherical cavity whose dome is oriented downwards and maintained isothermal. Its disk is subjected to a constant heat flux and inclined at an angle varying between 90° (vertical position) and 180° (disk horizontal with dome oriented downwards). The numerical approach is performed in transient regime by means of the finite volume method for Rayleigh numbers in the range of 10⁴ − 5 × 10⁸. These results are confirmed at steady state by measurements done for some configurations in a previous study for the same Rayleigh and inclination ranges. Otherwise, they complete other surveys considering inclination angles varying between 0° (horizontal cavity with dome oriented upwards) and 90° (vertical cavity) for a wider range of Rayleigh numbers. The correlations allow thermal control of devices submitted to natural convection in hemispherical cavities during the time preceding the steady state after their switch on.     

Role of finite element based grids and simulations on evaluation of Nusselt numbers for heatfunctions within square and triangular cavities involving multiple discrete heaters

Nusselt number is an important non-dimensional parameter which quantifies the heat transfer rate. Local Nusselt number is useful in predicting the heat transfer rate along the various hot and cold sections of the side walls in a discretely heated enclosed cavity. In addition, the overall heat balance in an enclosed cavity (total heat delivered by the hot isothermal walls should be equal to the total heat gained by the cold isothermal walls) can be validated via the average Nusselt numbers. Current finite element based simulations and post-processing have been carried out in order to analyze the influence of the multiple heaters on the Nusselt number along various sections (hot and cold) of the side walls in discretely heated square and triangular (design 1 and design 2) cavities. The working fluid is considered to be air (Pr = 0.7) and the numerical studies have been carried out for a large range of Rayleigh number (Ra = 10³–10⁵) for four different biquadratic elements (24 × 24, 28 × 28, 32 × 32 and 34 × 34). The current work also estimates the fractional error in the heat balance (ϵ) and it is clearly inferred that ϵ is comparatively lower for 34 × 34 biquadratic elements. Current work also reveals that the fractional error (ϵ) is mainly induced due to the sharp variations in the Nusselt number at the cold-hot junctions along the side walls. The present study also involves the detailed evaluation of the heatfunction (Π) expressions along the cold-hot junctions of the side walls. The computations of the heatfunctions are intrinsically related to the Nusselt numbers of the hot-cold junctions.     

Numerical investigation for natural convection of a nanofluid in an inclined L-shaped cavity in the presence of an inclined magnetic field

The problem of natural convection in an inclined L-shaped enclosure filled with Cu/water nanofluid that operates under differentially heated walls in the presence of an inclined magnetic field is presented in this paper. The fully implicit finite difference method is used to solve the governing equations. A comparison with previously published results in special case of the present study is performed and a very good agreement is found. Heat transfer and fluid flow are examined for parameters of the Hartmann number (0 ≤ Ha ≤ 100), the nanoparticles volume fraction (0% ≤ ϕ ≤ 20%), the cavity inclination angle (0° ≤ ϑ ≤ 300°), the magnetic field inclination angle (0° ≤ γ ≤ 270°), the cavity aspect ratio (0.25 ≤ AR ≤ 0.6) and the Rayleigh number (10³ ≤ Ra ≤ 10⁶). It is found that, the presence of the magnetic field in the fluid region causes a significant reduction in the fluid flow and heat transfer characteristics. Also, a good enhancement in the heat transfer rate can be obtained by adding the copper nanoparticles to the base fluid.     

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.     

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.     

Numerical investigation on heat transfer and friction factor characteristics of laminar and turbulent flow in an elliptic annulus utilizing nanofluid

In this paper, a numerical investigation on heat transfer performance and flow fields of different nanofluids flows through elliptic annulus in a laminar and turbulent flow regimes. The three-dimensional continuity, Navier–Stokes and energy equations are solved by using finite volume method (FVM) and the SIMPLE algorithm scheme is applied to examine the effects of laminar and turbulent flow on heat transfer characteristics. This study evaluates the effects of four different types of nanoparticles, Al₂O₃, CuO, SiO₂ and ZnO, with different volume fractions (0.5–4%) and diameters (25–80 nm) under constant heat flux boundary condition using water as a base fluid were used. The Reynolds number of laminar flow was in the range of 200 ≤ Re ≤ 1500, while for turbulent flow it was in the range of 4000 ≤ Re ≤ 10,000. The results have shown that SiO₂–water nanofluid has the highest Nusselt number, followed by ZnO–water, CuO–water, Al₂O₃–water, and lastly pure water. The Nusselt number for all cases increases with the volume fraction but it decreases with the rise in the diameter of nanoparticles. In all configurations, the Nusselt number increases with Reynolds number. It is found that the glycerine–SiO₂ shows the best heat transfer enhancement compared with other tested base fluids.     

Effect of solid volume fraction and tilt angle in a quarter circular solar thermal collectors filled with CNT–water nanofluid

Solar thermal collectors have significant importance due to its wide use in solar thermal technology. Augmentation of heat transfer is a key challenge for solar thermal technology. A quarter circular solar thermal collectors is investigated throughout the paper introducing carbon nanotube (CNT)–water nanofluid in the cavity. Tilt angle of this type of collector plays a vital role and heat transfer can be maximized for a particular tilt angle and solid volume fraction of the nanofluid. Galerkin weighted residual of FEM has been applied for the numerical solution of the problem. Grid independency test and code validation have been assessed for the accuracy of numerical solution. In this paper a wide range of solid volume fraction (δ = 0 to δ = 0.12) and tilt angle (ϕ = 0 to ϕ = 60°) has been investigated for Rayleigh number (Ra = 10⁵–10⁸) with varying dimensionless times. It has been found that both solid volume fraction and tilt angle play vital roles for the augmentation of heat transfer and a good heat transfer characteristic can be obtained by compromising between these two parameters. The results are shown using streamline, isotherm contour and related graph and chart.