截斷式脈沖射流流場結(jié)構(gòu)模擬與沖蝕硬巖能力分析

陸朝暉，盧義玉，Michael Hood，潘林華，賀 培

(1.重慶地質(zhì)礦產(chǎn)研究院 國土資源部頁巖氣資源勘查重點實驗室，重慶 400042；2.重慶大學(xué) 煤礦災(zāi)害動力學(xué)與控制國家重點實驗室，重慶 400044；3.昆士蘭大學(xué) 采礦合作研究中心，澳大利亞 布里斯班 4072)

截斷式脈沖射流流場結(jié)構(gòu)模擬與沖蝕硬巖能力分析

陸朝暉1,2,3，盧義玉2，Michael Hood3，潘林華1，賀 培1

(1.重慶地質(zhì)礦產(chǎn)研究院 國土資源部頁巖氣資源勘查重點實驗室，重慶 400042；2.重慶大學(xué) 煤礦災(zāi)害動力學(xué)與控制國家重點實驗室，重慶 400044；3.昆士蘭大學(xué) 采礦合作研究中心，澳大利亞 布里斯班 4072)

孔盤截斷式脈沖射流能夠充分利用水錘壓力、高速側(cè)向流、循環(huán)載荷等效應(yīng)，破巖效率高，在硬巖破碎中具有較高的開發(fā)潛力?；赩OF模型和動網(wǎng)格技術(shù)，建立與截斷式脈沖射流生成裝置幾何及運動過程一致的兩相流瞬態(tài)計算模型，結(jié)合高速攝像技術(shù)，研究了截斷式脈沖水射流的流體結(jié)構(gòu)的動態(tài)演變動力特征，測試了截斷式脈沖水射流的破巖能力，并分別從流暢結(jié)構(gòu)和破巖能力上與常規(guī)圓柱水射流進行了對比分析。研究結(jié)果表明：截斷式脈沖射流有助于在射流末端形成偏轉(zhuǎn)的瓦狀尖體結(jié)構(gòu)，可有效減弱或避免激波阻力的影響；截斷式脈沖水射流在破巖過程中，巖石靶物沖蝕孔口呈現(xiàn)非軸對稱形態(tài)，與常規(guī)圓柱射流相比，截斷式脈沖射流能夠大幅度提升破巖速度，增大破巖體積。

孔盤截斷式脈沖射流；流場結(jié)構(gòu)；動網(wǎng)格技術(shù)；硬巖沖蝕

脈沖水射流能夠有效利用其產(chǎn)生的水錘壓力、高速側(cè)向應(yīng)力和高頻動載荷等特性，在硬巖中形成剪切、拉伸和疲勞破壞，比連續(xù)射流更易形成階躍式破碎，被廣泛應(yīng)用到石油鉆井、采礦、隧道掘進等工業(yè)領(lǐng)域[1]。采用不同裝置不同動力源可產(chǎn)生不同類型的脈沖射流[2]，但在眾多脈沖射流類型中，孔盤截斷法是其中一種較為簡單易控的方法?？妆P截斷法能生成整潔脈沖，即相鄰單個脈沖完全斷開，有利于消除水墊影響，充分利用脈沖射流的水錘壓力，在破碎硬巖方面具有較高的開發(fā)潛力[3-4]。

高速脈沖射流對沖擊靶物的作用效果取決于水錘壓力的持續(xù)作用時間，而持續(xù)時間的長短又主要取決于脈沖射流流場，尤其是射流頭部流場結(jié)構(gòu)[5-6]。因此，與靶物碰撞初期的脈沖射流頭部流場結(jié)構(gòu)決定了射流的沖蝕能力，掌握射流流場結(jié)構(gòu)信息是研究其破硬巖機理的必要前提。高壓脈沖射流流場具有高速湍動的瞬變特性，截斷式脈沖射流流場結(jié)構(gòu)更為特殊，依靠純理論分析計算非常復(fù)雜。高速攝像技術(shù)能對射流的外部真實結(jié)構(gòu)及生成過程進行可視化研究[7]，但難以全面體現(xiàn)射流流場的變化發(fā)展。近年來，CFD(Computational Fluid Dynamics)技術(shù)的發(fā)展使流場細(xì)節(jié)的可視化成為可能，許多學(xué)者采用此方法對水射流進行了數(shù)值模擬研究[8-9]。

為了獲得截斷式脈沖射流的流場細(xì)節(jié)信息，筆者采用Fluent軟件，基于VOF(Volume of Fluid)模型和3D動網(wǎng)格技術(shù)，建立與截斷式脈沖射流生成裝置幾何及運動過程一致的計算模型，采用高速攝像法獲取截斷式脈沖射流的外部流場結(jié)構(gòu)加以驗證。應(yīng)用驗證后模型，模擬單個射流脈沖在100 mm經(jīng)驗最優(yōu)靶距范圍內(nèi)頭部結(jié)構(gòu)的動態(tài)演變及沖蝕靶體的動力特征，用于對其水錘壓力特性研究提供流場信息。為突出截斷式脈沖射流頭部流場結(jié)構(gòu)的特點，對比分析了與普通圓柱射流流場結(jié)構(gòu)的區(qū)別，研究了尖體流場演變機理。分別采用射流速度、尺度相同的兩種類型射流對花崗巖立方體進行了沖蝕實驗，對比分析了截斷式脈沖水射流的破巖機制。截斷式脈沖射流流場結(jié)構(gòu)與破巖機制的研究將對該類型射流在硬巖破碎應(yīng)用方面的優(yōu)化設(shè)計提供有效的的理論參考。

1 實驗系統(tǒng)設(shè)置

1.1截斷式脈沖射流生成原理

孔盤截斷式脈沖射流生成裝置如圖1所示，主要由高壓泵、噴嘴、液壓馬達(dá)、孔盤等組成。在高壓泵的作用下，高壓水流經(jīng)噴嘴高速噴出形成連續(xù)射流。液壓馬達(dá)帶動外置帶孔圓盤高速旋轉(zhuǎn)，可使連續(xù)射流周期性通過或遮斷，從而形成不連續(xù)的間斷射流，即脈沖射流。該裝置通過改變孔盤上切槽的數(shù)量、寬度或孔盤轉(zhuǎn)速即可改變單脈沖長度及脈沖頻率，對射流參數(shù)具有較高的可控性。

圖1 孔盤截斷式脈沖射流生成裝置示意圖Fig.1 Device for producing interrupted water jet by slotted disc

1.2高速攝像實驗設(shè)置

設(shè)計孔盤直徑為220 mm，厚度為10 mm，均勻布設(shè)16個盤孔，噴嘴軸線距孔盤軸線100 mm，盤孔寬度為10.2 mm。噴嘴直徑為3.429 mm，垂直于孔盤布置，間距9 mm。系統(tǒng)壓力設(shè)定為80 MPa，孔盤由液壓馬達(dá)帶動旋轉(zhuǎn)，轉(zhuǎn)速控制在780 r/min。拍攝時選用德國HiSpec 2G高速攝像機，將拍攝速度設(shè)置為7 437 fps，即連續(xù)兩張照片間的時間間隔為134.5 μs；曝光時間設(shè)置為2 μs；分辨率設(shè)置為592×122。假設(shè)孔盤旋轉(zhuǎn)時產(chǎn)生的每個單脈沖射流具有相似的非軸對稱結(jié)構(gòu)，為拍攝到截斷式脈沖射流的特有結(jié)構(gòu)，將高速攝像機鏡頭垂直于射流軸線方向，鏡頭視角中心對準(zhǔn)射流軸線和拍攝區(qū)域中點。圖2(a)為高速攝像實驗系統(tǒng)的實景圖。

1.3沖蝕花崗巖實驗設(shè)置

截斷式脈沖射流生成系統(tǒng)的設(shè)置與高速攝像實驗系統(tǒng)設(shè)置一致，獲得脈沖射流的頻率為208 Hz，單脈沖長度為634.05 mm，脈沖間隔距離為1 192.87 mm，功率為 72.9 kW。為測試沖蝕物理量，采用氣動的擋板開關(guān)執(zhí)行射流的通斷，共中斷實驗4次，時間間隔為 5 s，5 s，5 s，15 s，30 s和4 min，累計總沖蝕時間為5 min。其中沖蝕體積采用注射劑向沖蝕孔注射液體的方法獲得。實驗選用邊長為250 mm的花崗巖立方體，其物理力學(xué)屬性如表1，沖蝕靶距為100 mm(噴嘴出口至巖樣接觸面垂直距離)。圖2(b)為沖蝕巖石實驗實景圖。

(a) (b)圖2 高速攝像與沖蝕實驗實景圖Fig.2 The settings of the high speed camera and granite cube breaking experiment under interrupted water jet

表1 花崗巖物理力學(xué)參數(shù)Tab.1 Mechanical parameters of the granite sample

2 數(shù)值模擬

高速連續(xù)射流被孔盤干擾后，形成一系列射流脈沖。由于孔槽在孔盤均勻分布，且孔盤以勻速旋轉(zhuǎn)，因此理論上在一定泵壓的條件下，每個射流脈沖應(yīng)具有相同的射流結(jié)構(gòu)。故本數(shù)值模擬重點研究單個截斷式射流脈沖，模擬計算從孔盤阻擋射流至射流逐步從盤孔射出的整個射流頭部結(jié)構(gòu)的形成過程及射流沖蝕靶面的特殊流動狀態(tài)。在一定水動力條件下，噴嘴結(jié)構(gòu)是射流流場的主要影響因素[10]，但與通過盤孔的過程相比對截斷式射流結(jié)構(gòu)影響較小，本模擬研究主要針對后者。

2.1控制方程

脈沖射流噴入大氣介質(zhì)的流動為液、氣兩相流，對于這種非混合流體流動問題，確定相界面的位置和形狀后即可參照單相流體的處理方法來處理[11-12]。VOF模型通過求解一套動量方程和跟蹤處理穿過計算域的每一種流體的體積分?jǐn)?shù)來模擬兩相不能混合的流體，在追蹤液氣自由界面模擬中獲得了廣泛應(yīng)用，故選用VOF模型來追蹤流場中的兩相界面。

相間界面的追蹤是通過求解液氣兩相的體積分?jǐn)?shù)連續(xù)方程來實現(xiàn)的，第q相的體積分?jǐn)?shù)連續(xù)方程

(1)

其中，液、氣兩相還應(yīng)滿足約束條件

(2)

式中，αq為第q相的體積分?jǐn)?shù)，q=1為液相，q=2為氣相。

連續(xù)方程

(3)

動量方程

ρg+F

(4)

式中：ρ為密度；v為速度矢量；p為靜壓；μ為流體動力黏性系數(shù)；ρg為重力體積力；F為外部體積力。

因涉及射流的發(fā)散、解裂與霧化，沿兩相界面的表面張力效應(yīng)不可忽略[13]。采用Brackbill等提出的連續(xù)表面力模型(Continuum Surface Force，CSF)，在該模型下，考慮表面張力為常數(shù)的情況，且僅考慮兩相界面的法向力。根據(jù)散度定理，作用于相界面的表面張力可以表述為體積力，作為源項并入動量方程中

(5)

由表面張力導(dǎo)致的源項只在界面的某一側(cè)出現(xiàn)，文中指液相所在的一側(cè)。

(6)

Realizablek-ε湍流模型

Gb-ρε-YM+Sk

(7)

(8)

(9)

輸運方程中出現(xiàn)的物性參數(shù)取各相體積分?jǐn)?shù)的加權(quán)平均

(10)

(11)

2.2模型建立及網(wǎng)格劃分

在所研究的問題中，孔盤作為運動的實體，隨著孔盤的旋轉(zhuǎn)，射流的流場形態(tài)發(fā)生變化，這種流場形狀由于邊界運動而隨時間改變的問題可以采用動網(wǎng)格模型進行模擬。邊界的運動方式采用UDF(User Defined Functions)定義在網(wǎng)格面或網(wǎng)格區(qū)域上。由于流場中包含運動與不運動兩種區(qū)域，故在模型設(shè)置中用Fluent軟件提供的滑動界面功能將各區(qū)域連接起來，將它們組合在初始網(wǎng)格中以進行識別。

在具有移動邊界的任意控制體積V內(nèi)，流場用于描述流體通用標(biāo)量的積分形式的守恒方程可表示為

(12)

式中：ρ為流體密度；u為流動速度矢量；ug為網(wǎng)格移動速度；Γ為擴散系數(shù)；Sφ為φ的源項。

所研究問題的計算區(qū)域劃分為空氣間隔段、盤孔段和空氣段，各段軸向長度分別為9 mm、10 mm和81 mm，從射流入口至靶板處的軸向長度為100 mm。噴嘴直徑、盤孔尺寸及相對位置均參照射流結(jié)構(gòu)實驗的設(shè)置，其中噴嘴直徑為3.429 mm，盤孔寬度為10.2 mm，繞設(shè)定的孔盤軸線旋轉(zhuǎn)線速度為8.16 m/s。三個計算區(qū)域的運動狀態(tài)分別是：空氣間隔段與空氣段為靜止?fàn)顟B(tài)，盤孔段相對空氣間隔段和盤孔段繞軸線滑移旋轉(zhuǎn)，計算一個單脈沖的生成過程。在ICEM軟件應(yīng)用分區(qū)和分塊的劃分網(wǎng)格技術(shù)，對計算區(qū)域進行六面體網(wǎng)格劃分，網(wǎng)格數(shù)量總計200萬個。計算區(qū)域及網(wǎng)格劃分如圖3所示。

圖3 3D動態(tài)計算區(qū)域及六面體網(wǎng)格劃分Fig.3 3D dynamic computing domain and hexahedral grid

2.3離散方法

采用有限體積法將基本控制方程進行離散，生成離散方程組。為使VOF多相流模型求解精確，在對空間域進行離散化處理時采用修正的HRIC格式。在時間域離散上，采用隱式時間積分方案。為獲得VOF模型的精確瞬態(tài)特性，體積分?jǐn)?shù)方程采用顯式時間離散格式，位于兩相界面附近的單元通量采用幾何重建格式進行插值。對離散后的控制方程組求解時采用PISO(Pressure Implicit Split Operator)算法。

2.4邊界條件及收斂判據(jù)

射流入口采用壓力進口邊界條件，以湍動強度和水力直徑來設(shè)定入口的湍動狀態(tài)，空氣間隔段及空氣段區(qū)域的側(cè)面采用壓力出口，距射流入口100 mm處的靶面設(shè)定為無滑移光滑壁面。為了便于與實驗采集的數(shù)據(jù)相對比，計算過程中采用固定時間步長，且因所模擬的是微秒量級內(nèi)的高速湍動的瞬態(tài)問題，故選用小時間步長，設(shè)定為1×10-6。同時為了在運算的每一時間步獲得收斂解，將每個時間步長的最大迭代次數(shù)設(shè)定為50次，采用計算殘差曲線來評測整個計算區(qū)域內(nèi)所有守恒方程是否滿足要求來判斷解的收斂與否。

3 數(shù)值模擬結(jié)果及破巖效果分析

3.1數(shù)值模擬與高速攝像實驗圖片對比

圖4(b)為截斷式脈沖射流頭部結(jié)構(gòu)3D數(shù)值模擬結(jié)果。射流頭部形成一個前寬后窄、前薄后厚的瓦片狀尖體結(jié)構(gòu)，最前端射流向與盤孔旋轉(zhuǎn)方向相反的發(fā)生偏轉(zhuǎn)。模擬獲得的射流頭部偏轉(zhuǎn)特性與相同實驗設(shè)置條件下拍攝的高速攝像圖片(圖4(a))在形態(tài)與偏轉(zhuǎn)趨勢上具有較高的一致性。由高壓泵供水的壓力脈動、噴嘴與孔盤結(jié)構(gòu)、空氣淹沒環(huán)境的速度場狀態(tài)等未量化因素的影響，實際生成的射流結(jié)構(gòu)存在一定霧化擴散，但數(shù)值模擬結(jié)果應(yīng)對深入了解實際射流宏觀流場結(jié)構(gòu)與變化規(guī)律上具有一定參考意義。數(shù)值模擬所展示的在理論最優(yōu)靶距范圍內(nèi)的瓦狀尖體結(jié)構(gòu)，一定程度上再現(xiàn)了實際射流包裹在霧化液滴內(nèi)部的主體射流結(jié)構(gòu)，該部分是沖擊破碎巖石的最有效區(qū)域，深入探討其流場結(jié)構(gòu)與演化特點，對研究截斷式脈沖射流破巖機理和優(yōu)化射流結(jié)構(gòu)具有重要意義。

(a)截斷射流高速圖片(b)CFD數(shù)值模擬3D形態(tài)顯示

3.2射流結(jié)構(gòu)的動力特性分析

為深入了解截斷式脈沖射流頭部結(jié)構(gòu)的流場演化特點，將截斷式射流與未截斷射流的頭部流場結(jié)構(gòu)與速度場變化進行了對比分析。圖5(a)為未截斷圓柱射流噴出時的體積分?jǐn)?shù)模擬圖，與真實拍攝圖像(圖4(a)下)結(jié)構(gòu)上具有較高的一致性。該類型射流頭部呈現(xiàn)雨傘狀，射流前段直徑明顯大于從噴嘴噴出時的射流直徑，且射流頭部被向后滑移傘狀體包裹。截斷射流頭部結(jié)構(gòu)以射流軸線沿截斷切線方向平面呈面對稱結(jié)構(gòu)，圖5(a)為截斷射流頭部結(jié)構(gòu)以對稱面為剖面的體積分?jǐn)?shù)模擬圖，與未截斷射流具有較大的差異，相同時刻形成了長度明顯大于未截斷射流的尖體結(jié)構(gòu)，尖體完整長度約80 mm，最前端朝盤孔旋轉(zhuǎn)相反方向偏轉(zhuǎn)。圖6為兩種類型射流的速度場結(jié)構(gòu)對比，結(jié)合圖6分析可以發(fā)現(xiàn)，未截斷射流最前端直徑較大頭部結(jié)構(gòu)的速度較低，與其接觸的空氣區(qū)域在圓柱狀射流的高速推動下形成對稱的激波球面。根據(jù)流體動力學(xué)理論，當(dāng)射流馬赫數(shù)大于1時，射流前方氣體受壓縮形成激波，產(chǎn)生的音障效應(yīng)對射流阻力急劇增加。與圓柱射流相比，截斷式脈沖射流前端的尖體結(jié)構(gòu)與超音速戰(zhàn)斗機外形設(shè)計類似，在使射流在雷諾數(shù)明顯降低的情況下，可有效減弱或避免激波阻力的影響。圖中可以看出，截斷式脈沖射流前段向后滑移積聚現(xiàn)象不明顯，僅在射流偏轉(zhuǎn)一側(cè)形成了旋流，從而保證了尖體結(jié)構(gòu)300 m/s以上的有效打擊速度。

(a)圓柱射流(b)截斷射流

(a)圓柱射流(b)截斷射流

3.3截斷式脈沖射流破巖機制

脈沖射流破巖的主要影響因素有射流直徑、壓力、頻率、振幅等，本研究僅通過設(shè)計相同射流壓力和噴嘴直徑的一組連續(xù)圓柱射流與截斷式脈沖射流沖蝕花崗巖的對比實驗，實驗結(jié)果從一定程度上反映了截斷式射流特殊射流結(jié)構(gòu)與形式所形成的破巖機制。圖7是截斷式射流沖蝕孔形狀與連續(xù)射流的對比，連續(xù)射流沖蝕孔形狀為較為規(guī)則的圓柱形，直徑為13 mm，截斷式脈沖射流沖蝕孔截面為不規(guī)則下玄月型，弦長為16 mm，弦中點至圓弧的最短距離約為7 mm。不同的沖蝕孔形狀表明，射流頭部作為水錘壓力主要作用部分，是射流對巖石進行局部破碎的主體，決定了沖蝕孔的形狀。圖8為兩種射流在相同沖蝕時間內(nèi)沖蝕深度和沖蝕量的情況，其中沖蝕深度A和沖蝕體積A代表截斷式射流，沖蝕深度B和沖蝕體積B代表圓柱射流。截斷式脈沖射流與連續(xù)圓柱射流在累計沖蝕時間50 s之前均具有較高的遞增趨勢，最初5 s內(nèi)兩種射流的沖蝕深度與沖蝕體積近似相等，但截斷式脈沖射流的沖蝕深度與沖蝕體積的增加速率約為連續(xù)射流的2.8倍。累計沖蝕時間約50 s時，截斷式脈沖射流的沖蝕深度和沖蝕體積約為67 mm和13 ml，分別為連續(xù)射流的2.2倍和2.8倍。50 s后，連續(xù)射流作用下的沖蝕孔深度與沖蝕量趨于穩(wěn)定，而截斷式脈沖射流直到總沖蝕時間5 min時，沖蝕深度和沖蝕量隨沖蝕時間的增加而近似線性增加，截至5 min的沖蝕時間，截斷式是脈沖射流的總沖蝕深度為130 mm，沖蝕量為24 ml。以上沖蝕實驗表明：截斷式脈沖與連續(xù)圓柱射流雖沖蝕孔洞形狀不同，但截面積相當(dāng)；截斷式脈沖射流可降低空氣激波阻力對前端射流速度的影響，同時可避免在沖蝕孔內(nèi)所形成的水墊效應(yīng)，沖蝕能力較連續(xù)射流相比沖蝕能力大大提高。

圖7 連續(xù)射流與截斷式脈沖射流沖蝕花崗巖對比圖Fig.7 The erosion cavity on granite cube impacted by continuous cylindrical jet and interrupted pulse jet

圖8 未截斷與截斷射流破巖尺度對比Fig.8 The scale of erosion cavity on granite cube impacted by continuous cylindrical jet and interrupted pulse jet

4 結(jié) 論

基于VOF模型與3D動網(wǎng)格技術(shù)，模擬了單個截斷式射流脈沖速度場變化特性。利用室內(nèi)實驗，研究了截斷式脈沖射流的破巖效果。通過研究可以發(fā)現(xiàn)：

(1) 截斷式脈沖射流有助于在射流末端形成偏轉(zhuǎn)的尖體結(jié)構(gòu)，偏轉(zhuǎn)方向與孔盤轉(zhuǎn)動方向相反，偏轉(zhuǎn)的尖體結(jié)構(gòu)有助于射流破巖作用。

(2) 截斷式脈沖射流前端的尖體結(jié)構(gòu)與超音速戰(zhàn)斗機外形設(shè)計類似，在使射流在雷諾數(shù)明顯降低的情況下，可有效減弱或避免激波阻力的影響。

(3) 截斷式脈沖射流在形成過程中產(chǎn)生了瓦狀立體結(jié)構(gòu)，該瓦狀形態(tài)對射流的沖蝕效力產(chǎn)生一定的影響。因此，為了充分發(fā)揮截斷式脈沖射流的沖蝕破碎潛能，應(yīng)確定最優(yōu)靶距范圍。

(4) 在沖擊靶物時呈現(xiàn)出非軸對稱流動形態(tài)，這與破硬巖實驗中得到的沖蝕坑形態(tài)相吻合。

(5) 與常規(guī)圓柱射流相比，截斷式脈沖射流能夠大幅度提升破巖速度，增大破巖體積。

[1] 徐小荷, 余靜. 巖石破碎學(xué)[M]. 北京：煤炭工業(yè)出版社, 1984.

[2] 李曉紅，盧義玉，向文英．水射流技術(shù)及在礦業(yè)工程中的應(yīng)用[M]．重慶：重慶大學(xué)出版社，2007．

[3] DEHKHODA S, HOOD M. The internal failure of rock samples subjected to pulsed water jet impacts[J]. International Journal of Rock Mechanics and Mining Sciences, 2014, 66: 91-96.

[4] DEHKHODA S, HOOD M. An experimental study of surface and sub-surface damage in pulsed water-jet breakage of rocks[J]. International Journal of Rock Mechanics and Mining Sciences, 2013, 63: 138-147.

[5] 司鵠,薛永志. 基于SPH算法的脈沖射流破巖應(yīng)力波效應(yīng)數(shù)值分析[J].振動與沖擊,2016,35(5):146-152.

SI Hu, XUE Yongzhi. Numerical analysis for stress wave effects of rock broken under pulse jet[J].Journal of Vibration and Shock, 2016,35(5):146-152.

[6] 王樂勤,焦磊,徐如良,等.脈沖射流作用下駐點壓力特性的試驗研究[J].工程熱物理學(xué)報,2005, 26(1):69-71.

WANG Leqin, JIAO Lei, XU Ruliang, et al. Experimental study on stagnation pressure of pulse jet[J]. Journal of Engineering Thermophysics, 2005, 26(1): 69-71.

[7] MA Y, ZHU D Z, RAJARATNAM N, et al. Experimental study of the breakup of a free-falling turbulent water jet in air[J]. Journal of Hydraulic Engineering, 2016, 142(10): 06016014.

[8] 陸朝暉,盧義玉,夏彬偉,等.沖擊擠壓式脈沖射流動力特性數(shù)值模擬[J].中國石油大學(xué)學(xué)報(自然科學(xué)版),2013,37(4):104-108.

LU Zhaohui, LU Yiyu, XIA Binwei, et al. Numerical simulation on hydrodynamic characteristics of percussion pulsed jet[J]. Journal of China University of Petroleum, 2013,37(4):104-108.

[9] ANGLANI F, BARRY J, DEKKERS W, et al. CFD modelling of a water-jet cleaning process for concentrated solar thermal (CST) systems[C]//3rd Southern African Solar Energy Conference.[S.l.]:SASEC, 2015.

[10] LI D, KANG Y, DING X, et al. Effects of area discontinuity at nozzle inlet on the characteristics of high speed self-excited oscillation pulsed waterjets[J]. Experimental Thermal and Fluid Science, 2016, 29(4): 1-2.

[11] SRINIVASAN V, SALAZAR A J, SAITO K. Modeling the disintegration of modulated liquid jets using volume-of-fluid (VOF) methodology[J]. Applied Mathematical Modelling, 2011, 35(8):3710-3730.

[12] 汪劍鋒. 水下超聲速氣體射流噴射過程及周圍復(fù)雜流場的數(shù)值計算[D].杭州:浙江理工大學(xué),2014.

[13] PAN Z, WEIBEL J A, GARIMELLA S V. Spurious current suppression in VOF-CSF simulation of slug flow through small channels[J]. Numerical Heat Transfer, Part A: Applications, 2015, 67(1): 1-12.

Numericalsimulationandanalysisontheflowfieldstructureandhardrockerosionpotentialofadisc-slottedpulsewaterjet

LU Zhaohui1,2,3，LU Yiyu2，Michael Hood3，PAN Linhua1，HE Pei1

(1.Ministry of Land and Resources Key Laboratory of Shale Gas Exploration, Chongqing Institute of Geology and Mineral Resources, Chongqing 400042, China; 2.State Key Laboratory of Coal Mine Disaster Dynamics and Control, Chongqing University, Chongqing 400044, China; 3.Cooperative Research Centre of Mining, The University of Queensland, Brisbane 4072, Australia)

A disc-slotted pulse water jet is a potential tool to break hard rock due to its special loading styles, such as the water hammer pressure，ultra-speed lateral jetting and pulse dynamic load etc. Aiming at investigating the flow filed of the jet, a two-phase-flow transient computational model, matching with the geometry and motion of the interrupted water jet generating device, was established to simulate the dynamic evolution and characteristics of a single pulsation within 100 mm standoff based on the volume of fluid (VOF) model and dynamic mesh theory. The results show that at the head of the pulsed jet forms a deflective slug structure which is consistent with the result from high-speed photography experiments. The slug head velocity is lower than that at the jet outlet and the jet turbulence is mainly distributed over boundary layers between jet and air and at the deflective side. Tile-shaped stereo-structure is yielded during the pulse formation process. It presents non-axisymmetric flow pattern when impacting target, which consists with the irregular erosion cavity obtained by hard rock fragmentation experiments.

disc-slotted pulse water jet; flow field structure; dynamic mesh; hard rock erosion

TH212；TH213.3

A

10.13465/j.cnki.jvs.2017.19.028

國家自然科學(xué)基金(51304258)；重慶市自然科學(xué)基金重點項目(CSTC2013JJB90005)

2016-05-09 修改稿收到日期：2016-07-29

陸朝暉 男，博士，高級工程師，1984年生