納米流體強(qiáng)化沸騰換熱研究進(jìn)展

王東民，全曉軍，李金京

(上海交通大學(xué)機(jī)械與動力工程學(xué)院，上海 200240)

1 前 言

隨著工業(yè)技術(shù)的發(fā)展，散熱系統(tǒng)逐漸微型化并且迫切需要散熱能力更強(qiáng)的散熱系統(tǒng)，這需要散熱性能更好的換熱工質(zhì)和更先進(jìn)的冷卻系統(tǒng)。為此，科學(xué)家們發(fā)現(xiàn)了一種新型換熱工質(zhì)——納米流體。納米流體是指將納米尺寸的固體顆粒(通常在100 nm以下)均勻分散在液體中形成的穩(wěn)定膠體分散液，所采用的納米顆粒材料通常為金屬、金屬氧化物或非金屬材料，分散的基本流體通常為水、油、丙酮或乙二醇等。最初納米流體因?yàn)樵诤艿偷臐舛认戮惋@示出比傳統(tǒng)流體具有更高的熱導(dǎo)率而備受關(guān)注[1, 2]。除此之外，相比添加了微米顆?；蚝撩最w粒的懸浮液，納米流體具有更長久的穩(wěn)定性和更高的熱導(dǎo)率，并且在微通道內(nèi)流動所需的壓降也更小，同時也不會腐蝕或堵塞微通道。

納米流體因?yàn)楸葌鹘y(tǒng)流體具有眾多優(yōu)異的物理性質(zhì)而被廣泛應(yīng)用于各種工業(yè)領(lǐng)域。在核反應(yīng)堆散熱領(lǐng)域，如果采用納米流體作為冷卻水，那么在沸騰散熱過程中燃料棒上將覆蓋納米顆粒，使得核反應(yīng)裝置的壓力水冷反應(yīng)堆內(nèi)的燃料棒與冷卻水之間的沸騰換熱的臨界熱流密度提高[3]，從而可以提高核反應(yīng)堆的能量密度、可靠性和安全性。在電子器件散熱領(lǐng)域，因?yàn)榧{米流體具有比傳統(tǒng)流體更高的熱導(dǎo)率以及其強(qiáng)化沸騰換熱的特性，納米流體在微型高功率散熱組件如高速計算微芯片的冷卻中也具有重要應(yīng)用前景[4]。采用納米流體作為熱管的換熱工質(zhì)(如微型槽道熱管、網(wǎng)格吸液芯熱管、振蕩熱管以及兩相虹吸熱管)，可以極大增強(qiáng)換熱速率、減小總體熱阻、提高熱管換熱能力[5]。在汽車領(lǐng)域，因?yàn)榧{米流體的熱導(dǎo)率高于傳統(tǒng)流體和其優(yōu)良的潤滑性能，采用納米流體作為冷卻液代替?zhèn)鹘y(tǒng)發(fā)動機(jī)冷卻液(如：水和乙二醇的混合液)，能極大提高汽車發(fā)動機(jī)的散熱性能[6]。

沸騰換熱作為一種利用換熱工質(zhì)發(fā)生相變來散熱的方式，具有比對流換熱和導(dǎo)熱更高的散熱速率和更大的散熱能力，因此對納米流體強(qiáng)化沸騰換熱的研究成為近年來的研究熱點(diǎn)。下面介紹納米流體的主要制備方法、物理性質(zhì)、納米流體對沸騰換熱的影響，并對未來納米流體強(qiáng)化沸騰換熱的研究進(jìn)行了展望。

2 納米流體的制備及其物理性質(zhì)

2.1 納米流體的制備

納米流體的制備方法通常分為兩類：單步法和兩步法。單步法指納米顆粒的制備和納米流體的制備是聯(lián)合進(jìn)行的，通常指通過物理氣相沉積法或液相化學(xué)法來制備納米顆粒及納米流體。如Eastman等[7]利用流動的乙二醇液體冷卻銅蒸汽來制備銅納米流體，Liu 等[8]利用化學(xué)還原法制備水基銅納米流體等。單步法的優(yōu)點(diǎn)是不需要干燥和存儲納米顆粒，也不需要將固體納米顆粒分散到流體中去，因此避免了這些步驟所帶來的納米顆粒的團(tuán)聚問題，提高了納米流體的穩(wěn)定性。但該方法的缺點(diǎn)是只有蒸汽壓力較低的流體才適合采用單步法來制備納米流體。

兩步法指在第一步中先制備納米顆粒，第二步將納米顆粒分散在基本流體中。如Hong等[9]先通過化學(xué)蒸汽冷凝的方法制備鐵納米顆粒，然后將鐵納米顆粒分散在乙二醇中；Xuan等[10]利用兩步法制備了水基銅納米流體以及油基銅納米流體等。兩步法的優(yōu)點(diǎn)是目前納米顆粒的制備已經(jīng)實(shí)現(xiàn)商業(yè)化，可以很方便地購置。但由于納米顆粒具有很高的表面能，在利用兩步法制備納米流體時，納米顆粒在兩個步驟中都會發(fā)生不同程度的團(tuán)聚[11]。納米顆粒的團(tuán)聚會影響納米流體的換熱性能，甚至?xí)?dǎo)致納米顆粒從納米流體中析出、發(fā)生沉降或者堵塞換熱通道。

為了抑制納米流體中納米顆粒的團(tuán)聚，提高納米流體的穩(wěn)定性。通?？梢圆扇〕曁幚?、改變pH值、向納米流體中添加分散劑、修飾納米顆粒表面等手段來使納米流體保持長期的穩(wěn)定[12]。Li 等[13]發(fā)現(xiàn)在pH=9.5時，向水基銅納米流體中添加分散劑能使納米流體的穩(wěn)定性達(dá)到最佳。Murshed等[14]采用油酸和溴化十六烷基三甲銨(CATB)作為分散劑加入到水基TiO2納米流體中以達(dá)到提高納米流體穩(wěn)定性的目的。Yu等[15]采用等離子體聚合的方法修飾金剛石納米顆粒表面，提高納米顆粒表面的親水性，從而提高納米流體的穩(wěn)定性。楊雪飛等[16]通過在二氧化硅納米顆粒表面嫁接硅烷，利用空間位阻穩(wěn)定性原理來提高二氧化硅納米流體的穩(wěn)定性。

2.2 納米流體的物理性質(zhì)

納米顆粒分散到基本流體中可以極大增強(qiáng)基本流體的熱導(dǎo)率。Choi 等[1]認(rèn)為這是由于納米顆粒的熱導(dǎo)率高于基本流體；Xuan等[17]認(rèn)為除此之外，納米顆粒在納米流體中的布朗運(yùn)動會引起熱分散效應(yīng)，也是導(dǎo)致納米流體熱導(dǎo)率高于傳統(tǒng)流體的原因。通常納米流體熱導(dǎo)率增強(qiáng)的原因主要有：基本流體中的導(dǎo)熱、納米顆粒內(nèi)部的導(dǎo)熱、納米顆粒的布朗運(yùn)動、納米顆粒表面吸附的液體層以及納米顆粒的團(tuán)聚效應(yīng)[18]。

低濃度的納米顆粒分散到基本流體中會引起基本流體粘度的輕微變化，而對納米流體粘度的研究相比納米流體熱導(dǎo)率的研究還很少。在低濃度下，往往利用Einstein公式[19]估算納米流體的粘度，基于Einstein公式發(fā)展的Brinkman 公式[20]則被廣泛應(yīng)用于納米流體的數(shù)值模擬計算時針對納米流體粘度的估算。

Vafaei等[21]分析了碲化鉍納米顆粒的粒徑和濃度對氣-液表面張力的影響，發(fā)現(xiàn)隨著濃度的增大，納米流體氣-液表面張力減小，超過某一濃度后表面張力反而增大。在同一濃度下，2.5 nm比10.4 nm的納米顆粒對應(yīng)的納米流體具有更小的氣-液表面張力。綜上，納米顆粒的濃度和尺度對氣-液表面張力具有重要影響。

3 納米流體對沸騰換熱的影響

3.1 納米流體的沸騰換熱系數(shù)

沸騰換熱系數(shù)是評價一種換熱工質(zhì)在沸騰換熱過程中換熱速率的重要參數(shù)。當(dāng)加熱表面在沸騰換熱過程中以恒定散熱功率向不同換熱工質(zhì)中散熱時，如果加熱表面的散熱溫度越低，則表示對應(yīng)的換熱工質(zhì)沸騰換熱速率越大，其沸騰換熱系數(shù)越高。

目前國內(nèi)外學(xué)者在對納米流體能否增強(qiáng)沸騰換熱系數(shù)這一問題上還存在廣泛爭論。早期Sonltani等[22]認(rèn)為沸騰過程中氣泡的底部和加熱表面之間存在極薄的液體層，而納米顆粒在這層液體層內(nèi)能增強(qiáng)該液層內(nèi)部的導(dǎo)熱，提高加熱表面向沸騰氣泡內(nèi)的熱量傳遞速率，從而使得納米流體的沸騰換熱系數(shù)比純流體的更高。Wen和Ding[23]則認(rèn)為納米流體強(qiáng)化沸騰換熱并不僅僅只受到該種因素的影響。Kim等[24]提出在納米流體沸騰過程中，納米顆粒會沉積在加熱表面，增大加熱表面的親水性，這反而會使得加熱表面上的沸騰換熱系數(shù)降低。Li等[25]也認(rèn)為納米顆粒會在氣泡底部富集(圖1)。Narayan等[26]則認(rèn)為納米流體能否增強(qiáng)沸騰換熱系數(shù)，取決于納米顆粒的平均粒徑和加熱表面粗糙度的相對大小。如果納米顆粒的平均粒徑大于加熱表面的粗糙度，那么沉積在加熱表面上的納米顆粒會增大加熱表面的粗糙度，提高加熱表面的活化成核點(diǎn)密度，增強(qiáng)沸騰換熱系數(shù)(圖2[27]，圖3a[28])。而如果納米顆粒的平均粒徑小于加熱表面的粗糙度，那么沉積在加熱表面上的納米顆粒會填平加熱表面上的空穴，減小加熱表面的粗糙度，從而減小加熱表面的活化成核點(diǎn)密度，弱化沸騰換熱系數(shù)(圖3b[28])。但納米顆粒過度沉積到加熱表面上會導(dǎo)致加熱表面和換熱流體間導(dǎo)熱熱阻增大，使沸騰換熱系數(shù)減小。而且納米顆粒常常會以團(tuán)聚體的形式沉積到加熱表面上[29]。

圖1 隨著氣泡的生長納米顆粒在氣泡底部富集的示意圖[25]Fig.1 Schematic of enrichment and deposition of nanoparticles at the bottom of vapor bubbles as vapor bubbles grow[25]

圖2 納米流體沸騰前(a)和沸騰后(b)沉積納米顆粒的加熱表面的SEM照片[27]Fig.2 SEM images of heating surface before (a) and after (b) boiling of nanofulids with nanoparticles deposition[27]

圖3 納米顆粒沉積在光滑加熱表面(a)和粗糙加熱表面(b)對沸騰換熱造成不同影響的示意圖[28]Fig.3 Schematic of different effects of nanoparticles deposition on smooth heating surface (a) and rough heating surface (b) on boiling heat transfer[28]

最近，Quan等[30]發(fā)現(xiàn)中等親水性的納米顆粒(接觸角為37°左右)沉積在加熱表面上會增大加熱表面的粗糙度，進(jìn)而使加熱表面上的沸騰換熱系數(shù)增大；而強(qiáng)親水性的納米顆粒(接觸角為22°左右)沉積在加熱表面上并不會使得加熱表面粗糙度明顯增大，且會形成較厚的沉積層，增大納米流體和加熱表面之間的熱阻，惡化沸騰換熱系數(shù)(圖4)。

圖4 中等親水性納米顆粒 (a) 和強(qiáng)親水納米顆粒(b) 在加熱表面上沉積的照片[30]Fig.4 Photos of moderately hydrophilic nanoparticles (a) and strongly hydrophilic nanoparticles (b) deposition on heating surfaces[30]

除了沉積在加熱表面上的納米顆粒會影響沸騰換熱系數(shù)外，懸浮于納米流體中的納米顆粒也會通過影響沸騰氣泡對沸騰換熱系數(shù)造成重要影響。Vafaei等[21]就發(fā)現(xiàn)懸浮的納米顆粒會減小納米流體的表面張力，從而使納米流體中產(chǎn)生的沸騰氣泡尺寸更小，增強(qiáng)沸騰換熱系數(shù)。Tu等[31]也發(fā)現(xiàn)在水基Al2O3納米流體沸騰換熱實(shí)驗(yàn)中產(chǎn)生的沸騰氣泡小于純水中產(chǎn)生的沸騰氣泡。Xu等[32]通過開展Al2O3納米流體的流動沸騰實(shí)驗(yàn)，發(fā)現(xiàn)納米流體中產(chǎn)生的沸騰氣泡比純水中的更小，并提出這是因?yàn)榧{米顆粒分散在氣泡底部極薄的液體層內(nèi)，增大了氣泡底部液體層內(nèi)的壓力，促進(jìn)了氣泡的脫離。該壓力即為Wen[33]提出的納米顆粒聚集在氣泡底部引起的結(jié)構(gòu)分離壓力(圖5)。李啟明等[34]利用熱力學(xué)理論分析了納米流體沸騰換熱的成核過程，指出納米顆粒會在界面聚集，降低氣液界面自由能，因此納米流體內(nèi)更容易成核。

Quan等[30]發(fā)現(xiàn)中等親水性的納米顆粒(接觸角為37°左右)會吸附到沸騰氣泡的界面上，抑制氣泡之間的合并，促進(jìn)沸騰換熱系數(shù)的提高；而強(qiáng)親水性的納米顆粒(接觸角為22°左右)不會吸附到沸騰氣泡的界面上，也不會對沸騰氣泡的合并脫離等造成影響，因此不會促進(jìn)沸騰換熱系數(shù)的提高(圖6)。Hu等[35]則發(fā)現(xiàn)隨著熱流密度的增加，石墨烯-乙二醇水納米流體的沸騰換熱系數(shù)增加，其存在納米流體最優(yōu)濃度，使得沸騰換熱系數(shù)最大(圖7)。

圖5 納米顆粒在氣泡底部聚集引起結(jié)構(gòu)分離壓力[33]Fig.5 Concentrated nanoparticles under vapor bubbles resulting structure disjoint pressure[33]

圖6 中等/強(qiáng)親水性納米流體以及純水中沸騰氣泡的圖片[30]Fig.6 Photo of vapor bubbles in strongly/moderately hydrophilic nanoparticles suspension and pure water [30]

圖7 不同石墨烯納米顆粒濃度下石墨烯-乙二醇水納米流體沸騰換熱系數(shù)隨加熱熱流密度的變化曲線[35]Fig.7 Variation curves of boiling heat transfer of graphene-ethylene glycol/water nanofluids with heating heat flux density under different concentrations of graphene nanoparticle[35]

3.2 納米流體沸騰換熱的臨界熱流密度

沸騰換熱的臨界熱流密度是衡量換熱工質(zhì)最大換熱能力的重要參數(shù)。沸騰換熱的臨界熱流密度指當(dāng)散熱表面上的散熱功率不斷提高直至達(dá)到某一臨界值時，加熱表面上迅速產(chǎn)生的大量沸騰氣泡來不及脫離加熱表面，合并成一層氣膜覆蓋在加熱表面上。這時加熱表面與換熱工質(zhì)之間的沸騰換熱系數(shù)急劇下降，加熱表面溫度迅速上升，沸騰換熱狀態(tài)由核態(tài)沸騰變?yōu)槟B(tài)沸騰。因此，換熱工質(zhì)的臨界熱流密度越高，表示該換熱工質(zhì)能夠承受的散熱功率越高。

目前國內(nèi)外學(xué)者的研究結(jié)果一致表明納米流體的臨界熱流密度高于傳統(tǒng)流體。影響納米流體臨界熱流密度的因素主要包括納米顆粒的粒徑、納米顆粒濃度、納米流體表面張力、沉積到加熱表面的納米顆粒沉積層。Moreno等[36]發(fā)現(xiàn)3 nm的納米顆粒對應(yīng)的納米流體的臨界熱流密度相比純流體增大了31%，而隨著納米顆粒粒徑的增大，納米流體的臨界熱流密度相比純流體的增長幅度逐漸減小。You等[37]通過實(shí)驗(yàn)研究發(fā)現(xiàn)隨著水基Al2O3納米流體的濃度從0.001增長到0.025 g/L，納米流體的臨界熱流密度不斷增大，最大增幅可達(dá)200%。Jeong等[38]發(fā)現(xiàn)隨著磁性納米流體濃度從1×10-6增大到1×10-4，納米流體的臨界熱流密度增幅從160%增至240%。納米流體的表面張力也是納米流體臨界熱流密度增長的重要影響因素。Kumar和Milanova[39]通過實(shí)驗(yàn)研究發(fā)現(xiàn)納米流體和基本流體表面張力的不同是延遲碳納米管納米流體臨界熱流密度到來的原因。

沉積到加熱表面上的納米顆粒沉積層主要通過改變加熱表面的親水性和毛細(xì)作用影響納米流體的臨界熱流密度。Kim等[40]指出納米顆粒沉積到加熱面上后增大了加熱面的親水性，從而使得納米流體的臨界熱流密度比純流體的更高。Zhang和Jacobiy[41]也得到了類似的結(jié)論(圖8)。Park等[42]研究了水基Al2O3納米流體和石墨烯納米流體的臨界熱流密度，并指出與臨界熱流密度相關(guān)的臨界波長的改變是臨界熱流密度增大的主要原因。僅僅靠加熱表面親水性的增強(qiáng)還解釋不了納米流體臨界熱流密度比純流體大大增強(qiáng)的現(xiàn)象。Kim 等[43]指出納米顆粒在加熱面上的沉積層具有毛細(xì)吸液作用，能將更多的液體吸引到氣泡底部的干斑區(qū)域，抑制氣泡底部干斑的生長，延緩臨界熱流密度的到來(圖9)。近年來，Quan等[30]提出納米顆粒的接觸角也能對納米流體的臨界熱流密度產(chǎn)生重要影響。中等親水性的納米顆粒(接觸角為37°左右)能吸附在氣泡界面上，抑制氣泡合并成覆蓋在加熱表面上的氣膜，延緩納米流體膜態(tài)沸騰換熱的到來；而強(qiáng)親水性的納米顆粒(接觸角為22°左右)卻因?yàn)椴荒芪降綒馀萁缑嫔?，不能抑制氣泡合并，不能延緩納米流體膜態(tài)沸騰換熱的到來。

圖8 Al2O3納米顆粒在加熱表面上沸騰沉積的SEM照片及其與加熱表面接觸角[41]Fig.8 SEM images of Al2O3 nanoparticles deposition on heating surfaces after boiling heat transfer and the contact angle between Al2O3 nanoparticles and the heating surface[41]

圖9 納米顆粒沉積層的毛細(xì)作用對氣泡底部干斑的影響[43]Fig.9 Capillary effects of nanoparticles deposition on dry spot under a vapor bubble[43]

3.3 納米流體沸騰換熱的理論和數(shù)值研究

由于納米流體沸騰換熱的復(fù)雜性，目前對納米流體強(qiáng)化換熱的理論和數(shù)值研究主要集中于納米流體的單相對流或?qū)岬难芯?。對納米流體強(qiáng)化沸騰換熱的理論和數(shù)值模擬研究還很少，主要是通過修改傳統(tǒng)流體的沸騰換熱理論模型或數(shù)值模型中流體的物性參數(shù)來實(shí)現(xiàn)。如Behroyan等[44]、Mohammed等[45]以及Salehi和Hormozi等[46]利用計算流體力學(xué)方法(CFD)將納米流體處理成純流體或者混合流體來研究納米流體的沸騰換熱性能。Abedini等[47]采用有限體積法，將納米流體處理成純流體研究了水基Al2O3流體的沸騰換熱性能。Wang和Wu[48]采用移動顆粒半隱式方法研究了納米顆粒濃度對氣泡脫離特性的影響。

4 結(jié) 語

納米流體作為一種新型換熱工質(zhì)，與傳統(tǒng)流體相比除了易于制備、高熱導(dǎo)和良好的穩(wěn)定性之外，也能對沸騰換熱起到強(qiáng)化作用，特別是納米流體在低濃度下就能大幅提高沸騰換熱的臨界熱流密度。納米流體對沸騰換熱的影響主要是通過流體中的納米顆粒和沉積到加熱表面上的納米顆粒沉積層兩方面產(chǎn)生的。流體中的納米顆粒影響沸騰氣泡的脫離特性、氣泡的形態(tài)以及氣泡之間的合并，從而影響沸騰氣泡從加熱表面帶走的熱量和臨界熱流密度下加熱表面氣膜的形成，最終會影響沸騰換熱的換熱系數(shù)和臨界熱流密度。而沉積在加熱表面上的納米顆粒會改變加熱表面的潤濕性、粗糙度和毛細(xì)作用，從而影響加熱表面上的活化成核點(diǎn)密度、加熱表面與流體之間的導(dǎo)熱熱阻、氣泡脫離特性以及氣泡底部干斑的生長，最終影響沸騰換熱系數(shù)和臨界熱流密度。

但目前納米流體沸騰換熱的實(shí)驗(yàn)研究結(jié)果可重復(fù)性不高，導(dǎo)致不同研究者對同一納米流體研究得到的實(shí)驗(yàn)結(jié)果不盡相同，甚至互相矛盾，這主要是以下原因引起的：① 由于納米流體制備方法的不同或者即使采用同一種制備方法，往往會導(dǎo)致納米顆粒在納米流體內(nèi)的團(tuán)聚程度不一致，納米顆粒在納米流體內(nèi)并不是單分散的。② 為了保持納米流體的穩(wěn)定性，常常需要改變流體的pH值、或修飾納米顆粒表面或向納米流體中添加分散劑，這為實(shí)驗(yàn)研究納米流體對沸騰換熱的影響帶來了其他因素的影響。③ 即使同一納米流體在沸騰換熱過程中，納米顆粒會隨機(jī)沉積在加熱面上，導(dǎo)致不同研究者在實(shí)驗(yàn)過程中得到的納米顆粒沸騰沉積層也不盡一致，這增加了實(shí)驗(yàn)的不確定性，從而使得到的納米流體對沸騰換熱影響的結(jié)論互相矛盾。解決在實(shí)驗(yàn)研究納米流體對沸騰換熱的影響中出現(xiàn)的以上問題將能更清楚地揭示納米流體影響沸騰換熱的影響機(jī)理。另外，在實(shí)際應(yīng)用方面，納米流體強(qiáng)化沸騰換熱的持久性較差。這是因?yàn)樵诜序v換熱過程中，納米顆粒不斷沉積到加熱表面上，使流體中的納米顆粒濃度不斷減小，同時導(dǎo)致加熱表面的導(dǎo)熱熱阻不斷增大。

目前對納米流體強(qiáng)化沸騰換熱的研究手段大部分為實(shí)驗(yàn)研究。建立系統(tǒng)有效的納米流體沸騰換熱理論模型和納米流體沸騰換熱的數(shù)值模擬研究亟待開展，以指導(dǎo)納米流體沸騰換熱的實(shí)驗(yàn)研究和納米流體強(qiáng)化沸騰換熱的工程應(yīng)用。