Document Type : Full Lenght Research Article
Authors
^{1} Faculty of Mechanical Engineering, Semnan University, Semnan, Iran
^{2} Materials and Energy Research Center (MERC), Karaj, Iran
Abstract
Keywords
Journal of Heat and Mass Transfer Research
Journal homepage: http://jhmtr.journals.semnan.ac.ir 

A B S T R A C T
Rheological characteristics of Al_{2}O_{3}, CuO and TiO_{2} nano particles were investigated in oil as the base fluid at 1 and 2 wt.%. Constitutive relations for nonNewtonian fluid were discussed based on the powerlaw model. Measured viscosities of each nanofluid were used to evaluate the powerlaw and consistency index. Results indicated that the nanofluid viscosity decreased by increasing the concentration. Oil showed shear thickening behavior while nanofluids showed shear thinning behavior. An increase in nanoparticle concentration caused a decrease in the powerlaw index beside an increase in the consistency index. Moreover, the present study showed that the effective viscosity of fluids would be decreased by nanoparticle addition at some wt.% and some shear rates. Furthermore, results showed that the classic models for nanofluid viscosity couldn’t predict their real values of nano fluid viscosity, as the measured values are less than the predicted ones.
© 2014 Published by Semnan University Press. All rights reserved. 
An experimental investigation of rheological characteristics of nonNewtonian nanofluids
Milad Tajik JamalAbad^{1}, Maziar Dehghan^{1}, Seyfolah Saedodin^{1}, Mohammad Sadegh Valipour^{1}, Amirhossein Zamzamian*^{2}
^{1} Faculty of Mechanical Engineering, Semnan University, Semnan, Iran ^{2}Materials and Energy Research Center (MERC), Karaj, Iran

History: Received 4 August 2013 Received in revised form 15 October 2013 Accepted 16 November 2013
Keywords: Nanofluids Viscosity Rheological characteristics Powerlaw index Consistency index


1. Introduction
Numerous researchers have been investigating better ways to enhance the thermal performance of heat transfer fluids. One of the methods used is to add nanosized particles of high thermal conductive materials like carbon, metals, and metal oxides into the heat transfer fluid to improve the overall thermal conductivity of the fluid.

Viscosity of a heat transfer fluid is important with respect to studying its convective heat transfer and the pumping power required for practical applications. In industry, an optimization is required between heat transfer capabilities and the frictional characteristics arising from the viscosity of the fluid. Experimental data for the effective viscosity of aqueous nanofluids is limited to certain nanoparticles, such as Al_{2}O_{3} [14], CuO [4, 5], TiO_{2}[1] and MWCNT [6]. Most of these works have been directed towards metal oxide nanoparticles. The parameters against which viscosity was studied were particle volumetric concentration, temperature, and shear rate. Most of the works have mainly focused on spherical nanoparticles of metal oxides, and on the basis of Einstein theory [7]. Viscoelasticfluidbased nanofluids with dispersion of copper (Cu) nanoparticles in viscoelastic surfactant solution were investigated by Yang et al. [8]. They found that the viscoelasticfluidbased Cu nanofluid shows a nonNewtonian behavior in its viscosity, and the viscosity increases with a reduction in the temperature. Ravi Prasher et al. [9] performed some experiments and gave results on the viscosity of aluminabased nanofluids reported for various shear rates, temperature, nanoparticle diameter, and nanoparticle volume fraction. From their data it seems that the increase in the nanofluid viscosity is considerable in comparison with the enhancement in the thermal conductivity reported in the literature. Lee et al. [10] studied the viscosity of SiCwater nanofluids and showed that the viscosity of nanofluids increases along with increased concentration and decreased temperature. Utomo et al. [11] investigated thermal conductivity, viscosity, and heat transfer coefficient of titania and alumina nanofluids. Their results illustrated that the viscosity of alumina and titania nanofluids was higher than the prediction of Einstein–Batchelor model. The rheological properties of titania nanofluids with different base fluids were studied by Chen et al. [12, 13]. They found that the titaniaethylene glycol nanofluids showed a Newtonian behavior, while the titaniawater showed a nonNewtonian behavior. These finding show that the base fluid has a strong effect on the rheological behavior. Experimental results on the viscosity of the nanofluid prepared by dispersing alumina nanoparticles (<50 nm) in commercial car coolant were reported by Kole and Dey [14]. Their results showed that the pure base fluid displayed Newtonian behavior over the measured temperature, but it transformed to a nonNewtonian fluid with addition of a small amount of alumina nanoparticles. Viscosity of novel copper oxide coconut oil nanofluid was studied by Rashin and Hemalatha [15]. Experimental viscosity studies have shown that most of nanofluids have a nonNewtonian behavior. Rheological characteristics of nonNewtonian nanofluids were investigated by Hojjat et al. [16]. Results showed that nanofluids as well as the base fluid exhibit pseudoplastic (shear thinning) behavior. They found that the apparent viscosity of nanofluids decreased by increasing shear rate and increasing volume fraction.
The objective of the present study is to investigate the rheological behavior of some nanofluids. In few studies it was found that the viscosity of the fluid could be decreased by nanoparticle addition at some vol% and shear rates [17]. The present study was conducted to prove the idea. The powerlaw model was used to model the nonNewtonian behavior of the fluids. Two fitting parameters, the powerlaw index and the consistency index, have been calculated for different nanofluids in several concentrations. The base fluid used throughout this investigation was the heat transfer oil. Al2O3, TiO2 and CuO nanoparticles were used in the nanofluids at 1 and 2 wt.% concentrations. The results are compared with classical models, and their discrepancy is discussed.
Nanofluids are made form Al_{2}O_{3}, TiO_{2} and CuO nanoparticles of weight fractions 1% and 2% using the twostep method. Transmission electron microscope (TEM) image of CuO nanoparticles is shown in Fig. 1. The heat transfer oil used in this study was purchased from Behran Oil Company. According to ASTM D1290 standard the density of the above mentioned oil at the temperature of 15.6 ^{o}C is 870 Kg/m^{3} .Since these suspensions are not stable systems, stabilizing methods should be utilized to gain more stable suspensions. In the present study an appropriate amount of SDBS surfactant, Sodium Dodecyl Benzene Sulfonate, recommended by [18] was added to the mixture. Then, the solid particles deagglomerated by intensive ultrasonic device (Ms Hielscher model UP400S) for 3 h after mixing with the base fluid by magnetic forced agitation. In order to
Figure 1. Transmission electron microscopic image of CuO particles with a nominal size of 40 nm 
Table 1 Physical property of nanoparticles.
Shape of particles 
Nanoparticle’s Purity 
mean diameters (nm) 
Particle 
Spherical particles 
>99% 
40 
CuO 
Spherical particles 
>99% 
15 
Al_{2}O_{3} 
Spherical particles 
>99.9% 
20 
TiO_{2} 
prevent fluid from vaporizing, samples remained in mixture of ice and water during sonication. Nanofluids produced were found to be very stable and there was no visually observable sedimentation after more than a week. Table 1 contains other properties of the nanoparticles.
Viscosity of the nanofluids was measured by Anton Paar rheometer (model: MCR301, figure 2). The viscometer drives the spindle immersed into the sample holder containing the test fluid sample and provides a rotational speed that can be controlled to vary from 0.01 to 250 rpm. The viscometer has been tested and calibrated with the calibration fluid provided by Brookfield Engineering Laboratories. Meanwhile, all the measurements are performed under steady state conditions. Each experiment was repeated three times. The overall standard deviation of the results was below
Figure 2. Anton paar rheometer 
6%.
3. Results and discussions
Figure 3 shows the variation of apparent viscosity with respect to the strain rate (s^{1}). It can be seen that the viscosity of the base fluid (heat transfer oil) increases with an increase in the strain rate. On the contrary, the viscosity of nanofluid decreases with an increase in the strain rate by adding the nanoparticles (Al_{2}O_{3}, TiO_{2}, and CuO) to the base fluid. So, the base fluid (heat transfer oil) is a shear thickening (dilatant) fluid while the nanofluid, based on the heat transfer oil, is a shear thinning (pseudoplastic) fluid. In other words, the behavior of fluid changes by adding the nanoparticles.
The slope and yintercept of profiles in figure 3 indicate the Powerlaw index and consistency index of powerlaw model according to equation 1 [16, 1923]:
(1) 
μ is apparent viscosity and is shear rate. r and n are the constants of powerlaw model. Figure 4 shows the linear behavior of the measured viscosity of the base fluid. From the slope and yintercept of the fitted line in figure 4, the Powerlaw index and consistency index of the base fluid can be obtained which are 1.032 and 0.053, respectively. This procedure could be done for the Al_{2}O_{3}–oil nanofluid at 1 and 2 wt% which is presented in figure 5.
Figure 3. Logarithmic variation of the measured viscosities of the fluids 
Figure 4. Variation of viscosity of the pure base fluid (heat transfer oil) with respect to the imposed strain rate

An increase in the consistency index of the fluids with the presence of nanoparticles is illustrated in figure 7. The maximum increase in the consistency index is for the Al_{2}O_{3} nanoparticles the same as what could be seen for the powerlaw index in figure 6. The consistency index increases 4, 14 and 20% at 2 wt% concentration for CuOoil, TiO_{2}oil and Al_{2}O_{3}oil nanofluids, respectively.
The ratio of dynamic viscosity of nanofluid to the base fluid is shown in figures 8 and 9 at 1 and 2 wt% concentration of nanoparticles. This ratio decreases with an increase in the strain rate. Al_{2}O_{3}oil nanofluid has the maximum tendency in reduction in the viscosity. The dynamic viscosity ratio decreases about 12% at 1 wt% of Al_{2}O_{3}oil nanofluid while this reduction is about 8% at 1 wt% of TiO_{2}oil and CuOoil nanofluids, respectively. Moreover, the slope of reduction in the dynamic viscosity ratio for the Al_{2}O_{3}oil is the highest. The reduction in the
Figure 6. Effects of nanoparticle concentration on the powerlaw index 
Figure 7. Effects of nanoparticle concentration on the consistency index 
viscosity ratio is 10% for the CuOoil nanofluid at 2 wt% concentration. This reduction for the TiO_{2}oil and Al_{2}O_{3}oil nanofluids is only about 7%. Reduction in the viscosity of the fluid with an increase in the strain rate shows a shear thinning behavior. This shearthinning behavior is due to the nanoparticle agglomeration clusters which break down when the shear rate is increased and, as a result, the apparent viscosity decreases [16]. Another reason in the reduction of viscosity is the abruption in the base fluid structure. Nanoparticles could break off the connection between the base fluid molecules and act as a connector between fluid layers. Also, Brownian nature of particles and its differences form the Brownian characteristics of the base fluid could be the other reason. These differences force fluid and nanoparticles to move relatively. As a result of this relative motion, the connection between fluid layers and also between the fluid and nanoparticle is continuously constructed and broken off.
The correlations of Brinkman [24] and Einstein are used for comparison. Einstein proposed correlation (2) for the suspensions of solid particles in fluids. In 1999 Drew and Passman [25] showed that the Einstein’s relation can be used for the mixtures of micrometer particles with low (<2) volumetric concentrations.
(2) 
Brinkman [24] modified the Einstein relation to be applicable in the denser mixtures:
Figure 8. Effects of nanoparticle addition and shear rate on the viscosity of the nanofluids (at 1 wt.% concentration) 
Figure 9. Effects of nanoparticle addition and shear rate on the viscosity of the nanofluids (at 2 wt.% concentration) 

(3) 
None of the correlations can predict the experimental results in figures 8 and 9. At 1 wt%, all of the nanofluids have a lower value for the dynamic viscosity ratio than the correlations (2) and (3). At 2 wt%, the Al_{2}O_{3}oil nanofluid has values the same as the correlations at low shear rates. But, at high shear rates, the experimental results fall below the prediction of the correlations. Moreover, the discrepancies of the experimental results and the correlations are a function of shear rate. It seems that further studies on the viscosity of nanofluids are necessary to develop a correlation for the behavior of nonNewtonian nanofluids and maybe for the reduction in the viscosity. Also, this study shows that the viscosity of nanofluids could be lower than the base fluid at some concentrations and shear rates, on the contrary of all available correlations. The reduction in the viscosity is not a new finding. This behavior was not seen in the previous literatures except for some few investigations like Kole and Dey [17]. But, it wasn’t emphasized that the reduction in the fluid viscosity at some vol% of nanoparticle presence occurred in the study of Kole and Dey [17]. Also, the available correlations proposed for the effective viscosity cannot predict the reduction of the viscosity in presence of nanoparticles. However, the reduction in the effective viscosity has very important applications, especially in lubrications and pumping power. The lower the value of viscosity, the lower would be the loss of energy.
4. Conclusion
The aim of the present work was the study of viscosity of some nanofluids based on the nonNewtonian powerlaw model. The heat transfer oil was selected as the base fluid. Al_{2}O_{3}, TiO_{2}, and CuO nanoparticles were added to the base fluid at 1% and 2% weightedaverage concentrations. The highlights of the study are stated as follows:
The results showed that the effective viscosity of nanofluid could be reduced at some wt% and strain rates. The reduction in the viscosity by the use of nanoparticle would have important applications in lubrication and pumping. By reducing the viscosity of the fluid, the energy loss would decrease.
Appendix
The constitutive relation for the stress tensor of the fluids is presented in this section. The balance of momentum gives [19]:
(1) 
v is velocity vector, ρ is density, S is the Cauchy stress tensor, M is body force, and g is gravitational force. Based on the nature of the fluids, the Cauchy stress tensor is divided into two parts; deviatoric and spherical parts [20  26]:
(2) 
I is the unit tensor, τ is the shear stress tensor, and is the pressure of the fluid. The ReinerRivlin theory is used in defining the Cauchy stress tensor. Form the point, the Cauchy stress tensor is a function of the rate of deformation tensor (D) [23, 24]:
(3) 
and are material functions (or material coefficients). The rate of deformation tensor (D) is the symmetric part of the gradient of velocity vector and is defined as follows:
(4) 
has a constant value and , , and are the invariants of the D [23, 24]:
(5) 
The undefined coefficient ( and ) should be defined. A comparison between relations 2 and 3 gives that the is the pressure of the fluid. Through experimental investigations it was found that the dependency of to the third invariant of the rate of deformation tensor ( ) is weak. In other words, is a function of the second invariant of D ( ). Form the second law of thermodynamic it holds that should have positive values [27, 28]. So, one can write:
(6) 
The second invariant of D has a negative value, on the other hand is a positive function, and regarded to the experimental investigations the equation (6) could be rewritten in the form of equation (7). Fluid which follows this relation is called generalized Newtonian fluids (incompressible).
(7) 
is a tensor function and modeling should be used to determine it. The common model is the powerlaw model which is used in the present study. In the powerlaw model, has the following form [2330]:
(8) 
in which n is powerlaw index and r is consistency index. Therefore, the Cauchy stress tensor will be [23, 24, 30]:
(9) 
The values of n in the constitutive relation (9) define the behavior of the fluid. n<1 is the case of shear thinning (pseudoplastic) fluids. The resistance of the fluid against the shear rate (viscosity of the fluid) decreases with increased shear rate. n=1 is the case of Newtonian fluids. In this case the viscosity of the fluid is independent of the imposed deformation rate. Most of the fluids are categorized in these two cases (Newtonian and shear thinning fluids). However, some of the fluids have a different behavior and are called shear thickening (dilatant) fluids. For shear thickening fluids n is greater than one. The viscosity of this type of fluids increases with an increase in the shear rate [23, 28].
References
[1]. B.C. Pak, Y.I. Cho, “Hydrodynamic and Heat Transfer Study of Dispersed Fluids with Submicron Metallic Oxide Particles,” Experimental Heat Transfer, 11 (2), 151170 (1998).
[2]. S.K. Das, N. Putra, W. Roetzel, “Pool Boiling Characteristics of Nano fluids,” International Journal of Heat and Mass Transfer, 46 (5), 851862 (2003).
[3]. M. Mirzaei, M. Dehghan, “Investigation of flow and heat transfer of nanofluid in microchannel with variable property approach”, Heat Mass Transfer, 49, 18031811 (2013).
[4]. S.Z. Heris, S.G. Etemad, M.N. Esfahany, “Experimental Investigation of Oxide Nanofluids Laminar Flow Convective Heat Transfer,” Int. Comm. Heat Mass Trans., 33(4), 529535 (2006).
[5]. D.P. Kulkarni, D.K. Das, G.A. Chukwu, “Temperature Dependent Rheological Property of Copper Oxide Nanoparticles Suspension (Nanofluid),” J. Nanoscience Nanotechnology, 6(4), 11501154 (2006).
[6]. T. Phuoc, M. Massoudi, R. Chen., “Viscosity and thermal conductivity of nanofluids containing multiwalled carbon nanotubes stabilized by chitosan,” International Journal of Thermal Sciences, 50, 1218 (2011).
[7]. A. Einstein, “Eine neue Bestimmung der Molekuldimension,” Annalen der Physik, 19, 289306 (1906).
[8]. J. Yang, F. Li, W. Zhou, Y. He, B. Jiang , “Experimental investigation on the thermal conductivity and shear viscosity of viscoelasticfluidbased nanofluids”, International Journal of Heat and Mass Transfer, 55, 31603166 (2012).
[9]. P. Ravi, S. David, W. Jinlin, “Measurement of nanofluid viscosity and its implications for thermal applications”. Applied Phys. Lett., 89 (13), 133108 1331083 (2006).
[10]. S. Lee, S. Park, S. Kang, I. Bang, J. Kim , “Investigation of viscosity and thermal conductivity of SiC nanofluids for heat transfer applications,” International Journal of Heat and Mass Transfer, 54, 433438 (2011).
[11]. A.Utomo, H. Poth, P. Robbins, A. Pacek , “Experimental and theoretical studies of thermal conductivity, viscosity and heat transfer coefficient of titania and alumina nanofluids” ,International Journal of Heat and Mass Transfer, 55, 7772–7781 (2012).
[12]. H. Chen, Y. Ding, Y. He, C. Tan, “Rheological behavior of ethylene glycol based titania nanofluids”, Chem. Phys. Lett. 444, 333–337 (2007).
[13]. H. Chen, S. Witharana, Y. Jin, C. Kimd, Y. Ding, “Predicting thermal conductivity of liquid suspensions of nanoparticles based on rheology”, Particuology 7, 151–157 (2009).
[14]. M. Kole, T.K. Dey, “Viscosity of alumina nanoparticles dispersed in car engine coolant,” Experimental Thermal and Fluid Science, 34, 677–683 (2010).
[15]. M. Nabeel Rashin, J. Hemalatha, “Viscosity studies on novel copper oxide–coconut oil nanofluid,” Experimental Thermal and Fluid Science, 48, 67–72 (2013).
[16]. M. Hojjat, S.Gh. Etemad, R. Bagheri, J. Thibault, “Rheological characteristics of nonNewtonian nanofluids: Experimental investigation”, International Communications in Heat and Mass Transfer, 38, 144–148 (2011).
[17]. M. Kole, T.K. Dey, “Thermophysical and pool boiling characteristics of ZnOethylene glycol nanofluids”, International Journal of Thermal Sciences 62, 6170 (2012).
[18]. M.T. JamalAbad, A. Zamzamian, M. Dehghan, Experimental studies on the heat transfer and pressure drop characteristics of Cuwater and Alwater nanofluids in a spiral coil, Experimental Thermal and Fluid Science, 47, 206–212 (2013).
[19]. F.M. White, “Viscous Fluid Flow”, third ed. McGrawHill, New York, 2006.
[20]. M. Dehghan, H. Basirat Tabrizi, On nearwall behavior of particles in a dilute turbulent gas–solid flow using kinetic theory of granular flows, Powder Technology, 224, 273–280 (2012).
[21]. M. Dehghan, H. Basirat Tabrizi, Turbulence effects on the granular model of particle motion in a boundary layer flow, Canadian Journal of Chemical Engineering, 92, 189–195 (2014).
[22]. M. Dehghan, H. B. Tabrizi, Effects of coupling on turbulent gasparticle boundary layer flows at borderline volume fractions using kinetic theory, Journal of Heat and Mass Transfer Research, 1, 18 (2014).
[23]. M. Dehghan, M. Mirzaei, A. Mohammadzadeh, Numerical formulation and simulation of a nonNewtonian magnetic fluid flow in the boundary layer of a stretching sheet, Journal of Modeling in Engineering, 11 (34), 7382 (2013).
[24]. G. Astarita and G. Marrucci, “Principles of NonNewtonian Fluid Mechanics”, McGrawHill (UK), (1974).
[25]. D.C. Leigh, “NonNewtonian fluids and the second law of thermodynamics”, Physics of Fluids, 5, 501–502 (1962).
[26]. B.D. Coleman, H. Markovitz, W. Noll, “Viscometric Flows of NonNewtonian Fluids”, SpringerVerlag, (1966).
[27]. H.C. Brinkman, “The Viscosity of Concentrated Suspensions and Solution”, J. Chem. Phys. 20, 571581(1952).
[28]. D.A. Drew, D.A. Passman, “Theory of Multicomponent Fluids”, Springer, Berlin, (1999).
[29]. H.I. Andersson, B.S. Dandapat, “Flows of a power law fluid over a stretching sheet”. Stability Appl Anal Continuous Media,1, 339–347 (1991).
[30]. M. Dehghan, M. Mirzaei, M.S. Valipour, S. Saedodin, Flow of a nonNewtonian fluid over a linearly moving sheet at a transient state; new similarity variable and numerical solution scheme, Journal of Modeling in Engineering, (2014) (accepted in press).