雙輻板渦輪盤盤腔流動(dòng)換熱分析

韓玉琪，賈志剛，劉 紅，朱大明

(中國(guó)航空發(fā)動(dòng)機(jī)研究院，北京101304)

雙輻板渦輪盤盤腔流動(dòng)換熱分析

韓玉琪，賈志剛，劉 紅，朱大明

(中國(guó)航空發(fā)動(dòng)機(jī)研究院，北京101304)

基于UG軟件建立了雙輻板渦輪盤的三維模型，基于CFX軟件對(duì)雙輻板渦輪盤盤腔內(nèi)的流動(dòng)換熱現(xiàn)象進(jìn)行了數(shù)值模擬，獲得了典型工況下旋轉(zhuǎn)盤腔內(nèi)的流場(chǎng)和溫度場(chǎng)分布，初步總結(jié)了該盤腔結(jié)構(gòu)內(nèi)的流動(dòng)換熱規(guī)律。結(jié)果表明：該盤腔結(jié)構(gòu)內(nèi)的流動(dòng)形式為典型的源-匯流動(dòng)，冷氣的總壓和溫度沿徑向均逐漸增大，總壓增和溫度差隨轉(zhuǎn)速的增加而增大、隨冷氣流量的增加而減小，對(duì)流換熱系數(shù)在不同流場(chǎng)區(qū)域表現(xiàn)出不同的分布特征。為雙輻板渦輪盤的進(jìn)一步結(jié)構(gòu)優(yōu)化設(shè)計(jì)提供了參考。

航空發(fā)動(dòng)機(jī)；雙輻板渦輪盤；旋轉(zhuǎn)盤腔；冷卻；流動(dòng)；換熱；數(shù)值模擬

numericalsimulation

1 引言

據(jù)統(tǒng)計(jì)，在航空發(fā)動(dòng)機(jī)所有非包容事故中，大約一半是由輪盤損壞引起。渦輪盤工作過程中承受離心力載荷、熱應(yīng)力載荷和輪盤外載荷，是典型的壽命限制件[1-2]。對(duì)于渦輪發(fā)動(dòng)機(jī)，高推重比的目標(biāo)要求渦輪前溫度越來越高，在盤體材料耐溫極限增長(zhǎng)極其有限的現(xiàn)狀下，渦輪部件的高效冷卻成為航空發(fā)動(dòng)機(jī)的核心技術(shù)之一。雙輻板渦輪盤是針對(duì)下一代高推重比渦扇發(fā)動(dòng)機(jī)的高壓渦輪盤設(shè)計(jì)，在減輕質(zhì)量、提高冷卻效率方面具有很大的發(fā)展前景。由于雙輻板渦輪盤采用了空心結(jié)構(gòu)，冷卻空氣可以直接進(jìn)入兩個(gè)輻板形成的空腔進(jìn)行冷卻，與傳統(tǒng)的實(shí)心渦輪盤相比，其散熱面積增大、冷卻效果更好。

Cairo[3-4]申請(qǐng)了雙輻板渦輪盤結(jié)構(gòu)的專利，并對(duì)該雙輻板渦輪盤進(jìn)行了分析，實(shí)施了不同程度的模擬試驗(yàn)，關(guān)注了焊接區(qū)域的應(yīng)力和疲勞壽命。為防止兩個(gè)輻板在高轉(zhuǎn)速工況下盤心由于軸向變形而接觸，Harding[5]設(shè)計(jì)了一個(gè)盤轂尺寸較大的雙輻板渦輪盤。Burge[6]申請(qǐng)了應(yīng)用于高壓壓氣機(jī)最后一級(jí)輪盤的雙輻板輪盤專利。江和甫[7]指出，雙輻板渦輪盤的技術(shù)關(guān)鍵為焊接工藝、冷卻設(shè)計(jì)、強(qiáng)度結(jié)構(gòu)設(shè)計(jì)。陸山[8]提出并建立了雙輻板渦輪盤結(jié)構(gòu)優(yōu)化設(shè)計(jì)數(shù)學(xué)模型及方法，篩選了子午面形狀設(shè)計(jì)參數(shù)，并針對(duì)典型高負(fù)荷渦輪盤結(jié)構(gòu)優(yōu)化問題進(jìn)行了雙輻板盤結(jié)構(gòu)優(yōu)化設(shè)計(jì)。欒永先[9-10]對(duì)雙輻板渦輪盤的結(jié)構(gòu)特點(diǎn)和工藝難點(diǎn)進(jìn)行了介紹，通過與傳統(tǒng)渦輪盤進(jìn)行對(duì)比闡述了雙輻板結(jié)構(gòu)的先進(jìn)性，并應(yīng)用有限元分析軟件對(duì)傳統(tǒng)渦輪盤和雙輻板渦輪盤進(jìn)行了強(qiáng)度分析，提出其未來的研制設(shè)想。陸山[11]采用盤/榫結(jié)構(gòu)分部快速優(yōu)化/整體精細(xì)優(yōu)化的高效優(yōu)化策略，進(jìn)行了雙輻板渦輪盤/榫的結(jié)構(gòu)優(yōu)化設(shè)計(jì)。董少靜[12]采用漸進(jìn)結(jié)構(gòu)優(yōu)化算法確定出雙輻板渦輪盤的結(jié)構(gòu)形式，并對(duì)其進(jìn)行了有限元分析和尺寸優(yōu)化，使得同等應(yīng)力水平下的雙輻板渦輪盤比傳統(tǒng)渦輪盤的質(zhì)量下降了23.6%，還通過三維旋轉(zhuǎn)光彈試驗(yàn)驗(yàn)證了所提出的雙輻板渦輪盤結(jié)構(gòu)的合理性和相關(guān)計(jì)算的正確性。丁水汀[13]提出，對(duì)于渦輪盤腔的冷卻問題，可以采用阻力評(píng)價(jià)、換熱評(píng)價(jià)和溫度分布評(píng)價(jià)三種評(píng)價(jià)指標(biāo)構(gòu)成工程評(píng)價(jià)體系，前者決定了全系統(tǒng)的阻力損失的大小，后兩者表征了相對(duì)冷卻效果和熱應(yīng)力水平的高低。趙熙[14]的研究表明，在與實(shí)心盤等質(zhì)量的空心盤結(jié)構(gòu)中，中心進(jìn)氣轉(zhuǎn)靜系+高位進(jìn)氣旋轉(zhuǎn)系、高位進(jìn)氣轉(zhuǎn)靜系+高位進(jìn)氣旋轉(zhuǎn)系結(jié)構(gòu)具有較好的換熱特性。金琰[15]通過劃分各個(gè)區(qū)域的換熱模型，計(jì)算了雙輻板渦輪盤的溫度分布，并通過間接熱-結(jié)構(gòu)耦合計(jì)算了盤體上的應(yīng)力分布。

目前，在雙輻板渦輪盤領(lǐng)域公開的文獻(xiàn)很少，且集中于盤體強(qiáng)度分析和結(jié)構(gòu)優(yōu)化方面，所采用的模型簡(jiǎn)化程度大都很高，本文基于更接近工程實(shí)際的三維雙輻板渦輪盤的模型開展了冷卻效果分析。以雙輻板渦輪盤盤腔為研究對(duì)象，采用CFX軟件對(duì)盤腔內(nèi)的流動(dòng)換熱現(xiàn)象進(jìn)行數(shù)值模擬，研究該結(jié)構(gòu)盤腔中壓力場(chǎng)、流場(chǎng)和溫度場(chǎng)的分布，分析總結(jié)其流動(dòng)換熱規(guī)律，以期為雙輻板渦輪盤盤腔冷卻設(shè)計(jì)提供技術(shù)支撐。

2 物理模型與數(shù)值方法

本文分析的雙輻板渦輪盤模型如圖1所示，沿盤腔周向分布了78個(gè)出氣孔。

圖1 雙輻板渦輪盤模型Fig.1 Modelof twin-web rotor disk

為減小計(jì)算量，利用盤腔的周期對(duì)稱性，取全盤的1/78進(jìn)行網(wǎng)格劃分。網(wǎng)格劃分由商業(yè)軟件ICEM完成。為提高計(jì)算精度，流體區(qū)域和固體區(qū)域均采用六面體網(wǎng)格劃分方法。為保證數(shù)值模擬精度，對(duì)邊界層網(wǎng)格進(jìn)行了加密處理，并保證壁面第一層網(wǎng)格處y+＜2，邊界層厚度增長(zhǎng)率為1.2，以滿足SST模型的要求。經(jīng)網(wǎng)格無關(guān)解驗(yàn)證后，流體域包含58萬個(gè)網(wǎng)格單元，固體域包含54萬個(gè)網(wǎng)格單元，共計(jì)112萬個(gè)網(wǎng)格單元。

使用CFX軟件對(duì)雙輻板渦輪盤的流動(dòng)換熱問題進(jìn)行穩(wěn)態(tài)分析。流體控制方程N(yùn)-S方程組的離散格式為高精度格式，湍流模型選為SST兩方程模型，離散格式為高精度格式。流體入口給定質(zhì)量流量和溫度T=800 K，出口給定壓力p=1.5MPa，流固交界面上使用無滑移邊界條件，周期面上使用周期對(duì)稱邊界條件；固體盤緣給定熱流密度q=2 000W/m2，周期面上使用周期對(duì)稱邊界條件，流固交界面耦合求解，其余固體壁面使用絕熱條件，收斂條件為最大殘差小于10-5。數(shù)值模擬中使用的氣體為理想氣體，定性溫度取為入口氣流溫度(即800 K)，參考?jí)毫θ?.5 MPa。固體區(qū)域材料選擇鋼。分析了6種不同冷氣流量和3種不同轉(zhuǎn)速下盤腔的流動(dòng)換熱工況，共計(jì)18種工況。6種入口冷氣流量(m)分別為0.078、0.117、0.156、0.195、0.234、0.273 kg/s，對(duì)應(yīng)的無量綱流量系數(shù)Cw在104量級(jí)；3種輪盤轉(zhuǎn)速(ω)分別為6 000、8 000、10 000 r/min，對(duì)應(yīng)的無量綱旋轉(zhuǎn)雷諾數(shù)Rew在106至107量級(jí)。

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

3.1 流動(dòng)結(jié)構(gòu)

文中計(jì)算的各種工況下流動(dòng)結(jié)構(gòu)相似，下面以冷氣流量0.195 kg/s、盤腔轉(zhuǎn)速10 000 r/min工況為例進(jìn)行分析。圖2為對(duì)稱面上的流線分布圖，圖中清晰地展示了徑向內(nèi)流旋轉(zhuǎn)盤腔的流動(dòng)結(jié)構(gòu)。流體進(jìn)入盤腔后，首先是源區(qū)，隨后逐漸被卷吸進(jìn)入雙輻板內(nèi)側(cè)的埃克曼層內(nèi)，中間為核區(qū)，最后通過匯區(qū)流入出氣孔，呈現(xiàn)出典型的源-匯流動(dòng)形式。

為對(duì)比不同半徑處流體切向速度與當(dāng)?shù)匦D(zhuǎn)速度的相對(duì)大小，圖3給出了不同轉(zhuǎn)速下中心線上旋流系數(shù)β=w/ωr的分布曲線。其中，w為當(dāng)?shù)貧怏w的周向速度，r為該點(diǎn)的徑向坐標(biāo)，R為輪盤外半徑。從圖中可以看到，各種工況下旋流系數(shù)均處于0到1之間，即流體的切向速度小于當(dāng)?shù)乇P腔的旋轉(zhuǎn)速度。流體進(jìn)入盤腔后剛開始旋流系數(shù)為0，此時(shí)位于源區(qū)；隨后逐漸被卷吸進(jìn)入兩側(cè)的埃克曼層，此時(shí)中心核區(qū)流體的旋流系數(shù)逐漸增大，即圖中曲線0到1的增長(zhǎng)階段；然后通過匯區(qū)進(jìn)入出氣孔(圖中各曲線的最后兩個(gè)數(shù)據(jù)點(diǎn)位于出氣孔內(nèi))，旋流系數(shù)達(dá)到1。觀察各曲線中旋流系數(shù)為0的初始段可以發(fā)現(xiàn)，在相同轉(zhuǎn)速下，隨著冷氣流量的不斷加大，冷氣的慣性力逐漸變大，源區(qū)的長(zhǎng)度尺度也逐漸增大，核心區(qū)的旋流系數(shù)變小；在相同冷氣流量下，隨著盤腔旋轉(zhuǎn)速度的增加，源區(qū)的長(zhǎng)度尺度逐漸減小，核心區(qū)的旋流系數(shù)逐漸增大。

圖2 冷氣流量0.195 kg/s、盤腔轉(zhuǎn)速10 000 r/m in工況下對(duì)稱面上的流線分布Fig.2 Stream line distribution on symetric plane at condition m=0.195 kg/s,ω=10 000 r/min

3.2 總壓損失

圖3 不同轉(zhuǎn)速下中心線上的旋流系數(shù)分布Fig.3 Distribution of relative rotational speedβon centerline atdifferent rotational speeds

3.3 溫度分布

以冷氣流量0.195 kg/s、盤腔轉(zhuǎn)速10 000 r/min工況為例進(jìn)行分析。圖5為該工況下對(duì)稱面上的溫度分布云圖，可見沿半徑增大方向，流體域和固體域的溫度都逐漸升高，溫度梯度也逐漸變大，出氣孔內(nèi)流體與固體的接觸面積小、換熱不充分，相應(yīng)部位的固體域的溫度梯度較大。圖6示出了不同工況下盤體的最高溫度，可見：相同轉(zhuǎn)速下，隨著冷氣流量的增大，盤體最高溫度逐漸降低；相同冷氣流量下，隨著旋轉(zhuǎn)速度的增大，盤體的最高溫度逐漸升高。

圖4 不同工況下的總壓損失分布Fig.4 Distribution linesof the total pressure lossat different conditions

圖5 冷氣流量0.195 kg/s、盤腔轉(zhuǎn)速10 000 r/min工況下的溫度分布Fig.5 Distribution of temperature atcondition m=0.195 kg/s,ω=10 000 r/min

圖6 不同工況下盤體的最高溫度Fig.6 Themax temperature on disk atdifferent conditions

3.4 對(duì)流換熱系數(shù)分布

圖7為不同轉(zhuǎn)速下對(duì)流換熱系數(shù)h沿徑向的分布。可見，相同盤腔轉(zhuǎn)速下，初始的源區(qū)內(nèi)的對(duì)流換熱系數(shù)較小，且沿半徑增大方向略有減小；然后在埃克曼層區(qū)域內(nèi)，對(duì)流換熱系數(shù)隨半徑的增大而增大；最后在出氣孔內(nèi)，對(duì)流換熱系數(shù)隨半徑的增大而減小(各曲線最后兩個(gè)數(shù)據(jù)點(diǎn)位于出氣孔內(nèi))。

圖7 不同轉(zhuǎn)速下的對(duì)流換熱系數(shù)分布Fig.7 Distribution ofheat transfer coefficientatdifferent rotationalspeeds

4 結(jié)論

詳細(xì)研究了雙輻板渦輪盤腔內(nèi)的流動(dòng)和換熱現(xiàn)象，分析了該盤腔結(jié)構(gòu)內(nèi)流動(dòng)結(jié)構(gòu)、損失特性、溫度變化和對(duì)流換熱系數(shù)變化，主要結(jié)論如下：

(1)本文中的雙輻板盤腔內(nèi)的流動(dòng)是典型的源-匯流動(dòng)，盤腔中心線上的旋流系數(shù)隨冷氣流量的增大而減小，隨旋轉(zhuǎn)速度的增大而增大。

(2)受離心增壓作用，流體出口總壓大于進(jìn)口總壓，總壓差隨冷氣流量的增大而減小，隨旋轉(zhuǎn)速度的增大而增大。

(3)沿半徑增大方向，盤體溫度逐漸升高，在盤緣處達(dá)到最高；盤體最高溫度隨冷氣流量的增大而減小，隨旋轉(zhuǎn)速度的增大而增大。

(4)沿半徑增大方向，對(duì)流換熱系數(shù)在源區(qū)內(nèi)逐漸減小，在埃克曼層內(nèi)逐漸增大，進(jìn)入出氣孔后逐漸減小。

[1]FAA.Advisory Circular 33.75-1A:Guidancematerial for 14 CFR 33.75,safety analysis[S].

[2]Owen J M，Wilson M.Some current research in rotat?ing-disc systems[C]//.Turbine 2000 International Sympo?sium on Heat Transfer in Gas Turbine Systems.2001.

[3]Cairo R R.Composite ring reinforced turbine program fi?nal report[R].AFRL-PR-WP-TR-1999-2050，1999.

[4]Cairo R R，Sargent K A.Twin web disk:a step beyond con?vention[J].Journal of Engineering for Gas Turbines and Power，2002，124(2)：298—302.

[5]Harding B R，Curtiss D H.Contoured disk bore：US，7241111 B2[P].2007.

[6]Burge J C.Gas turbine compressor spool with structural and thermalupgrades：US，6267553B1[P].2001.

[7]江和甫.對(duì)渦輪盤材料的需求及展望[J].燃?xì)鉁u輪試驗(yàn)與研究，2002，15(4)：1—6.

[8]陸 山，李倫未.航空發(fā)動(dòng)機(jī)高負(fù)荷渦輪盤雙輻板結(jié)構(gòu)優(yōu)化設(shè)計(jì)[J].推進(jìn)技術(shù)，2011，32(5)：631—636.

[9]欒永先.雙輻板渦輪盤結(jié)構(gòu)強(qiáng)度分析[J].航空發(fā)動(dòng)機(jī)，2012，38(4)：38—41.

[10]欒永先.先進(jìn)渦輪盤結(jié)構(gòu)強(qiáng)度對(duì)比分析[J].航空發(fā)動(dòng)機(jī)，2013，39(3)：41—45.

[11]陸 山，趙 磊.雙輻板渦輪盤/榫結(jié)構(gòu)優(yōu)化設(shè)計(jì)方法[J].航空動(dòng)力學(xué)報(bào)，2014，29(4)：875—880.

[12]董少靜，申秀麗，康濱鵬，等.高推質(zhì)比雙輻板渦輪盤結(jié)構(gòu)研究及光彈試驗(yàn)驗(yàn)證[J].航空動(dòng)力學(xué)報(bào)，2015，30(1)：114—120.

[13]丁水汀，陶 智，徐國(guó)強(qiáng)，等.旋轉(zhuǎn)盤腔冷卻問題的工程評(píng)價(jià)[J].航空動(dòng)力學(xué)報(bào)，1999，14(1)：83—86.

[14]趙 熙，徐國(guó)強(qiáng)，羅 翔，等.等重量空心盤冷氣布置方案[J].北京航空航天大學(xué)學(xué)報(bào)，2009，35(5)：527—531.

[15]金 琰，郝艷華，黃致建.航空發(fā)動(dòng)機(jī)雙輻板渦輪盤溫度場(chǎng)與應(yīng)力場(chǎng)分析[J].鄭州輕工業(yè)學(xué)院學(xué)報(bào)，2015，30 (1)：81—84.

Flow and heat transfer for rotating cavity of tw in-web rotor disk

HAN Yu-qi，JIA Zhi-gang，LIU Hong，ZHU Da-ming
(Aero Engine Academy ofChina，Beijing 101304，China)

The three dimensionalmodelwas established with UG,the flow and heat transferwere simulated with CFX,the flow field and temperature field in typicalworking status of the rotating cavity for twin-web rotor disk were calculated and the conclusions of flow and heat transfer of rotating cavity were summarized. The results show that the flow structure in the rotating cavity is source-sink flow,total pressure and temper?ature of the cooling gas gradually increase along the radial direction;their increment increase with rotating speed and decrease with mass flow rate and heat transfer coefficient has different distributional features in different flow regions,which supports the furtheroptimized design of structure for twin-web rotor disk.

aero-engine；twin-web rotor disk；rotating cavity；cooling；flow；heat transfer；

V231.1；V232.3

A

1672-2620(2017)03-0048-05

2016-10-13；

2017-04-10

韓玉琪(1987-)，男，河南洛陽(yáng)人，工程師，博士，主要研究方向?yàn)楸P腔流動(dòng)換熱。