試論航空發(fā)動(dòng)機(jī)渦輪機(jī)匣傳熱分析技術(shù)

權(quán)少寧

（西安航空職業(yè)技術(shù)學(xué)院，陜西 西安 710089）

航空發(fā)動(dòng)機(jī)渦輪機(jī)匣能夠有效控制葉尖間隙，在進(jìn)行發(fā)動(dòng)機(jī)設(shè)計(jì)時(shí)，準(zhǔn)確預(yù)測(cè)運(yùn)轉(zhuǎn)溫度，獲取渦輪機(jī)匣在各種狀態(tài)下的變形量，能夠有效幫助葉尖間隙控制系統(tǒng)開(kāi)發(fā)。渦輪機(jī)匣是熱端部件，結(jié)構(gòu)十分復(fù)雜，包括各種各樣的流動(dòng)形式，即平板、凹槽、沖擊、多孔陣列等等，換熱規(guī)律也比較復(fù)雜。而因?yàn)闇u輪機(jī)匣的工作環(huán)境溫度比較高，流動(dòng)形式太過(guò)多元化，在溫度分析時(shí)，換熱模型實(shí)際上與具體結(jié)構(gòu)存在明顯差異，從而使得難以預(yù)測(cè)機(jī)匣溫度與變形，進(jìn)而不能確保葉尖間隙控制。因此，為了渦輪機(jī)匣溫度分析的準(zhǔn)確性得到提升，應(yīng)進(jìn)一步深入探究渦輪機(jī)匣表面換熱規(guī)律[1]。

1 航空發(fā)動(dòng)機(jī)渦輪機(jī)匣傳熱分析技術(shù)環(huán)節(jié)

基于高壓渦輪機(jī)匣進(jìn)行傳熱分析，其流動(dòng)傳熱分析過(guò)程比較繁雜，受各種因素影響，但是本質(zhì)上是流體動(dòng)力、熱交換理論、發(fā)動(dòng)機(jī)部件的有機(jī)結(jié)合。流動(dòng)特性是對(duì)發(fā)動(dòng)機(jī)二次流空氣的運(yùn)行狀態(tài)進(jìn)行探究，換熱特性是對(duì)發(fā)動(dòng)機(jī)零件表面與氣流熱交換原理進(jìn)行探究。流動(dòng)換熱的關(guān)鍵在于求解流場(chǎng)，現(xiàn)階段N-S 方程的準(zhǔn)確求解依舊是一大難題，需要以實(shí)驗(yàn)為載體加以彌補(bǔ)[2]。

充分了解高壓渦輪機(jī)匣和主動(dòng)間隙控制系統(tǒng)結(jié)構(gòu)特性，熟練掌握流動(dòng)傳熱的影響因素，有助于流動(dòng)與傳熱分析模型構(gòu)建。切實(shí)結(jié)合二次流空氣系統(tǒng)設(shè)計(jì)，進(jìn)行一維分析，從而獲取渦輪機(jī)匣氣體流動(dòng)環(huán)境狀態(tài)。就發(fā)動(dòng)機(jī)具體結(jié)構(gòu)與流動(dòng)狀況，就各流動(dòng)單元進(jìn)行準(zhǔn)則模塊化處理，以此構(gòu)成不同類(lèi)型節(jié)流單元串聯(lián)、并聯(lián)相組成的復(fù)雜空氣系統(tǒng)一維流體計(jì)算網(wǎng)絡(luò)，以進(jìn)一步做計(jì)算分析，進(jìn)而獲取高壓渦輪機(jī)匣流動(dòng)特性。在此基礎(chǔ)上，對(duì)渦輪機(jī)匣進(jìn)行CFD數(shù)值仿真研究，探究渦輪機(jī)匣流動(dòng)換熱分布的規(guī)律性，基于數(shù)值仿真驗(yàn)證換熱特性。數(shù)值仿真有利于獲取相關(guān)分布規(guī)律，但是并未全面實(shí)現(xiàn)對(duì)機(jī)匣流動(dòng)換熱特性的準(zhǔn)確分析，其特性的構(gòu)建需要進(jìn)行實(shí)驗(yàn)驗(yàn)證[3]。就熱物理模型與計(jì)算獲取邊界條件而言，基于有限元計(jì)算，獲得發(fā)動(dòng)機(jī)運(yùn)轉(zhuǎn)時(shí)機(jī)匣的瞬態(tài)溫度，有助于葉尖間隙分析和機(jī)匣強(qiáng)度設(shè)計(jì)。

總之，想要獲取渦輪機(jī)匣溫度，切實(shí)解決相關(guān)問(wèn)題，關(guān)鍵就在于計(jì)算換熱邊界條件。由于機(jī)匣的復(fù)雜性，導(dǎo)致難以預(yù)測(cè)周?chē)鲃?dòng)換熱環(huán)境。而換熱邊界準(zhǔn)確性在很大程度上影響著機(jī)匣溫度分析的精確性，因此換熱邊界條件分析是機(jī)匣溫度分析的關(guān)鍵，對(duì)此，必須健全既有計(jì)算方法，以數(shù)值仿真與實(shí)驗(yàn)分析保障換熱邊界條件的準(zhǔn)確性。

2 渦輪機(jī)匣數(shù)值仿真分析

高壓渦輪機(jī)匣內(nèi)部的構(gòu)造十分特殊，是一個(gè)腔體，其換熱備受氣流流動(dòng)影響，因?yàn)榻Y(jié)構(gòu)復(fù)雜導(dǎo)致無(wú)法基于既定準(zhǔn)則進(jìn)行計(jì)算。以一維空氣系統(tǒng)流動(dòng)計(jì)算的進(jìn)口壓力與質(zhì)量流量作為邊界條件，通過(guò)CFD 軟件進(jìn)行數(shù)值計(jì)算。渦輪機(jī)匣流場(chǎng)結(jié)構(gòu)在幾何結(jié)構(gòu)影響下，其腔體結(jié)構(gòu)強(qiáng)制內(nèi)部氣流改變方向，以構(gòu)成氣流渦，對(duì)壁面產(chǎn)生相應(yīng)作用，對(duì)氣流和壁面之間的換熱造成直接性影響。通過(guò)進(jìn)氣孔流入機(jī)匣內(nèi)部的氣流速率非常高，直接沖擊著壁面，在沖擊區(qū)的周?chē)刂吘壷饾u流動(dòng)，速率也會(huì)隨之漸漸降低[4]。沖擊孔與渦輪機(jī)匣的外壁面區(qū)域相對(duì)應(yīng)，存在十分顯著的斑狀沖擊強(qiáng)換熱區(qū)域，其特性與沖擊孔相同，是在劇烈沖擊下導(dǎo)致的，也是受空氣冷卻機(jī)匣的作用造成的。

在機(jī)匣的內(nèi)壁上，前部分，受入流孔沖擊，換熱系數(shù)較大，但是沖擊區(qū)域以外，換熱系數(shù)會(huì)不斷縮小，但是減小幅度不斷變?nèi)酢Ｖ虚g部分，氣流和壁面之間的相互作用處于均衡狀態(tài)下，換熱系數(shù)的分布也相應(yīng)的非常均衡，但是比較小。后部分，沖擊位置的氣流和壁面之間的作用十分強(qiáng)烈，換熱比較劇烈，系數(shù)也很大，但是沖擊區(qū)域外，相互作用減小，換熱系數(shù)也比較小。就數(shù)值分析結(jié)果可知，腔體內(nèi)部的氣流速度始終處于平衡狀態(tài)下，壁面換熱也比較均衡[5]。總之，孔出流對(duì)避免的沖擊作用很強(qiáng)，兩者之間的換熱也比較劇烈，系數(shù)較大。渦輪機(jī)匣結(jié)構(gòu)多元化，使得氣流方向逐漸變化，造成漩渦結(jié)構(gòu)復(fù)雜化。氣流渦所在位置上，氣流和壁面的相互作用比較強(qiáng)烈，換熱系數(shù)也較大。

3 渦輪機(jī)匣換熱實(shí)驗(yàn)

3.1 裝 置

為了降低實(shí)驗(yàn)難度，在不對(duì)機(jī)匣流動(dòng)換熱特性造成影響的基礎(chǔ)上，對(duì)既有結(jié)構(gòu)進(jìn)行簡(jiǎn)化，其中機(jī)匣屬于環(huán)腔結(jié)構(gòu)，受對(duì)稱(chēng)性的影響，選取環(huán)狀機(jī)匣具有周期性的1/27 進(jìn)行深入探究，并簡(jiǎn)化其內(nèi)部通道的圓弧面為平面。具體實(shí)驗(yàn)裝置如圖1所示。

圖1 實(shí)驗(yàn)裝置圖

3.2 原 理

基于瞬態(tài)法進(jìn)行機(jī)匣傳熱系數(shù)測(cè)量，通過(guò)加熱器進(jìn)行來(lái)流加熱，物體與環(huán)境、以及熱量傳遞尚未處于穩(wěn)定狀態(tài)的時(shí)候，應(yīng)就實(shí)驗(yàn)可以獲取的溫度數(shù)據(jù)與時(shí)間等相關(guān)參數(shù)，求解傳熱系數(shù)。所以，應(yīng)充分了解溫度數(shù)據(jù)、時(shí)間、傳熱系數(shù)之間的關(guān)系。因?yàn)閷?shí)驗(yàn)的時(shí)間比較短，加熱氣流和換熱之間導(dǎo)致的實(shí)驗(yàn)?zāi)Ｐ蜏囟入S之發(fā)生了一定變化，并未受到任何的其他影響，這時(shí)實(shí)驗(yàn)?zāi)Ｐ蛯?dǎo)熱問(wèn)題可以通過(guò)一維半無(wú)限大導(dǎo)熱加以解決。半無(wú)限大物體模型具體如圖2所示，其從0界面逐漸向正x向與其他坐標(biāo)延伸。這種物體并不存在，但在非穩(wěn)態(tài)導(dǎo)熱的初始階段，很有可能會(huì)把實(shí)物看作是半無(wú)限大物體，然后進(jìn)行適當(dāng)處理。在進(jìn)行實(shí)驗(yàn)時(shí)，把加熱器接通電源之后，氣流的溫度會(huì)漸漸上升，此時(shí)擾流柱通道端壁受熱氣流加熱的影響，其材料是有機(jī)玻璃，熱擴(kuò)散率比較低，測(cè)量時(shí)間較短，溫度變化只是停留在表面附近，并沒(méi)有深入到端壁內(nèi)部，所以，可以把端壁當(dāng)作半無(wú)限大物體，導(dǎo)熱過(guò)程為一維，只需要沿厚度傳熱即可[6-7]。

圖2 一維半無(wú)限大模型

在初始階段，模型溫度之中處于均衡的狀態(tài)下，并與氣流、環(huán)境溫度之間等同，然后在氣流溫度上升的時(shí)候，氣流和壁面進(jìn)行對(duì)流換熱，壁面溫度逐漸上升，同時(shí)朝向模型內(nèi)進(jìn)行傳熱。此過(guò)程為

壁面溫度為

式中：t0為模型與氣流初始溫度，℃；tw為壁面溫度，℃；tg為溫度上升之后氣流溫度，℃；t（x，τ）為τ時(shí)模型深度為x的溫度，℃；ρ為模型密度，kg/m3；c為比熱容，kJ/(kg·℃)；λ為導(dǎo)熱系數(shù)，W/(m·℃)。

在實(shí)驗(yàn)中，端壁非穩(wěn)態(tài)導(dǎo)熱邊界條件是第三類(lèi)條件，傳熱系數(shù)是有限值，此時(shí)半無(wú)限大導(dǎo)熱最大時(shí)間明顯要長(zhǎng)于計(jì)算時(shí)間。就有力玻璃的特性，可知最大時(shí)間為195 s。所以，在比最大時(shí)間小的狀態(tài)下，進(jìn)行端壁溫度場(chǎng)測(cè)試，可以利用半無(wú)限大模型。

3.3 測(cè) 量

雷諾數(shù)為

式中：為通過(guò)實(shí)驗(yàn)段質(zhì)量流量，kg/s；μ為動(dòng)力黏度，N·s/m2；Ain代表入口射流孔面積，m2；d為孔徑，mm。在實(shí)驗(yàn)中，雷諾數(shù)為：12×104～20×104。努賽爾數(shù)為

式中：λg為氣體導(dǎo)熱系數(shù)，W/(m·℃)；h為傳熱系數(shù)，W/(m·℃)。

在進(jìn)行實(shí)驗(yàn)的時(shí)候，電加熱器接通電源之后，氣流開(kāi)始加熱，溫度逐漸上升，此時(shí)應(yīng)開(kāi)啟旁通電磁閥，并將直通電磁閥及時(shí)關(guān)閉。通過(guò)閥門(mén)進(jìn)行工況調(diào)節(jié)，從而保證其穩(wěn)定性和可靠性。在滿足工況要求后，及時(shí)將電磁閥切換掉，加熱氣流基于實(shí)驗(yàn)件，攝像機(jī)也開(kāi)始運(yùn)行，詳細(xì)記錄熱帶偶對(duì)氣流具體溫度的測(cè)量。在結(jié)束并變色之后，實(shí)驗(yàn)結(jié)束[8]。攝像機(jī)全程記錄端壁壁面溫度變化中所產(chǎn)生的顏色變化，以及熱電偶在氣流溫度測(cè)量中產(chǎn)生的變化。在進(jìn)行實(shí)驗(yàn)時(shí)，加熱氣流基于實(shí)驗(yàn)件，壁面受熱氣流加熱，測(cè)量時(shí)間保持幾十秒，可以將端壁作為半無(wú)限大物體。測(cè)量端壁壁面溫度、加熱時(shí)間、氣流溫度等在時(shí)間推移下呈現(xiàn)上升狀態(tài)，結(jié)合模型物性參數(shù)與初始溫度，進(jìn)行傳熱系數(shù)計(jì)算。

3.4 結(jié)果分析

3.4.1 換熱分布

就實(shí)驗(yàn)件結(jié)構(gòu)特性，傳熱系數(shù)測(cè)量主要分為三大區(qū)域，如圖3所示，面a、面b、面c。

圖3 換熱測(cè)量表面

面a位于腔體前部分，氣流通過(guò)入口射流孔緩緩進(jìn)入，對(duì)壁面造成一定沖擊，在受到阻礙之后，及時(shí)改變方向，然后漸漸向下形成回流，流動(dòng)在射流沖擊的限制下進(jìn)行。就努賽爾數(shù)分布可知，面a對(duì)應(yīng)入口所構(gòu)成的高換熱區(qū)域，按照從中間向外的順序，換熱呈現(xiàn)減弱趨勢(shì)。距離沖擊區(qū)域較遠(yuǎn)的傳熱系數(shù)較小，尤其是最左端局部壁面換熱系數(shù)非常小。面b對(duì)應(yīng)腔體后部分，在面b左端所對(duì)應(yīng)的沖擊孔傳熱系數(shù)最大。在墻體凹槽位置的換熱速度下降非常快，這主要是因?yàn)闅饬飨掠瘟鲃?dòng)到凹槽位置的時(shí)候，會(huì)隨之形成順時(shí)針的漩渦，速度快速下降，傳熱系數(shù)也隨之降低。而在入口窄縫的傳熱系數(shù)卻有了一定程度的升高，因此經(jīng)過(guò)窄縫的時(shí)候，流動(dòng)速度提高，傳熱系數(shù)提高。氣流進(jìn)入腔體之后，換熱處于均衡狀態(tài)。面c對(duì)應(yīng)的射流孔沖擊區(qū)域換熱比較強(qiáng)，其中最右側(cè)最弱。面c換熱分布受熱口射流沖擊影響，從射流孔進(jìn)入沖擊斜壁面，周?chē)鷵Q熱顯著較高。但是不同沖擊區(qū)域換熱存在明顯差異，其中中間區(qū)域在沖擊之后反卷的氣流發(fā)生碰撞，能夠強(qiáng)化換熱，因此，此區(qū)域換熱明顯比兩端區(qū)域要高[9]。但是，最中間區(qū)域換熱與相鄰區(qū)域換熱相比卻較低，這是因?yàn)槿肟诹魉佥^小。

3.4.2 平均換熱

在雷諾數(shù)增加的趨勢(shì)下，各面平均努賽爾數(shù)都有所增加，其中，面a換熱>面b換熱>面c換熱。整面努賽爾數(shù)在雷諾數(shù)增加的趨勢(shì)下呈現(xiàn)線性增加的趨勢(shì)。而且在雷諾數(shù)增加的影響下，面c 的努賽爾數(shù)在雷諾數(shù)為200 000 工況的時(shí)候，突然增大。這是由于在此工況下，沖擊孔所對(duì)應(yīng)高換熱區(qū)相鄰區(qū)域內(nèi)不存在明顯的低換熱區(qū)，沖擊孔與中間孔對(duì)應(yīng)的高換熱區(qū)域幾乎呈現(xiàn)相連接的狀態(tài)[10]。

4 結(jié) 論

總之，基于詳細(xì)分析航空發(fā)動(dòng)機(jī)渦輪機(jī)匣流動(dòng)換熱與溫度，掌握機(jī)匣流動(dòng)換熱規(guī)律，構(gòu)建渦輪機(jī)匣傳熱分析技術(shù)體系。通過(guò)進(jìn)行流動(dòng)換熱CFD 數(shù)值仿真與換熱實(shí)驗(yàn)，獲得換熱分布特征。而且，沖擊孔出流強(qiáng)烈沖擊了壁面，在沖擊區(qū)域換熱很高，部分系數(shù)也比較大。由于渦輪機(jī)匣的內(nèi)部結(jié)構(gòu)十分復(fù)雜，氣流方向一直在變化，從而構(gòu)成了繁雜的漩渦結(jié)構(gòu)。在存在氣流渦時(shí)，氣流和壁面之間的相互作用比較強(qiáng)大，壁面部分區(qū)域的換熱系數(shù)很大。另外，在雷諾數(shù)不斷增加的趨勢(shì)下，壁面的平均換熱系數(shù)也按照一定順序在逐漸增大。通過(guò)以上研究成果，不僅可以準(zhǔn)確分析機(jī)匣溫度，還能夠進(jìn)一步健全優(yōu)化傳熱設(shè)計(jì)。