張 沖，焦光偉，祝 韜，張少波

(1.陸軍勤務(wù)學(xué)院 油料系， 重慶 401331；2. 61716部隊， 福州 350003； 3. 78355部隊， 云南 楚雄 651224)

圓錐長直管段噴嘴內(nèi)流場仿真研究

張 沖1，焦光偉1，祝 韜2，張少波3

(1.陸軍勤務(wù)學(xué)院 油料系， 重慶 401331；2. 61716部隊， 福州 350003； 3. 78355部隊， 云南 楚雄 651224)

通過ANSYS FLUENT軟件，采用歐拉多相流模型對前混式磨料水射流內(nèi)流場進行建模與仿真。觀察不同參數(shù)下的仿真結(jié)果并進行比較，結(jié)果表明：壓力升高繼而導(dǎo)致噴嘴出口速度提高，通過簡化公式計算所得噴嘴出口速度不一定準確；噴嘴出口直徑變大會導(dǎo)致出口速度減小，且小于理論計算值；提供能量一定時，由于磨料體積分數(shù)增大會使得噴嘴出口速度減小，因此磨料體積分數(shù)不宜過大。

磨料水射流；內(nèi)流場；數(shù)值模擬

水射流技術(shù)是一種新型的加工技術(shù)，該技術(shù)通過泵提供動力為水加壓，經(jīng)過噴嘴收縮段加速后噴出，從而將泵提供的能量轉(zhuǎn)化為射流沖擊力作用于工件上。由于水射流技術(shù)屬于冷態(tài)加工，不同于傳統(tǒng)加工技術(shù)會產(chǎn)生明火或局部高溫，具有安全性好、效率高、節(jié)能環(huán)保等特點，因而目前在材料清洗、剝離、切割等工業(yè)作業(yè)中應(yīng)用廣泛[1-2]。

水射流根據(jù)射流形式可分為連續(xù)射流、脈沖射流和空化射流。其中連續(xù)射流又分液體射流、固-液射流和液體-氣體-固體射流[3]。其中固-液射流即磨料水射流。磨料水射流是水通過高壓泵增壓后，又經(jīng)過磨料罐與磨料顆粒進行混合，再通過噴嘴噴出。由于加入磨料，射流的切割能力進一步增強，在切割鋼材、破碎巖石及剝離涂層等作業(yè)上有著很好的效果。

不同于傳統(tǒng)流體力學(xué)的計算方法，磨料水射流涉及面更廣，計算更為復(fù)雜；因而需要借助計算機軟件結(jié)合計算流體力學(xué)(CFD)理論，通過多相流模型模擬固-液射流在空氣中的運動狀態(tài)，從而預(yù)測射流流場的各項參數(shù)，研究其特性。目前計算流體力學(xué)領(lǐng)域的商業(yè)軟件FLUENT已被ANSYS收購并集成于ANSYS 17.0軟件中。該軟件功能強大，在仿真計算射流方面具有一定優(yōu)勢，能較準確地預(yù)測射流流場特性。

本研究的前混式磨料射流系統(tǒng)不同于傳統(tǒng)后混式磨料射流系統(tǒng)，磨料是在高壓管路中與水混合，混合更加充分，只經(jīng)歷一次能量轉(zhuǎn)化，射流所需工作壓力較以前低很多，能量利用率更高[4]。此外，由于系統(tǒng)中磨料罐采用均勻供料，能保證在不同壓力下射流的磨料體積分數(shù)保持在相同水平。

1 物理建模和網(wǎng)格劃分

考慮到計算的準確性采用三維建模。在Workbench中通過Geometry組件建立三維模型，并通過Mesh組件中ICEM CFD軟件劃分網(wǎng)格，采用非結(jié)構(gòu)化網(wǎng)格進行劃分，控制網(wǎng)格質(zhì)量。如圖1所示，模型包括完整的噴嘴部分。另外，在建模前作如下假設(shè)：

1) 固液兩相均看作連續(xù)介質(zhì)，水、磨料、空氣不存在分層，可相互混合；

2) 固液兩相均為不可壓縮流體；

3) 磨料顆粒流簡化為一種固體流體。

圖1 仿真模型

圖1為模型總體和噴嘴結(jié)構(gòu)示意圖，其中外流場模型為長200 mm、底面直徑為30 mm的圓柱形區(qū)域。噴嘴模型總長76.2 mm，收縮段長5.39 mm，收縮角為40°，噴嘴出口直徑為0.8 mm，后續(xù)根據(jù)噴嘴直徑變化調(diào)節(jié)對應(yīng)參數(shù)。圖2為模型的網(wǎng)格劃分情況。設(shè)置CFD偏好和FLUENT求解器偏好，適當(dāng)調(diào)節(jié)網(wǎng)格密度和其他參數(shù)，確保網(wǎng)格質(zhì)量。

圖2 模型網(wǎng)格劃分

2 求解設(shè)置

2.1 模型選擇

求解中多相流模型選擇歐拉模型，湍流模型選擇RNGk-ε模型，其中控制方程[5]為：

質(zhì)量守恒方程：

?αkρk/?t+▽·(αkρkνk)=0

(1)

動量守恒方程：

?αkρkνk/?t+▽·(αkρkνkνk)=

(2)

能量守恒方程：

?αkρkhk/?t+▽·(αkρkνkhk)=

αk(?Pk/?t)+αkτk:▽νk-▽·αkqk+Hkl

(3)

RNGk-ε模型[6]使用以下的輸運方程來描述湍動能k和湍流耗散率ε：

k方程：

Gk+Gb-ρε-YM+Sk

(4)

ε方程：

(5)

2.2 求解設(shè)置

在FLUENT中設(shè)置求解參數(shù)。設(shè)置多相流模型為歐拉多相流模型(Eulerian multiphase flow model)，湍流模型選用RNGk-ε模型。設(shè)置材料(包括水、空氣和磨料)參數(shù)，其中水為主項，磨料和空氣分別為第2相和第3相。磨料密度為3 900 kg/m3，運動黏度為1.85×10-5m2/s，顆粒直徑為0.000 18 m；水密度998.2 kg/m3，20 ℃時運動黏度約為1.01×10-6m2/s。

進行邊界條件設(shè)置。因已知條件只有管路壓力，故噴嘴入口處設(shè)置為壓力入口，分別設(shè)置水、磨料和空氣的體積分數(shù)。噴嘴出口為壓力出口，出口表壓為0。對于多相流模擬，相間影響應(yīng)考慮[9-10]，因此設(shè)置固液間拖曳力為syamlal-obrien模型，并考慮虛擬質(zhì)量力對磨料顆粒的影響。設(shè)置求解方法為PC-SIMPLE，采用1階迎風(fēng)格式進行計算，適當(dāng)調(diào)節(jié)亞松弛因子。初始化流場并將流場空氣體積分數(shù)初始化為1，設(shè)置迭代時間步長為10-5s，迭代5 000步進行計算。

3 結(jié)果與分析

3.1 參數(shù)組合

根據(jù)文獻[11-13]，仿真模擬在設(shè)定了射流壓力、噴嘴直徑及磨料體積分數(shù)3個參數(shù)的不同組合條件下，共進行了8個不同參數(shù)組合的計算并進行對比，其中理論噴嘴出口速度是根據(jù)文獻[14]中式(7)計算所得，僅供仿真計算時參考。參數(shù)設(shè)定如表1所示。

表1 仿真參數(shù)組合

3.2 不同壓力下內(nèi)流場的仿真結(jié)果

取1～4組內(nèi)流場X-Y平面上的仿真結(jié)果進行分析，所得水和磨料X軸的軸向速度變化曲線如圖3所示。

從圖3可以看出：噴嘴內(nèi)的流體速度在通過收縮段后急劇升高，其大小和管路內(nèi)噴嘴入口壓力有關(guān)。顯然，入口壓力越大，噴嘴最后出口速度越大，這是由于入口處提供的壓能經(jīng)過噴嘴收縮段轉(zhuǎn)換為流體的動能，使水和磨料在收縮段有了一個加速的過程。水和磨料通過收縮段之后加速過程并不會馬上停止，兩相在直管段繼續(xù)混合加速。然而由于磨料和水的密度不同且黏度有差異，因而產(chǎn)生了速度滑移，水的速度較大，因此磨料在噴嘴中的加速過程是一個被動加速的過程，且磨料速度向水速度接近。與此同時，由于湍流作用，磨料顆粒在噴嘴內(nèi)不斷碰撞，造成了能量的部分損失，其速度并不能完全達到水的速度，因此將仿真計算中水的出口速度和理論計算的流體出口速度相比較，發(fā)現(xiàn)仿真計算的結(jié)果要比理論值低很多。經(jīng)過計算，1～4組中水的噴嘴出口速度分別是理論值的73.75%、73.06%、72.48%、73.77%，因此有些文獻中計算噴嘴出口速度時采用簡化的理論公式，但其準確性有待商榷。

3.3 不同出口直徑下內(nèi)流場仿真結(jié)果

取1、5、6組內(nèi)流場X-Y平面上的仿真結(jié)果進行分析，所得水和磨料的速度云圖和軸線速度變化曲線如圖4～6所示。

圖3 水和磨料內(nèi)流場軸線速度變化曲線

圖4 不同噴嘴直徑的噴嘴內(nèi)流場速度云圖

圖5 水和磨料內(nèi)流場軸線速度變化曲線

圖6 不同噴嘴直徑下內(nèi)流場流體的速度變化

通過內(nèi)流場速度云圖對比可發(fā)現(xiàn)：相同壓力時，保持噴嘴入口直徑和收縮角不變，由于噴嘴直徑的不同，因此射流流體的加速效果也不同；當(dāng)噴嘴出口直徑增大，射流加速變緩，最終達到的最大速度也會變小。

根據(jù)連續(xù)性方程：ρ1v1A1=ρ2v2A2知：管路內(nèi)入口參數(shù)不變時，出口速度和噴嘴出口直徑的平方成反比。實際仿真計算的出口速度數(shù)值和根據(jù)連續(xù)性方程計算的數(shù)值有一定誤差，這是由于兩相流動的相間影響、湍流運動以及仿真計算本來就可能存在的誤差導(dǎo)致的。尤其是當(dāng)噴嘴出口直徑增大后，會導(dǎo)致顆粒間碰撞幅度提高，影響射流集束性。從噴嘴直徑改變后的速度衰減情況看，噴嘴直徑不宜過大。

3.4 不同磨料體積分數(shù)下內(nèi)流場仿真結(jié)果

取1、7、8組內(nèi)流場X-Y平面上的仿真結(jié)果進行分析，所得水和磨料軸線速度變化曲線如圖7所示。

由圖7變化曲線可知：磨料體積分數(shù)的增大反而使得噴嘴中流體的速度減小。其主要原因是當(dāng)管路提供能量相同時，根據(jù)能量守恒定理，磨料體積濃度增大使得流體質(zhì)量增大，繼而導(dǎo)致在獲得相同動能的情況下流體的速度有所減小，而且由于顆粒體積濃度的增大導(dǎo)致顆粒間碰撞加劇，使得動能損失也會增大[15]。

圖7 水和磨料內(nèi)流場軸線速度變化曲線

4 結(jié)論

1) 隨著壓力的增大，噴嘴的出口速度相應(yīng)增大，且影響顯著，但通過簡化的理論公式計算得到的噴嘴出口速度不一定準確。

2) 對于圓錐長直線型噴嘴，增大噴嘴出口直徑會降低出口速度，且小于理論計算值。為了不影響射流的沖擊效果，噴嘴直徑不宜過大。

3) 在提供的能量一定的情況下，磨料體積分數(shù)的增大會導(dǎo)致流體出口速度減小，同時考慮到相間影響和能量損失，磨料體積分數(shù)不宜過大。

[1] 陳玉凡.高壓水射流清洗技術(shù)現(xiàn)狀及發(fā)展前景[J].中國設(shè)備工程,2013(2):6-8.

[2] 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(4):610-618.

[3] 沈忠厚.水射流理論與技術(shù)[M].濟南:石油大學(xué)出版社,1998:331-333.

[4] 管金發(fā),鄧松圣,伍建林,等.前混合磨料水射流噴嘴磨損規(guī)律的數(shù)值模擬[J].重慶理工大學(xué)學(xué)報(自然科學(xué)),2016,30(12):48-54.

[5] 陳雁,張鋼,何少煒,等.基于歐拉多相流模型的空氣噴涂噴霧流場模擬[J].后勤工程學(xué)院學(xué)報,2016,32(5):84-88.

[6] NIE B,HUI W,LEI L I,et al.Numerical investigation of the flow field inside and outside high-pressure abrasive waterjet nozzle[J].Procedia Engineering,2011,26(4):48-55.

[7] 曹國強,肖詠,李東明,等.磨料水射流噴嘴結(jié)構(gòu)參數(shù)對流場性能影響的仿真與研究[J].機床與液壓,2015,43(14):91-93.

[8] 董宗正,付必偉,郭燦,等.基于CFD的高壓水射流噴嘴流場仿真分析[J].石油和化工設(shè)備,2016,19(7):20-23.

[9] 王明波,王瑞和.湍流對磨料水射流中磨料顆粒受力的影響[J].中國機械工程,2012,23(8):972-974.

[10] 陸國勝,龔烈航,王強,等.前混合磨料水射流磨料顆粒加速機理分析[J].解放軍理工大學(xué)學(xué)報(自然科學(xué)版),2006,7(3):275-280.

[11] 張沖,焦光偉,管金發(fā),等.磨料水射流剝離3PE防腐層試驗研究[J].后勤工程學(xué)院學(xué)報,2017,33(2):40-43.

[12] 郭聯(lián)歡.磨料水射流切割預(yù)測與模擬技術(shù)研究[D].重慶：中國人民解放軍后勤工程學(xué)院，2016.

[13] 劉會霞,丁圣銀,王霄,等.水射流切割模型及其性能分析[J].農(nóng)業(yè)機械學(xué)報,2006,37(11):122-124，136.

[14] 屈長龍,王喜順.基于FLUENT的高壓水射流除銹的流場仿真及射流參數(shù)優(yōu)化[J].機械與電子,2016(2):24-27.

[15] LI W Y,ZHU H T,WANG J,et al.An investigation into the radial-mode abrasive waterjet turning process on high tensile steels[J].International Journal of Mechanical Sciences,2013,77(4):365-376.

InvestigationontheSimulationAbouttheExternalFlowFieldofConicalLongStraightPipeNozzle

ZHANG Chong1, JIAO Guangwei1, ZHU Tao2, Zhang Shaobo3

(1.Department of Petroleum Supply, Army Logistics University of PLA, Chongqing 401331, China; 2.The No.61716thTroop of PLA, Fuzhou 350003,China; 3.The No.78355thTroop of PLA, Chuxiong 402260, China)

A simulation about the internal flow field of pre-mixed abrasive water jet was conducted by ANSYS FLUENT and Eulerian multiphase flow model. Results in different parameters were compared and indicated that: with the increase of the pressure,velocity in nozzle outlet will increase accordingly.Velocity in nozzle outlet which calculated by simplified formula was proved not exact.With the increase of the diameter of nozzle outlet,velocity in nozzle outlet will decrease correspondingly and below the theoretical calculating value.When the energy was constant,velocity in nozzle outlet will decrease due to the increasing of the abrasive volume fraction.Thus,the abrasive volume fraction should not be in a high level.

abrasive water jet; internal flow field ; simulation

2017-08-17

中石化項目“川氣東送管道應(yīng)急搶修關(guān)鍵設(shè)備研制”(G8603-12-2S-0057*)

張沖(1993—)，男，河南洛陽人，碩士研究生，主要從事油氣儲存技術(shù)與裝備研究，E-mail:516484430@qq.com。

張沖，焦光偉，祝韜，等.圓錐長直管段噴嘴內(nèi)流場仿真研究[J].重慶理工大學(xué)學(xué)報(自然科學(xué))，2017(12):105-110.

formatZHANG Chong, JIAO Guangwei, ZHU Tao, et al.Investigation on the Simulation About the External Flow Field of Conical Long Straight Pipe Nozzle[J].Journal of Chongqing University of Technology(Natural Science)，2017(12):105-110.

10.3969/j.issn.1674-8425(z).2017.12.018

TP602

A

1674-8425(2017)12-0105-06

(責(zé)任編輯林 芳)