王浩 徐進(jìn)良2)?

1) (華北電力大學(xué),低品位能源多相流與傳熱北京市重點(diǎn)實(shí)驗(yàn)室,北京 102206)

2) (華北電力大學(xué),電站能量傳遞轉(zhuǎn)化與系統(tǒng)教育部重點(diǎn)實(shí)驗(yàn)室,北京 102206)

熱油面液滴蒸發(fā)是自然現(xiàn)象,已有研究側(cè)重于單滴蒸發(fā),對(duì)于熱油面上多滴蒸發(fā)的認(rèn)識(shí)較少.本文研究了熱硅油面兩個(gè)等直徑FC-72 液滴的Leidenfrost 蒸發(fā),油溫為74.0—130.0 ℃,液滴初始直徑為1.5 mm,采用紅外熱成像及高速攝影測(cè)量,發(fā)現(xiàn)熱油面液滴蒸發(fā)存在非聚合、彈跳、分離3 個(gè)階段.本文理論分析了液滴在水平方向的受力,包括非均勻液滴溫度產(chǎn)生的Marangoni 力、重力水平分量、潤(rùn)滑推動(dòng)力、黏性力.尺度分析表明Marangoni力和重力水平分量起關(guān)鍵作用,Marangoni 力趨向于液滴分離,重力水平分量趨向于液滴聚合.在非聚合蒸發(fā)階段,重力水平分量克服Marangoni 力,但兩液滴間存在氣膜夾層,解釋了兩個(gè)液滴看似接觸但不聚合的現(xiàn)象.隨液滴尺寸減小,重力水平分量減小,不足以克服Marangoni 力,這是導(dǎo)致蒸發(fā)后期兩滴分離的主要原因.最后通過(guò)將模型得到的不同階段間的轉(zhuǎn)換時(shí)間同測(cè)量值進(jìn)行對(duì)比,證實(shí)了上述解釋.本文研究有助于理解復(fù)雜的Leidenfrost 液滴動(dòng)力學(xué)現(xiàn)象和機(jī)理.

1 引言

當(dāng)液滴滴落在熾熱的壁面上,其底部快速產(chǎn)生的蒸氣層使液滴懸浮,這種現(xiàn)象被稱為L(zhǎng)eidenfrost效應(yīng)[1],其發(fā)生的臨界溫度稱為L(zhǎng)eidenfrost 溫度(TL)[2].蒸氣膜的低導(dǎo)熱減慢了液滴蒸發(fā),延長(zhǎng)了液滴壽命.懸浮的液滴避免與基底接觸,減小了移動(dòng)的阻力.鑒于此,Leidenfrost 液滴在非接觸式微流控[3,4],化學(xué)和生物微反應(yīng)器[5,6],藥物無(wú)污染輸運(yùn)[7]等領(lǐng)域有著更為廣闊的應(yīng)用.

懸浮于平板表面上的液滴,蒸發(fā)過(guò)程會(huì)出現(xiàn)移動(dòng)[8]、彈跳[9]和振蕩[10]的隨機(jī)現(xiàn)象,這給液滴控制造成難度.在平板表面上制造規(guī)則鋸齒結(jié)構(gòu),可以有效地引導(dǎo)液滴底部的蒸氣流.使其在液滴兩側(cè)非對(duì)稱排放,從而打破液滴底部黏性力的平衡,驅(qū)動(dòng)液滴定向運(yùn)動(dòng)[11].此外,不均勻的壁溫分布影響了液滴局部氣化速度以及液滴內(nèi)部Marangoni 流,在液滴底部形成非對(duì)稱的蒸氣膜,引導(dǎo)液滴向冷的方向移動(dòng)[12].相比于剛性壁面,液滴在可變形表面的運(yùn)動(dòng)更為復(fù)雜.這是由于液滴和液池之間蒸氣層形態(tài)不同于平坦剛板上,氣膜厚度呈現(xiàn)不同的標(biāo)度律[13].高溫液池上懸浮的液滴,頸部蒸氣膜出口位置局部振蕩極易引起氣膜整體不穩(wěn)定,從而破壞蒸氣膜的對(duì)稱性,使液滴在液池表面自發(fā)移動(dòng)[13,14].在移動(dòng)過(guò)程中,液滴前緣出口氣膜厚度總是大于尾緣,液滴遵循蒸氣主流的方向[14,15].鑒于此,Gauthier等 [16]分析液滴與彎月面碰撞動(dòng)力學(xué),以探測(cè)界面形狀并鑒別液滴尺寸.由于液池表面光滑,液滴可在較低過(guò)熱度下進(jìn)行Leidenfrost 蒸發(fā).因此,與液滴在固體壁面上的蒸發(fā)相比,液滴在液池表面更易形成Leidenfrost 狀態(tài)[1,17].

2021 年,Pacheco 等 [18]通過(guò)觀察兩個(gè)液滴在Leidenfrost 狀態(tài)下的碰撞-合并行為,發(fā)現(xiàn)沉積在熱的凹板上不同液體的液滴碰撞時(shí),接觸區(qū)域內(nèi)較高飽和溫度的液滴加熱較低溫度液滴而發(fā)生的Leidenfrost 效應(yīng)推遲了液滴的合并.綜合液滴固-液界面發(fā)生的兩處Leidenfrost 蒸發(fā),他們將這一行為稱為三重Leidenfrost 效應(yīng) (triple Leidenfrost effect).相比于處在Leidenfrost 狀態(tài)不同液體的液滴間碰撞-合并行為,了解相同液體液滴的蒸發(fā)動(dòng)力學(xué)同樣重要,其有助于解析噴霧冷卻的相變機(jī)制,然而卻鮮有報(bào)道.對(duì)于單個(gè)液滴,蒸發(fā)速率受到蒸氣擴(kuò)散速率和液滴大小等因素的影響;然而,對(duì)于簡(jiǎn)單排列的幾個(gè)液滴或者有復(fù)雜尺寸分布的液滴,蒸發(fā)過(guò)程液滴間協(xié)同作用顯著[19,20].這種相互影響導(dǎo)致液滴內(nèi)出現(xiàn)非對(duì)稱流場(chǎng),從而影響液滴蒸發(fā)[21].因此,對(duì)于多個(gè)Leidenfrost 液滴,除了非對(duì)稱氣膜產(chǎn)生的驅(qū)動(dòng)機(jī)制外,還需要考慮液滴間相互作用對(duì)運(yùn)動(dòng)的影響,其中涉及的動(dòng)力學(xué)較為復(fù)雜.

與以往對(duì)常溫或振動(dòng)液池表面多滴作用機(jī)理研究[22-25]和高溫固體表面不同液體液滴的Leidenrost 蒸發(fā)機(jī)制研究 [18]不同,本文報(bào)道了相同液體且等大的相鄰液滴在熱油池上的Leidenfrost 蒸發(fā)行為,旨在探究Leidenfrost 態(tài)雙滴之間干涉效應(yīng)對(duì)其液滴運(yùn)動(dòng)的影響.研究發(fā)現(xiàn),雙滴與單滴相比,在液滴壽命、表面溫度和運(yùn)動(dòng)方式方面存在顯著差異;液滴表面的非對(duì)稱蒸發(fā)打破了液滴原有受力平衡,使液滴蒸發(fā)過(guò)程分為非聚合、彈跳以及最終液滴徹底分離3 個(gè)階段.從力學(xué)進(jìn)行分析,獲得3 個(gè)蒸發(fā)階段的轉(zhuǎn)換時(shí)間,它們偏離實(shí)驗(yàn)最大值僅占液滴壽命的3%,因此認(rèn)為可以較好地闡釋Leidenfrost雙液滴的蒸發(fā)機(jī)制.

2 實(shí)驗(yàn)方法

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

實(shí)驗(yàn)在溫濕度可控且安靜的房間中進(jìn)行.如圖1 所示,實(shí)驗(yàn)系統(tǒng)包含紫銅加熱塊、注射器、高速相機(jī)、紅外高速相機(jī)、同步器和控溫裝置.黏度μo=0.048 Pa·s 的硅油 (XIAMETER PMX-200,熱物性見(jiàn)表1) 盛放于銅塊頂部的矩形油槽中,槽體長(zhǎng)寬深分別為10,6 和2 mm.兩根功率為30 W的加熱棒嵌入銅塊底部.PID 控制器同電壓調(diào)節(jié)器共同確保銅塊的溫度波動(dòng)小于0.5 ℃,實(shí)驗(yàn)過(guò)程中油溫由此總能維持于給定溫度.與深液池相比,深度較淺的油槽保證了更大的銅油接觸面積與槽體積的比值,利于維持槽內(nèi)油溫的相對(duì)均勻.

圖1 (a) 實(shí)驗(yàn)裝置圖;(b),(c)放大的帶有薄液池的加熱銅塊(1-高速攝像機(jī),2-紅外高速 相機(jī),3-微量注射器,4-位移調(diào)節(jié)平臺(tái),5-電源變壓器,6-PID 溫度控制器,7-帶薄液池的銅塊,8-光源,9-用于釋放液滴的冷卻針頭)Fig.1.Photograph of experimental setup (a) and enlarged copper block with thin liquid pool (b),(c) (1-high speed camera,2-infrared radiation image camera,3-micro-syringe pump,4-displacement adjustment platform,5-voltage transformer,6-PID temperature controller,7-copper block with thin liquid pool,8-light source,9-cooled dual-needles for droplet release).

兩個(gè)FC-72 液滴 (3 MTMFluorinertTMElectronic Liquid,物性見(jiàn)表1)由微型注射泵(RSP01-BDG)產(chǎn)生,并靠自身重力脫離針尖.采用Y 形針頭,保證兩液滴初始大小相同,液滴從針頭分離同步,針頭間距5 mm.FC-72 液滴初始半徑Rini=0.75 mm,小于其毛細(xì)長(zhǎng)度l=[σd/(ρdg)]0.5=0.92 mm (σd為液滴表面張力,ρd為液滴密度,g為重力加速度).因此,液滴近似球形,避免因液滴的復(fù)雜形狀建模而帶來(lái)的計(jì)算難題.液滴從2.5 mm 的高度(≈3.3Rini)釋放,其韋伯?dāng)?shù)We=ρdv2Dini/σd=7.2,避免了因液滴沖撞池面而影響Leidenfrost轉(zhuǎn)換溫度.其中,v為液滴沖擊液面速度(與液滴釋放高度有 關(guān)),Dini為液滴初始直徑.FC-72 密度為1680 kg/m3,是硅油密度960 kg/m3的1.75 倍.在常溫條件下,當(dāng)兩種液體接觸,FC-72 液滴迅速沉沒(méi)于油池.由于兩者互不相溶,硅油的物性參數(shù)不受沉沒(méi)其中的液滴影響.池面加熱液滴導(dǎo)致的蒸發(fā)減小了液滴的尺寸,加劇了液滴脫離針頭難度.因此,針頭上裝有微型冷卻裝置,抑制液滴懸掛針頭期間的溫升.油池表面的過(guò)熱度ΔTo表示油面溫度To與液滴飽和溫度Tsat之差,本研究ΔTo介于17—70 ℃之間.

表1 在1 atm (1 atm=1.013×105 Pa) 的壓力下,FC-72 和硅油的物性參數(shù)Table 1.Physical properties of FC-72 and silicon oil at 1 atm (1 atm=1.013×105 Pa).

帶有微焦鏡頭的高速相機(jī)從水平方向捕捉液滴運(yùn)動(dòng)軌跡,空間分辨率為4 μm,拍攝頻率為2000 Hz.實(shí)驗(yàn)時(shí),油的深度略低于槽體深度,輕微下凹的油面使液滴處于油池中心位置,約束液滴的隨機(jī)運(yùn)動(dòng).油面寬度方向產(chǎn)生的彎液面確保蒸發(fā)的液滴總是位于相機(jī)焦平面上,這對(duì)清晰記錄液滴合并-碰撞行為至關(guān)重要.紅外相機(jī)從上方記錄液滴表面溫度的演化過(guò)程,在3—5 μm 的光譜范圍內(nèi)溫度分別率為0.02 ℃,空間分辨率為12.4 μm,記錄頻率為180 Hz.

2.2 溫度校驗(yàn)及液滴溫度邊界的確定

在正式實(shí)驗(yàn)之前,通過(guò)調(diào)節(jié)高精密x-y位移平臺(tái),完成熱電偶測(cè)量池內(nèi)溫度.以池面中心為基準(zhǔn)點(diǎn),熱電偶沿水平[坐標(biāo)為(x,0,0.2)]和豎直[坐標(biāo)為(0,0,z)]兩個(gè)方向的位移步長(zhǎng)分別為1.0 mm和0.2 mm.測(cè)量結(jié)果見(jiàn)圖2(a),水平方向x=0—4 mm 范圍內(nèi)最大溫差為0.5 ℃,豎直方向z=0—0.8 mm 深度內(nèi)最大溫差為0.7 ℃.因此,認(rèn)為實(shí)驗(yàn)過(guò)程中油面溫度恒定.用插入液面0.5 mm 熱電偶測(cè)得溫度與紅外測(cè)溫?cái)M合,完成校正紅外相機(jī).為了確保熱電偶測(cè)量溫度即為液面溫度,將裝有待測(cè)液體的試管放于恒溫油浴內(nèi),管內(nèi)液面比試管外油面低2—3 cm.如圖2(b)所示,兩種方法測(cè)得的溫度線性相關(guān),FC-72 和硅油的相關(guān)系數(shù)分別為0.999 和0.964.實(shí)驗(yàn)所用熱電偶精度為0.2 ℃.

圖2(c)說(shuō)明了液滴表面溫度的識(shí)別方法,紅外測(cè)溫覆蓋液滴區(qū)域及液滴外背景區(qū)域.僅當(dāng)液滴豎直方向上的厚度大于臨界光程長(zhǎng)度時(shí),所得液滴溫度是真實(shí)的.此時(shí),來(lái)自液滴背景的溫度信息不會(huì)干擾液滴溫度測(cè)量,因此液滴溫度近乎均勻.沿法線方向繼續(xù)延伸,背景信息的干擾使溫度突然上升.通常將溫度出現(xiàn)轉(zhuǎn)折的位置Tr定義為溫度邊界,它是液滴主導(dǎo)溫度區(qū)與和背景主導(dǎo)溫度區(qū)的分界線[26,27].由此可以獲得液滴表面溫度的分布.

圖2 (a) 油池在水平方向和深層方向的溫度分布;(b) 硅油和FC-72 的表面溫度測(cè)量的校準(zhǔn);(c) 通過(guò)紅外測(cè)量定位液滴界面的原理Fig.2.(a) Temperature dispersion in the oil bath’s horizontal and deep directions;(b) calibration of surface temperature measurement for silicon oil and FC-72;(c) the principle to locate the drop interface by IR measurement.

3 實(shí)驗(yàn)結(jié)果與討論

3.1 液滴的大小變化

本工作借助高速相機(jī)對(duì)高溫油面的單滴和雙滴蒸發(fā)行為進(jìn)行觀測(cè).硅油的黏度μo=0.048 Pa·s,油面溫度To包含50—130 ℃的區(qū)間,并且每個(gè)油溫下的液滴實(shí)驗(yàn)至少重復(fù)5 次.如圖3(a) 所示,當(dāng)液滴緩慢地釋放至不同溫度的液池表面時(shí),存在明顯的臨界油溫,即Leidenfrost 溫度TL.當(dāng)油溫To＜TL,沉積于油面的液滴在短暫懸停之后快速沉沒(méi)至油中 (sinking).這是因?yàn)榈蜏赜兔嬉旱握舭l(fā)緩慢,液滴底部的蒸氣膜變薄.當(dāng)液滴或油面擾動(dòng)引起的波高與蒸氣層厚度相近時(shí),液滴與油面局部區(qū)域接觸就會(huì)直接誘導(dǎo)蒸氣層坍塌.Zawala 等[28]將此過(guò)程稱為“排膜”.當(dāng)To＞TL,液滴快速蒸發(fā)增大了蒸氣膜厚度,使液滴可以長(zhǎng)期懸浮于液池之上(Leidenfrost evaporation).在本研究中,這兩種狀態(tài)間轉(zhuǎn)換的臨界溫度TL約為74.0 ℃,相應(yīng)油池表面過(guò)熱度ΔTo=17.4 ℃.

圖3(b)顯示了Leidenfrost 狀態(tài)下液滴大小隨時(shí)間的變化.對(duì)不同條件下液滴的直徑分別擬合,結(jié)果顯示液滴蒸發(fā)服從于液滴大小與時(shí)間的冪次律關(guān)系式[1]:

圖3 (a) 不同油面溫度To 下液滴壽命的兩種分區(qū);(b) 油面溫度為88.2 ℃和128.0 ℃時(shí)液滴的直徑與時(shí)間的關(guān)系Fig.3.(a) Droplet diameters versus time at oil surface temperature of 88.2 ℃ and 128.0 ℃;(b) two regimes distribution of droplet life time at different oil surface temperatures.

式中,Dini是液滴初始直徑,n是指數(shù)系數(shù),τ是特征時(shí)間,其通過(guò)外推液滴直徑D至 0 mm 處得到.所有的擬合曲線的相關(guān)性系數(shù)均大于0.994.在溫度較低To=88.2 ℃和溫度較高To=128 ℃的兩類油面上,雙滴相比單個(gè)液滴蒸發(fā)的速率減小,壽命普遍延長(zhǎng).這與Khilifi 等 [29]在常溫情況下觀察多滴蒸發(fā)得到的結(jié)論類似,他們發(fā)現(xiàn)當(dāng)液滴距離足夠近,與單個(gè)液滴相比,蒸發(fā)速率降低約55%.冪指數(shù)n是擬合的關(guān)鍵參數(shù),該值隨著油面溫度To的升高而減小,這在最近的研究中已被報(bào)道[30].本文關(guān)注于單個(gè)液滴和兩個(gè)液滴蒸發(fā)的差異.如圖3(b)所示,當(dāng)油溫To=88.2 ℃,單個(gè)液滴和兩個(gè)液滴蒸發(fā)分別對(duì)應(yīng)n=0.81,n=0.84.當(dāng)To=128.0 ℃時(shí),單個(gè)液滴和兩個(gè)液滴蒸發(fā)分別對(duì)應(yīng)n=0.64,n=0.72.結(jié)果顯示,在相同油溫To情況下,兩個(gè)液滴蒸發(fā)過(guò)程系數(shù)n總是大于單個(gè)液滴蒸發(fā)情況.圖3(b)插圖比較了不同初始直徑的單個(gè)和兩個(gè)液滴蒸發(fā)過(guò)程n值大小.顯然,圖中單個(gè)液滴數(shù)據(jù)點(diǎn)相比兩個(gè)液滴整體向上平移,驗(yàn)證了上述結(jié)論的普適性.這種差異取決于相鄰液滴之間存在的干涉蒸發(fā)效應(yīng),它使得相鄰雙滴內(nèi)側(cè)比外側(cè)有更高的蒸氣濃度,進(jìn)而抑制液滴總蒸發(fā)速率,將在后文討論.

3.2 相鄰Leidenfrost 液滴蒸發(fā)的分區(qū)

圖4 展示了相鄰的兩個(gè)Leidenfrost 液滴蒸發(fā)的動(dòng)態(tài)行為.通常情況下,當(dāng)兩個(gè)常溫的液滴接觸時(shí)(非劇烈碰撞),液滴在經(jīng)過(guò)短暫的變形后在最小表面自由能作用下最終合并[31,32].然而,研究發(fā)現(xiàn)在Leidenfrost 狀態(tài)下,油面上兩個(gè)同等大小的FC-72液滴在蒸發(fā)的過(guò)程中相互吸引但永不合并,這一新穎現(xiàn)象被稱為液滴的非聚合態(tài)蒸發(fā).為了探索液滴間蒸發(fā)的機(jī)理,對(duì)To=88.2 ℃的油池表面兩個(gè)液滴的運(yùn)動(dòng)行為進(jìn)行量化,顯示兩滴最小間隙δ隨時(shí)間的變化.圖4(b) 顯示了蒸發(fā)的液滴依次經(jīng)歷了3 個(gè)階段: 接觸非合并蒸發(fā) (non-coalescence)、彈跳蒸發(fā)(bouncing)和相互分離蒸發(fā) (separation).輕輕放在油面的兩個(gè)液滴在重力勢(shì)能的作用下,先相互吸引隨后開(kāi)啟液滴間的衰減式彈跳過(guò)程.持續(xù)180 ms 后,彈跳停止,液滴轉(zhuǎn)為安靜的非合并蒸發(fā)階段.兩個(gè)液滴間的氣膜夾層阻礙了它們聚合互溶.然而,受限于液滴“接觸”界面擾動(dòng)影響,高速相機(jī)無(wú)法直接捕獲這層氣膜.所以,這一階段所測(cè)量到的液滴最小間隙δ=0.直到t=2.8 s (約為蒸發(fā)壽命的29%),這種安靜的蒸發(fā)狀態(tài)被打破,液滴間的彈跳過(guò)程自動(dòng)開(kāi)啟.隨著液滴變小,彈跳發(fā)生的頻率從低到高,逐漸趨于穩(wěn)定.圖4(c)展示了液滴在5217—5342 ms 時(shí)間段內(nèi)穩(wěn)定的彈跳行為,周期和振幅分別為32 ms 和28 μm.液滴彈跳階段的時(shí)長(zhǎng)為3.11 s,約占蒸發(fā)壽命的32%.最后,兩個(gè)液滴在經(jīng)歷了短暫的平靜蒸發(fā)后,緩慢地分離開(kāi).在t=7.675 s 時(shí)刻,液滴間的最小間隙δ達(dá)到711 μm,遠(yuǎn)高出彈跳蒸發(fā)階段2 個(gè)數(shù)量級(jí).分離后的液滴各自蒸發(fā),直至壽命結(jié)束.

圖4 (a) 從a—i 的9 個(gè)特定時(shí)間的液滴動(dòng)態(tài)圖;(b),(c) 在To=88.2 ℃時(shí),液滴整個(gè)壽命期的3 種蒸發(fā)行為Fig.4.(a) Droplet dynamics at nine specific time from a to i;(b),(c) three-regimes behavior of droplet dynamics during the whole droplet lifetime at To=88.2 ℃.

3.3 Leidenfrost 液滴的蒸發(fā)干涉效應(yīng)

液滴的蒸發(fā)速率受液滴數(shù)量、陣列形式和液滴間距的影響[33-36].相鄰的液滴互相影響蒸發(fā)速度,且隨著液滴距離的減小,液滴間的干涉作用加劇[37].對(duì)于液滴準(zhǔn)穩(wěn)態(tài)蒸發(fā)的情況,傳質(zhì)通量的解析表達(dá)式為 [38]

式中,?RH為相對(duì)濕度,ρv為蒸氣密度,R為液滴半徑,f為與液滴和基底接觸角相關(guān)的形狀因子,dv為質(zhì)量擴(kuò)散系數(shù).當(dāng)f和dv一定的情況下,液滴所處環(huán)境濕度的增大會(huì)減緩液滴蒸發(fā).反之,則加速液滴蒸發(fā).干涉蒸發(fā)體現(xiàn)在鄰近液滴之間局部增大的蒸氣濃度/濕度,它使蒸發(fā)的驅(qū)動(dòng)力減小[39,40].對(duì)于陣列分布的液滴,干涉效應(yīng)使最外層液滴的蒸發(fā)速率普遍大于陣列內(nèi)部液滴[41].因此,Leidenfrost狀態(tài)下,液滴相互干涉的蒸發(fā)機(jī)制不能忽視.

圖5 所示為液滴表面的溫度等高線分布,Leidenfrost 狀態(tài)相鄰的兩個(gè)液滴外側(cè)區(qū)域相比二者之間區(qū)域溫度有明顯的下降,且整體來(lái)看兩個(gè)液滴具有相似的溫度分布特點(diǎn).以中心線左側(cè)液滴為例,t=5248 ms 時(shí)刻,液滴遠(yuǎn)離中心線一側(cè)的表面溫度梯度? T/?n=5.6 ℃/pix,明顯大于靠近中心線側(cè)的? T/?n=2.2 ℃/pix,n為矢量且垂直于液滴表面.事實(shí)上,液相和氣相之間的相變涉及質(zhì)量和熱量的傳遞.兩種傳遞方式通過(guò)路易斯數(shù)Le=a/dv關(guān)聯(lián),其中a為熱擴(kuò)散系數(shù).該標(biāo)度律顯示了蒸氣濃度梯度? Cv/?n和邊界處溫度梯度? T/?n的正相關(guān)特性[42].換言之,? T/?n越大,蒸發(fā)的驅(qū)動(dòng)力越大,液滴蒸發(fā)越快.因此,綜合前述討論,液滴兩側(cè)溫度分布的不對(duì)性易于被解釋.相鄰液滴蒸發(fā)干涉效應(yīng)導(dǎo)致液滴局部區(qū)域的溫度梯度/濃度梯度改變,由此產(chǎn)生了液滴表面蒸發(fā)速率不同,這決定了表面溫度的非對(duì)稱性分布.具體說(shuō),蒸發(fā)速率越快,相變帶走的潛熱越多,液滴該區(qū)域的表面溫度越低.需要指出,由于獲得溫度梯度是相對(duì)值,因此上述傳熱與傳質(zhì)的類比是可行的,它同Yu 等[43]所采用的方法一致.

圖5 雙滴接觸時(shí)刻,液滴區(qū)域和背景區(qū)域的溫度分布Fig.5.Temperature distribution on the droplet-dominated and background-dominated regions at contact time.

提取9 個(gè)特殊時(shí)刻的液滴真實(shí)溫度輪廓線,如圖6 所示.首先,在液滴赤道平面建立極坐標(biāo)系,A表示液滴距離最近點(diǎn),B代表距離最遠(yuǎn)點(diǎn),φ為以A為起始點(diǎn)的水平坐標(biāo)偏轉(zhuǎn)角.當(dāng)液滴接觸,表面變形使得φA＞ 0.對(duì)于不同蒸發(fā)階段,液滴赤道面B點(diǎn)兩側(cè)均表現(xiàn)出對(duì)稱的周向溫度分布.如前所述,局部強(qiáng)烈的蒸發(fā)冷卻作用使B點(diǎn)處溫度始終最低.當(dāng)液滴與油池接觸,過(guò)冷的液滴使圓周溫度均勻分布,此時(shí)? T/?φ ≈0 .隨著液滴溫度升高,非對(duì)稱性蒸發(fā)現(xiàn)象出現(xiàn),周向溫度梯度? T/?φ逐漸被建立.此時(shí),相鄰液滴的內(nèi)側(cè)溫度也開(kāi)始低于外側(cè).當(dāng)液滴進(jìn)入蒸發(fā)的第二階段—彈跳蒸發(fā),相鄰液滴干涉引起的非對(duì)稱蒸發(fā)加強(qiáng).t=5248 ms時(shí)刻,液滴內(nèi)外兩側(cè)的最大溫差升高到 5 ℃.進(jìn)入彈跳蒸發(fā)階段后期,液滴尺寸減小又會(huì)引起內(nèi)部Marangoni 流反向影響液滴溫度,使其變得均勻.因此,該時(shí)期液滴表面溫度的差異減小,直到液滴徹底分離,? T/?φ近似為0.

圖6 紅外相機(jī)俯拍得到液滴表面的溫度輪廓線 (溫度隨著圓周角變化)Fig.6.Droplet surface temperatures along top view circle(temperatures are plotted versus circumference angles).

如圖6 所示,需要關(guān)注一個(gè)重要現(xiàn)象,即3 個(gè)代表時(shí)刻,液滴 “接觸” (5248 ms)、分離 (5257 ms)和“再接觸” (5265 ms),周向溫度最大溫差升高-下降-再升高的趨勢(shì),分別對(duì)應(yīng)溫差為5.29,1.56 和5.06 ℃.總的來(lái)說(shuō),液滴表面溫差與液滴水平彈跳行為存在著明顯的協(xié)同關(guān)系.彈跳蒸發(fā)方式隨著液滴表面溫差增大而出現(xiàn),隨溫差的消失而結(jié)束.因此,在Leidenfrost 狀態(tài)下,液滴干涉蒸發(fā)效應(yīng)由于改變液滴表面溫度分布,同液滴非聚合態(tài)蒸發(fā)的動(dòng)力學(xué)特性存在著聯(lián)系.液滴表面溫度的變化會(huì)影響表面張力大小,從而可能打破液滴的受力平衡.

4 干涉作用下液滴的受力分析

為進(jìn)一步探索加熱池上兩個(gè)Leidenfrost 液滴3 種蒸發(fā)狀態(tài)轉(zhuǎn)變的物理機(jī)制,對(duì)液滴的運(yùn)動(dòng)建立了一個(gè)簡(jiǎn)單的模型.如圖7(a)所示,模型包括4 種作用力,即潤(rùn)滑推動(dòng)力Fp、Marangoni 水平分力Fσ,x、重力水平分力Fg,x和黏性阻力Fd.需強(qiáng)調(diào)的是,上述沿x軸的水平分力為廣義方向上力的分量,實(shí)際上受液滴重力作用,凹陷的油面使x坐標(biāo)軸與油面水平基線(level line)存在夾角α.油面水平基線取液滴背景的油面等高線(見(jiàn)圖7(a)),這同文獻(xiàn) [44]中方法一致.下文中“水平”均為廣義上稱呼,不再解釋.顯然,Fp和Fσ,x利于液滴分離,而Fg,x和Fd則阻礙液滴分離(見(jiàn)圖7(a)).

液滴底部蒸氣層內(nèi)自發(fā)的和恒定的對(duì)稱性破壞不僅是維持Leidenfrost 液滴懸浮于液池之上的關(guān)鍵,同時(shí)也為液滴水平方向運(yùn)動(dòng)提供動(dòng)力.這種不對(duì)稱的氣膜結(jié)構(gòu)部分改變了蒸氣的流通方向,引起液滴前后黏滯力的差異,從而使液滴自發(fā)推進(jìn)成為可能.如圖7(b)所示,黏滯力差可以通過(guò)乘以非對(duì)稱結(jié)構(gòu)因子Δε/ε得到:其中Δε是氣膜出口理論厚度差值,ε是氣膜的特征厚度,μv是蒸氣黏度,uv是氣膜內(nèi)泊肅葉流的典型流速.此外,氣膜潤(rùn)滑理論給出氣膜內(nèi)壓差力與液滴重力的守恒:(ρd是液滴密度).聯(lián)合兩個(gè)式子,得出潤(rùn)滑推動(dòng)力[14,15]:

圖7 油面雙滴非聚合蒸發(fā)的力學(xué)分析 (a) 施加在傾斜液滴上的各種力;(b) 潤(rùn)滑推動(dòng)力;(c) Leidenfrost 蒸氣層與空氣交界處的Marangoni 力;(d) 氣膜出口位置的液滴表面張力圓周分布;(e) 沿x 方向分布的液滴表面溫度;(f) 作用在傾斜油面上的液滴重力;(g) 油面傾斜角α 與液滴邦德數(shù)Bo 之間的關(guān)系Fig.7.Force analysis explaining the non-coalescence phenomenon of dual-droplets over oil surface: (a) Various forces exerted on inclined droplets;(b) lubrication-propulsion force;(c) Marangoni force along circumference direction at the junction between Leidenfrost vapor layer and air;(d) distribution of droplet surface tension along the x-direction at the vapor outlet;(e) distribution of droplet surface temperature along the x-direction;(f) droplet gravity on the inclined oil surface;(g) the relationship between the oil surface’s inclination angle α and the droplet’s bond number.

氣膜兩側(cè)理論厚度差值 Δε ≈Dcβ=2Rsin(θ)β,Dc為液滴在蒸氣出口處圓截面的直徑,θ表示氣膜出口位置偏離豎直方向的角度,β是Δε存在所產(chǎn)生的角度,β=4μvVmR2/(σdε2)[8].Vm表示液滴內(nèi)部流速,對(duì)于Leidenfrost 液滴,Vm約為1 cm/s[8].平衡液面彎曲力 2 πσoRsin2θ與液滴重力mg=4/3πR3ρdg并聯(lián)合Δε表達(dá)式后,可以得到:

其中,ρo是硅油密度,σd是液滴表面張力,Ts為液面溫度,lo為硅油的毛細(xì)長(zhǎng)度,lo=[σo/(ρog)]1/2=0.15 mm.本文應(yīng)用潤(rùn)滑理論預(yù)測(cè)氣膜特征厚度ε表達(dá)式為[14]

其中,L為液滴相變潛熱,λv為蒸氣的導(dǎo)熱系數(shù).To—Td表示油面與液滴間溫差,Td取液滴平均溫度.由此可計(jì)算氣膜兩側(cè)理論厚度差值Δε,將計(jì)算結(jié)果與Matsumoto 等 [15]報(bào)道比較,重合的液滴半徑范圍內(nèi)0.75—0.50 mm,本文獲得Δε對(duì)應(yīng)0.7—0.9 μm,符合文獻(xiàn)預(yù)測(cè).因此,驗(yàn)證了關(guān)于氣膜厚度的建模方法.聯(lián)合(3)式—(5)式,Leidenfrost液滴潤(rùn)滑推動(dòng)力最終表示為

表面張力是溫度的函數(shù).氣液界面存在的溫度梯度會(huì)在界面引起張力梯度,從而產(chǎn)生Marangoni力.該力是氣-液界面的宏觀作用力,并驅(qū)使液滴或氣泡的運(yùn)動(dòng)[45].本文借鑒文獻(xiàn)[14-16,46]對(duì)Leidenfrost 液滴自驅(qū)動(dòng)的建模計(jì)算過(guò)程中處理液滴內(nèi)部流場(chǎng)的方法,即由流場(chǎng)不均勻而產(chǎn)生的力認(rèn)為是液滴內(nèi)力,忽略其對(duì)液滴運(yùn)動(dòng)的影響.如圖7(c)所示,對(duì)左側(cè)的液滴選取圓心角為dφ的微元Δs,則它在水平方向受到的Marangoni 力為

其中σφ是指以A′為起點(diǎn)逆時(shí)針旋轉(zhuǎn)至φ處氣-液界面的張力,可參考圖7(d),σφ表達(dá)式為

式中,Ts表示FC-72 表面溫度,σΑ′表示液滴在A′點(diǎn)的表面張力.dσ/dTs代表單位溫升的FC-72表面張力變化量,它通過(guò)實(shí)驗(yàn)前擬合不同溫度的FC-72 液體表面張力給出.σ=(17.096-0.098Ts)×10-3,相關(guān)性系數(shù)為0.9935.dTs/dx為液滴在氣膜出口位置水平方向的溫度梯度,見(jiàn)圖7(c).注意到液滴溫度沿水平方向呈線性分布,如圖7(e)所示,一個(gè)代表時(shí)刻t=5428 ms,其中x?=x/D=0.5(1-cosφ),代表歸一化后的x坐標(biāo).假設(shè)來(lái)自鄰側(cè)液滴的干涉效應(yīng)僅在x方向變化,而不受液高 (A′ 到A方向)的影響,則dTs/dx僅為時(shí)間的函數(shù),不隨x變化.因此 dTs/dx可以通過(guò)擬合液滴赤道的圓周溫度得出.聯(lián)同(7)式和(8)式,并對(duì)?在 [0,2π]的范圍積分得到

在蒸發(fā)期間,液滴兩側(cè)的油面表現(xiàn)出明顯的高度差別,即外側(cè)高而內(nèi)側(cè)低(見(jiàn)圖7(a)).因此對(duì)于左側(cè)液滴來(lái)說(shuō),存在重力的水平分力,

其中,α表示雙滴重力作用下傾斜油面產(chǎn)生的夾角,α ?β～0.1°,β表示液滴與油池表面角度偏移量的相對(duì)值.圖7(g)顯示了液滴大小對(duì)油面絕對(duì)傾角α的影響,顯然,α的是液滴邦德數(shù)Bo=ρdR2g/σd的函數(shù),α=3.08Bo0.53,相關(guān)性系數(shù)為 0.981.

當(dāng)液滴運(yùn)動(dòng)時(shí),由流體力學(xué)的經(jīng)典理論可知,氣膜內(nèi)產(chǎn)生的二次庫(kù)埃特流將阻礙液滴移動(dòng),這就產(chǎn)生了黏性阻力Fd.該力的大小與蒸氣層與液滴接觸面在垂直方向上的投影截面積、液滴運(yùn)動(dòng)速度與蒸氣膜厚度之比u/ε有關(guān)[14,15],表示為Fd～液滴速度u取其在液滴彈跳階段的平均速度～0.1 cm/s.

圖8(a)顯示了液滴在“接觸”和“分離”時(shí)表面溫差 ΔTs=TA-TB隨時(shí)間的變化.對(duì)于“接觸”的液滴,tc≈5.6 s 時(shí)刻,ΔTs達(dá)到最大值6 ℃.在此之前和之后,ΔTs分別呈增大和減小趨勢(shì).對(duì)其隨時(shí)間關(guān)系擬合得到: ΔTs=0.88+4.9 exp[-0.5(t-tc)2/w2],相關(guān)性系數(shù)為0.970.當(dāng)t ＞tc時(shí),w=2.74.t ＜tc時(shí),w=0.81;在彈跳蒸發(fā)階段,ΔTs隨著液滴短暫的“分離”降低,當(dāng)ΔTs約等于1.03 ℃時(shí),分離過(guò)程結(jié)束,吸引過(guò)程開(kāi)始.綜合液滴的受力計(jì)算,得到液滴蒸發(fā)過(guò)程中4 種作用力瞬時(shí)值.圖8(b)對(duì)比了這4 種作用力的大小,以判定主要作用力.結(jié)果表明,在0—8 s 的蒸發(fā)時(shí)間內(nèi)Fσ,x,Fg,x,Fp以及Fd對(duì)應(yīng)的量級(jí)范圍分別為 10—7—10—8N,10—6—10—9N,10—8—10—12N 和10—10—10—12N.顯然,Fσ,x與Fg,x間的抗衡主導(dǎo)液滴蒸發(fā)在3 種狀態(tài)間轉(zhuǎn)變.重力Fσ,x趨向于將液滴維持在非聚合蒸發(fā)階段.Marangoni 力Fσ,x則趨向使液滴進(jìn)入水平彈跳和分離蒸發(fā)階段.相比而言,不再考慮潤(rùn)滑推動(dòng)力Fp以及黏性阻力Fd對(duì)液滴運(yùn)動(dòng)的影響.

Fσ,x和Fg,x在圖8(c)中顯示,縱坐標(biāo)取各種力的對(duì)數(shù)值.由于Fσ,x和Fg,x隨液滴大小的變化速率不同,因此兩個(gè)力間存在兩個(gè)交點(diǎn),分別對(duì)應(yīng)時(shí)間tA和tB.其中,tA代表油池上Leidenfrost 液滴由非聚合蒸發(fā)向彈跳蒸發(fā)階段轉(zhuǎn)變的時(shí)間,在此之前Fg,x＞Fσ,x,重力主導(dǎo)液滴運(yùn)動(dòng),液滴被相互吸引.TB代表液滴開(kāi)始進(jìn)入完全分離的階段,此后Fg,x＜Fσ,x,Marangoni 力主導(dǎo)液滴運(yùn)動(dòng),使液滴脫離吸引反向移動(dòng).在tA和tB之間時(shí)間段內(nèi),兩種力間相互競(jìng)爭(zhēng).兩個(gè)插圖有助于理解這種競(jìng)爭(zhēng)所主導(dǎo)的液滴間彈跳蒸發(fā)的行為.當(dāng)液滴靠近時(shí),雙液滴干涉引起非對(duì)稱性蒸發(fā)產(chǎn)生的Fσ,x大于油面傾斜所產(chǎn)生的Fg,x,凈力使液滴分離.但隨液滴互相分離,雙滴間減弱的干涉作用使蒸發(fā)逐漸對(duì)稱.當(dāng)Fσ,x無(wú)法繼續(xù)抗衡Fg,x時(shí),液滴重新被吸引,開(kāi)始相向移動(dòng).這種周期彈跳的行為直到tB時(shí)刻結(jié)束.得益于小液滴作用的油面曲率減小,在tB時(shí)刻,減小的油面傾角α使Leidenfrost 液滴更像在平板表面蒸發(fā)(見(jiàn)圖7(f)).tB之后,Fg,x永遠(yuǎn)無(wú)法與之抗衡Fσ,x,進(jìn)而進(jìn)入下一蒸發(fā)階段.需要強(qiáng)調(diào),這里液滴重力和Marangoni 力的競(jìng)爭(zhēng)機(jī)制主導(dǎo)液滴的運(yùn)動(dòng)方式,有別于毛細(xì)長(zhǎng)度ld的物理意義,即流體所受表面張力與重力的比值,它是判斷液滴是否容易形變的重要依據(jù).一個(gè)簡(jiǎn)單的例子,毛細(xì)長(zhǎng)度僅能約束高空下落液滴的形狀,而無(wú)法改變其向下運(yùn)動(dòng)的趨勢(shì),即使液滴足夠小.

圖8 液滴溫度和各種力的變化 (a) 測(cè)量的液滴溫度與時(shí)間的關(guān)系,以及用簡(jiǎn)單的擬合得到的兩條曲線;(b) 各種力大小的比較;(c) 競(jìng)爭(zhēng)的重力Fg,x 和Marangoni 力Fσ,x 主導(dǎo)了液滴動(dòng)力學(xué)的三態(tài)行為Fig.8.Variation of droplet temperatures and various forces: (a) The measured droplet temperature versus time and two curve obtained with simple fitting;(b) comparison of various forces magnitudes;(c) competing gravity force Fg,x and Marangoni force Fσ,x dominate the three-regimes behavior of droplet dynamics.

由此得到的液滴3 個(gè)蒸發(fā)階段的轉(zhuǎn)化時(shí)間tA=3.245 s 和tB=6.110 s,分別對(duì)應(yīng)實(shí)驗(yàn)值2.979 s和6.087 s.最大時(shí)間偏差為0.266 s,它約為液滴壽命的3%.較小的實(shí)驗(yàn)與計(jì)算值之間的偏差,證實(shí)了Leidenfrost 液滴重力與Marangoni 力競(jìng)爭(zhēng)的機(jī)制主導(dǎo)兩個(gè)液滴非聚合態(tài)蒸發(fā).

5 結(jié)論

本工作對(duì)高溫油面上相鄰液滴的Leidenfrost蒸發(fā)機(jī)制進(jìn)行了研究.與常溫液滴相互接觸時(shí)呈現(xiàn)立即或延遲互溶不同,高溫液面上蒸發(fā)的兩個(gè)等大的FC-72 液滴在宏觀尺度下表現(xiàn)為互相吸引但永不融合.這種現(xiàn)象被稱為L(zhǎng)eidenfrost 液滴的非聚合蒸發(fā).該蒸發(fā)方式下液滴半徑仍然遵循R(t) — (1 —t/τ)n定律.但指數(shù)系數(shù)n較液滴單獨(dú)蒸發(fā)有所增大,相鄰液滴的蒸發(fā)速率變慢.鑒于此,在液滴Leidenfrost 蒸發(fā)過(guò)程中多滴之間蒸發(fā)干涉效應(yīng)依然存在.液滴的非對(duì)稱性蒸發(fā)使液滴表面溫度周向不均勻性逐漸建立,最大溫差達(dá)6 ℃.

從微觀尺度上研究發(fā)現(xiàn),蒸發(fā)雙滴間依次經(jīng)歷彼此相互吸引的非聚合蒸發(fā)階段、彈跳蒸發(fā)階段和分離階段.對(duì)液滴所受的各種力分析,結(jié)果表明液滴蒸發(fā)過(guò)程中存在著兩個(gè)互相關(guān)聯(lián)的動(dòng)態(tài)機(jī)制: 一方面是蒸發(fā)干涉引起Marangoni 力誘導(dǎo)液滴間分離;另一方面是池面傾斜產(chǎn)生液滴重量分力誘導(dǎo)液滴間吸引.且兩種力均與液滴大小有關(guān),它們共同影響著雙滴蒸發(fā).所得液滴3 個(gè)蒸發(fā)階段間轉(zhuǎn)化時(shí)間tA=3.245 s 和tB=6.110 s,分別與實(shí)驗(yàn)值2.979 s和6.087 s 相對(duì)應(yīng).最大時(shí)間偏差0.266 s,約為蒸發(fā)壽命的3%,相鄰Leidenfrost 液滴的相互作用及運(yùn)動(dòng)機(jī)制得到了解釋.