SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

This article is part of the series Nanofluids.

Open Access Open Badges Nano Express

Heat transfer augmentation in nanofluids via nanofins

Peter Vadasz

Author Affiliations

Department of Mechanical Engineering, Northern Arizona University, P. O. Box 15600, Flagstaff, AZ 86011-5600, USA

Faculty of Engineering, University of KZ Natal, Durban 4041, South Africa

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

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

Received:11 September 2010
Accepted:18 February 2011
Published:18 February 2011

© 2011 Vadasz; 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.


Theoretical results derived in this article are combined with experimental data to conclude that, while there is no improvement in the effective thermal conductivity of nanofluids beyond the Maxwell's effective medium theory (J.C. Maxwell, Treatise on Electricity and Magnetism, 1891), there is substantial heat transfer augmentation via nanofins. The latter are formed as attachments on the hot wire surface by yet an unknown mechanism, which could be related to electrophoresis, but there is no conclusive evidence yet to prove this proposed mechanism.


The impressive heat transfer enhancement revealed experimentally in nanofluid suspensions by Eastman et al. [1], Lee et al. [2], and Choi et al. [3] conflicts apparently with Maxwell's [4] classical theory of estimating the effective thermal conductivity of suspensions, including higher-order corrections and other than spherical particle geometries developed by Hamilton and Crosser [5], Jeffrey [6], Davis [7], Lu and Lin [8], Bonnecaze and Brady [9,10]. Further attempts for independent confirmation of the experimental results showed conflicting outcomes with some experiments, such as Das et al. [11] and Li and Peterson [12], confirming at least partially the results presented by Eastman et al. [1], Lee et al. [2], and Choi et al. [3], while others, such as Buongiorno and Venerus [13], Buongiorno et al. [14], show in contrast results that are in agreement with Maxwell's [4] effective medium theory. All these experiments were performed using the Transient-Hot-Wire (THW) experimental method. On the other hand, most experimental results that used optical methods, such as the "optical beam deflection" [15], "all-optical thermal lensing method" [16], and "forced Rayleigh scattering" [17] did not reveal any thermal conductivity enhancement beyond what is predicted by the effective medium theory. A variety of possible reasons for the excessive values of the effective thermal conductivity obtained in some experiments have been investigated, but only few succeeded to show a viable explanation. Jang and Choi [18] and Prasher et al. [19] show that convection due to Brownian motion may explain the enhancement of the effective thermal conductivity. However, if indeed this is the case then it is difficult to explain why this enhancement of the effective thermal conductivity is selective and is not obtained in all the nanofluid experiments. Alternatively, Vadasz et al. [20] showed that hyperbolic heat conduction also provides a viable explanation for the latter, although their further research and comparison with later-published experimental data presented by Vadasz and Govender [21] led them to discard this possibility.

Vadasz [22] derived theoretically a model for the heat conduction mechanisms of nanofluid suspensions including the effect of the surface area-to-volume ratio of the suspended nanoparticles/nanotubes on the heat transfer. The theoretical model was shown to provide a viable explanation for the excessive values of the effective thermal conductivity obtained experimentally [1-3]. The explanation is based on the fact that the THW experimental method used in all the nanofluid suspensions experiments listed above needs a major correction factor when applied to non-homogeneous systems. This time-dependent correction factor is of the same order of magnitude as the claimed enhancement of the effective thermal conductivity. However, no direct comparison to experiments was possible because the authors [1-3] did not report so far their temperature readings as a function of time, the base upon which the effective thermal conductivity is being evaluated. Nevertheless, in their article, Liu et al. [23] reveal three important new results that allow the comparison of Vadasz's [22] theoretical model with experiments. The first important new result presented by Liu et al. [23] is reflected in the fact that the value of "effective thermal conductivity" revealed experimentally using the THW method is time dependent. The second new result is that those authors present graphically their time-dependent "effective thermal conductivity" for three specimens and therefore allow the comparison of their results with the theoretical predictions of this study showing a very good fit as presented in this article. The third new result is that their time dependent "effective thermal conductivity" converges at steady state to values that according to our calculations confirm the validity of the classical Maxwell's theory [4] and its extensions [5-10].

The objective of this article is to provide an explanation that settles the conflict between the apparent enhancement of the effective thermal conductivity in some experiments and the lack of enhancement in other experiments. It is demonstrated that the transient heat conduction process in nanofluid suspensions produces results that fit well with the experimental data [23] and validates Maxwell's [4] method of estimating the effective thermal conductivity of suspensions. The theoretical results derived in this article are combined with experimental data [23] to conclude that, while there is no improvement in the effective thermal conductivity of nanofluids beyond the Maxwell's effective medium theory [4], there is nevertheless substantial heat transfer augmentation via nanofins. The latter are formed as attachments on the hot wire surface by a mechanism that could be related to electrophoresis and therefore such attachments depend on the electrical current passing through the wire, and varies therefore amongst different experiments. Also since the effective thermal conductivity does not increase beyond the Maxwell's [4] effective medium theory, the experiments using optical methods, such as Putnam et al. [15], Rusconi et al. [16] and Venerus et al. [17], are also consistent with the conclusion of this study.

In this article, a contextual notation is introduced to distinguish between dimensional and dimensionless variables and parameters. The contextual notation implies that an asterisk subscript is used to identify dimensional variables and parameters only when ambiguity may arise when the asterisk subscript is not used. For example t* is the dimensional time, while t is its corresponding dimensionless counterpart. However, kf is the effective fluid phase thermal conductivity, a dimensional parameter that appears without an asterisk subscript without causing ambiguity.

Problem formulation

The theoretical model derived by Vadasz [22] to investigate the transient heat conduction in a fluid containing suspended solid particles by considering phase-averaged equations will be presented only briefly without including the details that can be obtained from [22]. The phase-averaged equations are

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M1">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M2">View MathML</a>


where t* is time, Tf (r*,t*), and Ts (r*,t*) are temperature values for the fluid and solid phases, respectively, averaged over a representative elementary volume (REV) that is large enough to be statistically valid but sufficiently small compared to the size of the domain, and where r* are the coordinates of the centroid of the REV. In Equations (1) and (2), γs = ερscs and γf = (1 - ε)ρfcp represent the effective heat capacity of the solid and fluid phases, respectively; with ρs and ρf are the densities of the solid and fluid phases, respectively; cs and cp are the specific heats of the solid and fluid phases, respectively; and ε is the volumetric solid fraction of the suspension. Similarly, kf is the effective thermal conductivity of the fluid that may be defined in the form <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M3">View MathML</a>, where <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M4">View MathML</a> is the thermal conductivity of the fluid, <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M5">View MathML</a> is the thermal conductivity ratio, and ε is the solid fraction of suspended particles in the suspension. In Equations (1) and (2), the parameter h, carrying units of W m-3 K-1, represents an integral heat transfer coefficient for the contribution of the heat conduction at the solid-fluid interface as a volumetric heat source/sink within an REV. It is assumed to be independent of time, and its general relationship to the surface-area-to-volume ratio (specific area) was derived in [22]. Note that Ts(r*,t*) is a function of the space variables represented by the position vector <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M6">View MathML</a>, in addition to its dependence on time, because Ts(r*,t*) depends on Tf(r*,t*) as explicitly stated in Equation (1), although no spatial derivatives appear in Equation (1). There is a lack of macroscopic level conduction mechanism in Equation (1) representing the heat transfer within the solid phase because the solid particles represent the dispersed phase in the fluid suspension, and therefore the solid particles can conduct heat between themselves only via the neighbouring fluid. When steady state is accomplished ∂Ts/∂t* = ∂Tf/∂t* = 0, leading to local thermal equilibrium between the solid and fluid phases, i.e. Ts(r) = Tf(r).

For the case of a thin hot wire embedded in a cylindrical container insulated on its top and bottom one can assume that the heat is transferred in the radial direction only, r*, rendering Equation (2) into

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M7">View MathML</a>


In a homogeneous medium without solid-suspended particles, Equation (1) is not relevant and the last term in Equation (3) can also be omitted. The boundary and initial conditions applicable are an initial ambient constant temperature, TC, within the whole domain, an ambient constant temperature, TC, at the outer radius of the container and a constant heat flux, q0, over the fluid-wire interface that is related to the Joule heating of the wire in the form q0 = iV/(πdw*l*), where dw* and l* are the diameter and the length of the wire respectively, i is the electric current and V is the voltage drop across the wire. Vadasz [22] showed that the problem formulated by Equations (1) and (3) subject to appropriate initial and boundary conditions represents a particular case of Dual-Phase-Lagging heat conduction (see also [24-28]).

An essential component in the application of the THW method for estimating experimentally the effective thermal conductivity of the nanofluid suspension is the assumption that the nanofluid suspension behaves basically like a homogeneous material following Fourier law for the bulk. The THW method is well established as the most accurate, reliable and robust technique [29] for evaluating the thermal conductivity of fluids. A very thin (5-80 μm in diameter) platinum (alternatively tantalum) wire is embedded vertically in the selected fluid and serves as a heat source as well as a thermometer (see [22] for details). Because of the very small diameter and high thermal conductivity of the platinum wire, it can be regarded as a line heat source in an otherwise infinite cylindrical medium. The rate of heat generated per unit length (l*) of platinum wire due to a step change in voltage is therefore <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M8">View MathML</a> W m-1. Solving for the radial heat conduction due to this line heat source leads to an approximated temperature solution in the wire's neighbourhood in the form

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M9">View MathML</a>


provided a validity condition for the approximation is enforced, i.e. <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M10">View MathML</a>, where rw* is the radius of the platinum wire, <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M11">View MathML</a> is the fluid's thermal diffusivity, and γ0 = 0.5772156649 is Euler's constant. Equation (4) reveals a linear relationship, on a logarithmic time scale, between the temperature and time. Therefore, one way of evaluating the thermal conductivity is from the slope of this relationship evaluated at r* = rw*. For any two readings of temperature, T1 and T2, recorded at times t1* and t2* respectively, the thermal conductivity can be approximated using Equation (4) in the form:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M12">View MathML</a>


Equation (5) is a very accurate way of estimating the thermal conductivity as long as the validity condition is fulfilled. The validity condition implies the application of Equation (5) for long times only. However, when evaluating this condition to data used in the nanofluid suspensions experiments, one obtains that t0* ~ 6 ms, and the time beyond which the solution (5) can be used reliably is therefore of the order of hundreds of milliseconds, not so long in the actual practical sense.

Two methods of solution

While the THW method is well established for homogeneous fluids, its applicability to two-phase systems such as fluid suspensions is still under development, and no reliable validity conditions for the latter exist so far (see Vadasz [30] for a discussion and initial study on the latter). As a result, one needs to refer to the two-equation model presented by Equations (1) and (3), instead of the one Fourier type equation that is applicable to homogeneous media.

Two methods of solution are in principle available to solve the system of Equations (1) and (3). The first is the elimination method while the second is the eigenvectors method. By means of the elimination method, one may eliminate Tf from Equation (1) in the form:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M13">View MathML</a>


and substitute it into Equation (3) hence rendering the two Equations (1) and (3), each of which depends on both Ts and Tf, into separate equations for Ts and Tf, respectively, in the form:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M14">View MathML</a>


where the index i takes the values i = s for the solid phase and i = f for the fluid phase, and the following notation was used:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M15">View MathML</a>


In Equation (8), τq and τT are the heat flux and temperature-related time lags linked to Dual-Phase-Lagging [22,24-27,31], while αe is the effective thermal diffusivity of the suspension. The resulting Equation (7) is identical for both fluid and solid phases. Vadasz [22] used this equation in providing the solution. The initial conditions applicable to the problem at hand are identical for both phases, i.e. both phases' temperatures are set to be equal to the ambient temperature TC

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M16">View MathML</a>


The boundary conditions are

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M17">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M18">View MathML</a>


where r0* is the radius of the cylindrical container. Equation (7) is second-order in time and second-order in space. The initial conditions (9) provide one such condition for each phase while the second-order Equation (7) requires two such conditions. To obtain the additional initial conditions, one may use Equations (1) and (3) in combination with (9). From (9), it is evident that both phases' initial temperatures at t* = 0 are identical and constant. Therefore, <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M104','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M104">View MathML</a>, leading to <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M105','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M105">View MathML</a> and <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M106','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M106">View MathML</a> to be substituted in (1) and (3), which in turn leads to the following additional initial conditions for each phase:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M19">View MathML</a>


The two boundary conditions (10) and (11) are sufficient to uniquely define the problem for the fluid phase; however, there are no boundary conditions set for the solid phase as the original Equation (1) for the solid phase had no spatial derivatives and did not require boundary conditions. To obtain the corresponding boundary conditions for the solid phase, which are required for the solution of Equation (7) corresponding to i = s, one may use first the fact that at r* = r0* both phases are exposed to the ambient temperature and therefore one may set

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M20">View MathML</a>


Second, one may use Equation (6) and taking its derivative with respect to r* yields

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M21">View MathML</a>


In Equation (14), the spatial variable r* plays no active role; it may therefore be regarded as a parameter. As a result, one may present Equation (14) for any specified value of r*. Choosing r* = rw* where the value of <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M107','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M107">View MathML</a> is known from the boundary condition (11), yields from (14) the following ordinary differential equation:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M22">View MathML</a>


At steady state, Equation (15) produces the solution

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M23">View MathML</a>


where Ts,st is the steady-state solution. The transient solution Ts,tr = Ts - Ts,st satisfies then the equation:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M24">View MathML</a>


subject to the initial condition

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M25">View MathML</a>


because <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M108','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M108">View MathML</a> for all values of <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M109','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M109">View MathML</a> given that according to (9) at t* = 0: <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M110','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M110">View MathML</a>. Equation (17) can be integrated to yield

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M26">View MathML</a>


which combined with the initial condition (18) produces the value of the integration constant A = 0 and therefore the transient solution becomes

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M27">View MathML</a>


The complete solution for the solid temperature gradient at the wire is therefore obtained by combining (20) with (16) leading to

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M28','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M28">View MathML</a>


producing the second boundary condition for the solid phase, which is identical to the corresponding boundary condition for the fluid phase. One may therefore conclude that the solution to the problem formulated in terms of Equation (7) that is identical to both phases, subject to initial conditions (9) and (12) that are identical to both phases, and boundary conditions (10), (11), and (13), (21) that are also identical to both phases, should be also identical to both phases, i.e. Ts (t*,r*) = Tf (t*,r*). This, however, may not happen because then Tf - Ts = 0 leads to conflicting results when substituted into (1) and (3). The result obtained here is identical to Vadasz [32] who demonstrated that a paradox revealed by Vadasz [33] can be avoided only by refraining from using this method of solution. While the paradox is revealed in the corresponding problem of a porous medium subject to a combination of Dirichlet and insulation boundary conditions, the latter may be applicable to fluids suspensions by setting the effective thermal conductivity of the solid phase to be zero. The fact that in the present case the boundary conditions differ, i.e. a constant heat flux is applied on one of the boundaries (such a boundary condition would have eliminated the paradox in porous media), does not eliminate the paradox in fluid suspensions mainly because in the latter case the steady-state solution is identical for both phases. In the porous media problem, the constant heat flux boundary condition leads to different solutions at steady state, and therefore the solutions for each phase even during the transient conditions differ.

The elimination method yields the same identical equation with identical boundary and initial conditions for both phases apparently leading to the wrong conclusion that the temperature of both phases should therefore be the same. A closer inspection shows that the discontinuity occurring on the boundaries' temperatures at t = 0, when a "ramp-type" of boundary condition is used, is the reason behind the occurring problem and the apparent paradox. The question that still remains is which phase temperature corresponds to the solution presented by Vadasz [22]; the fluid or the solid phase temperature?

By applying the eigenvectors method as presented by Vadasz [32], one may avoid the paradoxical solution and obtain both phases temperatures. The analytical solution to the problem using the eigenvectors method is obtained following the transformation of the equations into a dimensionless form by introducing the following dimensionless variables:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M29','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M29">View MathML</a>


where the following two dimensionless groups emerged:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M30','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M30">View MathML</a>


representing a heat flux Fourier number and a temperature Fourier number, respectively. The ratio between them is identical to the ratio between the time lags, i.e.

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M31','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M31">View MathML</a>


Equations (1) and (3) expressed in a dimensionless form using the transformation listed above are

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M32','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M32">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M33','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M33">View MathML</a>


where the following additional dimensionless groups emerged:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M34','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M34">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M35','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M35">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M36">View MathML</a>


where Nif is the fluid phase Nield number. The solutions to Equations (25) and (26) are subject to the following initial and boundary conditions obtained from (9), (10) and (11) transformed in a dimensionless form:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M37','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M37">View MathML</a>


The boundary conditions are

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M38','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M38">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M39','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M39">View MathML</a>


No boundary conditions are required for θs. The solution to the system of Equations (25)-(26) is obtained by a superposition of steady and transient solutions θi,st(r) and θi,tr (t,r), respectively, in the form:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M40','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M40">View MathML</a>


Substituting (33) into (25)-(26) yields to the following equations for the steady state:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M41">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M42">View MathML</a>


leading to the following steady solutions which satisfy the boundary conditions (31) and (32):

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M43','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M43">View MathML</a>


The transient part of the solutions θi,tr (t,r) can be obtained by using separation of variables leading to the following form of the complete solution:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M44">View MathML</a>


Substituting (37) into (25)-(26) yields, due to the separation of variables, the following equation for the unknown functions Ron (r):

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M45','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M45">View MathML</a>


subject to the boundary conditions

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M46','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M46">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M47','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M47">View MathML</a>


and the following system of equations for the unknown functions Sin (t), (i = s,f), i.e.

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M48','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M48">View MathML</a>



<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M49','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M49">View MathML</a>


and where the separation constant <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M111','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M111">View MathML</a> represents the eigenvalues in space.

Equation (38) is the Bessel equation of order 0 producing solutions in the form of Bessel functions

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M50','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M50">View MathML</a>


Where J0(κnr) and Y0(κnr) are the order 0 Bessel functions of the first and second kind, respectively. The solution (43) satisfies the boundary condition (39) as can easily be observed by substituting r = 1 in (43). Imposing the second boundary condition (40) yields a transcendental equation for the eigenvalues κn in the form:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M51','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M51">View MathML</a>


where J1(κnrw) and Y1(κnrw) are the order 1 Bessel functions of the first and second kind, respectively, evaluated at r = rw. The compete solution is obtained by substituting (43) into (37) and imposing the initial conditions (30) in the form

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M52','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M52">View MathML</a>


At t = 0, both phases' temperatures are the same leading to the conclusion that

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M53','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M53">View MathML</a>


Multiplying (45) by the orthogonal eigenfunction Rom (κm ,r) with respect to the weight function r and integrating the result over the domain [rw,1], i.e. <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M54','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M54">View MathML</a> yield

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M55','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M55">View MathML</a>


The integral on the right-hand side of (47) produces the following result due to the orthogonality conditions for Bessel functions:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M56','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M56">View MathML</a>


where the norm N(κn) is evaluated in the form:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M57','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M57">View MathML</a>


The integral on the left-hand side of (47) can be evaluated using integration by parts and the equation for the eigenvalues (44) to yield

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M58','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M58">View MathML</a>


Substituting (48) and (50) into (47) yields the values of Sin at t = 0, i.e. Sno = Ssn(0) = Sfn(0)

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M59','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M59">View MathML</a>

that need to be used as initial conditions for the solution of system (41)

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M60','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M60">View MathML</a>


to produce the explicit solutions in time. With the initial conditions for Sin evaluated (i = s,f), one may turn to solving system (41) that can be presented in the following vector form:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M61','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M61">View MathML</a>


where the matrix A is explicitly defined by

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M62','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M62">View MathML</a>


with the values of a,c and dn given by Equation (42), and the vector Sn defined in the form Sn = [Ssn,Sfn]T. The eigenvalues λn corresponding to (52) are obtained as the roots of the following quadratic algebraic equation:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M63','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M63">View MathML</a>


leading to

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M64','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M64">View MathML</a>

which upon substituting a,c and dn from Equation (42) yields

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M65','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M65">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M66','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M66">View MathML</a>


The following useful relationship is obtained from (55) and (56):

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M67','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M67">View MathML</a>


The corresponding eigenvectors υ1n and υ2n are evaluated in the form:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M68','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M68">View MathML</a>


leading to the following solution:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M69','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M69">View MathML</a>


and explicitly following the substitution of (58) and the initial conditions Sin (i = s,f), at t = 0, i.e. Ssn(0) = Sfn(0) = Sno with the values of Sno given by Equation (51)

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M70','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M70">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M71','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M71">View MathML</a>


Substituting (57) into (60) and (61) and the latter into the complete solution (37) yields

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M72','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M72">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M73','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M73">View MathML</a>


where Bn is

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M74','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M74">View MathML</a>


Comparing the solutions obtained above with the solution obtained by Vadasz [22] via the elimination method, one may conclude that the latter corresponds to the solid phase temperature θs.

The Fourier solution is presented now to compare the solution obtained from the Dual-Phase-Lagging model to the former. The Fourier solution is the result obtained by solving the thermal diffusion equation

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M75','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M75">View MathML</a>


subject to the boundary and initial conditions

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M76','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M76">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M77','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M77">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M78','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M78">View MathML</a>


where the same scaling as in Equation (22) was applied in transforming the equation into its dimensionless form, hence the reason for the coefficient 1/β in the equation. The Fourier solution for this problem has then the form [34]

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M79','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M79">View MathML</a>



<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M80','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M80">View MathML</a>


and the eigenvalues κn are the solution of the same transcendental Equation (44) and the eigenfunctions Ron(r) are also identical to the ones presented in Equation (43). The relationship between the Fourier coefficient Cn and the Dual-Phase-Lagging model's coefficient Bn is

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M81','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M81">View MathML</a>


Correction of the THW results

When evaluating the thermal conductivity by applying the THW method and using Fourier law, one obtains for the effective thermal conductivity the following relationship [22]:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M82','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M82">View MathML</a>


where the temperature difference [Tw(t) - TC] is represented by the recorded experimental data, and the value of the heat flux at the fluid-platinum-wire interface q0 is evaluated from the Joule heating of the hot wire. In Equation (72) <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M83','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M83">View MathML</a>, where the coefficient Cn is defined by (70) and the eigenvalues κn are defined by Equation (44). Note that the definition of Cn here is different than in [22]. The results obtained from the application of Equation (72) fit extremely well the approximation used by the THW method via Equation (5) within the validity limits of the approximation (5). Therefore, the THW method is extremely accurate for homogeneous materials.

On the other hand, for non-homogeneous materials, by means of the solutions (62) and (63) applicable to fluid suspensions evaluated at r = rw, one obtains

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M84','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M84">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M85','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M85">View MathML</a>


where kf,act is the actual effective thermal conductivity, Tsw (t) and Tfw (t) are the solid and fluid phases temperatures "felt" by the wire at the points of contact with each phase, respectively, and the functions gs (t) and gf (t) obtained from the solutions (62) and (63) evaluated at r = rw take the form

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M86','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M86">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M87','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M87">View MathML</a>


When the wire is exposed partly to the fluid phase and partly to the solid phase, there is no justification in assuming that the wire temperature is uniform: on the contrary the wire temperature will vary between the regions exposed to the fluid and solid phases. Assuming that some solid nanoparticles are in contact with the wire in a way that they form approximately "solid rings" around the wire, then the "effective" wire temperature can be evaluated as electrical resistances in series. By defining the relative wire area covered by the solid nanoparticles as as = As/Atot = As/2πrw*l* its corresponding wire area covered by the fluid is af = Af/Atot = 1 - as, then from the relationship between the electrical resistance and temperature accounting for electrical resistances connected in series, one obtains an expression for the effective wire temperature (i.e. the temperature that is evaluated using the wire's lumped electrical resistance in the THW Wheatstone bridge) Tw in the form:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M88','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M88">View MathML</a>


Substituting (73) and (74) into (77) yields

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M89','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M89">View MathML</a>


One may then use (78) to evaluate the actual nanofluid's effective thermal conductivity kf,act from (78) in the form

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M90','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M90">View MathML</a>


When using the single phase Fourier solution (72) applicable for homogeneous materials to evaluate the effective thermal conductivity of non-homogeneous materials like nanofluid suspensions instead of using Equation (79), one obtains a value that differs from the actual one by a factor of

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M91','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M91">View MathML</a>


where kf,app is the apparent effective thermal conductivity obtained from the single phase Fourier conduction solution while kf,act is the actual effective thermal conductivity that corresponds to data that follow a Dual-Phase-Lagging conduction according to the derivations presented above. The ratio between the two provides a correction factor for the deviation of the apparent effective thermal conductivity from the actual one. This correction factor when multiplied by the ratio <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M92','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M92">View MathML</a> produces the results for <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M93','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M93">View MathML</a>, where <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M94','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M94">View MathML</a> is the thermal conductivity of the base fluid without the suspended particles, and kf,act is the effective thermal conductivity evaluated using Maxwell's [4] theory, which for spherical particles can be expressed in the form:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M95','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M95">View MathML</a>


where kf,act is Maxwell's effective thermal conductivity, <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M96','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M96">View MathML</a> is the ratio between the thermal conductivity of the solid phase and the thermal conductivity of the base fluid, and ε is the volumetric solid fraction of the suspension. Then, these results of <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M97','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M97">View MathML</a> can be compared with the experimental results presented by Liu et al. [23].

Results and discussion

The results for the solid and fluid phases' temperature at r = rw as a function of time obtained from the solutions (62) and (63) are presented in Figures 1, 2 and 3 in comparison with the single-phase Fourier solution (69) for three different combinations of values of Foq and as, and plotted on a logarithmic time scale. While the quantitative results differ amongst the three figures, there are some similar qualitative features that are important to mention. First, it is evident from these figures that the fluid phase temperature is almost the same as the temperature obtained from the single-phase Fourier solution. Second, it is also evident that the solid phase temperature lags behind the fluid phase temperature by a substantial difference. They become closer as steady-state conditions approach. It is therefore imperative to conclude that the only way, an excessively higher effective thermal conductivity of the nanofluid suspension as obtained by Eastman et al. [1], Lee et al. [2] and Choi et al. [3] could have been obtained even in an apparent form, is if the wire was excessively exposed to the solid phase temperature. The latter could have occurred if the electric current passing through the wire created electric fields that activated a possible mechanism of electrophoresis that attracted the suspended nanoparticles towards the wire. Note that such a mechanism does not cause agglomeration in the usual sense of the word, because as soon as the electric field ceases, the agglomeration does not have to persist and the particles can move freely from the wire's surface. Therefore, testing the wire's surface after such an experiment for evidence of agglomeration on the wire's surface may not necessarily produce the required evidence for the latter.

thumbnailFigure 1. Dimensionless wire temperature. Comparison between the Fourier and Dual-Phase-Lagging solutions for the following dimensionless parameters values Foq = 1.45 × 10-2 and as = 0.45.

thumbnailFigure 2. Dimensionless wire temperature. Comparison between the Fourier and Dual-Phase-Lagging solutions for the following dimensionless parameters values Foq = 1.1 × 10-2 and as = 0.55.

thumbnailFigure 3. Dimensionless wire temperature. Comparison between the Fourier and Dual-Phase-Lagging solutions for the following dimensionless parameters values Foq = 6 × 10-3 and as = 0.35.

Liu et al. [23] used a very similar THW experimental method as the one used by Eastman et al. [1], Lee et al. [2] and Choi et al. [3] with the major distinction being in the method of producing the nanoparticles and a cylindrical container of different dimensions. They used water as the base fluid and Cu nanoparticles as the suspended elements at volumetric solid fractions of 0.1 and 0.2%. Their data that are relevant to the present discussion were digitized from their Figure 3[23] and used in the following presentation to compare our theoretical results. Three specimen data are presented in Figure 3[23] resulting in extensive overlap of the various curves, and therefore in some digitizing error which is difficult to estimate when using only this figure to capture the data.

The comparison between the theoretical results presented in this article with the experimental data [23] is presented in Figures 4, 5 and 6. The separation of these results into three different figures aims to better distinguish between the different curves and avoid overlapping as well as presenting the results on their appropriate scales. Figure 4 presents the results that are applicable to specimen No. 4 in Liu et al. [23] and corresponding to values of Foq = 1.45 × 10-2 and as = 0.45 in the theoretical model. Evaluating Maxwell's [4] effective thermal conductivity for specimen No. 4 leads to a value of 0.6018 W/mK, which is higher by 0.3% than that of the base fluid (water), i.e. <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M98','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M98">View MathML</a>. From the figure, it is evident that the theoretical results match very well with the digitized experimental data. Furthermore, the steady-state result for the ratio between the effective thermal conductivity and that of the base fluid was estimated from the digitized data to be <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M99','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M99">View MathML</a> clearly validating Maxwell's [4] predicted value. The results applicable to specimen No. 5 in Liu et al. [23] and corresponding to values of Foq = 1.1 × 10-2 and as = 0.55 in the theoretical model are presented in Figure 5. The very good match between the theory and the digitized experimental data is again evident. In addition, the ratio between the effective thermal conductivity and that of the base fluid was estimated from the digitized data to be <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M100','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M100">View MathML</a> again validating Maxwell's [4] predicted value of <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M101','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M101">View MathML</a>. The last result is presented in Figure 6, which corresponds to specimen No. 9 in Liu et al. [23] and to values of Foq = 6 × 10-3 and as = 0.35 in the theoretical model. The results are presented on an appropriately scaled vertical axis and show again a very good match between the theory presented in this article, and the experimental data as digitized from Liu et al. [23]. Since the volumetric solid fraction for this specimen was 0.2%, its corresponding Maxwell's [4] effective thermal conductivity for this specimen leads to a value of 0.6036 W/mK, which is higher by 0.6% than that of the base fluid (water), i.e. <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M102','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M102">View MathML</a>. The steady-state result for the ratio between the effective thermal conductivity and that of the base fluid was estimated from the digitized data to be <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/154/mathml/M103','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/154/mathml/M103">View MathML</a> validating again Maxwell's [4] predicted value.

thumbnailFigure 4. Comparison of the present theory with experimental data of Liu et al. [23] (here redrawn from published data) of the effective thermal conductivity ratio for conditions compatible with specimen No. 4, leading to a Fourier number of Foq = 1.45 × 10-2 and a solid particles to total wire area ratio of as = 0.45.

thumbnailFigure 5. Comparison of the present theory with experimental data of Liu et al. [23] (here redrawn from published data) of the effective thermal conductivity ratio for conditions compatible with specimen No. 5, leading to a Fourier number of Foq = 1.1 × 10-2 and a solid particles to total wire area ratio of as = 0.55.

thumbnailFigure 6. Comparison of the present theory with experimental data of Liu et al. [12] (here redrawn from published data) of the effective thermal conductivity ratio for conditions compatible with specimen No. 9, leading to a Fourier number of Foq = 6 × 10-3 and a solid particles to total wire area ratio of as = 0.35.

It should be mentioned that Liu et al. [23] explain their time-dependent effective thermal conductivity by claiming that it was caused by nanoparticle agglomeration, a conclusion that is consistent with the theoretical results of this study.


The theoretical results derived in this article combined with experimental data [23] lead to the conclusion that, while there is no improvement in the effective thermal conductivity of nanofluids beyond the Maxwell's effective medium theory [4], there is nevertheless the possibility of substantial heat transfer augmentation via nanofins. Nanoparticles attaching to the hot wire by a mechanism that could be related to electrophoresis depending on the strength of the electrical current passing through the wire suggests that such attachments can be deliberately designed and produced on any heat transfer surface to yield an agglomeration of nanofins that exchange heat effectively because of the extremely high heat transfer area as well as the flexibility of such nanofins to bend in the fluid's direction when fluid motion is present, hence extending its applicability to include a new, and what appears to be a very effective, type of heat convection. A quantitative estimate of the effectiveness of nanofins requires, however, an extension of the model presented in this article to include heat conduction within the nanofins.


REV: representative elementary volume; THW: transient-hot-wire.

Competing interests

The author declares that they have no competing interests.

Authors' contributions

PV conceived and carried out all work reported in this paper.


  1. Eastman JA, Choi SUS, Li S, Yu W, Thompson LJ: "Anomalously increased effective thermal conductivities of ethylene glycol-based nanofluids containing copper nanoparticles".

    Appl Phys Lett 2001, 78:718-720. Publisher Full Text OpenURL

  2. Lee S, Choi SUS, Li S, Eastman JA: "Measuring thermal conductivity of fluids containing oxide nanoparticles".

    J Heat Transfer 1999, 121:280-289. Publisher Full Text OpenURL

  3. Choi SUS, Zhang ZG, Yu W, Lockwood FE, Grulke EA: "Anomalous thermal conductivity enhancement in nanotube suspensions".

    Appl Phys Lett 2001, 79:2252-2254. Publisher Full Text OpenURL

  4. Maxwell JC: A Treatise on Electricity and Magnetism. 3rd edition. Dover, New York: Clarendon Press; 1891:435-441.

    1954 (reprint)


  5. Hamilton RL, Crosser OK: "Thermal conductivity of heterogeneous two-component systems".

    I&EC Fundam 1962, 1:187-191. OpenURL

  6. Jeffrey DJ: "Conduction through a random suspension of spheres".

    Proc R Soc Lond A 1973, 335:355-367. Publisher Full Text OpenURL

  7. Davis RH: "The effective thermal conductivity of a composite material with spherical inclusions".

    Int J Thermophys 1986, 7:609-620. Publisher Full Text OpenURL

  8. Lu S, Lin H: "Effective conductivity of composites containing aligned spheroidal inclusions of finite conductivity".

    J Appl Phys 1996, 79:6761-6769. Publisher Full Text OpenURL

  9. Bonnecaze RT, Brady JF: "A method for determining the effective conductivity of dispersions of particles".

    Proc R Soc Lond A 1990, 430:285-313. Publisher Full Text OpenURL

  10. Bonnecaze RT, Brady JF: "The effective conductivity of random suspensions of spherical particles".

    Proc R Soc Lond A 1991, 432:445-465. Publisher Full Text OpenURL

  11. Das KS, Putra N, Thiesen P, Roetzel W: "Temperature dependence of thermal conductivity enhancement for nanofluids".

    J Heat Transfer 2003, 125:567. Publisher Full Text OpenURL

  12. Li CH, Peterson GP: "Experimental investigation of temperature and volume fraction variation on the effective thermal conductivity of nanoparticles suspensions (nanofluids)".

    J Appl Phys 2006, 99:084314. Publisher Full Text OpenURL

  13. Buongiorno J, Venerus DC: "Letter to the Editor".

    Int J Heat Mass Transfer 2010, 53:2939-2940. Publisher Full Text OpenURL

  14. Buongiorno J, Venerus DC, Prabhat N, McKrell T, Townsend J, Christianson R, Tolmachev YV, Keblinski P, Hu LH, Alvarado JL, Bang IC, Bishnoi SW, Bonetti M, Botz F, Cecere A, et al.: "A benchmark study on the thermal conductivity of nanofluids".

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

  15. Putnam SA, Cahill DG, Braun PV, Ge Z, Shimmin RG: "Thermal conductivity of nanoparticle suspensions".

    J Appl Phys 2006, 99:084308. Publisher Full Text OpenURL

  16. Rusconi R, Rodari E, Piazza R: "Optical measurements of the thermal properties of nanofluids".

    Appl Phys Lett 2006, 89:261916. Publisher Full Text OpenURL

  17. Venerus DC, Kabadi MS, Lee S, Perez-Luna V: "Study of thermal transport in nanoparticle suspensions using forced Rayleigh scattering".

    J Appl Phys 2006, 100:094310. Publisher Full Text OpenURL

  18. Jang SP, Choi SU: "Role of Brownian motion in the enhanced thermal conductivity of nanofluids".

    Appl Phys Lett 2004, 84:4316-4318. Publisher Full Text OpenURL

  19. Prasher R, Bhattacharya P, Phelan PE: "Thermal conductivity of nanoscale colloidal solutions (nanofluids)".

    Phys Rev Lett 2005, 94:025901. PubMed Abstract | Publisher Full Text OpenURL

  20. Vadasz JJ, Govender S, Vadasz P: "Heat Transfer Enhancement in Nanofluids Suspensions: Possible Mechanisms and Explanations".

    Int J Heat Mass Transfer 2005, 48:2673-2683. Publisher Full Text OpenURL

  21. Vadasz JJ, Govender S: "Thermal wave effects on heat transfer enhancement in nanofluids suspensions".

    Int J Therm Sci 2010, 49:235-242. Publisher Full Text OpenURL

  22. Vadasz P: "Heat conduction in nanofluid suspensions".

    J Heat Transfer 2006, 128:465-477. Publisher Full Text OpenURL

  23. Liu MS, Lin MCC, Tsai CY, Wang CC: "Enhancement of thermal conductivity with Cu for nanofluids using chemical reduction method".

    Int J Heat Mass Transfer 2006, 49:3028-3033. Publisher Full Text OpenURL

  24. Wang L: "Solution structure of hyperbolic heat-conduction equation".

    Int J Heat Mass Transfer 2000, 43:365-373. Publisher Full Text OpenURL

  25. Wang L, Xu M, Zhou X: "Well-posedness and solution structure of dual-phase-lagging heat conduction".

    Int J Heat Mass Transfer 2001, 44:1659-1669. Publisher Full Text OpenURL

  26. Wang L, Xu M: "Well-posed problem of dual-phase-lagging heat conduction equation in 2D and 3D regions".

    Int J Heat Mass Transfer 2002, 45:1065-1071. OpenURL

  27. Xu M, Wang L: "Thermal oscillation and resonance in dual-phase-lagging heat conduction".

    Int J Heat Mass Transfer 2002, 45:1055-1061. Publisher Full Text OpenURL

  28. Cheng L, Xu M, Wang L: "From Boltzmann transport equation to single-phase-lagging heat conduction".

    Int J Heat Mass Transfer 2008, 51:6018-6023. Publisher Full Text OpenURL

  29. Hammerschmidt U, Sabuga W: "Transient Hot Wire (THW) method: Uncertainty assessment".

    Int J Thermophys 2000, 21:1255-1278. Publisher Full Text OpenURL

  30. Vadasz P: "Rendering the Transient Hot Wire Experimental Method for Thermal Conductivity Estimation to Two-Phase Systems - Theoretical Leading Order Results".

    J Heat Transfer 2010, 132:081601. Publisher Full Text OpenURL

  31. Tzou DY: Macro-to-Microscale Heat Transfer: The Lagging Behavior. Washington, DC: Taylor&Francis; 1997. OpenURL

  32. Vadasz P: "On the Paradox of Heat Conduction in Porous Media Subject to Lack of Local Thermal Equilibrium".

    Int J Heat Mass Transfer 2007, 50:4131-4140. Publisher Full Text OpenURL

  33. Vadasz P: "Explicit Conditions for Local Thermal Equilibrium in Porous Media Heat Conduction".

    Trans Porous Media 2005, 59:341-355. Publisher Full Text OpenURL

  34. Özisik MN: Heat Conduction. 2nd edition. New York: John Wiley & Sons, Inc; 1993. OpenURL