SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

Open Access Nano Review

Rheological and volumetric properties of TiO2-ethylene glycol nanofluids

David Cabaleiro, María J Pastoriza-Gallego, Carlos Gracia-Fernández, Manuel M Piñeiro and Luis Lugo*

Author Affiliations

Departamento de Física Aplicada, Facultade de Ciencias, Universidade de Vigo, Vigo 36310, Spain

For all author emails, please log on.

Nanoscale Research Letters 2013, 8:286  doi:10.1186/1556-276X-8-286


The electronic version of this article is the complete one and can be found online at: http://www.nanoscalereslett.com/content/8/1/286


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

© 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), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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

Background

Nowadays, 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 enhancement 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 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 titanium 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 s−1 and also that the pressure drop of these nanofluids is very close to that of the base 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 s−1. 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 nanoparticle 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 s−1, confirming a Newtonian behavior and the expected increase of viscosity with nanoparticle concentration. Chen et al. [14] have also found a Newtonian behavior for TiO2/EG nanofluids containing 0.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 nanocrystalline 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 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.

Table 1. Material description

thumbnailFigure 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).

thumbnailFigure 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.

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 × 10−4 g cm−3[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 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. Non-linear 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 properties

The 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 from 1.0627 g cm−3 for pure EG, at 0.1 MPa and 363.15 K, up to 1.1800 g cm−3 for A-TiO2/EG nanofluids and 1.1838 g cm−3 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 nanocrystalline 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.

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

thumbnailFigure 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.

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 (v = 1/ρ), using the following expression [34], were adjusted to the experimental data:

v T , p = v T , p ref · 1 p p ref B T , p , (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 , p ref = v ref T ref , p ref × e a . θ + b θ 2 with θ = T T ref , (2)

where a, b, and vref(Tref,pref) are the adjustable parameters, vref(Tref,pref) being the specific volume at the reference 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 than 2.8 cm3 g−1. The bulk modulus B(T, p) was adjusted as a function of pressure and temperature with the following polynomial:

B T , p = B p ref , T ref + c θ + d θ 2 + e Δp + f Δ p 2 with Δp = p p ref . (3)

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

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 × 10−4 cm3 g−1, 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 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]:

ρ 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 cm−3[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.

thumbnailFigure 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.

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 = 1 v · v T , p T p = v T , p ref T v T , p ref + B T , p T p ·Δp B T , p · B T , p Δp , (5)

κ T = 1 v · v T , p p T = B T , p B T , p p T ·Δp B T , p · B T , p Δ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 the αp values of nanofluids are very similar to or lower than those of EG, achieving decreases up to 1.0% and 1.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 m E , 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 mol−1, 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.

thumbnailFigure 5. Excess molar volumes of nanofluids V m E . 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].

Rheological behavior

As 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 s−1. 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 containing 0.5 to 8.0 wt.% spherical TiO2 nanoparticles. Chen et al. [14] affirmed that a Newtonian behavior is found at a shear rate higher than 0.05 s−1. 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.

thumbnailFigure 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.

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]:

η = γ ˙ n 1 , (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 R-TiO2/EG nanofluids are higher than those for A-TiO2/EG, as shown in Figure 6c. These n values range from 0.27 to 0.72 for A-TiO2/EG and from 0.33 to 0.83 for R-TiO2/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 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 s−1 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].

thumbnailFigure 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 temperatures: 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.

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

η = e E a / R·T , (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 s−1 for A-TiO2/EG and around 8 s−1 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 formation [55]. First, with the aim to identify the linear viscoelastic region, strain sweep tests (for strains between 0.01% and 1,000%) were carried out at 10 rad s−1 (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].

thumbnailFigure 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 s−1 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.%.

Frequency sweep tests (for angular frequencies between 0.1 and 600 rad s−1) 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 viscoelastic 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 s−1 at a concentration of 10 wt.% to 15 rad s−1 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 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 s−1. 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], η * ω η γ ˙ ] γ ˙ = ω , η* being the complex viscosity η* ≡ ( + 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 rutile 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 s−1 for A-TiO2/EG and 8 s−1 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 interests

The authors declare that they have no competing interests.

Authors' contributions

DC 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.

Acknowledgements

This work was supported by the Ministerio de Economía y Competitividad (Spain) and the FEDER program through the project ENE2012-32908. The authors also acknowledge the financial support from Fundación Iberdrola and Universidad de Vigo. DC and LL acknowledge the financial support under the FPU and Ramón y Cajal Program provided by the ‘Ministerio de Educación, Cultura y Deporte’ and ‘Ministerio de Ciencia e Innovación’ (Spain), respectively.

References

  1. Choi SUS: Nanofluids: from vision to reality through research.

    J Heat Transfer 2009, 131:033106/1. OpenURL

  2. Huminic G, Huminic A: Application of nanofluids in heat exchangers: a review.

    Renew Sustain Energy Rev 2012, 16:5625-5638. Publisher Full Text OpenURL

  3. Fan X, Chen H, Ding Y, Plucinski PK, Lapkin AA: Potential of nanofluids’ to further intensify microreactors.

    Green Chem 2008, 10:670-677. Publisher Full Text OpenURL

  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:1283-1290. Publisher Full Text OpenURL

  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:3-54. OpenURL

  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. Publisher Full Text OpenURL

  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:218-228. Publisher Full Text OpenURL

  8. Nayak AK, Singh RK, Kulkarni PP: Measurement of volumetric thermal expansion coefficient of various nanofluids.

    Tech Phys Lett 2010, 36:696-698. Publisher Full Text OpenURL

  9. Nayak AK, Singh RK, Kulkarni PP: Thermal expansion characteristics of Al2O3 nanofluids: more to understand than understood.

    Appl Phys Lett 2009, 94:094102. Publisher Full Text OpenURL

  10. Cabaleiro D, Pastoriza-Gallego MJ, Piñeiro MM, Lugo L: Characterization and measurements of thermal conductivity, density and rheological properties of zinc oxide nanoparticles dispersed in (ethane-1,2-diol + water) mixture.

    J Chem Thermodyn 2013, 58:405-415. OpenURL

  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:184-189. Publisher Full Text OpenURL

  12. Grassian VH, O’Shaughnessy 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:397-402. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  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:2272-2281. Publisher Full Text OpenURL

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

    New J Phys 2007, 9:367. Publisher Full Text OpenURL

  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:126-129. OpenURL

  16. Grant FA: Properties of rutile (titanium dioxide).

    Rev Mod Phys 1959, 31:646-674. Publisher Full Text OpenURL

  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:65-71. OpenURL

  18. Penkavova V, Tihon J, Wein O: Stability and rheology of dilute TiO2-water nanofluids.

    Nanoscale Res Lett 2011, 6:273. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  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:105-109. Publisher Full Text OpenURL

  20. Setia H, Gupta R, Wanchoo RK: Thermophysical properties of TiO2-Water based nanofluids.

    AIP Conf Proc 2011, 1393:267-268. OpenURL

  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. Publisher Full Text OpenURL

  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:463-472. Publisher Full Text OpenURL

  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:1213-1226. Publisher Full Text OpenURL

  24. Tseng WJ, Lin K-C: Rheology and colloidal structure of aqueous TiO2 nanoparticle suspensions.

    Mater Sci Eng, A 2003, 355:186-192. Publisher Full Text OpenURL

  25. Pak BC, Cho YI: Hydrodynamic and heat transfer study of dispersed fluids with submicron metallic oxide particles.

    Exp Heat Transfer 1998, 11:151-170. Publisher Full Text OpenURL

  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:2040-2042. Publisher Full Text OpenURL

  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. PubMed Abstract | Publisher Full Text OpenURL

  28. Pastoriza-Gallego MJ, Casanova C, Páramo R, Barbés B, Legido JL, Piñeiro MM: A study on stability and thermophysical properties (density and viscosity) of Al2O3 in water nanofluid.

    J App Phys 2009, 106:064301. Publisher Full Text OpenURL

  29. Segovia JJ, Fandiño O, López ER, Lugo L, Martín MC, Fernández J: Automated densimetric system: measurements and uncertainties for compressed fluids.

    J Chem Thermodyn 2009, 41:632-638. Publisher Full Text OpenURL

  30. Cabaleiro D, Pastoriza-Gallego MJ, Piñeiro MM, Legido JL, Lugo L: Thermophysical properties of (diphenyl ether + biphenyl) mixtures for their use as heat transfer fluids.

    J Chem Thermodyn 2012, 50:80-88. OpenURL

  31. Piñeiro MM, Bessières 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:125-134. Publisher Full Text OpenURL

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

    Nanoscale Res Lett 2011, 6:560. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  33. Mariano A, Pastoriza-Gallego MJ, Lugo L, Camacho A, Canzonieri S, Piñeiro MM: Thermal conductivity, rheological behaviour and density of non-Newtonian ethylene glycol-based SnO2 nanofluids.

    Fluid Phase Equilib 2013, 337:119-124. OpenURL

  34. Fine RA, Millero FJ: Compressibility of water as a function of temperature and pressure.

    J Chem Phys 1973, 59:5529-5536. Publisher Full Text OpenURL

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

    High Pressure Res 2009, 29:38-45. Publisher Full Text OpenURL

  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:1067-1072. Publisher Full Text OpenURL

  37. Diebold U: The surface science of titanium dioxide.

    Surf Sci Rep 2003, 48:53-229. Publisher Full Text OpenURL

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

    J Appl Phys 2011, 110:014309. Publisher Full Text OpenURL

  39. Pastoriza-Gallego MJ, Lugo L, Cabaleiro D, Legido JL, Piñeiro MM: Thermophysical profile of ethylene glycol-based ZnO nanofluids.

    J Chem Thermodyn 2013. Publisher Full Text OpenURL

  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:240-250. Publisher Full Text OpenURL

  41. Kwak K, Kim C: Viscosity and thermal conductivity of copper oxide nanofluid dispersed in ethylene glycol.

    Korea-Aust Rheol J 2005, 17:35-40. OpenURL

  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. Publisher Full Text OpenURL

  43. Chen H, Ding Y, Tan C: Rheological behaviour of nanofluids.

    New J Phys 2007, 9:367. Publisher Full Text OpenURL

  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:397-402. Publisher Full Text OpenURL

  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:135-177. OpenURL

  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. OpenURL

  47. Lindner A, Bonn D, Meunier J: Viscous fingering in a shear-thinning fluid.

    Phys Fluids 2000, 12:256-261. Publisher Full Text OpenURL

  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:391-400. Publisher Full Text OpenURL

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

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

  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:134-140. Publisher Full Text OpenURL

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

  53. Timofeeva EV, Routbort JL, Singh D: Particle shape effects on thermophysical properties of alumina nanofluids.

    J App Phys 2009, 106:014304. Publisher Full Text OpenURL

  54. Abdelhalim MAK, Mady MM, Ghannam MM: Rheological and dielectric properties of different gold nanoparticle sizes.

    Lipids Health Dis 2011, 10:208. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

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

    Europhys Lett 2006, 75:624-630. Publisher Full Text OpenURL

  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. OpenURL

  57. Cox WP, Merz EH: Correlation of dynamic and steady-flow viscosities.

    J Polym Sci 1958, 28:619-622. Publisher Full Text OpenURL

  58. Haleem BA, Nott PR: Rheology of particle-loaded semi-dilute polymer solutions.

    J Rheol 2009, 53:383-400. Publisher Full Text OpenURL