旋轉(zhuǎn)狀態(tài)平板收斂縫形氣膜孔改善冷卻效果的數(shù)值研究

姚 玉，韓 星，盛志強(qiáng)，吳逸飛

(南昌航空大學(xué) 飛行器工程學(xué)院，南昌 330063)

在航空燃?xì)鉁u輪發(fā)動機(jī)的發(fā)展過程中，為提高發(fā)動機(jī)的熱效率和推重比，壓氣機(jī)的增壓比和渦輪進(jìn)口溫度不斷被提高，導(dǎo)致渦輪部件的工作環(huán)境溫度已遠(yuǎn)超材料的耐熱極限。氣膜冷卻因固有的高效冷卻特性已被廣泛地應(yīng)用于渦輪葉片上，成為航空燃?xì)鉁u輪發(fā)動機(jī)高溫部件的主要冷卻措施之一。氣膜的噴入，雖然減少了高溫燃?xì)馀c葉片表面之間的熱量傳遞，但要實(shí)現(xiàn)氣膜冷卻，就得抽取一部分壓縮空氣，旁路過燃燒室直接冷卻熱端部件，會引起熱力學(xué)損失。因此，盡管國內(nèi)外針對氣膜冷卻的流動與換熱特性開展了大量的研究工作，但是如何進(jìn)一步提高氣膜的冷卻效率仍然是一個富于挑戰(zhàn)和創(chuàng)新的科學(xué)問題[1-5]。

氣膜冷卻的核心問題在于降低冷氣出流向主流的穿透率，避免主流繞過冷氣直接接觸壁面造成局部高溫以及增強(qiáng)氣膜出流向下游的延伸能力，增強(qiáng)冷卻效果等。在影響氣膜冷卻效率的眾多因素中，氣膜孔幾何結(jié)構(gòu)對冷卻效率和流動特性具有重要的影響，優(yōu)化氣膜孔的幾何結(jié)構(gòu)一直是研究的重點(diǎn)，各國學(xué)者在這一領(lǐng)域開展了大量的研究工作[6-11]，目的在于通過改善氣膜孔冷卻結(jié)構(gòu)，實(shí)現(xiàn)使用相對少的冷卻空氣量而獲得相對好的冷卻效果。

在氣膜孔結(jié)構(gòu)優(yōu)化方案中，具有出口擴(kuò)展型面的氣膜孔結(jié)構(gòu)得到了國內(nèi)外學(xué)者的普遍關(guān)注，如扇形孔[12]、帶復(fù)合傾角的扇形孔[13]、箭頭形孔[14]等。冷氣通道在出口附近橫截面面積的增大，導(dǎo)致了氣膜出流平均速度的下降和擴(kuò)散能力的增強(qiáng)，降低了氣膜向主流的穿透率，并增加了氣膜的展向覆蓋區(qū)域，各種擴(kuò)展型面結(jié)構(gòu)均不同程度地抑制了反向渦對的產(chǎn)生，降低了旋渦的強(qiáng)度，提高了射流的貼壁性，從而增強(qiáng)了壁面的冷卻效果，但是單純的擴(kuò)張形孔并沒有將冷氣充分利用，而且會帶來較大的氣動損失。

Sargison[15-16]等人提出了收斂縫形氣膜孔，此結(jié)構(gòu)的冷卻效果與縫形氣膜孔類似，冷卻效率較圓形孔高，并且氣膜出流后下游邊界層的流動情況與無氣膜時類似，大大降低了氣動損失；Azzi[17]等人指出收斂縫形氣膜孔下游形成渦對的旋向與圓形孔方向相反，已有的研究結(jié)果[18-24]表明此渦對可將低溫冷氣由旋渦中心旋出，隔離主流避免高溫氣體與壁面直接接觸，達(dá)到降低壁溫的目的，且渦對的尺寸隨著吹風(fēng)比的增加而增大，能夠隔離主流的冷氣越多，冷卻效率也隨之提高。但對于旋轉(zhuǎn)條件下的收斂縫形氣膜孔研究不夠充分，為此，本文通過數(shù)值模擬的方法，對收斂縫形氣膜孔旋轉(zhuǎn)時的冷卻特性進(jìn)行研究，著重與圓形孔對比。在相同氣動參數(shù)條件下，通過分析下游截面溫度場和局部、平均的冷卻效率，驗(yàn)證收斂縫形氣膜孔在旋轉(zhuǎn)條件下同樣具有改善冷卻效果的能力。

1 計算模型

旋轉(zhuǎn)狀態(tài)下平板收斂縫形氣膜孔模型采用簡化的矩形通道，其幾何結(jié)構(gòu)見圖1。計算模型區(qū)域由兩個子區(qū)域組成：主流區(qū)和氣膜孔區(qū)。主流區(qū)為簡化的矩形通道，氣膜孔區(qū)為圓形孔或收斂縫形氣膜孔。矩形通道左側(cè)為主流入口，二次流以一定入射角度由主流通道下方的傾斜氣膜孔進(jìn)入，主、次流相互摻混后由右側(cè)的壓力出口流出。氣膜孔板被簡化為平板結(jié)構(gòu)，且不考慮孔板壁厚的影響。主流矩形通道沿流動方向長為80 mm，高為12 mm，橫向?yàn)?4 mm；氣膜孔設(shè)置在主流入口下游10d(d為氣膜孔入口圓形直徑)處，計算模型坐標(biāo)原點(diǎn)取為氣膜孔出口中心，流動方向?yàn)閤軸(也為氣膜冷卻結(jié)構(gòu)的旋轉(zhuǎn)軸)，高度為y軸，展向?yàn)閦軸(即離心力的方向)；旋轉(zhuǎn)半徑為100d，正向旋轉(zhuǎn)的情況下，被冷卻面為壓力面。

圖1 計算模型示意圖

收斂縫形氣膜孔是一種進(jìn)口為圓形、出口為矩形的氣膜孔結(jié)構(gòu)(圖2)。在圓轉(zhuǎn)矩形的過渡過程中，氣膜孔沿出流方向逐漸收斂，而在展向則呈兩側(cè)對稱擴(kuò)張狀。由于其收斂的幅度大于擴(kuò)張，從進(jìn)口到出口其橫截面的面積逐漸縮小，因而冷氣從入口到出口逐漸加速，在出口狹縫處冷氣速度達(dá)到最大。氣膜孔的主要結(jié)構(gòu)尺寸為：氣膜孔入口圓形直徑d為2 mm，出口縫寬s為0.2d，出口縫長l為2d，氣膜孔高h(yuǎn)為1.2d，氣膜孔軸線與流動方向的夾角α為30°。

(a) 主視圖

(b) 俯視圖圖2 收斂縫形氣膜孔冷卻模型示意圖

2 計算方法和驗(yàn)證

2.1 邊界條件

計算域的邊界條件設(shè)置：主流和二次流進(jìn)口為流量入口；混合出口為壓力出口；氣膜孔出口為內(nèi)部面；被冷卻壁面采用絕熱壁面，其余壁面同樣采用絕熱壁面條件，并且均采用無滑移壁面條件。

主流入口參數(shù)：雷諾數(shù)Re為6 000，溫度為600 K。氣膜孔入口參數(shù)：根據(jù)吹風(fēng)比和主流雷諾數(shù)對二次流入口流量進(jìn)行調(diào)節(jié)，吹風(fēng)比范圍為0.8～2.4，溫度為300 K?；旌蠚怏w出口參數(shù)：壓力為101 325 Pa，并滿足質(zhì)量守恒條件。冷卻模型轉(zhuǎn)動參數(shù)：轉(zhuǎn)速Ω范圍為0～750 r/min。主、次流的馬赫數(shù)都不高，均選用不可壓縮的理想氣體。

2.2 參數(shù)定義

吹風(fēng)比定義為：

(1)

式中：ρ2、u2分別表示二次流的密度和速度；ρ∞、u∞分別表示主流的密度和速度；m2、m∞分別表示二次流和主流的質(zhì)量流量；A2、A∞分別表示二次流和主流通道的橫截面積。

吹風(fēng)比的定義中，為保證冷氣入口流量的一致，采用與收斂縫形氣膜孔有相同入口形狀的圓形孔作為比較的基準(zhǔn)。

保持主流和二次流的溫度不變，即保持主、次流的密度比一定，通過改變主流與二次流的質(zhì)量流量比來調(diào)節(jié)吹風(fēng)比的變化。

絕熱冷卻效率定義為：

(2)

式中：T∞、T2分別表示主流和二次流的入口溫度；Taw表示為有氣膜冷卻情況下沿下游某處絕熱壁面上的恢復(fù)溫度，也就是絕熱壁面的溫度。

通過在絕熱壁面上進(jìn)行展向積分，可得到平均冷卻效率ηav。

2.3 計算方法驗(yàn)證

為了對計算模型的精確性進(jìn)行驗(yàn)證，選用文獻(xiàn)[20]的實(shí)驗(yàn)結(jié)果作為比較基準(zhǔn)。圖3是采用與參考文獻(xiàn)[20]相同的圓形氣膜冷卻結(jié)構(gòu)和工況時，本文模型計算出沿流向的平均絕熱冷卻效率與實(shí)驗(yàn)結(jié)果的比較。可見，所選湍流模型的計算結(jié)果與實(shí)驗(yàn)相差較小，較為合理。

圖3 模型計算驗(yàn)證

本文采用Fluent分離隱式求解器進(jìn)行穩(wěn)態(tài)求解，選用Realizablek-ε雙方程湍流模型加非平衡的壁面函數(shù)；各物理量的離散格式均為二階迎風(fēng)；壓力-速度耦合采用Simplec算法；解收斂的標(biāo)準(zhǔn)是各項(xiàng)殘差精度均小于10-5、殘差曲線趨于平直且被冷卻壁面的平均溫度變化不大。

2.4 網(wǎng)格劃分

網(wǎng)格的疏密很大程度上影響了計算結(jié)果的準(zhǔn)確性及精確性，故采用五種網(wǎng)格劃分?jǐn)?shù)目進(jìn)行網(wǎng)格獨(dú)立性驗(yàn)證，根據(jù)氣膜孔板壁面第一層網(wǎng)格節(jié)點(diǎn)的距離來調(diào)整網(wǎng)格數(shù)量，驗(yàn)證工況見表1。

表1 網(wǎng)格獨(dú)立性驗(yàn)證工況

通過比較各種網(wǎng)格數(shù)目的壁面平均冷卻效率，來驗(yàn)證網(wǎng)格的獨(dú)立性(見圖4)。可以看出，五種網(wǎng)格數(shù)目的平均冷卻效率都很接近，但在50萬～84萬范圍內(nèi)，隨著網(wǎng)格數(shù)目的增加，平均冷卻效率小幅度降低，當(dāng)超過84萬后幾條曲線基本重合在一起，對冷卻效率的影響不大，此時認(rèn)為進(jìn)入了網(wǎng)格穩(wěn)定區(qū)，因而可認(rèn)為此范圍內(nèi)的網(wǎng)格滿足獨(dú)立性驗(yàn)證條件。綜合考慮計算時間，數(shù)值模擬中最終選擇的網(wǎng)格數(shù)目為84萬。

圖4 網(wǎng)格獨(dú)立性驗(yàn)證

研究的主要目的是獲取被冷卻壁面的冷卻特性，而由于氣膜的注入，使得靠近被冷卻壁面附近的主流的流動發(fā)生了較大的變化，兩股氣流之間會產(chǎn)生較強(qiáng)的剪切作用。為了準(zhǔn)確模擬這一流動過程，劃分了附面層網(wǎng)格。沿著被冷卻壁面的法線方向，第一層網(wǎng)格布置在距壁面0.02 mm的位置，此后每層網(wǎng)格間距以1.2的比率增加，到第15層之后，改為1.03的比率增加。為了節(jié)約網(wǎng)格數(shù)量，以氣膜孔出口狹縫為中心，沿上、下游網(wǎng)格逐漸稀疏，而展向上網(wǎng)格均勻分布。整個主流區(qū)域采用非均勻結(jié)構(gòu)化網(wǎng)格，共有網(wǎng)格約84萬，網(wǎng)格分布見圖5。

圖5 網(wǎng)格分布示意圖

3 結(jié)果與分析

圖6是吹風(fēng)比M=2.0時，收斂縫形氣膜孔和圓形氣膜孔在下游截面x/d=2處溫度的對比。可以看出，轉(zhuǎn)速Ω對兩種孔型氣膜的分布均有影響，都會朝著離心方向發(fā)生偏移。隨著轉(zhuǎn)速的增加,兩種孔型對主流的穿透作用均有增強(qiáng)，也就是冷氣與主流的摻混增多，這對冷卻是不利的；但是收斂縫形氣膜孔的穿透率依然要低于圓形氣膜孔，意味著旋轉(zhuǎn)條件下，收斂縫形氣膜孔的貼壁性依然較好，可以較好地形成氣膜。由于收斂縫形氣膜孔出口在展向范圍的擴(kuò)大，使得冷氣核心的溫度稍高于圓形氣膜孔，降低了冷卻能力，但圓形氣膜孔的冷氣核心未緊貼壁面，并沒有形成較好的冷卻效果。同時在旋轉(zhuǎn)的影響下，收斂縫形氣膜孔在展向冷氣覆蓋的區(qū)域得到了增大，也就是擴(kuò)大了冷卻的有效區(qū)，因而展向的冷卻有效范圍要優(yōu)于圓形孔。大吹風(fēng)比工況條件下，x/d=2截面處圓形孔的冷氣已經(jīng)出現(xiàn)脫離壁面的趨勢，而收斂縫形氣膜孔則始終貼覆壁面形成較好的保護(hù)氣膜。

圖6 M=2.0氣膜孔下游x/d=2截面的溫度分布

圖7是吹風(fēng)比M=2.0時，收斂縫形氣膜孔和圓形氣膜孔在下游截面x/d=5處溫度的對比。在x/d=5截面處，圓形氣膜孔同樣對主流的穿透高于收斂縫形氣膜孔，并且圓形氣膜孔冷氣核心的溫度與收斂縫形氣膜孔已相差無幾，說明局部冷卻效果較好的優(yōu)勢已經(jīng)不復(fù)存在，而收斂縫形氣膜孔展向范圍的冷氣覆蓋依然好于圓形氣膜孔。大吹風(fēng)比的工況下，雖然高轉(zhuǎn)速條件下圓形氣膜孔對主流的穿透得到抑制，但依然出現(xiàn)脫離壁面的現(xiàn)象，而收斂縫形氣膜孔的氣膜覆蓋依然較好，因而此截面處收斂縫形氣膜孔的冷卻效果明顯較優(yōu)。

圖7 M=2.0氣膜孔下游x/d=5截面的溫度分布

圖8是收斂縫形氣膜孔(console)和圓形氣膜孔(cylindrical)下游各截面展向冷卻效率的對比。由圖可以看出，旋轉(zhuǎn)條件下，與圓形氣膜孔在小吹風(fēng)比M=0.8工況相比，收斂縫形氣膜孔的優(yōu)勢并不明顯，甚至在氣膜孔中心線附近靠近氣膜孔出口區(qū)域，其冷卻效果低于圓形氣膜孔，僅在展向范圍稍優(yōu)于圓形氣膜孔；而在下游遠(yuǎn)處，兩種孔型的冷卻效果相差不大，收斂縫形氣膜孔僅有些許優(yōu)勢。雖然吹風(fēng)比M=1.2時，在x/d=2截面處有較小區(qū)域的冷卻效率低于圓形氣膜孔，但x/d=2截面處其余區(qū)域以及其他下游截面的冷卻效果無論是均勻性還是效率均好于圓形孔。隨著吹風(fēng)比的進(jìn)一步增加，收斂縫形氣膜孔的優(yōu)勢更加明顯，在冷卻范圍、冷卻效率以及冷卻均勻性等方面均高于圓形氣膜孔，特別是當(dāng)圓形氣膜孔出現(xiàn)脫離壁面的情況時，收斂縫形氣膜孔由于其固有的通道特點(diǎn)使冷氣依然緊貼壁面，因而其冷卻效果大大好于圓形氣膜孔。

圖8 收斂縫形氣膜孔與圓形氣膜孔下游截面展向冷卻效率的對比

圖9是轉(zhuǎn)速為500 r/min時，圓形氣膜孔和收斂縫形氣膜孔冷卻效率的對比分布圖?？梢钥闯?，旋轉(zhuǎn)條件下，氣膜均朝著離心的方向偏移，并且隨著吹風(fēng)比的增大，偏移越明顯。由于收斂縫形氣膜孔的出口在展向的擴(kuò)張以及氣膜出流對主流的穿透率較低，故旋轉(zhuǎn)條件下，雖然氣膜存在偏移，但依然不能改變收斂縫形氣膜孔優(yōu)異性的表現(xiàn)，無論是在展向上氣膜覆蓋的均勻度，還是沿流向上氣膜的延展度，收斂縫形氣膜孔的表現(xiàn)都要優(yōu)于圓形氣膜孔。特別是大吹風(fēng)比的工況，圓形氣膜孔的氣膜呈現(xiàn)出脫離壁面的特點(diǎn)，而收斂縫形氣膜孔則無此現(xiàn)象，反而由于冷氣流量的增多，可以更好地形成氣膜保護(hù)壁面，增強(qiáng)冷卻效果。因此，旋轉(zhuǎn)條件下，收斂縫形氣膜孔的冷卻效果依然要優(yōu)于圓形氣膜孔，起到改善冷卻效果的作用。

圖9 Ω=500 r/min 收斂縫形氣膜孔與圓形氣膜孔冷卻效率的對比

圖10是不同轉(zhuǎn)速條件下，收斂縫形氣膜孔(console)和圓形氣膜孔(cylindrical)平均冷卻效率的對比。由于平均冷卻效率曲線考慮了氣膜孔出口處的冷氣，而收斂縫形氣膜孔出口處的縫寬較小，因而氣膜孔出口區(qū)域的平均冷卻效率出現(xiàn)圓形氣膜孔高于收斂縫形氣膜孔的現(xiàn)象。而在氣膜孔出口下游x/d≥1的區(qū)域，無論何種工況，收斂縫形氣膜孔的冷卻效率均高于圓形氣膜孔，特別是大吹風(fēng)比的情況，較好的貼壁性使得收斂縫形氣膜孔的冷卻效率遠(yuǎn)高于圓形氣膜孔，呈現(xiàn)出優(yōu)異的冷卻效果。轉(zhuǎn)速在500～750 r/min范圍時，雖然轉(zhuǎn)速的增大提高了圓形氣膜孔的氣膜在高吹風(fēng)比時的貼壁性，增強(qiáng)了冷卻效果，但同時轉(zhuǎn)速對收斂縫形氣膜孔也產(chǎn)生了較大的影響，在x/d<15范圍內(nèi)冷卻效率也得到了明顯的提高，因此，高轉(zhuǎn)速高吹風(fēng)比情況下，收斂縫形氣膜孔依然保持較優(yōu)的冷卻效果。

圖10 收斂縫形氣膜孔與圓形氣膜孔平均冷卻效率的對比

4 結(jié)論

本文針對平板單個收斂縫形氣膜孔的冷卻結(jié)構(gòu)，通過數(shù)值研究與圓形氣膜孔進(jìn)行對比，說明它在旋轉(zhuǎn)條件下，同樣具備改善氣膜冷卻效率的特性?？傻玫饺缦陆Y(jié)論：

(1) 旋轉(zhuǎn)條件下，氣膜均朝著離心的方向發(fā)生偏移，并且隨著吹風(fēng)比的增大，偏移越明顯，偏移使得展向冷氣覆蓋的區(qū)域增大，擴(kuò)大了冷卻的有效區(qū)，特別是收斂縫形氣膜孔，展向的冷卻有效范圍明顯優(yōu)于圓形氣膜孔。

(2) 旋轉(zhuǎn)條件下，兩種孔型對主流的穿透作用均有增強(qiáng)，但是收斂縫形氣膜孔的穿透率依然要低于圓形氣膜孔，說明旋轉(zhuǎn)工況下，收斂縫形氣膜孔的貼壁性依然較好，可以更好地形成氣膜。

(3) 考慮氣膜孔出口處冷氣的情況下，雖然局部區(qū)域出現(xiàn)了圓形氣膜孔的平均冷卻效率高于收斂縫形氣膜孔的現(xiàn)象，但在氣膜孔出口下游x/d≥1的區(qū)域，無論何種工況收斂縫形氣膜孔的冷卻效率均高于圓形氣膜孔，特別是大吹風(fēng)比的工況，較好的貼壁性使得收斂縫形氣膜孔的冷卻效率遠(yuǎn)高于圓形氣膜孔，呈現(xiàn)出優(yōu)異的冷卻效果。

(4) 隨著吹風(fēng)比的增加，收斂縫形氣膜孔的優(yōu)勢更加明顯，無論是在展向上氣膜覆蓋的均勻度，還是沿流向上氣膜的延展度，收斂縫形氣膜孔的表現(xiàn)都要優(yōu)于圓形氣膜孔。

(5) 各種轉(zhuǎn)速條件下，大吹風(fēng)比的工況，圓形孔的氣膜呈現(xiàn)出脫離壁面的特點(diǎn)，而收斂縫形氣膜孔則無此現(xiàn)象，反而由于冷氣流量的增多，可以更好的形成氣膜保護(hù)壁面，增強(qiáng)冷卻效果，因而其冷卻效果好于圓形氣膜孔。