胡金源，金志光，田維康

(南京航空航天大學(xué) 能源與動力學(xué)院，江蘇 南京 210016)

0 引言

航空發(fā)動機唇口結(jié)冰會減小進(jìn)氣流通面積并造成流場畸變，而且冰塊的掉落可能損傷風(fēng)扇，對飛行安全的影響較大，所以目前飛機上廣泛采用了熱氣防冰系統(tǒng)[1]。熱氣防冰系統(tǒng)所引入的熱空氣一般來源于高壓壓氣機，引氣流量可達(dá)發(fā)動機總空氣流量的10%，對發(fā)動機的性能會產(chǎn)生較大影響[2]。

為獲得換熱性能最優(yōu)的防冰腔，國內(nèi)外對各型防冰腔結(jié)構(gòu)參數(shù)開展了研究。彭瓏[3]和姚會舉[4]等對防冰腔的射流孔排列方式和排數(shù)等結(jié)構(gòu)參數(shù)進(jìn)行了研究，分析不同防冰腔內(nèi)壁面的對流傳熱系數(shù)分布。張靖周等[5]研究了笛形管結(jié)構(gòu)參數(shù)對防冰腔蒙皮凹腔表面溫度分布的影響。SAEED F等[6]基于對防冰腔結(jié)構(gòu)參數(shù)影響的分析，開發(fā)了優(yōu)化程序來獲取最佳的結(jié)構(gòu)參數(shù)。郭之強等[7]設(shè)計了表面凸起結(jié)構(gòu)，研究了射流角度對熱氣防冰效果的影響。蘆婭妮等[8]在開孔面積相同的情況下，應(yīng)用CFD技術(shù)研究了大孔形狀對換熱器傳熱系數(shù)和壓降的影響。

目前關(guān)于微引射式防冰腔的熱氣防冰性能研究還較少，本文將研究分析微引射式防冰腔主要結(jié)構(gòu)參數(shù)對防冰熱效率和蒙皮內(nèi)表面溫度分布的影響，為微引射式防冰腔的結(jié)構(gòu)優(yōu)化提供借鑒。

1 微引射式防冰腔模型

微引射式防冰腔的基準(zhǔn)模型如圖1所示。防冰腔分為前腔和后腔，其中前腔包括射流孔出口、引射式混合腔和波紋板通道，后腔為波紋板通道出口到笛形管和熱氣出口區(qū)域?！繱為唇口駐點處到后腔體擋板的弧長。高溫?zé)釟鈴纳淞骺讎姵?，引射后腔體余氣經(jīng)過混合腔通道以射流沖擊換熱的形式加熱發(fā)動機蒙皮內(nèi)表面。然后高溫氣體進(jìn)入波紋板通道，繼續(xù)加熱唇口兩側(cè)蒙皮內(nèi)表面，最后從波紋板通道出口流出，進(jìn)入后腔變?yōu)橛鄽?。一部分余氣再次進(jìn)入引射式混合腔，一部分余氣從排氣孔排出。

圖1 微引射式防冰腔基準(zhǔn)模型

設(shè)置微引射式防冰腔基準(zhǔn)模型展向?qū)挾葹?20mm，笛形管直徑為38mm，射流孔直徑和間距分別為1.6mm和15mm，蒙皮厚度為2mm，混合腔當(dāng)量直徑D和長度L分別為5mm和35mm。波紋板截面如圖2所示，通道中心間隔a為30mm，通道寬度b為10mm，通道截面高度h為4mm，波紋板與發(fā)動機蒙皮狹縫通道的高度為1 mm。

圖2 波紋板通道截面

本文分別改變射流孔直徑、射流孔間距、混合腔長徑比和波紋板通道出口高度，得到了如表1所示的12種微引射式防冰腔構(gòu)型。

表1 微引射式防冰腔模型的結(jié)構(gòu)參數(shù)

2 網(wǎng)格劃分和計算模型設(shè)置

設(shè)置矩形外流場長度為10倍模型長度，高度為5倍模型高度，寬度與模型保持一致。在計算域中選取防冰腔的M1、M2兩個截面進(jìn)行后續(xù)分析。其中M1截面位于射流孔與混合腔中心的z=0mm處；M2截面位于射流孔和波紋板通道的中心。采用FLUENT Meshing劃分計算網(wǎng)格，對射流孔、混合腔和波紋板通道處的流體域進(jìn)行加密，如圖3所示。在壁面處劃分邊界層網(wǎng)格，保證壁面y+在20左右，以六面體核心的網(wǎng)格生成方式形成非結(jié)構(gòu)網(wǎng)格。生成網(wǎng)格總量分別為180萬、260萬、330萬、460萬、550萬共5套網(wǎng)格進(jìn)行數(shù)值模擬，發(fā)現(xiàn)網(wǎng)格數(shù)量達(dá)到330萬后，內(nèi)蒙皮表面平均溫度結(jié)果基本一致，最終選取330萬網(wǎng)格總量進(jìn)行數(shù)值模擬。

圖3 模型局部網(wǎng)格劃分

數(shù)值計算采用FLUENT軟件，選擇基于密度耦合的求解器，湍流模型使用k-ε模型。流體域介質(zhì)選擇Sutherland黏性條件的理想氣體。笛形管壁面為463K恒溫邊界，流體域兩側(cè)為對稱邊界。外流場為104 360 Pa壓力進(jìn)口和101 325Pa壓力出口，溫度為265K。射流孔總?cè)肟谫|(zhì)量流量為0.00505kg/s，溫度為463K。監(jiān)控殘差下降3個數(shù)量級，且蒙皮內(nèi)表面平均溫度收斂穩(wěn)定后完成數(shù)值計算。

采用文獻(xiàn)[9]和文獻(xiàn)[10]的模型進(jìn)行驗證，本文與文獻(xiàn)中所得防冰腔模型特定截面的對流換熱系數(shù)h對比如圖4所示。結(jié)果表明本文使用的網(wǎng)格劃分方式和FLUENT設(shè)置可以獲得較為準(zhǔn)確的結(jié)果。

圖4 驗證模型壁面對流換熱系數(shù)對比

3 數(shù)值計算結(jié)果與分析

引入防冰熱效率η來分析不同結(jié)構(gòu)微引射式防冰腔的熱氣防冰性能，具體表達(dá)式如下：

(1)

其中：Thotin為進(jìn)口熱氣平均溫度；Thotout為排氣孔出口熱氣平均溫度；Twall為蒙皮表面平均溫度。

首先得到如圖5所示的G1模型蒙皮內(nèi)表面溫度分布。從圖中可看出，沖擊射流駐點區(qū)域的溫度最高，達(dá)到了334.52K。蒙皮溫度隨著熱氣流向射流交互區(qū)發(fā)展逐漸降低，波紋板通道與雙蒙皮狹縫通道對應(yīng)的表面之間存在約2K的溫差。

圖5 G1模型蒙皮內(nèi)表面溫度分布

圖6為G1模型在M2截面射流孔處的流場。由于高速氣流的卷吸效應(yīng)，后腔中的余氣被卷吸加速，和主流一起進(jìn)入混合腔，沖擊加熱蒙皮內(nèi)表面。

圖6 G1模型在M2截面射流孔處流場

3.1 射流孔直徑的影響

在其余結(jié)構(gòu)參數(shù)不變的情況下改變射流孔的直徑，得到如圖7所示的G1-G4模型在M2截面處壁面溫度分布。G1-G4模型的射流孔出口速度分別為397.61m/s、504.97m/s、307.52m/s和238.17m/s，對應(yīng)的蒙皮內(nèi)表面駐點溫度分別為333.87K、335.19K、332.45K和330.12K，可知射流沖擊速度與對應(yīng)的壁面駐點溫度成正相關(guān)關(guān)系。在進(jìn)口流量一致的情況下，G2模型的孔徑最小，射流孔入口速度和熱空氣進(jìn)入波紋板通道后的速度最大，所以G2模型的蒙皮內(nèi)表面平均溫度最高。

圖7 G1-G4模型M2截面壁面溫度分布

表2為根據(jù)式(1)計算出的不同孔徑防冰腔模型的防冰熱效率。隨著射流孔徑的減小，蒙皮表面平均溫度和防冰熱效率逐漸增大。

表2 不同射流孔直徑模型的防冰熱效率

3.2 射流孔間距的影響

以G1模型為基準(zhǔn)，保持其余結(jié)構(gòu)參數(shù)不變，改變射流孔的間距得到G5-G7模型。圖8、圖9所示為G5、G7模型的前腔M1截面流場圖。射流孔間距改變時，在展向120mm的防冰腔模型上射流孔的數(shù)量隨之改變，G1、G5-G7模型的射流孔數(shù)量分別為8個、6個、10個、12個。射流孔數(shù)量減小時每個射流孔的質(zhì)量流量增大，射流孔的出口流速也變大。G1、G5-G7對應(yīng)的射流孔出口最大速度分別為397.61m/s、485.67m/s、347.18m/s、299.14m/s。

圖8 G5模型M1截面處流場

圖9 G7模型M1截面處流場

不同射流孔間距的G1、G5-G7模型在M2截面蒙皮內(nèi)表面溫度分布如圖10所示。G5模型的射流孔出口速度最大，因此對應(yīng)駐點區(qū)域的溫度最高。但G5模型的射流孔間距過大，射流速度在M1截面的交互區(qū)耗散嚴(yán)重，熱空氣進(jìn)入波紋板通道時速度顯著下降，導(dǎo)致波紋板通道對應(yīng)的蒙皮表面溫度較低。

圖10 G1、G5-G7模型M2截面壁面溫度分布

不同孔間距防冰腔模型的防冰熱效率如表3所示?？煽闯?，隨著孔間距的增大，蒙皮表面平均溫度和防冰熱效率先增加后減小，4個模型中G1模型有最高的蒙皮表面平均溫度和防冰熱效率，表明存在一個最佳孔間距范圍使得防冰腔性能最好。

表3 不同射流孔間距模型的防冰熱效率

3.3 混合腔長徑比的影響

保持混合腔的當(dāng)量直徑為5mm，僅改變混合腔長度，得到G8-G10模型。G1、G8-G10模型在M2截面的蒙皮內(nèi)表面溫度分布如圖11所示。G8模型的射流孔與壁面之間的距離最近，因此有最高的駐點溫度為335.6K。G9模型的駐點溫度僅比G8模型減少1.19K，而波紋板通道對應(yīng)的蒙皮表面區(qū)域溫度比G8模型高2K～5K。

圖11 G1、G8-G10模型M2截面壁面溫度分布

表4所示為不同混合腔長徑比模型的防冰熱效率，隨著混合腔長徑比的增加，防冰熱效率先增加再減小。通過監(jiān)測混合腔截面的質(zhì)量流量發(fā)現(xiàn)，混合腔長徑比增大時射流的引射效應(yīng)增強，能卷吸更多的余氣，使混合腔截面的質(zhì)量流量增加[11]。但長徑比增大到8以后混合腔截面的流量增加減慢，且長徑比過大會使駐點溫度降低，因此混合腔長徑比也存在一個最佳取值使防冰腔性能最好。

表4 不同混合腔長徑比模型的防冰熱效率

3.4 波紋板通道高度的影響

G1模型的波紋板為等截面通道，在其他結(jié)構(gòu)參數(shù)不變時，設(shè)置G11和G12模型的波紋板為變截面通道，波紋板通道的進(jìn)口高度為4mm，出口高度分別為2mm和1mm，熱氣通道的截面積逐漸變小。

數(shù)值模擬后得到G1、G11和G12模型在M2截面處壁面溫度分布如圖12所示。G1、G11和G12的駐點溫度相差很小，3個模型在波紋板入口處的蒙皮表面溫度幾乎相同，而G12在波紋板通道對應(yīng)的蒙皮表面溫度下降最慢。因此G12模型有最高的蒙皮表面平均溫度316.69K。

圖12 G1、G11和G12模型M2截面壁面溫度分布

變截面波紋板通道模型的防冰熱效率如表5所示，G11模型的防冰熱效率較等截面波紋板通道模型G1提高了2.91%，表明波紋板通道截面變化越劇烈則防冰熱效率越高。

表5 變截面波紋板通道模型的防冰熱效率

4 結(jié)語

本文研究分析了不同結(jié)構(gòu)參數(shù)對微引射式防冰腔的蒙皮內(nèi)表面溫度分布、防冰熱效率的影響規(guī)律，結(jié)果表明：

1)微引射式防冰腔的射流孔孔徑和波紋板通道出口高度的減小能提高防冰熱效率和蒙皮內(nèi)表面平均溫度，其中射流孔徑從2.2mm減小到1.3mm時，防冰熱效率提高了4.81%；波紋板通道出口高度從4mm減小到1mm時，防冰熱效率提高了2.91%。

2)微引射式防冰腔的射流孔間距減小時，蒙皮內(nèi)表面平均溫度和防冰熱效率先增大后減小；混合腔長徑比增加時，射流能引射更多余氣，但長徑比增大導(dǎo)致射流孔與蒙皮內(nèi)表面之間距離增加，蒙皮內(nèi)表面駐點溫度降低。因此射流孔間距和混合腔長徑比都存在一個最優(yōu)值使防冰腔換熱性能最好。