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.     

Effect of magnetic field on convection heat transfer inside a tilted square enclosure

Steady, laminar, natural-convection flow in the presence of a magnetic field in a tilted enclosure heated from below and cooled from top is considered. The enclosure is filled with liquid gallium. In our formulation of governing equations, mass, momentum and energy are applied to the enclosure. To solve the nonlinear governing differential equations a finite volume code based on PATANKAR's SIMPLER method is utilized. It is shown that for a given inclination angle (φ), as the value of Hartmann number (Ha) increases, the convection heat transfer reduces. Furthermore it is found that at Ra = 10⁴, value of Nusselt number depends strongly upon the inclination angle for relatively small values of Hartmann number. At Ra = 10⁵, the Nusselt number increases up to about φ = 45^o and then decrease as φ increases.     

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.     

Conjugate natural convection in the roof cavity of heavy construction building during summer

This paper presents the results of a study of conjugate natural convection inside a building attic in the shape of a rectangular enclosure bounded by realistic walls made from composite construction materials under summer day boundary conditions. The effects of cavity aspect ratio, Rayleigh number (Ra), and orientation of the external surfaces on the flow and heat transfer characteristics were the main focus of the investigation. The problem was formulated in terms of the vorticity-stream function procedure, and the governing equations for steady, laminar, two-dimensional conjugate natural convection heat transfer were solved by employing the Alternating Direction Implicit (ADI) control volume method along with under-relaxation factors for temperature, vorticity, and stream functions. For Ra ranging from 10³ to 10¹⁰, steady state results of the streamline and temperature contours in addition to local and mean Nusselt number at all surfaces of the cavity were obtained. The results show that the values of Ra and the aspect ratio have significant effect on the temperature and stream function contours within the enclosure. Another important finding of the study is that heat flux into the room increases with the increase of both the aspect ratio and Rayleigh number.     

Heatlines based natural convection analysis in tilted isosceles triangular enclosures with linearly heated inclined walls: effect of various orientations

Natural convection in isosceles triangular enclosures with various configurations (case 1 — inverted, case 2 — straight and case 3 — tilted) is studied via heatline analysis for linear heating of inclined walls. Detailed analysis and comparison for various base angles (φ = 45°, 60°) of triangular enclosures have been carried out for a range of fluids (Pr = 0.015 − 1000) within Ra = 10³ − 10⁵ using Galerkin finite element method. The heat flow distributions indicate conduction dominant heat transfer at low Ra (Ra = 10³) for case 1 and case 2 whereas in case 3, convective heat flow is observed due to high buoyancy force. As Ra increases, enhanced thermal mixing is observed at the core of the cavity. Wall to wall heat transfer occurs at walls AB and AC due to linear heating boundary condition in all the cases. Although the distributions of fluid flow and heat flow are qualitatively similar for φ = 45° and 60°, the intensity of fluid flow and heat flow decreases as φ increases. Strength of fluid flow and heat flow circulation cells is found to be higher in case 3 for identical parameters. Results show that upper side wall (AC) for case 3 exhibits higher heat transfer rates whereas heat transfer rates for walls AB and AC are the same for case 1 and case 2. Also Nu_AB is higher for case 2 followed by case 1 and case 3 at the middle portion of wall AB. Thus to achieve high heat transfer from fluid to wall at the central region, case 2 and case 3 configurations may be recommended at high Ra (Ra = 10⁵) and Pr, irrespective of φ.     

Characteristics of heat transfer and fluid flow in a channel with single-row plates array oblique to flow direction for photovoltaic/thermal system

This study has been carried out to investigate the characteristics of convective heat transfer and fluid flow for a single row of oblique plates array to the flow direction inside a channel. The flow inside the channel is laminar and the plates array have spanwise distance between the plates and heated by radiation. This configuration has been designed to be used for Photovoltaic/Thermal system (PV/T) applications. The theoretical results are validated with measured values, and a good agreement prevailed. The results show that an increase in the plate oblique angle (γ) in the range from 0 to 15 degrees, leads to an increase in the Nusselt number (Nu) up to a maximum value and then decreases. The oblique angle at the maximum value of Nu depends on the flow Reynolds Number (Re), and (?_w/?_pl), where (?_w/?_pl) is defined as the ratio of the plates’ spacing at zero oblique angle to the plate length. Furthermore, increasing (?_w/?_pl) results in a significant increase in the heat transfer coefficient depending on the values of Re, and plate oblique angle (γ). In addition, increasing (γ) from 0 to 15 degrees results in a decrease in the friction factor up to a certain value, after which the friction value approaches a constant value depending on Re value and (?_w/?_pl). It was found that for any value of the plate oblique angle (γ), the friction factor decreases with the increase of the values of (?_w/?_pl) and Re, respectively.     

Bejan’s heatline analysis of natural convection in right-angled triangular enclosures: Effects of aspect-ratio and thermal boundary conditions

Natural convection in right-angled triangular enclosures with various top angles ($\phi$=15°, 30°, 45°) is studied in detail via heat flow analysis for various uniform isothermal and linear isothermal heating thermal boundary conditions. Detailed analysis on the effects of aspect-ratio and thermal boundary conditions on the fluid and heat flow inside the triangular enclosures have been carried out for a range of fluids (Pr = 7.2, 1000, 0.015) within Ra = 10³–10⁵. Interesting features of heat flow patterns under various thermal boundary conditions are ‘visualized’ by heatlines. The effect of increase in φ of triangular enclosures is such that the maximum heat flux at the top vertex decreases and the thermal mixing in cavity increases with the increase in $\phi$. It is found that, the fluid in the lower corners is adequately heated in presence of hot right wall compared to that in left wall heating cases. Further, the heat transfer characteristics, in terms of local and average Nusselt numbers, indicate that isothermal heating cases exhibit exponential decrease in Nu_l whereas linear heating cases interestingly show local intermediate maxima. Also, various qualitative and quantitative features of Nu and $\overline{Nu}$ are adequately explained based on heatlines. Finally, the correlations for $\overline{N{u}_l}$ and Ra are obtained for various fluid with all heating situations.     

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.     

Numerical study on effect of vent locations on natural convection in an enclosure with an internal heat source

Natural convection is a widely studied phenomenon because of the extensive applications in cooling of large scale electrical and electronic equipments. The current study involves study of effect of vent locations on natural convection in enclosures with partial openings having an internal heat source. It involves the numerical simulation of 2D steady state natural convection in enclosure of different aspect ratios (H/W = 1, 2 and 3) for lower Rayleigh numbers (Ra_h = 10³, 10⁴ and 10⁵). Four different configurations have been considered based on the number and position of vents — same side (SS), diagonal side (DS), one inlet two outlets (1I2O) and two inlets one inlet (2I1O). The mass flow rate driven through the enclosure and the average Nusselt number over the heater surface for all the four configurations have been compared. It is found that the 2I1O configuration yielded better heat transfer rates of the four considered. It was found that the mass flow rates and Nu increased with increase in Ra_h and decrease in the aspect ratio.     

Natural convection slip flow in a vertical microchannel heated at uniform heat flux

Numerical solutions for steady state developing natural convection flow in air, in vertical parallel-plate microchannels are accomplished. An asymmetric heating is considered and the walls are assumed to be at uniform heat flux. A first-order model is used for slip and jump boundary conditions and an analytical solution for the fully developed flow is also given. Results are performed for air, for the heat flux ratio in the 0.0–1.0 range, for Rayleigh, Ra, and Knudsen, Kn, numbers from 10^?1 to 8 × 10³ and from 0.0 to 0.10, respectively. The maximum mass flow rate is always obtained for the highest considered Kn value, whereas the average Nusselt number, Nu, increases for lower Ra (<10) and decreases for Ra > 100. Wall temperature profiles have the lowest values for highest considered Kn value at lower Ra, whereas for the developing flow, they present opposite trends. For developing flow, velocity profiles for asymmetric and symmetric heating are completely different. In developing flow velocity profiles along the wall present the highest increases for asymmetric heating and the highest considered Kn value.     

Mixed convection between horizontal plates—II. Fully developed flow

Fully developed velocity profiles of longitudinal convection rolls in mixed convection between horizontal plates were measured in nitrogen by laser Doppler anemometry for a range 2472 < Ra < 8300 and 15 < Re < 150. It is shown analytically and experimentally that the transverse velocities of the longitudinal convection rolls are independent of the forced flow. The experimentally and numerically obtained w-profiles (Pr = 0.71) are in good agreement with theoretical predictions (Pr → ∞) and other experimental results (Pr = 11.1 and 930) for Rayleigh-Benard convection. A detailed study of the longitudinal velocity modulation Δu[w_max(Ra), Re] is presented. Also, asymmetric roll patterns were found in spite of the small temperature differences used between the horizontal plates.     

Heat transfer correlations for PCM-based heat sinks with plate fins

Characterization of melting process in a Phase Change Material (PCM)-based heat sink with plate fin type thermal conductivity enhancers (TCEs) is numerically studied in this paper. Detailed parametric investigations are performed to find the effect of aspect ratio of enclosure and the applied heat flux on the thermal performance of the heat sinks. Various non-dimensional numbers, such as Nusselt number (Nu), Rayleigh number (Ra), Stefan number (Ste) and Fourier number (Fo) based on a characteristic length scale, are identified as important parameters. The half fin thickness and the fin height are varied to obtain a wide range of aspect ratios of an enclosure. It is found that a single correlation of Nu with Ra is not applicable for all aspect ratios of enclosure with melt convection taken into account. To find appropriate length scales, enclosures with different aspect ratios are divided into three categories, viz. (a) shallow enclosure, (b) rectangular enclosure and (c) tall enclosure. Accordingly, an appropriate characteristic length scale is identified for each type of enclosure and correlation of Nu with Ra based on that characteristic length scale is developed.     

Visualization of heat flow due to natural convection within triangular cavities using Bejan’s heatline concept

Natural convection and flow circulation within a cavity has received significant attention in recent times. The wide range of applicability of flow inside a cavity (food processing industries, molten metal industries, etc.) requires thorough understanding for cost efficient processes. This paper is based on comprehensive analysis of heat flow pattern using Bejan’s heatline concept. The key parameters for our study are the Prandtl number, Rayleigh number and Nusselt number. The values of Prandtl number (0.015, 0.026, 0.7 and 1000) have been chosen based on wide range of applicability. The Rayleigh number has been varied from 10² to 10⁵. Interesting results were obtained. For low Rayleigh number, it is found that the heatlines are smooth and perfectly normal to the isotherms indicating the dominance of conduction. But as Ra increases, flow slowly becomes convection dominant. It is also observed that multiple secondary circulations are formed for fluids with low Pr whereas these features are absent in higher Pr fluids. Multiple circulation cells for smaller Pr also correspond multiple cells of heatlines which illustrate less thermal transport from hot wall. On the other hand, the dense heatlines at bottom wall display enhanced heat transport for larger Pr. Further, local heat transfer (Nu_l, Nu_t) are explained based on heatlines. The comprehensive analysis is concluded with the average Nusselt number plots. A correlation for average heat transfer rate and Ra has been developed and the range of Rayleigh number is also found, to depict the conduction dominant heat transfer.     

Natural convection heat transfer in horizontal and vertical closed narrow enclosures with heated rectangular finned base plate

Natural convection heat transfer and fluid flow characteristics in horizontal and vertical narrow enclosures with heated rectangular finned base plate have been experimentally investigated at a wide range of Rayleigh number (Ra) for different fin spacings and fin lengths. Quantitative comparisons of finned surface effectiveness (ε) and heat transfer rate between horizontal and vertical enclosures have been reported. In comparison with enclosure of a bare base plate, insertion of heat conducting fins always enhances heat transfer rate. Optimization of fin-array geometry has been addressed. The results gave an optimum fin spacing at which Nusselt number (Nu_H) and finned surface effectiveness (ε) are maximum. It has been found that: (1) increasing fin length increases Nu_H and ε; (2) increasing Ra increasesNu_H for any fin-array geometries and (3) for any fin-array geometry and at Ra > 10,000, increasing Ra decreases ε while for fin-array geometries of large fin spacing and at Ra < 10,000, increasing Ra increases ε. Useful design guidelines have been suggested. Correlations of Nu_H have been developed for horizontal and vertical enclosures. Correlations predictions have been compared with previous data and good agreement was found.     

Numerical simulation of liquid metals in differentially heated enclosure undergoing orthogonal rotation

In this work a numerical investigation has been performed to examine the characteristics of mixed convective heat transfer in square enclosures undergoing orthogonal rotation i.e. rotation axis and gravity axis are orthogonal to each other. A semi implicit finite difference code on a collocated grid is used to solve the momentum and energy equations subject to Boussinesq approximation. The study is carried out for a wide range of operating parameters such as Rayleigh number (Ra), Taylor number (Ta), Rotational Rayleigh number (Ra_w) for a fixed Prandtl number (Pr). The numerical experiments have been carried out for a fixed Pr = 0.01, Ra varies from 10⁵ to 10⁷ while Ta and Ra_w vary from almost 0 to 10⁹. Results reveal that significant increase or decrease in heat transfer rates can be achieved by the rotational effects, mainly influenced by centrifugal force.     

Combined radiation and convection heat transfer in a porous channel bounded by isothermal parallel plates

Combined radiation and convection heat transfer in a porous medium confined between gray isothermal parallel plates is investigated. The medium is absorbing, emitting and scattering. Cases of boundaries at temperatures higher or lower than the medium are considered. In the porous medium, the boundary effect on the fully developed laminar velocity field as proposed by Kaviany is accounted for. For various values of the extinction coefficient, the scattering albedo, the conduction-radiation parameter and the boundary emissivity, Nusselt number, temperature and heat flux distributions are found for the range of values including the extreme limits of the porous medium shape parameter (PMSP), γ=(W²φ/K)^1/2, where W is the channel width, φ the porosity and K the permeability. For the lower limiting value of the PMSP γ, the effect of the porous medium is negligible and the situation approaches that of Poiseuille flow. For this limiting case, results from the present work are compared with those available in the literature. For medium to high values of the PMSP γ, for the purpose of comparison, some results are presented in tabular form. Radiation is found to have a significant effect on various parameters studied. The discrete transfer method was used for the solution of the radiative part of the energy equation. An iterative finite difference scheme was used to solve the energy equation.     

Passive Laminar Heat Transfer Across Porous Cavities Using the Thermal Non-Equilibrium Model

This work presents a study on laminar free convection within a square cavity filled with a fluid saturated porous medium. Macroscopic flow equations are obtained by volume-averaging local instantaneous continuity and momentum equations. The so-called “two-energy equation model” is used, in which distinct macroscopic equations are applied to the working fluid and the solid material. Transport equations are discretized using the control-volume method and the system of algebraic equations is relaxed via the SIMPLE (Semi-Implicit Method for Pressure-Linked Equations) algorithm. The effect of Ra_m on Nu_w correctly predicted the enhancement of passive heat transfer across the cavity for increasing Ra_m. Increasing k_s/k_f enhances the conduction transport through the solid material and, consequently, dampens the overall Nusselt number, defined here as the ratio between conduction and convection mechanisms over conduction transport only. Further, results indicate that by increasing the void space within the porous material the overall Nusselt number is reduced rather than increased. Individual contributions to the average Nusselt number indicate that, although convection is enhanced with increasing porosity, the reduction of conduction heat transfer through the solid material is the controlling mechanics for Nu_w as porosity increases. The results herein might contribute to design and optimization of passive heat transfer systems.     

Analysis of Bejan’s heatlines on visualization of heat flow and thermal mixing in tilted square cavities

This article analyzes the detailed heat transfer phenomena during natural convection within tilted square cavities with isothermally cooled walls (BC and DA) and hot wall AB is parallel to the insulated wall CD. A penalty finite element analysis with bi-quadratic elements has been used to investigate the results in terms of streamlines, isotherms and heatlines. The present numerical procedure is performed over a wide range of parameters (10³ ? Ra ? 10⁵,0.015 ? Pr ? 1000,0° ? φ ? 90°). Secondary circulations cells are observed near corner regions of cavity for all φ’s at Pr = 0.015 with Ra = 10⁵. Two asymmetric flow circulation cells are found to occupy the entire cavity for φ = 15° at Pr = 0.7 and Pr = 1000 with Ra = 10⁵. Heatlines indicate that the cavity with inclination angle φ = 15° corresponds to large convective heat transfer from the wall AB to wall DA whereas the heat transfer to wall BC is maximum for φ = 75°. Heat transfer rates along the walls are obtained in terms of local and average Nusselt numbers and they are explained based on gradients of heatfunctions. Average Nusselt number distributions show that heat transfer rate along wall DA is larger for lower inclination angle (φ = 15°) whereas maximum heat transfer rate along wall BC occur for higher inclination angle (φ = 75°).