- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents

Advances in Mechanical Engineering

Volume 2014 (2014), Article ID 392610, 15 pages

http://dx.doi.org/10.1155/2014/392610

## Numerical Study of Natural Convection in Vertical Enclosures Utilizing Nanofluid

^{1}University of Applied Science and Technology, P.O. Box 14155-1622, Tehran, Iran^{2}School of Mechanical Engineering, Babol University of Technology, P.O. Box 484, Babol, Iran^{3}Department of Mechanical Engineering, Payame Noor University, P.O. Box 19395-3697, Tehran, Iran

Received 23 June 2013; Accepted 24 October 2013; Published 27 January 2014

Academic Editor: Kambiz Vafai

Copyright © 2014 M. Alipanah et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Enhancement of buoyancy-driven convection heat transfer within vertical cavities containing nanofluids subjected to different side wall temperatures and various aspect ratios is investigated. The computations are based on an iterative, finitevolume numerical procedure (SIMPLE) that incorporates the Boussinesq approximation to simulate the buoyancy term. With the base fluid being water, three different nanoparticles (Cu, TiO_{2}, and Al_{2}O_{3}) are considered as the nanofluids. This study has been carried out for the pertinent parameters in the following ranges: the Rayleigh number, = 10^{5}–10^{7} and the volumetric fraction of nanoparticle between 0 and 5 percent. The results are presented for different length-to-height ratios varying from 0.1 to 1.0. The comparisons show that the mean Nusselt numbers and velocity magnitudes increase with volume fraction for the whole range of the Rayleigh numbers. The predictions show a noticeable heat transfer enhancement compared to pure fluid. It is also found that the heat transfer enhancement utilizing nanofluid is more pronounced at low aspect ratios than high aspect ratios. Moreover, the results depict that the addition of nanoparticles to the pure fluid has more effects at lower Rayleigh numbers.

#### 1. Introduction

Natural convection is frequently encountered in various engineering applications such as cooling systems for electronic devices [1, 2], chemical vapor deposition instruments (CVD) [3], furnace engineering [4], solar energy collectors and building energy systems [5], non-Newtonian chemical processes [6, 7], and domains affected by electromagnetic fields [8, 9]. Heat transfer enhancement in these systems is an essential topic from an energy saving perspective. An innovative technique to improve heat transfer is by suspending nanoscale particles in the base fluid. Fluids with nanoparticles suspended in them are called nanofluids. Nanofluids introduce a unique opportunity for realizing more effective heat removal in thermal-fluid systems. Adoption of these fluids in a variety of processes and industries is anticipated since materials with sizes of nanometers possess unique physical and chemical properties. To date, theoretical, numerical, and experimental studies on nanofluids including thermal conductivity modeling [10, 11], viscosity estimation [12, 13] and boiling heat transfer and natural convection [14, 15] have appeared. As a pioneer, Masuda et al. [16] reported on thermal conductivity improvement of dispersed ultrafine (nanosize) particles in liquids. Soon thereafter, Choi [17] was the first to coin the term “nanofluids” for this new class of fluids with superior thermal properties. Khanafer et al. [18] investigated heat transfer enhancement in a two-dimensional enclosure utilizing nanofluids for a range of Grashof numbers and volume fractions. It was found that the heat transfer across the enclosure was found to increase with the volumetric fraction of the copper nanoparticles in water at any given Grashof number. Ho et al. [19] investigated the influences of uncertainties due to adopting various relations for the effective thermal conductivity and dynamic viscosity of alumina/water nanofluid on the heat transfer characteristics. It was found that the uncertainties associated with different formulas adopted for the effective thermal conductivity and dynamic viscosity of the nanofluid have a strong bearing on natural convection heat transfer characteristics in the enclosure. Nnanna [20] performed an experimental investigation on the heat transfer behavior of buoyancy-driven Al_{2}O_{3}-water nanofluid in an enclosure. It was found that natural convection heat transfer was enhanced at low volume fractions in the range of . However, for volume fractions above 0.2, natural convective heat transfer was decreased due to reduction in the Rayleigh number caused by an increase in kinematic viscosity. Ho and Lin [21] performed an experimental investigation of natural convection heat transfer of Al_{2}O_{3}-water nanofluid in vertical square enclosures of three different sizes. It was found that the systematic heat transfer degradation for the nanofluids containing nanoparticles of over the entire range of the Rayleigh number. They also reported that for the nanofluid containing much lower particle fraction of 0.001, heat transfer enhancement of around 18% compared with that of water was found to arise in the largest enclosure at sufficiently high Rayleigh number. Abu-Nada et al. [22] investigated the influences of nanoparticle on natural convection heat transfer enhancement in horizontal annuli with various nanoparticles and volume fractions, reporting an enhancement of heat transfer. Wang and Mujumdar [23] covered fluid flow and heat transfer characteristics of nanofluids in forced and free convection flows and potential applications of nanofluids. Khodadadi and Hosseinizadeh [24] investigated nanoparticles within conventional phase change materials such as water. Their findings show that nanoparticle-enhanced phase change materials (NEPCM) have a great potential for demanding thermalenergy storage applications. Although heat transfer predictions for pure fluids in closed cavities have been widely studied in the past, there has been little attempt to report on cases of thin cavities using nanofluids. Wang et al. [25] investigated free convection heat transfer of water/Al_{2}O_{3} nanofluids in horizontal and vertical rectangular enclosures. They reported that the ratio of heat transfer coefficient of nanofluids to that of base fluid decreased as the size of nanoparticles increases. Hwang et al. [26] presented the numerical solution of natural convection in a Al_{2}O_{3}-water mixture in a rectangular cavity heated from below. They used various models to obtain the effective thermal conductivity and viscosity. Their results show that the ratio of heat transfer coefficient of nanofluids to that of base fluid decreased as the average temperature of nanofluids was lowered. Karimipour et al. [27] numerically investigated the mixed convection of a water-copper nanofluid inside a rectangular cavity. They observed that when Reynolds number is less than one, heat transfer rate is much greater than when Reynolds number is more than one. Moreover, they found that increasing the volume fraction of the nanoparticles increases the heat transfer rate.

Oztop and Abu-Nada [28] simulated the natural convection flow in a rectangular cavity by adding a heater placed at the right-hand side of the cavity. Their findings show that the Cu-water mixture has a better heat transfer enhancement compared to Al_{2}O_{3}-water mixture. Jahanshahi et al. [29] presented the numerical simulation of free convection based on experimental measured conductivity in a square cavity using water/SiO2 nanofluid. They investigated the influences of uncertainties due to adopting various formulas for the effective thermal conductivity of silica-water nanofluid on the heat transfer characteristics. They reported that the increase of heat transfer due to experimental formula is more than numerical formula. Also their results have shown that the heat transfer due to numerical formula decreases with increase in volume fraction. Manca et al. [30] numerically investigated on laminar mixed convection in a water-Al_{2}O_{3} nanofluid, flowing in a triangular cross-sectioned duct. They survey the effects of different values of Richardson number and nanoparticle volume fractions on the convective heat transfer of nanofluid. They found that the average convective heat transfer increases by increasing values of Richardson number and nanoparticle volume fraction. Özerinç et al. [31] numerically analyzed the laminar forced convection heat transfer with temperature-dependent thermal conductivity of nanofluids and thermal dispersion inside a straight circular tube. They applied some recent correlations based on a thermal dispersion model for both constant wall heat flux and constant wall temperature boundary conditions. The correlations are according to single-phase approach assumption. The results show that the single-phase assumption is an accurate way of heat transfer enhancement analysis of nanofluids in convective heat transfer. Celli [32] applied a nonhomogenous model for investigating the spatial distribution of the nanoparticles dispersed inside a square cavity subject to different side wall temperatures using nanofluid for natural convection flow. The Brownian motion and the thermophoresis are considered as the leading physical transport mechanisms for the nanoparticles. They reported that for low Rayleigh number nonhomogenous method is appropriate for description the nanofluid systems and for high Rayleigh numbers homogenous method becomes reliable.

The main aim of the present study is the investigation of natural convection heat transfer utilizing nanofluids in vertical cavities. Three different nanoparticles were selected to compare the heat transfer enhancement variations due to the change of nanoparticles. Three different Rayleigh numbers up to the limit of laminar flow regime have been studied in order to elucidate the effect of buoyancy terms. The investigation covered low aspect ratios where conduction heat transfer is marked, up to an aspect ratios close to unity where the convection heat transfer is dominant. A very fine mesh distribution has been used in order to obtain the benchmark solutions for all aspect ratios. The results have been validated with those available in the literature for both pure and nanobased fluids. Finally, the average and maximum Nusselt numbers, streamlines and temperature fields for different values of volume fraction, Rayleigh number, and aspect ratio are illustrated.

#### 2. Problem Statement and Boundary Conditions

Consider a two-dimensional enclosure of height and width with impermeable walls that is filled with nanofluid as shown in Figure 1. The top and the bottom walls are assumed to be insulated, whereas the two vertical walls are maintained constant but with different temperatures. Gravity acts parallel to the active vertical walls pointing toward the bottom wall. The nanofluid is treated as an incompressible and Newtonian fluid. Thermophysical/transport properties of the nanofluid are assumed to be constant, whereas the density variation in the buoyancy force term is handled by the Boussinesq approximation.

The pertinent thermophysical/transport properties are given in Table 1.

#### 3. Governing Equations

Considering the nanofluid as a continuous media with thermal equilibrium between the base fluid and the solid nanoparticles, the governing equations are as follows.

Continuity [29]: X-momentum equation [29]: Y-momentum equation [29]: Energy equation [29]: The density of the nanofluid is given by [29]: whereas the heat capacitance of the nanofluid and part of the Boussinesq term are [29]: with being the volume fraction of the solid particles and subscripts , , and standing for base fluid, nanofluid, and solid, respectively. The viscosity of the nanofluid containing a dilute suspension of small rigid spherical particles is given by [29]: whereas the thermal conductivity of the stagnant nanofluid is [29]: The effective thermal conductivity of the nanofluid is [29]: and the thermal conductivity enhancement term due to thermal dispersion is given by [29]: The empirically determined constant is evaluated following the work of Wakao and Kaguei [34].

The boundary conditions areIn order to write the nondimensional form of the governing equations, we have to introduce a number of references quantities. For example, dimensionless coordinates are defined as where is the dimension of the chosen enclosure (Figure 1). Similarly, fluid properties and the flow variables can be nondimensionalized with respect to the reference quantities. This consideration yields [29] The reference velocity is defined as [29]

Using the above nondimensional variables/parameters, the nondimensional form of the governing equations can be written as
The dimensionless groupings, that is, the Rayleigh and Prandtl numbers and the Boussinesq term, are given by
The variations of the Nusselt numbers along vertical walls and averaged Nusselt number are derived as follows [29]:
Results of a grid independency check for the average Nusselt number are given in Table 2 and Figure 2. Since the difference between the numerical results for grid densities 201 × 201 and 301 × 301 at ^{5} and 10^{6} and grid densities 301 × 301 and 401 × 401 at is very small, 201 × 201 and 301 × 301 were deemed adequate for the case. Basically, these grid densities produce square cells inside the cavity. As for the other aspect ratios, the grid resolution in the -direction was properly reduced to maintain the same square grid resolution inside the domain.

Finally, the following criterion was invoked to assure convergence of the solution:

This definition guarantees a balance between total energy in and out of the cavity.

#### 4. Results and Discussion

The nanofluid in the cavity is chosen as suspension of Al_{2}O_{3}, Cu, and TiO_{2} in the water. The thermophysical/transport properties of the base fluid and nanoparticles are given in Table 1. Many test cases were investigated at different moderate-to-high Rayleigh numbers of 10^{5}, 10^{6}, and 10^{7}, six length-to-height aspect ratios of , 0.2, 0.25, 0.5, 0.75, and 1.0, and three volume fractions , 0.025 and 0.05. The computations are based on an iterative, finitevolume numerical procedure (SIMPLE) using a staggered grid arrangement in combination with the QUICK differencing scheme. The present code for nanofluid has been validated for a water/SiO_{2} nanofluid mixture [29] in a square cavity. Also, the predicted maximum axial and vertical velocities and the average and maximum Nusselt numbers are compared with those of de Vahl Davis [33] in Table 3. Moreover, the present numerical code was also validated against the results of Khanafer et al. [18] for natural convection in an enclosure filled with Cu-water nanofluid at and 10^{5}, and , as shown in Figure 3. Based on these comparisons of the present results to well-established prior work, the accuracy of the computer code was validated.

Figure 4 illustrates comparison of the streamline patterns between the case of Al_{2}O_{3}-water nanofluid at and base water fluid for three Rayleigh numbers. The results show that regardless of the Rayleigh number and the type of fluid, streamlines are characterized by a recirculatory pattern for all the ratios. The heated fluid next to the left wall rises vertically, thus replacing the fluid that travels horizontally toward the cold wall on the right side, and then sinks along the right wall. For a given ratio, as the Rayleigh number is increased, the intensity of convection increases as evidenced by packing of the streamlines. The streamlines at show that the central vortex breaks up into two, three, and four vortexes in the case of a base fluid at , 10^{6}, and 10^{7}, respectively. In the case of a nanofluid at the same streamfunction, the central vortex does not break up at the three Rayleigh numbers. The streamlines also show that the central vortex of a nanofluid occupies a larger zone than that for pure fluid at higher aspect ratios and is smaller at lower aspect ratios. As a result of these discussions, presence of nanoparticles enhances convection heat transfer at higher aspect ratios and enhances conduction heat transfer at lower aspect ratios compared to corresponding cases with a base fluid.

Figure 5 illustrates comparison of isotherm contours between the nanofluid at and base fluid for various Rayleigh numbers. For the highest aspect ratio being the case of a square cavity, the isotherms at the center of the cavity are nearly horizontal and become vertical within the thermal boundary layers next to the vertical active walls. However, the isotherms are vertical within the central part of the thinnest cavity; that is, . This indicates that the role of conduction is dominant in the central region of the thin cavities. However, convection remains important at the top and bottom ends of the thin cavities. The results also show that the isotherms of nanofluid are vertical more than base fluid at whole aspect ratios. This indicates that the effect of conduction in nanofluid is more than base fluid.

Figure 6 presents a comparison of the values of dimensionless determined at against the aspect ratio for the base fluid and nanofluid at for various Rayleigh numbers. For the base fluid, the magnitudes of for are 35.634, 80.832, and 181.52 for , 10^{6}, and 10^{7}, respectively. To make the figure more informative, the magnitude of for each Rayleigh number is nondimensionalized with the magnitude of for the base fluid () for and at the same Rayleigh number; that is, for ). The data presented in Figure 6 show that the magnitude of increases with increasing of the ratio for different Rayleigh numbers and subsequently decreases with increasing of the ratio. As it can be seen, when the Rayleigh number increases, the peaks of distributions shift from middle ratio to lower ratios. The figure indicates that becomes relatively independent of the ratio for and and 10^{7}. As it can be seen, in a nanoparticle suspension both of viscous terms and buoyant terms growth more compare to base fluid. The first one slows down the fluid flow and the second one speeds it up. At lower aspect ratios the first parameter is more dominant and for this we see the maximum velocity decreases when the nanoparticles are used. At higher aspect ratios, the maximum velocity increases because the buoyant terms are dominated.

Similarly, Figure 7 demonstrates the comparison of determined at in terms of the ratio between the base fluid and nanofluid with for three different Rayleigh numbers. To make the figure more informative, the magnitude of at each Rayleigh number is nondimensionalized with the magnitude of for base fluid () at at the same Rayleigh number; that is, at ). The magnitudes of for the case are 73.825, 236.09, and 737.7 for , 10^{6}, and 10^{7}, respectively. The figure shows that the magnitudes of and initially increase and then level off, becoming almost invariant at higher aspect ratios. The present results predict very low magnitude for at lower aspect ratios for and 10^{6}. The results also show that the magnitudes of at each aspect ratio is more than . This may happen by an increase in kinematic viscosity of nanofluid compared to base fluid.

Figure 8 presents the variation of in terms of ratio for different Rayleigh numbers and various volume fractions. The magnitudes of at each Rayleigh number are nondimensionalized with the magnitude of at at the same Rayleigh number; that is, at . The magnitudes of at are 8.495, 19.696, and 45.163 for , 10^{6}, and 10^{7}, respectively. For given geometries, the maximum values of Nusselt number are observed to increase consistently as the volume fraction of the nanoparticles is increased.

The variations of with the ratio at different Rayleigh numbers for the base fluid and nanofluid at are compared in Figure 9. As is observed, the maximum Nusselt number becomes more sensitive to ratio when it decreases for whole range of Rayleigh numbers. The details show that the values of the maximum Nusselt number of nanofluid increase remarkably as compared with the base fluid at three Rayleigh numbers except from to at where values of the maximum Nusselt number of nanofluid decrease as compared with those of base fluid. The details show that the values of maximum Nusselt number at increase after their initial decrease and finally decrease with increasing of the aspect ratio. But the values of maximum Nusselt number at decrease after their initial increase with increasing of the aspect ratio, and finally for , we see a fine decrease in whole aspect ratios.

Figure 10 shows the variation of Nu_{avg} in terms of ratio at different Rayleigh numbers and various volume fractions. The magnitudes of Nu_{avg} at each Rayleigh number are nondimensionalized with the magnitude of at at the same Rayleigh number; that is, at . The magnitudes of at are 4.725, 9.247, and 17.386 for , 10^{6}, and 10^{7}, respectively. As is observed, the values of average Nusselt number increase with increase in volume fraction. When the volume fraction increases, random movement of nanoparticles increases the thermal dispersion in the flow of nanofluid and consequently enhances the heat transfer rates in the enclosure. Also the values of Nu_{avg} decrease with increasing of the aspect ratio. The results show that the maximum and minimum values of Nu_{avg} at each volume fraction for and 10^{6} are seen at and and for at and , respectively. Moreover, the results depict that more increase in Rayleigh number occurred at lower Rayleigh number. It shows that the addition of nanoparticles to the pure fluid has more effects at lower Rayleigh numbers.

Figure 11 displays the Nusselt number distributions on the hot and cold walls for different aspect ratios (, 0.75, 0.5, 0.25, 0.2, and 0.1) and different volume fractions at . As is observed, the figures are entirely symmetric in all six parts of the vertical cavities at each volume fraction. The results show that the value of the Nusselt number increases with increase in volume fraction.

Figure 12 shows the variation of in terms of ratio using three different nanoparticles and Rayleigh numbers at volume fraction equal to 0.05 (). The magnitudes of Nu_{avg} at each Rayleigh number are nondimensionalized with the magnitude of at at the same Rayleigh number; that is, at ). As is observed, the lowest Nu_{avg} was obtained for TiO_{2} due to domination of conduction mode of heat transfer since TiO_{2} has the lowest value of thermal conductivity compared to Cu and Al_{2}O_{3}. However, the thermal conductivity of Al_{2}O_{3} is approximately one-tenth of Cu, as given in Table 1, the values of Nu_{avg} for Al_{2}O_{3} and Cu are close to each other specially at lowest Rayleigh number. However, a unique property of Al_{2}O_{3} is its low thermal diffusivity compared with Cu. The reduced value of thermal diffusivity leads to higher temperature gradients and, therefore, higher enhancements in heat transfer. The Cu nanoparticles have high values of thermal diffusivity and, therefore, this reduces temperature gradients which will affect the performance of Cu nanoparticles.

#### 5. Conclusion

Heat transfer enhancement in a wide range of thin-to-thick vertical cavities subject to different side wall temperatures using nanofluid is studied numerically. The results are presented at different Rayleigh numbers, a wide range of vertical cavity aspect ratios, different volume fractions, and different types of nanoparticles. The present results illustrate that the suspended nanoparticles substantially increase the heat transfer rate at any given Rayleigh number and aspect ratio. In addition, the results illustrate that the average and maximum Nusselt number increase with an increase in volume fraction of nanoparticles. As is observed, the most enhancement of average Nusselt number is seen at and at any volume fraction and also the results show that the average Nusselt number increases with decrease of aspect ratio except from to for . The type of nanofluid is a key factor for heat transfer enhancement. The results illustrate that the highest values of Nusselt number are obtained when using Cu nanoparticles.

#### Nomenclature

AR: | Aspect ratio, defined as |

Specific heat at constant pressure (J/kg K) | |

Nanoparticle diameter (m) | |

EAN: | Enhancement in average Nusselt number |

Gr: | Grashof number, defined as Ra/Pr |

Gravitational acceleration components | |

Cavity height (m) | |

Fluid thermal conductivity (W/m K) | |

Solid thermal conductivity (W/m K) | |

Cavity width (m) | |

Nu: | Nusselt number |

Pr: | Prandtl number, defined as |

Pressure (kg/m s^{2}) | |

Ra: | Rayleigh number, defined as |

Temperature (K) | |

Velocity components (m/s) | |

Nondimensional velocities, defined as | |

VF: | Particle volume fraction |

Cartesian coordinates | |

Dimensionless coordinates. |

*Greek Symbol*

Thermal diffusivity | |

Thermal expansion coefficient | |

Kinematic viscosity | |

Density | |

Dynamic viscosity | |

Particle volume fraction. |

*Subscripts and Superscripts*

avg: | Average magnitude |

: | Cold |

eff: | Effective |

f: | Fluid |

h: | Hot |

max: | Maximum magnitude |

nf: | Nanofluid |

Reference value | |

: | Solid |

*: | Nondimensionless parameter. |

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgment

The authors wish to thank the support and encouragement of their good friend Dr. Seyed Farid Hosseinizadeh, God bless his soul.

#### References

- L. A. Florio and A. Harnoy, “Combination technique for improving natural convection cooling in electronics,”
*International Journal of Thermal Sciences*, vol. 46, no. 1, pp. 76–92, 2007. View at Publisher · View at Google Scholar · View at Scopus - M. Alipanah, P. Hasannasab, S. F. Hosseinizadeh, and M. Darbandi, “Entropy generation for compressible natural convection with high gradient temperature in a square cavity,”
*International Communications in Heat and Mass Transfer*, vol. 37, no. 9, pp. 1388–1395, 2010. View at Publisher · View at Google Scholar - P. O. Iwanik and W. K. S. Chiu, “Temperature distribution of an optical fiber traversing through a chemical vapor deposition reactor,”
*Numerical Heat Transfer A*, vol. 43, no. 3, pp. 221–237, 2003. View at Publisher · View at Google Scholar · View at Scopus - K. O. Lim, K. S. Lee, and T. H. Song, “Primary and secondary instabilities in a glass-melting surface,”
*Numerical Heat Transfer A*, vol. 36, no. 3, pp. 309–325, 1999. View at Scopus - S. Singh and M. A. R. Sharif, “Mixed convective cooling of a rectangular cavity with inlet and exit openings on differentially heated side walls,”
*Numerical Heat Transfer A*, vol. 44, no. 3, pp. 233–253, 2003. View at Scopus - P. M. Patil and P. S. Kulkarni, “Effects of chemical reaction on free convective flow of a polar fluid through a porous medium in the presence of internal heat generation,”
*International Journal of Thermal Sciences*, vol. 47, no. 8, pp. 1043–1054, 2008. View at Publisher · View at Google Scholar - M. Lamsaadi, M. Naimi, M. Hasnaoui, and M. Mamou, “Natural convection in a vertical rectangular cavity filled with a non-newtonian power law fluid and subjected to a horizontal temperature gradient,”
*Numerical Heat Transfer A*, vol. 49, no. 10, pp. 969–990, 2006. View at Publisher · View at Google Scholar · View at Scopus - L. K. Saha, M. A. Hossain, and R. S. R. Gorla, “Effect of Hall current on the MHD laminar natural convection flow from a vertical permeable flat plate with uniform surface temperature,”
*International Journal of Thermal Sciences*, vol. 46, no. 8, pp. 790–801, 2007. View at Publisher · View at Google Scholar · View at Scopus - Y. Y. Yan, H. B. Zhang, and J. B. Hull, “Numerical modeling of electrohydrodynamic (EHD) effect on natural convection in an enclosure,”
*Numerical Heat Transfer A*, vol. 46, no. 5, pp. 453–471, 2004. View at Publisher · View at Google Scholar · View at Scopus - J. C. A. Maxwell,
*Treatise on Electricity and Magnetism*, vol. 1, Clarendon Press, Oxford, UK, 2nd edition, 1881. - R. L. Hamilton and O. K. Crosser, “Thermal conductivity of heterogeneous two-component systems,”
*Industrial and Engineering Chemistry Fundamentals*, vol. 1, no. 3, pp. 187–191, 1962. View at Scopus - A. Einstein, “Eine neue Bestimmung der Molekuldimensionen,”
*Annalen der Physik*, vol. 19, pp. 289–306, 1906. - H. C. Brinkman, “The viscosity of concentrated suspensions and solutions,”
*The Journal of Chemical Physics*, vol. 20, no. 4, p. 571, 1952. View at Scopus - Y. Xuan and Q. Li, “Heat transfer enhancement of nanofluids,”
*International Journal of Heat and Fluid Flow*, vol. 21, no. 1, pp. 58–64, 2000. View at Publisher · View at Google Scholar · View at Scopus - I. C. Bang and S. H. Chang, “Boiling heat transfer performance and phenomena of Al
_{2}O_{3}-water nano-fluids from a plain surface in a pool,”*International Journal of Heat and Mass Transfer*, vol. 48, no. 12, pp. 2407–2419, 2005. View at Publisher · View at Google Scholar · View at Scopus - H. Masuda, A. Ebata, K. Teramae, and N. Hishinuma, “Alteration of thermal conductivity and viscosity of liquid by dispersing ultra-fine particles,”
*Netsu Bussei*, vol. 7, pp. 227–233, 1993. - U. S. Choi, “Enhancing thermal conductivity of fluids with nanoparticles,” in
*Developments and Application of Non-Newtonian Flows*, pp. 99–105, ASME, 1995. - K. Khanafer, K. Vafai, and M. Lightstone, “Buoyancy-driven heat transfer enhancement in a two-dimensional enclosure utilizing nanofluids,”
*International Journal of Heat and Mass Transfer*, vol. 46, no. 19, pp. 3639–3653, 2003. View at Publisher · View at Google Scholar · View at Scopus - C. J. Ho, M. W. Chen, and Z. W. Li, “Numerical simulation of natural convection of nanofluid in a square enclosure: effects due to uncertainties of viscosity and thermal conductivity,”
*International Journal of Heat and Mass Transfer*, vol. 51, no. 17-18, pp. 4506–4516, 2008. View at Publisher · View at Google Scholar · View at Scopus - A. G. A. Nnanna, “Experimental model of temperature-driven nanofluid,”
*Journal of Heat Transfer*, vol. 129, no. 6, pp. 697–704, 2007. View at Publisher · View at Google Scholar · View at Scopus - C. J. Ho and C. C. Lin, “Experiments on natural convection heat transfer of a nanofluid in a square enclosure,” in
*Proceedings of the 13th International Heat Transfer Conference*, Sydney, Australia, 2006. - E. Abu-Nada, Z. Masoud, and A. Hijazi, “Natural convection heat transfer enhancement in horizontal concentric annuli using nanofluids,”
*International Communications in Heat and Mass Transfer*, vol. 35, no. 5, pp. 657–665, 2008. View at Publisher · View at Google Scholar · View at Scopus - X. Q. Wang and A. S. Mujumdar, “Heat transfer characteristics of nanofluids: a review,”
*International Journal of Thermal Sciences*, vol. 46, no. 1, pp. 1–19, 2007. View at Publisher · View at Google Scholar · View at Scopus - J. M. Khodadadi and S. F. Hosseinizadeh, “Nanoparticle-enhanced phase change materials (NEPCM) with great potential for improved thermal energy storage,”
*International Communications in Heat and Mass Transfer*, vol. 34, no. 5, pp. 534–543, 2007. View at Publisher · View at Google Scholar · View at Scopus - X. Q. Wang, A. S. Mujumdar, and C. Yap, “Free convection heat transfer in horizontal and vertical rectangular cavities filled with nanofluids,” in
*Proceedings of the International Heat Transfer Conference*, Sydney, Australia, 2006. - K. S. Hwang, J. H. Lee, and S. P. Jang, “Buoyancy-driven heat transfer of water-based Al
_{2}O_{3}nanofluids in a rectangular cavity,”*International Journal of Heat and Mass Transfer*, vol. 50, no. 19-20, pp. 4003–4010, 2007. View at Publisher · View at Google Scholar · View at Scopus - A. Karimipour, A. H. Nezhad, A. Behzadmehr, S. Alikhani, and E. Abedini, “Periodic mixed convection of a nanofluid in a cavity with top lid sinusoidal motion,”
*Proceedings of the Institution of Mechanical Engineers C*, vol. 225, no. 9, pp. 2149–2160, 2011. View at Publisher · View at Google Scholar · View at Scopus - H. F. Oztop and E. Abu-Nada, “Numerical study of natural convection in partially heated rectangular enclosures filled with nanofluids,”
*International Journal of Heat and Fluid Flow*, vol. 29, no. 5, pp. 1326–1336, 2008. View at Publisher · View at Google Scholar · View at Scopus - M. Jahanshahi, S. F. Hosseinizadeh, M. Alipanah, A. Dehghani, and G. R. Vakilinejad, “Numerical simulation of free convection based on experimental measured conductivity in a square cavity using water/SiO
_{2}nanofluid,”*International Communications in Heat and Mass Transfer*, vol. 37, no. 6, pp. 687–694, 2010. View at Publisher · View at Google Scholar · View at Scopus - O. Manca, S. Nardini, D. Ricci, and S. Tamburrino, “Numerical investigation on mixed convection in triangular cross-section ducts with nanofluids,”
*Advances in Mechanical Engineering*, vol. 2012, Article ID 139370, 13 pages, 2012. View at Publisher · View at Google Scholar · View at Scopus - Özerinç, A. G. Yazicioğlu, and S. Kakaç, “Numerical analysis of laminar forced convection with temperature-dependent thermal conductivity of nanofluids and thermal dispersion,”
*International Journal of Thermal Sciences*, vol. 62, pp. 138–148, 2012. View at Publisher · View at Google Scholar - M. Celli, “Non-homogeneous model for a side heated square cavity filled with a nanofluid,”
*International Journal of Heat and Fluid Flow*, vol. 44, pp. 327–335, 2013. View at Publisher · View at Google Scholar - G. de Vahl Davis, “Natural convection of air in square cavity: a benchmark solution,”
*International Journal for Numerical Methods in Fluids*, vol. 3, no. 3, pp. 249–264, 1983. View at Scopus - N. Wakao and S. Kaguei,
*Heat and Mass Transfer in Packed Beds*, Gordon and Breach Science Publishers, New York, NY, USA, 1982.