高偉龍, 葉 芳, 郭 航, 趙建福, 陳 浩

(1.北京工業(yè)大學(xué) 環(huán)境與生命學(xué)部能源與動(dòng)力工程學(xué)院,北京100124;2.北京工業(yè)大學(xué)傳熱強(qiáng)化與過(guò)程節(jié)能教育部重點(diǎn)實(shí)驗(yàn)室及傳熱與能源利用北京市重點(diǎn)實(shí)驗(yàn)室,北京100124;3.中國(guó)科學(xué)院力學(xué)研究所微重力重點(diǎn)實(shí)驗(yàn)室,北京100190;4.中國(guó)科學(xué)院大學(xué)工程科學(xué)學(xué)院,北京100049;5.中國(guó)科學(xué)院力學(xué)研究所高溫氣體動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京100190)

1 概述

隨著科學(xué)研究的深入和制造技術(shù)水平的提高,各領(lǐng)域的設(shè)備儀器都在逐漸小型化、集成化、緊湊化,功率密度不斷升高[1]。散熱問(wèn)題成為制約許多領(lǐng)域發(fā)展的瓶頸,使設(shè)備工作穩(wěn)定性、可靠性、壽命降低[2]。因此,高效散熱方式對(duì)大功率密度設(shè)備的散熱十分必要。常見(jiàn)的散熱方式有自然風(fēng)冷、強(qiáng)制風(fēng)冷、自然水冷、強(qiáng)制水冷、相變傳熱等,與其他散熱方式相比,相變傳熱利用傳熱工質(zhì)的顯熱和相變潛熱,傳熱能力更強(qiáng)。近年來(lái),相變傳熱一直被學(xué)者們廣泛研究,也被廣泛應(yīng)用于核電站、火力發(fā)電廠、電子元件、計(jì)算機(jī)數(shù)據(jù)中心、電動(dòng)汽車(chē)等領(lǐng)域[3-8]。

池沸騰散熱方式結(jié)構(gòu)簡(jiǎn)單、不需要復(fù)雜的機(jī)械設(shè)備,其中的核態(tài)沸騰能夠在較低的溫度梯度下散發(fā)大量熱量,因此成為最有效的散熱方式之一。核態(tài)沸騰中,液體工質(zhì)從受熱面吸熱,發(fā)生相變產(chǎn)生氣泡,熱量從受熱面排出。氣泡在受熱面生長(zhǎng)及脫離過(guò)程均受到不同力的作用,氣泡形成、生長(zhǎng)、聚并、脫離過(guò)程對(duì)傳熱有很大影響[9]。為進(jìn)一步理解核態(tài)沸騰的復(fù)雜過(guò)程,必須對(duì)核態(tài)沸騰中氣泡動(dòng)力學(xué)參數(shù)的影響因素進(jìn)行深入研究。

本文綜述影響因素對(duì)核態(tài)沸騰中氣泡動(dòng)力學(xué)參數(shù)的影響,整理氣泡動(dòng)力學(xué)參數(shù)的經(jīng)驗(yàn)關(guān)聯(lián)式。

2 氣泡動(dòng)力學(xué)參數(shù)及影響因素

2.1 氣泡動(dòng)力學(xué)參數(shù)

氣泡動(dòng)力學(xué)參數(shù)主要有:氣泡脫離直徑(半徑)、氣泡脫離頻率、氣泡等待時(shí)間、氣泡生長(zhǎng)時(shí)間、氣泡生長(zhǎng)周期、活性成核點(diǎn)密度。

氣泡脫離直徑(半徑)為氣泡脫離受熱面時(shí)的最終直徑(半徑),由于氣泡形狀不是規(guī)則的球形,因此一般用脫離時(shí)刻的氣泡等效直徑(半徑)來(lái)代替最終直徑(半徑),氣泡等效直徑(半徑)指與氣泡體積同等的球體積的球直徑(半徑)。氣泡脫離頻率為氣泡生長(zhǎng)周期的倒數(shù),氣泡生長(zhǎng)周期為氣泡生成等待時(shí)間(簡(jiǎn)稱(chēng)氣泡等待時(shí)間)和氣泡生成后生長(zhǎng)時(shí)間(簡(jiǎn)稱(chēng)氣泡生長(zhǎng)時(shí)間)之和。氣泡等待時(shí)間為從上一個(gè)氣泡脫離到下一個(gè)氣泡生成所需要的時(shí)間,氣泡生長(zhǎng)時(shí)間為從氣泡生成到氣泡脫離所需要的時(shí)間?；钚猿珊它c(diǎn)密度為一定面積受熱面上可以生成氣泡的成核點(diǎn)數(shù)量。

2.2 影響因素

核態(tài)沸騰中氣泡生長(zhǎng)、聚并、脫離過(guò)程均復(fù)雜多變,氣泡動(dòng)力學(xué)參數(shù)受很多因素影響。根據(jù)影響因素來(lái)源,可將影響因素分為內(nèi)在、外在兩類(lèi)。其中,內(nèi)在因素主要包括壓力、微重力、受熱面熱流密度、受熱面過(guò)熱度與工質(zhì)過(guò)冷度、受熱面傾角、受熱面表面粗糙度、接觸角、汽化核心間距等。外在因素主要為外加場(chǎng),包括電場(chǎng)、磁場(chǎng)、超聲波等。

① 壓力

一些學(xué)者研究了壓力對(duì)核態(tài)沸騰中氣泡動(dòng)力學(xué)參數(shù)的影響。Sakashita[10]研究了壓力為2.23、3.17、4.47 MPa時(shí),在水平受熱面的水池沸騰中,隨著壓力增加,氣泡生長(zhǎng)速率減小,氣泡脫離半徑減小。Hutter等人[11]用實(shí)驗(yàn)研究了FC-72在壓力分別為0.152、0.127、0.101 MPa時(shí),壓力對(duì)氣泡脫離頻率的影響,結(jié)果表明,氣泡脫離頻率隨壓力增大略微下降。Miglani等人[12]研究發(fā)現(xiàn),R134a池沸騰中氣泡脫離頻率隨壓力增大而增大,氣泡脫離直徑隨壓力增大而減小。

Surtaev等人[13]研究了壓力8.8～103.0 kPa對(duì)核態(tài)沸騰的影響,隨著壓力降低,活性成核點(diǎn)密度降低,而氣泡生長(zhǎng)速率和氣泡脫離直徑均增大。Michaie等人[14]研究了從4.2 kPa到當(dāng)?shù)卮髿鈮旱牟煌瑝毫?duì)氣泡動(dòng)力學(xué)參數(shù)的影響,結(jié)果表明,隨著壓力降低,氣泡脫離直徑和氣泡脫離頻率均增加。壓力不同時(shí),氣泡分別呈現(xiàn)蘑菇形或蒸汽柱、扁球形、球形。

此外,Lamas等人[15]研究了壓力為0.01～1.00 MPa時(shí),水、R134a、氨水等工質(zhì)的池沸騰。研究發(fā)現(xiàn),隨著壓力增加,氣泡脫離直徑減小,但變化不大,見(jiàn)圖1[15]。

圖1 不同工質(zhì)氣泡脫離直徑隨壓力的變化[15]

② 微重力

微重力條件下核態(tài)沸騰中氣泡的受力情況和常重力下有所不同,導(dǎo)致氣泡動(dòng)力學(xué)參數(shù)不同。因此,一些學(xué)者研究了重力對(duì)核態(tài)沸騰中氣泡動(dòng)力學(xué)參數(shù)的影響,并總結(jié)得出氣泡動(dòng)力學(xué)參數(shù)與重力加速度(符號(hào)為g,單位為m/s2)的關(guān)系。

Hazi等人[16]使用格子玻爾茲曼方法(Lattice Boltzmann Method,LBM)研究了水在水平受熱面上的核態(tài)沸騰中重力對(duì)氣泡脫離直徑和氣泡脫離頻率的影響。發(fā)現(xiàn)氣泡脫離直徑與g-0.5成正相關(guān),氣泡脫離頻率與g-0.75成正相關(guān)。

Gong等人[17]及曾建邦等人[18]也用LBM方法研究了水的核態(tài)沸騰過(guò)程,得到了類(lèi)似的關(guān)系式。易天浩等人[19]研究了在微重力條件下水池沸騰中單氣泡的動(dòng)力學(xué)行為及傳熱特性,模擬結(jié)果表明,氣泡脫離直徑與g-0.488成正相關(guān),氣泡生長(zhǎng)周期與g-1.113成正相關(guān)。趙建福等人[20]模擬了重力條件對(duì)水池沸騰的影響,發(fā)現(xiàn)重力對(duì)氣泡脫離直徑和氣泡生長(zhǎng)時(shí)間影響較大,氣泡脫離直徑與g-1/3成正相關(guān),氣泡生長(zhǎng)時(shí)間與g-4/5成正相關(guān)。

此外,還有學(xué)者定性研究了微重力條件對(duì)核態(tài)沸騰中氣泡脫離直徑和氣泡脫離頻率等氣泡動(dòng)力學(xué)參數(shù)的影響。Ma等人[21]采用數(shù)值方法研究了微重力對(duì)恒定壁溫水平親水受熱面氣泡動(dòng)力學(xué)參數(shù)的影響,研究發(fā)現(xiàn),重力水平的降低會(huì)導(dǎo)致氣泡脫離直徑和氣泡生長(zhǎng)周期增加。Yi等人[22]、金效興[23]、趙銳等人[24]也發(fā)現(xiàn)了氣泡脫離直徑及氣泡脫離時(shí)間隨重力變化的相同研究結(jié)果。楊燕等人[25]采用數(shù)值模擬方法研究了微重力條件對(duì)水池沸騰中單氣泡生長(zhǎng)特性的影響,模擬結(jié)果表明,隨著重力加速度減小,氣泡脫離直徑增大。

③ 受熱面熱流密度

受熱面熱流密度對(duì)核態(tài)沸騰中活性成核點(diǎn)密度、氣泡脫離直徑、氣泡脫離頻率有很大影響。

Gong等人[26]在純物質(zhì)0.1、0.3、0.5 MPa以及混合物0.3 MPa的壓力下,20～150 kW/m2的受熱面熱流密度范圍內(nèi),可視化觀測(cè)了乙烷(R170)、異丁烷(R600a)以及兩者混合物的核態(tài)沸騰。發(fā)現(xiàn),受熱面熱流密度增加時(shí),純物質(zhì)和混合物的氣泡脫離直徑、氣泡脫離頻率均顯著增加。Hamzekhani等人[27]在研究了純水、100%純度乙醇和二元混合物(乙醇/水、NaCl/水、Na2SO4/水)在大氣壓下的池沸騰過(guò)程中氣泡行為特性,結(jié)果表明,幾種工質(zhì)的氣泡脫離直徑均隨受熱面熱流密度增加而增加,且氣泡脫離直徑波動(dòng)幅度隨受熱面熱流密度增加而增大,硫酸鈉溶液質(zhì)量濃度為300 kg/m3時(shí)氣泡脫離直徑波動(dòng)最大,分別見(jiàn)圖2、3[27]。

圖2 純水、100%純度乙醇、乙醇溶液氣泡脫離直徑隨受熱面熱流密度的變化[27]

圖3 氯化鈉溶液、硫酸鈉溶液氣泡脫離直徑隨受熱面熱流密度的變化[27]

Hetsroni等人[28]在大氣壓條件下研究了受熱面熱流密度為10、50 kW/m2時(shí)對(duì)水和表面活性劑溶液池沸騰中氣泡生長(zhǎng)過(guò)程的影響,結(jié)果發(fā)現(xiàn),低熱流密度(10 kW/m2)時(shí),水和表面活性劑溶液中氣泡體積、生長(zhǎng)速率等差別不大。但在高熱流密度(50 kW/m2)時(shí),表面活性劑溶液中沸騰劇烈。隨著受熱面熱流密度增加,水中氣泡脫離直徑增大,而表面活性劑溶液中氣泡脫離直徑卻減小。

Chien等人[29]研究發(fā)現(xiàn),隨著受熱面熱流密度增大,R-123核態(tài)沸騰中氣泡等待時(shí)間和氣泡脫離直徑均減小,而活性成核點(diǎn)密度和氣泡脫離頻率均增加。鐘達(dá)文等人[30]研究了熱流密度對(duì)受熱面朝下的水池沸騰氣泡動(dòng)力學(xué)參數(shù)的影響因素,隨著受熱面熱流密度增加,氣泡生長(zhǎng)時(shí)間先減小后維持穩(wěn)定。姚遠(yuǎn)等人[31]用實(shí)驗(yàn)研究了壓力為0.2 MPa時(shí),在水平銅受熱面上熱流密度(14.65～80.79 kW/m2)對(duì)乙烷池沸騰的影響,實(shí)驗(yàn)發(fā)現(xiàn),隨著受熱面熱流密度增大,氣泡脫離直徑和脫離頻率均增大,脫離形狀分別為球形、橢球形和不規(guī)則形狀。Mchale等人[32]在大氣壓條件下研究了FC-77在光滑和粗糙兩種表面上池沸騰時(shí)氣泡成核情況,氣泡脫離頻率和上升速度均隨著熱流密度增加而增加。刁彥華等人[33]研究了R113在0.1 MPa時(shí)池沸騰中熱流密度對(duì)氣泡生長(zhǎng)、聚并、脫離過(guò)程的影響,研究發(fā)現(xiàn),隨著熱流密度增加,氣泡脫離時(shí)間縮短,活性成核點(diǎn)密度增大。

④ 受熱面過(guò)熱度和工質(zhì)過(guò)冷度

受熱面過(guò)熱度對(duì)氣泡動(dòng)力學(xué)參數(shù)有一定影響,一些學(xué)者針對(duì)這一問(wèn)題開(kāi)展了深入研究。Hutter等人[11]研究了FC-72在池沸騰中硅表面過(guò)熱度對(duì)氣泡脫離頻率、脫離直徑及等待時(shí)間等的影響,結(jié)果表明,高過(guò)熱度時(shí),氣泡脫離頻率幾乎不變,氣泡脫離直徑隨過(guò)熱度增加近似線(xiàn)性增加,氣泡等待時(shí)間隨過(guò)熱度增加而迅速縮短。Gong等人[17]也研究發(fā)現(xiàn),氣泡脫離直徑隨受熱面過(guò)熱度增加線(xiàn)性增加。趙建福等人[20]用Level Set數(shù)值模擬方法研究了常壓下水池沸騰過(guò)程中單氣泡的成長(zhǎng)和脫離過(guò)程,隨著受熱面過(guò)熱度增加,氣泡脫離直徑增加,氣泡生長(zhǎng)時(shí)間縮短。Mchale等人[32]研究了常壓下FC-77在粗糙受熱面和光滑受熱面上的池沸騰過(guò)程。研究發(fā)現(xiàn),隨著受熱面過(guò)熱度增加,氣泡脫離直徑、氣泡脫離頻率、活性成核點(diǎn)密度均增加。Zhang等人[34]研究了常壓下液氮池沸騰,發(fā)現(xiàn)氣泡脫離直徑隨受熱面過(guò)熱度增加而增大。Gao等人[35]研究了水在亞大氣壓(3.6～22.0 kPa)下不銹鋼受熱面池沸騰中氣泡動(dòng)力學(xué)參數(shù),研究表明,隨受熱面過(guò)熱度增大,熱邊界層變厚,氣泡脫離直徑變大。壓力為5.6 kPa時(shí)不同受熱面過(guò)熱度的水池沸騰中氣泡等效直徑隨時(shí)間的變化見(jiàn)圖4[35]。

此外,池沸騰中工質(zhì)溫度對(duì)氣泡動(dòng)力學(xué)參數(shù)和傳熱特性也有很大影響。因此,一些學(xué)者研究了工質(zhì)過(guò)冷度對(duì)氣泡動(dòng)力學(xué)參數(shù)的影響。畢景良等人[36]研究FC-72池沸騰發(fā)現(xiàn),氣泡脫離頻率和氣泡脫離直徑均隨工質(zhì)過(guò)冷度增加而減小,由于不同過(guò)冷度下氣泡的生長(zhǎng)速度不同,因此氣泡的生長(zhǎng)狀況不同,也導(dǎo)致氣泡生長(zhǎng)周期內(nèi)不同時(shí)刻受熱面熱流密度的變化。Demiray等人[37]采用微加熱器陣列研究了工質(zhì)過(guò)冷度對(duì)FC-72池沸騰的影響。實(shí)驗(yàn)表明,在低過(guò)冷度下,單氣泡脫離直徑和能量傳遞較大;高過(guò)冷度時(shí),氣泡脫離頻率增加導(dǎo)致傳熱效率提高。Kim等人[38]用實(shí)驗(yàn)方法研究了水池沸騰中工質(zhì)過(guò)冷度對(duì)氣泡動(dòng)力學(xué)參數(shù)和周?chē)黧w流動(dòng)的影響,實(shí)驗(yàn)結(jié)果表明,隨著工質(zhì)過(guò)冷度增加,氣泡生長(zhǎng)時(shí)間呈指數(shù)減小,氣泡脫離直徑也隨之減小。Coulibaly等人[39]研究了受熱面恒溫條件下過(guò)冷核態(tài)沸騰中氣泡的聚并過(guò)程,實(shí)驗(yàn)結(jié)果表明,氣泡脫離直徑和氣泡脫離頻率隨工質(zhì)過(guò)冷度增加而減小,且與單氣泡相比,氣泡聚并導(dǎo)致受熱面熱流密度波動(dòng)較大。

⑤ 受熱面傾角

對(duì)核態(tài)沸騰的早期研究主要在水平向上的受熱面開(kāi)展,觀察氣泡動(dòng)力學(xué)參數(shù)和對(duì)傳熱特性的影響。然而水平受熱面易較早出現(xiàn)壁面干燥,從而降低傳熱性能,導(dǎo)致過(guò)早達(dá)到臨界熱流密度,甚至將設(shè)備燒壞。加之一些設(shè)備的受熱面也確實(shí)存在一定傾角,因此部分學(xué)者研究了受熱面傾角對(duì)池沸騰的影響,池沸騰中受熱面傾角見(jiàn)圖5[40]。由圖5可知,0°為受熱面水平向上,90°為受熱面豎直向右,180°為受熱面水平向下。

圖5 池沸騰中受熱面傾角[40]

Tanjung等人[41]對(duì)去離子水池沸騰中印刷電路板(Printed Circuit Board,PCB)不同傾角對(duì)氣泡行為的影響開(kāi)展了實(shí)驗(yàn)研究,受熱面傾角為0°和45°時(shí),氣泡在壁面的垂直方向生成和脫離,但45°時(shí)氣泡會(huì)在脫離受熱面前滑移一段距離。傾角為90°和135°時(shí),氣泡生成后沿受熱面向上滑移。傾角為180°時(shí),氣泡產(chǎn)生后在壁面聚集。當(dāng)傾角為90°時(shí),核態(tài)沸騰開(kāi)始的受熱面熱流密度、氣泡等待時(shí)間、氣泡脫離直徑、活性成核點(diǎn)密度最大,氣泡脫離頻率最小。

Chuang等人[42]研究了不同受熱面傾角對(duì)水池沸騰中氣泡動(dòng)力學(xué)參數(shù)的影響,根據(jù)實(shí)驗(yàn)結(jié)果得知,隨著受熱面傾角增加,氣泡脫離直徑和氣泡脫離頻率增加,觀察發(fā)現(xiàn)活性成核點(diǎn)密度與傾角無(wú)關(guān)。Wang等人[43]實(shí)驗(yàn)研究了去離子水在常壓下池沸騰中氣泡動(dòng)力學(xué)參數(shù)受加熱面傾角的影響,研究發(fā)現(xiàn),傾角大于90°時(shí),隨受熱面傾角增大,氣泡脫離直徑增大,氣泡脫離頻率減小。

⑥ 受熱面表面粗糙度

核態(tài)沸騰中受熱面表面粗糙度對(duì)氣泡動(dòng)力學(xué)參數(shù)和傳熱過(guò)程有重要影響,因?yàn)椴煌砻娲植诙纫讓?dǎo)致活性成核點(diǎn)密度、氣泡脫離頻率及氣泡脫離直徑不同,從而導(dǎo)致傳熱性能不同。

Suszko等人[44]研究了PF-5600電介質(zhì)溶液在銅受熱面上核態(tài)沸騰過(guò)程,實(shí)驗(yàn)結(jié)果表明,隨著受熱面表面粗糙度增加,活性成核點(diǎn)密度和氣泡脫離頻率增大,但氣泡脫離直徑減小。Bovard等人[45]實(shí)驗(yàn)研究了純工質(zhì)(純水、乙醇、丙酮)的池沸騰氣泡動(dòng)力學(xué)參數(shù),結(jié)果發(fā)現(xiàn),隨表面粗糙度增加,氣泡脫離直徑減小。Mchale等人[32]研究了FC-77在不同表面粗糙度的受熱面池沸騰中的氣泡動(dòng)力學(xué)參數(shù),結(jié)果表明,表面粗糙度對(duì)氣泡脫離直徑有影響,活性成核點(diǎn)密度隨受熱面表面粗糙度增大而增大。Wang等人[43]實(shí)驗(yàn)研究了去離子水在常壓下氣泡動(dòng)力學(xué)參數(shù)受受熱面表面粗糙度的影響,研究發(fā)現(xiàn),隨著受熱面表面粗糙度增大,沸騰傳熱系數(shù)增大,受熱面傾角90°～180°時(shí),光滑受熱面上的氣泡脫離直徑小,氣泡脫離頻率大。

⑦ 接觸角

受熱面潤(rùn)濕性是影響核態(tài)沸騰的重要因素,它可以通過(guò)接觸角來(lái)表示。接觸角表征液體和固體壁面接觸的能力,是固、液和氣三相交界處形成的角,池沸騰中一般分為動(dòng)態(tài)接觸角、靜態(tài)接觸角。實(shí)際的氣泡生長(zhǎng)過(guò)程中隨著氣泡基底的收縮和擴(kuò)張,接觸角也在動(dòng)態(tài)變化,即動(dòng)態(tài)接觸角有利于準(zhǔn)確分析氣泡生長(zhǎng)過(guò)程。靜態(tài)接觸角為動(dòng)態(tài)接觸角范圍內(nèi)的某一確定值,一些研究中為了簡(jiǎn)化模型,將氣泡生長(zhǎng)過(guò)程中的接觸角設(shè)為恒定值,即靜態(tài)接觸角。未特殊指明時(shí),本文的接觸角指動(dòng)態(tài)接觸角。

Wang等人[46]研究了接觸角對(duì)水池沸騰中氣泡動(dòng)力學(xué)參數(shù)和傳熱特性的影響,結(jié)果表明,氣泡脫離時(shí)間隨接觸角增加而增加。Hazi等人[16]根據(jù)數(shù)值模擬結(jié)果得到水池沸騰中,靜態(tài)接觸角與氣泡脫離直徑相關(guān)性不大,但氣泡脫離頻率隨靜態(tài)接觸角增大而增大。Gong等人[17]采用LBM方法研究了水池沸騰中氣泡的生成、生長(zhǎng)、脫離過(guò)程,結(jié)果發(fā)現(xiàn),氣泡生長(zhǎng)周期隨著接觸角增加而增加,且接觸角對(duì)氣泡脫離直徑的影響不大。Ardron等人[47]對(duì)池沸騰中的接觸角進(jìn)行了建模,結(jié)果表明,接觸角對(duì)高壓下氣泡脫離直徑有很大影響。

⑧ 汽化核心間距

汽化核心間距是影響氣泡動(dòng)力學(xué)參數(shù)的一個(gè)重要因素,學(xué)者們對(duì)汽化核心間距進(jìn)行了一定研究。其中一些研究在受熱面上設(shè)置人工空腔,當(dāng)作汽化核心,空腔間距即為汽化核心間距,空腔排列方式即汽化核心排列方式。Zhang等人[48]研究了有2個(gè)人工空腔的受熱面上蒸餾水池沸騰中人工空腔的相互作用,發(fā)現(xiàn)影響人工空腔的3個(gè)因素:氣泡間流體動(dòng)力相互作用、人工空腔間熱相互作用、水平和傾斜氣泡聚并。當(dāng)汽化核心間距S與氣泡脫離直徑Dd的比S/Dd>3時(shí),人工空腔相互獨(dú)立,不受3個(gè)因素的影響。2

Nitesh等人[49]研究了大氣壓下FC-72池沸騰中空腔間距對(duì)氣泡動(dòng)力學(xué)參數(shù)的影響,研究發(fā)現(xiàn),空腔間距對(duì)氣泡脫離頻率和脫離直徑無(wú)影響,但空腔間距對(duì)活性成核點(diǎn)密度有很大影響,空腔間距為0.75 mm時(shí),活性成核點(diǎn)密度最大,空腔間距為0.5 mm時(shí)活性成核點(diǎn)密度最小。Hutter等人[50]研究了FC-72的人工空腔成核位點(diǎn)相互作用的池沸騰過(guò)程,結(jié)果表明,氣泡脫離直徑隨汽化核心間距增加而增加。Golobic等人[51]用實(shí)驗(yàn)研究了水池沸騰中受熱面上2～4個(gè)人工空腔間不同排列方式的相互作用,空腔間距為2.6～4.1 mm,結(jié)果表明,人工空腔為2個(gè)時(shí),人工空腔之間的相互作用,使1個(gè)人工空腔活性增加,另1個(gè)人工空腔活性降低甚至直接失去活性。Bon等人[52]研究了FC-72和正己烷核態(tài)沸騰過(guò)程,空腔間距為75～600 μm,結(jié)果表明,與光滑受熱面相比,有圓柱形空腔的受熱面?zhèn)鳠嵯禂?shù)提高,隨著空腔深度增加,活性成核點(diǎn)密度增加。Bi等人[53]用微陣列加熱器研究了FC-72池沸騰中汽化核心間距對(duì)雙氣泡聚并的影響,S/Dd為0.49時(shí),由于距離太小,生成的氣泡立即聚并為1個(gè)大氣泡,與單氣泡相比氣泡脫離頻率減小。S/Dd為0.80時(shí),由于氣泡發(fā)生連續(xù)聚并現(xiàn)象,雙氣泡先發(fā)生水平聚并,然后與新生成的氣泡再發(fā)生聚并,最后離開(kāi)受熱面。S/Dd為0.98時(shí),只發(fā)生雙氣泡間的水平聚并,在短時(shí)間內(nèi)完成聚并與脫離過(guò)程,受熱面熱流密度增大不是很明顯。S/Dd為0.94時(shí),氣泡在上升過(guò)程中聚并。Chatpun等人[54]在低熱流密度范圍內(nèi),采用實(shí)驗(yàn)方法研究了不同空腔距離和不同空腔排列方式對(duì)蒸餾水池沸騰氣泡動(dòng)力學(xué)參數(shù)的影響,結(jié)果表明,氣泡脫離直徑、氣泡脫離頻率均受空腔距離和空腔排列方式的影響。

綜合以上文獻(xiàn)可知,雙氣泡及多氣泡的生長(zhǎng)、聚并、脫離及氣泡動(dòng)力學(xué)參數(shù)、傳熱特性均受到汽化核心間距和排列方式的影響,不同汽化核心間距易導(dǎo)致氣泡聚并的不同結(jié)果,同時(shí)氣泡的聚并頻率、聚并后脫離直徑及脫離頻率和傳熱效果也不同。

⑨ 外加場(chǎng)

外加電場(chǎng)對(duì)核態(tài)沸騰的氣泡動(dòng)力學(xué)參數(shù)有很大影響,當(dāng)池沸騰中外加電場(chǎng)作用時(shí)會(huì)導(dǎo)致電流體動(dòng)力學(xué)效應(yīng)(EHD)。Gao等人[55]研究了外加電場(chǎng)對(duì)R113池沸騰中氣泡形成、生長(zhǎng)、脫離過(guò)程的影響,結(jié)果表明,隨著外加電場(chǎng)電壓增大,氣泡脫離直徑及氣泡脫離頻率均減小,見(jiàn)圖6、7[55]。

圖6 不同受熱面熱流密度下外加電場(chǎng)電壓對(duì)氣泡脫離直徑的影響[55]

圖7 不同受熱面熱流密度下外加電場(chǎng)電壓對(duì)氣泡脫離頻率的影響[55]

陳帥[56]研究了不同重力條件下均勻電場(chǎng)和非均勻電場(chǎng)對(duì)氣泡生長(zhǎng)、脫離過(guò)程的影響。常重力下,均勻電場(chǎng)作用使氣泡被橫向壓縮,縮短了收縮與脫離時(shí)間。隨重力降低,電場(chǎng)對(duì)加速氣泡脫離的作用減弱。相同重力下,均勻電場(chǎng)作用增大了氣泡上升速度。與均勻電場(chǎng)相比,非均勻電場(chǎng)在常重力下加速了氣泡脫離頻率,非均勻電場(chǎng)中,電場(chǎng)強(qiáng)度越大,氣泡上升速度越快。

外加磁場(chǎng)對(duì)池沸騰中氣泡動(dòng)力學(xué)參數(shù)及傳熱速率有一定影響,但在這方面的研究較少。目前的一些研究主要為在工質(zhì)中添加磁性物質(zhì),然后研究磁場(chǎng)的影響。Rahmati等人[57]研究了水在可變磁場(chǎng)中的池沸騰過(guò)程。實(shí)驗(yàn)通過(guò)在沸騰池中放置鐵磁珠,并改變電壓來(lái)控制磁場(chǎng),結(jié)果發(fā)現(xiàn),無(wú)磁場(chǎng)時(shí),沸騰傳熱系數(shù)變化不明顯,而有磁場(chǎng)時(shí),沸騰傳熱系數(shù)增加,氣泡脫離速率加快。Liu等人[58]實(shí)驗(yàn)研究了有無(wú)磁場(chǎng)時(shí),水基磁性流體的池沸騰過(guò)程。結(jié)果表明,施加非均勻磁場(chǎng)時(shí),氣泡脫離直徑減小。在磁場(chǎng)作用下,氣泡的形狀發(fā)生很大改變,氣泡底部面積變大,在中間形成氣泡頸,較大的底面積可以吸收更多的熱量,加速了氣泡生長(zhǎng),導(dǎo)致氣泡更快脫離,最終強(qiáng)化沸騰傳熱。

超聲波是池沸騰中氣泡動(dòng)力學(xué)參數(shù)和傳熱特性的影響因素之一,主要通過(guò)對(duì)池沸騰中的液體施加超聲波,研究影響作用。Khooshechin等人[59]通過(guò)添加表面活性劑和施加超聲波兩種方法對(duì)納米流體池沸騰進(jìn)行了實(shí)驗(yàn)研究,他們發(fā)現(xiàn),超聲波導(dǎo)致不同受熱面熱流密度下的沸騰傳熱系數(shù)均有所增加,同時(shí)超聲波導(dǎo)致氣泡脫離直徑增大。

3 氣泡動(dòng)力學(xué)參數(shù)經(jīng)驗(yàn)關(guān)聯(lián)式

3.1 氣泡脫離直徑

氣泡脫離直徑是重要的氣泡動(dòng)力學(xué)參數(shù),通?？梢杂脤?shí)驗(yàn)或力平衡的方法確定。采用實(shí)驗(yàn)方法時(shí),使用高速攝像機(jī)記錄沸騰過(guò)程,然后用圖像處理軟件處理拍攝結(jié)果,最終得到氣泡脫離直徑。另一種方法是通過(guò)分析計(jì)算氣泡脫離過(guò)程中所受各項(xiàng)力的平衡,獲得氣泡脫離直徑[60]。氣泡脫離過(guò)程中受力可以分為促進(jìn)氣泡脫離的力和阻礙氣泡脫離的力[61]。浮力和壓力為促進(jìn)氣泡脫離的力,液體慣性力、阻力和表面張力則是阻礙氣泡脫離的力?；谏鲜龇治?學(xué)者們開(kāi)發(fā)了關(guān)于氣泡脫離直徑的關(guān)聯(lián)式。

Fritz[62]考慮了使氣泡從受熱面脫離的浮力和使氣泡保持在受熱面的表面張力的平衡,得到了水和溶液的氣泡脫離直徑Dd關(guān)聯(lián)式(簡(jiǎn)稱(chēng)Fritz關(guān)聯(lián)式):

式中Dd——?dú)馀菝撾x直徑,m

θ——接觸角,(°)

σ——表面張力,N/m

g——重力加速度,m/s2

ρL——工質(zhì)液相密度,kg/m3

ρv——工質(zhì)氣相密度,kg/m3

Fritz關(guān)聯(lián)式中,水的接觸角為45°,溶液的接觸角為30°。由Fritz關(guān)聯(lián)式可知,氣泡脫離直徑和接觸角、表面張力以及工質(zhì)物理性質(zhì)有關(guān)。Fritz關(guān)聯(lián)式作為氣泡脫離直徑的基本形式,僅考慮了浮力和表面張力,接下來(lái)關(guān)于氣泡脫離直徑的部分研究是關(guān)于Fritz關(guān)聯(lián)式的修正。

Cole等人[63]實(shí)驗(yàn)觀察發(fā)現(xiàn),氣泡生長(zhǎng)速率和氣泡脫離直徑成正相關(guān),因此考慮氣泡生長(zhǎng)速率,修正了Fritz關(guān)聯(lián)式。Ruckenstein[64]、Stralen等人[65]以及Gorenflo等人[66]修正了Fritz關(guān)聯(lián)式,包含了工質(zhì)熱擴(kuò)散系數(shù)和雅各布數(shù)。Kocamustafaogullari[67]通過(guò)包括工質(zhì)液相和氣相的密度比來(lái)修正Fritz關(guān)聯(lián)式。Wenzel[68]通過(guò)包括雅各布數(shù)、普朗特?cái)?shù)和阿基米德數(shù)修正了Fritz關(guān)聯(lián)式。Phan等人[69]在Fritz關(guān)聯(lián)式基礎(chǔ)上進(jìn)行了修正,研究了接觸角對(duì)氣泡脫離的影響,將氣泡脫離直徑與接觸角和工質(zhì)物理性質(zhì)關(guān)聯(lián)起來(lái)。Nam等人[70]在具有隔離微腔的硅襯底上形成CuO納米結(jié)構(gòu)的超親水表面進(jìn)行了單氣泡動(dòng)力學(xué)特性研究,考慮了作用于氣泡的浮力和表面張力的平衡,得到了氣泡脫離直徑的模型。

Cole等人的關(guān)聯(lián)式為[63]:

式中g(shù)c——重力修正系數(shù)

D——?dú)馀葜睆?m

Ruckenstein的關(guān)聯(lián)式為[64]:

式中aL——工質(zhì)液相熱擴(kuò)散率,m2/s

Ja——雅各布數(shù)

Stralen等人的關(guān)聯(lián)式為[65]:

Gorenflo等人的關(guān)聯(lián)式為[66]:

式中c——系數(shù)[66]

Kocamustafaogullari的關(guān)聯(lián)式為[67]:

Wenzel的關(guān)聯(lián)式為[68]:

式中Pr——普朗特?cái)?shù)

Ar——阿基米德數(shù)

Phan等人的關(guān)聯(lián)式為[69]:

Nam等人的關(guān)聯(lián)式為[70]:

以上關(guān)于氣泡脫離直徑的關(guān)聯(lián)式及模型是基于Fritz關(guān)聯(lián)式的相關(guān)性展開(kāi)。此外,其他學(xué)者也通過(guò)實(shí)驗(yàn)研究或數(shù)值模擬研究方法提出了關(guān)于氣泡脫離直徑的關(guān)聯(lián)式。Jamialahmadi等人[71]進(jìn)行了圓柱形不銹鋼加熱器外部NaCl溶液、KNO3溶液和Na2SO4溶液的池沸騰實(shí)驗(yàn),測(cè)量了氣泡脫離直徑,提出了考慮受熱面熱流密度的電解質(zhì)溶液氣泡脫離直徑的關(guān)聯(lián)式。Kim等人[72]得到了亞大氣壓條件下氣泡脫離直徑的關(guān)聯(lián)式,將氣泡脫離直徑與雅各布數(shù)和工質(zhì)物理性質(zhì)聯(lián)系起來(lái),相關(guān)性預(yù)測(cè)實(shí)驗(yàn)數(shù)據(jù)的相對(duì)誤差在±20%以?xún)?nèi)。Suszko等人[44]在受熱面熱流密度0.5 W/cm2下,對(duì)PF-5060池沸騰中光滑和粗糙的銅受熱面進(jìn)行了實(shí)驗(yàn),分別得到了光滑和粗糙受熱面氣泡脫離直徑與氣泡生長(zhǎng)時(shí)間的關(guān)聯(lián)式。

Jamialahmadi等人的關(guān)聯(lián)式為[71]:

式中q——受熱面熱流密度,W/m2

Kim等人的關(guān)聯(lián)式為[72]:

Suszko的關(guān)聯(lián)式為[44]:

光滑表面:

粗糙表面:

式中tg——?dú)馀萆L(zhǎng)時(shí)間,s

以上研究表明,氣泡脫離直徑與工質(zhì)物理性質(zhì)、表面張力以及接觸角、受熱面熱流密度等參數(shù)有關(guān)。

3.2 氣泡脫離頻率

氣泡脫離頻率是研究核態(tài)沸騰過(guò)程的另一個(gè)重要的氣泡動(dòng)力學(xué)參數(shù)[60],取決于氣泡等待時(shí)間和氣泡生長(zhǎng)時(shí)間,是兩者之和的倒數(shù):

式中f——?dú)馀菝撾x頻率,Hz

tw——?dú)馀莸却龝r(shí)間,s

大量研究發(fā)現(xiàn),氣泡脫離頻率與受熱面過(guò)熱度、工質(zhì)物理性質(zhì)、接觸角、空腔距離及相鄰氣泡相互作用等有關(guān),受熱面表面粗糙度也有明顯影響[31]。由于早期的研究單獨(dú)確定氣泡生長(zhǎng)時(shí)間和等待時(shí)間比較困難,因此常把氣泡脫離頻率與其他參數(shù)結(jié)合起來(lái)計(jì)算。Jakob等人[73]研究發(fā)現(xiàn),在液氮和水的核態(tài)沸騰過(guò)程中,氣泡脫離直徑和脫離頻率的乘積是恒定的:

fDd=0.078

此外,Hatton等人[74]提出了一個(gè)相關(guān)性,該相關(guān)性顯示了氣泡脫離直徑與氣泡脫離頻率二次方的乘積和工質(zhì)熱擴(kuò)散率的關(guān)系:

Cole[75]給出了氣泡脫離頻率和氣泡脫離直徑的表達(dá)式:

根據(jù)實(shí)驗(yàn)數(shù)據(jù),Ivey[76]針對(duì)水核態(tài)沸騰提出了3個(gè)不同區(qū)域:動(dòng)力區(qū)域、過(guò)渡區(qū)域、熱力學(xué)區(qū)域。

動(dòng)力區(qū)域:浮力和阻力作用于不斷增長(zhǎng)的氣泡,有關(guān)聯(lián)式:

過(guò)渡區(qū)域:浮力、阻力和表面張力為主要力,作用于氣泡,有關(guān)聯(lián)式:

熱力學(xué)區(qū)域:氣泡增長(zhǎng)占主導(dǎo)地位,有關(guān)聯(lián)式:

式中C——常數(shù)

Kumada等人[77]研究了高熱流密度和大氣壓下圓形水平圓盤(pán)上氣泡脫離頻率,通過(guò)對(duì)力平衡方程進(jìn)行分析,提出了氣泡脫離頻率的半經(jīng)驗(yàn)方程,方程包括密度比、運(yùn)動(dòng)黏度和受熱面尺寸。Sakashita等人[78]在大氣壓至7 MPa的壓力下,研究了水池沸騰中高熱流密度下氣泡脫離頻率,提出了氣泡脫離頻率關(guān)聯(lián)式,氣泡脫離頻率不受受熱面熱流密度和壓力的影響,關(guān)聯(lián)式中包含工質(zhì)密度比、運(yùn)動(dòng)黏度和表面張力。Hamzekhani等人[79]利用實(shí)驗(yàn)方法研究了水和氯化鈉溶液在常壓、不同熱流密度下的池沸騰過(guò)程,提出了預(yù)測(cè)氣泡脫離頻率的關(guān)聯(lián)式,關(guān)聯(lián)式中包含氣泡脫離直徑、表面張力、重力加速度、氣液密度差和熱流密度等參數(shù)。

Kumada等人的關(guān)聯(lián)式為[77]:

式中νL——液相運(yùn)動(dòng)黏度,m2/s

Ds——受熱面直徑,mm

Sakashita等人的關(guān)聯(lián)式為[78]:

Hamzekhani等人的關(guān)聯(lián)式為[79]:

式中 Δρ——液相和氣相的密度差,kg/m3

綜合以上研究,大多數(shù)關(guān)聯(lián)式中氣泡脫離頻率與氣泡脫離直徑相關(guān)。氣泡脫離頻率受各種狀態(tài)和氣泡特性相關(guān)參數(shù)的影響,如氣泡脫離直徑、氣泡等待時(shí)間、氣泡生長(zhǎng)時(shí)間、表面張力、受熱面熱流密度和工質(zhì)熱物理性質(zhì),無(wú)法全面準(zhǔn)確預(yù)測(cè)所有情況的氣泡脫離頻率。因此,應(yīng)綜合考慮受熱面表面粗糙度、空腔尺寸、壓力等,從而建立適用性更強(qiáng)的氣泡脫離頻率關(guān)聯(lián)式。

3.3 氣泡等待時(shí)間

氣泡等待時(shí)間是一個(gè)重要的氣泡動(dòng)力學(xué)參數(shù),對(duì)池沸騰中受熱面的傳熱速率有重要影響,學(xué)者們提出了許多關(guān)聯(lián)式。

Han等人[80]提出了氣泡等待時(shí)間的關(guān)聯(lián)式,結(jié)果表明,隨著空腔尺寸增大,氣泡等待時(shí)間先增大后減小,隨著受熱面過(guò)熱度增加,氣泡等待時(shí)間減小:

式中Rc——空腔半徑,μm

Tw——受熱面溫度,K

TL——流體溫度,K

Tsat——工質(zhì)飽和溫度,K

r——汽化潛熱,kJ/kg

Stralen等人[81]建立了純液體池沸騰中氣泡等待時(shí)間和氣泡生長(zhǎng)時(shí)間之間的關(guān)系,結(jié)果發(fā)現(xiàn),同一成核位點(diǎn),氣泡等待時(shí)間是氣泡生長(zhǎng)時(shí)間的3倍:

tw=3tg

Jeremy等人[9]考慮了氣泡生長(zhǎng)階段由微液層蒸發(fā)引起的過(guò)熱降低的效應(yīng),提出了近似的氣泡等待時(shí)間關(guān)聯(lián)式:

式中γL——流體液相體膨脹系數(shù),K-1

3.4 氣泡生長(zhǎng)時(shí)間

氣泡生長(zhǎng)時(shí)間內(nèi)初始?xì)馀莶粩辔諢崃恐钡綒馀菝撾x,對(duì)氣泡動(dòng)力學(xué)特征有顯著的影響。Zuber[82]研究了非均勻溫度場(chǎng)中氣泡動(dòng)力學(xué)特性,提出了氣泡生長(zhǎng)時(shí)間與工質(zhì)熱擴(kuò)散率、氣泡脫離直徑和雅各布數(shù)有關(guān):

Lee等人[83]在不同邊界條件下,提出了氣泡生長(zhǎng)時(shí)間的預(yù)測(cè)關(guān)聯(lián)式,氣泡生長(zhǎng)時(shí)間與氣泡脫離直徑和雅各布數(shù)成正相關(guān):

以上對(duì)氣泡等待時(shí)間和生長(zhǎng)時(shí)間的總結(jié)表明,空腔尺寸、氣泡脫離直徑、受熱面過(guò)熱度及工質(zhì)物理性質(zhì)等對(duì)氣泡等待時(shí)間和氣泡生長(zhǎng)時(shí)間有重要影響。

3.5 活性成核點(diǎn)密度

活性成核點(diǎn)密度也是池沸騰中氣泡動(dòng)力學(xué)參數(shù)之一,隨著活性成核點(diǎn)密度增加,氣泡數(shù)量也隨之增加,最終導(dǎo)致壁面?zhèn)鳠嵯禂?shù)增加。許多研究中通常采用拋光、刻蝕、化學(xué)提拉法以及在受熱面使用各種類(lèi)型的涂層等方法來(lái)增加活性成核點(diǎn)密度[84]。研究人員開(kāi)發(fā)了許多確定成核點(diǎn)密度的關(guān)聯(lián)式。

Zou等人[85]實(shí)驗(yàn)研究了R134a過(guò)冷沸騰時(shí)不銹鋼和銅受熱面上的成核點(diǎn)分布。根據(jù)研究結(jié)果,與不銹鋼受熱面相比,銅受熱面的成核點(diǎn)分布更加均勻,這是不同熱導(dǎo)率金屬受熱面的成核點(diǎn)間熱相互作用的結(jié)果。Hibiki等人[86]根據(jù)實(shí)際受熱面,對(duì)活性成核點(diǎn)密度建立模型,模擬結(jié)果表明,活性成核點(diǎn)密度是空腔臨界尺寸和接觸角的函數(shù),隨著受熱面過(guò)熱度增加,活性成核點(diǎn)密度增加。Paul等人[87]在大氣壓條件下,研究了鉑絲上水池沸騰實(shí)驗(yàn)的氣泡動(dòng)力學(xué)參數(shù),結(jié)果表明,隨著受熱面熱流密度增加,活性成核點(diǎn)密度增加。Qi等人[88]研究了蒸餾水和乙醇在黃銅和不銹鋼受熱面的池沸騰,蒸餾水成核點(diǎn)密度與空腔臨界半徑的倒數(shù)趨勢(shì)相似,乙醇在粗糙和鏡面拋光黃銅受熱面上的活性成核點(diǎn)密度沒(méi)有顯著差異。Li等人[89]開(kāi)發(fā)了適用于壓力范圍0.101～19.8 MPa,包含受熱面過(guò)熱度、壓力、接觸角3個(gè)變量的池沸騰中活性成核點(diǎn)密度模型,研究發(fā)現(xiàn),受熱面過(guò)熱度對(duì)活性成核點(diǎn)密度的增加有促進(jìn)作用。

Benjamin等人[90]實(shí)驗(yàn)研究了中低熱流密度下蒸餾水、四氯化碳、正己烷、丙酮在不同表面粗糙度的不銹鋼和鋁壁受熱面上池沸騰中的活性成核點(diǎn)密度。結(jié)果發(fā)現(xiàn),活性成核點(diǎn)密度與受熱面表面粗糙度、表面張力、受熱面過(guò)熱度、受熱面熱物性、工質(zhì)物理性質(zhì)有關(guān),提出了活性成核點(diǎn)密度關(guān)聯(lián)式:

式中ns——活性成核點(diǎn)密度,m-2

λw——受熱面熱導(dǎo)率,W/(m·K)

ρw——受熱面密度,kg/m3

cw——受熱面比熱容,J/(kg·K)

λL——工質(zhì)液相熱導(dǎo)率,W/(m·K)

cp,L——工質(zhì)液相比定壓熱容,J/(kg·K)

R——受熱面相對(duì)表面粗糙度

Ra——受熱面平均表面粗糙度,μm

p——壓力,Pa

Sarafraz等人[91]提出了預(yù)測(cè)活性成核點(diǎn)密度的關(guān)聯(lián)式,發(fā)現(xiàn)活性成核點(diǎn)密度和受熱面熱流密度成指數(shù)函數(shù)關(guān)系:

ns=0.888 6q0.950 31

上述研究表明,活性成核點(diǎn)密度的確定涉及工質(zhì)物理性質(zhì)、表面張力、接觸角、受熱面熱流密度、受熱面過(guò)熱度、空腔尺寸等參數(shù)。此外,這些關(guān)聯(lián)式仍不能準(zhǔn)確預(yù)測(cè)實(shí)驗(yàn)數(shù)據(jù)。

4 結(jié)束語(yǔ)

綜述了不同影響因素對(duì)核態(tài)沸騰中氣泡動(dòng)力學(xué)參數(shù)的影響,整理了氣泡動(dòng)力學(xué)參數(shù)的經(jīng)驗(yàn)關(guān)聯(lián)式。雖然關(guān)于核態(tài)沸騰中氣泡動(dòng)力學(xué)特性的研究已有很多,但所提出的氣泡動(dòng)力學(xué)參數(shù)關(guān)聯(lián)式僅能準(zhǔn)確預(yù)測(cè)各自的實(shí)驗(yàn)結(jié)果,而不能準(zhǔn)確預(yù)測(cè)他人的實(shí)驗(yàn)結(jié)果。這是由于現(xiàn)有關(guān)聯(lián)式中考慮的影響因素有限,因此需要綜合考慮更多的影響因素來(lái)準(zhǔn)確預(yù)測(cè)氣泡動(dòng)力學(xué)參數(shù)。