康 燦，劉海霞，周 亮，楊敏官

(1.江蘇大學(xué)能源與動(dòng)力工程學(xué)院，江蘇鎮(zhèn)江212013;2.江蘇大學(xué)材料科學(xué)與工程學(xué)院，江蘇鎮(zhèn)江212013)

水射流是一種重要的流動(dòng)行為，也是多個(gè)重要工程領(lǐng)域內(nèi)關(guān)注的對(duì)象之一.在不同的射流壓力條件下，水射流流束的形態(tài)之間存在著很大的差別，不同形態(tài)射流束的內(nèi)部能量分布特征及破碎模式也不相同［1］.涉及多種非線性因素的射流束穩(wěn)定性問(wèn)題仍是目前具有挑戰(zhàn)性的難題之一，尤其是在較高的射流壓力條件下，環(huán)境流體對(duì)射流束的作用更加難以捕捉與?；?-3］.當(dāng)射流壓力超過(guò)200 MPa時(shí)，常規(guī)的流動(dòng)測(cè)量方法難以準(zhǔn)確獲取流動(dòng)參數(shù)分布.盡管一些先進(jìn)的流動(dòng)測(cè)量或顯示方法如紅外熱成像法等不斷在高壓射流的研究中得到實(shí)踐與發(fā)展，相關(guān)的定量研究結(jié)果仍然極為缺乏［4］.

近年來(lái)，對(duì)高壓水射流流動(dòng)開展的研究多集中在數(shù)值模擬方面，研究結(jié)論偏重定性分析.也有學(xué)者在900 m·s-1的高速純水射流沖蝕固體壁面的研究中，運(yùn)用有限元方法與光滑粒子流體動(dòng)力學(xué)方法相結(jié)合的手段，模擬分析了流場(chǎng)和射流的沖擊效應(yīng)，綜合考慮了流體和固體的相互影響［5-6］.在試驗(yàn)方法的探索中，有的學(xué)者為獲得射流沖擊固體壁面的沖擊力，采用在被沖擊的固體表面埋設(shè)傳感器的方法對(duì)射流破壞巖石的機(jī)理進(jìn)行了研究［7］.而通常的試驗(yàn)中多根據(jù)射流沖擊過(guò)的試樣的質(zhì)量損失來(lái)估計(jì)射流的作用力［8］.隨著材料表面檢測(cè)技術(shù)的發(fā)展，對(duì)試樣表面輪廓進(jìn)行光學(xué)掃描可以更為合理地評(píng)估水射流的打擊效果［9］.從工程應(yīng)用的角度來(lái)看，運(yùn)行參數(shù)如噴嘴移動(dòng)速度、靶距、射流壓力和流量等對(duì)水射流效果的影響倍受關(guān)注，針對(duì)水射流的特定應(yīng)用而開展的大多數(shù)試驗(yàn)研究也以尋找最優(yōu)化的運(yùn)行參數(shù)組合為目標(biāo)［5］.

水射流流束的內(nèi)部速度分布和流束波動(dòng)與被沖擊過(guò)的固體壁面的形貌特征之間存在著密切的關(guān)聯(lián).更為具體地講，受射流作用過(guò)的固體表面的形貌特征應(yīng)該與射流流束的局部特征相呼應(yīng).本研究遵循這一思路，以壓力為280 MPa的射流為研究對(duì)象，采用流場(chǎng)測(cè)量的方法，獲得射流束內(nèi)的能量集聚與脈動(dòng)特性;拍攝流束輪廓的瞬態(tài)波動(dòng)，并運(yùn)用圖像處理的方法對(duì)瞬時(shí)圖像進(jìn)行分析;進(jìn)而對(duì)磨料射流作用過(guò)的Al合金試樣表面進(jìn)行掃描電鏡觀察與三維輪廓分析，以期通過(guò)特征形貌來(lái)間接描述射流的本質(zhì)特性與沖蝕能力.

1 純水射流噴嘴結(jié)構(gòu)

試驗(yàn)中所采用的純水射流噴嘴結(jié)構(gòu)示意圖如圖1所示，圖中尺寸單位為mm，內(nèi)流道直徑為0.33 mm的人造藍(lán)寶石為主要磨損件，也是使水流加速的關(guān)鍵過(guò)流部件.噴嘴體外圍部分的材料為不銹鋼.射流自噴嘴噴出前會(huì)突然擴(kuò)散至直徑為1.00 mm的直段.試驗(yàn)過(guò)程中側(cè)重考慮射流流束的形態(tài)與能量特征，并不考慮噴嘴內(nèi)流道內(nèi)的能量損失與可能發(fā)生的空化現(xiàn)象.

圖1 純水射流噴嘴剖面圖

2 流動(dòng)測(cè)量試驗(yàn)

2.1 試驗(yàn)系統(tǒng)與測(cè)量體布置

目前在流動(dòng)速度的光學(xué)測(cè)量方面常采用的2種方法為粒子圖像速度場(chǎng)(PIV)法和激光多普勒粒子分析儀(PDPA)法.一般的PIV系統(tǒng)不能滿足湍流脈動(dòng)的測(cè)量，其在高壓射流測(cè)量中的應(yīng)用也鮮見報(bào)道.本試驗(yàn)中采用PDPA對(duì)水射流流場(chǎng)進(jìn)行流速測(cè)量，PDPA技術(shù)為單點(diǎn)測(cè)量技術(shù)，可以同時(shí)獲得示蹤粒子的流動(dòng)速度和粒徑，是定量測(cè)量射流流場(chǎng)的重要手段之一［10］.試驗(yàn)在自由射流條件下進(jìn)行，采用如圖2所示的前向接收模式.

圖2 試驗(yàn)系統(tǒng)示意圖

試驗(yàn)中未添加專門的示蹤粒子，借助介質(zhì)內(nèi)部的雜質(zhì)和散布在射流束中的液滴進(jìn)行速度追蹤.由于試驗(yàn)中不涉及液滴粒徑值的統(tǒng)計(jì)問(wèn)題，且測(cè)量的開展僅局限于近噴嘴區(qū)，所以認(rèn)為該示蹤方法可行.建立以噴嘴出口斷面中心為坐標(biāo)原點(diǎn)、z軸正向與射流主流方向一致的三維直角坐標(biāo)系.流場(chǎng)測(cè)量在z=3.0 mm和z=7.0 mm射流束截面上進(jìn)行.在每一個(gè)截面上，借助控制軟件實(shí)現(xiàn)PDPA發(fā)射和接收探頭的同步移動(dòng)，從而遍歷如圖3所示的所有測(cè)量體.

圖3 測(cè)量截面的測(cè)量體布置

2.2 流動(dòng)試驗(yàn)結(jié)果與分析

運(yùn)用PDPA技術(shù)測(cè)量得到的z=3.0，7.0 mm截面上的速度分布分別如圖4a，b所示，圖中平均速度單位為m·s-1，從速度分布形態(tài)來(lái)看，z=3.0 mm截面的速度分布呈現(xiàn)“平頂帽”形態(tài)，該區(qū)域的速度最大值達(dá)到了580 m·s-1.在高速區(qū)域的邊緣存在著明顯的速度梯度，該速度梯度預(yù)示著連續(xù)射流與強(qiáng)烈霧化區(qū)之間的分界位置.射流外層與外界空氣間的速度差是擾動(dòng)局部射流束甚至導(dǎo)致射流霧化的關(guān)鍵因素.在z=7.0 mm截面上，“平頂帽”狀的高速區(qū)依然存在，該高速區(qū)內(nèi)的總體速度量值較z=3.0 mm斷面小，且該高速區(qū)的面積較大.同時(shí)，z=7.0 mm截面上的高速區(qū)與低速區(qū)之間的速度梯度相對(duì)較小.

圖4 射流束截面上的平均速度分布

流體的流動(dòng)速度常被分解為2部分:①平均速度;②脈動(dòng)速度.脈動(dòng)速度是影響射流流體內(nèi)部動(dòng)量傳遞和射流束破碎的重要因素.

脈動(dòng)速度的均方根值，也稱為均方根速度，定義如下:

式中:N為同一測(cè)量體內(nèi)的采樣點(diǎn)數(shù)目;ui為單個(gè)粒子的速度;ua為基于同一測(cè)量體內(nèi)的所有粒子的瞬時(shí)速度而獲得的平均速度.

水射流攜帶的動(dòng)能是射流沖蝕固體壁面的重要因素.從2個(gè)截面的速度分布形態(tài)已可以推斷在不同的靶距處射流沖擊能力的明顯差異.脈動(dòng)速度是反映湍流物理本質(zhì)的重要參數(shù)，與射流中的湍動(dòng)能、湍動(dòng)能耗散率及耗散能量的最小渦尺度密切相關(guān).在高雷諾數(shù)Re流動(dòng)條件下，流場(chǎng)中的湍流脈動(dòng)是不可忽視的因素，因此，根據(jù)PDPA的測(cè)量結(jié)果構(gòu)建了脈動(dòng)速度在2個(gè)截面的分布，如圖5所示，圖中均方根速度單位為m·s-1.

圖5 射流束截面上的均方根速度分布

從圖5可以看出:2個(gè)截面的均方根速度分布與平均速度分布相呼應(yīng)，在z=3.0 mm截面上，流束截面中心處出現(xiàn)明顯的均方根速度低谷，該低谷的直徑約為0.8 mm.可以認(rèn)為:在該低谷所代表的區(qū)域內(nèi)，射流流束處于穩(wěn)定狀態(tài)，外界空氣對(duì)該區(qū)域的干擾可以忽略.同時(shí)，在z=7.0 mm截面，均方根速度分布的形態(tài)與z=3.0 mm截面相似，而均方根速度低谷覆蓋的面積增大.值得注意的是，在該低谷邊緣出現(xiàn)了明顯的周向不對(duì)稱性，射流有向不穩(wěn)定狀態(tài)過(guò)渡的趨勢(shì).從均方根速度的量級(jí)來(lái)看，在此高Re射流中，射流核心區(qū)的湍流強(qiáng)度相對(duì)于中低Re射流較小.

3 流束形態(tài)波動(dòng)試驗(yàn)

PDPA系統(tǒng)采用對(duì)測(cè)量體內(nèi)所通過(guò)的粒子信息進(jìn)行統(tǒng)計(jì)的方法，無(wú)法同時(shí)獲取流束形態(tài)的波動(dòng).高速數(shù)碼攝像技術(shù)是近年來(lái)在流動(dòng)顯示方面得到成功應(yīng)用的一種手段，其對(duì)于瞬態(tài)流動(dòng)現(xiàn)象的捕捉具有很強(qiáng)的針對(duì)性.為了記錄流束輪廓的波動(dòng)過(guò)程，應(yīng)用高速數(shù)碼攝像機(jī)，輔以LED光源，對(duì)近噴嘴區(qū)流束段進(jìn)行了拍攝.拍攝過(guò)程中采用的拍攝頻率為5 000幀·s-1，即相鄰的2幅圖像之間的時(shí)間間隔為0.2 ms.

圖6為對(duì)應(yīng)于某一拍攝時(shí)刻的流束形態(tài)原始圖像，該圖像中明暗區(qū)別的流束區(qū)域直觀地表示出了流束的集聚與徑向擴(kuò)散特征.為定量分析流束的波動(dòng)幅度沿射流方向的發(fā)展及瞬態(tài)變化，采用圖像處理軟件對(duì)原始圖像進(jìn)行處理，獲得了射流流束中心線上的圖像灰度值隨時(shí)間的變化.

圖6 射流束形態(tài)的原始圖像

圖7為隨機(jī)選取的6幅連續(xù)圖像(t1-t6)中的流束中心線灰度值分布，距離噴嘴出口1.0～6.0 mm的流束段，灰度值較高，且波動(dòng)幅度較小，該段的集束性較強(qiáng).而在距離噴嘴出口6.0～11.0 mm處，灰度值下降明顯，說(shuō)明流束的徑向擴(kuò)散嚴(yán)重.工程上應(yīng)用的破壞性射流的打擊靶距一般為2.0～5.0 mm，這一點(diǎn)從流束波動(dòng)特征上也可以反映出來(lái).從靶距1.0～11.0 mm的范圍內(nèi)的波動(dòng)差異也可以得出結(jié)論，射流束的連續(xù)性沿著射流前進(jìn)方向不斷被破壞，射流束對(duì)固體壁面的作用能力也會(huì)隨之減弱.

圖7 用圖像灰度值表示的射流束輪廓波動(dòng)特征

4 固壁表面形貌試驗(yàn)

為了借助射流作用過(guò)的固體表面形貌來(lái)描述射流流束的能量分布及其對(duì)固體壁面的沖蝕作用，選用厚度為8.0 mm的Al合金板作為射流沖擊對(duì)象，以其表面形貌特征作為描述射流能量分布的依據(jù).試驗(yàn)中采用與文獻(xiàn)［11］中相同結(jié)構(gòu)的噴嘴，選用粒徑為177 μm的石英砂作為磨料，混合管的出口直徑為0.8 mm.由于水流的流速較高，可近似認(rèn)為在所關(guān)注的z=3.0，7.0 mm射流束截面上，磨料的z向速度分量與位于同一位置的流體質(zhì)點(diǎn)的z向速度分量大小相同.試驗(yàn)中采取的射流壓力同樣為280 MPa，水的流量為 3.7 L·min-1，噴嘴支撐機(jī)構(gòu)與機(jī)床的懸臂相連，通過(guò)機(jī)床控制系統(tǒng)可實(shí)現(xiàn)噴嘴在3個(gè)直角坐標(biāo)方向的移動(dòng).本試驗(yàn)中，射流的主流方向始終與Al合金板垂直，噴嘴僅沿同一水平方向平移，平移速度為200 mm·min-1.噴嘴的平移速度決定了射流在試樣表面同一位置停留的時(shí)間，對(duì)射流作用效果有很大影響.

在靶距分別為3.0，7.0 mm 條件下，射流均能將Al合金板切穿，從而得到了被射流切割而形成的試樣斷面.通常認(rèn)為該切割過(guò)程的主導(dǎo)因素在于水楔作用和石英砂粒磨削作用.盡管利用有限元方法可以對(duì)受到砂粒沖擊的試樣進(jìn)行瞬態(tài)應(yīng)力模擬，卻無(wú)法獲得試樣斷面的形貌［11］.為觀察試樣斷面的沖蝕形態(tài)，首先采用掃描電鏡SEM對(duì)在不同靶距條件下獲得的試樣進(jìn)行了表面觀測(cè).由于水射流的作用機(jī)理是冷態(tài)作用，不存在組織形貌的改變，采用SEM得到的結(jié)果也僅僅反映試樣表面的一般特征.試樣表面的SEM照片如圖8所示.靶距為3.0 mm時(shí)的SEM觀測(cè)結(jié)果如圖8a所示，射流方向?yàn)榇怪毕蛳拢纬闪嗣黠@的沖蝕條紋，且條紋分布較為貫通、均勻.在該圖像中突出顯示了磨削痕和截?cái)嗪圻@2個(gè)明顯的射流作用痕跡.在相同的放大倍數(shù)下，靶距為7.0 mm時(shí)的沖蝕結(jié)果如圖8b所示，射流方向同樣為垂直向下，射流沖蝕所致的形貌較為雜亂，幾乎不存在上下貫通的條紋.同時(shí)，磨削痕相對(duì)傾斜，截?cái)嗪鄣漠a(chǎn)生更顯得吃力.在圖8b中還存在1個(gè)明顯的撞擊痕特征，該特征的出現(xiàn)一方面是由于靶距增加、流束徑向擴(kuò)散而引起砂粒的徑向速度值增加，另一方面，砂粒動(dòng)能隨靶距增加的衰減也是不可忽視的因素.可以預(yù)測(cè)，在砂粒與試樣發(fā)生碰撞的過(guò)程中，砂粒的破碎和砂粒嵌入試樣表面2種現(xiàn)象同樣存在著較大的發(fā)生概率，該作用過(guò)程到目前為止尚未獲得全面且準(zhǔn)確的解釋.

圖8 試樣表面的SEM照片

為了定量評(píng)估試樣表面的形貌特征，采用光學(xué)輪廓掃描技術(shù)對(duì)2種靶距條件下獲得的切割表面進(jìn)行了粗糙程度的三維測(cè)量.在圖9a，b中分別突出了磨削痕和截?cái)嗪圻@2個(gè)特征.在圖9中，縱坐標(biāo)z的原點(diǎn)并不在噴嘴出口斷面的中心，該坐標(biāo)僅表示所選取的采樣面沿射流主流方向的尺寸大小.同時(shí)，用色譜分布的方式來(lái)反映試樣表面上的采樣點(diǎn)高度分布.可以看出，受到磨料射流作用的試樣表面存在著形態(tài)各異的痕跡，這些痕跡的尺度與砂粒的單次或多次沖擊有關(guān)，更與射流能量的衰減密切關(guān)聯(lián).相對(duì)于圖9a，圖9b所示的形貌特征更加不規(guī)則.三維輪廓顯示方法是一種定量的方法，基于這種方法，可以獲得某一局部區(qū)域上的統(tǒng)計(jì)粗糙度.在2種靶距條件下，各隨機(jī)選取了3個(gè)觀測(cè)面，且所有觀測(cè)面的尺寸與圖9所示相同.在每一個(gè)觀測(cè)面上均取353 280個(gè)測(cè)點(diǎn)，記錄每個(gè)測(cè)點(diǎn)的高度值，獲得了與測(cè)點(diǎn)高度值對(duì)應(yīng)的概率密度，如圖10所示.

圖9 試樣表面的局部三維形貌

圖10 試樣表面采樣點(diǎn)高度的概率密度分布

測(cè)量得到的相對(duì)高度值均在-140～40 μm.在相對(duì)高度值小于-20 μm和大于16 μm時(shí)，其對(duì)應(yīng)的概率密度均小于0.5%，故在圖10中略去了這2個(gè)區(qū)間的概率密度顯示.由圖10可以看出:在z=3.0 mm時(shí)，3條概率密度曲線均向相對(duì)高度0值(參考水平面)收縮，且均呈現(xiàn)明顯的單峰分布.相比之下，z=7.0 mm時(shí)，對(duì)區(qū)域整體粗糙度有著顯著影響的采樣點(diǎn)高度值范圍較寬，且所隨機(jī)抽取的3個(gè)區(qū)域中均出現(xiàn)了雙峰分布特征，這也與圖8所示的不同靶距條件下的形貌特征差異相對(duì)應(yīng).沿著圖10所示的思路可以將各區(qū)域的采樣點(diǎn)高度在整個(gè)試樣表面進(jìn)行平均化，從而獲得評(píng)價(jià)射流對(duì)試樣表面作用的整體效果的關(guān)鍵度量參數(shù).

5 結(jié)論

1)高速水射流流束截面中心存在著高速核心區(qū)，其提供了破壞固體壁面的主要能量.該高速核心區(qū)同時(shí)以湍流均方根速度的谷狀分布為重要特征.在高速核心區(qū)外圍存在著明顯的速度梯度.隨著射流向下游延伸，高速核心區(qū)不斷沿徑向擴(kuò)散，且射流所攜帶的動(dòng)能不斷衰減.

2)由射流束中心線的圖像灰度值沿射流方向的分布可以判斷射流束輪廓的分段波動(dòng)特征，由此結(jié)果也可以映射出在高Re條件下，環(huán)境空氣對(duì)射流束的擾動(dòng)是一種小時(shí)間尺度現(xiàn)象，沿射流方向，射流束抵御這種擾動(dòng)的能力存在著明顯差異.

3)在受磨料水射流作用的Al合金試樣表面上，磨削痕、截?cái)嗪?、撞擊痕等明顯的形貌特征與復(fù)雜的射流作用密切相關(guān).盡管目前能夠從這些形貌特征中提取表面粗糙度、采樣點(diǎn)高度的概率密度等重要信息，但這些形貌特征的產(chǎn)生機(jī)制仍有待于進(jìn)一步探索.

References)

［1］Delteil J，Vincent S，Erriguible A，et al.Numerical investigations in Rayleigh breakup of round liquid jets with VOF methods［J］.Computers＆Fluids，2011，50:10-23.

［2］Shinjo J，Umemura A.Detailed simulation of primary atomization mechanisms in diesel jet sprays(isolated identification of liquid jet tip effects)［J］.Proceedings of the Combustion Institute，2011，33:2089-2097.

［3］萬(wàn)云霞，黃 勇，朱 英.液體圓柱射流破碎過(guò)程的實(shí)驗(yàn)［J］.航空動(dòng)力學(xué)報(bào)，2008，23(2):208-214.Wan Yunxia，Huang Yong，Zhu Ying.Experiment on the breakup process of free round liquid jet［J］.Journal of Aerospace Power，2008，23(2):208-214.(in Chinese)

［4］Cui Longlian，An Liqian，Mao Lingtao，et al.Application of infrared thermal testing and mathematical models for studying the temperature distributions of the highspeed waterjet［J］.Journal of Materials Processing Technology，2009，209:4360-4365.

［5］Hascalik A，?ayda? U，Gürün H.Effect of traverse speed on abrasive waterjet machining of Ti-6Al-4V alloy［J］.Materials and Design，2007，28:1953-1957.

［6］Ma Li，Bao Ronghao，Guo Yimu.Waterjet penetration simulation by hybrid code of SPH and FEA［J］.International Journal of Impact Engineering，2008，35:1035-1042.

［7］Guha A，Barron R M，Balachandar R.An experimental and numerical study of water jet cleaning process［J］.Journal of Materials Processing Technology，2011，211:610-618.

［8］Kong M C，Axinte D，Voice W.Aspects of material removal mechanism in plain waterjet milling on gamma titanium aluminide［J］.Journal of Materials Processing Technology，2010，210:573-584.

［9］Osta A R，Lee J，Sallam K A，et al.Study of the effects of the injector length/diameter ratio on the surface properties of turbulent liquid jets in still air using X-ray imaging［J］.International Journal of Multiphase Flow，2012，38:87-98.

［10］康 燦，張 峰，楊敏官，等.環(huán)形水射流流場(chǎng)的實(shí)驗(yàn)研究與統(tǒng)計(jì)分析［J］.實(shí)驗(yàn)流體力學(xué)，2011，25(1):7-12.Kang Can，Zhang Feng，Yang Minguan，et al.Experimental and statistical analysis of water jet discharged from an annular nozzle［J］.Journal of Experiments in Fluid Mechanics，2011，25(1):7-12.(in Chinese)

［11］康 燦，周 亮，王育立，等.磨料水射流沖擊材料壁面的模擬與試驗(yàn)［J］.江蘇大學(xué)學(xué)報(bào):自然科學(xué)版，2013，34(3):276-280.Kang Can，Zhou Liang，Wang Yuli，et al.Simulation and experiment of impinged material surface with abrasive water jet［J］.Journal of Jiangsu University:Natural Science Edition，2013，34(3):276-280.(in Chinese)