發(fā)動(dòng)機(jī)旋轉(zhuǎn)表面水滴撞擊特性數(shù)值研究

王治國(guó)，婁德倉(cāng)，郭文

(中國(guó)燃?xì)鉁u輪研究院，四川成都610500)

發(fā)動(dòng)機(jī)旋轉(zhuǎn)表面水滴撞擊特性數(shù)值研究

王治國(guó)，婁德倉(cāng)，郭文

(中國(guó)燃?xì)鉁u輪研究院，四川成都610500)

介紹了發(fā)動(dòng)機(jī)旋轉(zhuǎn)機(jī)械內(nèi)部過冷水滴軌跡的三維數(shù)值計(jì)算方法，闡述了旋轉(zhuǎn)坐標(biāo)系下氣流及粒子的運(yùn)動(dòng)規(guī)律。采用ANSYS-CFX軟件及其粒子輸運(yùn)模型，對(duì)某發(fā)動(dòng)機(jī)風(fēng)扇轉(zhuǎn)子葉片外圍空氣及水滴流場(chǎng)進(jìn)行了數(shù)值模擬。利用水滴的速度矢量、撞擊區(qū)域等參數(shù)表征水滴的撞擊特性，獲得了轉(zhuǎn)速、水滴直徑等對(duì)風(fēng)扇葉片表面水滴撞擊特性的影響：水滴撞擊區(qū)域集中在風(fēng)扇葉片迎風(fēng)面葉盆側(cè)，且水滴撞擊區(qū)域隨著轉(zhuǎn)速的增加而減??；水滴在葉片表面的撞擊范圍隨著水滴直徑的增大而減小。

航空發(fā)動(dòng)機(jī)；防冰；旋轉(zhuǎn)機(jī)械；水滴撞擊特性；粒子輸運(yùn)模型

particle transport model

1 引言

發(fā)動(dòng)機(jī)結(jié)冰對(duì)飛行器影響巨大，輕者造成發(fā)動(dòng)機(jī)輸出功率降低，嚴(yán)重時(shí)表面積冰脫落會(huì)損傷發(fā)動(dòng)機(jī)，造成機(jī)毀人亡[1，2]，因此世界各國(guó)的適航條例均對(duì)發(fā)動(dòng)機(jī)防冰能力提出了明確要求。飛行器表面結(jié)冰是由于大氣中含有的過冷水滴撞擊到迎風(fēng)表面后凝結(jié)而成，形成的冰晶等對(duì)翼型的氣動(dòng)性能有著重大影響[3]，過冷水滴對(duì)發(fā)動(dòng)機(jī)進(jìn)口部件迎風(fēng)表面的撞擊特性是防冰系統(tǒng)設(shè)計(jì)的基礎(chǔ)。文獻(xiàn)[4～6]對(duì)靜止部件表面水滴撞擊特性進(jìn)行了計(jì)算研究；劉華、朱劍鋆[7，8]等在此基礎(chǔ)上進(jìn)行了靜止部件表面防冰的換熱計(jì)算研究；Hamed等[9]采用三維計(jì)算方法分析了Rotor67轉(zhuǎn)子葉片表面的水滴撞擊特性，揭示了旋轉(zhuǎn)部件表面的流動(dòng)及撞擊特性。國(guó)內(nèi)在旋轉(zhuǎn)表面水滴撞擊特性方面研究較少。本文依據(jù)旋轉(zhuǎn)坐標(biāo)系下氣流及粒子的運(yùn)動(dòng)規(guī)律，采用ANSYS-CFX軟件及其粒子輸運(yùn)模型，對(duì)某發(fā)動(dòng)機(jī)風(fēng)扇轉(zhuǎn)子葉片外圍流場(chǎng)及其表面的水滴撞擊特性進(jìn)行計(jì)算分析。

2 水滴運(yùn)動(dòng)軌跡

為模擬空氣流場(chǎng)中水滴運(yùn)動(dòng)軌跡，需建立水滴運(yùn)動(dòng)方程。為便于分析，常做如下假設(shè)：①水滴直徑通常介于10～50 μm，水滴的存在不足以影響周圍空氣的流動(dòng)；②水滴密度為常數(shù)，即運(yùn)動(dòng)中沒有相變發(fā)生；③水滴為球形，且有統(tǒng)一的直徑；④初始速度等于自由流速度。

在歐拉坐標(biāo)系下建立控制方程?？諝庀嗟目刂品匠炭珊?jiǎn)化為：

過冷水滴在旋轉(zhuǎn)坐標(biāo)系下的運(yùn)動(dòng)方程簡(jiǎn)化為：

式中：rp、θp、zp定義了粒子在圓柱坐標(biāo)系下的位置，ωˉ代表葉片圓周速度。式(3)、式(4)中最后一項(xiàng)分別代表離心力和科氏力；三個(gè)方程中第一項(xiàng)代表因相對(duì)速度而產(chǎn)生的相間氣動(dòng)力Fˉ。重力、浮升力與氣動(dòng)力和離心力相比可忽略，對(duì)于小的過冷水滴，水滴粒子間的相互作用及壓力梯度可忽略。因此，F(xiàn)ˉ可用阻力系數(shù)和水滴滑移速度來表示：

式中：d表示水滴直徑；Vˉa、Vˉp分別表示空氣相和水滴的速度矢量；CD為相間阻力系數(shù)[10]，且有：

3 水滴撞擊特性

水滴撞擊率是運(yùn)動(dòng)氣流中水滴撞擊固體表面的速率，取決于來流、水滴特性及被撞擊物體的結(jié)構(gòu)。圖1為水滴撞擊示意圖，在接近和繞過葉型表面含水滴的空氣流場(chǎng)中，水滴慣性比空氣大，其繞過葉型的運(yùn)動(dòng)軌跡曲率比空氣流線小，在前緣區(qū)造成水滴撞擊現(xiàn)象[5]。圖中，SU、SL代表葉型表面計(jì)自駐點(diǎn)的弧長(zhǎng)，同時(shí)也分別為上、下表面水滴撞擊極限位置。θ為撞擊點(diǎn)水滴速度方向與微元表面切線方向的夾角。在水滴直徑較大情形下，可假設(shè)水滴不改變運(yùn)動(dòng)方向，此時(shí)θ等于物體表面切線方向與來流運(yùn)動(dòng)方向的夾角。葉型沿流線方向的投影高度(即支板迎風(fēng)寬度)為H。文中轉(zhuǎn)子葉片水滴撞擊極限Sm，指水滴與葉片表面左、右兩條相切軌跡所包圍的表面長(zhǎng)度S與H之比，即：

圖1 水滴撞擊示意圖Fig.1 Scheme of water droplet impingement

水滴撞擊效率β為實(shí)際水撞擊率與最大可能水撞擊率的比值。據(jù)此有

在采用CFX粒子離散相模型模擬水滴流場(chǎng)時(shí)，跟蹤粒子軌跡，可獲得粒子撞擊極限位置、水滴運(yùn)行參數(shù)等。通過這些參數(shù)，可表征葉型表面的水撞擊特性。

4 旋轉(zhuǎn)表面水滴撞擊特性模擬

為計(jì)算旋轉(zhuǎn)部件表面水滴撞擊特性，選取某發(fā)動(dòng)機(jī)進(jìn)氣風(fēng)扇葉片流道進(jìn)行計(jì)算。該風(fēng)扇有19個(gè)葉片，外徑154 mm，機(jī)匣與葉尖間隙0.25 mm。其具體性能參數(shù)見表1。

表1 某發(fā)動(dòng)機(jī)風(fēng)扇設(shè)計(jì)參數(shù)Table 1 Design parameters of an aero-engine fan

4.1 計(jì)算模型及網(wǎng)格

選取轉(zhuǎn)子葉片及進(jìn)出口流道段進(jìn)行計(jì)算。風(fēng)扇葉柵流道流場(chǎng)采用TURBOGRID進(jìn)行分網(wǎng)。葉片尖部采用局部加密，葉片表面網(wǎng)格及葉柵流道表面網(wǎng)格分別如圖2、圖3所示，網(wǎng)格單元數(shù)約10萬(wàn)。

圖2 葉片表面網(wǎng)格Fig.2 Grids of the blade surface

圖3 葉柵流道表面網(wǎng)格Fig.3 Grids of the cascade passage

4.2 水滴撞擊特性計(jì)算分析

4.2.1 水滴運(yùn)動(dòng)軌跡

風(fēng)扇流道內(nèi)部流場(chǎng)采用CFX軟件求解三維可壓縮RANS方程，對(duì)流項(xiàng)采用迎風(fēng)差分格式。湍流模擬采用k-ε兩方程模型，該流場(chǎng)在旋轉(zhuǎn)坐標(biāo)系下求解。邊界條件根據(jù)發(fā)動(dòng)機(jī)風(fēng)扇設(shè)計(jì)參數(shù)選取，水滴初始速度采用與空氣相同的進(jìn)口速度。

圖4給出了設(shè)計(jì)狀態(tài)下，不同相對(duì)葉高截面處的風(fēng)扇葉柵通道流場(chǎng)馬赫數(shù)分布。圖5給出了不同相對(duì)葉高截面在轉(zhuǎn)速30 000 r/min、d=20 μm下風(fēng)扇葉片流道內(nèi)水滴速度矢量。從圖中可看到，水滴在葉片表面不同截面位置的撞擊區(qū)域。

圖6給出了三維葉片表面水滴撞擊速度，其中藍(lán)色為無(wú)水滴撞擊區(qū)域，紅色則為水滴撞擊集中區(qū)域。從圖中可看到，葉片表面的水滴撞擊區(qū)域主要集中在葉盆側(cè)，該側(cè)相對(duì)于水滴的運(yùn)動(dòng)方向?yàn)橛L(fēng)面。葉盆表面下側(cè)尾緣處沒有水滴撞擊，這是由于水滴的慣性所致。圖6(a)還顯示出葉柵流道根部緣板表面的水滴撞擊區(qū)域。根部緣板由于沿流道方向變化，其迎風(fēng)表面存在水滴撞擊現(xiàn)象。

4.2.2 轉(zhuǎn)速對(duì)水滴撞擊特性的影響

圖4 風(fēng)扇通道內(nèi)馬赫數(shù)分布Fig.4 Mach number distribution inside the fan passage

圖5 水滴速度矢量分布Fig.5 Droplet velocity vector distribution

在轉(zhuǎn)速20 000～40 000 r/min范圍內(nèi)，數(shù)值模擬了轉(zhuǎn)速對(duì)水滴撞擊特性的影響。圖7給出了不同轉(zhuǎn)速下葉盆側(cè)水滴撞擊區(qū)域。從圖中看，水滴撞擊區(qū)域隨著轉(zhuǎn)速的增大逐漸減小，葉片根部無(wú)水滴區(qū)域逐漸擴(kuò)大。這是由于隨著離心力作用的增大，水滴徑向向外運(yùn)動(dòng)趨勢(shì)增強(qiáng)所致。從式(3)中可看到水滴運(yùn)動(dòng)受相間阻力、旋轉(zhuǎn)離心力及科氏力的共同作用，但在高轉(zhuǎn)速下，離心力影響起主導(dǎo)作用(105量級(jí))。

圖6 水滴撞擊區(qū)域Fig.6 Droplet impingement location

圖7 轉(zhuǎn)速對(duì)水滴撞擊特性的影響Fig.7 Effect of speed on the water droplet impingement characteristic

4.2.3 水滴直徑對(duì)撞擊特性的影響

圖8 不同直徑水滴速度矢量分布(10%相對(duì)葉高截面)Fig.8 Velocity vector distribution of droplet with different diameter(10%span from shroud)

在研究水滴直徑對(duì)葉片表面撞擊特性的影響時(shí)，主要計(jì)算了水滴直徑10～1 000 μm范圍水滴的運(yùn)動(dòng)軌跡。圖8給出了10%相對(duì)葉高截面上不同直徑水滴的速度矢量分布。從水滴的撞擊極限位置可看出，隨著水滴直徑的增大，水滴撞擊區(qū)域逐漸減小，但直徑超過200 μm后這種影響逐漸減小。這是由于相同密度下，水滴慣性隨尺寸的增大而增加。式(5)表明，水滴在軸向的運(yùn)動(dòng)加速度完全取決于相間阻力，隨著水滴直徑的增大，水滴運(yùn)動(dòng)慣性增強(qiáng)，阻力加速度相對(duì)減弱，即水滴保持其運(yùn)動(dòng)方向的能力增強(qiáng)。因此葉片表面撞擊范圍減小。

5 結(jié)論

(1)水滴撞擊區(qū)域集中在風(fēng)扇葉片迎風(fēng)面葉盆側(cè)，根部緣板表面迎風(fēng)位置也存在水滴撞擊現(xiàn)象。

(2)轉(zhuǎn)速增加導(dǎo)致離心力影響增強(qiáng)，使得水滴的撞擊區(qū)域隨著轉(zhuǎn)速的增加而減小。

(3)水滴在葉片表面的撞擊范圍隨著水滴直徑的增大而減小，但直徑超過200 μm后影響逐漸減小。

[1]裘燮綱，韓鳳華.飛機(jī)防冰系統(tǒng)[M].北京：航空專業(yè)教材編審組，1985：50—53.

[2]常士楠，艾素霄，畢文明，等.飛機(jī)發(fā)動(dòng)機(jī)進(jìn)氣道防冰系統(tǒng)的設(shè)計(jì)計(jì)算[J].北京航空航天大學(xué)學(xué)報(bào)，2007，33(6)： 649—652.

[3]王治國(guó)，楊軍，婁德倉(cāng).結(jié)冰對(duì)翼型流場(chǎng)影響的研究[J].燃?xì)鉁u輪試驗(yàn)與研究，2010，23(1)：31—36.

[4]Brette C，F(xiàn)regeau M，Trifu O，et al.A Three-Dimensional Water Droplet Trajectory and Impingement Analysis Pro?gram[R].AIAA 2005-4838，2005.

[5]楊軍，郭文，婁德倉(cāng).發(fā)動(dòng)機(jī)防冰支板水滴撞擊特性的數(shù)值研究[J].燃?xì)鉁u輪試驗(yàn)與研究，2011，24(1)，19—24.

[6]楊勝華，林貴平，申曉斌.三維復(fù)雜表面水滴撞擊特性計(jì)算[J].航空動(dòng)力學(xué)報(bào)，2010，25(2)：284—290.

[7]劉華，郭文，楊軍，等.發(fā)動(dòng)機(jī)進(jìn)氣帽罩防冰熱載荷的數(shù)值模擬研究[J].燃?xì)鉁u輪試驗(yàn)與研究，2012，25 (1)：44—48.

[8]朱劍鋆，董葳.防冰部件表面流動(dòng)換熱與溫度計(jì)算分析[J].燃?xì)鉁u輪試驗(yàn)與研究，2011，24(1)：15—18.

[9]Hamed A，Das K，Basu D.Numerical Simulations of Ice DropletTrajectoriesandCollectionEfficiencyon Aero-Engine Rotating Machinery[R].AIAA 2005-1248，2005.

[10]Morency F，Beaugendre H，Baruzzi G S，et al.FEN?SAP-ICE：A Comprehensive 3D Simulation System for In-Flight Icing[R].AIAA 2001-2566，2001.

Numerical Simulation of Water Droplet Impingement Characteristic on Aero-Engine Rotating Machinery

WANG Zhi-guo，LOU De-cang，GUO Wen
(China Gas Turbine Establishment，Chengdu 610500，China)

A methodology for three-dimensional numerical simulation of super-cooled water droplet trajec?tory through aero-engine rotating machinery is presented.Both flow and droplets'governing equations are described in the reference frame of rotating blades.By employing ANSYS-CFX and its particle transport model,the flow fields of both air and water droplet inside the fan blade passage are simulated.The water droplet impingement characteristic is displayed by its velocity vector and impingement limit.Simulation re?sults discover the impact of the rotational speed and droplet diameter on the impingement characteristic on the fan blade surface.The droplet impingement region is located on the pressure surface of the fan blade. The impingement area decreases with the increasing of the rotational speed and water droplet diameter.

aero-engine；anti-icing；rotating machinery；droplet impingement characteristic；

V231.3

A

1672-2620(2013)01-0035-05

2012-07-02；

2012-10-22

王治國(guó)(1981-)，男，山西運(yùn)城人，工程師，碩士，主要從事航空發(fā)動(dòng)機(jī)空氣系統(tǒng)與熱分析研究。