Rheological and volumetric properties of TiO2-ethylene glycol nanofluids

Abstract

Homogeneous stable suspensions obtained by dispersing dry TiO2 nanoparticles in pure ethylene glycol were prepared and studied. Two types of nanocrystalline structure were analyzed, namely anatase and rutile phases, which have been characterized by scanning electron microscopy. The rheological behavior was determined for both nanofluids at nanoparticle mass concentrations up to 25%, including flow curves and frequency-dependent storage and loss moduli, using a cone-plate rotational rheometer. The effect of temperature over these flow curve tests at the highest concentration was also analyzed from 283.15 to 323.15 K. Furthermore, the influence of temperature, pressure, nanocrystalline structure, and concentration on the volumetric properties, including densities and isobaric thermal expansivities, were also analyzed.

Cabaleiro et al. Nanoscale Research Letters 2013, 8:286 http://www.nanoscalereslett.com/content/8/1/286

NANO REVIEW Open Access

Rheological and volumetric properties of TiO2-ethylene glycol nanofluids

David Cabaleiro, Mara J Pastoriza-Gallego, Carlos Gracia-Fernndez, Manuel M Pieiro and Luis Lugo*

Abstract

Homogeneous stable suspensions obtained by dispersing dry TiO2 nanoparticles in pure ethylene glycol were prepared and studied. Two types of nanocrystalline structure were analyzed, namely anatase and rutile phases, which have been characterized by scanning electron microscopy. The rheological behavior was determined for both nanofluids at nanoparticle mass concentrations up to 25%, including flow curves and frequency-dependent storage and loss moduli, using a cone-plate rotational rheometer. The effect of temperature over these flow curve tests at the highest concentration was also analyzed from 283.15 to 323.15 K. Furthermore, the influence of temperature, pressure, nanocrystalline structure, and concentration on the volumetric properties, including densities and isobaric thermal expansivities, were also analyzed.

Keywords: Nanofluid; TiO2/EG; Anatase; Rutile; Non-Newtonian; Nanoparticles; Density; Thermal expansivity coefficient

Review

BackgroundNowadays, it is well known that the novel proposal of nanofluids represents a valuable way for the development of the heat transfer fluids currently available. Thus, nanofluids have recently emerged with new potential applications in heat exchangers or cooling devices, being widely used in many engineering applications as electronics cooling, vehicle engines, nuclear reactors, energy efficiency enhancers, food industry, air conditioning, refrigeration, and biomedicine [1-4]. As an example, it has been shown that by using nanofluids in radiators, pumps, or compressors in cars, the aerodynamic charge could be reduced, producing fuel savings up to 6% [5].Therefore, with the aim to improve the heat transfer properties of nanofluids, a considerable amount of research efforts are being devoted to the analysis of their thermal conductivity and convective heat transfer properties. Although it is possible to tailor nanofluids exhibiting negative thermal conductivity enhancement, or a decrease in the effective thermal conductivity of the dispersion if compared with that of the base liquid [6], in most cases, nanofluids exhibit a significant enhance-

ment in thermal conductivity. Therefore, nanofluids are expected to provide optimized convective heat transfer coefficients. However, this type of nanocolloidal dispersion affects also other thermophysical properties than thermal conductivity. Concerning the concentration dependence of nanofluids, a revision of the literature shows, besides the increase in thermal conductivity, decreases of heat capacity and a noticeable increase of density and viscosity, including the possibility of a non-Newtonian behavior. All these properties affect significantly the convective heat transfer coefficient. In addition, as the relation between this coefficient and the involved thermophysical properties could not follow classical laws, it is essentially required to determine accurately their trend with concentration, temperature, and/or pressure.

Recently, Huminic and Huminic [2] have reported a review on the application of nanofluids in various types of heat exchangers as plate, shell and tube, compact, and double pipe heat exchangers. The authors concluded that both the thermophysical properties and type of flow inside the heat exchanger played important roles in the efficiency of the nanofluid as a coolant. Moreover, in most practical applications, the heat transfer fluid is not stationary [3], and consequently, the analysis of the rheological properties is also essential to appropriately

* Correspondence: mailto:[email protected]

Web End [email protected] Departamento de Fsica Aplicada, Facultade de Ciencias, Universidade de Vigo, Vigo 36310, Spain

2013 Cabaleiro et al.; licensee Springer. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0

Web End =http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Cabaleiro et al. Nanoscale Research Letters 2013, 8:286 Page 2 of 13 http://www.nanoscalereslett.com/content/8/1/286

determine the increments on the average heat transfer coefficient of the flowing system, which generally increases with the concentration of nanoparticles as well as with the Reynolds number [2]. Numerical results [7] indicate that high-concentration nanofluids of TiO2 or

Al2O3 in water exhibit higher heat transfer enhancements and also higher pressure drops. On the other hand, Peyghambarzadeh et al. [4] have experimentally demonstrated, using water- and ethylene glycol (EG)-based nanofluids as cooling agents inside flat aluminum tubes of a car radiator, that the heat transfer behaviors of the nanofluids were highly dependent on particle concentration and flow conditions and otherwise weakly temperature dependent. From the results of Huminic and Huminic [2], it can be concluded that homogeneously dispersed and stabilized nanoparticles enhance the forced convective heat transfer coefficient of the base fluid in a range of 3% to 49%, observing a greater increase with increasing temperature and nanoparticle concentration. Therefore, a proper balance between the heat transfer enhancement and the pressure drop penalty, together with viscosity behavior, should be taken into account when seeking an appropriate nanofluid for a given application.

In addition to the knowledge of the cited rheological behavior, the volumetric properties including the isobaric thermal expansivity coefficient play as well an important role in many heat removal systems involving natural convection. The thermal expansivity coefficient is needed to apply nanofluids in engineering-scale systems [8,9], and this property is usually negligible for metallic oxide particles if compared to that of the base fluids as EG or water. Hence, it is often presumed that this coefficient should decrease with rising concentration of nanoparticles as we have previously reported [10]. Nevertheless, some works [8,9] have found the opposite behavior of the one resulting from considering the fluids to behave separately in the mixture for the case of water-based Al2O3 nanofluids. This is one of the singular properties of nanofluids that would find a remarkable application in many heat extraction systems using natural convection as a heat removal method [11]. Therefore, more attention should be paid to this magnitude with the goal to understand the complex interaction of nanoparticles with the base fluid molecules, and it could be also a powerful additional tool to characterize nanofluids.

In this work, we focus our attention on the volumetric and rheological behaviors of the suspension of two nanocrystalline forms of TiO2 nanoparticles, anatase and rutile, dispersed in pure EG as the base fluid. The influence of the nanocrystalline phase, temperature, pressure, and concentration on the isobaric thermal expansivity coefficient is also analyzed, looking for a verification of the surprising results for different nanofluids found by

Nayak et al. [8,9]. In addition to the reasons cited, the selection here of TiO2/EG nanofluids is inspired also on several other arguments. First, EG can be used over a wide temperature range. Then, an enhancement in the overall heat transfer coefficient of up to 35% in a compact reactor-heat exchanger, with a limited penalty of increase in pressure drop due to the introduction of nanoparticles, has been reported for TiO2/EG nanofluids [3]. Moreover,

TiO2 is a safe and harmless material for human and animals if compared with other nanomaterials [12]. From a production perspective, these nanoparticles are easily obtained because they are readily produced in large industrial scales. Concerning their physicochemical profile, they have an excellent stability when dispersed in a fluid even without stabilizer addition, and metal oxide nanoparticles are chemically more stable than their metallic counterparts [13]. Finally, remarkably few works are found in the literature [3,14,15] devoted to the study of thermal or rheological properties of TiO2/EG nanofluids, and up to our knowledge, their volumetric and viscoelastic properties have never been reported.

The experimental density of stable and homogeneous TiO2/EG nanofluids at percent mass concentrations (wt.%) of 1.00, 1.75, 2.50, 3.25, and 5.00, which correspond in percent volume (vol.%), respectively, of 0.29,0.51, 0.74, 1.04, and 1.51 for anatase and 0.26, 0.47, 0.67,0.94, and 1.36 for rutile, in wide pressure (from 0.1 to 45 MPa) and temperature (from 283.15 to 343.15 K) ranges was analyzed. From these density data for anatase titan-ium dioxide-EG nanofluids (A-TiO2/EG, from now on, for the sake of brevity) and rutile titanium dioxide-EG nanofluids (viz. R-TiO2/EG) [16], the derived thermal expansion and thermal compressibility coefficients were studied. Moreover, we have carried out a rheological study on samples of A-TiO2/EG and R-TiO2/EG nanofluids at mass concentrations of 5.00, 10.00, 15.00, 20.00, and 25.00 wt.%, which correspond to 1.51, 3.13, 4.88, 6.77, and 8.83 vol.% for A-TiO2/EG and to 1.36, 2.83, 4.43, 6.16, and 8.08 vol.% for R-TiO2/EG, respectively. The effect of the structure of nanoparticles, rutile and anatase, on linear and non-linear tests was analyzed on these samples, and the influence of the temperature was carried out over a temperature range of 283.15 to 333.15 K for the 25 wt.% concentration in both structures.

Several works in the literature have focused on water-or water + EG-based TiO2 nanofluids [13,17-24]. Bobbo et al. [17] and Penkavova et al. [18] studied the viscosity of TiO2/water nanofluids observing a Newtonian behavior for all compositions, while He et al. [13] concluded that aqueous TiO2 nanofluids, with anatase phase and a small amount of rutile phase, show a shear thinning behavior where the shear viscosity tends to be constant at shear rates above 100 s1 and also that the pressure drop of these nanofluids is very close to that of the base

Cabaleiro et al. Nanoscale Research Letters 2013, 8:286 Page 3 of 13 http://www.nanoscalereslett.com/content/8/1/286

liquid. Nevertheless, Tseng and Lin [24] have investigated the rheological behavior of suspensions of anatase TiO2 nanoparticles in water (0.05 to 0.12 vol.%), reporting a pseudoplastic flow for most of the shear rates examined, from 10 to 1,000 s1. Moreover, their tests suggest a time-dependent phenomenon, attributing to these suspensions a thixotropic response [24]. Several authors [19-23] have studied thermal conductivity enhancements, higher than 20% [21], increasing the nano-particle concentration. Concerning volumetric studies in TiO2/water nanofluids, only the work by Setia et al. [20]

can be cited, where the specific volume for several concentrations up to 0.75 vol.% of TiO2 nanoparticles for several temperatures is reported, finding significant deviations from the additive rule [25] for the samples with volume fractions higher than 0.5 vol.%.

Nevertheless, as pointed out above, few studies were focused on the thermophysical or rheological behavior of TiO2/EG nanofluids [3,14,15]. Fan et al. [3] determined the thermal conductivity at 303 K for the concentrations 0.5, 2.0, and 4.0 wt.% (corresponding respectively to 0.10, 0.43, and 0.86 vol.%) for TiO2/EG nanofluids and their corresponding viscosity in the shear rate range of 1 to 3,000 s1, confirming a Newtonian behavior and the expected increase of viscosity with nano-particle concentration. Chen et al. [14] have also found a Newtonian behavior for TiO2/EG nanofluids containing0.5, 1.0, 2.0, 4.0, and 8.0 wt.% spherical nanoparticles at 293.15 to 333.15 K and a relative viscosity dependent on particle concentration in a non-linear manner without temperature dependence. On the other hand, Lee et al.[15] have determined temperature-independent thermal conductivity enhancements up to 16% for 5.5 vol.% TiO2/EG nanofluids constituted by nanoparticles with rutile and anatase phases. On the other hand, to our knowledge, no evidence on non-Newtonian behavior for TiO2/EG nanofluids, or studies about their volumetric behavior, including densities, isothermal compressibility, and isobaric thermal expansivity coefficients, have been reported so far in the literature. Hence, there is a key need to address this issue.

Methods

Homogeneous and stable suspensions were prepared by dispersing dry TiO2 nanoparticles in pure EG. Two types of TiO2 powder, corresponding to the pure nanocrystal-line anatase and rutile phases, whose descriptions are shown in Table 1, were employed. Although rutile is the stable phase for bulk TiO2, the colloidal phase preparation methods for TiO2 generally favor the anatase structure [26,27]. Both types of nanoparticles were supplied by SkySpring Nanomaterials, Inc. (Houston, TX, USA) with a reported average size of 10 to 30 nm for rutile and 10 to 25 nm for anatase, with a chemical purity of 99.5% for both cases, while ethylene glycol with a mass purity of 99.5% was supplied by Sigma-Aldrich (St. Louis, MO, USA). With the aim to characterize the morphology of these nanomaterials, both types of TiO2 nanoparticles were characterized using the scanning electron microscopy (SEM) technique, obtaining the images with a JEOL JSM-6700 F field emission gun-SEM (Akishima-shi, Japan) operating at an acceleration voltage of 20 kV in a backscattering electron image (yttrium aluminum garnet-type detector). This device incorporates an energy-dispersive X-ray (EDS) spectrometer that was used to chemically characterize the samples. SEM samples were prepared by the deposition of a drop of nanopowder dispersed in methanol (analytical grade) on top of a silica film, dried under atmospheric conditions. SEM images for both types of nanocrystalline structures are shown in Figure 1. The magnification of close agglomerates in micrometers (Figure 1b,d) allows identifying the individual nanoscale globular or nearly spherical particles for anatase and rutile. Average particle sizes were estimated using SEM micrographs by counting a minimum of 100 particles, obtaining values of 35 17 nm for anatase and 47 18 nm for rutile. Figure 2 shows the chemical composition of the samples, obtained from the EDS spectra, determined from the area displayed in Figure 2a,c and represented in Figure 2b,d. The analysis of anatase nanoparticles shows that only Ti and O elements are detectable (Figure 2b), while for rutile, an amount inferior to 1% by mass of

Table 1 Material description

Material Supplier Mass purity (%) Medium size (nm) Crystalline structure Anatase titanium

dioxide (A-TiO2)

SkySpring

Nanomaterials

99.5 35 17 Tetragonal

Rutile titanium

dioxide (R-TiO2)

SkySpring

Nanomaterials

99.5 47 18 Tetragonal

Cabaleiro et al. Nanoscale Research Letters 2013, 8:286 Page 4 of 13 http://www.nanoscalereslett.com/content/8/1/286

a) b)

c) d)

Figure 1 SEM images of dry A-TiO2 and R-TiO2 nanoparticles. SEM images of anatase nanoparticles at two magnifications: 50,000 (a) and

200,000 (b), and rutile nanoparticles at two magnifications: 50,000 (c) and 200,000 (d).

Si is present, as shown in Figure 2d, probably due to the silica support. No relevant amounts of other compounds were identified for the samples studied.

The preparation of the nanofluid was carried out using the two-step method at the mass concentrations of 1.00,1.75, 2.50, 3.25, and 5.00 wt.% for volumetric measurements, whereas 5.00, 10.00, 15.00, 20.00, and 25.00 in wt.% concentrations were used for rheological tests, without adding any surfactant, in order to study the effect of nanoparticle aggregation. The uncertainty in the mass compositions for the different studied nanofluids ranges from 0.003% to 0.02%, increasing with the mass concentration. Subsequently, the nanofluids were dispersed by ultrasonic homogenization using a Bandelin Sonoplus HD 2200 (Bandelin Electronic, Berlin, Germany) for 16 min to prevent aggregation. More details about sonication methods have been previously published [28].

Concerning the characterization of the volumetric behavior of the cited R-TiO2/EG and A-TiO2/EG nanofluids, density measurements were experimentally carried out at concentrations up to 5% in mass fraction from atmospheric pressure up to 45 MPa and from

278.15 to 363.15 K. Temperature and pressure were measured within uncertainties of 0.02 MPa and 0.02 K for pressure and temperature, respectively. Density data were obtained from the period values measured using a commercially available vibrating tube densimeter (Anton Paar DMA 512P, Graz, Austria) with an estimated uncertainty of 5 104 g cm3 [29] over the whole pressure and temperature range. More details about the procedure, calibration, temperature, and pressure control can be found in our previous works [10,30,31].

Rheological properties of R-TiO2/EG and A-TiO2/ EG nanofluids were determined using a rotational Physica MCR 101 rheometer (Anton Paar, Graz, Austria), equipped with a cone-plate geometry with a cone diameter of 25 mm and a cone angle of 1. The cone went down to an imposed gap of 0.048 mm from the plate and covered the whole sample for all tests. The measurement consists of imposing the shear stress to the sample and recording the related shear rate. Temperature is controlled with a Peltier P-PTD 200 (Anton Paar, Graz, Austria), placed at the lower plate, with a diameter of 56 mm without

Cabaleiro et al. Nanoscale Research Letters 2013, 8:286 Page 5 of 13 http://www.nanoscalereslett.com/content/8/1/286

a) b)

c) d)

Figure 2 EDS images and microanalysis of TiO2 nanoparticles. EDS images of A-TiO2 (a) and R-TiO2 (c) nanoparticles, and microanalysis from

the area within the rectangle shown in EDS images for A-TiO2 (b) and R-TiO2 (d) nanoparticles.

groove. The linear and non-linear tests were developed from torques of 0.1 Nm in the temperature range of 283.15 to 323.15 K, each 10 K. A constant amount of 110 l of sample was considered [32] for the analysis and was placed on the Peltier plate. Nonlinear and linear viscoelastic experiments were carried out with the objective to analyze both relatively large deformations and small-amplitude oscillatory shear. Thus, the flow curves of the samples studied and the frequency-dependent storage (G) and loss (G) moduli were determined. More details about the experimental setup and operating conditions can be found in our previous papers [10,32,33].

Results and discussion

Volumetric propertiesThe density values of both sets of nanofluids, A-TiO2/

EG and R-TiO2/EG, at mass fractions up to 5 wt.% were experimentally measured at pressure up to 45 MPa in a wide temperature range of 278.15 to 363.15 K along eight isotherms. Table 2 reports the experimental density data for both nanofluids. The density values range from1.0627 g cm3 for pure EG, at 0.1 MPa and 363.15 K, up to 1.1800 g cm3 for A-TiO2/EG nanofluids and 1.1838 g cm3 for R-TiO2/EG nanofluids at 5 wt.%, p = 45 MPa, and T = 278.15 K. At equivalent temperature, pressure

and concentration, the density values of the A-TiO2/EG are lower than those of R-TiO2/EG, excepting the 1 wt.% sample, for which they agree to within the experimental uncertainty. Density values increase with nanoparticle concentration as expected, as shown in Figure 3a where the increments in relation to the base fluid reference value at different concentrations are shown, with higher increments also for the rutile nano-crystalline structure, reaching values of 3.8%. We have found that these increments with concentration are almost temperature and pressure independent. For a given concentration, density data show pressure and temperature dependences similar to the base fluid, increasing with pressure and decreasing with temperature. The average percentage density increments with pressure range between 1.5% at the lowest temperature and 2% at the highest temperature. On the other hand, the average percentage density diminutions with temperature at different pressures are gathered in Figure 3b, showing decreases between 5% and 5.4%. These temperature variations are very similar for both nanocrystalline structures and the base fluid, as can be appreciated in this figure.

With the aim to report a generalized temperature and pressure correlation of the volumetric behavior of the measured base fluid and nanofluids, the specific volumes

Cabaleiro et al. Nanoscale Research Letters 2013, 8:286 Page 6 of 13 http://www.nanoscalereslett.com/content/8/1/286

Table 2 Density (), isobaric thermal expansivity (p), and isothermal compressibility (T) of A-TiO2/EG and R-TiO2/EG nanofluids

p (MPa)

(gcm3) 104p (K1) 104T (MPa1) T = 283.15

K

T = 313.15

K

T = 343.15

K

T = 283.15

K

T = 313.15

K

T = 343.15

K

T = 283.15

K

T = 313.15

K

T = 343.15

K

Base fluid (EG) 0.10 1.1202 1.0989 1.0772 6.31 6.52 6.731.00 1.1206 1.0993 1.0776 6.30 6.51 6.72 3.52 3.89 4.3420.00 1.1279 1.1073 1.0861 6.09 6.27 6.43 3.34 3.69 4.0840.00 1.1353 1.1152 1.0950 5.89 6.03 6.14 3.33 3.66 4.0545.00 1.1373 1.1174 1.0973 5.84 5.97 6.07 A-TiO2/EG (1.75

wt.%)

0.10 1.1327 1.1117 1.0901 6.20 6.43 6.661.00 1.1332 1.1121 1.0905 6.20 6.42 6.65 3.35 3.61 3.9720.00 1.1407 1.1200 1.0988 6.06 6.23 6.37 3.38 3.63 4.0040.00 1.1482 1.1280 1.1076 5.92 6.03 6.09 3.27 3.51 3.8545.00 1.1503 1.1300 1.1100 5.89 5.99 6.03 A-TiO2/EG (5.00

wt.%)

0.10 1.1584 1.1366 1.1147 6.42 6.51 6.591.00 1.1589 1.1370 1.1150 6.41 6.50 6.58 3.61 3.96 4.3320.00 1.1667 1.1450 1.1239 6.21 6.29 6.36 3.35 3.65 3.9740.00 1.1745 1.1535 1.1324 6.02 6.08 6.15 3.39 3.70 4.0245.00 1.1766 1.1558 1.1349 5.97 6.03 6.10 R-TiO2/EG (1.75

wt.%)

0.10 1.1339 1.1126 1.0910 6.15 6.41 6.671.00 1.1343 1.1129 1.0914 6.14 6.40 6.66 3.62 0.03 4.5020.00 1.1414 1.1209 1.1001 5.93 6.16 6.39 3.28 3.61 3.9840.00 1.1491 1.1290 1.1093 5.71 5.92 6.12 3.45 3.82 4.2445.00 1.1513 1.1314 1.1113 5.65 5.85 6.04 R-TiO2/EG (5.00

wt.%)

0.10 1.1622 1.1405 1.1184 6.24 6.43 6.631.00 1.1626 1.1409 1.1188 6.23 6.42 6.62 3.52 3.75 4.0720.00 1.1706 1.1489 1.1271 6.10 6.26 6.40 3.41 3.63 3.9340.00 1.1779 1.1570 1.1362 5.98 6.09 6.18 3.34 3.55 3.8345.00 1.1802 1.1592 1.1382 5.95 6.05 6.12

(v = 1/), using the following expression [34], were adjusted to the experimental data:

v T; p

v T; pref

[notdef][notdef]

1 ppref B T; p

[notdef][notdef]; 1

where the reference pressure, pref, was taken as 0.1 MPa.

The dependence of the specific volume at this pressure was described by using the following expression [35,36].

v T; pref

[notdef][notdef]

vref Tref ; pref [notdef][notdef]

e a:b

2

ref

[notdef][notdef]

with TTref ; 2

where a, b, and vref(Tref,pref) are the adjustable parame

ters, vref(Tref,pref) being the specific volume at the refer

ence temperature Tref = 278.15 K and pressure pref = 0.1 MPa. These coefficients, whose values are given in Table 3, were fitted for the base fluid and the different nanofluids with standard deviations lower or equal than2.8 cm3 g1. The bulk modulus B(T, p) was adjusted as a

function of pressure and temperature with the following polynomial:

B T; p

B pref ; T

c d2 ep

f p2 with p ppref : 3

The values of B(pref,Tref), c, d, e, and f were determined

by fitting Equation 1 to all the experimental data at pressures different than pref by a least squares method using

a Marquardt-Levenberg-type algorithm. For the base fluid and all the studied nanofluids, the standard deviations obtained with this correlation are lower than or equal to 1.4 104 cm3 g1, and the coefficients are given in Table 3.

Although viscosity, heat capacity, and thermal conductivity are the main parameters involved in the calculation of the heat transfer rate of a nanofluid, the precise determination of density is also relevant because, as commented above, these properties may be quite

Cabaleiro et al. Nanoscale Research Letters 2013, 8:286 Page 7 of 13 http://www.nanoscalereslett.com/content/8/1/286

4.0

nf 1

0; 4

where is the volumetric fraction of nanoparticles and the subscripts np, 0, and nf refer to the nanoparticles, base liquid, and nanofluids, respectively. The densities of anatase and rutile titanium oxide are, respectively, 3.830 and 4.240 g cm3 [37]. With the aim to evaluate the goodness of this estimation, our experimental values were compared with those predicted using this equation. It was found that this equation overpredicts the density of the nanofluids studied in this work with deviations that it can reach 0.5% for A-TiO2/EG and 0.3% for R

TiO2/EG at the highest concentrations and temperatures, as shown in Figure 4.

Isobaric thermal expansivity, p, and isothermal compressibility, T, coefficients can be determined from specific volume correlations using their respective thermodynamic definitions according the following expressions:

p

a)

100( nf- 0)/ 0

3.0

2.0

1.0

0.0

0.0 1.0 2.0 3.0 4.0 5.0 6.0

wt%

5.6

b)

5.4

100(/)

5.2

5.0

4.8

1 v

v T; p

T

[notdef][notdef]p

v T;pref

T

v T; pref

[notdef][notdef]

B T;p

T

[notdef][notdef]p

p

B T; p

B T; p

p

;

5

0 10 20 30 40 50

p/MPa

Figure 3 Average density variations with nanoparticle concentration and pressure. (a) Enhancement of the density in relation to the base fluid (100 (nf 0)/0) vs. concentration (wt.%) for both nanocrystalline structures. (b) Density diminutions with temperature (100 ()/) vs. pressure (p) for the base fluid and both nanocrystalline structures. Cross mark, base fluid; diamond, A-TiO2/EG; triangle, R-TiO2/EG.

T

1 v

v T; p

p

[notdef][notdef]T

B T; p

B T;p

p

[notdef][notdef]T

p

6

In Table 2, the values calculated for p and T are reported for some temperatures and pressures for the base fluid (EG) and both nanofluids at two different concentrations (1.75 and 5.00 wt.%). The estimated uncertainties for p and T are 4% and 2%, respectively. The p values for both the base fluid and R-TiO2/EG and A-TiO2/EG nanofluids decrease when pressure rises (up to 9.8% for the base fluid) and increase with temperature (up to 6.6% for the base fluid). Concerning the concentration dependence, first, we have found that

B T; p

B T; p

p

:

different from those of the original pure fluid, and it can lead to erroneous mass balances. As we have pointed out, significant variations in density can be achieved when temperature, pressure, concentration, or the type of nanocrystalline structure are analyzed in detail. In order to check some conventional assumptions [3,20], we have determined the ideal nanofluid density from the nanoparticle and base fluid densities according to [25]:

Table 3 Density correlation coefficients and standard deviations () for the base fluid (EG) and the nanofluids

Base fluid

A-TiO2/EG (wt.%) R-TiO2/EG (wt.%)1.00 1.75 2.50 3.25 5.00 1.00 1.75 2.50 3.25 5.00 103a (C1) 0.62714 0.62327 0.61646 0.62116 0.63558 0.64060 0.61708 0.61084 0.62243 0.62955 0.62042 106b (C2) 0.35343 0.30347 0.38267 0.25865 0.17013 0.14365 0.38319 0.43431 0.24473 0.23998 0.32687 104 (cm3 g1) 1.1 1.2 1.2 1.9 1.4 2.8 1.6 1.4 1.8 1.3 1.1 B(pref,Tref) (MPa) 2,875.23 2,813.30 3,016.52 2,732.87 2,840.25 2,798.17 2,796.391 2,782.86 2,744.918 2,619.262 2,865.778 c (MPa C1) 9.1949 8.8432 6.1026 7.7217 10.4348 8.8384 9.8265 9.8347 10.4074 8.6823 5.4028 102d (MPa C2) 0.3779 0.4173 0.2270 0.5231 2.44 1.61 1.61 1.23 2.45 0.89114 1.48e 5.123 5.727 1.559 11.030 7.262 9.430 8.211 13.951 10.066 17.127 3.220 103f (MPa1) 57.3 12.3 49 103.1 50.9 108.5 50.8 190.2 71.4 187.5 12.3 104* (cm3 g1) 0.7 0.8 1.4 0.9 0.9 1.4 0.9 1.0 1.0 1.3 1.2

Cabaleiro et al. Nanoscale Research Letters 2013, 8:286 Page 8 of 13 http://www.nanoscalereslett.com/content/8/1/286

0.2

0.2

a)

b)

100( exp- pred )/ exp

0.0

100( exp- pred)/ exp

0.0

-0.2

-0.2

-0.4

-0.6

-0.4

0 1 2 3 4 5 6

0 1 2 3 4 5 6

wt%

wt%

Figure 4 Percentage deviations between experimental and predicted densities. Deviations between experimental density data (exp) and predicted values (pred) by Equation 4 vs. mass concentration (wt.%) for (a) A-TiO2/EG and (b) R-TiO2/EG nanofluids.

the p values of nanofluids are very similar to or lower than those of EG, achieving decreases up to 1.0% and1.9% for A-TiO2/EG and R-TiO2/EG nanofluids, respectively. These results are opposite to those previously found by Nayak et al. [8,9], reporting a significant increase in this property compared to the base fluid for water-based Al2O3, CuO, SiO2, and TiO2 nanofluids. It should be mentioned that Nayak et al. have determined the isobaric thermal expansivities by measuring the bulk variation with temperature for the samples in a glass flask with a long calibrate stem. Consequently, further studies about this property are still needed on EG- or water-based nanofluids. On the other hand, the T values of the studied samples do not exhibit evident concentration or nanocrystalline structure dependence (or these differences are within the uncertainty). The T values decrease when the pressure rises and increase with the temperature along the isobars for both the base fluid and nanofluid samples, as can be seen in Table 2.

In order to compare the volumetric behavior of the nanofluids with the ideal fluid behavior, excess molar volumes, V Em , were calculated [10,38]. Figure 5 shows an expansive volumetric behavior for both A-TiO2/EG and R-TiO2/EG. This behavior has also been found for other pure EG-based nanofluids, and it is contrary to that presented by nanofluids which use water or EG + water as the base fluid [28]. Excess molar volumes for A-TiO2/EG increase slightly with nanoparticle concentration ranging from 0.03 up to 0.11 cm3 mol1, which correspond to a variation in the molar volume between 3.3% and 14.3%. Concerning R-TiO2/EG, its behavior is closer to ideal, and it is almost concentration independent with a maximum variation in volume of 4.6%. No significant temperature or pressure dependences for this property were found.

Rheological behaviorAs pointed out, only a reduced number of studies about the rheological behavior of nanofluids can be found in the literature, and there are inconsistencies such as

Newtonian and non-Newtonian behaviors reported for the same nanofluid as well as discrepancies in the effects of temperature, particle size, and shape, and high shear viscosity values [40-44]. In this context, a key issue is to obtain nanofluid structural information, and one of the feasible methods is through detailed rheological analyses[45]. In this work, two types of studies have been carried out. Viscosity as a function of shear rate, the so-called flow curve, was determined for both nanofluids at 303.15 K and at five different mass concentrations (5, 10, 15, 20, and 25 wt.%). The applied torques started from 0.1 Nm, covering shear rate ranges from 0.1 to 1,000 s1. Figure 6a,b shows these flow curves for both nanofluids at different concentrations. Unlike the base fluid, both sets of nanofluids present a clear shear thinning (pseudoplastic) non-Newtonian behavior. In the lowest shear rate region, Newtonian plateaus are easily identified as the concentration rises. This non-Newtonian behavior opposes that reported previously by Chen et al. [14] that studied EG-based nanofluids containing0.5 to 8.0 wt.% spherical TiO2 nanoparticles. Chen et al.

0.80

0.60

V mE/cm3mol-1

0.40

0.20

0.00

-0.20

0 0.01 0.02 0.03 0.04

x

Figure 5 Excess molar volumes of nanofluids VEm. Molar excess

volumes vs. molar fraction for different EG nanofluids at 303.15 K and 20 MPa. Filled circle, A-TiO2/EG; filled triangle, R-TiO2/EG; empty

triangle, Fe3O4/EG [38]; empty diamond, Fe2O3/EG [38]; empty circle,

(48-nm ZnO)/EG [39]; empty square, (4.6-nm ZnO)/EG [39].

Cabaleiro et al. Nanoscale Research Letters 2013, 8:286 Page 9 of 13 http://www.nanoscalereslett.com/content/8/1/286

[14] affirmed that a Newtonian behavior is found at a shear rate higher than 0.05 s1. It should be taken into account that our viscosity results for Newtonian EG agree with those of Chen et al. [14] within an average deviation of 1.5% [32]. The controversies found in the literature on rheological studies indicate that the specific properties of the nanoparticles such as shape, structure, and size, and the interaction between the base liquid and nanoparticles can play an essential role in determining the rheological behavior of nanofluids. However, in this case, the main reasons of the different rheological behavior on TiO2/EG nanofluids may be attributed to the following: (1) the range of nanoparticle concentration studied by Chen et al. [14] (<8 wt.%) is lower than those analyzed in this work (<25 wt.%), (2) the range of shear stress studied in this work covers a wider area, and it is here where shear thinning appears, (3) the minimum shear rate which the equipment can reach is decisive to determine the first Newtonian plateau, especially at low nanoparticle concentration, and finally (4) the different stability and aggregation of particles affect flow conditions because the effective mass concentration can be higher than the actual solid mass.

The Ostwald-de Waele model (Power law) was used to describe the experimental shear dynamic viscosity data, , as a function of the shear rate, , in the shear thinning region for each concentration of both sets of nanofluids by using the following expression [46-48]:

K _

n1; 7

where the adjustable parameters K and n are the flow consistency factor and the flow behavior index, respectively. Good adjustments are obtained for all studied nanofluid samples, reaching percentage deviations in shear dynamic viscosity around 3%. At the same mass concentration, the flow behavior index values for RTiO2/EG nanofluids are higher than those for A-TiO2/

EG, as shown in Figure 6c. These n values range from0.27 to 0.72 for A-TiO2/EG and from 0.33 to 0.83 for RTiO2/EG, decreasing near-exponentially when the volume fraction increases, which evidences that the shear thinning behavior is more noticeable when the nanoparticle concentration increases. The n values are similar to those typically obtained for common thermoplastics[49]. It must also be pointed out that although this model offers a simple approximation of the shear thinning behavior, it does not predict the upper or lower Newtonian plateaus [47].

As a further test, the influence of temperature on the flow curves was studied for the highest mass concentration (25 wt.%) for both nanofluids between 283.15 and 323.15 K, as shown in Figure 7a,b, respectively. In these flow curves, we can observe the diminution of viscosity

10

a)

1

/ Pas

0.1

0.01

0.1 1 10 100 1000

.

/ s-1

10

b)

1

/ Pas

0.1

0.01

0.1 1 10 100 1000

.

/ s-1

1

0 0.00 0.02 0.04 0.06 0.08 0.10

c)

0.8

0.6

n

0.4

0.2

Figure 6 Viscosity () vs. shear rate ( _

) of EG/TiO2 nanofluids at

different concentrations. Flow curves for (a) A-TiO2/EG and (b)

R-TiO2/EG at 303.15 K and at different mass concentrations: cross

mark, EG; line, 5 wt.%; circle, 10 wt.%; square, 15 wt.%; diamond, 20 wt.%; triangle, 25 wt.%. (c) Flow behavior index (n) vs. volume fraction () for A-TiO2/EG (filled diamond) and R-TiO2/EG (empty

diamond) at 303.15 K.

Cabaleiro et al. Nanoscale Research Letters 2013, 8:286 Page 10 of 13 http://www.nanoscalereslett.com/content/8/1/286

when the temperature rises, as Chen et al [14] had found in their study between 293.15 and 333.15 K. Nevertheless, the shear viscosities reported in this work show a temperature dependence very influenced by the shear rate value. Moreover, we can observe that the shear viscosity is nearly independent of temperature at a shear rate around 10 s1 for both A-TiO2/EG and R-TiO2/EG nanofluids, which is not the case at a high or low shear rate. On the other hand, at the same concentration and temperature, A-TiO2/EG nanofluids present higher shear viscosities than R-TiO2/EG nanofluids for all shear rates. These viscosity differences increase with concentration. Applying the Ostwald-de Waele model on these flow curves at different temperatures, we have also obtained good results, finding that n values increase with temperature. This may be a result of the temperature effect on the better nanoparticle dispersion. Similar increases of the flow behavior index were also determined previously [50,51].

The influence of temperature, T, on the viscosity at each shear rate can be expressed in terms of an Arrhenius-type equation [52,53]:

AeEa= RT; 8

where R is the universal gas constant and A and Ea are the fitting parameters of the pre-exponential factor and energy of activation to fluid flow, respectively. This equation describes adequately the temperature dependence of the shear viscosity of the studied nanofluids. Figure 7c shows the obtained Ea values vs. shear rate for the 25 wt.% concentration of A-TiO2/EG and R-TiO2/ EG nanofluids. It is generally accepted that higher Ea values indicate a faster change in viscosity with temperature and high temperature dependency of viscosity [50]. Thus, lower Ea values found for A-TiO2/EG indicate an inferior temperature influence on viscosity for this nanofluid. Moreover, at shear rates around 6 s1

for A-TiO2/EG and around 8 s1 for R-TiO2/EG, a minimum of the energy of activation was detected, as can be observed in Figure 7c. The values obtained here for A-TiO2/EG and R-TiO2/EG are similar to those obtained by Abdelhalim et al. [54] for gold nanoparticles in an aqueous solution.

In addition, linear viscoelastic oscillatory experiments were performed for A-TiO2/EG in order to study their mechanical properties under small-amplitude oscillatory shear. The power of these tests is that stress can be separated into two terms and the elastic or storage modulus can be determined. Then, it can be established whether the nanofluid behaves as the base fluid without agglomerates or alternatively as a solid with a certain level of agglomerates due to the increase in the interactions and collisions among particles that lead to gel

10

a)

1

/ Pas

0.1

0.01

0.1 1 10 100 1000

.

/ s-1

10

b)

/ Pas

1

0.1

0.01

0.1 1 10 100 1000

.

/ s-1

400

c)

300

E A/Jmol-1

200

100

0

1 10 100 1000

.

/ s-1

Figure 7 Viscosity () vs. shear ( _

) rate of EG/TiO2 nanofluids at

different temperatures. Flow curves for (a) A-TiO2/EG and (b)

R-TiO2/EG at 25 wt.% mass concentration and at different tempera

tures: square, 283.15 K; circle, 293.15 K; triangle, 303.15 K; diamond, 313.15 K; cross mark, 323.15 K. (c) Energy of activation to fluid flow (Ea) vs. shear rate for A-TiO2/EG (filled diamond) and R-TiO2/EG

(empty diamond) 25 wt.% nanofluids.

Cabaleiro et al. Nanoscale Research Letters 2013, 8:286 Page 11 of 13 http://www.nanoscalereslett.com/content/8/1/286

formation [55]. First, with the aim to identify the linear viscoelastic region, strain sweep tests (for strains between0.01% and 1,000%) were carried out at 10 rad s1 (see Figure 8a,b). Smaller strain amplitudes were not considered due to equipment conditions as the strain waveform was not sinusoidal due to the presence of experimental noise. A linear regime was found, over which G and G remain constant at low strains with critical strains lower than 1%, which are weakly concentration dependent whereas the stress upper limit of the linear viscoelastic regime region increases with concentration. After this critical strain, G and G decrease as the strain increases in two steps, which may correspond to, first, the break of the structure and then the orientation of agglomerates aligned with the flow field at large deformations [55]. This two-step decrease presents two peaks, which become more evident at higher concentrations, that were previously described in the literature as an attractive gel structure [55,56]. Figure 8c shows the strain dependences of the shear stress for the deformation tests where the strains at these two peaks are identified with arrows. The loss modulus clearly decreases at a strain beyond 1%, and no overshoot trend is observed as found on other nanofluids [32].

Frequency sweep tests (for angular frequencies between0.1 and 600 rad s1) were performed for A-TiO2/EG

nanofluids, and the evolution of each modulus with the

oscillation frequency was obtained, as shown in Figure 8c,d. These experiments were carried out in the linear visco-elastic region using a constant strain value of 0.1% for all nanofluids. Both moduli increase with concentration at a given constant frequency which means that when the nanoparticle content is increased, the hydrodynamic interactions as well as the probability of collision become important, enhancing the aggregation processes. In all cases, the elastic modulus is higher than the viscous one at low frequencies, while the contrary occurs at high frequencies, where the suspensions behave like a liquid. Crossover frequencies, where G = G and a change in the viscoelastic behavior is detected, increase with the concentration of nanoparticles from around 4 rad s1 at a concentration of 10 wt.% to 15 rad s1 at 25 wt.%. That is in agreement with the fact that the degree of agglomeration of the particles is more important at the highest concentrations, but the alignment with the flow of the aggregates is achieved in a shorter time for higher concentrations. This analysis was not carried out for the lowest nanofluid concentration (5 wt.%) due to the availability of the minimum torque of the used device. Moreover, it should be taken into account that those data at elevated frequencies in which problems of inertia of equipment appear were not considered. This

c)

100

0.01 1 100

1

100

/ Pa

100

0.01

b)

10

0.0001

10

%

G' / Pa

1

G" / Pa

1

0.1

a)

0.01

0.1

0.01 0.1 1 10 100 1000

0.01 0.1 1 10 100 1000

%

%

1000

1000

d)

e)

100

100

G' / Pa

10

G" / Pa

10

1

1

0.1

0.1

0.1 1 10 100 1000

0.1 1 10 100 1000

/rads-1

/rads-1

Figure 8 Storage (G) and loss (G) moduli. (a) Storage modulus, (b) loss modulus, and (c) shear stress () as a function of strain () at an angular frequency of 10 rad s1 and a temperature of 303.15 K for different concentrations of A-TiO2/EG. (d) Storage and (e) loss moduli as a

function of frequency () at a strain of 0.1% and a temperature of 303.15 K for different concentrations of A-TiO2/EG. Line, 5 wt.%; circle, 10 wt.%;

square, 15 wt.%; diamond, 20 wt.%; triangle, 25 wt.%.

Cabaleiro et al. Nanoscale Research Letters 2013, 8:286 Page 12 of 13 http://www.nanoscalereslett.com/content/8/1/286

was done by taking into consideration the relationship between the complex viscosity and the frequency. The loss and storage moduli increase with frequency especially at frequencies higher than 10 rad s1. It can be also observed that the elastic modulus data fall on a straight line for the highest frequencies. Finally, we want to point out that the increase in nanoparticle concentration leads to an increase in the formation of agglomeration of the particle, but even the concentration of 5 wt.% for A-TiO2/EG nanofluid does not follow the conventional Cox-Merz rule [57],

j j

_

_

, * being the complex viscosity * (G + iG)/

, which is often valid for Newtonian or non-structured fluids. Our data demonstrate the Cox-Merz rule to become more inapplicable as the concentration of nanoparticles increases. Moreover, it is illustrated that the addition of nanoparticles makes the difference |*| increase as _

; 0 for all A-TiO2/EG concentrations. This behavior was previously found by Haleem and Nott [58] for suspensions of rigid spheres in semi-dilute polymer solutions. These authors attributed this anomalous behavior to the fact that an anisotropic particle microstructure can form at steady shear even in the limit _

0, while it cannot reach it for small-amplitude oscillatory shear. Up to our knowledge, no previous results were published on the Cox-Merz rule of nanofluids, and therefore, more studies exploring other nanofluid types are necessary.

Conclusions

The density values for R-TiO2/EG are higher than those for the A-TiO2/EG nanofluid at the same temperature, pressure, and mass concentration. The enhancement of density in relation to the base fluid is also higher for ru-tile nanofluids, reaching values of 3.8% at the highest concentration. These increments with the concentration are almost temperature and pressure independent. The isobaric thermal expansivity values of A-TiO2/EG and

R-TiO2/EG nanofluids decrease when the pressure rises and increase with temperature as the base fluid does, while we have found that these values for the nanofluids are very similar to or lower than those for the base fluid, achieving decreases up to 2%. The analyzed nanofluids present an expansive volumetric behavior which is more pronounced in A-TiO2/EG. This expansive behavior is also found for other EG-based nanofluids. Contrarily to what was previously published, a shear thinning non-Newtonian behavior was found for both sets of TiO2/EG nanofluids. As the concentration rises, Newtonian plateaus are found at the lowest shear rate and the shear thinning regions can be described using the Ostwald-de Waele model. At the same temperature and concentration conditions, A-TiO2 nanofluids show higher shear dynamic viscosity for all the shear rates. Minima in the energy of activation were found at shear rates around

6 s1 for A-TiO2/EG and 8 s1 for R-TiO2/EG when the shear dynamic viscosity data were fitted for the 25 wt.% concentrations using an Arrhenius-type equation. Finally, viscoelastic oscillatory experiments were performed for A-TiO2. The two-step decrease of the loss modulus present in the deformation tests evidence an attractive gel behavior of the studied nanofluids. Finally, the A-TiO2/EG nanofluid does not follow the conventional Cox-Merz rule. The differences between the viscosities determined in steady shear and the dynamic viscosities from the oscillatory rate are higher when the nanoparticle concentration increases.

Competing interestsThe authors declare that they have no competing interests.

Authors' contributionsDC performed the nanofluid sample characterization and experimental measurements and participated in the redaction and the critical evaluation of experimental results. MJPG contributed with the selection of the optimal experimental setting and type of tests to be performed. CGF participated in the critical evaluation of experimental and theoretical results. MMP analyzed the data and participated in the structuring of the work. LL conceived the study, developed its design, and coordinated the redaction of the manuscript. All authors read and approved the final manuscript.

AcknowledgementsThis work was supported by the Ministerio de Economa y Competitividad (Spain) and the FEDER program through the project ENE2012-32908. The authors also acknowledge the financial support from Fundacin Iberdrola and Universidad de Vigo. DC and LL acknowledge the financial support under the FPU and Ramn y Cajal Program provided by the Ministerio de Educacin, Cultura y Deporte and Ministerio de Ciencia e Innovacin (Spain), respectively.

Received: 15 April 2013 Accepted: 19 May 2013 Published: 13 June 2013

References1. Choi SUS: Nanofluids: from vision to reality through research. J Heat Transfer 2009, 131:033106/1.

2. Huminic G, Huminic A: Application of nanofluids in heat exchangers: a review. Renew Sustain Energy Rev 2012, 16:56255638.

3. Fan X, Chen H, Ding Y, Plucinski PK, Lapkin AA: Potential of nanofluids to further intensify microreactors. Green Chem 2008, 10:670677.

4. Peyghambarzadeh SM, Hashemabadi SH, Hoseini SM, Seifi JM: Experimental study of heat transfer enhancement using water/ethylene glycol based nanofluids as a new coolant for car radiators. Int Commun Heat Mass Transfer 2011, 38:12831290.

5. Mohseni M, Ramezanzadeh B, Yari H, Gudarzi MM: The role of nanotechnology in automotive industries. In New Advances in Vehicular Technology and Automotive Engineering. Edited by Carmo JP, Ribeiro JE. New York: Intech; 2012:354.

6. Teja AS, Beck MP, Yuan Y, Warrier P: The limiting behavior of the thermal conductivity of nanoparticles and nanofluids. J App Phys 2010, 107:114319.

7. Demir H, Dalkilic AS, Kurekci NA, Duangthongsuk W, Wongwises S: Numerical investigation on the single phase forced convection heat transfer characteristics of TiO2 nanofluids in a double-tube counter flow heat exchanger. Int Commun Heat Mass Transfer 2011, 38:218228.

8. Nayak AK, Singh RK, Kulkarni PP: Measurement of volumetric thermal expansion coefficient of various nanofluids. Tech Phys Lett 2010, 36:696698.

9. Nayak AK, Singh RK, Kulkarni PP: Thermal expansion characteristics of Al2O3 nanofluids: more to understand than understood. Appl Phys Lett 2009, 94:094102.

10. Cabaleiro D, Pastoriza-Gallego MJ, Pieiro MM, Lugo L: Characterization and measurements of thermal conductivity, density and rheological

Cabaleiro et al. Nanoscale Research Letters 2013, 8:286 Page 13 of 13 http://www.nanoscalereslett.com/content/8/1/286

properties of zinc oxide nanoparticles dispersed in (ethane-1,2-diol + water) mixture. J Chem Thermodyn 2013, 58:405415.11. Nayak AK, Gartia MR, Vijayan PK: An experimental investigation of single-phase natural circulation behavior in a rectangular loop with Al2O3

nanofluids. Exp Therm Fluid Sci 2008, 33:184189.12. Grassian VH, OShaughnessy PT, Adamcakova-Dodd A, Pettibone JM, Thorne PS: Inhalation exposure study of titanium dioxide nanoparticles with a primary particle size of 2 to 5 nm. Environ Health Perspect 2007, 115:397402.

13. He Y, Jin Y, Chen H, Ding Y, Cang D, Lu H: Heat transfer and flow behaviour of aqueous suspensions of TiO2 nanoparticles (nanofluids) flowing upward through a vertical pipe. Int J Heat Mass Transfer 2007, 50:22722281.

14. Chen H, Ding Y, Tan C: Rheological behaviour of nanofluids. New J Phys

2007, 9:367.

15. Lee G-J, Kim CK, Lee MK, Rhee CK: Characterization of ethylene glycol based TiO2 nanofluid prepared by pulsed wire evaporation (PWE)

method. Rev Adv Mater Sci 2011, 28:126129.16. Grant FA: Properties of rutile (titanium dioxide). Rev Mod Phys 1959, 31:646674.17. Bobbo S, Fedele L, Benetti A, Colla L, Fabrizio M, Pagura C, Barison S: Viscosity of water based SWCNH and TiO2 nanofluids. Exp Therm Fluid Sci 2012, 36:6571.

18. Penkavova V, Tihon J, Wein O: Stability and rheology of dilute TiO2-water nanofluids. Nanoscale Res Lett 2011, 6:273.

19. Reddy MCS, Rao VV, Reddy BCM, Sarada SN, Ramesh L: Thermal conductivity measurements of ethylene glycol water based TiO2 nanofluids. Nanosci Nanotech Let 2012, 4:105109.

20. Setia H, Gupta R, Wanchoo RK: Thermophysical properties of TiO2-Water based nanofluids. AIP Conf Proc 2011, 1393:267268.

21. Wamkam CT, Opoku MK, Hong H, Smith P: Effects of pH on heat transfer nanofluids containing ZrO2 and TiO2 nanoparticles. J Appl Phys 2011, 109:024305.

22. Xie H, Yu W, Chen W: MgO nanofluids: higher thermal conductivity and lower viscosity among ethylene glycol-based nanofluids containing oxide nanoparticles. J Exp Nanosci 2010, 5:463472.

23. Turgut A, Tavman I, Chirtoc M, Schuchmann HP, Sauter C, Tavman S: Thermal conductivity and viscosity measurements of water-based TiO2 nanofluids. Int J Thermophys 2009, 30:12131226.

24. Tseng WJ, Lin K-C: Rheology and colloidal structure of aqueous TiO2 nanoparticle suspensions. Mater Sci Eng, A 2003, 355:186192.

25. Pak BC, Cho YI: Hydrodynamic and heat transfer study of dispersed fluids with submicron metallic oxide particles. Exp Heat Transfer 1998, 11:151170.

26. Hu C, Duo S, Zhang R, Li M, Xiang J, Li W: Nanocrystalline anatase TiO2 prepared via a facile low temperature route. Mater Lett 2010, 64:20402042.

27. Reyes-Coronado D, Rodriguez-Gattorno G, Espinosa-Pesqueira ME, Cab C, de Coss R, Oskam G: Phase-pure TiO2 nanoparticles: anatase brookite and rutile. Nanotechnology 2008, 19:145605.

28. Pastoriza-Gallego MJ, Casanova C, Pramo R, Barbs B, Legido JL, Pieiro MM: A study on stability and thermophysical properties (density and viscosity) of Al2O3 in water nanofluid. J App Phys 2009, 106:064301.

29. Segovia JJ, Fandio O, Lpez ER, Lugo L, Martn MC, Fernndez J: Automated densimetric system: measurements and uncertainties for compressed fluids. J Chem Thermodyn 2009, 41:632638.

30. Cabaleiro D, Pastoriza-Gallego MJ, Pieiro MM, Legido JL, Lugo L: Thermophysical properties of (diphenyl ether + biphenyl) mixtures for their use as heat transfer fluids. J Chem Thermodyn 2012, 50:8088.

31. Pieiro MM, Bessires D, Gacio JM, Saint-Guirons H, Legido JL: Determination of high-pressure liquid density for n-perfluorohexane and n-perfluorononane. Fluid Phase Equilib 2004, 220:125134.

32. Pastoriza-Gallego MJ, Lugo L, Legido JL, Pieiro MM: Rheological non-Newtonian behaviour of ethylene glycol-based Fe2O3 nanofluids.

Nanoscale Res Lett 2011, 6:560.33. Mariano A, Pastoriza-Gallego MJ, Lugo L, Camacho A, Canzonieri S, Pieiro MM: Thermal conductivity, rheological behaviour and density of non-Newtonian ethylene glycol-based SnO2 nanofluids. Fluid Phase Equilib 2013, 337:119124.

34. Fine RA, Millero FJ: Compressibility of water as a function of temperature and pressure. J Chem Phys 1973, 59:55295536.

35. Guignon B, Aparicio C, Sanz PD: Volumetric properties of sunflower and olive oils at temperatures between 15 and 55C under pressures up to 350 MPa. High Pressure Res 2009, 29:3845.

36. Mikhailov GM, Mikhailov VG, Reva LS, Ryabchuk GV: Precision fitting of the temperature dependence of density and prediction of the thermal expansion coefficient of liquids. Russ J Appl Chem 2005, 78:10671072.

37. Diebold U: The surface science of titanium dioxide. Surf Sci Rep 2003,

48:53229.

38. Pastoriza-Gallego MJ, Lugo L, Legido JL, Pieiro MM: Enhancement of thermal conductivity and volumetric behaviour of FexOy nanofluids.

J Appl Phys 2011, 110:014309.39. Pastoriza-Gallego MJ, Lugo L, Cabaleiro D, Legido JL, Pieiro MM: Thermophysical profile of ethylene glycol-based ZnO nanofluids. J Chem Thermodyn 2013. doi:http://dx.doi.org/10.1016/j.jct.2013.07.002

Web End =10.1016/j.jct.2013.07.002 .

40. Ding Y, Alias H, Wen D, Williams RA: Heat transfer of aqueous suspensions of carbon nanotubes (CNT nanofluids). Int J Heat Mass Transfer 2006, 49:240250.

41. Kwak K, Kim C: Viscosity and thermal conductivity of copper oxide nanofluid dispersed in ethylene glycol. Korea-Aust Rheol J 2005, 17:3540.

42. Prasher R, Song D, Wang J, Phelan P: Measurements of nanofluid viscosity and its implications for thermal applications. App Phys Lett 2006, 89:133108.

43. Chen H, Ding Y, Tan C: Rheological behaviour of nanofluids. New J Phys 2007, 9:367.

44. Namburu PK, Kulkarni DP, Misra D, Das DK: Viscosity of copper oxide nanoparticles dispersed in ethylene glycol and water mixture. Exp Therm Fluid Sci 2007, 32:397402.

45. Chen H, Ding Y: Heat transfer and rheological behaviour of nanofluids a review. In Advances in Transport Phenomena. Edited by Wang L. Berlin: Springer; 2009:135177.

46. Haminiuk CWI, MacIel GM, Plata-Oviedo MSV, Quenehenn A, Scheer AP: Study of the rheological parameters of honey using the Mitschka method. Int J Food Eng 2009, 5:13.

47. Lindner A, Bonn D, Meunier J: Viscous fingering in a shear-thinning fluid. Phys Fluids 2000, 12:256261.

48. Santra AK, Sen S, Chakraborty N: Study of heat transfer due to laminar flow of copper-water nanofluid through two isothermally heated parallel plates. Int J Therm Sci 2009, 48:391400.

49. Alberto C, Naranjo TAO, Sierra JD: Plastics Testing and Characterization: Industrial Applications. Cincinnati: Hanser Gardner Publications; 2008.

50. Ah-Hen KS, Vega-Glvez A, Moraga NO, Lemus-Mondaca R: Modelling of rheological behaviour of pulps and pures from fresh and frozen-thawed murta (Ugni molinae Turcz) berries. Int: J Food Eng; 2012.

51. Chin NL, Chan SM, Yusof YA, Chuah TG, Talib RA: Modelling of rheological behaviour of pummelo juice concentrates using master-curve. J Food Eng 2009, 93:134140.

52. Larson RG: The Structure and Rheology of Complex Fluids. New York: Oxford University Press; 1999.

53. Timofeeva EV, Routbort JL, Singh D: Particle shape effects on thermophysical properties of alumina nanofluids. J App Phys 2009, 106:014304.

54. Abdelhalim MAK, Mady MM, Ghannam MM: Rheological and dielectric properties of different gold nanoparticle sizes. Lipids Health Dis 2011, 10:208.

55. Pham KN, Petekidis G, Vlassopoulos D, Egelhaaf SU, Pusey PN, Poon WCK: Yielding of colloidal glasses. Europhys Lett 2006, 75:624630.

56. Tanaka H, Meunier J, Bonn D: Nonergodic states of charged colloidal suspensions: repulsive and attractive glasses and gels. Phys Rev E Stat Nonlin 2004, 69:031404.

57. Cox WP, Merz EH: Correlation of dynamic and steady-flow viscosities. J Polym Sci 1958, 28:619622.

58. Haleem BA, Nott PR: Rheology of particle-loaded semi-dilute polymer solutions. J Rheol 2009, 53:383400.

doi:10.1186/1556-276X-8-286Cite this article as: Cabaleiro et al.: Rheological and volumetric properties of TiO2-ethylene glycol nanofluids. Nanoscale Research Letters 2013 8:286.