胡金鑫, 鄭光華, 孟玉航, 魏邦, 余柯鋒

(1.西北工業(yè)大學(xué) 動(dòng)力與能源學(xué)院, 陜西 西安 710129;2.中國(guó)航發(fā)湖南動(dòng)力機(jī)械研究所 測(cè)試技術(shù)部, 湖南 株洲 412002)

在小型發(fā)動(dòng)機(jī)的整機(jī)測(cè)試中，傳感器的大量布置會(huì)在一定程度上影響整機(jī)流場(chǎng)，測(cè)試結(jié)果也會(huì)受到流場(chǎng)的影響，進(jìn)而將影響發(fā)動(dòng)機(jī)的性能與調(diào)整[1]。探針作為一種接觸式測(cè)量手段，引起的堵塞效應(yīng)不可避免地對(duì)流場(chǎng)產(chǎn)生干擾，使性能下降甚至引起失速。圍繞探針上下游尾跡國(guó)內(nèi)外做過(guò)不少研究[2]。

Jan[3]采用實(shí)驗(yàn)的方法研究了壓氣機(jī)轉(zhuǎn)子中的旋轉(zhuǎn)探針引起的堵塞效應(yīng)，結(jié)果表明探針會(huì)降低葉柵通道中軸向流速。

Simon等[4]研究發(fā)現(xiàn)探針置于靜子葉片附近造成的堵塞效應(yīng)導(dǎo)致葉片和探頭之間的區(qū)域中壓力增加和速度降低并形成較大區(qū)域的尾跡，從而導(dǎo)致停滯點(diǎn)的偏移和相鄰探頭壓力端口的錯(cuò)誤讀數(shù)。

余柯鋒[5]研究了整流套對(duì)葉型探針測(cè)量結(jié)果的影響。結(jié)果表明在葉型總壓探針前端增加整流套，可有效增加探針的不敏感角以及高馬赫數(shù)的測(cè)量準(zhǔn)確度。

魏崇等[6]研究發(fā)現(xiàn)隨著馬赫數(shù)增加，探針對(duì)下游流場(chǎng)的影響程度越大，對(duì)下游流場(chǎng)流向的影響距離最大可達(dá)50倍探針支桿直徑，對(duì)周向的影響距離最大可達(dá)15°。

高杰等[7]研究分析了不同結(jié)構(gòu)的探針支桿對(duì)壓氣機(jī)葉柵性能的影響，發(fā)現(xiàn)總壓探針支桿會(huì)增加葉柵和對(duì)應(yīng)葉片的總壓損失，當(dāng)探針安裝在葉柵通道中部時(shí)，支桿對(duì)下游流場(chǎng)的影響最小。

付少林等[8]研究發(fā)現(xiàn)當(dāng)探針位于近壁面時(shí)，渦輪測(cè)試效率具有較高的精度，當(dāng)探針距離輪轂5%～10%葉高及距離機(jī)匣90%～95%葉高時(shí)，渦輪測(cè)試效率精度較高；探針位于渦輪出口下游3倍轉(zhuǎn)子弦長(zhǎng)以上時(shí)，在不同的工況下，渦輪修正效率與數(shù)值計(jì)算效率的誤差小于0.3%。

朱高平等[9]對(duì)比分析不同周向分布位置的探針耙對(duì)葉柵性能參數(shù)的影響，發(fā)現(xiàn)總壓探針耙會(huì)增加下游的流速，探針尾緣形成的激波與相鄰葉片尾緣的激波相互干涉，造成上游葉柵通道堵塞以及葉片的載荷分布不相同。

王肖等[10]研究發(fā)現(xiàn)總壓探針的插入使得被測(cè)流場(chǎng)流量減小，堵塞效應(yīng)受探針結(jié)構(gòu)的影響較小。

鄭光華等[11]基于數(shù)值計(jì)算研究發(fā)現(xiàn)相對(duì)近壁距離越小，總壓探針測(cè)量誤差越大，且當(dāng)近壁距離大于3倍探針直徑時(shí)，臨壁效應(yīng)消失。

本文采用數(shù)值模擬的方法，研究總壓探針在壓氣機(jī)進(jìn)口級(jí)出口流道中的測(cè)量結(jié)果以及探針對(duì)流場(chǎng)的影響。針對(duì)探針的測(cè)量精度，分析探針測(cè)量誤差隨發(fā)動(dòng)機(jī)工作狀態(tài)(巡航、最大連續(xù)、起飛)而變化的情況，重點(diǎn)分析不同安裝位置下，探針測(cè)量結(jié)果的變化趨勢(shì)，分析葉盆、葉背附近的壓力場(chǎng)與速度場(chǎng)；針對(duì)總壓損失系數(shù)，分析不同安裝位置下，流向曲線上總壓損失系數(shù)的變化情況；分析不同下游距離截面的速度均勻度。

1 數(shù)值模型

1.1 計(jì)算模型

計(jì)算域選取某渦軸發(fā)動(dòng)機(jī)壓氣機(jī)進(jìn)口級(jí)出口流道，如圖1所示，流道中流體做向心運(yùn)動(dòng)，帶有預(yù)旋進(jìn)入葉柵通道，平行中心軸線方向流出。入口截面為圓柱面，圓柱直徑為452 mm，入口截面寬度11 mm，出口截面為圓環(huán)截面，內(nèi)徑為44 mm，外徑為128 mm，整個(gè)流域在軸向上長(zhǎng)59 mm，在流域內(nèi)周向均勻分布20處葉片。軸向沿流向?yàn)閤軸正方向，中心軸線到進(jìn)口截面方向?yàn)閺较蚍较?，?duì)應(yīng)z軸正方向，另一方向切向，即為y軸。

圖1 計(jì)算模型

圖2為單點(diǎn)總壓探針模型?？倝禾结樃袎汗芸傞L(zhǎng)度為3 mm，壁厚為0.2 mm，外徑為1.6 mm，內(nèi)徑為1.2 mm。感壓管進(jìn)口處有45°的倒角，轉(zhuǎn)彎處的半徑為2 mm。

圖2 總壓探針模型

總壓探針置于葉柵通道出口，每個(gè)葉柵通道寬18°，葉片厚度為2.3°，因此葉柵通道的實(shí)際寬度為15.7°。再結(jié)合葉片尾緣與軸線的距離可知，葉柵通道實(shí)際寬度為15.7d(d為感壓管外徑)。如圖3所示，本文通過(guò)改變感壓管軸線與葉盆之間的距離D，研究不同安裝位置下，探針的測(cè)量準(zhǔn)確度以及探針對(duì)流場(chǎng)的影響。D的取值分別為1.5d,2.5d,3.5d,4.5d,5.5d,6.5d,7.5d,8.5d,9.5d,10.5d,12.5d,14.5d。

圖3 總壓探針在葉柵通道中的位置

1.2 網(wǎng)格劃分

利用ANSYS ICEM對(duì)計(jì)算模型進(jìn)行進(jìn)一步處理并劃分網(wǎng)格，整體采用非結(jié)構(gòu)網(wǎng)格，在探針附近進(jìn)行局部加密后設(shè)置1.2的增長(zhǎng)比，探針近壁面及計(jì)算域壁面處采用邊界層網(wǎng)格，保證近壁面y+為1左右，并且在近壁區(qū)內(nèi)設(shè)定10層邊界層網(wǎng)格。網(wǎng)格數(shù)量達(dá)到450萬(wàn)以后，進(jìn)出口總壓比基本穩(wěn)定，可認(rèn)為此時(shí)網(wǎng)格數(shù)量對(duì)計(jì)算結(jié)果的影響可忽略，將數(shù)量大于500萬(wàn)范圍的網(wǎng)格作為本文計(jì)算網(wǎng)格。數(shù)值計(jì)算的全部殘差曲線收斂至1×10-5以下。圖4表示網(wǎng)格劃分及邊界層處理。

圖4 網(wǎng)格劃分

1.3 邊界條件

數(shù)值計(jì)算采用Fluent18.0軟件完成，所有算例皆基于壓力求解器進(jìn)行穩(wěn)態(tài)流動(dòng)求解，計(jì)算域中流體選用理想氣體。入口選擇壓力入口，設(shè)定入口的總溫、總壓，流動(dòng)的方向；出口選擇壓力出口，流動(dòng)方向垂直于出口邊界，其余邊界皆設(shè)置為絕熱無(wú)滑移壁面。通過(guò)改變無(wú)探針時(shí)出口的背壓來(lái)實(shí)現(xiàn)有無(wú)探針時(shí)的進(jìn)口邊界條件一致。

表1 發(fā)動(dòng)機(jī)不同工作狀態(tài)氣動(dòng)參數(shù)表

表2 不同工作狀態(tài)進(jìn)口流量

1.4 參數(shù)定義

對(duì)于總壓探針，有總壓相對(duì)測(cè)量誤差w，表示同一位置，未放探針和放探針之后總壓值的對(duì)比，計(jì)算公式如下所示

(1)

除了對(duì)探針測(cè)量誤差的評(píng)價(jià)，本文還研究探針的放置對(duì)流場(chǎng)造成的影響，評(píng)價(jià)參數(shù)如下：

總壓損失系數(shù)

(2)

式中：pt,in表示進(jìn)口截面總壓；pt,i表示截面i的總壓平均值；ρ表示氣體密度；v表示氣流速度。

速度均勻度

(3)

2 結(jié)果與討論

2.1 不同安裝位置對(duì)探針測(cè)量結(jié)果的影響

本節(jié)具體分析不同質(zhì)量流量、不同安裝位置的探針對(duì)測(cè)量結(jié)果的影響。取感壓管彎管之前的1 mm直線上的總壓值平均作為總壓探針的測(cè)量值。取值線在感壓管的相對(duì)位置如圖5所示。

圖5 測(cè)量值取值線

圖6為總壓相對(duì)測(cè)量誤差隨安裝位置改變而變化的趨勢(shì)。圖中3條曲線分別為巡航、最大連續(xù)、起飛3個(gè)工作狀態(tài)下，探針總壓相對(duì)測(cè)量誤差隨安裝位置的變化規(guī)律。

圖6 不同安裝位置探針的總壓相對(duì)測(cè)量誤差

從圖中可以看出，當(dāng)D/d為6.5時(shí)，探針的總壓相對(duì)測(cè)量誤差最小。同時(shí)，比較不同工作狀態(tài)下，探針的總壓相對(duì)測(cè)量誤差，可以發(fā)現(xiàn)，隨著進(jìn)口流量增加，總壓相對(duì)測(cè)量誤差逐漸增加。這是因?yàn)殡S著流量增加，來(lái)流流速逐漸增加，探針對(duì)流場(chǎng)的擾動(dòng)相應(yīng)增加，在探針處的速度梯度逐漸增加，因此總壓相對(duì)測(cè)量誤差逐漸增加。當(dāng)D/d為6.5時(shí)，在巡航工作狀態(tài)下，總壓相對(duì)誤差為0.047%，轉(zhuǎn)變?yōu)樽畲筮B續(xù)工作狀態(tài)后，進(jìn)口流量相對(duì)增加14.6%，總壓相對(duì)測(cè)量誤差相對(duì)增加53%，流量進(jìn)一步增加，轉(zhuǎn)變?yōu)槠痫w工作狀態(tài)時(shí)，進(jìn)口流量相對(duì)增加4.6%，總壓相對(duì)測(cè)量誤差相對(duì)增加8.3%。分析總壓相對(duì)測(cè)量隨安裝位置變化而變化的規(guī)律可以發(fā)現(xiàn)，當(dāng)D/d較小時(shí)，即總壓探針位于葉盆尾緣附近，探針的總壓相對(duì)測(cè)量誤差非常大，隨著D/d逐漸增加，總壓相對(duì)測(cè)量誤差逐漸減小，直至D/d為6.5時(shí)，總壓相對(duì)測(cè)量誤差降到最低。當(dāng)安裝位置逐漸靠近葉背，可以發(fā)現(xiàn)總壓相對(duì)測(cè)量誤差又逐漸增加。

圖7是D/d為6.5時(shí)，探針感壓管速度云圖和流線圖。從圖7可知，在葉盆尾緣附近，存在一個(gè)速度梯度較大的區(qū)域，速度梯度對(duì)探針測(cè)量準(zhǔn)確度影響非常大，因此，當(dāng)總壓探針距離葉盆較近時(shí)，總壓相對(duì)測(cè)量誤差較大(圖6中D/d較小時(shí))，隨著探針遠(yuǎn)離葉盆附近的高壓力梯度區(qū)域，總壓相對(duì)誤差逐漸減小。若總壓探針與葉盆的距離進(jìn)一步增加，靠近葉背尾緣，則總壓探針將處于葉片尾跡中，因此總壓相對(duì)測(cè)量誤差同樣將逐漸增加，隨著總壓探針與葉盆距離再增加，探針則脫離葉片尾跡，因此總壓相對(duì)測(cè)量誤差減小(圖6中D/d較大時(shí))。

圖7 D/d=6.5時(shí)探針周圍速度分布(起飛狀態(tài))

圖8 不同安裝位置探針感壓管下游壓力分布(起飛狀態(tài))

圖8為起飛工作狀態(tài)下，當(dāng)D/d為2.5,6.5,14.5時(shí)，探針感壓管周圍總壓分布云圖。當(dāng)D/d為14.5時(shí)，如圖8a)所示，探針的感壓管正好在葉片的尾跡中，因此當(dāng)探針位于葉背下游時(shí)，相對(duì)測(cè)量誤差明顯比其他安裝位置探針的總壓測(cè)量誤差大；圖8b)表明，當(dāng)D/d為6.5時(shí)，總壓探針位于葉柵通道中間部位，葉柵通道中間部位壓力分布相對(duì)均勻，且速度同樣相對(duì)均勻，總壓探針沒(méi)有受到葉片尾跡的影響，此外，在感壓管外壁面附近壓力較低，氣流經(jīng)過(guò)感壓管之后，總壓損失也會(huì)相對(duì)較大；如圖8c)所示，當(dāng)D/d為2.5時(shí)，探針處于葉盆的出口的附近，相較于葉柵通道中間部位流場(chǎng)，此處氣流速度相對(duì)較小，但是感壓管距離葉片尾跡較近，因此總壓相對(duì)測(cè)量誤差也較大。

2.2 不同安裝位置的探針對(duì)流場(chǎng)的影響

為了研究探針下游流場(chǎng)的變化情況，在探針?biāo)诮孛嫒∫粭l沿流向的曲線，分析該曲線上總壓損失系數(shù)沿流向的變化趨勢(shì)。并在探針下游取了6個(gè)截面，分別距探針后壁面5,8,16,24,32 mm以及流域出口截面，同時(shí)考慮到探針對(duì)流場(chǎng)的影響有限，6個(gè)截面的周向角度覆蓋了探針安裝的葉柵通道(整體角度范圍為16°)，將這6個(gè)截面分別命名為y1～y6，如圖9所示。

圖9 探針下游曲線以及截面示意圖

圖10 總壓損失系數(shù)變化趨勢(shì)

圖10a)為最大連續(xù)工作狀態(tài)下，當(dāng)D/d為14.5時(shí)，總壓探針下游總壓損失系數(shù)沿流向變化趨勢(shì)。如圖10a)所示，總壓損失系數(shù)沿流程的整體變化趨勢(shì)：在總壓探針上游，總壓損失系數(shù)沿流程幾乎保持不變，數(shù)值接近0，至探針后壁面(流向156 mm)，總壓損失系數(shù)迅速上升至1.09，隨后，總壓損失系數(shù)沿流向迅速下降，直至探針下游28 mm，總壓損失系數(shù)約為6%，然后保持穩(wěn)定，直至探針下游出口。總壓損失系數(shù)保持在6%以下。這是因?yàn)榭倝禾结槍?duì)上游流場(chǎng)幾乎沒(méi)有影響，而對(duì)下游流場(chǎng)影響非常大。探針對(duì)流場(chǎng)的擾動(dòng)主要體現(xiàn)在下游，當(dāng)氣流經(jīng)過(guò)探針時(shí)，會(huì)造成局部總壓損失，探針正后方的總壓迅速減小，隨著流向距離逐漸增加，探針兩側(cè)氣流逐漸向探針后方匯集，因此氣流總壓逐漸增加，總壓損失系數(shù)逐漸減小，隨后，探針對(duì)流場(chǎng)的擾動(dòng)消失，總壓損失系數(shù)沿流向幾乎不變。圖10b)為最大連續(xù)工作狀態(tài)下，不同安裝位置的探針下游的總壓損失系數(shù)。由前文分析可知，探針對(duì)流場(chǎng)的擾動(dòng)主要體現(xiàn)在探針下游，因此圖10b)僅選取了探針下游的點(diǎn)，分析這些點(diǎn)的總壓損失系數(shù)變化規(guī)律。如圖10b)所示，在最大連續(xù)工作狀態(tài)下，當(dāng)D/d分別為6.5和2.5時(shí)，總壓探針下游總壓損失系數(shù)沿流程的變化曲線幾乎重合，而當(dāng)D/d為14.5時(shí)，在總壓探針下游4 mm(沿流程160 mm處)至探針下游28 mm，總壓損失系數(shù)略大于其他2處下游的總壓損失系數(shù)，但是，在探針下游28 mm之后，3個(gè)安裝位置探針下游的總壓損失系數(shù)相同，均在6%之下，直至流域出口。根據(jù)前文中葉柵通道速度云圖(見(jiàn)圖7)可知，當(dāng)安裝位置距離葉背很近的時(shí)候，探針將處于葉片的尾跡中，同時(shí)，此處的氣流速度相對(duì)較低，因此探針兩側(cè)氣流向探針后方匯集速度減緩，導(dǎo)致探針下游總壓損失系數(shù)在較長(zhǎng)的流向距離內(nèi)依然較大。

圖11為巡航工作狀態(tài)下，3個(gè)安裝位置探針下游的速度均勻度沿流向的變化趨勢(shì)。從圖中可以看出，不管是安裝探針還是未安裝探針的流場(chǎng)，在y1～y5截面，速度均勻度逐漸增加。這是因?yàn)閥1截面距離探針后壁面僅有5 mm，在y1截面由于受到葉片尾跡的影響，速度均勻度相對(duì)較低，而當(dāng)流向距離進(jìn)一步增加，葉片尾跡的影響消失，因此速度均勻度迅速增加。同時(shí)，在總壓探針下游24 mm之前，原始流場(chǎng)的速度均勻度略大于安裝探針流場(chǎng)相應(yīng)截面的速度均勻度，而在總壓探針下游24 mm之后，3個(gè)安裝位置探針下游的速度均勻度與原始流場(chǎng)相應(yīng)截面速度均勻度接近。從上文可知，在探針下游28 mm之前, 探針對(duì)流場(chǎng)的擾動(dòng)相對(duì)較大, 因此y1～y3截面的速度均勻度小于原始流場(chǎng)相應(yīng)截面的速度均勻度，在探針下游28 mm之后，探針對(duì)流場(chǎng)的影響逐漸消失，這里的現(xiàn)象再次驗(yàn)證了上文的結(jié)論。比較不同安裝位置探針下游y1截面到y(tǒng)3截面的速度均勻度，可以發(fā)現(xiàn)，當(dāng)安裝位置靠近葉柵通道中間位置時(shí)，探針下游的速度均勻度最大，當(dāng)安裝位置靠近葉盆，探針下游速度均勻度次之，當(dāng)安裝位置靠近葉背時(shí)葉片的尾跡和探針的尾跡疊加，對(duì)流場(chǎng)擾動(dòng)增加，因此探針下游速度均勻度最小。在探針下游5 mm，未放置探針時(shí)，流場(chǎng)速度均勻度為92.9%，當(dāng)安裝位置距離葉背6.5d時(shí)，速度均勻度減小0.14%，當(dāng)安裝位置在葉盆附近，D/d為2.5時(shí)，流場(chǎng)速度均勻度為92.6%，當(dāng)安裝位置在葉背附近，D/d為14.5時(shí)，流場(chǎng)速度均勻度減小了0.43%。

圖11 不同安裝位置探針下游速度均勻度

3 結(jié) 論

1) 在進(jìn)口級(jí)壓氣機(jī)出口截面，當(dāng)安裝位置距離葉盆5.5d到7.5d之間，探針的總壓相對(duì)測(cè)量誤差相對(duì)較小，且當(dāng)D/d為6.5時(shí)總壓相對(duì)測(cè)量誤差最小；

2) 隨著進(jìn)口流量增加，總壓探針的總壓相對(duì)測(cè)量誤差逐漸增加；在D/d為6.5的安裝位置下，探針的測(cè)量誤差與進(jìn)口流量呈線性關(guān)系；

3) 在3個(gè)典型的工作狀態(tài)下，總壓探針對(duì)流場(chǎng)的影響范圍有限，在探針下游28 mm之后，流場(chǎng)中心的總壓損失系數(shù)下降為6%；

4) 探針安裝位置越靠近葉背，探針下游的總壓損失系數(shù)下降速度越慢，探針在流向上的影響范圍更大。相較于葉柵通道中心的探針，靠近葉背位置的探針在流向上的影響范圍增加了16 mm；

5) 探針安裝位置越靠近葉背，探針下游的速度均勻度越小，相較于葉背位置，當(dāng)探針安裝位置在葉盆附近時(shí)，速度均勻度略微增加。