張 勃，吉洪湖，曹廣州，黃 偉

（南京航空航天大學(xué) 能源與動(dòng)力學(xué)院，南京 210016）

圓轉(zhuǎn)矩形噴管出口寬高比對(duì)射流雷諾剪應(yīng)力的影響

張 勃，吉洪湖，曹廣州，黃 偉

（南京航空航天大學(xué) 能源與動(dòng)力學(xué)院，南京 210016）

采用熱線(xiàn)風(fēng)速儀，利用單斜絲，對(duì)寬高比W／H分別為1，4，8，12，16的5個(gè)圓轉(zhuǎn)矩形收斂噴管和一個(gè)軸對(duì)稱(chēng)噴管的射流對(duì)稱(chēng)面上雷諾剪應(yīng)力分布特性進(jìn)行了實(shí)驗(yàn)研究。研究發(fā)現(xiàn)：在噴口下游不同截面上，射流寬、窄對(duì)稱(chēng)面上的雷諾剪應(yīng)力沿徑向均先緩慢增大，到達(dá)射流邊界后迅速減小，射流邊界逐漸沿徑向外移。矩形噴管射流相比軸對(duì)稱(chēng)射流具有較強(qiáng)的旋流，雷諾剪應(yīng)力較大，且隨著寬高比增大，旋流強(qiáng)度增大，剪應(yīng)力也逐漸提高，導(dǎo)致了射流與外流摻混增強(qiáng)。寬高比大于8以后，增大幅度逐漸減小。射流寬、窄對(duì)稱(chēng)面上的分布規(guī)律相同。

圓轉(zhuǎn)矩形噴管；射流；雷諾剪應(yīng)力；摻混；紅外抑制

0 引 言

飛行器發(fā)動(dòng)機(jī)的紅外輻射信號(hào)是敵方紅外制導(dǎo)導(dǎo)彈前向和側(cè)向攻擊時(shí)的主要搜索與跟蹤目標(biāo)。在眾多的輻射源中，噴氣發(fā)動(dòng)機(jī)的尾噴流紅外輻射是3～5μm波段上的主要輻射來(lái)源之一。二元非軸對(duì)稱(chēng)噴管可以強(qiáng)化尾噴流與外流的摻混，起到降低其紅外輻射的目的，在國(guó)外多種先進(jìn)飛機(jī)上得到實(shí)際應(yīng)用，取得了良好的效果［1］。

各國(guó)學(xué)者的研究顯示，非軸對(duì)稱(chēng)二元噴管加強(qiáng)摻混的機(jī)理在于其加強(qiáng)了射流剪切運(yùn)動(dòng)，具有較強(qiáng)的二次渦流，在其拐角區(qū)誘導(dǎo)出反向渦對(duì)，強(qiáng)化尾噴流與外界大氣的摻混，促進(jìn)其與中心區(qū)氣流的能量交換［2］。

射流強(qiáng)化摻混主要源于其與外流之間雷諾應(yīng)力的增大，A.Krothapalli［3］，Y.Tsuchiya［4］在低雷諾數(shù)條件下，分別對(duì)矩形與軸對(duì)稱(chēng)噴管射流流場(chǎng)進(jìn)行了實(shí)驗(yàn)研究，得到了射流雷諾剪應(yīng)力沿徑向先增大再減小，正應(yīng)力沿徑向逐漸減小的規(guī)律，D.R.Webster［5］則進(jìn)一步指出射流雷諾剪應(yīng)力大小取決于噴管形狀和氣流參數(shù)。Hinze［6］表明，射流與外流交界面附近的雷諾應(yīng)力大小是影響其摻混強(qiáng)弱的主要因素之一，且其噴口形狀變化對(duì)雷諾應(yīng)力具有明顯影響。

為了研究寬高比對(duì)矩形噴管射流摻混特性的影響規(guī)律，以一個(gè)軸對(duì)稱(chēng)噴管（AR0）為基礎(chǔ)，保持進(jìn)、出口面積不變，設(shè)計(jì)了出口寬高比為1，4，8，12，16（分別命名為AR1，AR4，AR8，AR12，AR16）的5個(gè)圓轉(zhuǎn)矩形收斂噴管，對(duì)其射流下游不同截面上的雷諾剪應(yīng)力分布特征進(jìn)行了實(shí)驗(yàn)研究。

1 實(shí)驗(yàn)設(shè)備與模型

采用的實(shí)驗(yàn)臺(tái)如圖1所示，主要包括收斂噴管、氣流穩(wěn)定段、風(fēng)機(jī)、變頻器和熱線(xiàn)風(fēng)速儀以及相關(guān)的數(shù)據(jù)采集設(shè)備。

圖1 實(shí)驗(yàn)臺(tái)架示意圖Fig.1 Sketch of the test rig

實(shí)驗(yàn)中，在風(fēng)機(jī)出口處布置了一段350mm的漸擴(kuò)段，在其與噴口之間布置一段與噴管進(jìn)口等徑的長(zhǎng)900mm的平直管段，使得流體到達(dá)噴管進(jìn)口時(shí)的速度分布達(dá)到管流充分發(fā)展時(shí)的速度分布。噴管進(jìn)、出口半徑分別為115mm、70mm，長(zhǎng)145mm。

圖2中以模型AR4為例，給出了噴管寬邊、窄邊、過(guò)渡面的形狀。

圖2 模型照片F(xiàn)ig.2 Photo of the model

圖3中給出了寬高比分別為1、4、8、12、16的5個(gè)矩形噴管模型與一個(gè)軸對(duì)稱(chēng)噴管模型照片，不同噴管的進(jìn)口面積均為0.04m2，出口面積均為0.02m2。

圖3 實(shí)驗(yàn)?zāi)Ｐ驼掌現(xiàn)ig.3 Photo of the test models

實(shí)驗(yàn)中使用變頻器對(duì)風(fēng)機(jī)轉(zhuǎn)速進(jìn)行調(diào)節(jié)以滿(mǎn)足實(shí)驗(yàn)雷諾數(shù)要求。氣流出口雷諾數(shù)均為2.5×105（Re＝UD／μ，U為氣流出口速度，D為試驗(yàn)中軸對(duì)稱(chēng)噴管AR0出口直徑），環(huán)境壓力為1.013×105Pa，環(huán)境溫度為303K。

2 測(cè)量參數(shù)與測(cè)量方法

為了研究射流的雷諾剪應(yīng)力ε＝ρuiuj，利用噴管出口平均速度的動(dòng)壓為參照，采用相對(duì)雷諾剪應(yīng)力ε為射流軸向X 速度脈動(dòng)，uj（j＝2，3）為側(cè)向Y，Z速度脈動(dòng)，U0是噴管出口平均速度），衡量射流的側(cè)向速度脈動(dòng)特征。

試驗(yàn)采用的斜絲為直徑3.8μm，長(zhǎng)3mm的45°鉑絲（如圖4中熱線(xiàn)支桿左端兩根金屬桿之間細(xì)絲），熱絲兩端鍍有合金，減少了桿繞流的干擾，有效提高測(cè)量精度，速度測(cè)量誤差約為1%。由于研究的矩形射流具有明顯的二維特征，采用將斜絲旋轉(zhuǎn)180°的方法測(cè)量求取側(cè)向速度脈動(dòng)，這種方法與采用叉絲測(cè)量相比，對(duì)流場(chǎng)影響較小，在流場(chǎng)相對(duì)穩(wěn)定的情況下，試驗(yàn)精度良好［6－9］。

圖4 45°斜絲Fig.4 Inclined hot wire（45°）

旋轉(zhuǎn)前后熱絲位置如圖5（a）和（b）所示。

在求取雷諾剪應(yīng)力時(shí)，對(duì)于圖5（a）中熱絲的45°放置位置，建立如下方程：

圖5 熱線(xiàn)的兩種放置位置Fig.5 Two locations of hot wire

式中，C′＝，U為平均速度，E為時(shí)均電壓，e則為脈動(dòng)電壓。

針對(duì)熱絲旋轉(zhuǎn)前后位置分別列出對(duì)應(yīng)的方程，對(duì)旋轉(zhuǎn)前后的方程進(jìn)一步取時(shí)均計(jì)算，可以得到：

實(shí)驗(yàn)選取射流的兩個(gè)特征面，即寬對(duì)稱(chēng)面與窄對(duì)稱(chēng)面，在其下游 X＝0.1D＊，0.5D＊，1D＊，2D＊，4D＊，8D＊（D＊＝4A／C，A 和C 分別為軸對(duì)稱(chēng)噴管AR0出口面積與濕邊周長(zhǎng)）6個(gè)不同軸向截面進(jìn)行測(cè)量，截面位置如圖6所示。

圖6 噴口流場(chǎng)下游測(cè)量平面圖Fig.6 Measurement planes of flow field at nozzle downstream

3 實(shí)驗(yàn)結(jié)果分析

圓到矩形過(guò)渡段的存在使得流體在流動(dòng)過(guò)程中存在一個(gè)從等徑的圓截面向不等徑的矩形出口流動(dòng)的過(guò)程，在寬邊與窄邊方向，噴管截面分別具有擴(kuò)張與收縮特征，導(dǎo)致了明顯旋流特征，射流在兩個(gè)方向上也分別表現(xiàn)出不同的流動(dòng)特征。

噴管寬高比的變化會(huì)使得射流出口旋流強(qiáng)度發(fā)生變化，從而引起速度脈動(dòng)的變化。下文針對(duì)不同寬高比條件下，噴管射流的雷諾剪應(yīng)力在寬、窄對(duì)稱(chēng)面上的分布規(guī)律進(jìn)行分析。

3.1 對(duì)稱(chēng)面上雷諾剪應(yīng)力分布

圖7中給出了軸對(duì)稱(chēng)噴管（AR0）射流在下游不同截面處雷諾剪應(yīng)力沿徑向分布。

圖7 AR0對(duì)稱(chēng)面上的雷諾剪應(yīng)力Fig.7 Reynolds shear stress at AR0symmetrical plane

可以看出，各個(gè)截面的雷諾剪應(yīng)力均從中心軸線(xiàn)沿徑向先增大后減小。以X／D＊＝0.1處的曲線(xiàn)為例，其中心線(xiàn)上的雷諾剪應(yīng)力最小，隨著徑向距離增大逐漸增大，在射流與外流交界面上達(dá)到最大，此時(shí)，射流與外流發(fā)生劇烈摻混，能量發(fā)生耗散，而后沿徑向應(yīng)力減小。射流在擴(kuò)張過(guò)程中，應(yīng)力沿軸向也逐漸減小，與董志勇［10］所述變化規(guī)律相同。

圖8 AR1對(duì)稱(chēng)面的雷諾剪應(yīng)力Fig.8 Reynolds shear stress at AR1symmetrical plane

不同軸向截面上雷諾剪應(yīng)力分布規(guī)律類(lèi)似，且隨著軸向距離的增大，其最大值均沿徑向逐漸外移，量值逐漸減小，說(shuō)明射流邊界逐漸擴(kuò)大。這個(gè)過(guò)程中，射流與外流發(fā)生能量交換，摻混能力沿程略有減低。在X／D＊＝0.1和8的兩個(gè)截面上，最大值出現(xiàn)的徑向位置從Y＝0.5D＊外移到Y(jié)＝2D＊附近，應(yīng)力值相應(yīng)從0.018減小為0.0038。

圖8給出了正方形噴管AR1對(duì)稱(chēng)面上的雷諾剪應(yīng)力分布，不同軸向距離處，各個(gè)截面上的雷諾剪應(yīng)力沿徑向先增大后減小。

圖9（a）與（b）分別給出了模型AR8寬、窄對(duì)稱(chēng)面上的雷諾剪應(yīng)力分布。

圖9 AR8對(duì)稱(chēng)面上的雷諾剪應(yīng)力Fig.9 Reynolds shear stress at AR8symmetrical plane

射流寬邊與窄邊上雷諾剪應(yīng)力量值相近，分布規(guī)律與AR1相同。

3.2 不同寬高比條件下，雷諾剪應(yīng)力的變化

為了研究寬高比變化對(duì)射流雷諾剪應(yīng)力的影響，將相同軸向截面上，不同噴管寬、窄對(duì)稱(chēng)面上的雷諾剪應(yīng)力分別繪制在兩張圖上，如圖10（a），（b）所示。

圖10 下游1D＊處對(duì)稱(chēng)面上雷諾剪應(yīng)力分布Fig.10 Reynolds shear stress distribution at symmetrical plane of 1D＊

圖10中給出噴口下游1D＊處，不同噴管射流雷諾剪應(yīng)力沿徑向的分布。寬對(duì)稱(chēng)面上（圖10（a）），噴管AR0中，雷諾剪應(yīng)力沿徑向從中心處的低值上升至射流邊界的0.001左右，正方形噴管AR1的變化規(guī)律與AR0相同，沿徑向先增大后減小。隨著噴管寬高比增大，矩形噴管角區(qū)渦對(duì)的發(fā)展逐漸增強(qiáng)，摻混能力強(qiáng)的特點(diǎn)逐漸顯示出來(lái)，AR4噴管中的雷諾剪應(yīng)力達(dá)到了0.12，寬高比為8的AR8噴管，雷諾剪應(yīng)力進(jìn)一步增大，這是由于矩形噴管出口存在不對(duì)稱(chēng)結(jié)構(gòu)，其對(duì)周?chē)黧w的卷吸作用引起的側(cè)向拉伸明顯強(qiáng)于軸對(duì)稱(chēng)噴管導(dǎo)致的。但是當(dāng)寬高比進(jìn)一步增大至12、16時(shí)，射流雷諾剪應(yīng)力增大速度略有下降。觀(guān)察發(fā)現(xiàn)，AR12與AR16噴管射流雷諾剪應(yīng)力沿徑向衰減較慢，距離較長(zhǎng)，這是由于射流湍動(dòng)能在射流邊界附近強(qiáng)度較大，而在其它位置較小導(dǎo)致的。

相比而言，窄對(duì)稱(chēng)面上（圖10（b）），不同寬高比噴管射流的衰減距離區(qū)別較小，這是由于在噴口面積一定的情況下，寬高比的增大意味著噴口高度的降低，一定程度上限制了出口流向渦在近噴口區(qū)域的發(fā)展。

圖11中給出了噴口下游4D＊處，寬、窄對(duì)稱(chēng)面上，不同噴管射流雷諾剪應(yīng)力沿徑向的分布。

圖11 下游4 D＊處對(duì)稱(chēng)面上雷諾剪應(yīng)力Fig.11 Reynolds shear stress distribution at symmetrical plane of 4D＊

與圖10中所示下游1D＊截面上的數(shù)據(jù)相比，不同寬高比噴管射流沿徑向的變化規(guī)律相同，均為先增大至射流邊界，而后迅速減小。但是不同噴管之間的區(qū)別明顯減小。這是由于圓轉(zhuǎn)矩形過(guò)渡段對(duì)射流的影響主要表現(xiàn)在近噴口區(qū)，此處，其誘導(dǎo)出的軸向與角區(qū)渦對(duì)強(qiáng)度較大。隨著軸向距離增大，噴口形狀變化的影響逐漸減弱，渦對(duì)強(qiáng)度減小，不同噴管射流之間應(yīng)力數(shù)值接近。

Rajaratnam［11］，Sforza［12］對(duì)直噴管射流的研究中也得到了軸向距離增大，噴管出口影響減弱的結(jié)論。

4 結(jié) 論

利用熱線(xiàn)風(fēng)速儀，在低雷諾數(shù)條件下，對(duì)圓轉(zhuǎn)矩形噴管射流下游的雷諾剪應(yīng)力分布特征以及寬高比的影響進(jìn)行了實(shí)驗(yàn)研究，得出如下結(jié)論：

（1）在射流下游不同軸向截面上，雷諾剪應(yīng)力沿徑向先緩慢增大，到射流邊界附近與外流摻混加劇，剪應(yīng)力迅速減小，射流邊界沿軸向逐漸外移；

（2）矩形噴管誘導(dǎo)出較強(qiáng)的出口旋流，雷諾剪應(yīng)力大于軸對(duì)稱(chēng)噴管，并且隨著寬高比增大逐漸提高，當(dāng)寬高比大于8以后，噴口高度對(duì)旋流形成壓制，其值則略有減?。?/p>

（3）在亞聲速射流中，噴管出口形狀、出口寬高比的不同對(duì)射流雷諾剪應(yīng)力影響較大，隨著射流軸向距離的增大，噴管不對(duì)稱(chēng)結(jié)構(gòu)誘導(dǎo)出的渦系逐漸減弱，使得對(duì)雷諾剪應(yīng)力影響逐漸減小。

［1］理查森D.現(xiàn)代隱身飛機(jī)［M］.北京：科學(xué)出版社，1991.

［2］HERMANN S.Boundary－layer theory［M］.NewYork：McGraw－Hill，1968.

［3］KROTHAPALLI A，BAGANOFF D.On the mixing of a rectangular jet［J］.J.Fluid Mech.，1981，107：201－220.

［4］TSUCHIYA Y，HORIKOSHI C.On the spreading of rectangular jets［J］.Experiments in Fluid，1986，16（4）：197－204.

［5］WEBSTER D R，LIU Y.Velocity measurement of turbulence callapse in a linearly satisfied jet［J］.Experiments in Fluid，2001，31：394－400.

［6］HINZE J O.Turbulence［M］.New York：McGraw－Hill，1975.

［7］田杰，陳旭，郝輝，等.旋轉(zhuǎn)單斜絲熱線(xiàn)測(cè)量葉輪機(jī)出口流場(chǎng)［J］.流體機(jī)械，2003，31（7）：25－28.

［8］MULLER O R.On the accuracy of turbtdenee measurements with inclined hot wires［J］.J.Fluid Mech.，1982，119：155－172.

［9］RUSS S，SIMON T W.On the rotation，slanted，hotwire technique［J］.Experiments in Fluids，1991，（12）：76－80.

［10］董志勇.射流力學(xué)［M］.北京：科學(xué)出版社，2005.

［11］RAJARATNAM N.Turbulent jets［M］.Amsterdam：Elsevier Scientific Publishing Company.

［12］SFORZA P M，STEIGER M H，TRENTACOSTE N.Studies on three－dimensional viscous jets［J］.AIAA J，1966，5（4）：800－806.

張 勃（1976－），男，江蘇南京人，副教授，博士。研究方向：強(qiáng)化傳熱，流動(dòng)密封。通訊地址：南京航空航天大學(xué)能源與動(dòng)力學(xué)院（210016），聯(lián)系電話(huà)：13601580069，E－mail：zhangbo＿pe＠nuaa.edu.cn

Reynolds shear stress characteristics of jets from circular－rectangular transition nozzles with different aspect ratios

ZHANG Bo，JI Hong－h(huán)u，CAO Guang－zhou，HUANG Wei
（Nanjing University of Aeronautics and Astronautics，Nanjing 210016，China）

The flow field and Reynolds shear stress characteristics of the jets from 5circularrectangular transition nozzles and one axisymmetric nozzle were investigated with inclined hot wire.The distributions of Reynolds shear stress in the wide and narrow symmetrical planes were measured on different downstream sections.Results show that in the radial direction，Reynolds shear stress increases slowly till the interface between jet and ambient gas，then decreases fast.Rectangular nozzles can induce stronger swirl flow than axisymmetric one and larger Reynolds shear stress，which increases with increasing aspect ratio.The Reynolds shear stress increment gradually decreases when the ratio exceeds 8.The distributions are the same in both wide and narrow symmetrical planes.

circular－rectangular transition nozzle；jet；Reynolds shear stress；mixing；infrared suppress

V231.1

A

1672－9897（2012）05－0042－05

2011－09－01；

2012－03－12