來流速度對(duì)防冰表面溢流水流動(dòng)換熱的影響

鄭梅，董威，朱劍鋆，郭之強(qiáng)

(上海交通大學(xué)機(jī)械與動(dòng)力工程學(xué)院，上海200240)

來流速度對(duì)防冰表面溢流水流動(dòng)換熱的影響

鄭梅，董威*，朱劍鋆，郭之強(qiáng)

(上海交通大學(xué)機(jī)械與動(dòng)力工程學(xué)院，上海200240)

為研究來流速度對(duì)防冰表面溢流水流動(dòng)形態(tài)及換熱的影響，基于空氣-水兩層相互作用的質(zhì)量、動(dòng)量和能量守恒，建立防冰表面溢流水水膜流動(dòng)換熱及破裂的數(shù)學(xué)模型，分析了防冰表面溢流水在不同來流條件下的流動(dòng)形態(tài)和表面換熱情況。計(jì)算分析表明:來流速度增加時(shí)，防冰表面相同位置處的連續(xù)水膜厚度減小，水膜破裂位置隨之延后;較高來流速度條件下，破裂處水膜厚度稍有增加，使得破裂后形成的溪流厚度和寬度增大;作為主要的表面散熱項(xiàng)，連續(xù)水膜表面蒸發(fā)及對(duì)流換熱熱流均隨來流速度的增加而增大。此外，由水膜破裂引起的表面溢流水流態(tài)變化對(duì)防冰表面蒸發(fā)熱流有一定影響。

防冰;溢流水;流動(dòng)換熱;水膜;來流速度

0 引言

結(jié)冰是影響飛行安全的一大隱患［1］。通常認(rèn)為，結(jié)冰發(fā)生的原因是云層中含有溫度低于冰點(diǎn)的亞穩(wěn)態(tài)液態(tài)水，當(dāng)飛機(jī)穿越由這些過冷水滴組成的云團(tuán)時(shí)，便會(huì)在迎風(fēng)部件表面發(fā)生結(jié)冰現(xiàn)象［2］。研究表明，飛機(jī)結(jié)冰會(huì)破壞飛機(jī)的氣動(dòng)外形，影響穩(wěn)定性和操縱性;嚴(yán)重的結(jié)冰現(xiàn)象極有可能會(huì)造成機(jī)毀人亡的重大飛行事故［3-4］。因此，有效的防/除冰系統(tǒng)對(duì)于結(jié)冰環(huán)境下的飛行安全具有十分重要的意義。

防冰系統(tǒng)能保證部件表面溫度維持在冰點(diǎn)以上。當(dāng)云層中的過冷水滴撞擊表面時(shí)，若防冰系統(tǒng)提供的熱量不足以完全蒸發(fā)表面撞擊水，未蒸發(fā)的液態(tài)水會(huì)在氣動(dòng)力作用下沿部件表面向后流動(dòng)，形成溢流水。溢流水的存在及其在表面呈現(xiàn)的不同流動(dòng)形態(tài)會(huì)直接影響部件表面防冰熱流的分布情況，尤其在防冰不足的情況下，溢流水會(huì)對(duì)表面結(jié)冰區(qū)域和積冰冰形產(chǎn)生決定性影響。Messinger模型［5］作為結(jié)冰與防冰表面流動(dòng)換熱計(jì)算分析的基礎(chǔ)模型，考慮部件表面的質(zhì)量守恒和能量守恒，但并未引入溢流水流動(dòng)的因素。此后的研究工作也大多圍繞結(jié)冰預(yù)測及防冰分析而展開［6-13］。針對(duì)液膜流動(dòng)及破裂問題，研究工作主要集中于降膜流動(dòng)及破裂問題［14-17］，考慮飛機(jī)防冰表面的水膜破裂及溪流形態(tài)的文獻(xiàn)還十分有限。Al-Khalil等人［18］曾對(duì)防冰表面的水膜流動(dòng)進(jìn)行數(shù)學(xué)建模，耦合了部件表面的能量平衡及溢流水的質(zhì)量和能量平衡，并提出采用最小能量法預(yù)測水膜破裂。Silva等人［19-21］同樣將熱力學(xué)第一定律應(yīng)用于翼型表面及溢流水的分析、建模，其重點(diǎn)在于對(duì)表面溫度分布的預(yù)測;針對(duì)水膜破裂問題，沿用了Al-Khalil的最小能量法，但提出了濕潤系數(shù)的不同定義?？傮w而言，已有模型的計(jì)算結(jié)果與實(shí)驗(yàn)數(shù)據(jù)相比還是存在一定的誤差，同時(shí)在預(yù)測液膜破裂內(nèi)在機(jī)理等方面仍存在一定的局限性。

本文主要關(guān)注溢流水在防冰表面的流動(dòng)換熱以及水膜破裂等問題。基于空氣-水兩層相互作用的質(zhì)量、動(dòng)量和能量守恒，建立防冰表面溢流水水膜流動(dòng)換熱及破裂的數(shù)學(xué)模型，計(jì)算分析了來流速度對(duì)表面溢流水分布、流態(tài)及表面換熱的影響，獲得了防冰表面連續(xù)水膜厚度、破裂處溪流厚度和寬度，以及連續(xù)水膜表面主要散熱熱流的分布情況，并分析了溢流水及其流態(tài)對(duì)表面換熱的影響。

1 物理描述

云層中的過冷水滴撞擊飛機(jī)防冰部件表面，受表面加熱作用蒸發(fā)，部分未蒸發(fā)的撞擊水在氣動(dòng)力作用下形成溢流水，沿部件表面向后流動(dòng)。溢流水在流動(dòng)過程中會(huì)形成連續(xù)水膜，但通常情況下，這種流動(dòng)是不穩(wěn)定的。

在撞擊區(qū)域內(nèi)，由于水滴的收集量相對(duì)較大，部件表面通常完全濕潤并形成連續(xù)流動(dòng)的水膜。但在撞擊區(qū)域之外，隨著水膜不斷向后流動(dòng)，在表面張力、氣動(dòng)剪切力、壁面粘滯力和壓差力等力的作用下以及蒸發(fā)、對(duì)流換熱等傳熱傳質(zhì)過程的共同作用下，水膜在流動(dòng)過程中會(huì)發(fā)生動(dòng)量和能量損耗。在此過程中，水膜厚度會(huì)不斷發(fā)生變化，而其流動(dòng)狀態(tài)也會(huì)變得不穩(wěn)定，易發(fā)生破裂形成溪流。此時(shí)，部件表面將分為干、濕兩個(gè)區(qū)域，如圖1所示。

圖1 防冰部件表面溢流水流動(dòng)形態(tài)示意圖Fig.1 Runback water flow on anti-icing surface

1.1 連續(xù)水膜流動(dòng)換熱

如圖2所示，部件表面形成連續(xù)水膜流動(dòng)，流動(dòng)過程中存在質(zhì)量、動(dòng)量和能量的交換。假設(shè)部件表面水膜流動(dòng)及空氣邊界層均為不可壓縮、定常、層流流動(dòng)。由于水膜很薄，且空氣中液態(tài)水含量(LWC)很小，因此僅考慮空氣對(duì)水膜的單向耦合作用。相對(duì)空氣剪切力的大小，可忽略水膜重力影響，而僅考慮空氣剪切力及表面壓力差對(duì)水膜流動(dòng)的驅(qū)動(dòng)效應(yīng)。當(dāng)防冰系統(tǒng)運(yùn)行足夠長時(shí)間后，部件表面溫度始終保持在0℃以上，因此本文僅考慮無結(jié)冰情況下的穩(wěn)態(tài)水膜流動(dòng)。

圖2 防冰部件表面水膜流動(dòng)示意圖Fig.2 W ater film and air boundary flow on anti-icing surface

圖2中，mevap為水膜表面蒸發(fā)質(zhì)量，mimp為水滴撞擊質(zhì)量，min和mout分別為流入和流出控制體積的質(zhì)量，Qevap為蒸發(fā)散熱，Qimp為水滴撞擊動(dòng)能，Qin和Qout分別為控制體積進(jìn)出口焓值，Qc為表面對(duì)流換熱，Qd為加熱撞擊水滴所損失的熱能，Qw為壁面導(dǎo)熱。

1.2 水膜破裂

連續(xù)水膜在流動(dòng)過程中易發(fā)生破裂，在表面形成若干條溪流，如圖3所示為假定的破裂處的溪流模型在y-z平面上的橫截面形狀。

圖3 溪流模型示意圖Fig.3 Rivulet model

圖3中，δ0為溪流厚度，W為溪流寬度，L表示單條溪流及其臨界干表面的總寬度。假設(shè)溪流截面為圓的一部分，其半徑為R，因此在氣-液界面上溪流的幾何外形可以表示為:

式中，θ為液相與氣相之間形成的角度，即接觸角。如圖3中所示，氣-液-固三相接觸點(diǎn)存在力平衡條件，即Young方程［22］:

其中，σsg、σsl、σlg分別表示固-氣界面、固-液界面、液-氣界面張力。

根據(jù)幾何關(guān)系，溪流的厚度和寬度可以表示為:

2 數(shù)學(xué)模型

2.1 連續(xù)水膜流動(dòng)換熱

取控制體積內(nèi)的水膜為研究對(duì)象，應(yīng)同時(shí)滿足質(zhì)量守恒、動(dòng)量定理及熱力學(xué)第一定律，得到水膜流動(dòng)的連續(xù)性方程、動(dòng)量方程及能量方程:

其中，δ1和δ2分別表示水膜厚度和空氣邊界層厚度。p為機(jī)翼表面靜壓力;τ0為壁面粘滯應(yīng)力，τair為氣-液界面空氣對(duì)水膜的剪切應(yīng)力，剪切應(yīng)力由τ= μ(d u/d y)確定，μ為動(dòng)力粘度;uδ1為氣-液界面水膜速度;u∞，s為水滴撞擊表面時(shí)沿s方向的速度;ρwater、λwater、Cp，water均為水的物性參數(shù)，分別表示水的密度、導(dǎo)熱系數(shù)及定壓比熱容;hc為對(duì)流換熱系數(shù);Tw-s為氣-液界面水膜表面溫度;T∞為自由來流溫度;Td為碰撞水滴的溫度，可以視為與來流溫度相等;Le為液態(tài)水的蒸發(fā)潛熱;uwater和Twater分別為水膜邊界層的速度分布和溫度分布。根據(jù)邊界條件可以表示為:

其中，Twall為防冰部件表面溫度;uδ1、δ1、τair、Tw-s均為關(guān)于s的函數(shù)。

考慮空氣對(duì)水膜的單向耦合作用，在連續(xù)水膜流動(dòng)換熱模型中須引入空氣邊界層的動(dòng)量方程。取空氣邊界層的控制體積為研究對(duì)象。與水膜控制體積的分析方法相類似，其動(dòng)量方程為:

式(10)中，τ'air為氣-液界面水膜對(duì)空氣的粘滯力，與τair互為作用力與反作用力;ue為空氣邊界層靠近主流位置的速度;ρa(bǔ)ir為自由來流空氣密度;uair為空氣邊界層的速度分布，同樣可以根據(jù)邊界條件求得:

2.2 水膜破裂

在水膜破裂處，水膜與溪流之間存在質(zhì)量守恒和能量守恒關(guān)系:

其中，下標(biāo)f和r分別指代水膜和溪流。

根據(jù)上述兩個(gè)守恒方程，可以化簡得到破裂點(diǎn)處無量綱控制方程:

其中，

上式中δf為破裂處臨界水膜厚度，δ+為對(duì)應(yīng)的無量綱值，R+為無量綱溪流截面半徑。

在給定破裂位置的條件下，根據(jù)式(14)可求得破裂位置溪流截面半徑的無量綱參數(shù)，進(jìn)而可以獲得溪流的幾何外形參數(shù)。

對(duì)于連續(xù)水膜流動(dòng)換熱及水膜破裂模型的詳細(xì)推導(dǎo)過程可參考文獻(xiàn)［23］。對(duì)于水膜破裂位置的確定，可以按照文獻(xiàn)［18］所提出的最小能量法得出破裂處的臨界液膜厚度。對(duì)于本文的算例而言，撞擊區(qū)域內(nèi)形成的連續(xù)水膜厚度均無法達(dá)到臨界液膜厚度，可以認(rèn)為在撞擊極限處水膜將發(fā)生破裂。目前看來最小能量法得到的破裂點(diǎn)與實(shí)驗(yàn)觀察到的水膜破裂位置還存在很大誤差，還需要進(jìn)一步開展研究工作。

2.3 防冰表面熱流

水膜破裂后，防冰部件表面被分為干表面、濕表面，其換熱特性與連續(xù)水膜流動(dòng)時(shí)的完全濕潤表面以及無溢流水時(shí)的干表面有較大區(qū)別。若不考慮溢流水的存在，即撞擊極限之外，完全為干表面，此時(shí)表面無蒸發(fā)熱流且對(duì)流換熱僅考慮部件表面與空氣之間的熱傳遞過程。若不考慮水膜破裂，即部件表面始終由連續(xù)水膜覆蓋，完全為濕表面，此時(shí)僅考慮水膜表面的蒸發(fā)熱流以及水膜與空氣之間的對(duì)流換熱。當(dāng)考慮水膜破裂時(shí)，翼型表面被分為干/濕兩個(gè)部分，此時(shí)僅濕表面存在蒸發(fā)熱流，而對(duì)流換熱熱流必須綜合考慮濕表面上溪流和空氣之間的對(duì)流換熱以及干表面上壁面與空氣之間的對(duì)流換熱。

因此，三種不同溢流水分布形態(tài)條件下，表面蒸發(fā)熱流及對(duì)流換熱熱流可以分別根據(jù)以下公式進(jìn)行計(jì)算:

(1)完全濕潤表面(即連續(xù)水膜表面):

(2)干/濕表面(即發(fā)生水膜破裂的表面):

(3)干表面(即無溢流水表面):

其中，F(xiàn)為水膜破裂處表面濕潤系數(shù)，可表示為

3 模型求解

求解表面水膜流動(dòng)換熱及水膜破裂模型的前提是部件表面空氣流場及水滴撞擊特性的計(jì)算結(jié)果。因此，先通過CFD軟件及水滴撞擊計(jì)算程序獲得部件表面空氣流場和表面水滴撞擊特性。在Fluent中計(jì)算空氣流場采用k-ε湍流模型;進(jìn)出口邊界分別采用速度入口和壓力出口;翼型壁面采用無滑移邊界條件?？諝饬鲌龅諗亢螅捎脷W拉法計(jì)算部件表面局部水收集系數(shù)。外流場和局部水收集系數(shù)的計(jì)算結(jié)果將作為本文模型求解必要的邊界條件。

水膜流動(dòng)換熱及破裂模型的控制方程均在Matlab中編程求解。首先迭代計(jì)算連續(xù)水膜控制方程;連續(xù)水膜模型迭代收斂后，可獲得水膜表面各項(xiàng)熱流值;其次，在給定破裂位置可對(duì)水膜破裂模型進(jìn)行求解，并獲得溪流的厚度和寬度，以及表面換熱。整個(gè)計(jì)算過程如圖4所示。

圖4 計(jì)算過程Fig.4 Computational procedure

4 計(jì)算結(jié)果及分析

Zhang K等人［24］利用DIP技術(shù)對(duì)NACA 0012翼型表面液態(tài)水的流動(dòng)形態(tài)和水膜/溪流厚度進(jìn)行了實(shí)驗(yàn)測量。實(shí)驗(yàn)中NACA 0012翼型弦長為0.101m，攻角0°。采用本文所述的數(shù)學(xué)模型針對(duì)該實(shí)驗(yàn)?zāi)Ｐ瓦M(jìn)行了連續(xù)水膜流動(dòng)換熱及水膜破裂的計(jì)算分析。在實(shí)驗(yàn)給定的低速流動(dòng)條件下，連續(xù)水膜厚度及溪流厚度和寬度的計(jì)算結(jié)果在數(shù)量級(jí)和變化趨勢上均與實(shí)驗(yàn)數(shù)據(jù)吻合較好。限于篇幅，本文不再詳述，具體對(duì)比結(jié)果參見文獻(xiàn)［23］。

防冰表面溢流水的流動(dòng)受氣動(dòng)力和表面張力的影響最大，氣流速度直接影響水膜及溪流所受氣動(dòng)力的大小。云層參數(shù)不變的條件下，防冰部件表面的流動(dòng)狀態(tài)由于飛行速度的影響會(huì)呈現(xiàn)不同的特性，從而影響表面換熱過程。本文對(duì)不同來流速度下水膜流動(dòng)及破裂情況進(jìn)行了研究，分析了破裂前后防冰表面蒸發(fā)及對(duì)流換熱的特點(diǎn)。

計(jì)算仍采用弦長為0.101 m，攻角0°的NACA0012翼型的幾何模型，具體的計(jì)算狀態(tài)參數(shù)如表1所示。

表1 不同來流速度條件下的計(jì)算狀態(tài)參數(shù)Table 1 Calculation state parameter at different free stream velocities

4.1 水滴撞擊特性

不同來流速度條件下翼型表面的局部水收集系數(shù)分布如圖5所示。從整體上看，隨著來流速度的不斷增加，翼型表面的局部水收集系數(shù)無論在峰值還是在撞擊區(qū)域上都有所增加，表明來流速度高時(shí)將有更多的液態(tài)水到達(dá)翼型表面。

圖5 不同來流速度條件下局部水收集系數(shù)分布Fig.5 Local collection efficiency at different free stream velocities

4.2 連續(xù)水膜厚度分布

如圖6所示為不同來流速度條件下的連續(xù)水膜厚度在翼型上表面的分布情況。從整體上看，在不同來流速度條件下，連續(xù)水膜沿翼型表面的分布均呈現(xiàn)逐漸增大的變化趨勢。這是由于在撞擊區(qū)域內(nèi)，水滴撞擊質(zhì)量遠(yuǎn)大于表面蒸發(fā)質(zhì)量，導(dǎo)致溢流水量逐漸變大，從而使得水膜在沿表面流動(dòng)過程中逐漸增厚。此外，從圖6中還可以看出，隨著來流速度的增大，表面同一位置處水膜厚度有明顯的減小。這是由于前緣處水膜受氣流法向作用力影響較大。當(dāng)來流速度較大時(shí)，作用在水膜表面的法向氣動(dòng)力也較大，使得水膜變薄。

由于本文將水膜破裂點(diǎn)設(shè)置在水滴撞擊極限位置，由局部水收集系數(shù)的分布圖(圖5)可以看出，來流速度越大，撞擊區(qū)域越大，即水膜破裂位置隨來流速度的增大而延后。從圖6中可以看出，當(dāng)來流速度較大時(shí)，連續(xù)水膜在表面上的分布范圍更大，破裂位置更加遠(yuǎn)離駐點(diǎn)，此時(shí)破裂處的連續(xù)水膜厚度相對(duì)較大。這一特性將對(duì)破裂后溪流形態(tài)產(chǎn)生重要影響。

圖6 不同來流速度條件下連續(xù)水膜厚度分布Fig.6 W ater film thickness at different free stream velocities

4.3 連續(xù)水膜主要熱流分布

圖7 不同來流速度條件下水膜表面蒸發(fā)散熱熱流量Fig.7 Evaporation heat fluxes on water film surface at different free stream velocities

圖8 不同來流速度條件下水膜表面對(duì)流換熱熱流量Fig.8 Convection heat fluxes on water film surface at different free stream velocities

蒸發(fā)換熱與對(duì)流換熱是水膜表面主要的熱量損失［25］。圖7和圖8分別為不同來流速度條件下水膜表面蒸發(fā)熱流與對(duì)流換熱熱流沿表面分布圖。從圖中可以看出，隨著來流速度的增大，兩項(xiàng)熱流均有明顯增加。這是由于當(dāng)來流速度增加時(shí)，表面對(duì)流換熱系數(shù)會(huì)隨之增大，進(jìn)而影響表面蒸發(fā)和對(duì)流換熱熱流。

4.4 水膜破裂及溪流形態(tài)

圖9和10所示為不同來流速度條件下，水膜破裂處(即水滴撞擊極限處)溪流的厚度及寬度的分布情況。從圖上可以看出，相比于連續(xù)水膜，在相同來流速度條件下，破裂處溪流的厚度和寬度都有所增加;同時(shí)，當(dāng)來流速度不斷增大時(shí)，溪流的厚度和寬度也都隨之增大。正如上文所述，來流速度增大，表面溢流水量增加，在破裂處連續(xù)水膜厚度相對(duì)較大，因此溪流的厚度會(huì)有所增加。同樣，溪流寬度也隨之增大。但從計(jì)算結(jié)果的數(shù)值來看，幾個(gè)不同來流速度條件下的溪流厚度和寬度差別并不是很大。

圖9 不同來流速度條件下破裂處溪流厚度分布Fig.9 Rivulet thickness at breakup point at different free stream velocities

圖10 不同來流速度條件下破裂處溪流寬度分布Fig.10 Rivulet w idth at breakup point at different free stream velocities

4.5 防冰溢流水表面換熱特性

溢流水的流動(dòng)形態(tài)會(huì)影響表面換熱特性。圖11和12所示為不同來流速度條件下，三種不同溢流水流動(dòng)形態(tài)表面蒸發(fā)和對(duì)流換熱熱流的分布情況。首先，在相同來流速度下，完全濕潤表面(連續(xù)水膜表面)的蒸發(fā)熱流總大于部分濕潤表面(破裂處溪流干/濕表面)，這是由于部分濕潤表面中的干表面無蒸發(fā)熱流。其次，這兩種狀態(tài)下的蒸發(fā)熱流隨著來流速度的增加均呈現(xiàn)明顯增大的變化趨勢。這是由于來流速度較大時(shí)，表面?zhèn)髻|(zhì)系數(shù)相應(yīng)增大，從而影響表面蒸發(fā)熱流。

表面對(duì)流換熱熱流隨來流速度增加的變化趨勢與蒸發(fā)熱流相同，而對(duì)于同一來流速度下三種不同溢流水形態(tài)條件下表面對(duì)流換熱熱流的區(qū)別并沒有蒸發(fā)熱流那么明顯。從圖12還可以看出，相同速度下，完全干表面(無溢流水)的對(duì)流換熱熱流相對(duì)較大，部分濕潤表面次之，完全濕潤表面對(duì)流換熱量最小。這說明，盡管表面溢流水形成的水膜/溪流厚度很小，仍會(huì)對(duì)對(duì)流換熱表面的溫度產(chǎn)生一定的影響，只是影響量有限。

圖11 不同來流速度條件下表面蒸發(fā)熱流量Fig.11 Evaporation heat fluxes on anti-icing surface at different free stream velocities

圖12 不同來流速度條件下不同表面對(duì)流換熱熱流量Fig.12 Convection heat fluxes on anti-icing surface at different free stream velocities

5 結(jié)論

基于空氣-水兩層相互作用的質(zhì)量、動(dòng)量和能量守恒建立的防冰表面溢流水水膜流動(dòng)換熱及破裂的數(shù)學(xué)模型，研究分析來流速度對(duì)表面連續(xù)水膜厚度及主要散熱項(xiàng)、破裂點(diǎn)溪流厚度與寬度的影響。同時(shí)，對(duì)三種不同溢流水流動(dòng)形態(tài)下的表面蒸發(fā)與對(duì)流換熱熱流進(jìn)行對(duì)比，分析溢流水對(duì)表面換熱特性的影響。主要結(jié)論可以歸納為:

1)撞擊區(qū)域內(nèi)，在撞擊區(qū)域內(nèi)的蒸發(fā)質(zhì)量遠(yuǎn)小于撞擊質(zhì)量，因此溢流水質(zhì)量會(huì)不斷積累，使得連續(xù)水膜厚度沿表面的分布呈現(xiàn)不斷增大的變化趨勢;

2)來流速度的增加使得表面撞擊區(qū)域以及撞擊水量增加，表現(xiàn)為連續(xù)水膜在表面的鋪展面積及厚度均隨之增加;

3)來流速度對(duì)蒸發(fā)和對(duì)流換熱熱流的影響主要取決于換熱系數(shù)，來流速度增大，表面換熱劇烈，從而使得表面蒸發(fā)熱流和對(duì)流換熱熱流增加;

4)由于來流速度的增加，連續(xù)水膜在破裂處的厚度也會(huì)變大，使得破裂點(diǎn)溪流厚度和寬度增大;

5)在撞擊區(qū)域之外，翼型表面考慮溢流水的存在使得表面換熱過程蒸發(fā)熱流增加，而當(dāng)溢流水流動(dòng)形態(tài)由連續(xù)水膜破裂成溪流時(shí)，表面蒸發(fā)熱流有所下降;對(duì)于表面對(duì)流換熱而言，無溢流水流動(dòng)的表面熱流量略有增加，而完全濕潤表面換熱熱流量最小。

致謝:論文得到了上海交通大學(xué)機(jī)械與動(dòng)力工程學(xué)院工程熱物理研究所陳勇老師、雷桂林博士的支持和幫助，這里對(duì)他們的幫助表示感謝。

［1］Melody J W.In-flight characterization of aircraft icing［D］.Illinois: University of Illinois at Urbana-Champaign，2004.

［2］Politovich M K.Aircraft icing［J］.Encyclopedia of Atmospheric Sciences，2003:68-75.

［3］Qiu X G，Han F H.Aircraft anti-icing system［M］.Beijing: Aviation Industry Press，1985:44-58.(in Chinese)裘燮綱，韓鳳華.飛機(jī)防冰系統(tǒng)［M］.北京:航空工業(yè)出版社，1985:44-58.

［4］Addy H E.Ice accretions and icing effects for modern airfoils［R］.NASA TP 2000-210031.

［5］Messinger B L.Equilibrium temperature of an unheated icing surface as a function of air speed［J］.Journal of the Aeronautical Sciences，1953，20(1):29-42.

［6］Mac A C.Numerical simulation of airfoil ice accretion［R］.AIAA 83-0112.

［7］Yi X，Zhu G L.Computation of glaze ice accretion on airfoil［J］.Acta Aerodynamica Sinica，2004，22(4):490-493.(in Chinese)易賢，朱國林.考慮傳質(zhì)傳熱效應(yīng)的翼型積冰計(jì)算［J］.空氣動(dòng)力學(xué)學(xué)報(bào)，2004，22(4):490-493.

［8］Hedde T，Guffond D.ONERA three-dimensional icing model［J］.AIAA Journal，1995，33(6):1038-1045.

［9］Honsek R，Habashi W G.FENSAP-ICE:eulerian modeling of droplet impingement in the SLD regime of aircraft icing［R］.AIAA 2006-465.

［10］Morency F，Tezok F，Paraschivoiu I.Anti-icing system simulation using CANICE［J］.Journal of Aircraft，1999，36(3):999-1006.

［11］Al-Khalil K M，Wright W B，Miller D R.Validation of NASA thermal ice protection computer codes:part 3-the validation of ANTICE［R］.AIAA 1997-0051.

［12］Dong W，Zhu J，Zheng M，et al.Thermal analysis and testing of non-rotating cone leading edge hot air anti-icing system［J］.Journal of Propulsion and Power，2015，31(1):896-902.

［13］Dong W，Zhu J，Zhou Z X，et al.Heat transfer and temperature analysis of an aero-engine strut under icing condition［J］.Journal of Aircraft，2015，52(1):216-225.

［14］El-Genk M S，Saber H H.Minimum thickness of a flowing down liquid film on a vertical surface［J］.International Journal of Heat and Mass Transfer，2001，44(15):2809-2825.

［15］El-Genk M S，Saber H H.An investigation of the breakup of an evaporating liquid film，falling down a vertical，uniformly heated wall［J］.Journal of Heat Transfer，2002，124(1):39-50.

［16］Zhou X Q，Chen P L.Performance of the water flow film along a vertical plate［J］.Journal of Tongji University，1995，23(6):621-626.(in Chinese)周孝清，陳沛霖.豎直平板上的水膜流動(dòng)特征［J］.同濟(jì)大學(xué)學(xué)報(bào):自然科學(xué)版1995，23(6):621-626.

［17］Ye X M，Yan J G，Li C X，et al.Breakdown model of locally heated and subcooled films driving by interfacial shear and gravity［J］.Journal of North China Electric Power University，2010，37 (2):68-73.(in Chinese)葉學(xué)明，閆俊剛，李春曦，等.界面切應(yīng)力和重力驅(qū)動(dòng)下受熱過冷液膜的破斷模型研究［J］.華北電力大學(xué)學(xué)報(bào)，2010，37 (2):68-73.

［18］Al-Khalil K M.Numerical simulation of an aircraft anti-icing system incorporating a rivulet model for the runback water［D］.Toledo: The University of Toledo，1991.

［19］Silva G A L，Silvares O，Zerbini G J，et al.Numerical simulation of airfoil thermal anti-ice operation，part 1:mathematical modelling［J］.Journal of Aircraft，2007，44(2):627-633.

［20］Silva G A L，Silvares O，Zerbini G J，et al.Numerical simulation of airfoil thermal anti-ice operation，part 2:implementation and results［J］.Journal of Aircraft，2007，44(2):634-641.

［21］Silva G A L，Silvares O，Zerbini G J，et al.Differential boundarylayer ananlysis and runback water flow model applied to flow around airfoils with thermal anti-ice［R］.AIAA 2009-3967.

［22］Thundat T，Oden P I，Warmack R J.Microcantilever sensors［J］.Microscale Thermophysical Engineering，1997，1(3):185-199.

［23］Dong W，Zheng M，Zhu J，et al.Calculation and analysis of runback water flow on anti-icing surface［R］.AIAA 2015-0538.

［24］Zhang K，Blake J，Rothmayer A，et al.An experimental investigation on wind-driven rivulet/film flows over a NACA0012 airfoil by using digital image projection technique［R］.AIAA 2014-0741.

［25］Zheng M，Dong W，Lei G L，et al.Study of the flow and heat transfer of water film on hot air anti-icing airfoil surface［R］.IHTC 15-8671.

Influence of free stream velocity on runback water flow and heat transfer on anti-icing surface

Zheng Mei，Dong Wei*，Zhu Jianjun，Guo Zhiqiang
(School of Mechanical Engineering，Shanghai Jiao Tong University，Shanghai 200240，China)

The purpose of this paper is to investigate the effect of free stream velocity on the runback water flow and heat transfer on the anti-icing surface.Based on the mass，momentum and energy conservations of the runback water flow and the air flow，a mathematical model of the runback water film flow and rivulet flow was developed to investigate the effect of the free stream velocity on the heat and mass transfer on the anti-icing surface.The computation analysis indicates that the water film thickness at the same position on the anti-icing surface decreases with the free stream velocity increasing，and the rivulet thickness and width at the breakup point increase due to larger water film thickness at higher free stream velocity.Meanwhile，as the main heat losses on the anti-icing surface，the evaporation heat flux and the convection heat flux on the water film surface increase with the free stream velocity increasing.In addition，the characteristics of the heat and mass transfer on the dry surface，the fully wet surface and the partially wet surface were also investigated.The results show that the patterns of the runback water have some influence on the heat and mass transfer on the anti-icing surface.

anti-icing;runback water flow;flow and heat transfer;water film;free stream velocity

V211.3

A

10.7638/kqdlxxb-2015.0219

0258-1825(2016)03-0295-07

2015-12-21;

2015-12-23

國家自然科學(xué)基金(11272212，11572195);國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(2015CB755800)

鄭梅(1991-)，女，浙江衢州人，博士研究生，研究方向:飛機(jī)及發(fā)動(dòng)機(jī)結(jié)冰與防冰.E-mail:may-zheng@sjtu.edu.cn

董威*(1970-)，男，教授，博士生導(dǎo)師，研究方向:飛機(jī)及發(fā)動(dòng)機(jī)結(jié)冰與防冰.E-mail:wdong@sjtu.edu.cn

鄭梅，董威，朱劍鋆，等.來流速度對(duì)防冰表面溢流水流動(dòng)換熱的影響［J］.空氣動(dòng)力學(xué)學(xué)報(bào)，2016，34(3):295-301.

10.7638/kqdlxxb-2015.0219 Zheng M，Dong W，Zhu J J，et al.Influence of free stream velocity on runback water flow and heat transfer on anti-icing surface［J］.Acta Aerodynamica Sinica，2016，34(3):295-301.