Journal of Applied MathematicsVolume 2012 (2012), Article ID 259284, 18 pagesdoi:10.1155/2012/259284

Research Article

Numerical Investigation of Nanofluid Forced Convection in Channels with Discrete Heat Sources

Payam Rahim Mashaei, Seyed Mostafa Hosseinalipour, and Mehdi Bahiraei

CAE Lab & CFD Center, Mechanical Engineering Department, Faculty of Engineering, Iran University of Science & Technology (IUST), Tehran, Iran

Received 6 January 2012; Revised 21 January 2012; Accepted 22 January 2012

Academic Editor: M. F. El-Amin

Copyright © 2012 Payam Rahim Mashaei 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.

Numerical simulation is performed to investigate the laminar force convection of Al₂O₃/water nanofluid in a flow channel with discrete heat sources. The heat sources are placed on the bottom wall of channel which produce much thermal energy that must be evacuated from the system. The remaining surfaces of channel are kept adiabatic to exchange energy between nanofluid and heat sources. In the present study the effects of Reynolds number (, and 1000), particle volume fraction ( (distilled water), 1 and 4%) on the average heat transfer coefficient (h), pressure drop (), and wall temperature () are evaluated. The use of nanofluid can produce an asymmetric velocity along the height of the channel. The results show a maximum value 38% increase in average heat transfer coefficient and 68% increase in pressure drop for all the considered cases when compared to basefluid (i.e., water). It is also observed that the wall temperature decreases remarkably as Re and ϕ increase. Finally, thermal-hydraulic performance (η) is evaluated and it is seen that best performance can be obtained for and %.

1. Introduction

Localized areas of high temperature on microprocessors and various electronic components produce hot spots that have an unfavorable effect on their performance and operating conditions. With increasing of power density of these electronic components, good attempts have been carried out to enhance the heat exchanger rate of them by active as well as passive methods. While the former usually offers higher augmentation, it requires additional external forces that can increase the capital and operating cost of the system. In contrast, passive heat transfer enhancement can be obtained by changing the geometry or modifying thermal properties of working fluid [1]. In recent years, the advances in manufacturing technology have caused the production of nanoparticles and created a special class of fluids, called “nanofluids”. The term “nanofluid” refers to a two-phase mixture where the continuous phase is usually a liquid and the dispersed phase is created of extremely fine metallic particles of size about 50 nm, called “nanoparticles” [2].

Many researchers experimentally showed nanofluids have higher thermal conductivity than those of the base fluids and a lot of correlations were reported. For example Lee et al. [3] evaluated Al₂O₃-water/ethylene glycol with particle diameters 24.4 and 38.4 nm as well as CuO-water/ethylene glycol with particle diameters 18.6 and 23.6 nm and showed that thermal conductivity increases to 20% as particle volume fraction increases from 0 to 4%. Chopkar el al. [4] measured thermal conductivity of Ag₂Al-water nanofluids and Al₂Cu-water nanofluids and found that it increases about 130% with a volume fraction less than 1%. Some researchers [5–7] investigated the thermal conductivity of nanofluids as a function of volume fraction and temperature. They showed that the thermal conductivity of nanofluids remarkably increases as temperature rises. Wang and Mujumdar [8] listed a large number of correlations, which are available in the literatures, for thermal conductivity of nanofluids.

Santra et al. [11] numerically studied the effect of CuO-water nanofluid as a cooling medium to simulate the heat transfer behavior in a two-dimensional (infinite depth) horizontal rectangular duct, where top and bottom walls were two isothermal symmetric heat sources. They considered the fluid Newtonian as well as non-Newtonian for a wide range of Reynolds numbers and solid volume fractions. These authors observed that the heat transfer enhancement is possible using nanofluid in comparison with conventional fluid for the both cases. Maiga et al. [12] numerically investigated the flow and thermal behavior of water-Al₂O₃ and Ethylene Glycol-Al₂O₃ mixtures for a system of parallel, coaxial, and heated disks. A remarkable augmentation of heat transfer coefficient has been observed with increasing of the volume fraction of nanoparticles for both nanofluids. They have reported that the rate of increase of heat transfer is more for Ethylene Glycol-Al₂O₃ nanofluid in comparison with the water-Al₂O₃ nanofluid. However, the wall shear stress also increases considerably with increasing of volume fraction of nanoparticles. Feng and Kleinstreuer [13] executed a numerical simulation for alumina-water nanofluid flow with heat transfer between parallel disks. They reported that nanofluid produce smoother flow fields and temperature distributions and heat transfer rate increases with higher volume fraction, smaller nanoparticle diameter, reduced disk spacing, and larger inlet Reynolds number.

With respect the problem under study, that is, heat transfer of discrete heat sources in channel flows, there are numerous works [14–17], both numerical and experimental, which consider such a problem.

Bhowmik et al. [14] performed steady-state experiments to investigate general heat transfer from an in-line four simulated electronic chips in a vertical rectangular channel employing water as the fluid working. They evaluated the effects of heat fluxes, flow rate, and geometrical parameters on heat transfer coefficient and illustrated heat transfer rate strongly depends on flow rate. da Silva et al. [15] employed two different analytical approaches to determine how to arrange discrete heat sources on wall cooled by force convection: (i) large number of small heat sources and (ii) small number of heat sources with finite length, which are mounted on a flat wall. Both analyses proved that heat sources should be placed nonuniformly on the wall. Arquis et al. [16] numerically examined the fluid flow and heat transfer characteristics associated with cooling an in-line array of discrete heated blocks in a channel by using a single laminar slot air jet. They studied the effects of various values such as channel height, slot weight, jet Reynolds number, spacing between blocks, block height, and block thermal conductivities. In general, the heat transfer rate increases with the increase of Reynolds number and the decrease of channel weight. The effective cooling of blocks is observed to increase for shorter and widely spaced heated blocks. Jassim and Muzychka [17] studied the problem how to allocate discrete heat source to the space on a wall of a convergent flow. They showed that the heat sources should be distributed nonuniformly. Furthermore, the optimal spacing between heat sources depends on both Reynolds number and channel shape factor.

In the present paper, the flow and heat transfer characteristics of channel flow with discrete heat sources for base fluid (distilled water) and a nanofluid that is composed of distilled water and Al₂O₃ nanoparticles are numerically investigated. The main aim of this study is how the nanofluid affects on the heat transfer rate and pressure drop of flow in a channel with hot spots.

2. Mathematical Formulation

2.1. Geometry Configuration

In this study the velocity and temperature fields are determined in a parallel plates channel with height and Length , as shown in Figure 1. The length of duct is 50 times more than height. Nineteen heat sources of of length are placed on the bottom wall of the channel. The remaining bottom wall and whole top wall are isolated. The spacing between heat sources is .

Figure 1: Parallel plates channel with discrete heat sources.

2.2. Governing Equations and Boundary Conditions

Since nanofluids are composed of extremely small particles, it is assumed that the nanoparticles and basefluid are in thermal equilibrium and they flow at same velocity. In the present work, the nanofluid is considered incompressible with temperature-dependent properties. The compression work and viscous dissipation terms were considered negligible in the energy equation. Under such assumptions, the general governing equations written are as the followings.

Conservation of mass:

Conservation of momentum:

Conservation of energy: Equations (2.1)–(2.3) must be solve by using appropriate boundary conditions. The flow boundary conditions are a uniform velocity at the channel inlet, equal to , no-slip at the channel walls and zero relative pressure at channel outlet. The thermal boundary conditions are a uniform temperature at the channel entrance, equal to and an adiabatic condition at channel walls, except along the heat sources, in which a uniform heat flux is applied.

2.3. Nanofluid Thermophysical Properties

The thermophysical properties of nanofluid are chiefly functions of particle volumetric concentration and temperature. In the absence of experimental data, nanofluid density and specific heat are defined only as a function of volume fraction as follow.

Density:

Specific heat: In above equations, subscripts “bf”, “” and “nf” refer to basefluid, nanoparticle and the nanofluid, respectively.

Viscosity:

Thermal conductivity: Equations (2.6) to (2.8) were presented by Roy et al. [18] for water-Al₂O₃ nanofluid based on available experimental results published by Putra et al. [19].

3. Numerical Method

The governing differential equations are solved using the control volume method. A second order upwind method is used for energy and momentum equations. The SIMPLE procedure is chosen to couple pressure and velocity. The solution converge was met when the normalized residuals for all equations reached to the 10⁻⁷. The algebraic discretized equations throughout the physical domain are solved by means of well-known TDMA techniques.

In order to assess the grid independent of numerical solution, three grid densities are checked. Figure 2 illustrates the effect of grid size on the predicted temperature at outlet channel for and . As it can be seen a grid density of provides satisfactory solution for the shown example.

Figure 2: Effect of grid density on the temperature of channel outlet (, , ).

4. Results and Discussion

The thermal performance of the channel is characterized in terms of average heat transfer coefficient along the heat sources, , defined as: where is the coordinate along the heat sources and is given by where and are heat flux of heat sources and wall temperature, respectively.

Thermalhydraulic performance factor is defined as: where and are, respectively, average heat transfer coefficient and pressure drop ratio, referred to values obtained for basefluid.

The heat transfer performance of cooling channel is discussed in term of the figure of merit, FoM, which is given by [1] Here is the pumping power, which is given as: where is pump efficiency which assumed to be 70%. The total heat transfer rate is calculated as: where, is mass flow rate and and are mean temperature at inlet and outlet, respectively. Mean temperature is obtained by

4.1. Code Validation

In order to show the validity and also accuracy of the model and numerical method, two comparisons with the available data are carried out. The first comparison is related to a parallel plates channel that all its walls are heated with a constant heat flux and the water is used as fluid working. In this case, the Nusselt number is compared which is given by the following definitions: where is bulk temperature of fluid. The comparison is depicted in Figure 3. It is seen that the results of the present study are in good agreement with previous data [9].

Figure 3: Comparison of present numerical results with those obtained by Shah and London [9].

The second comparison is concerned with experimental data of nanofluid heat transfer in a circular tube with diameter () equal to 4.75 mm and Length of 1.2 m, given by Anoop et al. [10]. The nanofluid is water-Al₂O₃ with at a constant Reynolds number 1588. Figure 4 shows remarkable agreement between the present results on local heat transfer coefficient and those provided by [10].

Figure 4: Comparison of heat transfer coefficient between present simulation and experimental data [10] for water-Al₂O₃ nanofluid with in the circular tube.

4.2. Velocity Field

The values of velocity at channel outlet for and are depicted in Figures 5(a) and 5(b), respectively, for considered volume fractions. It can be seen that trend of velocity profile is affected by adding the nanoparticle in basefluid for both Reynolds numbers. However, with the increase of nanoparticle loading, as it is expected, the peak of velocity profile will increase. This produces larger velocity gradient near the wall which can lead to increase wall shear stress and pressure drop. Another interesting issue, which one can perceive in Figure 5, is the asymmetric profile for and because of temperature-depended viscosity. The magnitude of velocity near the bottom wall is lower due to decreased viscosity. This is because of higher temperature near the bottom wall that decreases the viscosity according to (2.6). The velocity profile is more asymmetric for and . The values of velocity at and (i.e., near the bottom and top walls) are, respectively, 0.00083 m/s and 0.00069 m/s for and which indicate about 20% difference, while this value is equal to 0 (symmetric) and 16% for (distilled water) and , respectively.

(a)

(b)

Figure 5: Effect of particle loading parameter on velocity profile at channel outlet for: (a) , (b) .

4.3. Local Nanofluid Properties

Figures 6(a) and 6(b) show temperature effects on local nanofluid conductivity which is computed inside the numerical model with (2.7) and (2.8). Since the whole top wall is adiabatic and there is no heat transfer, the nanofluid conductivity is only illustrated on bottom wall. One can easily notice that the behavior of nanofluid conductivity is periodic along the channel. This is due the fact that heat sources increase wall temperature and nanofluid conductivity. Furthermore, the remaining bottom wall, which is adiabatic, has lower temperature that leads to lower nanofluid conductivity. With comparison of Figures 6(a) and 6(b), It can be seen that presents higher conductivity when compared to . The reason is explained as flows. Although the higher Reynolds number increases the heat transfer rate, but, it is not allowed fluid particles to residence for a long time in channel because of its higher mass flow rate. So, for the nanofluid temperature rises more in comparison with and presents higher conductivity.

(a)

(b)

Figure 6: Effect of particle loading parameter on nanofluid conductivity profile along the channel for: (a) , (b) .

Since the variations of temperature on the top wall is negligible (less than 2 K), the local viscosity is illustrated in Figures 7(a) and 7(b) for and , respectively, only for bottom wall. It can be observed that the local viscosity has a decreasing periodic behavior on bottom wall because of the growing temperature along the channel. As shown in Figure 7, similar to the conductivity profile, the amplitude of variations in a period of local viscosity is higher for in comparison with . This is due this fact that the lower Reynolds number increases the residence time distribution (RTD) of nanofluid in a period of channel and causes more temperature difference.

(a)

(b)

Figure 7: Effect of particle loading parameter on nanofluid viscosity along the channel for: (a) , (b) .

4.4. Effect of Particle Volume Fraction on Heat Transfer and Temperature Field

Results reveal that the presence of nanoparticles has a remarkable effect on heat transfer enhancement. The average heat transfer coefficient profiles on the heat sources as a function of Reynolds number are depicted in Figure 8(a) for , and 4. In general, the average heat transfer coefficient increases as Reynolds number is increased. It is also observed that as increases becomes higher for a fixed value of Re. Thus, passing from to %, the maximum value of about 2366 W/m² k for is found at and .

(a)

(b)

Figure 8: Effect of parameters and Re on: (a) average heat transfer coefficient and (b) average heat transfer coefficient ratio.

Figure 8(b) illustrates the average heat transfer coefficient ratio (), referred to the values obtained for basefluid, as a function of Reynolds number for and 4%. It can be seen that the is greater than one for all considered cases and rises slightly as Reynolds number increases from 200 to 1000 for both volume fractions and 4%. Furthermore, there is a considerable augmentation for higher volume fraction, in fact the highest value of 1.38 is detected at and .

As we are also interested to analysis the cooling benefits of nanofluid on bottom wall, Figures 9(a) and 9(b) show the wall temperature for various volume fractions for and 1000, respectively. Due to a nonuniform heat flux on bottom wall, an increasing periodic behavior is detected in Figure 9. One can see that the wall temperature is more for in comparison with . As mentioned earlier, this is due this fact that the fluid particles spend more time in channel. The peak values of bottom wall temperature are related to the heat sources, while the minimum values belong to adiabatic part of bottom wall. It is also observed that among the considered particle volume fractions, the lowest minimum and maximum values of wall temperature in each period belong to . It can also be noted the effect of volume fraction on wall cooling becomes more within the farther distance from inlet, that is, greater .

(a)

(b)

Figure 9: Effect of particle loading parameter on bottom wall temperature at: (a) , (b) .

The effect of particle volume fraction on temperature along the channel height can be seen on Figures 10(a) and 10(b), depicting the temperature profile at the first of the last heat source for and 1000, respectively. It is observed that the values of temperature are higher for . On the bottom wall, that is, , passing from to the temperature gradient increases which presents a better heat transfer enhancement. By comparison of Figures 10(a) and 10(b), it is noted that the effect of the use of nanofluid on the temperature along the height of channel is greater for lower Reynolds.

(a)

(b)

Figure 10: Effect of volume fraction on temperature profile along channel height at the first of last heat source for: (a) , (b) .

4.5. Effect of Particle Volume Fraction on Pressure Drop

It will be predictable that the use of nanoparticles in basefluid can have an adverse influence on pressure drop because of increased viscosity. The pressure drop profiles as a function of Reynolds number are depicted in Figure 11(a) for , and 4%. As it can be noted, the pressure drop increases as Reynolds number and particle volume fraction increase and the highest value is found at and . By comparison , and , it is observed that the difference between the values of pressure drop is less at lower Reynolds number, while it will be more considerable at higher Reynolds number. This can be explained by the temperature-dependent viscosity of nanofluid.

(a)

(b)

Figure 11: Effect of parameters and Re on: (a) pressure drop and (b) pressure drop ratio.

The pressure drop ratio (), referred to the base fluid is described in Figure 11(b). It is observed that the is greater than one for all considered cases, similar to , and rises slightly as Reynolds number increases from 200 to 1000 for both volume fractions and 4%. In fact the highest value of pressure drop ratio is equal to 1.68 that is dedicated at and .

4.6. Effect of Heat Flux of Heat Sources on Heat Transfer and Pressure Drop

Since we use temperature-dependent properties, it is interesting to determine the effect of heat sources heat flux () on heat transfer coefficient and pressure drop. The effects of () on the average heat transfer and pressure drop can be seen in Figures 12(a) and 12(b), depicting, respectively, the heat transfer coefficient and pressure drop for and 4%. It can easily be seen in Figure 12(a) that the average heat transfer coefficient increases as () is increased because of increased nanofluid conductivity. On the other hand, as shown in Figure 12(b), the pressure drop decreases as () increases because of decreased nanofluid viscosity. Thus, the augmentation of () has a good effect on thermal-hydraulic performance of channel.

(a)

(b)

Figure 12: Effect of parameters and on (a) averaged heat transfer coefficient and (b) pressure drop. ().

4.7. Thermal Hydraulic and Overall Heat Transfer Performance

As discussed in previous sections, the use of nanofluids increases the heat transfer rate as well as pressure drop. In order to investigate the order of magnitude of augmentation of heat transfer and pressure drop for various Reynolds number and volume fraction, the thermal hydraulic performance factor as a function of Reynolds number is depicted in Figure 13 for and 4%.

Figure 13: Effect of parameters and Re on .

It is observed that the thermal hydraulic performance has different behaviors for and 4%. For , thermal hydraulic performance decreases as Reynolds number increases from 50 to 200 and after it, the profile will be unchangeable. For , thermal hydraulic performance decreases as Reynolds number increases from 50 to 200, similar to %, but there is a remarkable augmentation as Reynolds number increases from 200 to 1000 and the best performance of 1.16 is found at and %.

In order to have a comparison between ratio of output (heat transfer rate) to the input (pumping power), the overal heat transfer performance of channel in term of figure of merit (FoM) as a function of Reynolds number is illustrated in Figure 14 for and 4%. It can be seen that FoM decreases as volume fraction increases. This is due this fact that specific heat decreases as volume fraction increases. Furthermore, the increase of volume fraction leads to increased inlet velocity, which decreases residence time distribution of nanofluid. So, although the heat transfer coefficient increases, but overall heat transfer performance decreases as volume fraction increases. As Figure 14 indicates, the FoM remarkably decreases as Reynolds number increases and the best value of 207347 is found at and %.

Figure 14: Effect of parameters and Re on FoM.

5. Conclusion

A numerical analysis of flow and heat transfer characteristics of nanofluid in a parallel plates channel with discrete heat sources has been presented. The heat sources are placed on bottom wall at a constant heat flux and remaining channel surfaces are considered adiabatic. The basefluid is water and three volume fractions of Al₂O₃ nanoparticles ( (distilled water) 1 and 4%) are taken into account with a single-phase model. Furthermore, different Reynolds numbers in the range 50–1000 are considered. The use of nanofluid can cause an asymmetric profile velocity along the height of channel. Results clearly show that the use of nanofluid can remarkably increase heat transfer rate and decrease the wall temperature. Furthermore, for increasing volume fraction, nanofluid can produce lower wall temperature and higher heat transfer rate. The average heat transfer coefficient and pressure profiles present that heat transfer rate increases as applied heat flux of heat sources is increased while the pressure drop decreases.

Nomenclature

Symbols

: | nabla operatoe (1/m)
: | specific heat of the fluid (J/kg K)
FoM: | figure of merit
: | channel height (m)
: | average heat transfer coefficient along heat sources (W/m² K)
: | thermal conductivity (W/m K)
: | mass flow rate (Kg/s)
Nu: | local Nusselt number ()
: | pressure (Pa)
: | heat flux of heat sources (W/m²)
Re: | reynolds number ()
: | coordinate along heat sources (m)
: | temperature (K)
: | velocity vector (m/s)
: | total heat transfer rate (J/s)
: | inlet velocity (m/s).

Greeks

: | thermal hydraulic performance factor
: | pump efficiency
: | dynamic viscosity (Pa·s)
: | density (kg/m³)
: | particle volume fraction (%).

Subscripts

: | refers to the wall conditions
bf: | refers to the base fluid
nf: | refers to the nanofluid
: | refers to particles
: | refers to a ratio
0: | refers to the inlet conditions.

Acknowledgments

First, the authors would like to thank Editor Dr. M. F. El-Amin and two anonymous referees for their useful suggestions, which greatly improve the paper. We also would like to express our appreciation to the Iran nanotechnology initiative council for their financial support for this project.

References

A.P. Saamito, J.C. Kurnia, and A. S. Mujumdar, “Numerical evaluation of laminar heat transfer enhancement in nanofluid flow in coiled square tube,” Nanoscale Research Letter, vol. 6, p. 376, 2011.

G. Roy, S. J. Palm, and C. T. Nguyen, “Heat transfer and fluid flow of nanofluids in laminar radial flow cooling systems,” Journal of Thermal Science, vol. 14, no. 4, pp. 362–367, 2005.

S. Lee, S. U. S. Choi, S. Li, and J. A. Eastman, “Measuring thermal conductivity of fluids containing oxide nanoparticles,” Journal of Heat Transfer, vol. 121, no. 2, pp. 280–288, 1999.

I. Chopkar, S. Sudarshan, P. K. Das, and I. Manna, “Effect of particle size on thermal conductivity of nanofluid,” Metallurgical and Materials Transactions A, vol. 39, no. 7, pp. 1535–1542, 2008.

C. H. Li and G. P. Peterson, “Experimental investigation of temperature and volume fraction variations on the effective thermal conductivity of nanoparticle suspensions (nanofluids),” Journal of Applied Physics, vol. 99, no. 8, Article ID 084314, 2006.

S. K. Das, N. Putra, P. Thiesen, and W. Roetzel, “Temperature dependence of thermal conductivity enhancement for nanofluids,” Journal of Heat Transfer, vol. 125, no. 4, pp. 567–574, 2003.

T. P. Teng, Y. H. Hung, T. C. Teng, and J. H. Chen, “Performance evaluation on an air-cooled heat exchanger for alumina nanofluid under laminar flow,” Nanoscale Research Letter, vol. 6, p. 488, 2011.

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.

R. K. Shah and A. L. London, Laminar Flow Forced Convection in Ducts, Academic Press, New York, NY, USA, 1978.

K. B. Anoop, T. Sundararajan, and S. K. Das, “Effect of particle size on the convective heat transfer in nanofluid in the developing region,” International Journal of Heat and Mass Transfer, vol. 52, no. 9-10, pp. 2189–2195, 2009.

A. K. Santra, S. Sen, and N. Chakraborty, “Study of heat transfer due to laminar flow of copper-water nanofluid through two isothermally heated parallel plates,” International Journal of Thermal Sciences, vol. 48, no. 2, pp. 391–400, 2009.

S. E. B. Maiga, S. J. Palm, C. T. Nguyen, G. Roy, and N. Galanis, “Heat transfer enhancement by using nanofluids in forced convection flows,” International Journal of Heat and Fluid Flow, vol. 26, no. 4, pp. 530–546, 2005.

Y. Feng and C. Kleinstreuer, “Nanofluid convective heat transfer in a parallel-disk system,” International Journal of Heat and Mass Transfer, vol. 53, no. 21-22, pp. 4619–4628, 2010.

H. Bhowmik, C. P. Tso, K. W. Tou, and F. L. Tan, “Convection heat transfer from discrete heat sources in a liquid cooled rectangular channel,” Applied Thermal Engineering, vol. 25, no. 16, pp. 2532–2542, 2005.

A. K. da Silva, S. Lorente, and A. Bejan, “Optimal distribution of discrete heat sources on a plate with laminar forced convection,” International Journal of Heat and Mass Transfer, vol. 47, no. 10-11, pp. 2139–2148, 2004.

E. Arquis, M. A. Rady, and S. A. Nada, “A numerical investigation and parametric study of cooling an array of multiple protruding heat sources by a laminar slot air jet,” International Journal of Heat and Fluid Flow, vol. 28, no. 4, pp. 787–805, 2007.

E. Jassim and Y. S. Muzychka, “Optimal distribution of heat sources in convergent channels cooled by laminar forced convection,” Journal of Heat Transfer, vol. 132, no. 1, pp. 1–8, 2010.

G. Roy, S. J. Palm, and C. T. Nguyen, “Heat transfer and fluid flow of nanofluids in laminar radial flow cooling systems,” Journal of Thermal Science, vol. 14, no. 4, pp. 362–367, 2005.

N. Putra, W. Roetzel, and S. K. Das, “Natural convection of nano-fluids,” Heat and Mass Transfe, vol. 39, no. 8-9, pp. 775–784, 2003.