SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

Open Access Nano Express

Rheological non-Newtonian behaviour of ethylene glycol-based Fe2O3 nanofluids

María Jose Pastoriza-Gallego, Luis Lugo, José Luis Legido and Manuel M Piñeiro*

Author affiliations

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

For all author emails, please log on.

Citation and License

Nanoscale Research Letters 2011, 6:560  doi:10.1186/1556-276X-6-560

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


Received:4 July 2011
Accepted:25 October 2011
Published:25 October 2011

© 2011 Pastoriza-Gallego 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

The rheological behaviour of ethylene glycol-based nanofluids containing hexagonal scalenohedral-shaped α-Fe2O3 (hematite) nanoparticles at 303.15 K and particle weight concentrations up to 25% has been carried out using a cone-plate Physica MCR rheometer. The tests performed show that the studied nanofluids present non-Newtonian shear-thinning behaviour. In addition, the viscosity at a given shear rate is time dependent, i.e. the fluid is thixotropic. Finally, using strain sweep and frequency sweep tests, the storage modulus G', loss modulus G″ and damping factor were determined as a function of the frequency showing viscoelastic behaviour for all samples.

Keywords:
nanofluid; rheology; hematite; ferrofluids; storage modulus; loss modulus; nanoparticles; shear thinning

Introduction

Research on nanofluids characterization has progressed remarkably in the last decade [1-8]. The first studies were performed at the US Argonne National Laboratories reporting anomalous thermal conductivity enhancements, beyond the prediction of classic models. Nowadays, it is well known that this effect depends on particle size, concentration, nature of base fluids, pH, temperature and nanoparticles clustering [5,7]. Moreover, it has been shown that other transport properties exhibit also unusual behaviour, including viscosity and rheological properties [9-12]. This means that the thermophysical profile of a nanofluid may be tuned to meet the requirements for a given industrial application. From a practical point of view, this offers numerous benefits [8,13,14] as improved heat transfer and stability, microchannel cooling without clogging, or reduction in required pumping power. Thus, nanofluids have emerged as suitable tailored working fluids in industrial, engineering and medical applications [5,13-15], but this requires a rigorous analysis of heat transfer and rheological properties. The effective viscosity of a nanofluid constitutes a key property as it governs the ease of flow, pressure drop and thus the pumping power involved during flow applications [16]. Concerning rheological behaviour of nanofluids, only a reduced number of studies can be found (see, e.g. Prasher et al. [17], Kwak and Kim [18], Chen et al. [2,9,19], Rao [20], Namburu et al. [21], Chevalier et al. [11]), evidencing a gap where further studies concerning rigorous characterization of their Newtonian behaviour limits and their viscoelastic trend are necessary. Moreover, recent studies have identified nanoparticle structuring/aggregation as a dominant mechanism for the thermal conductivity enhancement of nanofluids, and rheological analysis can provide a useful insight on their structure [2].

Following our previous research on nanofluids [22-25], we present experimental evidence of non-Newtonian behaviour of nanofluids obtained by dispersing hematite (Fe2O3) nanoparticles in ethylene glycol (EG). These ferrofluids are termed as smart functional fluids, due to some of its unique features, manifesting simultaneously fluid and magnetic properties, and have found applications in mechanical engineering, aerospace and bioengineering [26,27]. The selected base fluid, EG, constitutes an excellent benchmark to compare viscosity results with literature. As an example of rheological analysis of ferrofluids, Hong et al. [26] studied water-based Fe3O4 nanofluids, reporting shear-thinning behaviour.

Experimental

In this work, homogeneous and stable suspensions of commercial hexagonal scalenohedral-shaped α-Fe2O3 (hematite) nanoparticles in EG were prepared at concentrations up to 25% in mass fraction (6.6% in volume fraction). The average nanoparticle diameter value determined was 29 ± 18 nm. More details about nanofluid preparation, stability and characterization have been recently reported [25]. These nanofluids were subjected to rheological analyses using a Physica MCR 101 rheometer (Anton Paar, Graz, Austria). The equipment allows to control torques between 0.5 μN·m and 125 mN·m and normal force from 0.1 to 30 N. The cone-plate geometry with a cone diameter of 25 mm and a cone angle of 1° was used. All experiments are conducted at a constant gap of 0.048 mm, and an initial stabilization period of 100 s is given for achieving constant temperature (303.15 K) using a Peltier system. Three replicates at each experimental condition were carried out.

Results and discussion

With the aim to check the operation of this rheometer using a cone-plate geometry and at shear rates up to 1,000 s-1 in the flow curves, initial experiments based on flow curves at controlled shear stress were carried out for pure EG, diisodecyl phthalate (DiDP) and polyalpha olefin (PAO-40). DiDP and PAO-40 represent Newtonian reference materials [28] in the moderate- to high-viscosity region. If compared with literature [10,25,29-32], excellent agreement is obtained for viscosities, with average deviations of 1.5%, 1.1% and 0.8% for EG [10,25,32], DiDP [30,31] and PAO-40 [29], respectively.

The rheological studies were performed under two types of flow [33,34]. The first is a non-linear viscoelastic experiment, the flow curve, or measurement of shear viscosity (η) as a function of shear rate (<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/560/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/560/mathml/M1">View MathML</a>). The second is the linear viscoelastic oscillatory experiment, leading to the determination of frequency-dependent energy storage modulus G' (elastic) and loss modulus G″ (viscous), which reveal the mechanical properties of the material under small amplitude oscillatory shear. Oscillatory shear measurements within the linear viscoelastic domain, intended to measure G' and G″, represent a useful way of characterizing complex fluids.

Figure 1 shows the viscosity of EG vs. <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/560/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/560/mathml/M1">View MathML</a> at 303.15 K, obtained from controlled shear stress tests. The applied torques start from 0.1 μN·m, covering a wide range of <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/560/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/560/mathml/M1">View MathML</a> (3 to 1,000 s-1). Shear viscosity is independent of <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/560/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/560/mathml/M1">View MathML</a>, indicating Newtonian behaviour for EG. The flow curves of the EG/Fe2O3 nanofluids are also plotted, showing shear-thinning (pseudoplastic) non-Newtonian behaviour. As concentration rises, a Newtonian plateau with shear thinning appears in the lowest <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/560/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/560/mathml/M1">View MathML</a> region, except for the lower concentration, due to the minimum <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/560/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/560/mathml/M1">View MathML</a> threshold value of the rheometer used. The shear-thinning behaviour has also been observed by Hong et al. [26] for water/Fe3O4 ferrofluids, He et al. [35] and Tseng and Lin [36] for water/TiO2 nanofluids or Chen et al. [9] for EG/titanate nanotube nanofluids. However, this behaviour is completely different to the observed Newtonian one of EG/TiO2 [10] or propylene glycol/Al2O3 [17] nanofluids. While in some papers Newtonian behaviour has been supposed a priori, these results evidence that rheological tests must be always carried out to support such an affirmation. Most commercial viscometers work at fixed <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/560/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/560/mathml/M1">View MathML</a> (≈100 s-1), and in these cases, the existence of unnoticed shear thinning may lead to a trend in the experimental viscosity measurements that might be erroneously attributed to the appearance of an anomalous enhancement.

thumbnailFigure 1. Viscosity vs. shear rate (<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/560/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/560/mathml/M1">View MathML</a>) dependence of EG/Fe2O3 nanofluids. At 303.15 K and t = 500 s for different weight concentrations: multiplication sign, EG; solid circle, 5 wt.%; solid inverted triangle, 10 wt.%; solid square, 15 wt.%, solid diamond, 20 wt.% and solid triangle, 25 wt.%. The inset shows viscosity vs. time dependence for 25 wt.% nanofluids at controlled shear stress.

Shear thinning of well-dispersed suspensions can be linked to the modifications in the structure and arrangement of interacting particles [37]. Shearing may cause the particles to orient in the direction of flow and its gradient. This can break agglomerates and hence reduce the amount of solvent immobilized by the particles. The interaction forces may then decrease and cause a lowering in the flow resistance and the apparent viscosity of the system.

The inset in Figure 1 shows the time evolution of shear viscosity for the 25 wt.% EG/Fe2O3 sample, and its decrease evidences thixotropic behaviour or a structure loss under shear. For this reason, all flow curves were measured after the preliminary application of a constant stress during 500 s. This time evolution of viscosity had not been reported for nanofluids so far, but it must be considered when performing viscosity measurements of nanofluids because it may also produce spurious trends for the measured data.

The following step was performing oscillatory or dynamic experiments to determine the viscoelastic behaviour. The power of dynamic testing is that stress can be separated into its elastic and viscous contributions, and the elastic or storage modulus G' and the viscous or loss modulus G″ can be calculated. First, strain sweep tests at constant ω = 10 rad s-1 were carried out (cf. Figure 2) to identify the linear viscoelastic region in the strain range from 0.1% to 1,000%. The linear regime, where G' and G″ are constant regardless of strain amplitude, is clearly observed. G' decreases monotonically as strain increases (Figure 2, upper panel), while G″ goes through a maximum, exhibiting an overshoot phenomenon (Figure 2, lower panel). The interpretation is that when an external strain is imposed, the structure of nanofluids resists the deformation up to a certain strain, where G″ increases, and then the structure is lost by the disaggregation of nanoparticles due to large deformations over the critical strain, after which the nanoparticles align with the flow field, decreasing G' and G″. This effect, more important at higher concentrations [25,38,39], may be governed by the aggregates dimension and radius of gyration. In order to evidence the well-defined linear viscoelastic range of this test, carried out at 10 rad s-1, the stress-strain curves are also displayed in Figure 3, for all concentrations. The critical strain is shown to be independent of concentration, while the stress upper limit of the linear viscoelastic regime region increases linearly with concentration, as shown in the inset of the figure.

thumbnailFigure 2. Storage (G') and loss (G″) moduli. As a function of strain at 10 rad/s and 303.15 K for different concentrations: (solid circle, empty circle) 5 wt.%, (solid inverted triangle, empty inverted triangle) 10 wt.%, (solid square, empty square) 15 wt.%, (solid diamond, empty diamond) 20 wt.% and (solid triangle, empty triangle) 25 wt.%.

thumbnailFigure 3. Shear stress (σ) as function of strain (γ) at 10 rad s-1 for different concentrations. Solid circle, 5 wt.%; solid inverted triangle, 10 wt.%; solid square, 15 wt.%; solid diamond, 20 wt.% and solid triangle, 25 wt.%. The inset shows yield stress as a function of the nanoparticles concentrations.

Then, frequency sweep tests were carried out in the linear viscoelastic region, with angular frequencies ranging from 0.1 to 400 rad s-1, with a constant strain value of 1%. The experimental data of storage (elastic) and loss (viscous) moduli are shown in Figure 4. The storage modulus exceeds the loss modulus, G' > G″, especially for higher concentrations, and G' values are practically constant in the low frequency range, indicating a typical gel structure and the dominant elastic nature of the material under these conditions. However, for 5 wt.% concentration, a crossover frequency (G' = G″) appears, meaning that below 10 rad s-1, it shows liquid-like behaviour but changes beyond that value. Both moduli increase with concentration at a given constant frequency, they increase with frequency beyond an approximate value of 10 rad s-1, and almost all G″ values fall on a straight line for the highest frequencies. These results must be underlined because, starting from a base fluid that exhibits Newtonian behaviour, the addition of hematite nanoparticles produces, even at low concentrations and frequencies, a continuous transition towards elastic behaviour, which means that rheological studies of its viscoelastic nature become essential to determine its potential practical use for any technical application. Moreover, due to their magnetic nature, the influence of an external magnetic field in the viscoelastic behaviour of this nanofluid becomes an attractive topic.

thumbnailFigure 4. Storage (G') and loss (G″) moduli. As a function of frequency for EG/Fe2O3 nanofluids at 1% strain: (solid circle, empty circle) 5 wt.%; (solid inverted triangle, empty inverted triangle) 10 wt.%; (solid square, empty square) 15 wt.%, (solid diamond, empty diamond) 20 wt.% and (solid triangle, empty triangle) 25 wt.%.

Finally, the trend of G″/G' (tangent of the phase angle δ, denoted as damping factor) with frequency was determined (Figure 5). For metals, the damping factor is typically small (approximately 0.0005), whereas for viscoelastic materials, it may exceed the unity. This damping factor is shown to decrease as concentration rises, corresponding to a G' increase, or an evolution from fluid towards elastic behaviour, as stated. Another remarkable feature is the presence of a well-defined maximum in the damping factor, appearing at frequencies shifting to higher values with sample concentration, and its height decreases exponentially with concentration.

thumbnailFigure 5. Damping factor (tg δ). As a function of frequency at 1% strain for various concentrations of EG/Fe2O3 nanofluids: empty circle, 5 wt.%; empty inverted triangle, 10 wt.%; empty square, 15 wt.%; empty diamond, 20 wt.%; and empty triangle, 25 wt.%. The inset shows maximum damping factor as a function of the weight concentration of nanofluids.

Conclusions

This work evidences the non-Newtonian nature of EG/Fe2O3 nanofluids, showing shear thinning and thixotropy. All samples show viscoelastic nature, suggesting that a combination of particle aggregation and shape effects is the mechanism for its high-shear rheological behaviour, which is also supported by the thermal conductivity measurements [7,24,39]. G' decreases after a certain critical strain, and G″ presents an overshoot phenomenon. Finally, the results of the frequency sweep show that the damping factor presents a maximum against frequency, corresponding to a continuous evolution with concentration from viscous to elastic nature. This is an evidence of important aggregation and structural changes in the samples, a subject still poorly studied that deserves further attention.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

MJPG performed the nanofluid samples characterization and experimental measurements. LL contributed with the selection of the optimal experimental setting of the rheometer and type of tests to be performed, and coordinated the redaction of the manuscript. JLL participated in the critical evaluation of experimental results. MMP conceived the study, and participated in its design and coordination. All authors read and approved the final manuscript.

Acknowledgements

The authors acknowledge CACTI (Univ. de Vigo) for technical assistance, and Univ. de Vigo, Xunta de Galicia (grant PGIDIT07PXIB314181PR), Min. de Educación y Ciencia (grant CTQ2006-15537-C02/PPQ) and Min. de Ciencia e Innovación (Ramón y Cajal Program), all in Spain, for financial support.

References

  1. Buongiorno J, et al.: A benchmark study on the thermal conductivity of nanofluids.

    J Appl Phys 2009, 106:094312. Publisher Full Text OpenURL

  2. Chen H, Ding Y: Heat transfer and rheological behaviour of nanofluids - a review. In Advances in Transport Phenomena. Volume 1. Edited by Wang L. Berlin: Springer; 2009:135. OpenURL

  3. Das SK, Choi SUS, Yu W, Pradeep T: Nanofluids: Science and Technology. New York: Wiley; 2008. OpenURL

  4. Li Y, Zhou J, Tung S, Schneider E, Xi S: A review on development of nanofluid preparation and characterization.

    Powder Technology 2009, 196:89. Publisher Full Text OpenURL

  5. Murshed SMS, Leong KC, Yang C: Thermophysical and electrokinetic properties of nanofluids - a critical review.

    App Thermal Eng 2008, 28:2109. Publisher Full Text OpenURL

  6. Philip J, Shima PD, Raj B: Nanofluid with tunable thermal properties.

    Appl Phys Lett 2008, 92:043108. OpenURL

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

  8. Wang XQ, Mujumdar AS: Heat transfer characteristics of nanofluids: a review.

    Int J Therm Sci 2007, 46:1. Publisher Full Text OpenURL

  9. Chen H, Ding Y, Lapkin A, Fan X: Rheological behavior of ethylene glycol-titanate nanotube nanofluids.

    J Nanopart Res 2009, 11:1513. Publisher Full Text OpenURL

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

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

  11. Chevalier J, Tillement O, Ayela F: Rheological properties of nanofluids flowing through microchannels.

    Appl Phys Lett 2007, 91:233103. OpenURL

  12. Heine DR, Petersen MK, Grest GS: Effect of particle shape and charge on bulk rheology of nanoparticle suspensions.

    J Chem Phys 2010, 132:184509. OpenURL

  13. Choi SUS, Zhang ZG, Keblinski P: Nanofluids.

    Encyclopedia of Nanoscience and Nanotechnology 2004, 6:757-773. OpenURL

  14. Zussman S: New Nanofluids Increase Heat Transfer Capability. USA: Argonne National Laboratory; 1997:4. OpenURL

  15. Wang L, Fan J: Nanofluids research: key issues.

    Nanoscale Res Lett 2010, 5:1241. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. Anoop KB, Kabelac S, Sundararajan T, Das SK: Rheological and flow characteristics of nanofluids: Influence of electroviscous effects and particle agglomeration.

    J Appl Phys 2009, 106:034909. OpenURL

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

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

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

  19. Chen H, Ding Y, Lapkin A: Rheological behaviour of nanofluids containing tube/rod-like nanoparticles.

    Powder Technol 2009, 194:132. Publisher Full Text OpenURL

  20. Rao Y: Nanofluids: stability, phase diagram, rheology and applications.

    Particuology 2010, 8:549. Publisher Full Text OpenURL

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

  22. Pastoriza-Gallego MJ, Casanova C, Legido JL, Piñeiro MM: CuO in water nanofluid: influence of particle size and polydispersity on volumetric behaviour and viscosity.

    Fluid Phase Equilib 2011, 300:188. Publisher Full Text OpenURL

  23. 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 Al2 O3 in water nanofluid.

    J App Physics 2009, 106:064301. OpenURL

  24. Pastoriza-Gallego MJ, Lugo L, Legido JL, Piñeiro MM: Thermal conductivity and viscosity measurements of ethylene glycol-based Al2O3 nanofluids.

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

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

  26. Hong RY, Ren ZQ, Han YP, Li HZ, Zheng Y, Ding J: Rheological properties of water-based Fe3O4 ferrofluids.

    Chem Eng Sci 2007, 62:5912. Publisher Full Text OpenURL

  27. Odenbach S: Ferrofluids: Magnetically Controllable Fluids and Their Applications. Berlin: Springer-Verlag; 2002. OpenURL

  28. Caetano FJP, Fareleira JMNA, Fröba AP, Harris KR, Leipertz A, Oliveira CMBP, Trusler JPM, Wakeham WA: An industrial reference fluid for moderately high viscosity.

    J Chem Eng Data 2008, 53:2003. Publisher Full Text OpenURL

  29. Bair S: Pressure-viscosity behavior of lubricants to 1.4 GPa and its relation to EHD traction.

    Tribol Transactions 2000, 43:91. Publisher Full Text OpenURL

  30. Harris KR: Temperature and pressure dependence of the viscosities of 2-ethylhexyl benzoate, bis(2-ethylhexyl) phthalate, 2,6,10,15,19,23-hexamethyltetracosane (squalane), and diisodecyl phthalate.

    J Chem Eng Data 2009, 54:2729. Publisher Full Text OpenURL

  31. Paredes X, Fandiño O, Comuñas MJP, Pensado AS, Fernández J: Study of the effects of pressure on the viscosity and density of diisodecyl phthalate.

    J Chem Thermodyn 2009, 41:1007. Publisher Full Text OpenURL

  32. Sun T, Teja AS: Density, viscosity, and thermal conductivity of aqueous ethylene, diethylene, and triethylene glycol mixtures between 290 K and 450 K.

    J Chem Eng Data 2003, 48:198. Publisher Full Text OpenURL

  33. Tadros TF: Rheology of dispersions: principles and applications. Weinheim: Wiley-VCH; 2010.

  34. Hyun K, Kim SH, Ahn KH, Lee SJ: Large amplitude oscillatory shear as a way to classify the complex fluids.

    J Non-Newtonian Fluid Mech 2002, 107:51. Publisher Full Text OpenURL

  35. 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 and Mass Transfer 2007, 50:2272. Publisher Full Text OpenURL

  36. Tseng WJ, Lin KC: Rheology and colloidal structure of aqueous TiO2 nanoparticle suspensions.

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

  37. Larson RG: The rheology of dilute solutions of flexible polymers: progress and problems.

    J Rheol 2005, 49:1. Publisher Full Text OpenURL

  38. Gharagozloo PE, Goodson KE: Aggregate fractal dimensions and thermal conduction in nanofluids.

    J Appl Phys 2010, 108:074309/1. OpenURL

  39. Prasher R, Evans W, Meakin P, Fish J, Phelan P, Keblinski P: Effect of aggregation on thermal conduction in colloidal nanofluids.

    Applied Phys Lett 2006, 89:143119. Publisher Full Text OpenURL