多壁碳納米管-水納米流體導(dǎo)熱機(jī)理及重力熱管實(shí)驗(yàn)研究

吳晗，楊峻

?

多壁碳納米管-水納米流體導(dǎo)熱機(jī)理及重力熱管實(shí)驗(yàn)研究

吳晗，楊峻

（南京工業(yè)大學(xué)機(jī)械與動(dòng)力工程學(xué)院，江蘇南京211816）

在考慮了納米層的情況下對原有的Xue導(dǎo)熱模型進(jìn)行改進(jìn)，研究多壁碳納米管-水納米流體的熱傳導(dǎo)性能，推導(dǎo)出納米流體導(dǎo)熱公式并將納米流體運(yùn)用到碳鋼重力熱管中，在不同質(zhì)量分?jǐn)?shù)下對單管傳熱進(jìn)行實(shí)驗(yàn)研究。結(jié)果表明納米層的存在提高了有效熱導(dǎo)率；在相同條件下質(zhì)量分?jǐn)?shù)2%的多壁碳納米管-水納米流體重力熱管傳熱性能最好，傳熱系數(shù)比普通水重力熱管最大提高了40%。

納米粒子；納米流體；熱傳導(dǎo)；重力熱管；傳熱

引 言

納米流體[1]作為一種新型高導(dǎo)熱介質(zhì)而得到廣泛應(yīng)用。納米流體的較高的導(dǎo)熱能力引起研究者的關(guān)注并從不同角度探究導(dǎo)熱機(jī)理，截至目前研究角度主要為懸浮液的特性[2-4]、納米顆粒布朗運(yùn)動(dòng)[5-7]、納米顆粒的團(tuán)聚[8-9]、納米層[10-12]等，由于納米流體的導(dǎo)熱機(jī)理復(fù)雜，對其研究者眾多，到目前還沒有確切結(jié)論。近些年將納米流體運(yùn)用于熱管傳熱也成為研究方向之一，其中重力熱管由于制造方便、成本低廉、傳熱良好等優(yōu)點(diǎn)被廣泛應(yīng)用于余熱回收[13]、空調(diào)制冷[14]、電子散熱[15]等傳熱設(shè)備中，且通過在重力熱管中添加納米顆粒，有顯著發(fā)現(xiàn)。郭廣亮等[16]使用多壁碳納米管-水納米流體代替水后，熱管臨界熱通量比原來最大提高了120%。宮玉英等[17]使用SiO2-水納米流體改進(jìn)水重力熱管，發(fā)現(xiàn)改進(jìn)后的熱管傳熱性能比普通水重力熱管提高了1.35～1.70倍。Noie等[18]使用Al2O3-水納米流體在體積濃度為1%～3%時(shí)研究重力熱管傳熱性能，發(fā)現(xiàn)納米流體重力熱管相對于普通水重力熱管管壁溫度分布更加均勻，傳熱效率提高了14.7%。納米顆粒的加入一定程度上提高了熱管傳熱性能，效果可觀。然而目前將納米流體應(yīng)用于重力熱管研究報(bào)道相對較少，此方面仍處于初級階段[19]。

納米顆粒一般有金屬粒子、氧化物粒子、碳化物粒子、多壁碳納米管等，其中多壁碳納米管由于優(yōu)越的力學(xué)、電學(xué)和化學(xué)性能及較高的導(dǎo)熱能力而比其他納米顆粒具有更大的強(qiáng)化傳熱潛力[20-21]。

本文選取多壁碳納米管為研究對象，水為基液，模型簡化之后從納米層的角度研究納米流體導(dǎo)熱機(jī)理，并運(yùn)用到重力熱管中進(jìn)行單管傳熱性能實(shí)驗(yàn)研究，掌握換熱規(guī)律，為余熱回收及其應(yīng)用提供技術(shù)支持。

1 納米流體導(dǎo)熱研究

1.1 模型建立和公式推導(dǎo)

考慮到多壁碳納米管的形狀和納米層的存在，為研究方便將多壁碳納米管簡化為半徑p、高度為的細(xì)長圓柱，納米層為包裹在圓柱體外的同心圓柱，納米層厚度為（圖1），忽略各壁面間影響。

根據(jù)Xue[22]導(dǎo)熱模型，在不考慮納米層時(shí)納米流體有效熱導(dǎo)率為

由于納米層的存在引起體積分?jǐn)?shù)和固相熱導(dǎo)率的變化，使用等效體積分?jǐn)?shù)e和等效固相熱導(dǎo)率pe分別代替和p。

等效體積分?jǐn)?shù)e的計(jì)算

式中，為單位體積的顆粒數(shù)目；=1+/p。

為計(jì)算等效固相熱導(dǎo)率pe，作出如下假設(shè)：

①納米流體的濃度比較低，顆粒相互之間的作用效果忽略不計(jì)；

②溫度場連續(xù)；

③溫度場滿足拉普拉斯方程。

為方便求解建立二維平面的柱坐標(biāo)系，平面任意一點(diǎn)的坐標(biāo)為(,)，假設(shè)溫度場方向如圖2所示。

平面上點(diǎn)溫度的拉普拉斯方程為

利用(,)=(, ?)，(,)=(,+2π)得出該拉普拉斯方程溫度方程的通解為：

當(dāng)0≤

p=cos(4)

當(dāng)p≤

其中，、、為常數(shù)，在納米層和納米顆粒交界處溫度和熱通量分別對應(yīng)相等，則

(6)

聯(lián)立式(4)～式(7)求得

(8)

將式(8)、式(9)代入到式(4)、式(5)中求出溫度沿方向上的溫度梯度

(10)

式(11)中取同心圓環(huán)平均直徑，即

(12)

由熱導(dǎo)率定義得出

式中，〈q〉、〈〉分別為平均熱通量和平均溫度梯度，φ為納米層的體積分?jǐn)?shù)，由式(2)推出φ=(2?1)，將式(10)、式(12)代入式(13)推出

(14)

其中=1，1=+1，2=?1。

對于納米層熱導(dǎo)率的計(jì)算，假設(shè)納米層隨著半徑呈線性變化，利用等效熱阻法求出k

當(dāng)考慮納米層存在有效熱導(dǎo)率為

(16)

1.2 分析討論

經(jīng)研究發(fā)現(xiàn)納米層厚度約為1 nm[23-25]，在此基礎(chǔ)上研究納米層對有效熱導(dǎo)率的影響。

影響多壁碳納米管導(dǎo)熱的因素較多，研究計(jì)算方法也有很多種[26-29]，本研究在室溫=300 K情況下進(jìn)行討論，根據(jù)研究結(jié)果[27]選取p=3000 W·m-1·K-1，f=0.6 W·m?1·K?1，并做出在不同半徑和體積分?jǐn)?shù)下有效熱導(dǎo)率和基液的比值曲線，并和原Xue模型結(jié)果對比（圖3）。

由圖3可知改進(jìn)后的模型增加了有效熱導(dǎo)率，有效熱導(dǎo)率和體積分?jǐn)?shù)大致呈線性關(guān)系，多壁碳納米管半徑越小有效熱導(dǎo)率越大，圖中有效熱導(dǎo)率比Xue模型增加了0.3%～5%，由此可知納米層的存在對納米流體導(dǎo)熱性能提高有重要影響。

2 重力熱管傳熱實(shí)驗(yàn)

將多壁碳納米管-水納米流體運(yùn)用于重力熱管實(shí)驗(yàn)，在一定充液量下研究不同的質(zhì)量分?jǐn)?shù)的多壁碳納米管-水納米流體重力熱管傳熱性能并和普通水重力熱管對比，得出最佳濃度配比。

2.1 實(shí)驗(yàn)裝置和方法

實(shí)驗(yàn)裝置如圖4所示，重力熱管采用碳鋼制成，外徑25 mm，壁厚2.5 mm，蒸發(fā)段、絕熱段、冷凝段分別設(shè)定為500、150、600 mm，本實(shí)驗(yàn)采用強(qiáng)制對流冷卻并在冷凝段焊接有長500 mm，外徑45 mm，壁厚2.5 mm的水套，通過冷卻水在水套內(nèi)的流動(dòng)來冷卻熱管。溫度測量采用T型熱電偶，其中蒸發(fā)段布置3對，絕熱段布置2對，冷凝段布置2 對。熱電偶的測量點(diǎn)分布如圖5所示。

1—three-phase alternating current; 2—distribution box; 3—electric heating furnace; 4—heating wire; 5—insulation; 6—insulation cotton; 7—insulation brick; 8—cooling jacket; 9—exhaust valve; 10,14—sink; 11—rotor flow meter; 12—water pump; 13—thermometer; 15—thermocouple; 16—acquisition system

納米工質(zhì)使用直接購置的多壁碳納米管，具體規(guī)格見表1，分別用兩步法配制質(zhì)量分?jǐn)?shù)為0.5%、1%、1.5%、2%、2.5%的納米流體，并用超聲波破碎儀振蕩3 h，熱管充液量為120 ml并采用加熱排氣法排成真空制成重力熱管，蒸發(fā)段用電爐加熱，設(shè)定爐膛溫度300℃，進(jìn)出的冷卻水流量為0.2 L·min?1，為防止冷凝段空氣對流引起的散熱采用保溫棉將冷凝段包裹，同時(shí)為增強(qiáng)對比將一定質(zhì)量分?jǐn)?shù)的多壁碳納米管-水納米流體重力熱管和相同規(guī)格與充液量的水重力熱管同時(shí)放入爐膛里進(jìn)行實(shí)驗(yàn)，待系統(tǒng)穩(wěn)定后分別記錄進(jìn)出水的溫度和熱管管壁熱電偶讀數(shù)。

表1 多壁碳納米管規(guī)格

2.2 實(shí)驗(yàn)結(jié)果分析

利用傳熱公式得出單管傳熱系數(shù)為

其中

=(o-i) (18)

=p(20)

式中，為水的比熱容；為水的質(zhì)量流量；io為冷卻水進(jìn)出溫度；、為蒸發(fā)段和冷凝段的平均溫度；為熱管外徑；為熱管長度。并作出不同質(zhì)量分?jǐn)?shù)下重力熱管壁溫分布、傳輸功率、單管傳熱系數(shù)曲線，如圖6～圖8所示。

由圖6得出在相同的條件下，多壁碳納米管-水納米流體重力熱管和水重力熱管的絕熱段和冷凝段的溫度基本相等，在蒸發(fā)段質(zhì)量分?jǐn)?shù)0.5%、1%、1.5%、2%的多壁碳納米管-水重力熱管平均壁溫比水重力熱管分別降低了3、8、10、12℃，質(zhì)量分?jǐn)?shù)2.5%的平均壁溫比水重力熱管升高了2℃；圖7中質(zhì)量分?jǐn)?shù)0.5%、1%、1.5%、2%多壁碳納米管-水納米流體重力熱管傳輸功率相對水重力熱管分別提高了3%～6%、6%～10%、10%～13%、12%～17%，質(zhì)量分?jǐn)?shù)2.5%的降低了2%～5.5%；圖8中質(zhì)量分?jǐn)?shù)0.5%、1%、1.5%、2%多壁碳納米管-水納米流體重力熱管傳熱系數(shù)比水重力熱管分別提高了5%～11%、20%～25%、30%～35%、35%～40%，質(zhì)量分?jǐn)?shù)2.5%的傳熱系數(shù)比水重力熱管降低6%～12%。本實(shí)驗(yàn)單管傳熱性能=2%>=1.5%>=1%>=0.5%> 水重力熱管>=2.5%。

本實(shí)驗(yàn)對一定范圍濃度的納米流體強(qiáng)化重力熱管傳熱進(jìn)行分析，當(dāng)蒸發(fā)段加熱時(shí)，熱管工質(zhì)在管內(nèi)液池處于核態(tài)沸騰，大量的多壁碳納米顆粒吸附在管壁上。當(dāng)納米顆粒的粒徑大于液體均質(zhì)成核的臨界半徑時(shí)[30]，納米顆粒促進(jìn)了流體內(nèi)的氣泡成核，增加了氣化并產(chǎn)生了密度比水大很多的氣泡，又由于布朗運(yùn)動(dòng)增加了擾動(dòng)并迅速帶走氣泡，降低了液膜熱阻從而強(qiáng)化傳熱；當(dāng)濃度超過一定范圍時(shí)，納米顆粒發(fā)生團(tuán)聚降低了擾動(dòng)程度，過高的濃度又使得流動(dòng)阻力變大增加了蒸發(fā)段熱阻，從而出現(xiàn)傳熱系數(shù)迅速降低的現(xiàn)象。

3 結(jié) 論

（1）在考慮了納米層情況下，改進(jìn)后的模型增加了多壁碳納米管-水納米流體的熱導(dǎo)率；一定體積分?jǐn)?shù)下，多壁碳納米管半徑越小，有效熱導(dǎo)率越大。

（2）不同質(zhì)量分?jǐn)?shù)下的單管實(shí)驗(yàn)中，多壁碳納米管質(zhì)量分?jǐn)?shù)0.5%、1%、1.5%、2%時(shí)的重力熱管相對水重力熱管強(qiáng)化了傳熱，質(zhì)量分?jǐn)?shù)2.5%時(shí)傳熱性能變差。

（3）質(zhì)量分?jǐn)?shù)2%的多壁碳納米管-水納米流體重力熱管傳熱性能最好，和普通水重力熱管相比蒸發(fā)段平均壁溫降低了12℃；傳輸功率提高了12%～17%；單管傳熱系數(shù)提高了35%～40%。

符 號 說 明

D——熱管外徑，mm H——多壁碳納米管的高度，mm h——熱管傳熱系數(shù)，W·m?2·K?1 k——熱導(dǎo)率，W·m?1·K?1 L——熱管長度，mm Q——熱管傳輸功率，W r——多壁碳納米管半徑，nm T——溫度，K Tw——熱管管壁溫度，℃ t——時(shí)間，h δ——納米層厚度，nm φ——納米流體體積分?jǐn)?shù) ω——納米流體質(zhì)量分?jǐn)?shù) 下角標(biāo) e——等效體積 eff——有效 f——基液 l——納米層 p——納米顆粒 pe——等效固相

References

[1] CHOIS S. Enhancing Thermal Conductivity of Fluids with Nanoparticles[M]. New York: ASME, 1995: 99.

[2] GUPTE S K, ADVANI S G. Role of micro-convection due to non-acne motion of particles in a mono-disperse suspension[J]. International Journal of Heat and Mass Transfer, 1995, 38(16): 2945-2958.

[3] DAVIS R H. The effective thermal conductivity of a composite material with spherical inclusions[J]. International Journal of Thermophysics, 1986, 7(3): 609-620.

[4] LU S, LIN H. Effective conductivity of composites containing aligned spherical inclusions of finite conductivity[J]. Journal of Applied Physics, 1996, 79(9): 6761-6769.

[5] EINSTEIN A. Investigations on the Theory of the Brownian Movement[M]. New York: Dover, 1956.

[6] 宣益民. 納米流體能量傳遞理論與應(yīng)用[J]. 中國科學(xué): 技術(shù)科學(xué), 44(3): 269-279.XUAN Y M. An overview on nanofluids and applications[J]. Sci. Sin. Tech., 44(3): 269-279.

[7] WANG J J, ZHENG R T, GAO J W,. Heat conduction mechanisms in nanofluids and suspensions[J]. Nano Today, 2012, 7(2): 124-136.

[8] 宣益民, 胡衛(wèi)峰, 李強(qiáng). 納米流體的聚集結(jié)構(gòu)和熱導(dǎo)率模擬[J]. 工程熱物理學(xué)報(bào), 2002, 23(2): 206-208. XUAN Y M, HU W F, LI Q. Simulations of structure and thermal conductivity of nanofluids[J]. Journal of Engineering Thermophysics, 2002, 23(2): 206-208.

[9] 王補(bǔ)宣, 盛文彥. 納米流體熱導(dǎo)率的團(tuán)簇宏觀分析模型[J]. 自然科學(xué)進(jìn)展, 2007, 17(7): 984-988. WANG B X, SHENG W Y. Cluster macro analysis model of nanofluids conductivity[J]. Progress in Natural Science, 2007, 17(7): 984-988.

[10] YU W, CHOI S. The role of interfacial layers in the enhanced thermal conductivity of nanofluids: a renovated maxwell model[J]. Journal of Nanoparticle Research, 2003, 6(4): 355-361.

[11] XIE H, FU M, ZHANG X. Effect of interfacial nanolayer on effective thermal conductivity of naooparticle-fluid mixture[J]. International Journal of Heat and Mass Transfer, 2005, 48(14): 2926-2932.

[12] APLIPOUR R, JAFARI A, MIZAEE M,. The direct effect of interfacial nanolayers on thermal conductivity of nanofluids[J]. Heat Mass Transfer, 2014, 50(12): 1727-1735.

[13] 徐偉, 陳思嘉, 何燕, 等. 熱管技術(shù)在余熱回收中的應(yīng)用研究進(jìn)展[J]. 廣東化工, 2007, 34(2): 40-42. XU W, CHEN S J, HE Y,. Advances in applications research of heat pipe technique in waste heat recovery[J]. Guangdong Chemical Industry, 2007, 34(2): 40-42.

[14] 王愛輝, 羅高喬, 汪韓送. 重力式熱管空調(diào)機(jī)組運(yùn)行特性試驗(yàn)研究[J]. 制冷技術(shù), 2013, 33(2): 14-16. WANG A H, LUO G Q, WANG H S. Experimental investigation of operating characteristic of two-phase closed thermosyphon type air conditioning unit[J]. Refrigeration Technology, 2013, 33(2): 14-16.

[15] 趙安林. 重力熱管散熱器實(shí)驗(yàn)研究[D]. 重慶: 重慶大學(xué), 2002. ZHAO A L. Experimental study on gravity heat pipe radiator[D]. Chongqing: Chongqing University, 2002.

[16] 郭廣亮, 劉振華. 碳納米管懸浮液強(qiáng)化小型重力型熱管換熱特性[J]. 化工學(xué)報(bào), 2007, 58(12): 3006-3010. GUO G L, LIU Z H. Heat transfer enhancement of small thermosyphon using carbon nanotube suspensions[J]. Journal of Chemical Industry and Engineering (China), 2007, 58(12): 3006-3010.

[17] 宮玉英, 朱保杰, 劉宗明, 等. SiO2-水納米流體熱管傳熱性能的研究[J]. 化工機(jī)械, 2013, 40(3): 302-305.GONG Y Y, ZHU B J, LIU Z M,. Study on heat transfer performance of SiO2-DW nanofluids heat pipe[J]. Chemical Machinery, 2013, 40(3): 302-305.

[18] NOIE S, HERIS S Z, KAHANI M,Heat transfer enhancement using Al2O3-DW nanofluids in a two-phase closed thermosyphon[J]. International Journal of Heat and Fluid Flow, 2009, 30(4): 700-705.

[19] 周根明, 周少華, 趙忠超, 等. 納米流體重力熱管啟動(dòng)性能的試驗(yàn)研究[J]. 江蘇科技大學(xué)學(xué)報(bào), 2013, 27(4): 376-380. ZHOU G M, ZHOU S H, ZHAO Z C,. Experimental study of start-up gravity heat pipe filled with nanofluids[J]. Journal of Jiangsu University of Science and Technology, 2013, 27(4): 376-380.

[20] 周慶祥, 肖軍平, 汪衛(wèi)東, 等. 碳納米管應(yīng)用研究進(jìn)展[J]. 化工進(jìn)展, 2006, 25(7): 750-754.ZHOU Q X, XIAO J P, WANG W D,. Progress of the application of carbon nanotubes[J]. Chemical Industry and Engineering Progress, 2006, 25(7): 750-754.

[21] 薛懷生, 樊建人, 胡亞才, 等. 碳納米管懸浮液在重力熱管中的沸騰特性[J]. 化工學(xué)報(bào), 2006, 57(11): 2562-2567.XUE H S, FAN J R, HU Y C,. Boiling characteristics of carbon nanotube suspension in gravity-assisted thermosyphon[J]. Journal of Chemical Industry and Engineering (China), 2006, 57(11): 2562-2567.

[22] XUE Q Z. Model for thermal conductivity of carbon nanotube-based composites[J]. Physica B, 2005, 368(4): 302-307.

[23] HASHIMOTO T, FUJIMURA M, KAWAI H. Domain-boundary structure of styrene-isoprene block copolymer films cast from solutions[J]. Macromolecules, 1974, 13(3): 660-669.

[24] YU C, RICHTER A, DATTA A,. Molecular layering in a liquid on a solid substrate: an X-ray reflectivity study[J]. Physica B, 2000, 283(1/2/3): 27-31.

[25] XUE L, KEBLINSKI P, PHILLPOT S,. Effect of liquid layering at the liquid-solid interface on thermaltransport[J]. Int. J. Heat Mass Transf., 2004, 47(19/20): 4277-4284.

[26] CAO J X, YAN X H, XIAO Y,. Electronic transport properties of metallic carbon nanotubes[J]. Phys. Rev. B, 2003, 12(12):1440-1444.

[27] CHE J W, CAGIN T, GODDARD W. Thermal conductivity of carbon nanotubes[J]. Nanotechnology, 2000, 11(2): 2083-2085.

[28] ZHANG W, ZHU Z, WANG F,. Chirality dependence of the thermal conductivity of carbon nanotubes[J]. Nanotechnology, 2004, 15(8): 936-939.

[29] 李慶威. 碳納米管熱傳導(dǎo)研究[D]. 北京: 清華大學(xué), 2011. LI Q W. Study on heat conduction of carbon nanotubes[D]. Beijing: Tsinghua University, 2011.

[30] 向軍, 李菊香. 納米懸浮液熱虹吸管的傳熱性能試驗(yàn)研究[J]. 熱能動(dòng)力工程, 2010, 25(2): 190-196. XIANG J, LI J X. Experimental study on heat transfer performance of thermosyphon[J]. Journal of Engineering for Thermal Energy and Power, 2010, 25(2): 190-196.

Thermal conduction mechanism of multi-walled carbon nanotubes-deionized water nanofluids and experimental research in gravity heat pipe

WU Han, YANG Jun

(College of Mechanical and Power Engineering,Nanjing Tech University,Nanjing211816, Jiangsu,China)

An improved model of thermal conduction of multi-walled carbon nanotubes– deionized water(MWCNTs-DW) nanofluids is proposed based on the Xue model by considering interfacial nanolayer, analyzingheat conduction, derivating thermal conductivity of MWCNTs-DW nanofluids and applyingMWCNTs-DW nanofluids to the research of heat transfer in gravity heat pipe. Different mass fraction of nanofluids gravity heat pipe are prepared and investigated using the nanofluid as the work liquid and the carbon steel pipe as material. The result shows that the interface nanolayer increases the effective thermal conductivity and the mass fraction of 2% is best performance under the same conditions, the heat transfer coefficient of nanofluids gravity heat pipe increases by 40% to the max than that of the base fluid water.

nanoparticles; nanofluids; heat conduction; gravity heat pipe; heat transfer

10.11949/j.issn.0438-1157.20161492

TK 124

A

0438—1157（2017）06—2315—06

楊峻。

吳晗（1992—），男，碩士研究生。

2016-10-23收到初稿，2017-02-25收到修改稿。

2016-10-23.

YANG Jun, yangjun@njtech.edu.cn