環(huán)形水射流流場的實(shí)驗(yàn)研究與統(tǒng)計(jì)分析

康 燦,張 峰, 楊敏官, 肖勝男

(江蘇大學(xué)能源與動(dòng)力工程學(xué)院,江蘇鎮(zhèn)江 212013)

0 引 言

環(huán)形射流是一種新型射流,多被應(yīng)用在航空航天和動(dòng)力工程領(lǐng)域,近年來環(huán)形射流在工業(yè)清洗領(lǐng)域內(nèi)被嘗試并取得了較好的應(yīng)用效果。在三維空間,液流自環(huán)形噴嘴噴出后首先形成具有一定厚度的圓環(huán),在離開噴嘴一定距離后,圓環(huán)開始沿徑向擴(kuò)散,圓環(huán)的中心由空心變?yōu)槌錆M著液體[1]。環(huán)形液體射流中存在著射流穩(wěn)定性、射流束聚并、液滴撞擊和液滴破碎等難題[2-3]。從目前的數(shù)值模擬水平來看,高速環(huán)形射流場的高精度數(shù)值模擬和多尺度流動(dòng)結(jié)構(gòu)的求解均具有較大的難度[4]。而目前針對環(huán)形射流的定量實(shí)驗(yàn)研究報(bào)道也較少。

當(dāng)射流壓力達(dá)到10MPa時(shí),液體的速度一般會達(dá)到100m/s[5]。在這類較高流速的射流場中,能量分布是評價(jià)射流場品質(zhì)的重要指標(biāo),更是決定射流作業(yè)能力的關(guān)鍵因素,這一點(diǎn)不同于低壓圓柱射流和雙流體霧化[6]。同時(shí),流場中的液滴尺寸分布對射流流體施加到被作用物體表面的打擊作用有著重要的影響[7-8]。

對一環(huán)形噴嘴的自由射流流場進(jìn)行實(shí)驗(yàn)研究,設(shè)定3種射流壓力p=11,12和15MPa,對不同射流壓力下的射流流場進(jìn)行PDPA光學(xué)測量。測量流場中的單點(diǎn)平均速度和脈動(dòng)速度,以說明射流能量的沿程衰減和流場的湍動(dòng)性能。通過對液滴粒徑分布的測量來衡量該流場內(nèi)的液滴破碎機(jī)制。對通過特定位置控制體的單個(gè)液滴粒徑和速度進(jìn)行統(tǒng)計(jì)分析,以描述群體液滴表現(xiàn)出的統(tǒng)計(jì)規(guī)律。

1 實(shí)驗(yàn)系統(tǒng)

采用的實(shí)驗(yàn)系統(tǒng)如圖1所示,水的壓能由功率為9kW的三缸柱塞泵提供,柱塞泵出口和噴嘴通過高壓軟管相連,高壓軟管內(nèi)的能量損失由柱塞泵出口和噴嘴入口的壓差反映。實(shí)驗(yàn)中采用的環(huán)形噴嘴如圖2所示,其出口直徑為1.68mm,中心圓柱的直徑為1.0mm,中心圓柱的下游端距離噴嘴出口斷面的軸向距離為1mm,采用該縮進(jìn)形式的目的是為了縮短圓環(huán)中心充滿液流所需的軸向距離。射流出口的環(huán)境壓力為當(dāng)?shù)卮髿鈮?。在距噴嘴出口不同軸向距離處選取6個(gè)測量斷面以進(jìn)行射流場測量與分析。

圖1 實(shí)驗(yàn)系統(tǒng)組成Fig.1 Sketch of the experiment system

圖2 環(huán)形噴嘴型式及射流斷面選取Fig.2 Annular nozzle and transverse sections

2 PDPA參數(shù)

PDPA的基本原理是散射光干涉原理,即測量物體散射光線時(shí)產(chǎn)生的多普勒頻移以獲得物體的移動(dòng)速度,與激光多普勒測速(LDA)技術(shù)最大的不同在于其接收光學(xué)單元有3個(gè)探測器且其位置不能與發(fā)射光束的對稱軸同線,在測量多維速度分量的同時(shí),亦可通過測量信號的相位得到示蹤粒子的粒徑[9]。本實(shí)驗(yàn)中采用丹麥Dantec公司生產(chǎn)的PDPA系統(tǒng),該系統(tǒng)由最大功率為5W的氬離子激光光源、二維綠光及藍(lán)光發(fā)射探頭、接收探頭、光纖驅(qū)動(dòng)器、光電接收器、信號處理器、控制軟件BSA Flow以及智能三維位移機(jī)構(gòu)等組成。實(shí)驗(yàn)中在前側(cè)向布置接收器,接收散射角為35°。該系統(tǒng)的主要性能參數(shù)如表1所示。

表1 PDPA系統(tǒng)的性能參數(shù)Table 1 Primary parameters of PDPA system

測量時(shí),激光光源產(chǎn)生的單色激光束經(jīng)光纖驅(qū)動(dòng)器后分成波長分別為514.5和488.0nm的綠光和藍(lán)光,經(jīng)過耦合傳輸后,進(jìn)入焦距為500mm的發(fā)射探頭。上下分射的兩束綠光可測量液滴軸向速度分量及液滴索太爾平均粒徑,水平分射的兩束藍(lán)光可對徑向速度分量進(jìn)行測量。4束光聚集在流場中的某一點(diǎn),形成一個(gè)橢球測量體(如圖3所示)。測量時(shí),發(fā)射探頭和接收探頭固定在智能三維坐標(biāo)架上,其同步移動(dòng)由計(jì)算機(jī)控制,最小移動(dòng)步長為0.1mm。每次采樣過程中的樣本液滴數(shù)設(shè)定為2000個(gè),通過PDPA采樣驗(yàn)證,液滴數(shù)符合測量精度要求。

圖3 PDPA原理圖Fig.3 Schematic diagram of PDPA

3 速度分布

3.1 軸向速度分布

軸向速度是表征射流能量的最重要指標(biāo),圖4表示出了在11MPa射流壓力條件下,6個(gè)射流斷面上的軸向速度分布。本文采用有量綱的表達(dá)方法,一是為了說明所研究射流的速度量級,二是在本研究中,射流的半擴(kuò)展厚→寬度難以精確確定,因?yàn)樵趯α魇吘夁M(jìn)行PDPA采樣時(shí),存在著一定的采樣誤差。

由圖4可以看出,由截面1～6,速度分布曲線在徑向的跨度增加,同時(shí)流束邊緣的速度分布趨于平坦。由距離噴嘴較近的截面1和2來看,環(huán)形射流在圓環(huán)中心發(fā)生聚并后,迅速發(fā)展成為圓柱射流。同時(shí)在截面1和2的速度分布圖上可發(fā)現(xiàn)射流中心的核心區(qū)較寬,且與核心區(qū)外層間存在著較大的速度梯度。自截面1～6,斷面最大軸向速度衰減百分比僅為0.62%,射流能量得到很好地保持。

圖4 軸向速度分布(p=11MPa)Fig.4 Axial velocity distribution under jet pressure of 11MPa

3.2 射流壓力對軸向速度分布的影響

圖5表示出截面1和6的軸向速度分布隨射流壓力的變化情況。如圖5(a)所示,隨著射流壓力的提高,截面1的軸向速度分布曲線基本上向縱軸整體平移,流束中心部分平移的量更大一些。在圖5(b)中,在不同壓力下,截面6上的軸向速度分布呈現(xiàn)光滑的二次曲線形式。遠(yuǎn)離噴嘴,射流束在徑向的擴(kuò)散增強(qiáng)。隨著射流壓力的提高,流束中間的射流流體加速更為明顯。在截面6的邊緣,不同壓力下的分布曲線交織在一起,說明混合層內(nèi)液體與空氣的相互作用更為強(qiáng)烈。

3.3 徑向平均速度分布

在較高壓力下的射流場中,徑向速度通常被忽略,但該文研究的圓環(huán)射流流束的直徑較大,徑向速度分布在分析流束擴(kuò)散程度和流束外緣液滴運(yùn)動(dòng)特性時(shí)起著重要的作用,故將6個(gè)截面上的徑向速度分布表示在圖6中?？梢钥闯?每個(gè)截面上的徑向流速分布都基本保持對稱,且由截面1～6,徑向平均速度分布曲線趨于平坦。同時(shí),由截面1～6,截面邊緣處的徑向速度逐漸減小,一方面,水射流與周圍空氣的摩擦使射流流體能量不斷衰減,另一方面,在流束邊緣處,大尺度液滴的破碎使液滴表面和空氣的接觸面積增加,液滴的徑向運(yùn)動(dòng)受到進(jìn)一步阻滯。

圖6 徑向速度分布(p=11MPa)Fig.6 Radial velocity distribution under jet pressure of 11MPa

4 湍流脈動(dòng)分析

4.1 軸向均方根速度

均方根速度是表征湍流脈動(dòng)程度的重要參數(shù),圖7表示出了在11MPa射流壓力條件下,截面1、3和5上的軸向均方根速度分布情況。在每個(gè)截面的射流中心區(qū),軸向均方根速度分布曲線均存在著明顯的波谷,而該波谷對應(yīng)著每個(gè)截面上軸向速度最大的區(qū)域。從另一角度看,在射流中心區(qū)外圍,分布著兩個(gè)軸向均方根速度的波峰,該處對應(yīng)著液體與周圍空氣開始發(fā)生強(qiáng)烈剪切作用的位置。

圖7 軸向均方根速度分布(p=11MPa)Fig.7 Axial root mean square velocity distribution under jet pressure of 11MPa

4.2 徑向均方根速度

對應(yīng)于圖6的徑向速度分布,圖8給出了截面1、3和5上的徑向均方根速度分布。從徑向均方根速度的量值來看,其與徑向速度的大小相當(dāng),說明液滴的徑向湍流脈動(dòng)很強(qiáng)烈。而從分布曲線來看,在遠(yuǎn)離噴嘴的截面5上,徑向均方根速度分布趨于平坦。在截面1和3處,截面邊緣處的徑向均方根速度分布曲線梯度較大,且射流流束中心的徑向均方根速度較小。

圖8 徑向均方根速度分布(p=11MPa)Fig.8 Radial root mean square velocity distribution under jet pressure of 11MPa

5 液滴粒徑分布

液滴是射流能量的載體,其尺寸及尺寸分布直接影響著射流的作用效果,在PDPA采樣過程中,對通過不同位置控制體的每個(gè)液滴都進(jìn)行了記錄。液滴尺寸的表達(dá)主要采用兩個(gè)參數(shù)：液滴數(shù)平均粒徑和索太爾平均粒徑。液滴數(shù)平均粒徑的定義為

索太爾平均粒徑可定義為

由于不同的粒徑表示方法得到的結(jié)果可能存在較大的差異,故圖9中同時(shí)表示出了11MPa射流壓力下的液滴數(shù)平均粒徑和索太爾平均粒徑分布?？梢钥闯?由截面1～4,液滴數(shù)平均粒徑在流束中心呈現(xiàn)略微下凹的分布形態(tài),在接近射流流束外緣時(shí),液滴數(shù)平均粒徑分布較為平坦。在截面5和6上,在流束外緣,隨著距射流軸線的徑向距離增加,液滴數(shù)平均粒徑又出現(xiàn)增大的趨勢。而索太爾平均粒徑分布則基本上與圖7中的軸向均方根速度分布相對應(yīng),在射流流束中心區(qū),湍動(dòng)程度小,液滴的索太爾平均粒徑大;而在流束外緣,湍動(dòng)程度大,且液滴與周圍空氣相互作用劇烈,液滴破碎幾率增大,故小尺寸液滴增加。

接近流束中心,液滴數(shù)平均粒徑和索太爾平均粒徑之間存在著明顯的差異,而在截面1的流束中心,這種差異最為顯著,說明該處的粒徑譜較寬,液滴尺寸分布不均勻。由于環(huán)形噴嘴的中心為一圓柱體,高速液流經(jīng)過時(shí)在圓柱體尾部形成滯止區(qū),在滯止區(qū)的外層存在著較大的速度梯度,高速和低速流體之間不斷發(fā)生剪切作用,而環(huán)形射流流束又在徑向迅速擴(kuò)散,噴嘴附近出現(xiàn)滯止區(qū)外層剝落和流動(dòng)干涉,這是導(dǎo)致此處液滴尺寸不均勻分布的重要因素之一。

6 液滴統(tǒng)計(jì)分析

為進(jìn)一步討論流場中的液滴行為及尺寸特征,在更短的時(shí)間尺度和更小的空間尺度條件下對液滴參數(shù)進(jìn)行統(tǒng)計(jì)分析。對在截面1和6的中心點(diǎn)和邊緣點(diǎn)處設(shè)置的控制體進(jìn)行監(jiān)測,對不同時(shí)刻通過控制體的單個(gè)液滴的粒徑和軸向速度進(jìn)行統(tǒng)計(jì)分析。

6.1 粒徑百分比

液滴粒徑與液滴數(shù)百分比之間的對應(yīng)關(guān)系如圖10所示。對于截面1的中心點(diǎn)和邊緣點(diǎn)處的控制體,其通過的不同粒徑液滴所占的百分比關(guān)系呈現(xiàn)明顯的統(tǒng)計(jì)特征,通過中心點(diǎn)控制體的液滴覆蓋的尺寸范圍更大,8μ m左右的液滴數(shù)量最多,約占整個(gè)樣本數(shù)量的10%,且壓力升高對不同粒徑液滴所占的百分比影響不大。通過截面1邊緣點(diǎn)處設(shè)置的控制體的液滴數(shù)量百分比與液滴粒徑的關(guān)系呈尖峰分布,占最大比例的液滴尺寸約為18μ m,占樣本數(shù)量的24%左右,且射流壓力升高到15MPa時(shí),該尺寸液滴所占比例增加至29%左右。

在圖10(b)中,截面6的液滴數(shù)量百分比與液滴粒徑的關(guān)系受射流壓力的影響同樣不明顯,且與截面1表現(xiàn)出相似的分布特征。對于中心點(diǎn)處的控制體,其通過的液滴的數(shù)量與尺寸對應(yīng)關(guān)系更為平坦。而對于邊緣點(diǎn),其關(guān)系曲線的峰值對應(yīng)的液滴粒徑達(dá)到36μ m左右,數(shù)量百分比卻下降至11%左右。射流壓力增加到15MPa時(shí),占最大百分比的液滴的粒徑略有減小,百分比的峰值卻增加至13.5%左右。

6.2 粒徑與速度統(tǒng)計(jì)關(guān)系

圖10 不同粒徑液滴的數(shù)量百分比Fig.10 Percentage of droplets with different diameters

圖11 不同軸向速度對應(yīng)的液滴數(shù)百分比Fig.11 Relation between droplet percentage and axial velocity

液滴數(shù)百分比與液滴軸向速度間的對應(yīng)關(guān)系如圖11所示。相對于截面6的邊緣點(diǎn),在截面1的邊緣點(diǎn)處,液滴數(shù)百分比曲線覆蓋的液滴軸向速度范圍較大,且曲線峰值對應(yīng)的液滴軸向速度量值較大。射流壓力對邊緣點(diǎn)處液滴數(shù)百分比曲線的形態(tài)沒有明顯的影響。而對于兩個(gè)截面的中心點(diǎn),射流壓力的增加使液滴數(shù)百分比曲線向橫坐標(biāo)的正向平移,且在中心點(diǎn)處,最大液滴數(shù)百分比的值均在40%左右,當(dāng)射流壓力增加到15MPa時(shí),該峰值略有下降。

7 結(jié) 論

(1)實(shí)驗(yàn)壓力條件下,環(huán)形射流的能量沿射流方向衰減緩慢,且射流流束中心的高速區(qū)較寬,該高速區(qū)內(nèi)的軸向和徑向的湍流脈動(dòng)程度均較弱;

(2)接近流束中心,液滴數(shù)平均粒徑和索太爾平均粒徑存在著較大的差異,尤其是在距噴嘴最近的截面中心位置,該差異更為顯著,說明射流流束中心的液滴粒徑譜較寬;

(3)通過截面邊緣控制體的小尺寸液滴占有絕對的數(shù)量優(yōu)勢,隨著距噴嘴的軸向距離增加,占最大數(shù)量比例的液滴的尺寸有所增加。通過射流中心位置控制體的高速液滴的數(shù)量占絕對優(yōu)勢,而邊緣控制體處的液滴數(shù)百分比與液滴軸向速度關(guān)系呈現(xiàn)明顯的二次曲線分布,具有明顯的統(tǒng)計(jì)特征。研究結(jié)果為環(huán)形射流在更深層面的應(yīng)用提供了參考。

[1] GEORGIOU G.Annular liquid jets at high Reynolds number[J].International Journal for Numerical Methods in Fluids.2003,42：117-130.

[2] 劉寧,張相炎.高速環(huán)形射流破碎霧化特性及動(dòng)態(tài)噴霧模擬[J].彈道學(xué)報(bào),2008,20(1)：11-14.

[3] BROECKHOVEN T,BROUNS M,VANHERZEELE J,et al.PIV measurements of a double annular jet for validation of numerical simulations[C].13th Int Symp.on Applications of Laser Techniques to Fluid Mechanics,Portugal,26-29 June,2006.

[4] CHRISTIAN DEL T,LARS B,JURG G,et al.Numerical and experimental investigation of an annular jet flow with large blockage[J].Journal of Fluids Engineering,2004,126：375-384.

[5] 薛勝雄.高壓水射流技術(shù)工程[M].安徽合肥：合肥工業(yè)大學(xué)出版社,2006.

[6] 宮偉力,趙海燕,安里千.高壓水射流結(jié)構(gòu)的紅外熱像特征[J].實(shí)驗(yàn)流體力學(xué),2008,22(3)：31-35.

[7] 張玉榮,金玉明,余志剛.高壓下液體射流破碎霧化特性實(shí)驗(yàn)研究[J].航空動(dòng)力學(xué)報(bào),1999,14(4)：409-412.

[8] 萬云霞,黃勇,朱英.液體圓柱射流破碎過程的實(shí)驗(yàn)[J].航空動(dòng)力學(xué)報(bào),2008,23(2)：208-214.

[9] 沈熊,魏乃龍,彭濤.應(yīng)用激光相位-多普勒系統(tǒng)測量霧化液滴顆粒和流動(dòng)特性[J].流體力學(xué)實(shí)驗(yàn)與測量,2000,14(2)：54-60.