液滴撞擊固體壁面是一種常見(jiàn)的物理現象，廣泛應用于噴墨打印、等離子噴涂、噴霧冷卻、金屬淬火和抑冰等領(lǐng)域。液滴與固體壁面的相互作用是影響上述應用的重要因素，如墨滴撞擊紙張的鋪展特性會(huì )影響噴墨打印的分辨率，液滴與壁面的接觸時(shí)間會(huì )影響噴霧冷卻的效率，而抑冰又要求減少液滴與壁面的接觸時(shí)間。類(lèi)似的情況還有許多，盡管很多學(xué)者已經(jīng)開(kāi)展了大量研究，但液滴與壁面的相互作用仍因其復雜性與多樣性而未被全部知曉。

對于液滴撞擊固體壁面，通常從兩方面展開(kāi)研究，即液滴的鋪展特性與傳熱傳質(zhì)特性。學(xué)者們圍繞這兩方面進(jìn)行了一定研究。隨著(zhù)越來(lái)越多的動(dòng)力學(xué)及傳熱學(xué)特性被發(fā)現，研究工況的溫度也隨之升高，其中一個(gè)研究方向逐漸聚焦于高溫情況下顯著(zhù)的液滴撞擊特征，Leidenfrost現象。1756年，Leidenfrost提出了這一標志性的液滴撞擊高溫壁面特征，開(kāi)創(chuàng )了新的研究領(lǐng)域。Leidenfrost現象是液滴受到高溫壁面的影響，氣液快速交換熱量，液滴底部產(chǎn)生蒸汽膜的物理現象。蒸汽膜厚度受多種因素影響而不斷發(fā)生變化，并存在程度不一的振蕩情況。實(shí)驗結果表明，Leidenfrost現象中液滴底部蒸汽為徑向的層流流動(dòng)，并有學(xué)者對蒸汽層的擴散提出了相關(guān)預測模型。

除底部氣膜外，Leidenfrost現象下液滴的動(dòng)力學(xué)與傳熱特性也與其它工況不同。首先，達到Leidenfrost狀態(tài)需要有足夠高的溫度，這一溫度通常被稱(chēng)為L(cháng)eidenfrost點(diǎn)，其大小受到壁面粗糙度，重力和表面張力等因素的影響。其次，液滴撞擊壁面的We數較大時(shí)，相對于未加熱壁面與濕壁面，Leidenfrost狀態(tài)下液滴更容易發(fā)生破碎。當液滴開(kāi)始鋪展時(shí)，由于液滴與壁面間的蒸汽膜存在，壁面粗糙度與接觸角對液滴鋪展的影響被削弱，這使得Leidenfrost狀態(tài)下的液滴鋪展情況大于膜態(tài)蒸發(fā)和核態(tài)沸騰狀態(tài)下的液滴鋪展情況，且液滴的最大鋪展情況與脫離時(shí)間仍然由We數決定。

雖然，很多學(xué)者對Leidenfrost現象進(jìn)行了一定研究，但是大部分研究成果聚焦于高溫平面上的液滴動(dòng)力學(xué)與傳熱特性，而非平面上Leidenfrost現象的研究非常有限。但是，液滴撞擊高溫曲面廣泛存在于電子工程、機械工程、動(dòng)力工程、化學(xué)工程等領(lǐng)域。如液滴撞擊高溫球面，常見(jiàn)于流化催化裂化和流化焦化等工藝，這一過(guò)程影響著(zhù)多種過(guò)程系統的性能，尤其是對于生產(chǎn)更有價(jià)值的汽車(chē)燃料具有重要意義。鑒于此，本文進(jìn)一步探究液滴與高溫球面撞擊過(guò)程中液滴的Leidenfrost現象，采用數值模擬方法進(jìn)行研究，使用Fluent軟件構建了二維軸對稱(chēng)數值模型，采用VOF模型及Lee模型，分析液滴撞擊高溫球面的動(dòng)力學(xué)特性與傳熱特性，探討Leidenfrost現象的機理及各因素的影響規律。

1.模型及方法

1.1物理模型

液滴撞擊球面的物理模型如圖1所示，半徑為r的球狀液滴以初速度v撞擊半徑為rn的固體球面，并在球面上鋪展、收縮和反彈。液滴材料為水，考慮重力加速度的影響，液滴與壁面撞擊時(shí)t=0.

圖1物理模型

定義無(wú)量綱系數鋪展因子β

式中：rm為液滴鋪展半徑，mm.鋪展系數β可以直觀(guān)地反映液滴在壁面上鋪展變化的規律。鋪展系數越大，說(shuō)明液滴鋪展越充分。

與撞擊平面不同，液滴撞擊球面時(shí)會(huì )沿球面鋪展，因此鋪展半徑rm的定義與平面有所不同

（2）

式中，φ為液滴鋪展角度，rd為液滴鋪展最大處與球心間的距離。

1.2控制方程

為了模擬液滴撞擊熱壁面的動(dòng)力學(xué)行為，需要求解連續性方程、動(dòng)量方程與能量方程，假設氣液兩相均為不可壓縮流體，表面張力系數、密度、粘度、導熱系數和比熱等物性參數均為常數

式中，u為速度矢量；ρ為流體密度；S為相變引起的質(zhì)量源項；μ是流體黏度；g為重力加速度；Fσ為表面張力；T為流體溫度；Cp為定壓比熱；k為導熱系數；Se為相變引起的能量源項。

1.3 VOF模型

VOF模型是一種兩相流動(dòng)計算模型，其通過(guò)計算網(wǎng)格中相的體積分數以捕捉兩相的分界面，對于目標流體，流體體積函數C可表示為

流體體積函數滿(mǎn)足的輸運方程為

采用隱式格式求解流體體積函數方程，并包含隱式體積力，界面建模為sharp模型，采用的界面捕捉算法為Compressive方法。

兩相流體的密度ρ、黏度μ、導熱系數k的計算式為

式中，下標l、v分別表示液相和氣相。

采用CSF模型計算表面張力，表達式為

式中，σ為表面張力系數；n為界面法向；κ為界面曲率。

1.4 Lee模型

本文采用Lee模型模擬氣液相變過(guò)程，控制方程為

式中，αv為氣相體積分數，uv為氣相速度，

和為蒸發(fā)和冷凝的傳質(zhì)率。對于蒸發(fā)過(guò)程，傳質(zhì)表達式為

式中，αl為液相的體積分數；Tl為流體溫度；Tsat為飽和溫度；Acoeff為控制相變強度的因子，可以表示為

式中，db為液體蒸發(fā)時(shí)產(chǎn)生的氣泡直徑；χ表示液體表面吸附蒸氣分子的量，為調節系數；M為摩爾質(zhì)量；R為通用氣體常數，8.314 J/（mol·K）；L為蒸發(fā)潛熱；在本文中，Acoeff取值為500.

1.5網(wǎng)格無(wú)關(guān)性驗證

為了選擇合理的網(wǎng)格精度，首先開(kāi)展網(wǎng)格無(wú)關(guān)性驗證。同時(shí)，為了驗證數值模擬的準確性，也需與實(shí)驗結果進(jìn)行對比驗證。本文以文獻中的工況為基準，采用異丙醇作為液相，展開(kāi)網(wǎng)格無(wú)關(guān)性驗證與模型驗證。計算區域的尺寸為5 mm×8 mm,液滴初始溫度為355.65 K、初始直徑為1.61 mm、液滴撞擊速度為0.348 m/s、We數等于8.6、球體直徑為10 mm、壁面溫度為623.15 K、壁面接觸角為67°、重力加速度為?9.81 m/s2.其它流體物性如表1所示。

表1流體物性參數表

為了定量分析網(wǎng)格數量對數值模擬結果準確性的影響，以液滴撞擊壁面時(shí)的鋪展因子作為定量分析指標及網(wǎng)格無(wú)關(guān)性驗證的評判標準。網(wǎng)格數量從36 w開(kāi)始增加，每次增長(cháng)大約30 w,并繪制液滴鋪展因子β隨時(shí)間變化與網(wǎng)格數量的關(guān)系圖。如圖2所示，網(wǎng)格數在86 w以上時(shí)基本不會(huì )對模擬結果產(chǎn)生影響，為保證模擬精度并確保求解時(shí)間，本文采用86 w網(wǎng)格數作為模型的基準網(wǎng)格數進(jìn)行模擬。從圖中也可以看出，數值模擬結果與參考文獻中的實(shí)驗結果趨于一致，說(shuō)明本文采用的模型可以用于研究液滴撞擊球面的Leidenfrost現象。

圖2網(wǎng)格無(wú)關(guān)性驗證和實(shí)驗驗證

圖3進(jìn)一步給出了不同時(shí)刻液滴撞擊球面的形態(tài)變化，圖3（a）為二維模擬結果、圖3（b）為文獻的實(shí)驗結果、圖3 c為paraview處理后的三維化模擬結果。如圖所示，0 ms時(shí)刻，液滴以0.348 m/s的速度撞擊球面，隨后液滴由慣性力驅動(dòng)開(kāi)始鋪展，鋪展過(guò)程中動(dòng)能不斷轉化為表面能和耗散能。1.5 ms-4.5 ms時(shí)段，液滴由球狀逐漸轉變?yōu)閳A盤(pán)狀。達到最大鋪展系數后，表面張力促使液滴回縮，6.0 ms-7.5 ms時(shí)段，液滴逐漸回縮并開(kāi)始反彈，上下部分出現較明顯分界，整體呈葫蘆狀。10.0 ms時(shí)，液滴與壁面的接觸面積明顯減少。15.0 ms時(shí)，液滴與壁面完全分離?？梢钥闯?，模擬結果與實(shí)驗情況吻合良好，進(jìn)一步證明本文數值模擬的準確性。

圖3模擬結果，（a）二維模擬結果（b）實(shí)驗結果（c）三維處理