半封閉圓管沖擊射流流動特性PIV試驗研究

陳 明，黃海津，王多銀，呂森鵬，陳沿吉

(1.重慶交通大學(xué) 水利水運工程教育部重點實驗室，重慶 400074；2.重慶交通大學(xué) 國家內(nèi)河航道整治工程技術(shù)研究中心，重慶 400074)

長期以來沖擊射流被廣泛應(yīng)用于諸多工程領(lǐng)域[1-2]，實際上沖擊射流不僅包含自由沖擊射流，還常常遇到限域空間中的沖擊射流，如船閘明溝消能結(jié)構(gòu)、燃?xì)鈾C(jī)渦輪冷卻[3]、電子封裝[4]等。該類射流的限域空間主要由噴嘴出口處封閉板與沖擊板所形成，通常被稱為半封閉沖擊射流。如圖1所示，半封閉沖擊射流的流場結(jié)構(gòu)由四個具有明顯流動特征的區(qū)域組成：流體從孔口噴出形成自由射流，勢流核周圍自由剪切層卷吸環(huán)境流體，射流沿程擴(kuò)展；當(dāng)射流接近沖擊板時，軸向速度快速衰減，徑向速度增大，主流轉(zhuǎn)向；主流貼著沖擊板流出，并向軸向擴(kuò)展，形成壁面射流；壁面射流與封閉板相互作用形成回流。由于封閉板的存在，局部能量傳遞和耗散特性將受到影響，而所涉的能量傳遞和耗散特性又與流場特性密切相關(guān)。因此，研究半封閉沖擊射流流動特性具有重要意義。

圖1 半封閉沖擊射流示意圖

現(xiàn)有對半封閉沖擊射流的研究多側(cè)重于其傳熱特性[5]，對流動特性的關(guān)注還相對較少。根據(jù)噴嘴的類型，半封閉沖擊射流主要包含半封閉孔板沖擊射流和半封閉圓管沖擊射流。其中，孔板射流因受出口收縮效應(yīng)的影響，相較圓管射流將產(chǎn)生更薄的剪切層，從而增強(qiáng)方位角不穩(wěn)定性[6]，進(jìn)而導(dǎo)致半封閉孔板沖擊射流與半封閉圓管沖擊射流流動特性的差異。在半封閉孔板沖擊射流研究方面，基于流動可視化技術(shù)，Garimella等[7]發(fā)現(xiàn)存在一個由壁面射流和封閉板共同驅(qū)動的回流區(qū)；Fitzgerald等[8]采用激光多普勒測速技術(shù)(laser doppler velocimetry, LDV)對其流場分布進(jìn)行了詳細(xì)測量，研究了噴嘴直徑、雷諾數(shù)和沖擊距對回流區(qū)、速度分布及紊動強(qiáng)度的影響；Guo等采用粒子測速技術(shù)(particle image velocimetry, PIV)系統(tǒng)研究了不同雷諾數(shù)和沖擊距下回流區(qū)的運動規(guī)律、射流撞擊沖擊板前的摻混特性、壁面射流區(qū)的擴(kuò)展率、速度衰減及自相似性，并與自由沖擊射流進(jìn)行了比較。在半封閉圓管沖擊射流研究方面，Ashforth-Frost等[9]較早地采用LDV測量了其速度、紊動強(qiáng)度及紊動能的分布，但其主要目的是為了給數(shù)值計算提供驗證數(shù)據(jù)；徐驚雷等[10-11]采用熱線風(fēng)速儀對不同雷諾數(shù)和沖擊距下的流場進(jìn)行了測量，但只分析了雷諾數(shù)、沖擊距對主流速度分布和徑向紊動能分布的影響；Baydar等[12]研究了沖擊板上的壓力分布，并與自由沖擊射流進(jìn)行了比較，指出封閉板的影響導(dǎo)致在沖擊板上產(chǎn)生負(fù)壓區(qū)，且該負(fù)壓區(qū)與紊動強(qiáng)度、傳熱系數(shù)峰值存在一定聯(lián)系；Gao等[13]研究了雷諾數(shù)和沖擊距對傳熱特性的影響，并分析了封閉板對沖擊板上靜壓分布及脈動壓力變化的影響。上述兩項研究主要是針對封閉板對沖擊板上的壓力分布影響展開研究，尚未涉及封閉板對流場特性的影響。

綜上所述，與半封閉孔板沖擊射流相比，針對半封閉圓管沖擊射流流動特性的研究還存在些許滯后，主要表現(xiàn)在：① 僅停留在速度、紊動特性分布的描述，缺乏對流場勢流核長度、擴(kuò)展率、速度衰減指數(shù)等重要特征量的統(tǒng)計，而這些特征量通常可直接用來指導(dǎo)工程設(shè)計，如船閘的明溝寬度通常根據(jù)勢流核長度確定[14]；② 與自由沖擊射流的區(qū)別特征尚未進(jìn)行比較，而封閉板的影響將導(dǎo)致射流穩(wěn)定性及與周圍流體的摻混特征發(fā)生變化，從而引起流動特性的改變。因此，針對半封閉圓管沖擊射流的流動特性有待進(jìn)行更深入的研究。

鑒于此，本文采用PIV技術(shù)對不同雷諾數(shù)和沖擊距下的半封閉圓管沖擊射流進(jìn)行了詳細(xì)測量，基于時均流場著重分析了半封閉圓管沖擊射流的統(tǒng)計特性，包括回流區(qū)運動、射流核長度、射流沖擊前卷吸特性、壁面射流區(qū)擴(kuò)展率、速度衰減及自相似性。同時，將壁面射流區(qū)的擴(kuò)展率和速度衰減指數(shù)與半封閉孔板沖擊射流和自由沖擊射流比較，旨在揭示圓管沖擊射流在半封閉空間中的流動特性，為后續(xù)解釋該邊界條件下的射流能量傳遞和耗散機(jī)制奠定基礎(chǔ)。

1 試驗系統(tǒng)與測量方案

1.1 射流裝置

射流試驗在定制的水箱中進(jìn)行，如圖2所示，主要包括射流水箱、回水箱、變頻泵、電磁流量計及一些管路組成。射流水箱由超白玻璃制作而成，試驗段為230 mm×175 mm×300 mm(長×寬×高)，水箱四個角落的溢流板用于控制淹沒水深。射流由一根聚氯乙烯圓管中噴出，圓管內(nèi)徑d=16 mm，外徑20 mm，長度960 mm(60d)，該長徑比是為了保證射流出口充分發(fā)展。射流管布置在水箱的平分線上，其中心線距底板4d，并與底板平行，射流由變頻泵提供，出口流速由電磁流量計和變頻泵共同控制。為避免淺水效應(yīng)引起水面波動對PIV測量效果的影響，通過反復(fù)調(diào)試，確定整個試驗過程中淹沒水深恒定為270 mm。

(a)裝置示意圖

1.2 PIV測量系統(tǒng)

采用的PIV測量系統(tǒng)是由北京江宜科技有限公司研發(fā)的TR-PIV系統(tǒng)，配有大光圈定焦鏡頭(Canon EF50 mm f/1.0L)、CMOS相機(jī)(IDT NX5-S2)、8W連續(xù)激光器及PIV流場計算軟件。CMOS相機(jī)最高分辨率2 560×1 920像素，滿幅采樣頻率730 Hz。激光束經(jīng)3片柱面鏡和一片鮑威爾棱鏡后轉(zhuǎn)變?yōu)楹穸? mm的矩形片光[15]，由射流水箱側(cè)面打入，照亮射流孔中心水平面上的示蹤粒子。示蹤粒子采用平均直徑為10 μm、密度為1.03 g/mm3的聚酰胺微珠。PIV流場計算采用基于圖像變形的三重網(wǎng)格迭代法，第一級判讀窗口為64×64像素，第三級(最后一級)判讀窗口為16×16像素，相鄰判讀窗口之間重疊50%，最終得到空間分辨率為8×8像素的流場，相鄰矢量之間的間隔為0.55 mm。在圖像分析時，相關(guān)峰的位置采用三點高斯擬合進(jìn)行亞像素插值，計算精度可達(dá)0.1像素[16]。

1.3 試驗工況

試驗工況的選取，重點從以下幾個方面考慮：① 射流出口的速度不受沖擊板的影響，需保證H/d>1.5[8]；② 當(dāng)H/d>6時，封閉板對長直圓管沖擊射流的流動特性基本無影響；③ 沖擊距與雷諾數(shù)的組合是工程應(yīng)用中常遇流動型態(tài)的代表[17]。基于上述考慮，本文設(shè)計了3組初始速度、3組沖擊距共9組試驗工況(表1)。采用長時間低頻測量的方式采集流場數(shù)據(jù)，即每采集2幀粒子圖像(一個瞬時流場)后，先儲存到服務(wù)器硬盤，再采集下2幀粒子圖像，粒子圖像的采樣間隔為1.25 ms，瞬時流場間的平均間隔約為1 s。為滿足流場統(tǒng)計分析要求，各工況瞬時流場樣本容量均為5 000次[18]。

表1 試驗條件

2 結(jié)果與討論

2.1 回流區(qū)運動特性

回流區(qū)是半封閉沖擊射流區(qū)別自由沖擊射流最明顯的特征。為研究雷諾數(shù)對回流區(qū)的影響，圖3選取了典型沖擊距(H/d=2)條件下不同雷諾數(shù)時均流線分布圖，其中云圖為無量綱軸向速度。可見，當(dāng)Re=5 084、8 179時，形成了明顯的回流結(jié)構(gòu)，并隨著雷諾數(shù)增大，回流區(qū)中心由r/d=5.7下移到r/d=6.5(圖中虛線標(biāo)識位置)；當(dāng)Re=13 926時，回流區(qū)中心已不在測量范圍內(nèi)。由此可得，隨著雷諾數(shù)增加，回流區(qū)向出流方向移動，這與半封閉孔板沖擊射流的研究結(jié)果一致。

為研究沖擊距對回流區(qū)的影響，圖4選取了典型雷諾數(shù)(Re=5 084)條件下不同沖擊距時均流線分布圖。受激光寬度及相機(jī)拍攝范圍的限制，當(dāng)H/d=6時，僅給出了距沖擊板4.8d的拍攝范圍。結(jié)合圖3(a)，當(dāng)H/d=2時，測量范圍內(nèi)形成明顯回流結(jié)構(gòu)，當(dāng)H/d=4和H/d=6時，測量范圍內(nèi)均無明顯回流結(jié)構(gòu)。Guo等指出，隨著沖擊距增加，回流區(qū)將進(jìn)一步向出流方向移動，導(dǎo)致其位于測量區(qū)域以外。由此可以推測，隨著沖擊距增加，與半封閉孔板沖擊射流特征相同，回流區(qū)也向出流方向移動。

(a)Re=5 084

(a)H/d=4

2.2 射流沖擊前流動特性

封閉板的存在將對射流初期的勢流核長度[19]、卷吸特性產(chǎn)生影響，為了量化壁面射流形成的初始條件，利用時均流場分別計算了勢流核長度(l)和卷吸系數(shù)(c)。

勢流核長度按“95%準(zhǔn)則[20]”確定，即從孔口至um衰減到0.95uj點的軸向長度，uj為噴嘴孔口速度。圖5統(tǒng)計了各試驗工況勢流核長度值，圖中虛線標(biāo)識為沖擊板的位置。可見，相同雷諾數(shù)下，隨著沖擊距增大，勢流核長度增大，勢流核末端距沖擊板的距離增大；相同沖擊距下，雷諾數(shù)對勢流核長度影響不大。其中，H/d=2時，隨著雷諾數(shù)增大，勢流核長度幾乎保持恒定，長度值與Guo等的統(tǒng)計結(jié)果基本一致；但當(dāng)沖擊距增大至H/d=4時，本文統(tǒng)計的勢流核長度明顯大于Guo等的統(tǒng)計結(jié)果。分析原因，H/d=2時，狹窄的限制空間使射流剪切層還未觸及勢流核中心便開始快速衰減了，因此勢流核長度基本一致；而當(dāng)H/d=4時，勢流核長度與噴嘴出口速度分布有關(guān)，充分發(fā)展管道流速度分布(本文)產(chǎn)生的勢流核長度大于頂帽型速度分布(Guo等)產(chǎn)生的勢流核長度。

圖5 勢流核長度統(tǒng)計圖

卷吸系數(shù)采用Zhang等[21]提出的公式計算：

(1)

圖6統(tǒng)計了各試驗工況卷吸系數(shù)擴(kuò)展率α，統(tǒng)計范圍為0.2H～0.8H。可見，相同雷諾數(shù)下，隨著沖擊距增大，卷吸系數(shù)擴(kuò)展率增大，這是由于隨著射流的發(fā)展，射流斷面沿程擴(kuò)大，并持續(xù)卷吸周圍流體，使斷面質(zhì)量流量沿程增大，從而引起大沖擊距時更大的擴(kuò)展率。相同沖擊距下，隨著雷諾數(shù)增大，卷吸系數(shù)擴(kuò)展率基本上呈略微增大趨勢，這是由于隨著雷諾數(shù)增大，射流剪切層紊動加強(qiáng)，卷吸更多的周圍流體，從而引起卷吸系數(shù)擴(kuò)展率增大。然而，當(dāng)沖擊距增大至H/d=6時，雷諾數(shù)從5 084增大到8 179，卷吸系數(shù)擴(kuò)展率明顯減小。Guo等在分析半封閉孔板沖擊射流的擴(kuò)展角度時也發(fā)現(xiàn)了類似現(xiàn)象，他們發(fā)現(xiàn)H/d=8時，雷諾數(shù)從5 000增大到7 500，擴(kuò)展角度明顯減小。針對上述卷吸系數(shù)擴(kuò)展率或射流擴(kuò)展角隨雷諾數(shù)增大而減小的現(xiàn)象，分析認(rèn)為這可能與封閉板引起的開爾文-亥姆霍茲渦(Kelvin-Helmholtz vortices)不穩(wěn)定性有關(guān)。

圖6 卷吸系數(shù)擴(kuò)展率統(tǒng)計圖

2.3 壁面射流區(qū)流動特性

射流撞壁后會在壁面附近誘導(dǎo)形成壁面射流，壁面射流沿徑向發(fā)展將經(jīng)歷形成階段和充分發(fā)展階段。在充分發(fā)展階段壁面射流將達(dá)到自相似性，此時擴(kuò)展率按線性規(guī)律增長[22-23]、速度按冪指數(shù)規(guī)律衰減[24-26]。為研究封閉板對壁面射流的影響，本節(jié)將主要分析各雷諾數(shù)、沖擊距下壁面射流的擴(kuò)展率和速度衰減指數(shù)，并與Guo等半封閉孔板沖擊射流的研究結(jié)果和自由沖擊射流的統(tǒng)計結(jié)果進(jìn)行比較。

圖7(a)～圖7(c)分別給出了各試驗工況下，壁面射流半值寬z0.5隨徑向位置的變化，z0.5定義為徑向速度等于vm/2時距沖擊板的距離。結(jié)果顯示：① 從r/d=0開始，壁面射流先后經(jīng)歷了形成階段和充分發(fā)展階段，兩階段的半值寬具有各自的特征。② 在形成階段，半值寬隨徑向距離增加略微減小；在充分發(fā)展階段，半值寬隨徑向距離增加而線性增長，表明壁面射流進(jìn)入自相似區(qū)。③ 當(dāng)H/d=2時(圖7(a))，雷諾數(shù)從5 084增大到13 926，z0.5/d滿足線性增長的區(qū)域分別為1.9

(a)H/d=2

(2)

圖8統(tǒng)計了各試驗工況下壁面射流擴(kuò)展率β，可見，相同沖擊距下，隨著雷諾數(shù)增加，擴(kuò)展率減小，Van hout等[17]在研究自由沖擊射流時也發(fā)現(xiàn)相同的規(guī)律，并指出這可能與較低雷諾數(shù)時(Re<23 000)初生渦、次生渦的產(chǎn)生及相互作用有關(guān)。相同雷諾數(shù)下，隨著沖擊距增大，擴(kuò)展率整體上呈增長趨勢，與Cooper等、Knowles等在研究自由沖擊射流時得出的規(guī)律相同，這是由于隨著沖擊距增大，射流在撞擊沖擊板前產(chǎn)生更多的摻混，從而引起壁面射流區(qū)更大的擴(kuò)展率。但當(dāng)Re=5 084、8 179時，H/d=4的擴(kuò)展率明顯小于H/d=2的擴(kuò)展率，這是由于回流區(qū)的影響導(dǎo)致的，隨著回流區(qū)移出測量范圍(Re=13 926)，擴(kuò)展率恢復(fù)為隨沖擊距增大而增大。

圖8 壁面射流擴(kuò)展率統(tǒng)計圖

圖9(a)～圖9(c)分別給出了各試驗工況下，壁面射流最大徑向速度vm隨徑向位置的變化。結(jié)果顯示：① 最大徑向速度從r/d=0開始增大，到r/d=1.0左右達(dá)到最大值，而后隨著r/d的增大而減小，這表明最大徑向速度沿徑向位置存在一個由加速運動轉(zhuǎn)變?yōu)闇p速流動的過程，這與徑向的壓力梯度變化有關(guān)[10-11]。② 當(dāng)H/d=2，Re=13 926和H/d=4，Re=13 926時，最大徑向速度沿徑向位置存在兩個明顯的峰值，第二個峰值是由壁面流動分離引起的[27]。③ 當(dāng)H/d=2時(圖9(a))，最大徑向速度滿足冪指數(shù)衰減的區(qū)域與雷諾數(shù)有關(guān)，雷諾數(shù)從5 084增大到13 926，vm/uj滿足冪指數(shù)衰減的區(qū)域分別為1.9

(a)H/d=2

(3)

圖10統(tǒng)計了各試驗工況下壁面射流衰減指數(shù)γ，可見，相同沖擊距下，當(dāng)H/d=2時，隨著雷諾數(shù)增大，衰減指數(shù)增大；當(dāng)H/d≥4時，衰減指數(shù)與雷諾數(shù)無關(guān)，這與Xu等在研究自由沖擊射流時得出的規(guī)律相同。相同雷諾數(shù)下，當(dāng)H/d≥4時，隨著沖擊距增大，衰減指數(shù)降低。回流區(qū)也會對衰減指數(shù)隨沖擊距變化的一般規(guī)律產(chǎn)生影響，如當(dāng)Re=5 084、8 179時，H/d=2的衰減指數(shù)明顯小于H/d=6的衰減指數(shù)。圖中也給出了Guo等半封閉孔板沖擊射流的試驗結(jié)果，可見，① 兩者的絕對值存在差異；② 相較本文結(jié)果，Guo等試驗結(jié)果隨雷諾數(shù)或沖擊距變化的一致性均較差，這可能與孔板噴嘴產(chǎn)生的射流具有更強(qiáng)的不穩(wěn)定性及三維性有關(guān)。

圖10 壁面射流衰減指數(shù)統(tǒng)計圖

為進(jìn)一步證實回流區(qū)對壁面射流自相似性的破壞，分析了壁面射流無量綱徑向速度分布隨徑向位置的變化。試驗數(shù)據(jù)顯示，當(dāng)H/d=6時，不同雷諾數(shù)下壁面射流充分發(fā)展區(qū)徑向速度分布自相似性與H/d=4較為一致，為節(jié)省篇幅，圖11僅給出了沖擊距H/d=2、H/d=4時不同雷諾數(shù)下壁面射流無量綱徑向速度分布圖。由圖可見：① 各工況下徑向速度沿壁法向均呈現(xiàn)先增大再減小的壁面射流特性。② 當(dāng)H/d=2時，在半值寬線性增長區(qū)內(nèi)、最大徑向速度冪指數(shù)衰減區(qū)內(nèi)，無量綱速度分布均表現(xiàn)出較強(qiáng)的自相似性，而超出該區(qū)域后，無量綱速度分布開始不具備自相似性(圖11(a)，r/d=6.0；圖11(b)，r/d=6.5；圖11(c)，r/d=7.0)，由此可以證實回流區(qū)破壞了壁面射流的自相似性，從而導(dǎo)致半值寬線性增長規(guī)律和最大徑向速度冪指數(shù)衰減規(guī)律終止。③ 當(dāng)H/d=4時，壁面射流充分發(fā)展區(qū)各斷面無量綱速度分布均表現(xiàn)出較強(qiáng)的自相似性(圖11(d)～(f))，與擴(kuò)展率分析、速度衰減分析表現(xiàn)的規(guī)律一致。

(a)H/d=2，Re=5 084

表2統(tǒng)計了近年來自由沖擊射流在壁面射流區(qū)內(nèi)的擴(kuò)展率和衰減指數(shù)。同時，為了方便比較，用虛線將擴(kuò)展率、衰減指數(shù)的范圍分別描繪于圖8、圖10中。可見，半封閉圓管沖擊射流壁面射流區(qū)擴(kuò)展率基本位于自由沖擊射流統(tǒng)計的范圍內(nèi)，速度衰減指數(shù)顯著小于自由沖擊射流統(tǒng)計的范圍，這表明封閉板的存在將加速圓管沖擊射流壁面射流區(qū)速度衰減，而對圓管沖擊射流壁面射流區(qū)擴(kuò)展率的影響不顯著。該結(jié)果與Guo等的研究結(jié)果存在顯著差異，他們在研究半封閉孔板沖擊射流時得出，與自由沖擊射流相比，封閉板的存在將使孔板沖擊射流壁面射流區(qū)擴(kuò)展率增大，而對孔板沖擊射流壁面射流區(qū)速度衰減的影響不顯著。

表2 沖擊射流壁面射流區(qū)擴(kuò)展率、衰減指數(shù)統(tǒng)計

3 結(jié) 論

本文采用粒子圖像測速技術(shù)研究了半封閉圓管沖擊射流的時均流動特性，得到了不同雷諾數(shù)Re(5 084～13 926)和沖擊距下H/d(2～6)回流區(qū)的運動規(guī)律、射流沖擊前及壁面射流區(qū)的流動特性，并與Guo等半封閉孔板沖擊射流的研究結(jié)果和自由沖擊射流的統(tǒng)計結(jié)果進(jìn)行了比較。主要結(jié)論如下：

(1)對回流區(qū)的運動規(guī)律研究表明：雷諾數(shù)和沖擊距對回流區(qū)運動特征的影響規(guī)律與半封閉孔板沖擊射流一致，即回流區(qū)隨雷諾數(shù)或沖擊距的增大逐漸向出流方向移動。

(2)對射流沖擊前的流動特性研究表明：相同雷諾數(shù)下，隨著沖擊距增大，勢流核長度和卷吸系數(shù)擴(kuò)展率均增大；相同沖擊距下，雷諾數(shù)對勢流核長度影響不大，而隨著雷諾數(shù)增大，卷吸系數(shù)擴(kuò)展率基本上呈略微增大趨勢。

(3)對壁面射流區(qū)的流動特性研究表明：相同雷諾數(shù)下，隨著沖擊距增大，擴(kuò)展率增大(H/d≥4)，衰減指數(shù)減小(H/d≥4)；相同沖擊距下，隨著雷諾數(shù)增大，擴(kuò)展率減小，衰減指數(shù)與雷諾數(shù)無關(guān)(H/d≥4)；較小沖擊距(H/d=2)時，回流區(qū)將會破壞壁面射流自相似性，從而對擴(kuò)展率和衰減指數(shù)產(chǎn)生影響。

(4)與自由沖擊射流比較，封閉板的存在使得壁面射流區(qū)速度衰減加快，而對壁面射流區(qū)擴(kuò)展率影響不顯著，該結(jié)果與Guo等關(guān)于封閉板對孔板沖擊射流影響特性的研究成果存在顯著差異。

需要指出的是，本文僅給出了統(tǒng)計意義上半封閉圓管沖擊射流的時均流動特性。實際情況下，半封閉圓管沖擊射流的瞬態(tài)運動非常復(fù)雜，渦系演化迅速，現(xiàn)階段尚未捕捉到其瞬態(tài)運動特性，后續(xù)將在高頻采樣的條件下，結(jié)合數(shù)值模擬手段，系統(tǒng)研究半封閉圓管沖擊射流的瞬態(tài)運動規(guī)律。