粗糙度對邊界層流動及壓氣機(jī)氣動性能的影響

孫海鷗， 葉楠， 王松， 王萌， 王忠義

(哈爾濱工程大學(xué) 動力與能源工程學(xué)院，黑龍江 哈爾濱 150001)

?

粗糙度對邊界層流動及壓氣機(jī)氣動性能的影響

孫海鷗， 葉楠， 王松， 王萌， 王忠義

(哈爾濱工程大學(xué) 動力與能源工程學(xué)院，黑龍江 哈爾濱 150001)

數(shù)值模擬；軸流壓氣機(jī)；表面粗糙度；邊界層流動；等效沙粒模型；工作效率；壓氣機(jī)

燃?xì)廨啓C(jī)作為一種先進(jìn)的動力設(shè)備，在艦船、輸氣、發(fā)電等領(lǐng)域有著廣泛的應(yīng)用。對于航母艦載機(jī)和長時間在海洋環(huán)境中工作的船用燃?xì)廨啓C(jī)來說，海洋大氣中的鹽霧氣溶膠、塵埃顆粒等使燃?xì)廨啓C(jī)通流部件承受著腐蝕、積垢及磨損等風(fēng)險[1-6]。腐蝕、積垢及磨損等是造成燃?xì)廨啓C(jī)葉片偏離設(shè)計工況甚至失效的常見形式，而葉片表面粗糙度增大則是最直接的表現(xiàn)形式，并且隨著運行時間其表面粗糙度也相應(yīng)增加。David Linden[7]對服役37個月的壓氣機(jī)進(jìn)口導(dǎo)葉進(jìn)行觀察，由于點蝕等造成孔徑甚至達(dá)到3.2 mm，深度為0.8 mm。在國內(nèi)，霍武軍等[8]對海軍航空兵所用發(fā)動機(jī)壓氣機(jī)葉片的觀察測量也表明海洋環(huán)境下葉片表面的點蝕蝕孔孔徑可達(dá)20～30 μm，由于蝕孔內(nèi)聚集H+、Cl-、SO42-等，有的蝕坑深度可達(dá)0.1～0.2 mm。盧明亮[9]在實驗室條件下對不銹鋼基材進(jìn)行鹽霧試驗，在經(jīng)過7、14和21 d后分別測量表面粗糙度發(fā)現(xiàn)由原來的2 μm分別增加到6.3、31.7和86.7 μm。由此可見，進(jìn)行葉片表面粗糙度對壓氣機(jī)性能影響的相關(guān)研究具有非常重要的意義。

在近幾十年中，國內(nèi)外的一些學(xué)者開展了很多針對壓氣機(jī)葉片表面粗糙度的相關(guān)研究。在國外，Syverud等基于對GE J85-13的實驗，得到葉片表面粗糙度對軸流壓氣機(jī)性能衰退的影響[10]。Back和Song等通過壓氣機(jī)平面葉柵粗糙度實驗，考察不同粗糙度大小和分布位置對葉柵壓力損失系數(shù)和效率損失系數(shù)的影響[11-12]。Mirko Morini等以NASA Stage37為對象，分別考察了均勻粗糙度和非均勻粗糙度對壓氣機(jī)總體性能衰退和內(nèi)部流動特征的影響[13-14]。在國內(nèi)，王松濤以NASA Stage35為計算模型分別計算了積垢和非均勻葉片粗糙度下壓氣機(jī)性能衰退規(guī)律和內(nèi)部流場特征的變化規(guī)律[15-16]。李冬使用數(shù)值計算方法研究了某型壓氣機(jī)葉片光潔和粗糙下的性能對比[17]。盧明亮使用改變?nèi)~型尺寸的方法研究了葉片腐蝕后粗糙度對壓氣機(jī)性能變化的影響，結(jié)果表明葉片腐蝕后的壓氣機(jī)性能出現(xiàn)了一定程度的衰退[9]。這些研究少有采用等效沙粒模型方法來研究葉片表面粗糙度對壓氣機(jī)性能的影響。本文采用數(shù)值模擬的方法，以T3系列平板轉(zhuǎn)捩實驗為基礎(chǔ)，驗證了等效沙粒粗糙度計算方法的可行性，再以等效沙粒粗糙度模型為基礎(chǔ)，以NASA Stage35型壓氣機(jī)為對象，研究了葉片表面粗糙度變化對壓氣機(jī)性能的影響。

1 等效沙粒模型數(shù)值方法

本文采用商業(yè)軟件CFX對研究對象進(jìn)行全三維數(shù)值求解。CFX使用的壁面函數(shù)方法是Launder和Spalding所提方法的擴(kuò)展。在對數(shù)規(guī)律的區(qū)域，流體近壁切向速度與壁面切應(yīng)力呈對數(shù)的關(guān)系，并且使用經(jīng)驗公式連接平均流動的近壁邊界條件和湍流輸運方程[18]。近壁速度的對數(shù)關(guān)系為

(1)

其中

(2)

(3)

式中：u+為近壁速度，uτ為摩擦速度，Ut為已知的距離壁面Δy距離處的切向速度，y+為距離壁面的無量綱距離，τw為壁面剪切應(yīng)力，κ為VonKarman常數(shù)，C為與壁面粗糙度相關(guān)的常數(shù)。

壁面表面粗糙度能夠明顯增加近壁湍流產(chǎn)生項，這又反過來導(dǎo)致壁面切應(yīng)力和傳熱系數(shù)的增加。為了使表面粗糙度的影響與實驗數(shù)據(jù)更好的吻合，CFX使用如下關(guān)系式考慮粗糙度的影響。

(4)

式中：B取5.2，偏移量ΔB是無量綱粗糙度h+(h+=huτ/υ)的函數(shù)。

對于砂粒粗糙度，偏移量ΔB可以表達(dá)成如下形式：

(5)

2 數(shù)值方法驗證

2.1 T3系列平板轉(zhuǎn)捩實驗

T3系列平板試驗是90年代ERCOFTAC(European Research Community on Flow， Turbulence and Combustion)在羅-羅公司所做的實驗[19]，是評價和驗證各種轉(zhuǎn)捩模型的經(jīng)典模型。T3系列平板分為零壓力梯度的T3A、T3B實驗和帶有壓力梯度的T3C2、T3C3、T3C4、T3C5等，實驗平板相同，前緣有圓角0.75 mm，長度為1.5 m。選取零壓力梯度的T3A和T3B進(jìn)行數(shù)值模擬來研究粗糙度對轉(zhuǎn)捩的影響，具體的邊界條件如表1所示。

表1 T3平板實驗的邊界條件

依據(jù)照實驗中的幾何尺寸建立了數(shù)值模擬所需的幾何模型并進(jìn)行了網(wǎng)格劃分。根據(jù)經(jīng)驗值，邊界層厚度約為2～3 mm，網(wǎng)格劃分保證至少40個節(jié)點，平板壁面第一層網(wǎng)格點的y+<1，流向網(wǎng)格節(jié)點數(shù)為610，高度方向網(wǎng)格點數(shù)為140，展向取15個網(wǎng)格點，網(wǎng)格總數(shù)約為138萬，如圖1所示。

2.2 光滑平板的數(shù)值模擬

本文分別使用不同的湍流模型和和兩種轉(zhuǎn)捩模型模擬平板邊界層轉(zhuǎn)捩流動過程，圖2是T3A和T3B對應(yīng)光滑條件下的不同模型的模擬結(jié)果以及與實驗值的對比。

圖1 T3平板計算域網(wǎng)格劃分示意圖Fig.1 The tablet domain mesh of T3 plate

圖2 T3A和T3B平板表面摩擦系數(shù)Cf沿流向分布Fig.2 Surface friction coefficient Cf along flow direction of T3A and T3B plate

圖中的參數(shù)Re為使用SST模型的Gamma轉(zhuǎn)捩模型時必須給定的轉(zhuǎn)捩開始動量厚度雷諾數(shù)，橫坐標(biāo)Rex為基于邊界層外主流速度和距平板前緣距離的雷諾數(shù)，Rex及表面摩擦系數(shù)Cf的義分別如下

(6)

(7)

式中：Ue為主流速度，ν為運動粘度系數(shù)，ρ為氣體密度，τ為流向方向的壁面切應(yīng)力。

從圖2(a)可以看出，SST模型的轉(zhuǎn)捩模型Gamma Theta和Gamma的計算結(jié)果與實驗值吻合較好，可以較準(zhǔn)確地預(yù)測邊界層轉(zhuǎn)捩過程，而其他模型的計算結(jié)果與實驗值差別很大，只是在湍流區(qū)域的計算值與實驗值吻合較好。

2.3 粗糙度表面數(shù)值模擬方法驗證

為了考察粗糙度對邊界層流動的影響，CFX軟件對粗糙度的處理方法是在壁面函數(shù)的求解中加入粗糙度的敏感項(等效砂粒粗糙度模型)，本文通過使用CFX軟件自帶的等效砂粒粗糙度模型和構(gòu)建Gauss型粗糙表面兩種方法進(jìn)行計算，考察粗糙度對邊界層流動參數(shù)的影響以及兩種方法計算結(jié)果的差別。基于上節(jié)的討論，對邊界層求解計算使用的湍流模型選用能夠較準(zhǔn)確預(yù)測轉(zhuǎn)捩過程的SST Gamma模型。

為了能夠觀察到粗糙度特征對邊界層流動影響的具體情況和特征，本節(jié)通過構(gòu)造服從Gauss分布的粗糙平板進(jìn)行數(shù)值模擬，進(jìn)而研究粗糙表面上的流動情況。

在笛卡爾坐標(biāo)系下，建立Gauss型粗糙表面模型時，x和y方向上的增量取Δx=40 μm、Δy=40 μm，基準(zhǔn)面取Z=0，通過控制輪廓峰值高度R的取值，分別建立具有不同粗糙高度值的表面，這里取輪廓算術(shù)平均偏差Ra作為粗糙度的衡量參數(shù)，其定義如下

(8)

這里共建立了3個粗糙表面模型，在進(jìn)行數(shù)值模擬時為了進(jìn)行對比，另取光滑平板模型的情況，分別列入表2，其中等效粗糙度ks與輪廓算數(shù)平均偏差Ra的經(jīng)驗關(guān)系式，取為ks=6.2Ra。

表2 建立的粗糙表面模型參數(shù)

把T3平板實驗的平板部分替換為構(gòu)造的Gauss型粗糙表面，前緣部分圓弧仍采用光滑過渡，計算域其他邊界不變，進(jìn)行網(wǎng)格生成。平板寬度取0.4 mm，網(wǎng)格節(jié)點數(shù)為12，流向取12 000網(wǎng)格節(jié)點，并且在x為0.2～0.8 m加密處理，圖4為網(wǎng)格的局部展示。

以表1中的邊界條件對具有Gauss型粗糙表面的平板進(jìn)行數(shù)值計算，同時按表2在CFX軟件中設(shè)

置等效粗糙度進(jìn)行計算，將計算結(jié)果表示在同一圖中以便進(jìn)行對比。圖4為兩種方法計算得到的平板表面摩擦系數(shù)結(jié)果。從圖中可以看出，對具有Gauss型粗糙度的T3平板的計算結(jié)果是波動的，這是因為表面摩擦系數(shù)的統(tǒng)計位置可能落在粗糙單元的峰頂或谷底，而兩種不同位置的剪切應(yīng)力存在著相對較大的差別，導(dǎo)至Cf值在某一小范圍內(nèi)跳動。同時可以看出Cf值的均值總體變化趨勢和不同流動狀態(tài)下粗糙度大小對其的影響與使用壁面函數(shù)法得到的計算結(jié)果都是一致的。

圖3 Gauss型粗糙度表面的網(wǎng)格劃分示意圖Fig.3 Mesh of the Gauss rough surface

圖4 兩種方法對應(yīng)的T3A和T3B的Cf計算值Fig.4 Cf of T3A and T3B corresponding two methods

通過以上研究可知，采用壁面函數(shù)(等效砂粒模型)的研究方法雖然不能完全得到邊界層內(nèi)的流動特征，但是從整體來看，兩種方法的計算結(jié)果是一致的。在工程問題中，一般只關(guān)注由粗糙度引起的宏觀結(jié)果，而對附面層內(nèi)的流動不作詳細(xì)研究。因此，使用壁面函數(shù)的方法考察粗糙度的影響具有一定的工程應(yīng)用價值，可以應(yīng)用在研究粗糙葉片對壓氣機(jī)性能影響的研究中。

3 均勻表面粗糙度對壓氣機(jī)性能的影響

3.1 物理模型及網(wǎng)格劃分

為了研究表面粗糙度對壓氣機(jī)性能的影響，本文以NASA Stage35試驗壓氣機(jī)為研究對象，該型壓氣機(jī)由一排轉(zhuǎn)子與一排靜子組成，其中轉(zhuǎn)子葉片數(shù)目為36，靜子葉片數(shù)目為46。設(shè)計轉(zhuǎn)速為17 188.7 r/min，最大葉尖速度為454.456 m/s。整個計算域采用HOH型網(wǎng)格結(jié)構(gòu)，并在葉片近壁區(qū)采用O型網(wǎng)格剖分，動葉葉尖間隙部分采用“蝶形網(wǎng)格”的結(jié)構(gòu)進(jìn)行劃分。對葉片附近網(wǎng)格進(jìn)行加密，設(shè)置第一層網(wǎng)格高度為10 μm。網(wǎng)格示意圖如圖5所示。在網(wǎng)格劃分完成后，葉片粗糙度按照CFX軟件中對粗糙表面的壁面函數(shù)處理方法給定等效砂粒粗糙度ks，分別計算不同表面粗糙度對壓氣機(jī)氣動特性的影響。其他邊界條件設(shè)定如下：

進(jìn)口：給定氣流總壓101 325Pa，總溫288.15K；

出口：按徑向平衡方程給定平均靜壓，調(diào)整壓力值模擬壓氣機(jī)不同的工況點。

壁面：絕熱，速度滿足無滑移條件。

圖5 壓氣機(jī)計算域網(wǎng)格劃分示意圖Fig.5 Mesh of the compressor domain

3.2 均勻粗糙度數(shù)值模擬結(jié)果

本文采用等效砂粒粗糙度模型作為壁面函數(shù)的方法進(jìn)行有關(guān)壓氣機(jī)表面粗糙度的數(shù)值模擬研究。表3為在動葉葉片上附加的粗糙度，相對粗糙度ks/c是以葉頂弦長為參考值，取等效粗糙度與參考值之比作為無量綱化的粗糙度值。其中ks=1μm時Rek<90，可以認(rèn)為是水力光滑，用以驗證CFX軟件中粗糙度模型的敏感性。

圖6為100%轉(zhuǎn)速下對動葉附加不同粗糙度計算得到的壓氣機(jī)工作特性線。首先可以看出對于Stage35壓氣機(jī)，在其近堵塞工況下流量變化很小而在近失速工況下其壓比幾乎不變，粗糙度的增加并未改變壓氣機(jī)特性曲線的變化趨勢和形態(tài)。對應(yīng)ks=1μm的壓氣機(jī)特性線與光滑情況幾乎重合，可以認(rèn)為流動處于水力光滑區(qū)。

表3 數(shù)值計算時葉片附加的粗糙度

圖6 不同動葉表面粗糙度下的壓氣機(jī)級工作特性線Fig.6 Characteristic line of the compressor with different roughness

從圖6中可以看出，當(dāng)動葉表面粗糙度ks由1μm增大至100μm時，壓氣機(jī)峰值效率減小3.5%，對應(yīng)工況點的壓比減小0.8%，總溫比升高0.4%，同時粗糙度增加后的壓氣機(jī)最大壓比與光滑情況下相比減小4.5%，這對于單級壓氣機(jī)來說是很大程度的衰退。

3.3 粗糙度對壓氣機(jī)損失的影響

(9)

式中：上標(biāo)“′”表示相對動葉葉片的參數(shù)，下標(biāo)“l(fā)e”和“te”分別表示前緣和尾緣，P和p則分別代表總壓和靜壓。這里以前緣處的總壓近似尾緣處理想過程的總壓。

為了更好地描述葉片粗糙度對壓氣機(jī)總體性能的影響，使用效率損失系數(shù)ζ來衡量表面粗糙度對壓氣機(jī)效率參數(shù)影響的函數(shù)關(guān)系，效率損失系數(shù)的定義[20]為

ζ=1-ηis

(10)

表4 峰值效率點對應(yīng)的參數(shù)計算結(jié)果

Table 4 Results of different cases corresponding peak efficiency point

編號ks/μmk+sωζ區(qū)域R110.390.13060.1672光滑R252.290.13320.1707粗糙R3104.960.13530.1740粗糙R42513.910.13990.1801粗糙R54023.690.14330.1843粗糙R65534.060.14670.1877粗糙R77044.870.14920.1908粗糙R88556.070.15220.1934粗糙R910067.590.15440.1960粗糙

圖7 總損失系數(shù)相對值smooth與的關(guān)系Fig.7 Relationship between smooth and

圖8 效率損失系數(shù)相對值ζ/ζsmooth與的關(guān)系Fig.8 Relationship between ζ/ζsmooth and

在一般工程問題中，在已知葉片平均粗糙度的情況下，根據(jù)上式可直接估算出壓氣機(jī)設(shè)計工況的總損失和效率損失。

4 結(jié)論

1)采用壁面函數(shù)(等效砂粒模型)的研究方法雖然不能完全得到邊界層內(nèi)的流動特征，但是該方法具有一定的工程價值，可以應(yīng)用在研究粗糙葉片對壓氣機(jī)性能影響的研究中。

2)動葉葉片表面粗糙度對壓氣機(jī)工作性能影響很大。當(dāng)動葉表面粗糙度ks由1 μm增大至100 μm時，壓氣機(jī)峰值效率減小3.5%，對應(yīng)工況點的壓比減小0.8%，總溫比升高0.4%，同時粗糙度增加后的壓氣機(jī)最大壓比與光滑情況下相比減小4.5%，這對于單級壓氣機(jī)來說是很大程度的衰退。

[1]孔祥鑫， 徐可君. 航空發(fā)動機(jī)葉片的腐蝕與防護(hù)[J]. 航空科學(xué)技術(shù)， 1997 (2): 24-26. KONG Xiangxin， XU Kejun. Aeroengine blade corrosion and its control[J]. Aeronautical science and technology， 1997 (2): 24-26.

[2]KURZ R， BRUN K， WOLLIE M. Degradation effects on industrial gas turbines[J]. Journal of engineering for gas turbines and power， 2009 (6): 62401-62407.

[3]GRANOVSKIY A， MANAEV I， VASSILIEV V， et al. Effects of blade degradation on turbine performance[C]//ASME 2013 Turbine Blade Tip Symposium. American Society of Mechanical Engineers， 2013: V001T01A003.

[4]CAGUIAT D E， ZIPKIN D M， PATTERSON J S. Compressor fouling testing on rolls royce/allison 501-K17 and general electric LM2500 gas turbine engines[C]//ASME Turbo Expo 2002: Power for Land， Sea， and Air. American Society of Mechanical Engineers， 2002: 933-942.

[5]HAMED A A， TABAKOFF W， RIVIR R B， et al. Turbine blade surface deterioration by erosion[J]. Journal of turbomachinery， 2005， 127(3): 445-452.

[6]王松，王國輝，韓青，等.葉片積垢對壓氣機(jī)性能衰退的影響[J].哈爾濱工程大學(xué)學(xué)報， 2014， (12): 1524-1528. WANG Song， WANG Guohui， HAN Qing， et al. Compressor performance deterioration caused by blade fouling[J]. Journal of Harbin Engineering University， 2014， (12): 1524-1528.

[7]LINDEN D H. Long term operating experience with corrosion control in industrial axial flow compressors[C]∥Proceedings of the Fortieth Turbomachinery Symposium.[S.l.], 2011: 93-106.

[8]霍武軍， 孫護(hù)國. 海航發(fā)動機(jī)壓氣機(jī)葉片腐蝕與防護(hù)措施[J]. 航空工程與維修， 2002， (6): 39-41. HUO Wujun， SUN Huguo. Corrosion and precaution for blades of aero-engine compressor in the navy air arm[J]. Journal of aviation maintenance and engineering， 2002 (6): 39-41.

[9]盧明亮. 鹽霧腐蝕對壓氣機(jī)性能影響的研究[D]. 哈爾濱工程大學(xué)， 2013. LU Mingliang. Reaearch on the effects of salts spray corrosion on the compressor performance[D]. Harbin: Harbin Engineering University， 2013.

[10]SYVERUD E， BAKKEN L E. The impact of surface roughness on axial compressor performance deterioration[C]// ASME Turbo Expo 2006: Power for Land， Sea， and Air. American Society of Mechanical Engineers， 2006: 491-501.

[11]BACK S C， HOBSON G V， SONG S J， et al. Effects of reynolds number and surface roughness magnitude and location on compressor cascade performance[J]. Journal of turbomachinery， 2012， 134(5): 051013.

[12]BACK S C， JEONG I C， SOHN J L， et al. Influence of surface roughness on the performance of a compressor blade in a linear cascade: experiment and modeling[C]//ASME Turbo Expo 2009: Power for Land， Sea， and Air. American Society of Mechanical Engineers， 2009: 239-247.

[13]MORINI M， PINELLI M， SPINA P R， et al. Computational fluid dynamics simulation of fouling on axial compressor stages[J]. Journal of engineering for gas turbines and power， 2010， 132(7): 072401.

[14]AIDI N， MORINI M， PINELLI M， et al. Performance evaluation of nonuniformly fouled axial compressor stages by means of computational fluid dynamics analyses[J]. Journal of turbomachinery， 2013， 136(2).

[15]CHEN S， SUN S， XU H， et al. Influence of local surface roughness of rotor blade on performance of an axial compressor stage[C]//ASME Turbo Expo 2013: Turbine Technical Conference and Exposition. American Society of Mechanical Engineers， 2013: V06AT35A015-V06AT35A015.

[16]CHEN S， ZHANG C， SHI H， et al. Study on the impact of fouling on axial compressor stage[C]//ASME Turbo Expo 2012: Turbine Technical Conference and Exposition. American Society of Mechanical Engineers， 2012: 1-8.

[17]李冬， 樊照遠(yuǎn)， 張娟，等. 壓氣機(jī)葉片粗糙度對其性能衰退的影響研究[J]. 航空發(fā)動機(jī)， 2009， (5): 32-35. LI Dong， FAN Zhaoyuan， ZHANG Juan， et al. Influence of blade roughness on compressor performance derterioration[J]. Aeroengine， 2009， (5): 32-35.

[18]CFX-Solver Theory Guide， ANSYS 13.0 Help[EB/OL]. 2010. http://www.doc88.com/p-908539666887.html.

[19]SAVILL A M. A synthesis of T3 test case predictions[J]. Numerical simulation of unsteady flows and transition to turbulence， 1992: 404-442.

[20]BOYLE R J， SENYITKO R G. Measurements and predictions of surface roughness effects on the turbine vane aerodynamics[C]//ASME Turbo Expo 2003， collocated with the 2003 International Joint Power Generation Conference. American Society of Mechanical Engineers， 2003: 291-303.

Influence of surface roughness on boundary-layer flow and characteristics of an axial compressor

SUN Haiou， YE Nan， WANG Song，WANG Meng， WANG Zhongyi

(College of power and energy engineering， Harbin engineering university， Harbin 150001， China)

numerical simulation; axial compressor; surface roughness; boundary-layer flow; equivalent sand model; efficiency

2016-08-02.

日期：2017-03-17.

國家自然科學(xué)基金項目(51309063，51679051);高等學(xué)校博士學(xué)科點專項科研基金新教師類資助課題(20132304120012).

孫海鷗(1962-)，男，教授，博士生導(dǎo)師; 王忠義(1982-)，男，副教授.

王忠義，E-mail: b205030024@126.com.

10.11990/jheu.201608006

V232.4

A

1006-7043(2017)04-0554-07

孫海鷗，葉楠，王松，等.粗糙度對邊界層流動及壓氣機(jī)氣動性能的影響[J]. 哈爾濱工程大學(xué)學(xué)報， 2017， 38(4): 554-560.

SUN Haiou， YE Nan， WANG Song， et al. Influence of surface roughness on boundary-layer flow and characteristics of an axial compressor [J]. Journal of Harbin Engineering University， 2017， 38(4): 554-560.

網(wǎng)絡(luò)出版地址：http://kns.cnki.net/kcms/detail/23.1390.u.20170317.0858.010.html