馬昌友, 侯敏杰, 凌代軍, 幸曉龍

(中國航空工業(yè)集團公司 中國燃?xì)鉁u輪研究院, 四川 江油 621703)

0 引 言

平面葉柵風(fēng)洞在葉柵性能基礎(chǔ)研究中占據(jù)了極為重要的地位，因此，對葉柵流場的準(zhǔn)確測量是至關(guān)重要的。目前通常采用三孔或五孔尾跡測壓探針對葉柵流場性能進(jìn)行測試[1]，其特點是只能進(jìn)行單點測量，對流場影響大，無法對葉片間槽道流場進(jìn)行測量，測量時間長、效率低, 并且得到的是平均意義下的氣流參數(shù)。當(dāng)葉柵發(fā)生較嚴(yán)重的氣流分離時，柵后尾跡區(qū)的氣流角變化較大，并且存在氣流旋渦，真實氣流方向與探針軸線夾角遠(yuǎn)遠(yuǎn)大于探針校準(zhǔn)角度，超出了探針應(yīng)用范圍，這時測得的值準(zhǔn)確度降低，往往不能正確地反映實際流動[2]。

近年來，非接觸式瞬態(tài)流場測速技術(shù)PIV (Particle Image Velocimetry)，克服了接觸式單點測量設(shè)備的局限性,可以無干擾地獲得葉柵槽道瞬時速度場、渦量場等, 已成為研究葉柵流場的一種先進(jìn)測試手段[3-4]。國內(nèi)外在葉柵流場PIV測試技術(shù)研究方面已取得了卓著的成績[5-10]，但由于PIV 技術(shù)比較復(fù)雜，影響測試結(jié)果因素比較多，如示蹤粒子的選取與播撒[11-12]、激光片光與CCD的布局等，都會影響所測流場的準(zhǔn)確性，國內(nèi)在高速風(fēng)洞中開展的PIV應(yīng)用研究較少。

作者基于暫沖式平面葉柵風(fēng)洞，利用2D-PIV系統(tǒng)對某亞聲速擴壓葉柵流場進(jìn)行了測量，以解決示蹤粒子撒播、PIV光路布局及標(biāo)定等技術(shù)問題，并對零迎角不同進(jìn)口馬赫數(shù)狀態(tài)下的葉柵中截面槽道及尾跡速度場分布進(jìn)行分析研究。為驗證PIV測試結(jié)果的準(zhǔn)確性，還利用三孔尾跡探針的測試結(jié)果和數(shù)值計算結(jié)果對2D-PIV在葉柵流場測試結(jié)果的可靠性進(jìn)行驗證。

1 實驗設(shè)置

本實驗在中國燃?xì)鉁u輪研究院暫沖式超、跨聲速平面葉柵風(fēng)洞試驗器上進(jìn)行[13]。葉型為某亞聲速擴壓葉型，弦長B=65mm，柵距T=49.1mm，安裝角γ=60.65°，設(shè)計進(jìn)口氣流角β1=45.7°，設(shè)計進(jìn)口馬赫數(shù)M1=0.66。圖1為葉柵試驗件實物圖，與常規(guī)葉柵設(shè)計相比，該葉柵在設(shè)計允許的條件下，增大柵距，以便減小由于葉片彎曲帶來的光路遮擋；葉柵一側(cè)柵板材料為透明的航空有機玻璃，兼作CCD拍照視窗；葉柵另一側(cè)柵板為普通鋼柵板，以便在葉柵內(nèi)部形成暗室，降低PIV測量時背景光的影響。

圖1 葉柵試驗件實物圖

實驗所使用的粒子圖像測速系統(tǒng)從德國LaVision公司引進(jìn)。該系統(tǒng)是由雙脈沖激光器、CCD、同步控制器、圖像采集和裝有Davis分析軟件的計算機組成。激光器為雙諧振脈沖式Nd：YAG激光器，工作頻率為15Hz，單脈沖最大能量360mJ，激光波長為532nm，通過激光臂，可以方便地將激光打入到葉柵槽道中。CCD相機緩存為512MB，像素1600×1200，拍攝面積85mm×64mm。相機支架可對CCD在3個方向上進(jìn)行精確旋轉(zhuǎn)，再配合三維位移機構(gòu)，可全方位控制CCD的方向和位置。試驗時，Davis軟件通過同步器控制激光打光和CCD拍照,并進(jìn)行圖像處理，獲得葉柵流道速度場分布。

為了驗證PIV試驗數(shù)據(jù)的準(zhǔn)確性，采用三孔楔形尾跡探針對葉柵性能進(jìn)行了常規(guī)測量。三孔楔形尾跡探針安裝在位移機構(gòu)上，可以控制尾跡探針沿柵距方向移動。

2 實驗方案

2.1示蹤粒子投放

粒子發(fā)生器選用德國Lavision公司的Laskin PIV Part40壓力霧化式粒子發(fā)生器，需使用高壓氣源供氣，最大允許進(jìn)、出口壓差可達(dá)0.6MPa。示蹤劑采用Lavision公司研制的癸二酸(2-乙基已基)酯(簡稱為DEHS)，該示蹤粒子適用于風(fēng)洞、葉輪機械、發(fā)動機內(nèi)流等流場，粒子平均半徑為1μm左右，優(yōu)點是無毒，穩(wěn)定性好，可以在不同的情況下存留幾個小時保持尺寸形狀不變。文獻(xiàn)[6]利用PIV開展渦輪葉片尾緣噴氣研究中使用了該示蹤粒子，使用效果良好。

考慮到試驗器是暫沖式的，來自氣源站的高壓空氣需經(jīng)過凈化器、干燥器和儲氣罐后，送至試驗段，故示蹤粒子不能直接從氣源站的壓縮機進(jìn)口吸入，而是在試驗段上游管道中注入，如圖2所示。粒子發(fā)生器所需的高壓氣源來自于氣源站儲氣系統(tǒng)。試驗前,來自氣源站的壓縮空氣先為儲氣系統(tǒng)供氣充壓至0.6MPa，將小儲氣罐通過閘閥與大儲氣罐斷開；試驗時小儲氣罐獨立作為粒子發(fā)生器的氣源，大儲氣罐作為葉柵主流氣源；當(dāng)小儲氣罐壓力較低時，將閘閥打開，與大儲氣罐連通，氣源站即可為小儲氣罐再次打壓至0.6MPa。示蹤粒子流通過安裝在試驗器上游的調(diào)壓閥與穩(wěn)壓箱之間的撒播器投放到主流中。調(diào)壓閥后主流管道的壓力一般在200kPa以內(nèi)，足夠的壓差使得示蹤粒子流可以壓入主流中。另外，撒播器距離試驗段約8m，中間經(jīng)擴壓器、蜂窩整流器、絲網(wǎng)整流器，確保示蹤粒子在被主流輸送到試驗段前與主流充分混合均勻。

圖2 供氣流程圖

2.2光路布置與標(biāo)定

圖3為試驗時激光片光與CCD的布局方案。試驗采用從葉柵下游即試驗段出口處逆向打光，為減小葉柵出口氣流對光學(xué)部件的沖擊力，打光方向適當(dāng)避開氣流方向，并確保激光片光與待測截面重合。由于采用單CCD對葉柵中截面槽道及其下游尾跡速度場進(jìn)行二維速度場測量，因此要求CCD拍照方向與激光片光嚴(yán)格垂直。

圖3 激光片光與CCD布局方案

試驗前首先對測試區(qū)域進(jìn)行標(biāo)定。由于視窗為有機玻璃柵板，對粒子成像的畸變影響較小，加之葉柵通道不便于放置標(biāo)定靶，因此直接將標(biāo)有刻度的校準(zhǔn)尺置于觀測區(qū)域獲取圖像，如圖4所示，用于標(biāo)定CCD像素與實際測量范圍的換算關(guān)系，從而可實現(xiàn)對流場結(jié)構(gòu)定量測量和標(biāo)定。

3 數(shù)值模擬方法簡介

在進(jìn)行PIV測量結(jié)果與三孔楔形尾跡探針測量結(jié)果對比之前，首先進(jìn)行PIV測量結(jié)果與數(shù)值模擬結(jié)果對比分析，以考察PIV測量結(jié)果的可靠性。作者采用全三維粘性流場數(shù)值模擬的手段，計算獲得了該葉柵在相應(yīng)狀態(tài)下的流場信息。計算使用商用NUMECA軟件的FINE/Turbo模塊，采用時間推進(jìn)法求解三維雷諾平均的N-S(RANS)方程，采用二階精度的中心差分格式進(jìn)行空間離散，時間項采用四階Runge-Kutta法迭代求解，湍流模型為一方程的S-A(Spalart- Allmaras)模型。計算采用多重網(wǎng)格法，結(jié)合當(dāng)?shù)貢r間步長及隱式殘差光順等加速收斂技術(shù)。計算域的網(wǎng)格劃分，采用了結(jié)構(gòu)化網(wǎng)格的HOH型網(wǎng)格，計算的總網(wǎng)格量約30萬，壁面處理采用了絕熱無滑移條件。

圖4 二維葉柵流場PIV標(biāo)定方法

4 結(jié)果分析

4.1葉柵流場周期性

為保證試驗數(shù)據(jù)準(zhǔn)確可靠，PIV測試前需在葉柵進(jìn)、出口建立周期性流場。圖5繪出了進(jìn)口氣流角為β1=45.7°時，葉柵進(jìn)口馬赫數(shù)分別為M1=0.5、0.66、0.73和0.8時柵前馬赫數(shù)沿額線方向3個柵距的分布。由圖可見，柵前進(jìn)口馬赫數(shù)沿柵距分布周期性是比較好的。在出口流場方面，圖6繪出了柵后測量截面兩個柵距的尾跡總壓恢復(fù)系數(shù)分布圖，也說明試驗流場周期性是比較好的。基于此，本文只對一個柵距的葉柵流場進(jìn)行了PIV測量。另外通過圖6定性地比較和分析尾跡區(qū)寬度、深度等表征葉柵損失大小的特征，可以發(fā)現(xiàn)，隨著進(jìn)口馬赫數(shù)的提高，葉柵損失呈現(xiàn)快速上升的趨勢。

圖5 葉柵進(jìn)口速度沿額線分布

圖6 葉柵出口尾跡參數(shù)分布

4.2葉柵全流場PIV測量與數(shù)值計算結(jié)果對比

圖7給出了在設(shè)計進(jìn)口氣流角β1=45.7°，進(jìn)口馬赫數(shù)M1分別為0.6和0.73狀態(tài)下，利用PIV系統(tǒng)采集并經(jīng)粒子圖像計算得到的葉柵流場的瞬態(tài)速度矢量圖。PIV圖像分析查問域(Interrogating Window)為64pixel×64pixel，成像系數(shù)為0.05333，合3.4mm×3.4mm，重疊率為50%。由圖7可見，氣流在葉背處先加速至最大值后逐漸減速增壓。從圖中可以看到葉片尾緣處的速度值較小，該處為葉柵尾跡區(qū)，但由于本文關(guān)注葉柵槽道及尾跡流場的全域測量，未能獲得尾跡區(qū)的細(xì)微結(jié)構(gòu)。由于葉片的存在，使得靠近葉片表面區(qū)域沒有速度向量。主要因素有片光被彎曲的葉片遮擋、CCD被葉片端截面遮擋、葉片附近的粒子反光不規(guī)律或粒子濃度過低，這些因素不可能根本消除，只能盡可能減輕近壁面問題的影響。計算中通過創(chuàng)建遮罩Mask方法來屏蔽該區(qū)域。

圖8給出了在設(shè)計進(jìn)口氣流角β1=45.7°，進(jìn)口馬赫數(shù)M1分別為0.6和0.73狀態(tài)下，通過CFD數(shù)值模擬得到葉柵流場的速度矢量圖。對比PIV測量結(jié)果和CFD數(shù)值模擬結(jié)果可以看出，PIV測得的葉柵中截面二維速度矢量場合理地反映了葉片槽道及尾跡的流動結(jié)構(gòu)。從PIV試驗和CFD計算結(jié)果可以看出，在設(shè)計迎角下，葉背沒有發(fā)生氣流分離。PIV捕捉到了葉柵尾跡，從圖中可以看到葉片尾緣處的速度向量較小。但由于本文關(guān)注葉柵槽道及尾跡流場的全域測量，未能捕捉到尾跡區(qū)的細(xì)微結(jié)構(gòu)。

另外，從圖7可以看出，圖中局部小區(qū)域的速度矢量與圖8相應(yīng)位置存在較大差異，主要集中在葉背，這是因為該位置處的葉柵柵板受到液體示蹤劑的污染，并不斷聚集，導(dǎo)致該區(qū)域的粒子圖像沒能合理反映出葉柵中截面的流動結(jié)構(gòu)。因此今后開展葉柵流場測量，特別是在較高馬赫數(shù)或大迎角下，需要考慮壁面污染對測量的影響。

(a) β1=45.7°,M1=0.6

(b) β1=45.7°,M1=0.73

(a) β1=45.7°,M1=0.6

(b) β1=45.7°,M1=0.73

4.3葉柵出口流場PIV與尾跡探針測量結(jié)果對比

由于尾跡探針測量得到的是穩(wěn)態(tài)定常數(shù)據(jù)，而PIV測得的為瞬態(tài)數(shù)據(jù)。為了進(jìn)行尾跡探針和PIV測量結(jié)果對比， PIV在同一狀態(tài)進(jìn)行了30次速度場數(shù)據(jù)采集，對30組瞬態(tài)速度場數(shù)據(jù)進(jìn)行了平均，從而得到時均意義下穩(wěn)態(tài)速度場。

圖9為進(jìn)口馬赫數(shù)為0.60和0.73時采用PIV和尾跡探針測得的葉柵出口距離尾緣45%T位置處氣流速度W2沿一個柵距的分布。從圖中可以看出，PIV和尾跡探針兩者測得的出口氣流速度沿柵距分布規(guī)律非常接近，并且相同位置的速度值差異較小。圖10則為相應(yīng)的葉柵出口氣流角β2沿一個柵距的分布。從圖中可以看出，PIV和尾跡探針兩者測得的出口氣流角沿柵距分布規(guī)律在主流區(qū)比較一致，但在尾跡區(qū)差異略大。

圖9 葉柵出口速度沿柵距分布

對圖9和10中尾跡探針和PIV測得的葉柵出口速度值和氣流角分別在一個柵距內(nèi)做算術(shù)平均，即可得到葉柵出口遠(yuǎn)場的氣流速度和方向。圖11和12分別繪出了葉柵出口速度值和氣流角隨進(jìn)口馬赫數(shù)的變化趨勢。由圖11可見，從尾跡探針和PIV獲得的葉柵遠(yuǎn)場出口氣流速度隨著進(jìn)口馬赫數(shù)的提高均呈現(xiàn)線性增大趨勢，并且兩種測量手段獲得的速度值差別較小。由圖12可見，葉柵在設(shè)計進(jìn)口馬赫數(shù)以下工作時，從尾跡探針和PIV獲得的葉柵遠(yuǎn)場出口氣流角均不隨進(jìn)口馬赫數(shù)發(fā)生明顯的變化，并且兩者測得值也較接近；當(dāng)葉柵在設(shè)計進(jìn)口馬赫數(shù)以上工作時，此時葉柵進(jìn)口馬赫數(shù)超過臨界馬赫數(shù)，甚至接近最大馬赫數(shù)(此時葉柵堵塞)，葉片邊界層分離嚴(yán)重，導(dǎo)致葉柵出口氣流角隨著進(jìn)口馬赫數(shù)繼續(xù)升高呈現(xiàn)明顯下降的趨勢，此時尾跡探針和PIV獲得的葉柵遠(yuǎn)場出口氣流角誤差較大。

圖10 葉柵出口氣流角沿柵距分布

圖11 葉柵出口速度隨進(jìn)口馬赫數(shù)的變化趨勢

圖12 葉柵出口氣流角隨進(jìn)口馬赫數(shù)的變化趨勢

5 誤差分析

從試驗結(jié)果對比分析來看，PIV與三孔楔形探針對平面擴壓葉柵流場測量結(jié)果基本吻合，但仍存在誤差，特別是葉柵出口氣流角的測量誤差較大。經(jīng)分析，產(chǎn)生的誤差主要存在如下兩個方面:

一是尾跡區(qū)的氣流方向沿額線變化較大，而PIV圖像處理時查問域設(shè)置過大，導(dǎo)致未能捕捉到尾跡區(qū)的細(xì)微結(jié)構(gòu)。由于進(jìn)口馬赫數(shù)較大時，葉柵流場的速度變化范圍變大，若將查問域設(shè)置過小，會導(dǎo)致葉背處速度較大區(qū)域沒有有效的粒子對。如只對葉柵出口尾跡區(qū)進(jìn)行PIV測量，試驗時選取合適的兩束激光時間差Δt，圖像處理時采用較小的查問域，便可提高對尾跡區(qū)流場細(xì)微結(jié)構(gòu)的捕捉能力。

二是三孔楔形探針的方向特性線較差，特別是在尾跡區(qū)真實氣流方向與探針安裝方向偏角大于探針校準(zhǔn)偏角時，用探針數(shù)據(jù)外插所得到的出口氣流角偏差較大，從而進(jìn)一步加大與PIV測量結(jié)果的差異。這一點也顯示出了PIV在葉柵尾跡測量比三孔楔形尾跡探針更有優(yōu)勢。

6 結(jié) 論

基于暫沖式葉柵風(fēng)洞，采用數(shù)字PIV技術(shù)，并配合粒子投放和光路布置技術(shù)，對某擴壓葉型葉柵槽道及尾跡成功地進(jìn)行了可視化測量。試驗結(jié)果對比分析表明：采用的PIV可以應(yīng)用于葉柵流場的精確測量。

PIV與三孔楔形尾跡探針在葉柵出口尾跡的測量所獲得的氣流速度和主流區(qū)的出口氣流角重合性較好；尾跡分離區(qū)的出口氣流角重合性略差，主要原因是在變化較大的尾跡區(qū)氣流偏角超出了三孔尾跡探針校準(zhǔn)范圍，恰好凸出了PIV技術(shù)優(yōu)勢。但對于跨聲速葉柵流場，葉柵槽道及尾跡區(qū)流場速度變化較大，建議采用分塊開展PIV測量，后期采用拼接處理，可提高對尾跡區(qū)等流場細(xì)微結(jié)構(gòu)的捕捉能力。本文提出的PIV測量技術(shù)也可用于連續(xù)式葉柵風(fēng)洞中。

參考文獻(xiàn):

[1]李鋼, 聶超群, 朱俊強, 等. 介質(zhì)阻擋放電等離子體對壓氣機葉柵性能影響的實驗[J]. 航空動力學(xué)報, 2008, 23(3): 522-526.

Li Gang, Nie Chaoqun, Zhu Junqiang, et al. Experimental investigations of the effect of dielectric barrier discharge plasma on performance of compressor cascade[J]. Journal of Aerospace Power, 2008, 23(3): 522-526.

[2]陳武, 侯安平, 鄭新前, 等. 在大分離區(qū)內(nèi)使用三孔探針測試的探索[J]. 流體力學(xué)實驗與測量, 2003, 17(3): 84-87,92.

Chen Wu, Hou Anping, Zheng Xinqian, et al. Research of test by three-hole pressure probe in large-separated zone[J]. Experiments and Measurements in Fluid Mechanics, 2003, 17(3): 84-87, 92.

[3]李志平, 李秋實, 袁巍, 等. 葉輪機轉(zhuǎn)子葉排非定常旋渦脫落頻譜特性研究新方法[J]. 實驗流體力學(xué), 2006, (02): 7-11.

Li Zhiping, Li Qiushi, Yuan Wei, et al. A new method for studying the unsteady vortex frequency characteristics on rotor blades of turbine engine[J]. Journal of Experiments in Fluid Mechanics, 2006, (02): 7-11.

[4]梁華, 李應(yīng)紅, 宋慧敏, 等. 等離子體氣動激勵誘導(dǎo)空氣流動的PIV研究[J]. 實驗流體力學(xué), 2011, 25(04): 22-25,36.

Liang Hua, Li Yinghong, Song Huimin, et al. PIV investigation on flow induced by plasma aerodynamic actuation[J]. Journal of Experiments in Fluid Mechanics, 2011, 25(04): 22-25,36.

[5]Jordi Estevadeordal, Matthew D Langford, et al. PIV investigation of the flow field in a transonic compressor cascade with a moving shock wave[R]. AIAA 2007-506433, 2007.

[6]Chase Nessler, Christopher Marks. PIV investigation of periodic unsteady wakes over a highly loaded LPT blade[R]. AIAA 2009-5107, 2009.

[7]劉波, 項效镕, 南向誼, 等. 附面層抽吸對葉柵表面分離流動控制的實驗研究[J]. 推進(jìn)技術(shù), 2009, 30(6)：703-708.

Liu Bo, Xiang Xiaorong, Nan Xiangyi, et al. Experimental investigation for suppression boundary layer separation on cascade surface by BLS[J]. Journal of Propulsion Technology, 2009, 30(6)：703-708.

[8]孫衛(wèi)亮, 宋彥萍, 谷君, 等. 大迎角平面壓氣機葉柵內(nèi)的PIV 應(yīng)用[J]. 節(jié)能技術(shù), 2006, 24(4): 291-293.

Sun Weiliang, Song Yanpin, Gu Jun, et al. Application of PIV in planar compressor cascade with large positive attack angle[J]. Energy Conservation Technology, 2006, 24(4): 291-293.

[9]戴靜君, 姜義忠, 董守平. PIV技術(shù)在渦輪葉柵內(nèi)流場試驗中的應(yīng)用[J]. 流體力學(xué)實驗與測量, 2003, (04): 68-70,83.

Dai Jingjun, Jiang Yizhong, Dong Shouping. The application of particle image velocimetry to fluid experimentation in a turbine cascade[J]. Experiments and Measurements in Fluid Mechanics, 2003, (04): 68-70,83.

[10] Freund Oliver, Rehder Hans-juergen, Schaefer Philipp, et al. Experimental investigations on cooling air ejection at a straight turbine cascade using PIV and QIS[R]. ASME GT2011-45296.

[11] 榮臻, 陳方, 劉洪, 等. 超聲速PIV示蹤粒子布撒技術(shù)研究[J]. 實驗流體力學(xué), 2012, (02): 64-67,85.

Rong Zhen, Chen Fang, Liu Hong, et al. Research on seeding technique of PIV in supersonic flow[J]. Journal of Experiments in Fliud Mechanics, 2012, (02): 64-67,85.

[12] 陳瑩, 鄧學(xué)鎣, 王延奎, 等. 低速風(fēng)洞PIV實驗中的示蹤粒子投放技術(shù)[J]. 實驗流體力學(xué), 2007, 21(02): 78-81.

Chen Ying, Deng Xueying, Wang Yankui, et al. Testing technique of PIV in low-speed wind tunnel[J]. Journal of Experiments in Fluid Mechanics, 2007, 21(02): 78-81.

[13] 凌代軍, 姜正禮, 仲永興. 平面葉柵非定常流動試驗方法研究[J]. 燃?xì)鉁u輪試驗與研究, 2006, 19(3): 24-29,36.

Ling Daijun, Jiang Zhengli, Zhong Yongxing. Experimental procedure study on unsteady flow in plane cascade[J]. Gas Turbine Experiment and Research, 2006, 19(3): 24-29,36.

作者簡介:

馬昌友(1981-)，男，安徽天長人，高級工程師，在讀博士。研究方向：氣動葉片試驗技術(shù)研究。通訊地址：四川省江油市305信箱三室(621703)。E-mail：gtemacy@163. com