Abstract
In this paper, we investigate the effect of adding a stabilizer on the dynamic thermal properties of ZnO nanofluid (containing 5 to 10 nm diameter of ZnO nanocrystals) measured using a 3ω method. Addition of the stabilizer leads to the stabilization of the nanofluid and also substantial reduction of the enhancement of thermal transport compared to that seen in the bare ZnO nanofluid. This also alters the frequency dependence of the thermal transport and the characteristic time scale associated with it. It is suggested that the addition of the stabilizer inhibits the thermodiffusionassisted local aggregation thus leading to substantial reduction of the enhancement of thermal transport properties of the bare nanofluid as proposed in some recent models, and this also alters the characteristic time scales by altering the scale of aggregation.
Keywords:
Heat transport; Nanofluid; ZnO nanoparticles; StabilizerBackground
Nanofluids are dispersions of nanoparticles (typically sizes approximately 5 to 20 nm) in liquid medium. In recent years, they have attracted considerable attention due to enhanced heat transport properties as seen through enhanced thermal conductance [1,2]. In general, heat transport due to conducting metallic or solid inclusions in nonconducting fluids leads to an enhancement. However, in the nanofluids, which have solid inclusions of sizes in the range of few nanometers or few tens of nanometers, the enhancement in thermal conductivity was found to be much larger than that expected from Maxwell’s effective medium theories [3,4].
A number of mechanisms have been proposed that could be responsible for the enhancement of the thermal conductivity. They include the (a) Brownian motion of the nanoparticles [5,6], (b) molecularlevel layering of the liquid at the liquidparticle interface [7], (c) ballistic heat transport in nanoparticles [8], and (d) local clustering of nanoparticles [9,10]. The suggested mechanisms do provide some level of explanation of the enhancement. However, there is no accepted theory/mechanism that can explain all the observations adequately.
Recently reported experimental studies suggest that the formation of local nanoparticle aggregate can play a significant role in the thermal transport in nanofluids [9,10]. In the context of nanofluids containing Fe nanoparticles, it was demonstrated [11] that Fe nanoparticles in the nanofluids can locally assemble into aggregate of micronsize clusters. It was found in CuO nanofluids that large thermal conductivity enhancements are often accompanied by sharp viscosity that increases at low nanoparticle volume fractions, which has been inferred as an indicative of local aggregation effects [12]. The aggregation can be controlled by surface charge, and the critical importance of particle surface charge in nanofluid thermal conductivity has been demonstrated [13].
In this paper, we carry out an investigation on the effect of local aggregation on the thermal transport in nanofluids. This was done in nanofluids containing ZnO nanoparticles with and without stabilizer. The stabilizer can affect local aggregation which in turn can substantially change the enhancement of the thermal conduction in nanofluids. Importantly, we also show that this affects the characteristic frequency scales associated with the dynamical heat transport in such nanofluids. The characteristic frequency, as derived directly from the dynamic thermal property measurement, is a measure of the scale of local aggregation. The reduction of the scale of local aggregation can reduce the magnitude of the thermal transport enhancement, providing a direct link between the two.
The choice of ZnO nanofluid for the investigation originates from the fact that unlike many metallic nanofluids, ZnO nanofluids can be a stable suspension over hours even without added stabilizers. This stability arises due to surface charges on asprepared ZnO nanoparticles [14]. The stability over hours is long enough that it enables us to carry out the thermal measurements. The addition of polyvinylpyrrolidone (PVP) as a stabilizer enhances the stability even further to weeks and even months. Thus, the system chosen is a very suitable system where the measurements can be carried out in nanofluids with and without stabilizers and thus track the changes in thermal parameters in the addition of the stabilizer.
In our earlier work on ZnO nanofluids [15], which is carried out using a dynamic 3ω technique, it has shown that the parameter effusivity (C_{p}κ, C_{p} = heat capacity, κ = thermal conductivity) has a prominent frequency dependence. The measured effusivity shows appreciable enhancement at low frequency, but above a characteristic frequency, the enhancement is significantly reduced and it approaches the parameters of the base liquid. In this paper, we investigate what happens to the enhancement of C_{p}κ as well as its frequency dependence when a stabilizer is added to the system. We find that the presence of stabilizer, which reduces the local aggregation, actually leads to a significant decrease of the C_{p}κ. We also find that the frequency dependence of C_{p}κ in bare ZnO nanofluid gets quantitatively modified when the stabilizer is attached. In addition, we carry out an analysis of the frequency dependence of the temperature oscillation to separate out the contributions of C_{p} and κ components and find that the enhancement in C_{p}κ is primarily due to the enhancement of thermal conductivity κ.
Methods
Nanofluid synthesis
Stable ZnO nanofluid, which is a dispersion of ZnO nanocrystals in ethanol, is prepared by wet chemical method [16]. The nanocrystals of ZnO were synthesized at low temperature (<90°C) in an alkaline medium using Zn acetate. The nanocrystals of ZnO are crystalline with an average size of approximately 10 nm as seen using the transmission electron microscope (TEM). The typical TEM image of a nanoparticle is shown in Figure 1a. Two nanofluids were prepared. One is a pure dispersion of the ZnO nanocrystals in ethanol, and the other is made by adding PVP as a stabilizer. PVP binds to the polar surface of ZnO. The ZnO nanofluid, even without PVP, can be stabilized in scales of hours. The addition of PVP leads to substantial enhancement of the stability of the nanofluid. PVP has been used in the past to make stable metal colloids of Pd [17], Au, and Ag [18]. We have monitored the stability of the nanofluids by measuring the UV–vis spectra of the fluids as a function of time. Both the nanofluids show characteristic absorption around λ ≈ 360 nm, which is the absorption edge for ZnO. For the ZnO nanofluid without PVP, the absorption initially decreases with time as shown in Figure 1b. The decrease is rapid initially and then slows down considerably after 30 min, when the absorption decreases by about 3% in 1 h. This stability is long enough to carry out the thermal measurements over a period of 2 h. The addition of the PVP leads to a very stable nanofluid that is stable over few weeks as can be seen in Figure 1c where there is no perceptible change in the UVvisible absorption even after 2 weeks.
Figure 1. TEM image of ZnO nanocrystal used. (a) Time dependence of UV–vis spectra for ZnO nanofluids, (b) without and (c) with PVP stabilizer.
Thermal measurements using 3ω technique
The thermal measurements were done using a 3ω technique [1921], where we use a platinum film both as a thermometer and a heater. The method, as applied to nanofluids, is explained elsewhere [15]. Here, we provide a small gist for quick reference.
In this method, the Pt film (width of 300 μm, thickness of 50 nm, and length of 5 mm grown on a glass substrate by magnetron sputtering) carrying a current at frequency f is immersed in the liquid in which measurements have to be made [19]. The periodic heating of the film, due to the sinusoidal current, makes the temperature oscillate around the average with an amplitude δT_{2ω} at a frequency 2ω (ω = 2πf). This leads to resistance oscillations of amplitude δT_{2ω} at frequency 2ω around the mean, where δR_{2ω} = αR_{0}δT_{2ω,}α is the temperature coefficient of resistance (TCR) of the heater, and R_{0} is the average resistance of the heater. The resistance oscillation δR_{2ω} at frequency 2ω mixes with the current at frequency ω to produce a potential drop () with a component at 3ω (sum band). The experiment measures the complex voltage with its phase and amplitude, using a phasesensitive detection technique. The thermal properties of the heateronsubstrate (S) and surrounding liquid (L) are given by two parameters Z and the phase φ. These parameters are obtained experimentally from the observed 3ω signal , the area of the heater (A), the power dissipated (P), and the measured TCR (α) of the Pt film using the equation [19]
where the thermal parameter is the effusivity given as ξ ≡ C_{p}κ. L and S refer to the liquid and the substrate, respectively.
The Pt film has a resistance of ≈ 100 Ω and a measured temperature coefficient of resistivity α ≈ 3.5 × 10^{−3}/K. The relative size of the heater width and the thermodiffusion length (D = thermal diffusivity) determines the lowfrequency range of the experiment. In our case for the base liquid ethanol (D ≈ 9 × 10^{−8} m^{2}/s), the working frequency is for the width of the heater used (approximately 300 μm). At highfrequency range, the limit arises due to the low value of the signal. The experiment was carried out in a temperaturecontrolled bath stabilized to better than ±0.01 K which houses a cylindrical copper shell as the sample container. The typical datataking time for a given frequency scan over the full range is 30 min. After each scan, the suspension is shaken in an ultrasonic shaker before the next run begins.
Using relation and , we obtain the ξ_{NF} for the nanofluid given as [19]
In addition to the effusivity ξ_{NF}, we also find the thermal conductivity κ using the frequency dependence of the temperature oscillation δT_{2ω}. The δT_{2ω} for a line heater has a total width of 2b dissipating power P_{L}/unit length and immersed in a liquid [20]:
where K is the integration variable, , refer to the solid (substratecarrying heater) and the liquid, respectively. The value of the interfacial resistance is expressed as R_{interface} ≈ 6.1 × 10^{−7} m^{2} K/W [20]. From Equation 4, it can be shown that the frequency dependence of δT_{2ω} has a logarithmic dependence on f whose slope is given as [21]
We also determine the specific heat C_{p} of the base liquid and the nanofluids using a differential scanning calorimeter, operating in modulation mode (with frequency <10 mHz).
Results and discussions
Change in thermal effusivity in the addition of stabilizer
The representative data on the detected temperature oscillation δT_{2ω} as a function of frequency is shown in Figure 2. It shows the typical δT_{2ω} data for ZnOPVP nanofluids. From this data, we do the analysis of thermal conductivity of respective nanofluids.
Figure 2. Typical temperature oscillation δT_{2ω}as a function of frequency measured in PVPstabilized ZnO nanofluid.
In Figure 3, we show the effusivity ξ_{NF} = C_{p}κ of the base fluid ethanol along with two nanofluids: the bare ZnO nanofluid as well as the ZnO nanofluid with stabilizer PVP. The data for the base liquid ethanol are also shown. The parameters are obtained from Equations 2 and 4 using the measured data. Both the nanofluids have the same volume fraction of 1.5% and have similar average particle size.
Figure 3. Frequency dependence of effusivity of base liquid ethanol, bare ZnO nanofluid, and PVPstabilized ZnO nanofluid.
The enhancement of ξ_{NF} in the nanofluids, at low frequency, compared to that in ethanol is clearly seen. Importantly, it is observed that the enhancement in the bare nanofluid (without stabilizer) is much larger compared with that in the nanofluid with the PVP stabilizer. The results are summarized in Table 1, where we show the enhancement of the effusivity ξ = C_{p}κ as a ratio taken with respect to (wrt) the base fluid as determined from the analysis of the signal. The lowfrequencylimiting values for ξ were used for the parameters in Table 1.
Table 1. Comparison of thermal parameters for nanofluids as measured by two methods
The value of the measured specific heat C_{p} of the base fluid as well as the nanofluids are comparable (C_{p} ≈ 2.5 J/g K). It is thus clear that the enhancement of the effusivity in both the nanofluids is arising primarily due to the enhancement of the thermal conductivity κ.
To make an independent check on the enhancement of the thermal conductivity, we used the measured frequency dependence of the thermal oscillation δT_{2ω}. Equation 4 gives a limiting lowtemperature slope for δT_{2ω} wrt the frequency (log f) that is proportional to κ^{−1}. Using this information, we obtain the relative enhancement of the thermal conductivity wrt the base fluid ethanol. The data for both the nanofluids are shown in Table 1. It can be seen that this also gives nearly the same value for enhancement (within 15% to 20%), which confirms that there is indeed an enhancement in κ in the nanofluids. It is gratifying that the analysis from both the parameters δT_{2ω} and gives similar results.
It can be seen from Table 1 that the enhancement κ for the bare ZnO nanofluid is significantly larger than that seen in the PVPstabilized ZnO nanofluid. This gives us the first important result that there is indeed a significant reduction in the effusivity and thermal conductivity on stabilizing the ZnO nanofluid with stabilizer that inhibits the local aggregation significantly, which in turn leads to its longterm stability. This observation establishes a direct connection between the enhancement of κ and the local aggregate formation.
The frequency dependence of the enhancement and its analysis
The enhancement of the effusivity in nanofluids has a frequency dependence as shown in Figure 3, where the enhancement decreased at higher frequency, and for f > 30 Hz, the values of C_{p}κ for both the nanofluids approach that of the base fluid ethanol. This frequency dependence of the effusivity for bare ZnO nanofluid (without PVP) has been reported elsewhere [15]. It was proposed that the frequency dependence can arise from dynamic local aggregation. In this paper, we explore the proposed hypothesis whether the frequency dependence indeed has a connection to the local aggregation. At low frequency (f ≤ 10 Hz), the enhancement is large, and it reaches a frequencyindependent value.
The decrease in the effusivity at higher frequency in both the nanofluids can be fitted by the lowpass filter relation:
The corner frequency f_{c} and the order of the filter n can be obtained from the fit to the data. For the ZnO nanofluid without PVP, the data can be fitted by the firstorder filter function (n = 1). For fluid with PVP, we got a different higher order value, which is n = 5.
In Figure 4, we show the fit of the data to Equation 5. The data for both the nanofluids are shown.
Figure 4. Lowpass filter response fit for ZnO nanofluids and ZnOPVP nanofluid.
The data are summarized in Table 2. From the table, it can be seen that there is a significant reduction in the corner frequency f_{c} in the addition of the stabilizer.
Table 2. Corner frequency, relaxation time, and estimated length scale of local agglomeration obtained from the data
The thermally driven local aggregation, which would enhance the local thermal transport and hence the value of the thermal conductivity, would lead to solidlike aggregated region in the nanofluids. It is proposed that the response of the type shown in Equation 5 is a manifestation of this local aggregation. The local aggregates respond to an oscillating temperature field δT_{2ω} with a characteristic thermal relaxation time τ_{c}. This will be related to the characteristic length scales of the local aggregate L_{A} through the thermal diffusivity D by the relation τ_{c} ≈ D^{−1}L_{A}^{2}. The relaxation time will determine the corner frequency f_{c} ≈ (4πτ_{c})^{−1} (the extra factor of 2 arises because the temperature oscillation is at frequency 2f). For frequencies larger than 2f_{c}, the temperature oscillation is too fast for the aggregate to respond leading to a decrease in the enhancement of heat transport.
In Table 2, we show the characteristic time τ_{c} as well as the aggregation length L_{A} as derived from the data. We find that the addition of the stabilizer leads to the reduction of the aggregation length L_{A} by 25% to 30%. The corresponding reduction in effusivity or the thermal conductivity is around 40%. This agrees well with the hypothesis that the local aggregation can control the enhancement of the thermal transport as well as the frequency response.
Conclusions
We have investigated the dynamical thermal property (effusivity and thermal conductivity) of ZnO nanofluids containing ZnO nanocrystals with an average diameter of 10 nm with and without PVP stabilizer. This was done to investigate the role of the stabilizer in the enhancement of thermal transport properties of nanofluids. It had been suggested that thermodiffusionassisted ‘solidlike’ local aggregation of the nanoparticles in the nanofluids can be the origin of enhancement of thermal conductivity in nanofluids. The investigations carried out on bare ZnO nanofluids as well as PVPstabilized nanofluids show that addition of a stabilizer, which inhibits diffusionassisted local aggregation due to attached moiety, leads to reduction in the enhancement of thermal parameters that are observed in bare ZnO nanofluids. It has also been shown, from characteristic time scales of the dynamic thermal measurements, that the scale of aggregation gets reduced in the addition of stabilizers. The experimental results provide evidence that the origin of enhancement of thermal conductivity in nanofluids can arise from local aggregation that occurs by thermodiffusion.
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
RKN and AKR jointly did the planning of the experiment, analysis of the data, and writing the manuscript. RKN did the synthesis, characterization, and the measurements. Both authors read and approved the final manuscript.
Acknowledgments
The authors acknowledge the financial support from the Department of Science and Technology, Government of India as a Unit in Nanoscience and Nanotechnology (UNANST). Thanks to Dr. K. Das and Mr. Rajib Nath for their help and useful discussions.
References

Eastman JA, Phillpot SR, Choi SUS, Keblinski P: Thermal transport in nanofluids.
Annual Rev Mater Res 2004, 34:219246. Publisher Full Text

Fan J, Wang LQ: Review of heat conduction in nanofluids.
J Heat Transfer 2011, 133:040801. Publisher Full Text

Maxwell JC: A Treatise on Electricity and Magnetism. Oxford: Oxford University Press; 1873.

Hamilton RL, Crosser OK: Thermal conductivity of heterogeneous two components systems.
Ind Eng Chem Fundam 1962, 1:187191. Publisher Full Text

Prasher R, Bhattacharya P, Phelan PE: Thermal conductivity of nanoscale colloidal solution (nanofluid).

Bhattacharya P, Saha SK, Yadav A, Phelan PE, Prasher RS: Brownian dynamics simulation to determine the effective thermal conductivity of nanofluids.
J Appl Phys 2004, 95:64926494. Publisher Full Text

Yu W, Choi SUS: The role of interfacial layers in the enhanced thermal conductivity of nanofluids: a renovated Hamilton–Crosser model.
J Nanoparticle Res 2004, 6:355361. Publisher Full Text

Keblinski P, Phillpot SR, Choi SUS, Eastman JA: Mechanisms of heat flow in suspensions of nanosized particles (nanofluids).
Int J Heat Mass Tranfer 2002, 45:855863. Publisher Full Text

Timofeeva EV, Gavrilov AN, McCloskey JM, Tolmachev YV, Sprunt S, Lopatina LM, Selinger JV: Thermal conductivity and particle agglomeration in alumina nanofluids: experiment and theory.

Wu C, Cho TJ, Xu J, Lee D, Yang B, Zachariah MR: Effect of nanoparticle clustering on the effective thermal conductivity of concentrated silica colloids.

Hong TK, Yang HS, Choi CJ: Study of the enhanced thermal conductivity of Fe nanofluids.
J Appl Phys 2005, 97:064311. Publisher Full Text

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

Lee D, Kim JW, Kim BG: A new parameter to control heat transport in nanofluids: surface charge state of the particle in suspension.
J Phys Chem B 2006, 110:4323. PubMed Abstract  Publisher Full Text

Ghosh M, Raychaudhuri AK: Ionic environment control of visible photoluminescence from ZnO nanoparticles.
Appl Phys Letts 2008, 93:123113. Publisher Full Text

Neogy RK, Raychaudhuri AK: Frequency dependent enhancement of heat transport in a nanofluid with ZnO nanoparticles.
Nanotechnology 2009, 20:305706. PubMed Abstract  Publisher Full Text

Ghosh M, Raychaudhuri AK: Structural and optical properties of Zn_{1−x}Mg_{x}O nanocrystals obtained by low temperature method.
J Appl Phys 2006, 100:034315. Publisher Full Text

Durap F, Metin O, Aydemir M, Özkar S: New route to synthesis of PVPstabilized palladium(0) nanoclusters and their enhanced catalytic activity in Heck and Suzuki crosscoupling reactions.
Appl Organometal Chem 2009, 23:498503. Publisher Full Text

Shin HS, Yang HJ, Kim SB, Lee MS: Mechanism of growth of colloidal silver nanoparticles stabilized by polyvinyl pyrrolidone in gammairradiated silver nitrate solution.
J Colloid Interface Sci 2004, 274:8994. PubMed Abstract  Publisher Full Text

Chen F, Shulman J, Xue Y, Chu CW, Nolas GS: Thermal conductivity measurement under hydrostatic pressure using the 3ω method.
Rev Sci Instrum 2004, 75:4578. Publisher Full Text

Cahill DG: Thermal conductivity measurement from 30 to 750 K: the 3ω method.
Rev Sci Instrum 1990, 61:802. Publisher Full Text