Applications of carbon nanotube (CNT) and graphene in thermal management have recently attracted significant attention. However, the lack of efficient prediction formula for heat transfer coefficient between nanomaterials and gas environment limits the further development of this technique. In this work, a kinetic model has been established to predict the heat transfer coefficient of an individual CNT in gas environment. The heat dissipation around the CNT is governed by molecular collisions, and outside the collision layer, the heat conduction is dominant. At nanoscales, the natural convection can be neglected. In order to describe the intermolecular collisions around the CNT quantitatively, a correction factor 1/24 is introduced and agrees well with the experimental observation. The prediction of the present model is in good agreement with our experimental results in free molecular regime. Further, a maximum heat transfer coefficient occurs at a critical diameter of several nanometers, providing guidelines on the practical design of CNT-based heat spreaders.
1. Introduction
As the scale of electronic devices continuously decreases to nanometers, the power density increases to a very high level; for example, an individual microprocessor operating at a low GHz frequency has a typical power density of 100 Wcm−2 [1]. Meanwhile, the waste heat at this high generation rate has to be dissipated efficiently in a limited space; otherwise, the resulting local high temperature may cause thermal failure of the device and highly reduce its lifetime. Hence, the thermal management at nanoscales has become a critical issue and attracted much attention. Recently, the carbon nanotube (CNT), nanowires, and graphene have been found to be the promising materials to remove the waste heat from nanoelectronic devices due to their super high thermal conductivities [2–5]. At two-dimensional scales, the graphene has important applications as heat spreaders and thermal interfacial materials [6–8]. Meanwhile, the CNT can be designed into a one-dimensional heat spreader connecting the hotspot with high power density and the heat sinks. These new techniques have potential applications in many areas, such as nanoelectromechanical systems [9], gas sensors [10], field emitters [11], and nanowire-based lasers [12, 13]. Since all the nanodevices will be used in the air environment, the study of heat transfer between the nanomaterials and surrounding air forms the foundation of the practical design of nano heat spreaders.
Heat transfer across the interface between macroscale solid and gas environment is well understood, where many experimental data and empirical formulas are available. However, as the characteristic size scales down to nanometers, it becomes a challenging task to determine the heat transfer coefficient (HTC) of nanomaterials quantitatively. In the experimental study, the commonly used resistance thermometer or thermocouples at macroscales are no longer applicable, and a noncontact measurement should be performed at the nanometer level in order to avoid the uncertainty of temperature measurement at interfaces. Recently, a noncontact technique of micro-Raman spectroscopy has been utilized to measure the thermal properties of CNT and graphene [9, 10, 14–16]. In this technique, a focused laser is used to heat the sample locally and the temperature response can be measured simultaneously from the downshifted G-band frequency. Combining with electrical measurement, both the optical absorption and thermal conductivity of an individual CNT can be obtained [16]. While this technique has the potential to provide both thermal conductivity and HTC of an individual CNT, the presented experimental data are still too few to conclude a convincing variation trend of the thermal property on a single parameter of CNT, like its length, diameter, or chirality.
In the theoretical study, no quantitative theoretical model has been fully developed to predict the HTC between an individual CNT and gas environment. Some molecular dynamics (MD) simulations of CNT-air interfacial thermal conductance have been reported and give a result of about 1.0 × 105 Wm−2 K−1 [17]. The associated thermal resistance is equivalent to the resistance of about 250 nm thick layer of air. Because of the high surface to volume ratios, the heat dissipation from an individual CNT to the gas environment becomes increasingly important compared with that at macroscales [5]. The heat dissipated across the CNT-air interface is more than 50% for a 2.04 nm, 5 μm long suspended CNT [18], and this fraction will approach 100% for CNTs longer than 60 μm [19]. Hence, the heat dissipation to the air environment can be exploited as an effective way for cooling CNT associated devices [20]. It is particularly important to develop an accurate prediction model for the HTC of individual CNT or nanowire, benefiting the practical thermal design of efficient heat spreaders at nanoscales. On the other hand, from CNT to nanowire, the characteristic diameter of HTC changes from several nanometers to hundreds of nanometers, ranging from the free molecular regime to the transition regime. The development of a unified theoretical model also has important academic significance.
In this work, a theoretical model of heat transfer across the CNT-gas interface has been established based on the kinetic theory. It demonstrates that a critical diameter exists at several nanometers for a maximum HTC of 1.08 × 105 Wm−2 K−1. The theoretical prediction agrees well with our experimental results in the free molecular regime [19].
2. Two-Layer Model of Heat Transfer
Fuchs proposed a general physical model to estimate the heat transfer between the gas environment and isolated particle at submicroscales [21, 22]. In this model, the whole space is separated into two layers: the layer around the particle with a thickness comparable to the mean free path (MFP) of molecules and the outside layer where the continuum theory is applicable. A schematic diagram of this two-layer model is shown in Figure 1.
Schematic diagram of the two-layer model.
Figure 1 shows the two-layer model for a spherical particle. The black circle in the center represents the solid part, around which a thin layer is the noncontinuum layer with a thickness of Δ. Ts, Tnc, and T∞ are the temperatures at the solid surface and at the interface of noncontinuum layer and of surrounding gas, respectively. rs and r∞ are the radiuses of solid particle and outer continuum layer when T=T∞, respectively. In this model, the heat transfer from the solid surface to the outer gas environment is finished in two steps: (1) energy exchange between the solid surface and gas molecules by collisions in the noncontinuum layer; (2) heat dissipation in the outer continuum layer. In Fuchs’ model, the first step heat transfer is described by molecular kinetics and the second step is described by means of continuum theory. Following this model, a governing equation for the HTC between the nanoparticles and surrounding gas is given as [23]
(1)h=Q4πrs2(Ts-T∞)=λfncr(ξ+γ)l[fncrξ2+β(ξ+γ)],
where h, Q, λ, and l are the heat coefficient, heat dissipated into the gas environment, thermal conductivity of gas, and molecular MFP, respectively; ξ=rs/l is the dimensionless radius, and γ=Δ/l is the dimensionless thickness of inner layer. fncr=ξ[1+ξeξEi(-ξ)] is the correction factor in the noncontinuum layer, where Ei(x) is the exponential integral of x. β is the correction factor for molecular collisions.
Based on (1), the Nusselt number Nu can be summarized as
(2)Nu=2rshλ=2fncrξ(ξ+γ)fncrξ2+β(ξ+γ).
This two-layer model is capable of predicting the HTC of particle with a diameter ranging from nanometers to micrometers, and the corresponding Knudsen number Kn is from 0.01 to 10. However, some questions still remain to limit the prediction accuracy: (1) an exact equivalent particle radius is needed for a good prediction, which is usually hard to get for aerosolized nanoparticles; (2) the molecular motion inside the inner layer is considered to be collisionless and there is no transition layer between the two layers in Figure 1 [24]. Ignoring the interactions between molecules may cause prediction uncertainty. Reference [23] takes into account the intermolecular collisions by changing the inner layer thickness and this thickness is assumed to increase as the particle diameter increases. But the accurate inner layer thickness is difficult to evaluate in practical applications.
In this work, we developed a prediction formula for the HTC of an individual CNT based on the two-layer model. The physical model shown in Figure 1 is applicable using a cylindrical coordinate system instead. The noncontinuum layer is adjacent to the CNT surface, where the kinetic theory is applied. The outside continuum layer is concentric with the CNT and the noncontinuum layer, where the continuum theory is applied. In the noncontinuum layer, the average velocity of molecules is given as [25]
(3)v-=2(2kBTπm)1/2,
where kB and m are the Boltzmann constant and mass of a molecule, respectively. In a cylindrical coordinate system, the number of molecules moving towards the solid surface with an incident angle θ is givens as
(4)N=nrdθdrdxcosθds2πrL,
where n, r, L, and ds are the molecular number density, radius, CNT length, and area element, respectively. Using the classical Maxwell distribution function, the fraction of molecules with a velocity between v and dv is given as
(5)Ndv=4π(m2πkBT)3/2v2exp(-mv22kBT)dv.
The molecular impact flux on the solid surface is given as [23]
(6)ψ=2mnπ2kBTl∫0∞(rsrs+Δ)2exp(-Δl)dΔ×∫0∞v2exp(-mv22kBT)dv=v-n2πξ[1+ξeξEi(-ξ)]=v-n2πfncr.
The dissipated energy Q at the wire surface can be described as
(7)Q=ψSαfCvm(Ts-Tnc).
Substituting the integrated result of ψ into (7), one can get
(8)Q=2αfdLnCvmfncrkBT2πm(Ts-Tnc),
where S, αf, Cv, Ts, Tnc, and d are the wire surface, accommodation coefficient, specific heat capacity at constant volume, solid surface temperature, noncontinuum layer temperature, and diameter of CNT, respectively.
According to the kinetic theory of gas, the thermal conductivity is given as [26]
(9)λ=2β′πv-n4Cvml=βnCvmlkBT2πm,
where l is the molecular MFP and β is the correction factor representing the effect of intermolecular collisions. The constant 2/π is caused by the different geometric coordinates of cylinder and sphere. Combining (8) and (9), the heat flux from the CNT surface to the noncontinuum layer is derived as
(10)qs-nc=αfλfncrβl(Ts-Tnc).
Figure 2 shows a schematic diagram of molecular collisions around an individual CNT. Because the diameter of CNT is much smaller than the molecular MFP, the probability of molecules colliding with CNT decreases as the diameter decreases. In this case, the surrounding gas appears to be rarefied for the CNT. As shown in Figure 2, the CNT is “covered” by a layer of molecules with a thickness of molecular MFP, referred to as the noncontinuum layer. Out of this layer, the gas is continuous and the collision cross-section is given as σ=πl2. To consider the probability that a molecule appears in the volume of CNT with a diameter of d, the cross-section is given as σCNT=0.25πd2. σCNT/σ equals zero as the diameter approaches zero, and it equals unity as the radius of CNT approaches the molecular MFP. The integration average value of σCNT/σ from zero to l equals 1/3. On the other hand, the reflected molecules from the CNT surface will collide with the other incident molecules within the layer of MFP, and the molecular interaction will decrease the collision probability. If every reflected molecule collides with the incoming one, that gives a limiting value of 1/2 for each direction. In a three-dimensional space, this fraction is 1/8. As a result, the correction factor β can be estimated as about 1/24.
Schematic diagram of molecular collision model.
In the continuum layer, the hear flux caused by heat conduction from radius r1 to radius r2 is q=(λ/rs)ΔT/ln(r2/r1), where rs is the solid wire radius. The heat conductance q/ΔT is inversely proportional to the radius rs. When rw=2μm, the HTC of natural convection is about 3.4 × 103 Wm−2 K−1 [26, 27], while the heat conductance q/ΔT is about 1.4 × 104 Wm−2 K−1, and the most heat is dissipated by conduction. As the radius rs decreases to submicrometers, the Grashof number Gr decreases rapidly below 1 × 10−7 in the third power of radius, implying that the buoyancy force is negligible compared with the viscous force and the heat conduction becomes the only dominant way. Governed by the one-dimensional heat conduction equation in cylindrical coordinates, the heat flux from the noncontinuum layer to the gas environment is given as
(11)qnc-∞=λrs+Δ(Tnc-T∞ln(nr(rs+r0)/(rs+Δ))),
where T∞ is the environment temperature and r∞=nr(rw+r0) is the radius when T=T∞. Imposing the continuity condition on heat flux at the interface between noncontinuum layer and continuum layer, the temperature Tnc is given as
(12)Tnc=Ts-[1+fncrβl(rs+Δ)ln(r∞rs+Δ)]-1(Ts-T∞).
Hence, we derived the HTC h as
(13)h=QπdL(Ts-T∞)=(2/π)αfnCvmfncrkBT/2πm1+(fncr/βl)(rs+Δ)ln(nr((rs+r0)/(rs+Δ))).
3. Results and Discussion
In order to verify the accuracy of the present theoretical model, we measured the HTCs of several single-walled CNTs (SWCNTs) using a Raman spectroscopy technique [19]. The experimental data and theoretical predictions are compared in Figure 3.
Comparison between the experimental data and theoretical prediction. The red dash line is the maximum value given by the kinetic theory.
As shown in Figure 3, the experimental data (red circles) match quite well the present model (solid curve) in a diameter range from 0.97 nm to 1.47 nm. Meanwhile, the literature data of microwires agree well with the present model as well [26]. When the accommodation coefficient αf equals unity, the maximum value of h is given as [28]
(14)hmax=58nukB=1.1×105Wm-2K-1,
where u=(3kBT/m)0.5 is the root mean square velocity. The maximum value predicted by (13) is 1.08 × 105 Wm−2 K−1, consistent with the prediction of (14). The best fitted parameters for the calculation result of Figure 3 are αf=1, β=0.04, nr=100, and r0=0.4l. As mentioned previously, the correction factor β can be calculated as(15)β=18∫02lσσCNTdx=18∫02lπx24πl2dx=124≈0.04,
and the best fitted value is consistent with the theoretical prediction, showing a good accuracy of the present model.
Figure 4 shows the ratio between the thermal resistances of noncontinuum layer and continuum layer, that is, Rnc and Rc, respectively. Here, the previous thermal resistance is calculated as Rnc=l/λ, where the thermal conductivity λ is calculated by (9) containing the correction factor β. The latter thermal resistance is calculated as Rc=(r∞/λair)ln[r∞/(rs+l)], where the thermal conductivity of air λair is used. It is seen in Figure 4 that the continuum layer contributes to the most thermal resistance when the diameter is at micrometers; in this case, the noncontinuum layer around the thin wire can be totally neglected. However, as the diameter decreases, the resistance Rnc becomes more and more important and the difference between Rnc and Rc decreases rapidly below 100 nm. At a critical diameter of about 5 nm, Rnc equals Rc and a maximum HTC value occurs as shown in Figure 3. When the diameter drops below the critical value, Rnc becomes larger than Rc and the HTC decreases as well. As a result, the existence of a maximum HTC is caused by the competition between Rnc and Rc.
Ratio between the thermal resistances of noncontinuum and continuum layers plotted with respect to the diameter. The red dash line marks the critical diameter when Rnc=Rc.
Figure 5 shows the HTC predicted by (13) with a correction factor β changing from 0.04 to 0.32. For a given temperature difference Ts-Tnc, the dissipated heat Q is proportional to the molecular number density n based on (8). As the correction factor β increases, the probability of intermolecular collisions inside the noncontinuum layer increases and the effective molecular density decreases. Thus the dissipated heat Q decreases and HTC h decreases as well. Changing the factor β only affects the HTC with a diameter below 100 nm, and the results of thick wires remain unchanged.
The effect of molecular interactions on the HTC.
Based on (13), a dimensionless Nusselt number Nu is given as
(16)Nu=hdλ=4αffncrξ/πβ+fncr(ξ+Δ/l)ln(nr((ξ+r0/l)/(ξ+Δ/l))).
The calculation results of (16) are given in Figure 6 as follows.
Nusselt number plotted with respect to Knudsen’s number, and different values of factor β are used for calculation. The regimes I, II, and III are the continuum regime, transition regime, and free molecular regime, respectively.
In Figure 6, the black solid and dash-dot curves are the prediction results of the present model for CNTs. The red dash curve and blue dot curve are the prediction results for nanoparticles, where the previous one takes into account the intermolecular collisions while the latter one does not [23]. Equation (2) relates the molecular interaction with the thickness of noncontinuum layer, and the increasing thickness decreases the dissipated heat at the particle surface. Similarly in the present model, as the CNT diameter increases, more molecules are reflected from the solid surface and collide with other molecules and the correction factor β increases. Using gas dynamics, the exact value of β can be estimated.
Figure 6 can be divided into three regimes, that is, the continuum regime (Kn<10-2), transition regime (10-2<Kn<10), and free molecular regime (Kn>10), respectively. In the continuum regime, the continuum theory is applicable and the thermal resistance induced by the molecular collisions at the solid surface is negligible. As the diameter decreases to the transition regime, the molecular interactions become more significant and reduce the Nu number accordingly. Finally in the free molecular regime, the intermolecular interactions within the MFP layer play a dominant role and decide the HTC. The present model given in (13) and (16) has been proved valid by the experiments in both free molecular and continuum regimes, while the experimental data in the transition regime are still lacking. In our future work, the HTC in a wide Kn number from 0.01 to 100 will be measured.
4. Summary
So far, there is no efficient theoretical model proposed to predict the HTC between an individual CNT and gas environment. In this work, a quantitative model has been established to predict HTC, and the prediction of the present model is in good agreement with the experimental data in both free molecular and continuum regimes.
The present model demonstrates that the natural convection can be neglected in the continuum layer outside the CNT, where the heat conduction plays a dominant role.
In the previous models for nanoparticles, it is difficult to evaluate the effect of intermolecular collisions on HTC around the solid surface. In the present model, the intermolecular effect can be quantitatively calculated using a correction factor β, and the analytical value is consistent with the experimental observation.
It is predicted by the present model that a maximum HTC exists in the transition regime, corresponding to a critical diameter at several nanometers. Similar to the concept of critical diameter of thermal insulation layer, the critical diameter here accounts for the dividing line of relative importance between the thermal resistances of noncontinuum layer and continuum layer. In practical applications, the performance of CNT-based heat spreader can be improved by choosing the critical diameters.
Acknowledgments
The authors would like to acknowledge the great help and support from Professor Zeng-Yuan Guo. The main idea of this work is derived from the discussion with Guo. This work is financially supported by the National Natural Science Foundation of China (Grant nos. 51076080 and 51136001), China Postdoctoral Science Foundation funded project, and Tsinghua University Initiative Scientific Research Program.
PopE.Energy dissipation and transport in nanoscale devices2010331471692-s2.0-7795363677210.1007/s12274-010-1019-zChenS.MooreA. L.CaiW.SukJ. W.AnJ.MishraC.AmosC.MagnusonC. W.KangJ.ShiL.RuoffR. S.Raman measurements of thermal transport in suspended monolayer graphene of variable sizes in vacuum and gaseous environments2011513213282-s2.0-7995541712710.1021/nn102915xKimP.ShiL.MajumdarA.McEuenP. L.Thermal transport measurements of individual multiwalled nanotubes2001872142-s2.0-0035914983215502RotkinS. V.PerebeinosV.PetrovA. G.AvourisP.An essential mechanism of heat dissipation in carbon nanotube electronics200995185018552-s2.0-6644912788310.1021/nl803835zChengC.FanW.CaoJ.RyuS.JiJ.GrigoropoulosC. P.WuJ.Heat transfer across the interface between nanoscale solids and gas201151210102101072-s2.0-8455520288310.1021/nn204072nBalandinA. A.Thermal properties of graphene and nanostructured carbon materials20111085695812-s2.0-7996064463110.1038/nmat3064ShahilK. M. F.BalandinA. A.Thermal properties of graphene and multilayer graphene: applications in thermal interface materials201215215133113402-s2.0-8486054669810.1016/j.ssc.2012.04.034NikaD. L.BalandinA. A.Two-dimensional phonon transport in graphene2012242323320310.1088/0953-8984/24/23/233203BlickR. H.GrifoniM.Focus on nano-electromechanical systems2005713672630StrelcovE.LilachY.KolmakovA.Gas sensor based on metal-insulator transition in VO2 nanowire thermistor200996232223262-s2.0-6684908407910.1021/nl900676nSheJ.XiaoZ.YangY.DengS.ChenJ.YangG.XuN.Correlation between resistance and field emission performance of individual ZnO one-dimensional nanostructures2008210201520222-s2.0-5684912460110.1021/nn800283uHuangM. H.MaoS.FeickH.YanH.WuY.KindH.WeberE.RussoR.YangP.Room-temperature ultraviolet nanowire nanolasers20012925523189718992-s2.0-003582730410.1126/science.1060367van VugtL. K.RuhleS.VanmaekelberghD.Phase-correlated nondirectional laser emission from the end facets of a ZnO nanowire2006612270727112-s2.0-3384634359010.1021/nl0616227DeshpandeV. V.HsiehS.BushmakerA. W.BockrathM.CroninS. B.Spatially resolved temperature measurements of electrically heated carbon nanotubes20091021042-s2.0-6314919772010.1103/PhysRevLett.102.105501105501BalandinA. A.GhoshS.BaoW.CalizoI.TeweldebrhanD.MiaoF.LauC. N.Superior thermal conductivity of single-layer graphene2008839029072-s2.0-4234908722510.1021/nl0731872HsuI.PettesM. T.BushmakerA.AykolM.ShiL.CroninS. B.Optical absorption and thermal transport of individual suspended carbon nanotube bundles2009925905942-s2.0-6524909193110.1021/n1802737qHuM.ShenoginS.KeblinskiP.RaravikarN.Thermal energy exchange between carbon nanotube and air2007902332-s2.0-3425075655910.1063/1.2746954231905HsuI.PettesM. T.AykolM.ShiL.CroninS. B.The effect of gas environment on electrical heating in suspended carbon nanotubes2010108842-s2.0-7814944486410.1063/1.3499256084307WangH. D.HuiL. J.GuoZ. Y.ZhangX.ZhangR. F.WeiF.LiT. Y.Thermal transport across the interface between a suspended single-wall carbon nanotube and airNanoscale and Microscale Thermophysical Engineering. In pressBeskokA.KarniadakisG. E.A model for flows in channels, pipes, and ducts at micro and nano scales19993143772-s2.0-0033239234FuchsN. A.On the stationary charge distribution on aerosol particles in a bipolar ionic atmosphere19635611851932-s2.0-3425058442610.1007/BF01993343FuchsN. A.1989New York, NY, USADoverKleinH. H.KarniJ.Ben-ZviR.BertocchiR.Heat transfer in a directly irradiated solar receiver/reactor for solid-gas reactions20078110122712392-s2.0-3454884612710.1016/j.solener.2007.01.004FilippovA. V.RosnerD. E.Energy transfer between an aerosol particle and gas at high temperature ratios in the Knudsen transition regime20004311271382-s2.0-034263254510.1016/S0017-9310(99)00113-1JohnJ. E. A.1933Boston, Mass, USAAllyn and BaconFujiiT.FujiiM.HondaT.Theoretical and experimental studies of the free convection around a long horizontal thin wire in air2Proceedings of the 7th International Heat Transfer Conference1982Munich, Germany311316ZhangX.FujiwaraS.QiZ.FujiiM.Natural convection effect on transient short-hot-wire method199916129135HsuI. K.PettesM. T.AykolM.ChangC. C.HungW. H.TheissJ.ShiL.CroninS. B.Direct observation of heat dissipation in individual suspended carbon nanotubes using a two-laser technique2011110452-s2.0-8005240257310.1063/1.3627236044328