西北工業(yè)大學(xué)力學(xué)與土木建筑學(xué)院飛行器可靠性工程研究所 姚會(huì)舉 劉永壽 張 峰 南 華

飛機(jī)以亞音速在含過(guò)冷水滴的云層中飛行時(shí)，由于云層中過(guò)冷水滴對(duì)機(jī)翼的撞擊作用，可能發(fā)生機(jī)翼積冰現(xiàn)象，而機(jī)翼積冰會(huì)對(duì)飛機(jī)產(chǎn)生很大的危害，不僅增加飛機(jī)負(fù)重，而且會(huì)使飛機(jī)的空氣動(dòng)力學(xué)性能惡化，升力減小，阻力增大，輕者影響飛機(jī)的穩(wěn)定性和操作性，重者甚至?xí)l(fā)機(jī)毀人亡的飛行事故[1-2]。

防冰系統(tǒng)性能的好壞直接影響著飛機(jī)的防/除冰能力，進(jìn)而影響飛行安全，其設(shè)計(jì)一直是生命保障和環(huán)控系統(tǒng)設(shè)計(jì)的一個(gè)重要組成部分[3]。飛機(jī)防冰系統(tǒng)按防/除冰方式的不同可分為機(jī)械除冰、液體防冰、電熱防/除冰和熱氣防/除冰。由于穩(wěn)定性、簡(jiǎn)易性以及經(jīng)濟(jì)性等方面的優(yōu)勢(shì)，大部分現(xiàn)役客機(jī)機(jī)翼采用熱氣防冰系統(tǒng)。然而，熱氣防冰系統(tǒng)存在著防冰熱氣有效利用率不高，大量熱氣未經(jīng)充分利用即被排出的缺陷。因此，為提高飛機(jī)性能，研究防冰腔結(jié)構(gòu)對(duì)熱氣防冰系統(tǒng)性能的影響有重要意義。

有關(guān)機(jī)翼氣熱防/除冰的理論、試驗(yàn)以及仿真分析也得到了廣泛開(kāi)展。Planquart等[4]采用熱薄膜方法，借助紅外溫度記錄儀試驗(yàn)測(cè)得了叉排直噴式熱氣防冰腔內(nèi)蒙皮表面的換熱系數(shù)分布情況，并與Fluent仿真結(jié)果一致，得出防冰腔的換熱性能與射流孔的弦向和展向分布間距以及熱氣的沖擊雷諾數(shù)有關(guān)。Farooq等[5]利用一種基于基因算法的優(yōu)化程序確定了熱氣防冰腔的典型參數(shù)，并通過(guò)算例驗(yàn)證了算法的可行性。卜雪琴等[6]對(duì)熱氣防冰系統(tǒng)中短小銷(xiāo)釘對(duì)強(qiáng)化傳熱性能的影響進(jìn)行了研究。

結(jié)構(gòu)參數(shù)對(duì)直噴式雙蒙皮防冰腔防冰熱性能的影響目前在國(guó)內(nèi)的研究還不多見(jiàn)。本文通過(guò)數(shù)值模擬來(lái)評(píng)估射流孔到蒙皮沖擊面距離H和射流孔孔徑d對(duì)防冰熱性能的影響。采用計(jì)算流體力學(xué)CFD方法，對(duì)不同結(jié)構(gòu)尺寸的直噴式雙蒙皮防冰腔進(jìn)行了仿真，目的是獲得不同結(jié)構(gòu)尺寸因素對(duì)強(qiáng)化傳熱系數(shù)和防冰效率的影響，為熱氣防冰系統(tǒng)的設(shè)計(jì)和優(yōu)化提供指導(dǎo)。

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

由于防冰腔內(nèi)為射流孔高速?lài)姵龅臒釟猓瑹釟庠诜辣粌?nèi)流動(dòng)，受腔內(nèi)結(jié)構(gòu)、熱氣速度的影響，其流動(dòng)形式為湍流，為準(zhǔn)確模擬熱氣在防冰腔內(nèi)的流動(dòng)，仿真計(jì)算采用Spalart-Allaras(簡(jiǎn)稱(chēng)S-A)湍流模型，S-A湍流模型對(duì)沖擊射流曲面的流動(dòng)問(wèn)題具有計(jì)算簡(jiǎn)單、與試驗(yàn)結(jié)果吻合較好等優(yōu)點(diǎn)[7-8]。S-A湍流模型是渦粘性模型,其核心思想是通過(guò)求解中間變量 的輸運(yùn)方程來(lái)獲得湍流粘度μt的輸運(yùn)方程為[9]：

求得后，湍流粘度μt通過(guò)濕潤(rùn)系數(shù)獲得：

2 計(jì)算模型

本文針對(duì)直噴式雙蒙皮防冰腔結(jié)構(gòu)模型開(kāi)展分析，引氣系統(tǒng)從發(fā)動(dòng)機(jī)引來(lái)的熱氣經(jīng)過(guò)管路流到機(jī)翼前緣中的笛形管內(nèi)，熱氣從笛形管壁上的射流孔噴射到防冰腔前腔，加熱蒙皮表面，熱氣經(jīng)過(guò)雙蒙皮防冰通道后流到后腔，并從后腔的排氣口排出。為方便網(wǎng)格劃分，將防冰腔出口簡(jiǎn)化為方形出口，防冰腔示意圖如圖1所示。

圖1 防冰腔示意圖Fig.1 Diagram of anti-icing cavity

蒙皮厚度1.5mm，笛形管直徑40mm，上面分布雙排直徑為d的圓形射流孔，射流孔間距為25mm，將熱氣從射流孔噴出，沿軸線(xiàn)到達(dá)蒙皮距離記為H，d和H如圖2所示，其尺寸參數(shù)范圍如表1所示。射流孔按照排列方式可以分為平行排列和交叉排列兩種，如圖3所示。

為了研究射流孔噴射距離H、射流孔孔徑d以及射流孔的排列形式對(duì)強(qiáng)化防冰腔傳熱系數(shù)和防冰效率的影響，設(shè)計(jì)了3組幾何尺寸不同的防冰腔模型，分別為：(1)H相同但d與H比值不同，即射流孔的孔徑不同的平行排列； (2)兩者之比相同但H不同的平行陣列；(3)對(duì)應(yīng)相同d與H比值的交叉排列。具體參數(shù)如表2所示。

圖2 d、H 結(jié)構(gòu)示意圖Fig.2 Structure diagram of d, H

表1 幾何尺寸參數(shù)范圍

圖3 笛形管結(jié)構(gòu)尺寸示意圖Fig.3 Diagram of piccolo tube structure size

表2 模型具體參數(shù)

防冰腔流域計(jì)算模型采用結(jié)構(gòu)體網(wǎng)格，如圖4所示，在蒙皮、防冰腔前腔近壁面部分采用附面層網(wǎng)格，保障蒙皮網(wǎng)格的密度與質(zhì)量，并在射流口區(qū)域進(jìn)行網(wǎng)格局部加密，最終網(wǎng)格總數(shù)在63萬(wàn)左右。這樣既保證了關(guān)注地方的網(wǎng)格細(xì)化，同時(shí)網(wǎng)格總量不會(huì)太大，加快流場(chǎng)計(jì)算速度。

圖4 防冰腔網(wǎng)格Fig.4 Anti-icing cavity mesh

3 邊界條件

將射流孔處設(shè)置為壓力入口，射流孔處壓力0.2MPa，熱氣溫度200℃；出口處設(shè)置為與環(huán)境壓力相同的壓力出口，外部環(huán)境溫度設(shè)置為-10℃；兩側(cè)擋板設(shè)為對(duì)稱(chēng)邊界；防冰腔內(nèi)、外蒙皮和擋板設(shè)為壁面邊界條件。

防冰腔采用S-A湍流模型，比較適合沖擊射流曲面的流動(dòng)計(jì)算，開(kāi)啟能量方程，控制方程的離散采用控制容積法，對(duì)流項(xiàng)差分格式采用二階迎風(fēng)格式，采用Simple算法求解方程。監(jiān)控出口溫度變化以及殘差來(lái)判斷計(jì)算的收斂，收斂精度為105。

仿真計(jì)算中，由于流場(chǎng)結(jié)構(gòu)復(fù)雜，不易收斂，通過(guò)調(diào)整松弛因子大小來(lái)控制收斂。判定是否收斂一般有殘差監(jiān)視和表面溫度平衡兩種方式。計(jì)算完成后，通過(guò)檢驗(yàn)計(jì)算流域的質(zhì)量和能量是否守恒來(lái)進(jìn)一步驗(yàn)證收斂性。

4 計(jì)算結(jié)果分析與討論

4.1 數(shù)據(jù)處理

Fluent仿真計(jì)算收斂后得到防冰腔外蒙皮表面溫度分布受熱氣流對(duì)防水腔前緣內(nèi)表面的沖擊影響，在前緣駐點(diǎn)附近蒙皮溫度達(dá)到最高值，蒙皮溫度沿弦向遞減，防冰腔下表面末端達(dá)到最低值，溫差達(dá)60K左右。上表面受雙蒙皮通道影響，在進(jìn)入雙蒙皮通道之后的防冰區(qū)域蒙皮溫度有一個(gè)明顯上升趨勢(shì)，而后再沿弦向遞減，整體來(lái)看，上蒙皮表面溫度比下表面溫度高。由于氣流在雙蒙皮通道與蒙皮接觸更充分，換熱效果明顯優(yōu)于下表面區(qū)域，根據(jù)參考文獻(xiàn)[10]中相似工況下仿真結(jié)果，可知本文仿真具有一定的參考價(jià)值。

為了對(duì)防冰腔防冰性能進(jìn)行直觀描述，本文引用防冰效率η的概念[11-12]，防冰效率是指系統(tǒng)所提供的熱氣能量中用于加熱防冰表面的百分比，計(jì)算公式如下：

式中，Tinlet、Toutlet分別為防冰系統(tǒng)供氣溫度和排氣溫度，Twall為防冰區(qū)蒙皮平均溫度。

4.2 結(jié)果分析

4.2.1 幾何尺寸的影響

分析1～5號(hào)模型的結(jié)果數(shù)據(jù)來(lái)研究射流孔孔徑對(duì)流動(dòng)換熱的影響。圖5、圖6分別為蒙皮表面平均傳熱系數(shù)h和模型的熱氣有效利用率η與d/H的關(guān)系。從圖5中可以看出，隨著d/H的增大，模型強(qiáng)化蒙皮表面平均傳熱系數(shù)明顯增大，模型的熱氣有效利用率η不斷下降。也就是說(shuō)，在相同射流孔到蒙皮沖擊面距離的情況下，射流孔孔徑越大，強(qiáng)化傳熱效果越強(qiáng)，但是用于防冰的熱氣有效利用率也越低。這表明，雖然增大射流孔孔徑即增大熱氣量可以使蒙皮的傳熱效果得到強(qiáng)化，但用于防冰的熱氣比例卻在降低，過(guò)低的防冰熱效率會(huì)加大發(fā)動(dòng)機(jī)的負(fù)擔(dān)。

圖5 平均傳熱系數(shù)h隨d / H的變化Fig.5 Varitions of average heat transfer coefficient h with d / H

圖6 熱氣有效利用率η隨d / H的變化Fig.6 Varitions of thermal efficiency η with d / H

表3 1～5號(hào)模型對(duì)比3號(hào)模型換熱和熱效率的變化

表3中所示為以3號(hào)模型的數(shù)據(jù)為基準(zhǔn)，其他模型內(nèi)表面?zhèn)鳠嵯禂?shù)和熱氣效率相對(duì)變化值，負(fù)號(hào)表示相對(duì)減小?？梢钥闯?，5號(hào)模型表面?zhèn)鳠嵯禂?shù)增大了45.9%，但熱氣效率卻減小了12.58%，說(shuō)明為了增大表面換熱系數(shù)而采用如此大的d/H帶來(lái)熱氣效率的損失過(guò)大，該模型是不合適的。強(qiáng)化換熱和降低熱效率之間的平衡點(diǎn)應(yīng)該在實(shí)際應(yīng)用中根據(jù)經(jīng)驗(yàn)需求數(shù)據(jù)酌情而定，d/H要適中。

分析6～10號(hào)模型的結(jié)果數(shù)據(jù)來(lái)研究d/H相同下射流孔到蒙皮沖擊面距離的變化對(duì)流動(dòng)換熱的影響。表4為以3號(hào)模型的數(shù)據(jù)為基準(zhǔn)，其他模型蒙皮表面?zhèn)鳠嵯禂?shù)和熱氣效率的相對(duì)變化值?？梢钥闯?，在d/H相同下傳熱系數(shù)和熱氣效率相對(duì)變化不大，幅度都在6%以?xún)?nèi)，可知d與H的比值是影響防冰腔性能的一個(gè)重要因素。

表4 6～10號(hào)模型對(duì)比3號(hào)模型換熱和熱效率的變化

4.2.2 射流孔排列形式的影響

表5給出了對(duì)應(yīng)結(jié)構(gòu)尺寸相同的射流孔平行、交叉排列形式時(shí)的強(qiáng)化換熱和熱氣效率變化百分比?？梢钥吹?，與平行排列相比，交叉排列形式的模型的表面換熱系數(shù)與傳熱效率都略大，這是由于交叉排列射流孔噴出的熱氣與防冰腔內(nèi)表面的沖擊面積更大，熱交換更充分，因此同樣條件下交叉排列的防冰效果優(yōu)于平行排列。

表5 平行和交叉排列對(duì)強(qiáng)化換熱和效率的影響

5 結(jié)論

本文對(duì)熱氣防冰系統(tǒng)中笛形管位置與射流孔孔徑強(qiáng)化傳熱效果進(jìn)行了研究，定量分析了射流孔噴射距離、射流孔孔徑以及射流孔排列形式對(duì)換熱以及傳熱效率的影響，為防冰腔的改進(jìn)設(shè)計(jì)提供參考。得出結(jié)論如下：

(1)射流孔孔徑d和射流孔與蒙皮沖擊面距離H的比值能夠較大程度地影響換熱與傳熱效率。d/H越大，即H為定值時(shí)，射流孔孔徑d越大，蒙皮換熱系數(shù)越高。但是，對(duì)于較大的射流孔，防冰系統(tǒng)的傳熱效率也相應(yīng)的降低。

(2)交叉排列形式模型的表面換熱系數(shù)與傳熱效率都略大，防冰效果優(yōu)于同等條件下平行排列的模型。

[1] 管寧. 三維機(jī)翼防冰熱載荷的數(shù)值模擬[D]. 南京： 南京航空航天大學(xué), 2007.

[2] 卜雪琴, 林貴平, 郁嘉. 三維內(nèi)外熱耦合計(jì)算熱氣防冰系統(tǒng)表面溫度. 航空動(dòng)力學(xué)報(bào), 2009, 24(11)： 2495-2500.

[3] Pellissier M, Wagdi G, Habashi. Design optimization of hot-air anti-icing systems by FENSAP-ICE. Journal of Aircraft, 2010, 48(4)： 1238-1251.

[4] Planquart P, Vanden G B. Experimental and numerical optimization of a wing leading edge hot air anti-icing system. Journal of Aircraft, 2005, 43(4)： 1259-1277.

[5] Saeed F, Paraschivoiu I. Optimization of a hot-air anti-icing system//AIAA. Proceedings of 41st AIAA Aerospace Meeting & Exhibit,Reston：AIAA, 2003：0733.

[6] 卜雪琴, 林貴平. 基于CFD的水收集系數(shù)及防冰表面溫度預(yù)測(cè). 北京航空航天大學(xué)學(xué)報(bào), 2007, 33(10)： 1182-1185.

[7] Papadakis M, Wong S J. Parametric investigation of a bleed air ice protection system//AIAA. Proceedings of 44th AIAA Aerospace Sciences Meeting and Exhibit，Reston：AIAA, 2006：1013.

[8] Mathieu F, Saeed F. Numerical heat transfer correlation for array of hot-air jets impinging on 3-dimensional concave surface. Journal of Aircraft, 2005,42(3)： 665-670.

[9] 陶文銓. 數(shù)值傳熱學(xué).第2版.西安：西安交通大學(xué)出版社,2001.

[10] Farroq Saeed, Ahmad Z.Al-Garni. Numerical Simulation of Surface Heat Tranfer from an Array of Hot-Air Jets//AIAA. Proceedings of 25th AIAA Applied Aerodynamics Conference, Reston：AIAA, 2007：4287.

[11] 卜雪琴,郁嘉,林貴平,等.機(jī)翼氣熱防冰系統(tǒng)設(shè)計(jì). 北京航空航天大學(xué)學(xué)報(bào), 2010, 36(8)： 927-930.

[12] Hugh H, Liu T, Hua J.Three-dimensional intergrated thermodynamic simulation for wing anti-icing system. Journal of aircraft,2004, 41(6)：1291-1297.