A numerical study is carried out concerning mixed convection of the nanofluid in two-sided lid-driven square cavity with a pair of triangular heat sources. The upper and bottom moving walls are thermally insulated while the left and right walls are cooled at constant temperature. Two-dimensional Navier-Stokes and energy equations are solved using the finite volume discretization method with SIMPLE algorithm. The method used is validated against previous works. Two cases were considered depending on the direction of moving walls. Effects of various design parameters such as Richardson number (0.1≤Ri≤100), nanoparticle volume fraction (0≤φ≤0.05), and size (25nm≤dp≤145nm) and type (Cu,Al2O3,TiO2) of nanoparticles on the heat transfer rate are investigated. The results of this investigation illustrate that, by reducing the diameter of the nanoparticles and Ri, the heat transfer rate increases. Moreover, it is found that by changing horizontal direction of the moving walls the heat transfer rate variation is negligible.
1. Introduction
Mixed convection problem with lid-driven flows in enclosures occurs in many engineering applications such as heating and cooling nuclear systems of reactors, lubrication technologies, cooling of electronic devices, and heat exchangers. There are many studies concerning the heat transfer in rectangular or square cavities driven by buoyancy and shear which have been studied extensively in the literature. Work of Oztop et al. [1] can be mentioned as example of such studies, in which they numerically studied mixed convection in square cavities with two moving walls. Their results suggest that when the vertical walls move upwards in the same direction, the heat transfer decreases significantly compared to when the vertical walls move in opposite directions. Work of Khanafer and Aithal [2] was concentrated on the effect of mixed convection flow and heat transfer characteristics in a lid-driven cavity with a circular body inside. The results showed that the average Nusselt number increases with an increase in the radius of the cylinder for various Richardson numbers and the optimal heat transfer results are obtained when placing the cylinder near the bottom wall. Islam et al. [3] investigated a numerical study on a lid-driven cavity with a heated square blockage. They found that size, location, and Richardson number of the heater eccentricities affect the average Nusselt number of heater. Kalteh et al. [4] considered laminar mixed convection of nanofluid in a lid-driven square cavity with a triangular heat source. Their simulations indicate that increasing the volume fraction, the nanoparticles diameter, and Reynolds number leads to an increase in average Nusselt number.
El Abdallaoui et al. [5] performed a numerical simulation of natural convection between a decentered triangular heating cylinder and centered triangular heating cylinder [6] in a square outer cylinder filled with a pure fluid or a nanofluid using the lattice Boltzmann method. The results indicate that the horizontal displacement from the centered position to decentered position leads to a considerable increase of heat transfer at weak Rayleigh and the vertical displacement has most important effect on heat transfer at high values of Rayleigh. Oztop et al. [7] investigated mixed convection in a lid-driven cavity with a cylindrical blockage inside. The vertical left sidewall was moving up or down. The computation was carried out for wide ranges of Richardson number, inner cylinder diameter, and center and location of the inner cylinder. They concluded that the most effective parameter on flow field and temperature distribution was the direction of the moving lid.
In the present study, the main aim is to examine the mixed convection heat transfer of nanofluid in a square enclosure with a pair of triangular heat sources. The first case under investigation is characterized by the numerical models used in our study. The computational procedure elaborated in this study is validated against the numerical results of other investigations. We studied the effects of various design parameters such as Richardson number (0.1≤Ri≤100), nanoparticle volume fraction (0≤φ≤0.05), and size (25nm≤dp≤145nm) and type (Cu,Al2O3,TiO2) of nanoparticles on the heat transfer. The new models of the thermal conductivity and effective viscosity investigated by Corcione [8] are used to estimate thermophysical proprieties of the nanofluid. Our numerical results are presented in the form of plots of isotherms, streamlines, and average Nusselt numbers to show the influence of nanofluid and design parameters.
2. Problem Statement
The studied configurations and coordinate system of the considered enclosure in the present study are shown in Figure 1. It is a two-sided lid-driven square cavity with a pair of triangular heating cylinders, filled with nanofluids. The triangular heat source is maintained at a temperature Th=310 K. The vertical walls are cooled at constant temperatures Tc=290 K. Two different cases were considered as shown in Figure 1. In Case 1, the upper wall (adiabatic) is moving to the right while bottom wall (adiabatic) is moving to the left. In Case 2, the direction of the two moving walls is the same. It is assumed that the nanofluid is Newtonian, incompressible, and laminar and the base fluid and the nanoparticles are in a thermal equilibrium state. The thermophysical properties of the nanofluid used in this study are evaluated at the average fluid temperature (Tc+Th)/2 as listed in Table 1.
Thermophysical properties of water and nanoparticles at T = 300 K [14].
Cp (J/KgK)
ρ (Kg/m3)
K (W/mK)
β×105 (K−1)
μ×106 (kgm−1s−1)
dp (nm)
Copper (Cu)
385
8933
401
1.67
—
25, 85, 145
Al2O3
765
3970
36
0.85
—
25
TiO2
710
4157
8.4
0.9
—
25
Water (H2O)
4179
997.1
0.613
27.6
855
0.385
Schematic of the cavity with a pair of triangular heat sources and boundary conditions.
3. Mathematical Formulation
The governing equations including the two-dimensional transient equations of the continuity, momentum, and energy for an incompressible flow are expressed in the following format:(1)∂u∂x+∂v∂y=0,∂u∂t+u∂u∂x+v∂u∂y=-1ρnf∂p∂x+μnfρnf∂2u∂x2+∂2u∂y2,∂v∂t+u∂v∂x+v∂v∂y=-1ρnf∂p∂y+μnfρnf∂2v∂x2+∂2v∂y2+ρβnfρnfT-Tc,∂T∂t+u∂T∂x+v∂T∂y=αnf∂2T∂x2+∂2T∂y2,where the nanofluid effective density, heat capacity, thermal expansion coefficient, and thermal diffusivity are calculated from the following equations [9, 10]:(2)ρnf=1-φρf+φρs,ρCpnf=1-φρCpf+φρCps,ρβnf=1-φρβf+φρβs,αnf=knfρCpnf.Corcione model [8–10] for dynamic viscosity and thermal conductivity of the nanofluid is given by(3)μnf=μf1-34.87dp/df-0.3φ1.03,knfkf=1+4.4ReB0.4Pr0.66TTfr10kpkf0.03φ0.66,ReB=ρfuBdpμf,uB=2kbTπμfdp2.All terms are defined in Nomenclature.
The boundary conditions for mixed convection are written as(4)u=-U0orU0,v=0,∂T∂y=0,onbottomwallofthecavity,u=U0,v=0,∂T∂y=0,onupperwallofthecavity,v=0,v=0,T=Tc,onrightwallofthecavity,u=0,v=0,T=Tc,onleftwallofthecavity.The following dimensionless variables for mixed convection are defined based on properties of pure fluid: (5)τ=tH/Uref,X=xH,Y=yH,U=uUref,V=vUref,P=pρnfUref2,θ=T-TcTh-Tc,where Uref is considered to be U0 for mixed convection. Dimensionless numbers for the system are defined as(6)Re=urefHνnf,Ri=GrRe2,Gr=gβfTh-TcH3νf2,Ra=Gr·Pr=gβfTh-TcH3αfνf,Pr=νfαf.The governing equations (1) are written in the following dimensionless form:(7)∂U∂X+∂V∂Y=0,∂U∂τ+U∂U∂X+V∂U∂Y=-∂P∂X+1Reρfρnfμnfμf∂2U∂X2+∂2U∂Y2,∂V∂τ+U∂V∂X+V∂V∂Y=-∂p∂Y+1Reρfρnfμnfμf∂2V∂X2+∂2V∂Y2+Riρβnfρnfβfθ,∂θ∂τ+U∂θ∂X+V∂θ∂Y=1Re·Prαnfαf∂2θ∂X2+∂2θ∂Y2.Dimensionless form of the boundary conditions can be written as(8)U=1or-1,V=0,∂θ∂y=0,onbottomwallofthecavity,U=1,V=0,∂θ∂y=0,onupperwallofthecavity,V=0,V=0,θ=0,onrightwallofthecavity,U=0,V=0,θ=0,onleftwallofthecavity.The total mean Nusselt number of all cavity’s wall is defined as(9)Nu-tot=1H∫0Hknfφk∂θ∂Xleft+∂θ∂XrightdY.
4. Numerical Details
The discretization procedure of the governing equations (7) and boundary conditions described by (8) is based on a finite volume formulation, given by Patankar [11] on a staggered grid. SIMPLE (Semi-Implicit Method for Pressure Linked Equations) is used to solve the coupled pressure-velocity equation while Hybrid Differencing Scheme (HDS) of Spalding [12] is used for the convective terms. Line by line application of TDMA (Tridiagonal Matrix Algorithm) method [12] is applied on equation systems until sum of the residuals became less than 10-6. The developed algorithm was implemented in FORTRAN program.
4.1. Grid Independence Study
In order to determine a proper grid for the numerical simulation, a square cavity filled with Cu-water nanofluid φ=0.05 having a pair of triangular heat sources with size w=0.2 is analyzed in two extreme Richardson numbers (Ra=0.1and100). The mean Nusselt number obtained using different grid numbers for particular cases is presented in Table 2. As can be observed from the table, a uniform 103×103 grid is sufficiently fine for the numerical calculation.
Effect of the grid size on Nu-tot for the cavity filled with the Cu-water nanofluid φ=0.05 having a pair of triangular heat sources with size W=w/H=0.2 for case 1 (see Figure 1).
Ri
63×63
83×83
103×103
123×123
0.1
19.781
20.295
20.351
20.357
100
8.921
9.377
9.438
9.439
4.2. Validations
The present numerical scheme was validated against various numerical results available in the literature, two different heat convection problems are chosen. The first case is the numerical results of Iwatsu et al. [13] and Oztop et al. [7] for a top heated moving lid and bottom cooled square cavity filled with air (Pr=0.71). A 100×100 mesh was used and the computations were done for three different Richardson numbers. Table 3 demonstrates an excellent comparison of the average Nusselt number between the present results and the numerical results found in the literature [7–13] with a maximum discrepancy of about 1.6%. The second case is a mixed convection flow and heat transfer characteristics in a lid-driven cavity with a square heater inside. Figure 2 illustrates a comparison of the isotherms and streamlines between the present results and the results reported by Islam et al. [3] at Reynolds number (Re=100) and at Richardson numbers (Ri=0.1).
Comparison of Nu- at the hot lid between the present results and those reported in the literature.
Ri
Average Nusselt number at the hot lid
Iwatsu et al. [13]
Oztop et al. [7]
Present
1.0
1.34
1.30
1.36
0.0625
3.62
3.63
3.68
0.01
6.29
6.34
6.29
Comparison of the streamlines and isotherms between the present results and that of Islam et al. [3] (h/H=0.25,Re=100,andRi=0.1).
5. Results and Discussion
In this section, the numerical results for mixed convection heat transfer of the nanofluid in a square cavity with horizontal moving lids are discussed. Results pertinent to the effects of the moving wall direction of motion, volume fraction (0≤φ≤0.05), diameter (25nm≤dp≤145nm), and type Cu,Al2O3,andTiO2 of the nanoparticles on the fluid flow and heat transfer rate are presented for the square cavity with a pair of triangular heat sources placed horizontally (see Figure 1, Cases 1 and 2). Size of triangular heat sources is fixed at w=0.2H.
Figures 3(a) and 3(b) show the effects of the moving wall direction of motion on the streamlines and isotherms for various values of Ri ranging from 0.1 to 100 while Grashof number is fixed at Gr=104. The cavity is filled with Cu-water φ=5%, and for comparisons the streamlines and the isotherms for pure fluid and nanofluid are shown by dashed line and solid line, respectively. As shown in Figure 3, some differences are observed in streamlines and isotherms of pure fluid and nanofluid. We can explain this by the higher viscosity of nanofluid compared to that of the pure fluid which increases the diffusion of momentum in the nanofluid. Figure 3(a) shows that, at Ri=0.1 for Case 1, the flow field is characterized by a clockwise vortex created in the upper half of the enclosure, and another created in the lower half. These two vortexes are mostly generated by horizontal moving lids. At this Ri, we can see in Figure 3(b) that the concentration of the isotherms, close to the pair of triangular heat sources, increases, which indicates that the heat transfer is through forced convection. By increasing Ri, the effects of natural convection and buoyancy force become stronger leading to formation of a counterclockwise rotating eddy close to the left wall of the cavity. The effects of moving wall decrease by increasing Ri. The uniformly distributed isotherms at high Ri show that the main heat transfer mechanism is through the natural convection.
(a) Streamlines and (b) isotherms inside the cavity with a pair of triangular heat sources (W=w/H=0.2) filled with Cu-water pure fluid (dashed line) and Cu-water nanofluid (solid line) with φ=5% and at different Ris for Cases 1 and 2. Gr=104.
Figure 3(a) shows that, at low Ri for Case 2, the flow field is characterized also by a clockwise vortex created in the upper half of the enclosure, and a counterclockwise vortex created in the lower half. These two vortices move to the right side due to an increase in lids speed. Figure 3(b) shows that the concentration of the isotherms, close to the pair of triangular heat sources and left wall of the cavity, increases. The isotherms move to the left side and the diffusion of heat to the right side decreases with increasing the speed of the horizontal lids (decreasing Ri). This observation shows the higher forced convection effects in the right side of the cavity. When Ri increases to 100, the effects of moving wall decrease due to the enhancement in the buoyancy force. A secondary eddy appears on the left indicates that the natural convection becomes more dominant there.
Figure 4 shows values of Nu-tot for the illustrated cases in Figure 3. According to Figure 4(b), it is clear that the heat transfer rate increases with decreasing the nanoparticle diameter and the highest values of the heat transfer rate occur at 25nm diameter. It is also noticed that by increasing volume fraction of nanoparticles and decreasing Ri the heat transfer rate increases. At Ri>10 the effects of natural convection are dominant so there is an optimum volume fraction of nanoparticles which maximize the heat transfer which is about 3% in the most cases.
Variations Nu-tot with respect to the volume fraction of the nanoparticles at different Richardson numbers, diameters, and type of nanoparticles for (a) Case 1 and (b) Case 2. Gr=104.
Figure 4(a) shows the effects of using different types of nanofluids on Nu-tot. In all Richardson numbers nonsignificant variation on the heat transfer rate is observed when different kinds of nanofluids are used. It can be seen that among all nanoparticles those nanoparticles with higher thermal conductivity (such as Cu) produce slightly higher rate of heat transfer compared to nanoparticles with lower thermal conductivity (such as TiO2). It is clear from Figure 4 that, at all Ri, if the direction of the two moving walls is same, the heat transfer rate increases. From Cases 2 to 1, the achieved reduction of Nu-tot is approximately 4.67%, 0.28%, 0.09%, and 0.44%, at Ri=0.1,1,10,and100, respectively. Consequently, the heat transfer rate variation between these two cases is negligible.
6. Conclusion
This study investigates mixed convection of the nanofluid in two-sided lid-driven square cavity with a pair of triangular heating cylinders. The effects of various design parameters such as Richardson number (0.1≤Ri≤100), nanoparticle volume fraction (0≤φ≤0.05), and size (25nm≤dp≤145nm) and type (Cu,Al2O3,TiO2) of nanoparticles on the heat transfer rate are investigated. According to the presented results, the following conclusions are drawn:
By increasing volume fraction of nanoparticles (φ) and increasing Richardson number (Ri), the heat transfer rate increases.
The heat transfer rate increases with decreasing the nanoparticle diameter and the highest values of the heat transfer rate occur at 25 nm diameter.
In all Richardson numbers nonsignificant variation on the heat transfer rate is observed when different kinds of nanofluids are used. It can be seen that among all nanoparticles those nanoparticles with higher thermal conductivity (such as Cu) produce slightly higher rate of heat transfer compared to nanoparticles with lower thermal conductivity (such as TiO2).
At Ri>10 the effects of natural convection are dominant so there is an optimum volume fraction of nanoparticles which maximize the heat transfer which is about 3% in the most cases.
The effect of changing direction of the two horizontal moving walls on the heat transfer rate (Nu-tot) is negligible.
NomenclatureCp:
Specific heat, Jkg-1K-1
dp:
Diameter of the nanoparticle, m
df:
Diameter of the base fluid molecule, m
g:
Gravitational acceleration, ms-2
Gr:
Grashof number (=gβΔTH3/ν2)
H:
Enclosure height, m
k:
Thermal conductivity, Wm-1K-1
kb:
Boltzmann’s constant =1.38066×10-23
Nu-tot:
Heat transfer of cavity’s wall
p:
Dimensional pressure, Nm-2
P:
Dimensionless pressure
Pr:
Prandtl number =νf/αf
Ra:
Rayleigh number =gβfTh-TcH3/αfνf
ReB:
Brownian-motion Reynolds number
Re:
Reynolds number =U0H/ν
Ri:
Richardson number =Gr/Re2
T:
Temperature, K
t:
Time, s
τ:
Dimensionless time (t/H/Uref)
u,v:
Velocity components, ms-1
uB:
Brownian velocity of the nanoparticle, ms-1
U,V:
Dimensionless velocity components
w:
Dimensional width of the triangular heater
W:
Dimensionless width of the triangular heater
x,y:
Cartesian coordinates, m
X,Y:
Dimensionless Cartesian coordinates (x,y)/H.
Greek Symbolsα:
Thermal diffusivity, m2s-1
β:
Thermal expansion coefficient, K-1
θ:
Dimensionless temperature
μ:
Dynamic viscosity, kgm-1s-1
ν:
Kinematic viscosity, m2s-1
ρ:
Density, kgm-3
φ:
Volume fraction of the nanoparticles
ψ:
Stream function (=-∫YoYU∂Y+ψ(X,Y0)).
Subscriptsc:
Cold
f:
Fluid
h:
Hot
nf:
Nanofluid
p:
Solid nanoparticles.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
OztopH. F.DagtekinI.BahloulA.Comparison of position of a heated thin plate located in a cavity for natural convection200431112113210.1016/s0735-1933(03)00207-02-s2.0-0347051801KhanaferK.AithalS. M.Laminar mixed convection flow and heat transfer characteristics in a lid driven cavity with a circular cylinder20136620020910.1016/j.ijheatmasstransfer.2013.07.0232-s2.0-84881281929IslamA. W.SharifM. A. R.CarlsonE. S.Mixed convection in a lid driven square cavity with an isothermally heated square blockage inside20125519-205244525510.1016/j.ijheatmasstransfer.2012.05.0322-s2.0-84863526039KaltehM.JavaherdehK.AzarbarzinT.Numerical solution of nanofluid mixed convection heat transfer in a lid-driven square cavity with a triangular heat source201425378078810.1016/j.powtec.2013.12.0392-s2.0-84892663300El AbdallaouiM.HasnaouiM.AmahmidA.Numerical simulation of natural convection between a decentered triangular heating cylinder and a square outer cylinder filled with a pure fluid or a nanofluid using the lattice Boltzmann method201527193205El AbdallaouiM.HasnaouiM.AmahmidA.Lattice-boltzmann modeling of natural convection between a square outer cylinder and an inner isosceles triangular heating body20146691076109610.1080/10407782.2014.8943982-s2.0-84904469310OztopH. F.ZhaoZ.YuB.Fluid flow due to combined convection in lid-driven enclosure having a circular body200930588690110.1016/j.ijheatfluidflow.2009.04.0092-s2.0-70249120325CorcioneM.Empirical correlating equations for predicting the effective thermal conductivity and dynamic viscosity of nanofluids201152178979310.1016/j.enconman.2010.06.0722-s2.0-79961179051CorcioneM.Heat transfer features of buoyancy-driven nanofluids inside rectangular enclosures differentially heated at the sidewalls20104991536154610.1016/j.ijthermalsci.2010.05.0052-s2.0-77955412769HaddadZ.OztopH. F.Abu-NadaE.MataouiA.A review on natural convective heat transfer of nanofluids20121675363537810.1016/j.rser.2012.04.0032-s2.0-84863437157PatankarS. V.1980Washington, DC, USAMcGraw-HillSpaldingD. B.A novel finite difference formulation for differential expressions involving both first and second derivatives19724455155910.1002/nme.16200404092-s2.0-0015368687IwatsuR.HyunJ. M.KuwaharaK.Mixed convection in a driven cavity with a stable vertical temperature gradient19933661601160810.1016/s0017-9310(05)80069-92-s2.0-0027572792IncroperaF. P.DeWittD. P.2002New York, NY, USAJohn Wiley & Sons