Document Type : Full Lenght Research Article
Authors
Semnan University
Abstract
Keywords
Main Subjects
Suspending nanoparticles in a host fluid will result in a new fluid that is called nanofluid. Water, Ethylene glycol and propylene glycol are from famous host fluids. In 1995 Choi et al. introduced new age of fluids that now is known as nanofluid. Existence of high thermal conductivity metallic nanoparticles (e.g., copper, aluminum, silver and Titanium) increases the thermal conductivity of such mixtures, thus enhancing their overall heat transfer capability (Xuan et al. 2003). In recent years, Nanofluids have attracted attention as a new innovation of heat transfer fluids in building heating, in various heat exchangers, in plants and in automotive cooling applications, because of their excellent thermal performance. Various benefits of the application of nanofluids include: improved heat transfer, heat transfer system size reduction, minimal clogging, micro channel cooling and miniaturization of systems (Choi, 1995).
Numerous investigates have been conducted on the thermophysical properties of nanofluids (effective dynamic viscosity, thermal conductivity and etc.) and the energy transport in nanofluids. Study of thermophysical properties of nanofluids can be observed in several literatures such as Lee et al. (2000), Xie et al. (2002), Patel et al. (2005), and Chang et al. (2005), Also many theoretical, numerical simulations and empirical studies on influence of existence of nanoparticle in convective heat transfer have been reported.
On the other hand, fluid flow and heat transfer in a cavity filled by pure fluid which is driven by buoyancy and shear have been studied extensively in literature. Mixed Convection (a kind of convection Including natural convection heat transfer and forced convection heat transfer) has significant role in many applications in industry and engineering. lake and reservoirs (Imberger and Hamblin, 1982), food processing, crystal growth (Moallemi and Jang, 1992), electronic cooling devices, drying technologies, solar ponds (Cha and Jaluria, 1984) solar collectors(Ideriah, 1980) and float glass production (Pilkington, 1969), are among its current applications.
Effect of several parameters on mixed convection of various cavities are investigated in many studies like Talebi et al. (2010), Abunada et al.(2010), Mahmoodi(2011), Sadodin et al.(2011), Hemmat Esfe et al.(2012), Abbasian Arani et al.(2012), Hemmat Esfe (2012), Fereidoon et al.(2013), Saedodin et al.(2013) and Zarei et al.(2013). They studied effects of parameters like cavity geometry, base fluid and effect of boundary condition etc.
Nikfar and Mahmoodi (2012) simulated natural convection in a square cavity containing water based nanofluid. They stated that Nusselt number of cavities hot wall has direct relationship with solid volume fraction.
In another study Mahmoodi and Mazrouei(2012) studied the effect of putting adiabatic obstacles in cavities on fluid flow and natural heat transfer. The working fluid in their study was a water based fluid containing Cu nanoparticles. They stated that at low Rayleigh numbers, the rate of heat transfer reduces by increasing the size of adiabatic obstacles.
Mahmoodi (2012) studied mixed convection in a lid-driven cavity filled with Al2O3-water nanofluid. Having a moving top wall in their simulation is one of their novelties in their study. In another research Mazrouei and Sebdani et al.(2012) checked the variable properties of nanofluids effect on mixed convection in a rectangular cavity.
In present study an inclined square cavity is simulated in order to study the mixed convection of flow in a cavity filled with Al2O3 nanoparticles and also with a heated obstacle and upper moving walls.
Based on the authors' knowledge no research or paper has been done or publish in cavity with Simultaneous changes in parameters such as diameter of nanoparticles and inclination angles. The effects of changes in important parameters such as Richardson number, inclination angle, diameter of nanoparticles and solid volume concentration of the Al2O3 nanoparticles on the flow and thermal fields and heat transfer inside the cavity with heated obstacle has not been reported so far. Utilizing of new variable properties models for estimate the thermal conductivity and dynamic viscosity of nanofluid are other aspects that distinguished this survey from other numerical studies in this domain.
In Fig. 1 a schematic of studied geometry in present study can be seen. As it stated the cavity contains a water based nanofluid with Al2O3 nanopartilces. Adiabatic condition is considered for left, top and bottom sides of cavity. For the right wall of cavity alow temperature like Tc is considered. Also Th is determined as temperature of obstacle in order to induce the buoyancy effect.
The nanofluid behavior used in cavity is Newtonian, Incompressible and the flow is laminar. Also it is considered that the host fluid and nanoparticles as solid phase, flow with the same velocity in cavity. The density variation in the body force term of the momentum equation is satisfied by Boussinesq’s approximation. The thermophysical properties of nanoparticles and the water as the base fluid at T = 25˚C are presented in Table1.
The thermal conductivity and the viscosity of the nanofluid are taken into consideration as variable properties; both of them change with volume fraction and temperature of nanoparticles. Under the above assumptions, the system of governing equations is:
(1)
Table 1. Thermophysical properties of all ingredients of nanofluid at T =25˚C.
Physical properties |
Fluid phase (Water) |
Solid (Al2O3) |
Cp(J/kg k) |
4179 |
765 |
(kg/m3) |
997.1 |
3970 |
K (W m-1 K-1) |
0.6 |
25 |
β×10-5 (1/K) |
21. |
0.85 |
μ×10-4(Kg/ms) |
8.9 |
…….. |
dp (nanometer) |
--- |
47 |
Fig. 1. Schematic diagram of current study
(2)
(3)
and
(4)
The dimensionless parameters may be presented as
(5)
Hence,
(6)
The dimensionless form of the above governing equations (1) to (4) become
(7)
(8)
(9)
and
(10)
2.1 Thermal diffusivity and effective density
Thermal diffusivity and effective density of the nanofluid are
(11)
(12)
2.2 Heat capacity and thermal expansion coefficient
Heat capacity and thermal expansion coefficient of the nanofluid are therefore
(13)
(14)
2.3 Viscosity
The effective viscosity of nanofluid was calculated by:
(15)
This well-validated model is presented by Jang et al. (2007) for a fluid containing a dilute suspension of small rigid spherical particles and it accounts for the slip mechanism in nanofluids. The empirical constant and are -0.25 and 280 for Al2O3, respectively.
It is worth mentioning that the viscosity of the base fluid (water) is considered to vary with temperature and the flowing equation is used to evaluate the viscosity of water,
where (16)
2.4 Dimensionless stagnant thermal conductivity:
Hamilton and crosser (H-C model)(1962), is used for calculating effective thermal conductivity of suspended nanoparticles in a host fluid.
(17)
2.5 Total dimensionless thermal conductivity of nanofluids:
(18)
The presented model is achieved by results of Xu et al. (2006) study and has been chosen to reach the thermal conductivity of nanofluids. In presented formula “c” is an experimental constant that is a certain number for each host fluid. It is 85 for deionized water and is 280 for ethylene glycol. This constant is not related to type of nanoparticles.
Nup is the Nusselt number for liquid flowing around a spherical particle and equal to two for a single particle. The fluid molecular diameter is df =4.5*10-10(m) for water in present study and the fractal dimension Df is determined by:
(19)
where dp,max and dp,min are the maximum and minimum diameters of nanoparticles, respectively. Ratio of minimum to maximum nanoparticles dp,min/dp,max is R.
(20)
3. Numerical approach
Governing equations for continuity, momentum and energy equations associated with the boundary conditions in this investigation were calculated numerically based on the finite volume method and associated staggered grid system, using FORTRAN computer code. The SIMPLER algorithm is used to solve the coupled system of governing equations. The convection terms is approximated by a hybrid-scheme which is conducive to a stable solution. In addition, a second-order central differencing scheme is utilized for the diffusion terms. The algebraic system resulting from numerical discretization was calculated utilizing TDMA applied in a line going through all volumes in the computational domain. The solution process is repeated until an acceptable convergence criterion is reached. A FORTRAN computer code has been developed to solve the equations as described above. The process is repeated until the following convergence criterion is satisfied:
Fig. 2. Mesh valid at X=0.6L, φ=0.05, γ=90˚
(21)
In above equation, M and N are number of grid points in x and y directions, respectively. n shows number of iteration and is representative of any scalar transport quantity. Nine different mesh sizes were used to check and qualify the grid independency.
As can be observed, 101 × 101 uniform grid size yields the required accuracy and was hence applied for all simulation exercises in this work as presented in the following section.
To ensure the accuracy and validity of this new model, we analyze a square cavity filled with base fluid with Pr = 0.7 and different Ra numbers. Table 2 shows the comparison between the results obtained with our research and the values presented in the literature. The quantitative comparisons for the average Nusselt numbers indicate an excellent agreement between them.
Table 2. Validation analysis results |
|||||
Ha et al. |
Fusegi et al. |
Tiwari and Das |
Hadjisophocleous |
Present study |
|
|
|
|
|
|
|
|
|
|
|
|
Ra=103 |
1.072 |
1.085 |
1.0871 |
1.141 |
1.1366 |
Nu |
|
|
1.508 |
1.540 |
1.5738 |
Numax |
|
|
0.6901 |
0.727 |
0.6468 |
Numin |
|
|
|
|
|
|
|
|
|
|
|
Ra=104 |
2.070 |
2.1 |
2.195 |
2.29 |
2.3256 |
Nu |
|
|
3.5585 |
3.84 |
3.5456 |
Numax |
|
|
0.5809 |
0.670 |
0.7172 |
Numin |
|
|
|
|
|
|
|
|
|
|
|
Ra=105 |
4.464 |
4.361 |
4.450 |
4.964 |
4.4928 |
Nu |
|
|
7.9371 |
8.93 |
7.3037 |
Numax |
|
|
0.7173 |
1.01 |
0.9906 |
Numin |
|
|
|
|
|
|
|
|
|
|
|
Ra=106 |
|
|
8.803 |
10.39 |
8.6388 |
Nu |
|
|
19.2675 |
21.41 |
14.2521 |
Numax |
|
|
0.9420 |
1.58 |
1.6247 |
Numin |
4. Results and discussion
In this study, thermal characteristics and flow patterns inside an inclined cavity filled with nanofluid with hot barrier are investigated. Some parameters such as diameter of nanoparticles, cavity inclination angles, solid volume fraction and Richardson number are considered and their effects on streamlines, isotherms and total heat transfer are studied. Grashof number is assumed constant and equal to 104 while
Figure 3 shows flow and temperature behavior in different cavity inclination angles at. Flow pattern for horizontal situation of the cavity shows presence of a clockwise primary cell in upper parts of the cavity while one small vortex is formed at right side of hot barrier. Presence of vortexes results from two buoyancy force due to temperature difference and shear force due to upper lid movement. In horizontal situation, buoyancy and shear forces assist each other. With increasing slope of this cavity, buoyancy acts in a direction opposite to shear forces and formation of a small vortex near upper lid proves this fact. It is expected that increasing cavity inclination angles to 90oresults in relative neutralization of buoyancy forces and shear forces and consequently vortex forces and heat transfer inside the cavity decrease. Isotherm lines in horizontal position of the cavity also show density of lines near isotherm walls. With increasing slope of the cavity, density of isotherm lines and thereupon temperature gradient near walls decrease and this decrease is expected to reduce heat transfer inside the cavity.
Figure 4 portrays streamlines and isotherms to demonstrate the effect of nanoparticles diameter for Ri=1, T=300 and . 4 different diameters (20, 40, 60 and 80 nano meter) are shown in this figure. As it is observed, changes in nanoparticles diameter produces no substantial changes in flow patterns and temperature. Flow pattern in this parameter range shows formation of two clockwise and anti-clockwise vortexes which upper clockwise vortex is stronger than lower and occupies more spaces of the cavity. The strength of vortexes slightly decreases with increasing nanoparticles diameter but their core and primary form is maintained. Temperature lines also demonstrate formation of a thermal boundary layer inside the cavity and near the barrier. Isotherm lines are slightly separated from each other with increasing diameter of nanoparticles and therefore, temperature gradient decreases. According to isotherm lines it is expected that increasing diameter of nanoparticles causes total heat transfer inside the cavity to decrease slightly.
Variation of isotherms and streamlines versus volume fraction of nanoparticle inside the cavity are demonstrated in figure 5 at Ri=0.01, dp=40, T=300, . Streamlines show formation of a strong vortex in upper and middle parts of the cavity while two strong and weak vortexes can be observed at sides of hot rectangular barrier. No significant changes occur in flow pattern with increasing volume fraction of nanoparticles. Temperature lines in this range of parameters present intense density of isotherm lines near isothermal walls. Density of lines decreases very slightly with increasing volume fraction and temperature gradient reduces. Despite a reduction in temperature gradient with increasing volume fraction, no accurate prediction can be made for total heat transfer inside the cavity. Increasing solid concentration increases thermal conductivity of nanofluid while it slightly decreases temperature gradient.
|
||
Fig. 3. streamlines and isotherms in different inclination angles in Ri=100, φ=0.05%, dp=20 nm
|
dp=20 nanometer |
||
dp=40 nanometer |
||
dp=60 nanometer |
||
dp=80 nanometer |
||
Fig. 4. streamlines and isotherms in different diameter of nanoparticles at Ri=1, φ=0.05%, T=300, γ=30 |
Fig. 5. streamlines and isotherms in different solid volume concentration at Ri=0.01, dp=40, T=300, |
Figure 6 illustrates heat transfer inside the cavity versus changes of slope and Richardson number for dp=20, T=300, . As it is obvious in this diagram, with increasing Richardson number and (buoyancy forces dominate shear force) Nusselt number and consequently heat transfer inside the cavity decrease. On the other hand, with increasing cavity inclination angles from horizontal position, buoyancy force and shear force counteract each other and as a result total heat transfer inside the cavity decreases. This matter was discussed earlier when figure 2 was studied.
Fig. 6. Variation of Nusselt number versus cavity inclination angle and Richardson number for dp=20, T=300, .
Fig. 7. Variation of Nusselt number versus Richardson number for different diameters of nanoparticles dispersed in water for T=300o, and %
Fig. 8. Variation of Nusselt number versus Richardson number and volume fraction of nanoparticles at dp=40, T=300o and
Figure 7 shows variation of Nusselt number versus Richardson number for different diameters of nanoparticles dispersed in water for T=300o, and . At this condition, values of Nusselt number for all diameters of nanoparticles decrease with increasing Richardson number. With increasing nanoparticles diameter from 20 nm to 40 nm, a relatively significant change occurs in heat transfer while with increasing diameter more than 40 nm, no significant changes occur in heat transfer and Nusselt number.
Figure 8 illustrates variation of Nusselt number versus changes in Richardson number and volume fraction of nanoparticles at dp=40, T=300o and . As it is seen in this diagram, adding nanoparticles to base fluid causes a considerable increase in heat transfer inside the cavity and this increasing trend continues with increasing volume fraction. On the other side, in this range of parameters increasing Richardson number (i.e. buoyancy force dominates shear force) also results in a reduction in Nusselt number and consequently heat transfer inside the cavity.
In this paper thermal behavior of nanoparticles inside an inclined cavity with a hot barrier was investigated. Effects of some important parameters including nanoparticle diameter, cavity inclination angles, volume fraction and Richardson number were studied and below results were obtained:
1- With increasing Richardson number and predominance of buoyancy force over shear force, Nusselt number and heat transfer decrease.
2- In horizontal situation that buoyancy and shear forces have the same direction, the maximum values are obtained for Nusselt number and heat transfer.
3- With increasing slope of the cavity to 90 degree, heat transfer continuously decreases at all studied Richardson numbers.
4- With decreasing diameter of dispersed nanoparticles in water, heat transfer increases.
5- Adding nanoparticles to base fluid in addition to maintaining flow pattern causes a significant increase in heat transfer inside the cavity.
Greek symbols |
Nomenclature
|
|
thermal diffusivity, m2 s |
α |
specific heat, J kg-1 K-1 |
thermal expansion coefficient, K-1 |
β |
Grashof number |
dimensionless temperature |
θ |
gravitational acceleration, m s-2 |
dynamic viscosity, Kg m-1 s-1 |
μ |
heat transfer coefficient, W m-2 K-1 |
kinematic viscosity, m2 s-1 |
ν |
enclosure length, m |
density, kg m-3 |
ρ |
thermal conductivity, W m-1 K-1 |
volume fraction of the nanoparticles |
φ |
Nusselt number |
|
|
pressure, N m-2 |
|
|
dimensionless pressure |
Subscripts |
Prandtl number |
|
cold |
c |
heat flux, W m-2 |
effective |
eff |
Reynolds number |
fluid |
f |
Richardson number |
hot |
h |
dimensional temperature, K |
nanofluid |
nf |
dimensional velocities components in x and y direction, m s-1 |
solid particles |
s |
dimensionless velocities components in X and Y direction |
wall |
w |
lid velocity |
|
|
dimensional Cartesian coordinates, m |
|
|
dimensionless Cartesian coordinates |