劉 禎, 孔維梁, 劉 洪, 王福新

(上海交通大學(xué) 航空航天學(xué)院, 上海 200240)

0 引 言

冰附著力受到關(guān)注的主要原因在于：一是相對(duì)傳統(tǒng)的防除冰方法，現(xiàn)有研究已表明防冰材料可以減小冰附著力、推遲液滴結(jié)冰時(shí)間以及減少冰積聚，從而以其獨(dú)特的優(yōu)勢(shì)在防除冰中顯示出巨大的潛能[1-3]，而防除冰方法的開(kāi)發(fā)也需要對(duì)冰附著力的認(rèn)知；二是對(duì)飛機(jī)上非防護(hù)表面的結(jié)冰脫落的預(yù)測(cè)需要準(zhǔn)確地分析其融冰和脫冰特性，因其事關(guān)飛機(jī)的起降安全。

冰附著力的研究主要包含冰附著力的機(jī)理，冰附著力的影響因素2個(gè)方面。

對(duì)冰附著力的機(jī)理的研究一般都是從基底表面的濕潤(rùn)性能出發(fā)的，以Young等式[4]、Wenzel模型和Cassie-Baxter模型[5-6]為理論依據(jù)。冰與基底表面之間的相互作用會(huì)影響冰附著力的大小[2-3,7]。Soltis等人[3]通過(guò)實(shí)驗(yàn)得出，隨著表面粗糙度增加，過(guò)冷水與表面微觀(guān)結(jié)構(gòu)接觸越充分，因而冰附著力越大。而Saito 等人[7]發(fā)現(xiàn)當(dāng)基底的表面粗糙度增大到一定程度時(shí)，其冰附著力反而會(huì)減小，這是因?yàn)楸砻嫖⒂^(guān)結(jié)構(gòu)與冰之間的空氣泡的數(shù)量的增加導(dǎo)致冰與基底之間附著點(diǎn)減少。表面粗糙度相同時(shí)，基底表面微觀(guān)結(jié)構(gòu)不同也會(huì)導(dǎo)致冰附著力的不同，因?yàn)樘厥獾谋砻嫖⒂^(guān)結(jié)構(gòu)會(huì)使冰有不同數(shù)量的附著點(diǎn)[3]。綜上所述，研究冰附著力的機(jī)理必須考慮冰與基底表面之間的相互作用，冰與基底之間附著點(diǎn)會(huì)影響冰附著力的大小。

冰附著力的影響因素包括基底因素和環(huán)境因素[1-3]。其中影響冰附著力的環(huán)境因素主要包括：成冰(動(dòng)態(tài)或者靜態(tài))以及溫度。常見(jiàn)的成冰方式有2種：一種為靜態(tài)結(jié)冰，即靜止的水滴或者水層緩慢結(jié)冰；另一種為動(dòng)態(tài)結(jié)冰，即水滴以動(dòng)態(tài)的方式撞擊基底表面然后結(jié)冰[1-3]。靜態(tài)結(jié)冰中的水基于毛細(xì)作用力達(dá)到力學(xué)平衡，而動(dòng)態(tài)結(jié)冰涉及水滴撞擊基底的過(guò)程[1-3]。在飛機(jī)飛行中的飛機(jī)結(jié)冰問(wèn)題主要為動(dòng)態(tài)結(jié)冰問(wèn)題。在動(dòng)態(tài)結(jié)冰的條件下，影響冰附著力的環(huán)境因素包括：溫度(T)、風(fēng)速(v)、水含量(LWC)和平均粒徑(MVD)[1-3]。

文獻(xiàn)中因?yàn)榛缀铜h(huán)境因素以及測(cè)量方法的不同，實(shí)驗(yàn)結(jié)果也不同[2]，范圍從70～931kPa，這是因?yàn)楸街Φ挠绊懸蛩乇姸喽覐?fù)雜，不同的基底、不同的測(cè)試方法、不同的成冰方式以及不同的溫度均會(huì)導(dǎo)致冰附著力結(jié)果的不同。

有關(guān)溫度對(duì)冰附著力的影響的研究起始于1962年，Stallabrass和Price[8]的實(shí)驗(yàn)結(jié)果呈現(xiàn)出動(dòng)態(tài)結(jié)冰的冰附著力隨溫度降低而降低的趨勢(shì)。而在1983年，Itagaki等人[9]在分析其實(shí)驗(yàn)結(jié)果中卻得出動(dòng)態(tài)結(jié)冰的冰附著力受溫度影響很小。溫度是影響冰附著力的重要環(huán)境因素，而溫度對(duì)冰附著力的影響呈現(xiàn)出多樣的規(guī)律。這是因?yàn)楸街Φ挠绊懸蛩氐膹?fù)雜性：在最近的研究中，在Tarquini 等人[2]的實(shí)驗(yàn)中溫度范圍為-15～0℃，結(jié)果表明動(dòng)態(tài)結(jié)冰的冰附著力隨著溫度降低而增大；Soltis等人[3]的實(shí)驗(yàn)中溫度范圍為-16～-8℃，在其所選材料和溫度范圍內(nèi)，動(dòng)態(tài)結(jié)冰的冰附著力隨溫度線(xiàn)性下降；Guerin等人[10]在研究冰附著力時(shí)發(fā)現(xiàn)了動(dòng)態(tài)結(jié)冰的冰附著力隨溫度變化的3個(gè)區(qū)間(溫度范圍-25～0℃)，區(qū)間的范圍和基底材料有關(guān)，冰附著力在溫度從0℃下降中會(huì)依次經(jīng)過(guò)AZ(Adhesive Zone)、TZ(Transition Zone)、CZ(Cohesive Zone)，冰附著力先增加后減小。以上的實(shí)驗(yàn)研究不但溫度范圍不同，其他冰附著力的影響因素也不同。此外，Janjua[11]探索了靜態(tài)結(jié)冰的冰附著力受基底溫度以及環(huán)境溫度的影響，得出的結(jié)論是兩者均會(huì)影響冰附著力，但是環(huán)境溫度的影響更大，并從傳熱的角度解釋了這一現(xiàn)象，這表明在實(shí)驗(yàn)中需要考慮傳熱因素。

除此之外，其它環(huán)境因素對(duì)冰附著力的影響也在展開(kāi)，Tarquini[2]等人的實(shí)驗(yàn)中未發(fā)現(xiàn)冰附著力與LWC的關(guān)系，但是發(fā)現(xiàn)在高LWC情況下，冰附著力有所下降。Momen等人[12]在實(shí)驗(yàn)中發(fā)現(xiàn)在低LWC和小MVD情況下，冰附著力會(huì)略微減小,而且冰附著力隨著LWC的增加而減小，并指出MVD對(duì)冰附著力的影響需要進(jìn)一步的研究。但關(guān)于環(huán)境因素對(duì)冰附著力影響的研究，正如Tarquini等人[2]所指出的，大多數(shù)文獻(xiàn)未列出詳細(xì)具體的環(huán)境條件。例如：Stallabrass和Price[8]只提供了水霧的MVD；Itagaki[9]未提供任何水霧的參數(shù)；Tarquini等人[2]、Soltis等人[3]、Momen等人[12]的實(shí)驗(yàn)中提供了LWC以及MVD，但是只采用了常規(guī)小粒徑水霧(50μm以下)，未研究大粒徑水霧的情況。這種情況的存在使得各研究間缺乏可比性。其次，現(xiàn)可查閱的文獻(xiàn)中對(duì)冰附著力的研究只分別進(jìn)行了靜態(tài)結(jié)冰和動(dòng)態(tài)結(jié)冰的實(shí)驗(yàn)且有的進(jìn)行了相關(guān)機(jī)理的研究，但未進(jìn)行相應(yīng)的對(duì)比分析。

另外，值得注意的是，張辰等人[13]研究發(fā)現(xiàn)，過(guò)冷大粒徑水滴撞擊壁面與小粒徑水滴撞擊壁相比存在異常,表明大粒徑水霧的結(jié)冰特性和常規(guī)小粒徑水霧相比是存在較大差異的。

基于以上分析，本文在保持LWC不變的條件下，開(kāi)展了3種不同粒徑的水霧(包含大粒徑水霧和常規(guī)小粒徑水霧)分別在靜態(tài)結(jié)冰和動(dòng)態(tài)結(jié)冰下的冰附著力(剪切力)的對(duì)比性試驗(yàn)，并分析了兩者的差異。

1 實(shí)驗(yàn)方案

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

如圖1所示，實(shí)驗(yàn)是在溫控箱(高低溫交變濕熱試驗(yàn)箱LRHS-504B-LJS)中進(jìn)行的，實(shí)驗(yàn)系統(tǒng)主要分為4個(gè)系統(tǒng)：噴霧發(fā)生系統(tǒng)、溫度控制系統(tǒng)、實(shí)驗(yàn)參數(shù)測(cè)定系統(tǒng)和冰附著力測(cè)試系統(tǒng)。噴霧發(fā)生系統(tǒng)的主體是霧化噴嘴，霧化噴嘴兩頭通路分別為水路和氣路，水路(實(shí)驗(yàn)中采用的水是0℃的純凈水冰水混合物)上有水泵、過(guò)濾器和壓力傳感器以及流量傳感器，氣路上有氣泵、過(guò)濾器、壓力傳感器。噴霧系統(tǒng)通過(guò)調(diào)節(jié)水壓、流量和氣壓來(lái)調(diào)節(jié)水霧的LWC與MVD，其中MVD的范圍為40～280μm。

1 水霧噴頭; 2 氣路; 3 水路; 4 USB插口; 5 離心冰附著力測(cè)試系統(tǒng)

溫度控制系統(tǒng)主要通過(guò)溫控箱，可以在-50～0℃范圍內(nèi)控制溫度。同時(shí)，在溫控箱外的管路上使用絕熱材料，以減少熱交換。

1.2 參數(shù)的設(shè)置與測(cè)量

實(shí)驗(yàn)中環(huán)境溫度T的測(cè)量基于溫控箱中的溫度傳感器。

根據(jù)水含量的定義，水霧LWC可通過(guò)測(cè)量單位時(shí)間內(nèi)撞擊單位面積基底表面的水滴質(zhì)量來(lái)表征。單位時(shí)間內(nèi)撞擊單位面積基底表面的水滴質(zhì)量在實(shí)驗(yàn)中控制為不變量，其值不變，可認(rèn)為L(zhǎng)WC不變。水霧的MVD通過(guò)Winner319工業(yè)噴霧激光粒度分析儀來(lái)測(cè)定，如圖2所示。

圖2 Winner319工業(yè)噴霧激光粒度分析儀

1.3 冰附著力的測(cè)量

離心裝置是一種常用的測(cè)量冰附著力的方法[1-3,8-12,14]，在冰附著力的測(cè)量中具有簡(jiǎn)單、快速和重復(fù)性較高的特點(diǎn)。而且，因?yàn)楸奶厥庑再|(zhì)，動(dòng)態(tài)結(jié)冰所形成的冰難以直接施加力的作用，離心裝置適合用于動(dòng)態(tài)結(jié)冰的冰附著力的測(cè)量。

如圖3所示，冰附著力測(cè)試系統(tǒng)采用離心測(cè)力裝置，其中電機(jī)采用Maxon RE35，該電機(jī)可以在低溫環(huán)境下工作。懸臂梁臂長(zhǎng)620mm，基底位于懸臂梁的一端，材料為Al6061，平面尺寸為20mm×30mm。在懸臂梁的另一端設(shè)有配重，以保持整個(gè)離心測(cè)力裝置的平衡。冰的質(zhì)量m的測(cè)量在冰脫落發(fā)生后，采用電子秤測(cè)量。靜態(tài)結(jié)冰的冰，將0℃的水置入冰附著力測(cè)量裝置上的無(wú)底容器(底部為基底)生成；動(dòng)態(tài)結(jié)冰的冰，采用噴霧生成。

在實(shí)驗(yàn)中，在冰的附著完成后，起動(dòng)電機(jī)，控制懸臂梁以一定轉(zhuǎn)速旋轉(zhuǎn)，逐漸增加轉(zhuǎn)速，直到冰脫落發(fā)生。冰脫落發(fā)生時(shí)，冰附著力與離心力相等，即F=mrω2，其中m為冰的質(zhì)量，r為半懸臂長(zhǎng)，ω為轉(zhuǎn)速。相應(yīng)的剪切應(yīng)力τ=F/A，其中A為基底的面積。

圖3 離心測(cè)力裝置

2 實(shí)驗(yàn)結(jié)果與分析

2.1 實(shí)驗(yàn)結(jié)果

在實(shí)驗(yàn)中采用Al6061作為基底材料，保持不變。在不同溫度T=-5、 -10、-15、-20和-25℃時(shí)，分別測(cè)量不同形成方式冰的附著力，第一實(shí)驗(yàn)組采用靜態(tài)結(jié)冰的方式，第2、3、4實(shí)驗(yàn)組采用動(dòng)態(tài)結(jié)冰的方式，調(diào)節(jié)噴霧系統(tǒng)保持單位時(shí)間內(nèi)撞擊單位面積基底表面的水滴質(zhì)量為2g/min不變，MVD分別為40、80和250μm，測(cè)量冰附著力。為消除傳熱因素的影響，在實(shí)驗(yàn)中冰的附著完成后，靜置30min后再進(jìn)行冰附著力的測(cè)量。每種條件下，冰附著力測(cè)量值為測(cè)量3次后所取平均值。實(shí)驗(yàn)得到了如圖4所示的結(jié)果。

2.2 動(dòng)態(tài)結(jié)冰與靜態(tài)結(jié)冰的對(duì)比

在本實(shí)驗(yàn)中，冰附著力隨實(shí)驗(yàn)條件不同呈現(xiàn)出不同的規(guī)律。由圖4的實(shí)驗(yàn)結(jié)果可以看出，對(duì)于靜態(tài)結(jié)冰和MVD分別為80和250μm的動(dòng)態(tài)結(jié)冰，隨著溫度的降低，冰附著力逐漸增大，與Stallabrass、Price[8]和Tarquini等人[2]以及 Soltis等人[3]得出的結(jié)果相同。但是，對(duì)于MVD為40μm的情況，在溫度降低為-15℃時(shí)，動(dòng)態(tài)結(jié)冰的冰附著力反而減??；以及對(duì)于MVD為80μm的情況，在溫度降低為-20℃時(shí)，動(dòng)態(tài)結(jié)冰的冰附著力減小。Guerin等人[10]的實(shí)驗(yàn)結(jié)果所呈現(xiàn)的曲線(xiàn)同樣出現(xiàn)了轉(zhuǎn)折點(diǎn)，這里將該轉(zhuǎn)折點(diǎn)所對(duì)應(yīng)的溫度稱(chēng)為臨界溫度。不同的是，本實(shí)驗(yàn)測(cè)得動(dòng)態(tài)結(jié)冰的冰附著力隨溫度的變化在臨界溫度之后是無(wú)規(guī)律的，而非隨著溫度的降低而增加。

圖4 Al6061基底在T=-5、-10、-15、-20和-25℃時(shí)的冰附著力

不同粒徑下動(dòng)態(tài)結(jié)冰和靜態(tài)結(jié)冰的冰附著力的對(duì)比結(jié)果如圖5所示，可以發(fā)現(xiàn)動(dòng)態(tài)結(jié)冰和靜態(tài)結(jié)冰這2種成冰方式在到達(dá)臨界溫度前，對(duì)冰附著力的大小影響很小，動(dòng)態(tài)結(jié)冰與靜態(tài)結(jié)冰所成冰的冰附著力相近。但是在溫度到達(dá)臨界溫度之后，動(dòng)態(tài)結(jié)冰的冰附著力與靜態(tài)結(jié)冰相比明顯要小，而且其隨溫度的變化是無(wú)規(guī)律的。

為了確定臨界溫度的范圍，在臨界溫度附近增加了實(shí)驗(yàn)，對(duì)于MVD為40μm的情況，該臨界溫度在-15℃左右(-17.5～-12.5℃)；對(duì)于MVD為80μm的情況該臨界溫度在-17.5℃左右(-20～-15℃)；對(duì)于MVD為250μm的情況，在-25℃之前未發(fā)現(xiàn)臨界溫度。

在實(shí)驗(yàn)中同時(shí)觀(guān)察到了在臨界溫度附近冰的類(lèi)型的變化(見(jiàn)圖6)。冰的類(lèi)型有明冰、霜冰以及混合冰3種：明冰形成的環(huán)境條件是溫度相對(duì)較高、水含量較大，部分水滴在撞擊基底表面后會(huì)在結(jié)冰前流動(dòng)，形成的冰致密透明，密度較大；霜冰形成的環(huán)境條件是溫度相對(duì)較低、水含量較低、水滴平均粒徑較小，水滴在撞擊基地表面時(shí)迅速結(jié)冰，甚至在撞擊前結(jié)冰，形成的冰疏松多孔，乳白色，密度較小?；旌媳敲鞅退幕旌稀?/p>

動(dòng)態(tài)結(jié)冰的機(jī)理是水滴通過(guò)撞擊與基底表面接觸，在凝固前進(jìn)入基底表面微觀(guān)結(jié)構(gòu)，然后凍結(jié)[2-3,11]。在到達(dá)臨界溫度前，圖6(a)、(b)和(c)中動(dòng)態(tài)結(jié)冰所成冰為明冰，形成的冰透明致密，在形成最終冰的形態(tài)之前會(huì)有水膜的出現(xiàn)，與靜態(tài)結(jié)冰所成冰的類(lèi)型是相同的，所以此時(shí)靜態(tài)結(jié)冰和動(dòng)態(tài)結(jié)冰的冰附著力相近。而在到達(dá)臨界溫度后，圖6(d)和(e)中動(dòng)態(tài)結(jié)冰所成冰為混合冰或者霜冰，為乳白色的冰，在冰形成的過(guò)程中，可以觀(guān)察到水霧在撞擊基底表面后迅速結(jié)冰，而且圖6(d)和(e)中動(dòng)態(tài)結(jié)冰所成冰可以明顯觀(guān)察到疏松多孔的結(jié)構(gòu)，此時(shí)表面微觀(guān)結(jié)構(gòu)與冰之間空氣泡的數(shù)量較多，冰與基底之間的附著點(diǎn)較少，因此冰附著力較小。所以對(duì)于MVD為40μm，溫度為-15℃以及MVD為80μm，溫度為-20℃的情況，動(dòng)態(tài)結(jié)冰的冰附著力反而減小。

(a) Ice adhesion of freezer and impact(MVD=40μm)

(b) Ice adhesion of freezer and impact(MVD=80μm)

(c) Ice adhesion of freezer and impact(MVD=250μm)

(a) -5℃，MVD=250μm (b) -15℃，MVD=250μm (c) -25℃，MVD=250μm

綜上所述，與靜態(tài)結(jié)冰相比，在動(dòng)態(tài)結(jié)冰過(guò)程中，由于存在水滴對(duì)基底表面的撞擊過(guò)程，在不同的環(huán)境條件下，所成冰的類(lèi)型可能完全不同。在基底表面可能出現(xiàn)明冰、霜冰或者混合冰的冰層，嚴(yán)重影響冰的附著力的大小。這也可以解釋?zhuān)煌瑮l件下冰的附著力隨溫度的變化在到達(dá)臨界溫度后呈現(xiàn)出不同的趨勢(shì)，且不具有明顯的規(guī)律性。平均粒徑越大，越不容易產(chǎn)生混合冰或者霜冰，所以在實(shí)驗(yàn)結(jié)果中，其它環(huán)境參數(shù)不變的情況下，平均粒徑越大，臨界溫度越低。

2.3 平均粒徑對(duì)冰附著力的影響

從圖5可以初步分析得出平均粒徑對(duì)冰附著力的影響，在實(shí)驗(yàn)范圍內(nèi)，其他環(huán)境參數(shù)不變的情況下，平均粒徑越大，臨界溫度越低。臨界溫度的一個(gè)重要特點(diǎn)是冰的類(lèi)型發(fā)生了明顯的變化，在到達(dá)臨界溫度之前，所成冰為明冰，而在到達(dá)臨界溫度附近時(shí)所成冰為混合冰，到達(dá)臨界溫度之后所成冰則為混合冰或者霜冰。如圖6所示，在基底和LWC不變的情況下，大粒徑水霧形成的冰結(jié)構(gòu)較為致密，為明冰，而小粒徑水霧形成的冰結(jié)構(gòu)較為疏松，為混合冰或者霜冰，冰的類(lèi)型不同導(dǎo)致了冰附著力的差異。

大粒徑的環(huán)境更利于產(chǎn)生明冰，而小粒徑的環(huán)境更利于產(chǎn)生混合冰和霜冰。而且溫度越低，越利于產(chǎn)生混合冰和霜冰。因此，在其他環(huán)境參數(shù)不變的情況下，平均粒徑越大，臨界溫度越低。

3 結(jié) 論

本實(shí)驗(yàn)研究了以Al6061為基底在-25～-5℃之間的靜態(tài)結(jié)冰與動(dòng)態(tài)結(jié)冰(水霧LWC恒定，MVD=40、80、250μm)的冰附著力的變化情況，得出了以下結(jié)論：

(1) 存在一個(gè)臨界溫度，動(dòng)態(tài)結(jié)冰和靜態(tài)結(jié)冰這2種成冰方式在到達(dá)臨界溫度前，對(duì)冰附著力的大小影響較小，而且冰附著力的大小隨著溫度降低而增加；但是在溫度到達(dá)臨界溫度之后時(shí)，動(dòng)態(tài)結(jié)冰所成冰的冰附著力相對(duì)于靜態(tài)結(jié)冰明顯要小，且冰的類(lèi)型從明冰逐漸轉(zhuǎn)變?yōu)榛旌媳蛩?/p>

(2) 平均粒徑對(duì)冰附著力的影響主要體現(xiàn)在對(duì)臨界溫度的影響。在實(shí)驗(yàn)范圍內(nèi)，其它環(huán)境參數(shù)不變的情況下，平均粒徑越大，臨界溫度越低。

為更深入地探討在環(huán)境參數(shù)影響下的變化規(guī)律，在后續(xù)工作中，有必要擴(kuò)大可變化的參數(shù)類(lèi)型和變化范圍等，并著重研究大粒徑水霧動(dòng)態(tài)結(jié)冰的冰附著力的變化特性。

參考文獻(xiàn):

[1]Yeong Y H, Athanasios M, Eric L. Atmospheric ice adhesion on water-repellent coatings: wetting and surface topology effects[J]. Langmuir, 2015, 31(48): 13107-13116.

[2]Stefania T, Carlo A, Alidad A, et al. Investigation of ice shedding properties of superhydrophobic coatings on helicopter blades[J]. . Cold Regions Science and Technology, 2014, 100: 50-58.

[3]Jared S, Jose P, Timothy E, et al. Ice adhesion mechanisms of erosion-resistant coatings[J]. AIAA Journal, 2015, 53(3): 654-662.

[4]Young T. An essay on the cohesion of fluids[J]. Phil Trans R Soc London, 1805, 95: 65-87.

[5]Wenzel R N. Resistance of solid surfaces to wetting by water[J]. Ind Eng Chem, 1936, 28(8): 988-994.

[6]Cassie A B D, Baxter S. Wettability of porous surfaces[J]. Trans Faraday Soc, 1944, 40: 546-551.

[7]Saito H, Takai K, Yamauchi G. Water- and ice-repellent coatings[J]. Surf Coat Int, 1997, 80(4): 168-171.

[8]Stallabrass J R, Price R D. On the adhesion of ice to various materials[J]. National Research Laboratory, Report LR-350, 1962.

[9]Itagaki K. Mechanical ice release processes I. Self-shedding from high speed rotors[R]. CCREL Report 83-26, 1983.

[10]Frederic G, Caroline L, Marie-Isabelle F, et al. Analytical model based on experimental data of centrifuge ice adhesion tests with different substrates[J]. Cold Regions Science and Technology, 2016, 121: 93-99.

[11]Zaid A J. The influence of freezing and ambient temperature on the adhesion strength of ice[J]. Cold Regions Science and Technology, 2017, 140: 14-19.

[12]Momen G, Jafari R, Farzaneh M. Ice repellency behaviour of superhydrophobic surfaces: effects of atmospheric icing conditions and surface roughness[J]. Appl Surf Sci, 2015, 349: 211-218.

[13]Zhang C, Liu H. Effect of drop size on the impact thermodynamics for supercooled large droplet in aircraft freezing[J]. Physics of Fluids, 2016: 062107.

[14]Laforte C, Laforte J L. How a solid coating can reduce ice adhesion on structures[C]//Proceedings of the 10th International Workshop of Atmospheric Icing of Structures, Brno, Czech Republic, 2002.