黃平，卜雪琴,*，劉一鳴，林貴平，楊坤

1.北京航空航天大學(xué) 航空科學(xué)與工程學(xué)院，北京 100083

2.中國航發(fā)商用航空發(fā)動機有限責(zé)任公司, 上海 200241

隨著中國航空航天事業(yè)的高速發(fā)展，民用飛機和大型運輸機更加注重經(jīng)濟性、安全性、舒適性。飛機/發(fā)動機結(jié)冰及防除冰是飛機/發(fā)動機設(shè)計及適航驗證中面臨的關(guān)鍵問題。飛機在含有過冷水滴或冰晶的云層中飛行時可能面臨嚴重的結(jié)冰現(xiàn)象，使得飛機的安全性受到威脅，飛機需要安裝防除冰系統(tǒng)以保障飛行安全。飛機/發(fā)動機結(jié)冰和防除冰是一個復(fù)雜的多相流、傳熱、傳質(zhì)等多物理場耦合問題。隨著人們對飛機/發(fā)動機結(jié)冰事故發(fā)生原因的猜測進而確認，研究者對飛機結(jié)冰研究的熱點也在不斷地發(fā)展，以了解結(jié)冰的機理、影響和危害。

本文介紹了飛機結(jié)冰研究的發(fā)展，進而引申到目前國內(nèi)外的研究熱點混合相/冰晶結(jié)冰及其對飛機/發(fā)動機的影響；總結(jié)了混合相/冰晶結(jié)冰與過冷水結(jié)冰的差異以及對飛機的危害；從冰風(fēng)洞試驗和數(shù)值模擬方面總結(jié)了混合相/冰晶結(jié)冰的研究進展；最后分析了混合相/冰晶結(jié)冰研究的不足及展望。

1 結(jié)冰研究的發(fā)展

從現(xiàn)代飛機飛行開始，飛機結(jié)冰問題就逐漸引起了人們的重視。最早報道的結(jié)冰事故可以追溯到20世紀20年代。首先引起人們關(guān)注的是過冷水結(jié)冰。1994年，美國印第安納州發(fā)生一起嚴重的空難事故，一架載有68人的ATR-72飛機發(fā)生墜機，引起事故的原因歸結(jié)于過冷大水滴(直徑大于50 μm)結(jié)冰。自此,過冷大水滴結(jié)冰開始倍受關(guān)注。

20世紀90年代初,有關(guān)高空冰晶引發(fā)的商用飛機發(fā)動機失效事件引起了人們的關(guān)注。Mason等對1990—2006年的240起結(jié)冰事故進行了總結(jié)，發(fā)現(xiàn)其中有62起是因冰晶結(jié)冰引起渦扇發(fā)動機動力損失。自此,有關(guān)冰晶結(jié)冰的研究也開始逐漸展開。2011年NASA梳理了發(fā)動機高冰-水含量(High Ice-Water Content, HIWC)條件下的結(jié)冰研究框架并開展相關(guān)研究。2012年歐盟成立了高空冰晶(High Altitude Ice Crystals，HAIC)研究項目。2014年美國聯(lián)邦航空規(guī)章FAR25-140和33-34號修正案將過冷大水滴和冰晶適航納入其中。

NASA和歐盟對于冰晶結(jié)冰研究的框架圖分別如圖1和圖2所示。其中歐盟的HAIC項目聯(lián)合了11個歐洲國家的34家航空相關(guān)單位、澳大利亞的5家單位、加拿大以及美國。從NASA和歐盟的冰晶結(jié)冰框架圖來看，主要集中在冰晶結(jié)冰機理與危害、地面試驗與測量技術(shù)、空中探測與試驗技術(shù)、冰晶氣象研究及預(yù)報技術(shù)、冰晶結(jié)冰數(shù)值模擬技術(shù)等5個方面。

圖1 NASA冰晶積冰研究項目框架圖[11]Fig.1 NASA High Ice-Water Content program[11]

圖2 歐盟HAIC項目技術(shù)框架圖[12]Fig.2 European High Altitude Ice Crystals program[12]

國內(nèi)高等院校、航空工業(yè)部門等近幾年來也逐漸意識到冰晶結(jié)冰對航空發(fā)動機的危害。袁慶浩等針對航空發(fā)動機內(nèi)部冰晶結(jié)冰研究進行了調(diào)研和總結(jié)，從地面模擬試驗、計算模擬兩方面論述了對冰晶機理的研究進展。沈浩等從數(shù)值和試驗兩方面對國內(nèi)外已開展的壓氣機內(nèi)結(jié)冰研究進行了分析與總結(jié)，討論了當(dāng)前冰晶結(jié)冰的主要研究成果和存在的局限，提出了開展冰晶結(jié)冰研究需重點關(guān)注的方向。姜飛飛等采用拉格朗日方法分析了冰晶在發(fā)動機內(nèi)涵道的運動和融化過程，得到了冰晶在低壓壓氣機內(nèi)涵通道內(nèi)運動的軌跡及與葉片碰撞的特性。卜雪琴等將Messinger模型進行拓展，考慮冰晶黏附效應(yīng)，采用歐拉法計算了二維翼型混合相結(jié)冰。張麗芬等采用拉格朗日法計算粒子運動軌跡，考慮冰晶撞擊模型，計算了NACA0012翼型混合相結(jié)冰。總體來說，國內(nèi)混合相及冰晶的結(jié)冰研究剛剛起步。

2 混合相/冰晶結(jié)冰與過冷水結(jié)冰的差異及其危害

2.1 與過冷水結(jié)冰的差異

冰晶結(jié)冰和過冷大水滴結(jié)冰有很大不同，主要體現(xiàn)在結(jié)冰發(fā)生高度不同、結(jié)冰位置(環(huán)境)不同、結(jié)冰機理不同。

在結(jié)冰高度方面，根據(jù)Bravin等統(tǒng)計數(shù)據(jù)分析，冰晶結(jié)冰事故發(fā)生的平均溫度為-36 ℃，平均高度為10.7 km。事故發(fā)生高度比過冷水滴結(jié)冰發(fā)生高度(一般7 km以下)高。過冷水結(jié)冰通常發(fā)生在飛機起飛、爬升、進場、著陸等階段，而冰晶結(jié)冰主要發(fā)生在巡航階段。冰晶出現(xiàn)在強對流云層的上部，暖濕空氣急劇上升并凍結(jié)成微小冰晶，在風(fēng)的作用下冰晶向周圍擴散，如圖3所示。冰晶不易被肉眼觀測和機載雷達探測。飛機為躲避對流云層中易察覺的下部降雨區(qū)域，在上部冰晶區(qū)域巡航，使得發(fā)動機吸入大量冰晶，冰晶在發(fā)動機內(nèi)部熱環(huán)境下融化形成冰水混合物，進一步發(fā)生冰晶黏附產(chǎn)生積冰。

圖3 過冷大水滴和冰晶在對流云層中的分布[19]Fig.3 Distribution of supercooled large droplets and ice crystals in convective clouds[19]

在結(jié)冰位置(環(huán)境)方面，過冷水滴結(jié)冰一般發(fā)生在機翼、機身、進氣道短艙、帽罩等冷表面，而冰晶結(jié)冰主要發(fā)生在發(fā)動機低壓壓氣機級靜子葉片、傳感器等熱表面。冰晶在干冷表面撞擊后往往會彈開或破碎而不發(fā)生結(jié)冰。過冷水滴進入發(fā)動機內(nèi)部高溫環(huán)境會很快上升至冰點以上而不發(fā)生結(jié)冰，而冰晶進入后會部分融化形成冰水混合物，撞擊葉片，葉片表面形成液膜，隨后固態(tài)冰晶不斷黏附在葉片表面使葉片降溫累積成冰(見圖4)，冰晶不斷深入到壓氣機內(nèi)部，從而結(jié)冰區(qū)域也隨之延伸，如圖5所示。

圖4 撞擊表面?zhèn)鳠醾髻|(zhì)現(xiàn)象Fig.4 Heat and mass transfer on impingement surface

圖5 典型渦扇發(fā)動機及過冷水/冰晶結(jié)冰位置[10]Fig.5 Icing positions of supercooled droplets and ice crystals of typical turbofan engine[10]

結(jié)冰方式和結(jié)冰環(huán)境的不同必然會導(dǎo)致結(jié)冰機理不同，冰晶結(jié)冰機理的定性分析如前所述，主要影響因素之一是撞擊表面的液膜狀態(tài)，而液膜狀態(tài)與冰水含量、環(huán)境溫度、壓力、撞擊壁面溫度等因素有很大聯(lián)系，目前關(guān)于這些因素以及液膜狀態(tài)對冰晶結(jié)冰的影響機理研究還很缺乏?？偨Y(jié)冰晶結(jié)冰和過冷大水滴結(jié)冰的差異如表1所示。

表1 冰晶結(jié)冰和過冷大水滴結(jié)冰差異Table 1 Difference of icing between ice crystals and supercooled droplets

2.2 冰晶結(jié)冰危害

冰晶結(jié)冰在飛機各個迎風(fēng)部件的冷熱表面均有可能發(fā)生。在既含冰晶又含過冷水的混合氣象條件下，典型的冷表面如機翼、尾翼、風(fēng)擋等可能發(fā)生混合相結(jié)冰；典型的熱表面如發(fā)動機內(nèi)部核心處的靜子、各類傳感器熱表面等可能發(fā)生冰晶氣象條件下的結(jié)冰。不同部位的結(jié)冰會對飛機的飛行安全帶來不同程度的危害。

機翼、尾翼等升力表面的混合相結(jié)冰會降低升力，增大摩擦阻力和壓差阻力，從而影響到飛機的操控性和穩(wěn)定性；風(fēng)擋的結(jié)冰會影響飛行員的視線，降低飛行員對當(dāng)前飛行狀態(tài)以及外部環(huán)境的準確判斷；傳感器表面結(jié)冰會影響傳感器的正常工作，降低傳感器的精確度，獲得誤差較大的參數(shù)干擾飛行員的判斷；冰晶結(jié)冰通常伴隨著總溫傳感器異常(顯示0 ℃)。發(fā)動機核心部件結(jié)冰引起發(fā)動機推力下降，甚至引發(fā)喘振、熄火、冰脫落損傷葉片等嚴重危險。2016年1月一架使用GEnx發(fā)動機的波音787飛機發(fā)生高空冰晶結(jié)冰事件，發(fā)動機空中停車后無法重新啟動，幸運的是飛機依靠另一臺發(fā)動機著陸。

波音公司的統(tǒng)計數(shù)據(jù)表明自1990至今的航空發(fā)展歷程中至少發(fā)生了100起發(fā)動機不正常或失效事件，其中大多數(shù)與冰晶結(jié)冰有關(guān)。而Mason等經(jīng)過數(shù)據(jù)統(tǒng)計，從1990—2014年間有162起發(fā)動機冰晶結(jié)冰事故，涵蓋了12架不同的發(fā)動機。事故發(fā)生的地點主要集中在(亞)熱帶地區(qū)，文獻[10]顯示，中國東海岸是發(fā)動機冰晶結(jié)冰事故率相對較高的地區(qū)。冰晶結(jié)冰問題是飛機/發(fā)動機設(shè)計和適航驗證中面臨的關(guān)鍵問題,研究冰晶結(jié)冰問題對保障飛機的安全飛行十分重要。

3 混合相/冰晶結(jié)冰研究進展

3.1 冰風(fēng)洞地面試驗設(shè)備

冰風(fēng)洞(IWT)是地面試驗最先進的設(shè)備之一，能夠在一定范圍內(nèi)模擬結(jié)冰氣象條件，實現(xiàn)可控的結(jié)冰試驗。早在20世紀40～50年代，用于結(jié)冰試驗研究的初型冰風(fēng)洞就已經(jīng)建成，利用自然氣候制造冷環(huán)境。經(jīng)過七八十年的發(fā)展，世界很多國家如美國、加拿大、英國、意大利等都建立了冰風(fēng)洞，冰風(fēng)洞的類型越來越多，規(guī)模也越來越大，能夠?qū)崿F(xiàn)更加復(fù)雜的試驗條件。各國利用現(xiàn)有冰風(fēng)洞開展冰晶結(jié)冰的地面試驗技術(shù)研究，主要包括：在冰風(fēng)洞中實現(xiàn)冰晶生成、存儲和運輸?shù)龋倪M結(jié)冰試驗設(shè)備使其具備模擬冰晶結(jié)冰能力；冰晶粒子的撞擊特性研究；混合相條件下冰晶的黏附機理研究；冰晶/混合相結(jié)冰的影響因素研究等。

3.1.1 Cox冰風(fēng)洞

2002年在聯(lián)邦航空局(FAA)的資助下，威奇托州立大學(xué)、Cox公司和NASA Glenn研究中心合作，在Cox冰風(fēng)洞(見圖6)開展了混合相結(jié)冰氣象條件對熱防除冰系統(tǒng)性能影響的研究。

圖6 Cox IWT的整體布局[22]Fig.6 Overall layout of Cox IWT[22]

Cox研發(fā)了混合相和冰晶相氣象條件的生成方法。其冰晶粒子的產(chǎn)生方法有兩種，一是利用噴嘴噴出水霧然后凍結(jié)形成球形粒子，這樣生成的冰晶粒子接近球形；二是先在低溫下生成冰塊，然后利用旋轉(zhuǎn)機械切割出碎冰粒，這樣生成的粒子形狀較不規(guī)則且體積較大。同時噴射冰晶粒子和過冷水滴即可生成混合相條件。

2003年Cox的混合相結(jié)冰試驗是文獻記錄中有關(guān)混合相/冰晶條件下結(jié)冰展開的最早的研究，后續(xù)其他冰風(fēng)洞大多借鑒了Cox冰風(fēng)洞混合相的生成方式，并在此基礎(chǔ)上進行改進。

3.1.2 NRC高空試驗研究設(shè)備

加拿大國家研究委員會(National Research Council, NRC)燃氣輪機實驗室于2007年以來對發(fā)動機高空試驗研究設(shè)備(Research Altitude Test Facility, RATFac)進行冰晶結(jié)冰試驗設(shè)備的升級改造，增加了冰晶生成系統(tǒng)和運輸系統(tǒng)，見圖7。冰晶的生成利用旋轉(zhuǎn)電機切割冰塊的方法，調(diào)整不同刀片的角度和電機轉(zhuǎn)速可以控制冰晶粒子的大小，通過控制步進電機送入總冰塊的質(zhì)量來調(diào)節(jié)總冰晶的生成速率。后續(xù)利用此冰晶實驗裝置開展了很多小型樣件的冰晶結(jié)冰機理性研究，包括雙楔翼型(2012年、2016年)、半球形(2014年、2016年)、半圓柱形(2014年)、圓錐體(2014年)、NACA0012翼型(2018)。

圖7 NRC冰晶試驗系統(tǒng)[24]Fig.7 NRC ice crystals test system[24]

3.1.3 美國推進系統(tǒng)實驗室

美國Glenn研究中心于20世紀50年代初期建成了大型結(jié)冰研究風(fēng)洞(Icing Research Tunnel, IRT)，并成立了NASA的結(jié)冰研究中心。2015年以來NASA Glenn中心的推進系統(tǒng)實驗室(Propulsion Systems Laboratory, PSL)針對全尺寸高空發(fā)動機試驗裝置進行改造，利用噴霧系統(tǒng)噴出的水霧在冷氣流中快速凍結(jié)，形成冰晶。由于實驗系統(tǒng)(見圖8)龐大，升級改造復(fù)雜，目前開展了一些試驗參數(shù)如粒子直徑、液態(tài)水和冰晶云霧總溫的標(biāo)定工作。2018年開展了翼型標(biāo)模混合相結(jié)冰試驗，詳見3.2.3節(jié)。

圖8 NASA PSL冰晶結(jié)冰試驗系統(tǒng)[32]Fig.8 NASA PSL ice crystals accretion test system[32]

3.1.4 德國Braunschweig冰風(fēng)洞

德國Braunschweig冰風(fēng)洞是歐盟HAIC項目發(fā)展起來的,該冰風(fēng)洞采用了一套目前比較先進的試驗設(shè)備。Baumert等在此冰風(fēng)洞中針對NACA0012翼型和圓柱開展了混合相結(jié)冰試驗。整個系統(tǒng)由循環(huán)冰風(fēng)洞、冰晶制備和輸送系統(tǒng)(Ice crystal Generation and conveyance System, IGS)、制冷系統(tǒng)(Refrigeration System)3大部分組成生成，如圖9所示。

圖9 布倫瑞克冰風(fēng)洞試驗設(shè)備示意圖[38]Fig.9 Shematic of Braunschweig icing wind tunnel experiment equipment[38]

圖9中的冰晶制備和輸送系統(tǒng)(IGS)是整套設(shè)備的創(chuàng)新之處。首先由超聲霧化器產(chǎn)生一種具有均勻霧滴大小的準靜止霧，通過小型風(fēng)扇轉(zhuǎn)移到云室，在云室內(nèi)通過壓縮機的脈沖工作生成小冰晶，冰晶不斷產(chǎn)生并在云室的底部沉積，最后將收集到的冰晶放入-60 ℃的冰箱儲存。整套設(shè)備在冷卻室中運行。對于冰晶運輸專門設(shè)計了管道系統(tǒng)，采用一臺高壓徑向風(fēng)機進行空氣輸送。冰風(fēng)洞采用閉環(huán)設(shè)計，靜溫范圍在-20～30 ℃之內(nèi)。通過外部制冷系統(tǒng)控制溫度，盛有容量為4 000 L 冷卻介質(zhì)的容器提供冷凍液體之后進入換熱器，混合相生成系統(tǒng)連接的外部制冷系統(tǒng)可以提供80 kW的冷卻能量和36 kW的軸向風(fēng)扇輸入，風(fēng)扇可以使風(fēng)洞中氣流循環(huán)流動。冰風(fēng)洞噴口收縮比為10:1，軸向長度為3.5 m，相對較長的軸向延伸能夠使冰風(fēng)洞內(nèi)產(chǎn)生正確的氣流曲率，讓水滴和冰晶的流動特性得到增益。測試段為方形截面，側(cè)面長度為500 mm，流速可達到為40 m/s，在風(fēng)洞進口處安裝噴桿系統(tǒng)，使液滴能夠噴射進入風(fēng)道。液態(tài)水含量可在1.4～3.4 g/m范圍內(nèi)進行調(diào)節(jié)。噴管入口處的沉降室配備整流器和湍流濾網(wǎng)來提高流動質(zhì)量。

Cox的IWT、加拿大的RATFac、美國的PSL、德國布倫瑞克冰風(fēng)洞的相關(guān)參數(shù)總結(jié)如表2 所示。國內(nèi)有關(guān)混合相/冰晶結(jié)冰的試驗研究還在起步階段，建設(shè)具有模擬冰晶環(huán)境的大型冰風(fēng)洞是進行冰晶結(jié)冰試驗的關(guān)鍵。

表2 典型冰風(fēng)洞開展冰晶結(jié)冰研究的試驗?zāi)芰able 2 Experimental ability of typical icing wind tunnel on ice crystals accretion study

3.2 冰風(fēng)洞試驗現(xiàn)象

試驗發(fā)現(xiàn)，液態(tài)水的存在是冰晶結(jié)冰的主要影響因素。液態(tài)水使得冰晶可以黏附在暖表面如探測器或發(fā)動機靜子葉片上，黏附的冰晶在暖表面融化使得暖表面變冷而達到冰點，產(chǎn)生結(jié)冰，如圖4所示。為了定量分析實驗現(xiàn)象，理解液態(tài)水對冰晶結(jié)冰的影響，引入?yún)?shù)和c，定義為

(1)

(2)

式中：表示冰晶結(jié)冰氣象條件中所包含的所有液態(tài)水的比例；為冰晶融化比，表示冰晶粒子中融化的液態(tài)水所占的比例；LWC表示空氣中液態(tài)水滴的含量；LWC表示部分融化冰晶中液態(tài)水的含量；IWC表示空氣中的冰晶含量；TWC表示空氣中的總冰水含量。當(dāng)沒有LWC時，=，因此稱之為總?cè)诨剩瑸楸诨取?/p>

3.2.1 Cox冰風(fēng)洞試驗

Al-Khalil等在Cox冰風(fēng)洞中開展了低于0 ℃的混合相結(jié)冰試驗，其目的是為了驗證混合相環(huán)境對熱防冰系統(tǒng)性能的影響，同時研究冰晶粒子撞擊到機翼前緣后的物理現(xiàn)象。實驗對象是NACA0012；TWC相對較低，為0.3～1.4 g/m；較低，為0.165；溫度為-12～-7 ℃；為0.3～0.7。不同條件得到的冰型不同，見圖10。圖中:為速度;OAT為外界空氣溫度。在溫度相對較高的-5.6 ℃時，產(chǎn)生了不規(guī)則冰型，并且表面存在液態(tài)水；在溫度相對較低的-11.1 ℃ 時，產(chǎn)生的是光滑的干態(tài)冰。因此Al-Khalil等參考過冷水結(jié)冰將混合相結(jié)冰分為霜冰條件和明冰條件。

圖10 兩種典型結(jié)冰條件下的冰形[22]Fig.10 Ice shapes under two typical icing conditions[22]

實驗中觀察到了很明顯的侵蝕現(xiàn)象，特別是在明冰條件下，形成的積冰表面光滑且呈不透明的白色。

3.2.2 NRC冰風(fēng)洞試驗

2013年，NRC在試驗中摻混熱空氣，得到濕球溫度()高于0 ℃的混合氣流，使冰晶粒子運動中發(fā)生融化，融化的冰晶顆粒撞擊試驗件發(fā)生結(jié)冰，稱之為純冰晶融化相態(tài)結(jié)冰。試驗中不同的通過調(diào)整混合氣體的濕球溫度而獲得。NRC的研究表明，在10%～25%時，冰晶結(jié)冰比較明顯，保持相對穩(wěn)定的結(jié)冰強度水平，如圖11所示。出現(xiàn)這種現(xiàn)象的原因主要有：①比較小時，沒有足夠的液態(tài)水使得冰晶有效的黏附；②比較大時，液態(tài)水太多，使得表面結(jié)冰不夠致密以抵抗冰晶撞擊帶來的侵蝕，從而發(fā)生結(jié)冰質(zhì)量損失。另外，試驗發(fā)現(xiàn)增大粒子速度和直徑會減小“結(jié)冰平臺期”，這是因為冰晶侵蝕率的增大以及撞擊粒子動能耗散增大使得積冰融化，導(dǎo)致圖11中“結(jié)冰平臺期”在總?cè)诨矢偷臅r候就結(jié)束。需要注意的是，NRC的研究中，的增大是通過提高環(huán)境溫度來實現(xiàn)的，因此越高時，結(jié)冰表面的對流散熱會降低，這樣也會降低結(jié)冰率。以上說明，冰晶融化相態(tài)結(jié)冰存在結(jié)冰平臺期，對應(yīng)的冰晶融化率范圍會隨著冰晶參數(shù)的不同發(fā)生改變。

圖11 總?cè)诨屎徒Y(jié)冰強度的關(guān)系[26]Fig.11 Relationship between melting ratio and icing intensity[26]

NRC試驗發(fā)現(xiàn)，純冰晶融化相態(tài)結(jié)冰冰形通常為楔形冰。圖12是NRC在RATFac中的試驗結(jié)果，分別為0.25和0.4，總水含量為6 g/m。試驗發(fā)現(xiàn)，經(jīng)過一定時間后，結(jié)冰達到了穩(wěn)態(tài)，冰形不隨時間發(fā)生變化，說明冰晶結(jié)冰量和冰晶侵蝕損失量達到了平衡，稱之為穩(wěn)態(tài)結(jié)冰。對比圖12(a)和圖12(b)，=0.25的情況下，=0.14時的結(jié)冰量明顯高于=0.314的結(jié)冰量。原因如前所述，=0.14時達到了結(jié)冰平臺期，結(jié)冰速率相對較大，而=0.314時，結(jié)冰速率反而有所降低。對比圖12(a)和圖12(c)，均為0.14,=0.4時的結(jié)冰量明顯低于=0.25時的結(jié)冰量，說明速度對結(jié)冰量具有較強的反作用，速度越大，冰晶侵蝕導(dǎo)致的結(jié)冰損失量會顯著增大。

對比圖12(a)和圖13，表明在=0.25,=3 ℃,=0.14～0.2時，隨著TWC的增大，結(jié)冰量明顯增大。然而，當(dāng)速度增大(=0.4)時，NRC試驗發(fā)現(xiàn)，此時TWC從6 g/m增大到9.6 g/m時，結(jié)冰量并沒有明顯的增大，冰尖處的高度僅從大約6.3 mm增大到約6.9 mm。說明速度越大，TWC對結(jié)冰量的影響減弱,這和速度越大時結(jié)冰侵蝕量越大有關(guān)。

圖12 NRC試驗: 半圓柱形表面結(jié)冰冰形 (TWC=6 g/m3, 1 ℃≤Twb≤5.5 ℃)[26]Fig.12 NRC experiments: icing shape for cylinder at (TWC=6 g/m3, 1 ℃≤Twb≤5.5 ℃)[26]

圖13 NRC試驗: 半圓柱形表面結(jié)冰冰形 (Ma=0.25, TWC=12 g/m3)Fig.13 NRC experiments: icing shape for cylinder at Ma=0.25, TWC=12 g/m3

以上說明冰晶融化相態(tài)結(jié)冰時，速度越大，冰晶侵蝕效果越強，將大大削弱由于其他因素引起的結(jié)冰增量。

濕球溫度對冰層的增長有很大影響，如前所述，NRC的RATFac設(shè)備中云霧的液相是依靠粒子融化而獲得，而空氣的濕球溫度是影響粒子融化的一個重要特征參數(shù)。另外，濕球溫度較高時，冰層經(jīng)常發(fā)生脫落，這是因為溫度越高，多孔冰層(Porous Ice Layer)中殘留的液態(tài)水增多，使得冰層和表面的黏附力降低，易被氣流吹拂。在發(fā)動機內(nèi)部，氣流的濕球溫度是相對較高的，結(jié)冰易脫落而損傷發(fā)動機葉片，甚至導(dǎo)致發(fā)動機熄火。

上面討論的現(xiàn)象都是濕球溫度高于0 ℃的情況(比如在發(fā)動機內(nèi)部)，冰晶在運動過程中融化從而產(chǎn)生液態(tài)水，冰水混合粒子撞擊表面發(fā)生結(jié)冰。在濕球溫度低于0 ℃的混合相條件下，由于液態(tài)水滴的存在，飛機機翼表面也可能發(fā)生結(jié)冰；或者混合相中的冰晶撞擊到加熱表面上時發(fā)生融化黏附，從而降低表面溫度繼而產(chǎn)生結(jié)冰。在濕球溫度低于0 ℃的混合相條件下，圖11中的結(jié)冰包線的右側(cè)將不會下降到0，當(dāng)=1時，此時的結(jié)冰就是過冷水滴結(jié)冰過程。

2016年Currie和Fuleki在NRC的RATFac中開展了溫度低于0 ℃的混合相態(tài)結(jié)冰試驗。試驗發(fā)現(xiàn)此條件下的結(jié)冰與表面之間表現(xiàn)出更好的黏附性，同時，試驗在半球形表面得到了帶蘑菇狀底座的錐狀冰，如圖14(a)所示。而在大于0 ℃的純冰晶融化相態(tài)結(jié)冰時，試驗得到的冰形是和表面相切的錐狀冰，如圖14(b)所示。分析其原因，可能還是和冰層密度以及冰晶粒子的侵蝕效果有關(guān)。如前所述，濕球溫度越高，冰層中殘留液態(tài)水越多，因此冰層密度越低，導(dǎo)致冰角越容易被侵蝕而不易形成蘑菇狀基底。

圖14 濕球溫度對冰形的影響[28]Fig.14 Influence of wet bulb temperature on ice shape[28]

2018年Struk等在NRC的RATFac中針對NACA0012開展了混合相結(jié)冰試驗。目的是進一步驗證冰晶結(jié)冰機理，試驗條件同球形試件相似。當(dāng)濕球溫度大于0 ℃時，在試驗中觀察到了結(jié)冰和脫落過程。當(dāng)濕球溫度進一步降低到0 ℃后，積冰不再脫落，積冰厚度會隨著時間增長，冰形為楔形冰，并未出現(xiàn)半球形表面的蘑菇底座冰形。試驗中還觀察到了積冰冰形出現(xiàn)了不均勻的現(xiàn)象，原因是流向的擾動和混合相條件生成的不均勻性。

3.2.3 NASA冰風(fēng)洞試驗

2018年NASA的PSL實驗室開展了NACA0012翼型混合相態(tài)結(jié)冰試驗和分析，試驗狀態(tài)類似渦扇噴氣式發(fā)動機的壓氣機處，通過控制主流氣體的相對濕度來控制噴流水滴在運動過程中的蒸發(fā)冷卻以及凍結(jié)程度，從而達到一定凍結(jié)率(融化率+凍結(jié)率=1)的混合相條件。

試驗獲得了產(chǎn)生冰晶結(jié)冰的最小融化率，也觀測到了結(jié)冰“平臺期”現(xiàn)象，平臺期的最大黏附效率為0.2，比NRC試驗得出的黏附效率要低。但是NASA試驗沒有獲得“平臺期”的最大冰晶融化率，沒有觀測到穩(wěn)態(tài)結(jié)冰現(xiàn)象，如圖15所示，在融化率為0.91時，結(jié)冰量仍較多。原因可能是NASA的冰晶條件是通過水滴部分凍結(jié)產(chǎn)生的，試驗中還存在并未凍結(jié)的過冷液滴，試驗中的濕球溫度低于且接近0 ℃。濕球溫度低于0 ℃時沒有結(jié)冰“平臺期”和NRC的濕球溫度低于0 ℃ 的試驗結(jié)果一致。但是NASA和NRC試驗中的冰形有較大區(qū)別。

圖15 不同冰晶融化率不同時刻的冰形[36]Fig.15 Ice shapes at different melting ratios of ice crystal[36]

Bartkus等基于不同的表面能量平衡提出了兩種不同的積冰類型：融化主導(dǎo)式結(jié)冰和凍結(jié)主導(dǎo)式結(jié)冰。前者發(fā)生在濕球溫度高于0 ℃；后者發(fā)生在濕球溫度低于0 ℃。融化主導(dǎo)式結(jié)冰為固態(tài)冰晶直接在表面積聚，此時的積冰對于表面的黏附力較弱，結(jié)冰量主要取決于固態(tài)冰晶中未融化的固態(tài)相冰晶。試驗觀察到在這種情況下，冰晶撞擊過程中只有部分冰殘留在表面，大部分冰晶通過飛濺、反彈和侵蝕作用損失。凍結(jié)主導(dǎo)式積冰表現(xiàn)在冰晶撞擊表面的同時液態(tài)水部分凍結(jié)，此時的積冰對于表面有很強的黏附力，積冰量主要取決于液態(tài)水的凍結(jié)量。

3.2.4 Braunschwerg冰風(fēng)洞試驗

Baumert等在Braunschwerg冰風(fēng)洞對NACA0012 (弦長0.5 m) 翼型和圓柱開展了混合相結(jié)冰試驗，通過調(diào)整過冷水滴含量和冰晶含量來獲得不同的，稱之為混合相態(tài)結(jié)冰試驗。速度40 m/s； TWC為7～17 g/m；=0～0.6；溫度為0 ℃、-5 ℃、-15 ℃，濕度達到飽和，=，即試驗中濕球溫度小于或等于0 ℃。

為了定量評估駐點處黏附/凍結(jié)成冰的量，定義駐點積冰系數(shù)，表示前緣位置積冰量占總水收集量的比值：

(3)

Braunschwerg冰風(fēng)洞的試驗發(fā)現(xiàn)，冰晶粒子只有在駐點附近撞擊破碎彈開后會二次撞擊表面，在其他區(qū)域撞擊破碎后會被氣流帶走，見圖16，這說明二次粒子的撞擊對整個結(jié)冰的貢獻不大。

圖16 純冰晶條件下冰晶撞擊表面現(xiàn)象[39]Fig.16 Phenomenon of ice crystal impinging on surface under glaciated ice condition[39]

試驗發(fā)現(xiàn)在環(huán)境溫度低于0 ℃的混合相結(jié)冰時，溫度越低，表面結(jié)冰量越少，如圖17所示。溫度低對結(jié)冰會產(chǎn)生兩個相反的影響：① 溫度低使冰凍結(jié)更易發(fā)生；② 但會導(dǎo)致表面滯留和溢流水減小，進而導(dǎo)致黏附系數(shù)降低，固態(tài)冰晶黏附量減少。另外，可能在溫度比較低(如-10 ℃以下)時：溫度越低，表面結(jié)冰越干燥脆弱，容易被侵蝕，使結(jié)冰量減小。因此混合相結(jié)冰數(shù)值模擬中的黏附模型和侵蝕模型應(yīng)考慮溫度帶來的這些影響。

圖17 NACA0012混合相條件120 s時結(jié)冰情況 (TWC=12 g/m3, mr=0.28) [39]Fig.17 NACA0012 icing under mixed phase condition after 120 s at TWC=12 g/m3, mr=0.28[39]

試驗發(fā)現(xiàn)溫度越高時，結(jié)冰更容易脫落，NACA0012表面結(jié)冰為楔形冰，和NRC的結(jié)果一致。同時，在以及粒子直徑均相同的情況下，駐點處結(jié)冰厚度隨TWC呈線性關(guān)系，也就是駐點積冰系數(shù)基本不隨TWC變化，如圖18所示。但布倫瑞克冰風(fēng)洞試驗中的速度較低，為40 m/s，以上結(jié)論不一定適用于速度高的情況，有待進一步考證。

圖18 前緣結(jié)冰120 s情況與TWC的關(guān)系 (mr=0.2 s)[39]Fig.18 Leading edge ice accretion depending on total water content after 120 s at mr=0.2[39]

比較圓柱和NACA0012翼型的結(jié)冰，試驗發(fā)現(xiàn)在相同的結(jié)冰條件下，駐點處最終的結(jié)冰高度增長速率和冰角均基本相同，只不過圓柱表面達到穩(wěn)定的結(jié)冰增長速率和穩(wěn)定冰角所需的時間更長。原因主要是圓柱的直徑(60 mm)大于NACA0012的前緣直徑(15.8 mm)，這說明越薄的物體在混合相條件下結(jié)冰，越容易達到穩(wěn)定增長。

需要指出的是，冰晶融化相態(tài)結(jié)冰以及混合相態(tài)結(jié)冰的試驗數(shù)據(jù)仍然不是特別充足，上述NRC、Cox以及布倫瑞克冰風(fēng)洞試驗中的速度均較低，部分試驗結(jié)論可能在高速情況下并不適用，不同速度、溫度、融化率、TWC情況下的黏附效應(yīng)和侵蝕效應(yīng)對結(jié)冰的影響還有待開展更多的基礎(chǔ)性試驗。

3.3 混合相/冰晶結(jié)冰數(shù)值模擬

冰晶結(jié)冰的數(shù)值模擬主要集中在如下3個方面：冰晶運動的軌跡及融化相變；冰晶撞擊壁面時的黏附模型、侵蝕模型；結(jié)冰熱力學(xué)模型。Wright、Habashi、 Villedieu、 Trontin以及冰晶/混合相的結(jié)冰熱力學(xué)模型均基于Messinger模型進行擴展。由于在溫暖環(huán)境下面結(jié)冰含水量較大，結(jié)冰附著力不夠經(jīng)常導(dǎo)致冰脫落，Bennani和Kintea等將冰脫落因素考慮到了數(shù)值模擬中。

3.3.1 冰晶運動融化相變

冰晶運動軌跡及運動過程中的相變影響到冰晶的表面收集量以及到達部件表面時的狀態(tài)(融化率、溫度、粒徑等)。目前冰晶運動軌跡的模擬基本上都是沿用過冷水滴軌跡計算常用的拉格朗日法和歐拉法，但冰晶粒子在環(huán)境氣流中運動的控制方程與傳統(tǒng)上過冷水的控制方程主要有兩點區(qū)別，一是有融化相變現(xiàn)象，在壓氣機前部有一段總溫介于 0～10 ℃的區(qū)域，冰晶在溫暖環(huán)境下運動時，伴隨著與環(huán)境之間的熱量交換，有可能發(fā)生融化相變現(xiàn)象；二是冰晶粒子的形狀不規(guī)則，非球形冰晶顆粒在融化過程中隨著液態(tài)水的增多，顆粒形狀漸趨球形。因此在運動相變過程中顆粒形狀變化引起的對阻力系數(shù)、運動軌跡以及對流換熱系數(shù)的影響需考慮入內(nèi)。

Villedieu和Trontin等在現(xiàn)有的水滴結(jié)冰模型的基礎(chǔ)上，考慮了冰晶結(jié)冰相關(guān)的換熱相變模型、阻力模型和撞擊模型，提出了新的冰晶結(jié)冰計算模型。關(guān)于冰晶粒子運動過程中的融化相變模型，Villedieu做出了如下簡化：

1) 固態(tài)階段，冰晶粒子溫度低于融化溫度。顆粒為純固態(tài)冰晶，冰晶顆粒與外界環(huán)境進行對流換熱，冰晶溫度上升，冰晶升華，直至溫度達到融化溫度。

2) 融化階段，冰晶粒子溫度等于融化溫度。在進一步換熱作用下，冰晶由外側(cè)開始融化，水膜環(huán)繞內(nèi)側(cè)冰晶核，在此階段部分融化冰晶顆粒吸收的熱量全部用來供給融化潛熱和蒸發(fā)潛熱，溫度維持不變。

3) 液態(tài)階段，冰晶粒子溫度大于融化溫度。冰晶顆粒完全融化呈液滴狀態(tài)，繼續(xù)從外環(huán)境吸收熱量，液態(tài)水滴蒸發(fā)，溫度繼續(xù)上升。

Villedieu根據(jù)3個階段不同的熱力學(xué)特性分別建立控制方程，定義粒子球形度Ф表示粒子形狀的不規(guī)則程度，并基于粒子球形度改進了阻力系數(shù)模型和換熱模型。

(4)

式中：表示相同體積球體表面積；表示不規(guī)則粒子的真實表面積。

考慮到部分融化冰晶表面的蒸發(fā)效應(yīng)，其蒸發(fā)模型為

(5)

Norde對上述Villedieu相變模型進行了計算驗證，將得到的融化時間與Hauk做的懸浮冰晶粒子融化試驗結(jié)果進行對比。通過對比不同的阻力系數(shù)模型和換熱模型計算融化時間與試驗值的誤差，找到了與試驗結(jié)果吻合最好的阻力系數(shù)模型和換熱模型。

阻力模型采用了Ganser模型：

(6)

式中：為相對雷諾數(shù)；為斯托克斯形狀因子；為牛頓形狀因子。

換熱系數(shù)模型采用Villedieu等的模型：

(7)

式中：為普朗特數(shù)。

Ayan、Grift、Aouizerate等采用拉格朗日方法分別對翼型、發(fā)動機壓氣機葉片、三維發(fā)動機進氣道以及壓氣機級內(nèi)冰晶顆粒的運動和相變以及可能的結(jié)冰進行了模擬；Iuliano等采用歐拉方法分析了運動過程中不規(guī)則粒子形狀對運動軌跡的影響，對比不同的阻力模型的計算結(jié)果與試驗值之間誤差，同時考慮了粒子形狀對換熱性能的影響。

3.3.2 黏附模型

冰晶黏附特性是其結(jié)冰與過冷水結(jié)冰較大不同之處。2014年，Villedieu等給出了冰晶黏附系數(shù)計算式，見式(8)，認為黏附系數(shù)與表面水膜厚度、冰晶粒子(部分融化)表面水膜厚度、粒子直徑、法向動能恢復(fù)系數(shù)有關(guān)，并且能夠保證當(dāng)液膜厚度越大、法向動能恢復(fù)系數(shù)越小時，冰晶粒子黏附系數(shù)越大。

(8)

Villedieu根據(jù)NASA-NRC的2個混合相結(jié)冰工況的試驗數(shù)據(jù)，通過對比計算結(jié)果，得到系數(shù)=40時結(jié)果的誤差最小。

Villedieu的黏附概率模型雖然考慮了液膜厚度、融化冰晶液態(tài)水厚度、冰晶粒徑、法向動能(速度)的影響，但模型中的僅通過兩個試驗數(shù)據(jù)得到，經(jīng)過更多試驗數(shù)據(jù)的對比驗證，該模型存在一些不足之處：① 液膜厚度很難精確預(yù)測；② 從前緣到撞擊極限范圍，粒子法向恢復(fù)系數(shù)快速增長，使得除前緣駐點處的黏附系數(shù)預(yù)測偏低，導(dǎo)致除駐點處的結(jié)冰量預(yù)測偏低，而駐點處結(jié)冰量則看起來大很多，出現(xiàn)了和實際物理現(xiàn)象不符的結(jié)冰冰形。

為了避免上述問題，2016年Trontin等提出了新模型，這個模型使用了冰晶粒子的平均融化比、總液態(tài)水含量和總含水量。利用NRC冰晶運動融化后撞擊表面結(jié)冰的試驗數(shù)據(jù)以及NASA-NRC混合氣象條件結(jié)冰的試驗數(shù)據(jù)來校核。

=MAX(,)

(9)

式中：

(10)

=(+)

(11)

(12)

(13)

考慮到表面溫度越低時，冰晶黏附越少，特別是溫度較低的干態(tài)冰表面時，冰晶不會黏附。于是，黏附模型中的常數(shù)和表面溫度有關(guān)：

(14)

的取值根據(jù)NRC冰晶融化相態(tài)結(jié)冰的試驗數(shù)據(jù)得到。根據(jù)加拿大NRC的試驗結(jié)果，相同LWC/TWC(所有液態(tài)水包括過冷水滴和融化冰晶所含水比上總水含量)時，部分融化的純冰晶條件下的黏附系數(shù)要高于混合相條件下的黏附系數(shù)。

2016年Trontin的黏附模型相對來說更容易用已有試驗數(shù)據(jù)來校驗，不需要在計算過程中和液膜厚度進行耦合迭代計算。通過考慮LWC/TWC對黏附效率的貢獻，間接考慮了壁膜厚度的影響?？偟膩碚f，現(xiàn)有的黏附模型沒有考慮粒子速度、尺寸、形狀的影響，這意味著它僅僅是一個基于經(jīng)驗的“宏觀模型”，其滿足有效性的范圍有限，有待發(fā)展更加普適的冰晶粒子黏附模型。

3.3.3 侵蝕模型

在結(jié)冰表面形成一定的冰層后，冰晶對冰層的撞擊就會驅(qū)逐冰層表面的一部分冰，表現(xiàn)為冰晶的侵蝕作用。2003年Al-Khalil等的試驗清楚地證明了侵蝕對結(jié)冰過程的影響。侵蝕會導(dǎo)致冰層形狀的變化，結(jié)冰量的減少以及冰層表面粗糙度的變化等。因此要得到準確的結(jié)冰預(yù)測就需要考慮冰晶的侵蝕作用。2010年Wright等率先根據(jù)實驗數(shù)據(jù)建立了自己的有關(guān)侵蝕率的經(jīng)驗公式，認為侵蝕不僅與冰層的物理性質(zhì)有關(guān)，還與撞擊粒子的速度、角度以及慣性等有關(guān),后面的研究都是對其公式的改進。Trontin等基于與Currie試驗結(jié)果一致的假設(shè)，建立了侵蝕模型，認為：

1) 侵蝕程度受切向速度的影響甚于受法向速度的影響。

2) 液態(tài)水含量越高，冰晶侵蝕效果越明顯；反之，固態(tài)冰晶含量越高，侵蝕效果越弱。

3) 侵蝕效果隨冰晶直徑的增加而增強。

4) 冰層溫度越低，冰越堅硬，侵蝕效果越弱。

5) 冰晶溫度越低，侵蝕效果越強。

6) 侵蝕效果隨表面冰層曲率的增加而增強。

侵蝕效率表示侵蝕效應(yīng)對結(jié)冰量影響的嚴重程度，受到多個因素的影響，表達式為

(15)

式中：為侵蝕強度；MMD為冰晶的平均質(zhì)量直徑；為撞擊到表面冰晶的切向速度；為冰晶活化能，取值48.2 kJ/mol；是通用氣體常數(shù)，其值為8.314 J·mol·K；為壁面溫度；為融化溫度；為冰層的當(dāng)?shù)厍?下角標(biāo)0代表參考值。

結(jié)合Currie等的試驗數(shù)據(jù)，Trontin和Villedieu進一步得到如下侵蝕效率經(jīng)驗?zāi)Ｐ停?/p>

(16)

式中：指滯留在冰層中的液態(tài)水占冰層的質(zhì)量分數(shù)；=0.6；=0.015；=84.5 m/s。式(16) 不僅結(jié)合Currie有限的冰風(fēng)洞試驗數(shù)據(jù)得到了經(jīng)驗系數(shù)，并且將滯留在冰層中的液態(tài)水質(zhì)量分數(shù)的影響考慮進來,認為存在臨界的液態(tài)水質(zhì)量分數(shù)即0.6，此時侵蝕效率將達到無窮大，表面結(jié)冰為含水量很高的濕雪冰，冰層很容易被侵蝕造成質(zhì)量損失。這和冰晶融化率較高時結(jié)冰強度反而降低的試驗現(xiàn)象(圖11)相一致。另外，在駐點附近，由于切向速度很小，使得模型的侵蝕效果出現(xiàn)比實際小很多的情況，導(dǎo)致計算的冰形在駐點附近出現(xiàn)異常，因此Trontin引入了曲率的影響(駐點附近曲率大)來修正侵蝕模型。侵蝕的質(zhì)量流量由侵蝕效率決定，但肯定要小于表面的結(jié)冰量。

目前的侵蝕模型是基于有限冰風(fēng)洞試驗得到的經(jīng)驗公式。首先，試驗樣本少，導(dǎo)致此經(jīng)驗公式適用范圍有限；其次，冰風(fēng)洞中的部分參數(shù)無法精確標(biāo)定和測量，將引入誤差；再次，在高液態(tài)水含量時沒有考慮冰晶顆粒被液態(tài)水夾帶著流走的效果，使得侵蝕效果偏低。因此有必要開展專門的冰晶粒子侵蝕試驗研究，獲得大量試驗數(shù)據(jù)后總結(jié)適用面更廣的侵蝕模型。

3.3.4 結(jié)冰過程模擬

目前混合相/冰晶結(jié)冰的熱力學(xué)模型均是在Messinger模型的基礎(chǔ)上進行改進發(fā)展而來的。其中Wright、Trontin、Norde等均在過冷水結(jié)冰熱力學(xué)模型中增加冰晶相關(guān)的質(zhì)量和能量項，得到了冰晶/混合相結(jié)冰熱力學(xué)模型。Trontin、Norde等在結(jié)冰模型中不僅考慮了冰晶的黏附效應(yīng)、侵蝕效應(yīng)，還將溢流水在多孔冰層中的滯留現(xiàn)象考慮入內(nèi)。

混合相/冰晶結(jié)冰過程中控制體的質(zhì)量和能量守恒示意圖如圖19所示。

圖19 控制體內(nèi)的質(zhì)量守恒與能量守恒示意圖[56]Fig.19 Mass and energy conservation in control volume[56]

對于每個控制體，進入控制體和離開控制體的液態(tài)水質(zhì)量流量守恒：

(17)

對于每個控制體，進入控制體和離開控制體的固態(tài)冰質(zhì)量流量守恒：

(18)

對于每個控制體，進入控制體和離開控制體的能量守恒：

(19)

考慮冰晶黏附效應(yīng)后，撞擊冰晶中實際黏附在結(jié)冰表面參與結(jié)冰過程的冰晶質(zhì)量流量為

(20)

式中：為黏附系數(shù)；為入口邊界處冰晶速度；IWC為固態(tài)冰含量；為冰晶融化比；為壁面處冰晶的收集系數(shù)；為控制體的底面積。

(21)

因此考慮侵蝕效應(yīng)后的結(jié)冰質(zhì)量流量和溢流水質(zhì)量流量為

(22)

(23)

在冰晶結(jié)冰條件下，有試驗表明液態(tài)水會被滯留在結(jié)冰顆粒之間。Trontin等稱此時的結(jié)冰為濕雪冰(Slushy Ice)，濕雪冰的密度、強度將明顯低于過冷水結(jié)冰冰層的，從而增大侵蝕效果。

Trontin等在結(jié)冰熱力學(xué)模型中考慮了液態(tài)水在多孔冰層中的滯留。方法如下：IWC/TWC比值越高，液態(tài)水滯留在積冰顆粒之間的可能性越大。也就是說，混合相中的冰粒子越多，液態(tài)水越有可能被困在粒子之間的孔隙中，在冰層中積累；反之，液態(tài)水占壁面撞擊粒子的比例越大，液態(tài)水也不容易滯留，而是容易溢流出去。

另外，考慮冰晶結(jié)冰存在冰晶黏附、侵蝕效應(yīng)以及溢流水滯留等特殊現(xiàn)象，在時間離散化方面，采用多步算法來捕捉整個結(jié)冰過程中由于上述現(xiàn)象引起的結(jié)冰形狀的演變，這種演變甚至有可能使得結(jié)冰達到一種動態(tài)平衡，即結(jié)冰量將不再發(fā)生變化。也就是說，將結(jié)冰時間Δ分成個時間步，每個時間步長內(nèi)按照如下3個步驟開展結(jié)冰計算：

1) 計算結(jié)冰質(zhì)量和能量守恒(考慮黏附效應(yīng))。

2) 考慮侵蝕效應(yīng)更新溢流水和結(jié)冰量(質(zhì)量損失，溫度不變)。

3) 考慮液態(tài)水的滯留，得到液態(tài)水的滯留量和溢流量。

由于在步驟2)中溫度保持不變，這意味著存在如下假設(shè)：被侵蝕的冰粒子和飛濺的水滴在離開表面之前已經(jīng)與冰層達到熱平衡。

到目前為止，所有的混合相/冰晶結(jié)冰熱力學(xué)模型認為冰晶只能在撞擊處黏附或反彈離開表面，并不會被溢流水帶走至下游。然而在比較稀的濕雪冰時，很有可能存在流動的濕雪冰，即表面的結(jié)冰會和液態(tài)水一起流動至下游。

冰晶結(jié)冰過程中，當(dāng)濕球溫度高于0 ℃時，結(jié)冰后期經(jīng)常會發(fā)生脫落，這個現(xiàn)象在目前的結(jié)冰模擬中還未考慮。因此要想將混合相/冰晶結(jié)冰過程模擬更加準確，模型中尚有較多因素的影響有待考慮。

4 研究展望

綜上所述，混合相/冰晶結(jié)冰過程和過冷水結(jié)冰機理存在很大差別，有關(guān)混合相/冰晶結(jié)冰的試驗和數(shù)值模擬研究工作才剛剛起步。各種因素包括冰晶融化率、溫度、速度、撞擊角度等對混合相/冰晶結(jié)冰過程中的黏附效應(yīng)和侵蝕效應(yīng)有很大影響，現(xiàn)有的黏附模型和侵蝕模型存在很多不足之處，現(xiàn)有的結(jié)冰模型尚有溢流冰和冰脫落現(xiàn)象未考慮。因此針對混合相/冰晶結(jié)冰，重點要關(guān)注如下研究內(nèi)容：

1) 混合相/冰晶結(jié)冰的試驗?zāi)M和測量技術(shù)，包括冰晶粒子的生成、混合相的形成、冰水混合含量的測量、冰晶融化率的測量等。

2) 結(jié)冰風(fēng)洞中混合相/冰晶結(jié)冰機理試驗，以獲得豐富的試驗數(shù)據(jù)。

3) 冰晶粒子沿程運動融化相變試驗和數(shù)值模擬，以獲得能夠準確評估冰晶融化時間和融化率的模型和模擬方法。

4) 冰晶黏附試驗和黏附模型，定量評估表面液膜對冰晶黏附系數(shù)的影響。

5) 冰晶侵蝕試驗和侵蝕模型，定量評估溫度、冰晶速度、撞擊角度、冰層密度等因素對侵蝕系數(shù)的影響。

6) 濕雪冰的模擬以及剪切力作用下冰水混合物在表面的運動，從而修正現(xiàn)有的結(jié)冰熱力學(xué)模型。

7) 在濕球溫度高于0 ℃時冰晶結(jié)冰容易發(fā)生脫落，有必要開展冰脫落現(xiàn)象及新型防除冰方法的研究。

5 結(jié) 論

混合相/冰晶結(jié)冰是飛機和發(fā)動機面臨的和過冷水結(jié)冰有很大區(qū)別的結(jié)冰現(xiàn)象，深入掌握其結(jié)冰機理，理解結(jié)冰危害，對提高飛機/發(fā)動機的飛行安全具有重要意義。

針對混合相/冰晶結(jié)冰問題，從與過冷水結(jié)冰的區(qū)別、結(jié)冰危害,冰晶結(jié)冰的主要地面試驗設(shè)備,地面混合相/冰晶試驗現(xiàn)象,冰晶結(jié)冰數(shù)值模擬4個方面進行了詳細的論述。指出了目前有關(guān)混合相/冰晶結(jié)冰研究的不足之處和重點需要關(guān)注的方面。