Abstract
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 , nanoparticle volume fraction , and size and type of nanoparticles on the heat transfer rate are investigated. The results of this investigation illustrate that, by reducing the diameter of the nanoparticles and , 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 , nanoparticle volume fraction , and size and type 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 K. The vertical walls are cooled at constant temperatures 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 as listed in Table 1.
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:where the nanofluid effective density, heat capacity, thermal expansion coefficient, and thermal diffusivity are calculated from the following equations [9, 10]:Corcione model [8–10] for dynamic viscosity and thermal conductivity of the nanofluid is given byAll terms are defined in Nomenclature.
The boundary conditions for mixed convection are written asThe following dimensionless variables for mixed convection are defined based on properties of pure fluid: where is considered to be for mixed convection. Dimensionless numbers for the system are defined asThe governing equations (1) are written in the following dimensionless form:Dimensionless form of the boundary conditions can be written asThe total mean Nusselt number of all cavity’s wall is defined as
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 . 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 having a pair of triangular heat sources with size is analyzed in two extreme Richardson numbers . 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 grid is sufficiently fine for the numerical calculation.
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 . A 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 . 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 ( and at Richardson numbers (.
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 , diameter , and type 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 and ). Size of triangular heat sources is fixed at .
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 . The cavity is filled with Cu-water , 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 for Case , 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 , 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 , 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 . The uniformly distributed isotherms at high show that the main heat transfer mechanism is through the natural convection.
Figure 3(a) shows that, at low for Case , 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 ). This observation shows the higher forced convection effects in the right side of the cavity. When 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 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 diameter. It is also noticed that by increasing volume fraction of nanoparticles and decreasing the heat transfer rate increases. At 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.
(a)
(b)
Figure 4(a) shows the effects of using different types of nanofluids on . 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 ). It is clear from Figure 4 that, at all , if the direction of the two moving walls is same, the heat transfer rate increases. From Cases to 1, the achieved reduction of is approximately 4.67%, 0.28%, 0.09%, and 0.44%, at , 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 , nanoparticle volume fraction , and size and type of nanoparticles on the heat transfer rate are investigated. According to the presented results, the following conclusions are drawn:(i)By increasing volume fraction of nanoparticles and increasing Richardson number , the heat transfer rate increases.(ii)The heat transfer rate increases with decreasing the nanoparticle diameter and the highest values of the heat transfer rate occur at 25 nm diameter.(iii)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 ).(iv)At 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.(v)The effect of changing direction of the two horizontal moving walls on the heat transfer rate is negligible.
Nomenclature
| : | Specific heat, |
| : | Diameter of the nanoparticle, m |
| : | Diameter of the base fluid molecule, m |
| Gravitational acceleration, | |
| Grashof number | |
| Enclosure height, m | |
| Thermal conductivity, | |
| : | Boltzmann’s constant |
| : | Heat transfer of cavity’s wall |
| Dimensional pressure, | |
| Dimensionless pressure | |
| Prandtl number | |
| Rayleigh number | |
| : | Brownian-motion Reynolds number |
| Reynolds number | |
| Richardson number | |
| Temperature, K | |
| Time, s | |
| Dimensionless time | |
| Velocity components, | |
| : | Brownian velocity of the nanoparticle, |
| Dimensionless velocity components | |
| Dimensional width of the triangular heater | |
| Dimensionless width of the triangular heater | |
| Cartesian coordinates, m | |
| Dimensionless Cartesian coordinates . |
| Thermal diffusivity, | |
| Thermal expansion coefficient, | |
| Dimensionless temperature | |
| Dynamic viscosity, | |
| Kinematic viscosity, | |
| Density, | |
| Volume fraction of the nanoparticles | |
| Stream function . |
| Cold | |
| Fluid | |
| Hot | |
| Nanofluid | |
| Solid nanoparticles. |
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.