基于后混合式磨料水射流磨料顆粒運(yùn)動研究

潘崢正，萬慶豐，雷玉勇，蔣代君，陳忠敏

(西華大學(xué)機(jī)械工程與自動化學(xué)院，四川成都 610039)

基于后混合式磨料水射流磨料顆粒運(yùn)動研究

潘崢正，萬慶豐，雷玉勇，蔣代君，陳忠敏

(西華大學(xué)機(jī)械工程與自動化學(xué)院，四川成都 610039)

為了改善后混合式磨料水射流的切割性能，以固液兩相流理論為基礎(chǔ)，通過FLUENT模塊對噴嘴內(nèi)磨料的運(yùn)動進(jìn)行建模和數(shù)值分析。以聚焦管內(nèi)的單個(gè)顆粒為研究對象，優(yōu)化傳統(tǒng)的顆粒相控制模型，獲得磨料顆粒在噴嘴內(nèi)基本的運(yùn)動情況。研究表明，磨料顆粒吸入混合腔后通過各相之間的碰撞進(jìn)入高壓射流，形成液固兩相流。兩相流形成初期，磨料顆粒的速度在高壓水射流的攜帶作用下迅速上升。由于混合腔內(nèi)部結(jié)構(gòu)突變等原因，顆粒相速度浮動較大且處于非穩(wěn)態(tài)。當(dāng)顆粒通過收斂段進(jìn)入聚焦管后，顆粒相速度以指數(shù)形式逼近射流相速度并逐漸穩(wěn)定。但由于存在沿程能量損失等因素，顆粒相速度與射流相速度始終存在滑移，最終顆粒相速度仍略小射流相速度。

后混合式磨料水射流;磨料顆粒;數(shù)值模擬;FLUENT

磨料水射流 (Abrasive Water Jet，AWJ)，是將一定數(shù)量的磨料 (如碳化硅、石榴子石、金剛砂以及石英砂等砂類)與高壓水混合形成固液兩相射流［1－2］。高壓水的部分動能傳遞給磨料，改變了射流對靶體的作用方式，將高壓水射流對靶體的持續(xù)作用改變?yōu)槟チ蠈Π畜w的沖擊、磨削作用，提高了射流的品質(zhì)和工作效率。噴嘴作為后混合式磨料水射流系統(tǒng)的重要結(jié)構(gòu)部件，直接影響產(chǎn)品加工的質(zhì)量。通過研究混砂管內(nèi)磨料顆粒的運(yùn)動規(guī)律，可以有效地改善噴嘴的性能。文中針對后混合磨料水射流系統(tǒng)，以固液兩相流理論為基礎(chǔ)，對后混合式磨料水射流噴嘴內(nèi)磨料顆粒的運(yùn)動情況進(jìn)行建模，利用ANASYS軟件平臺下的FLUENT模塊進(jìn)行數(shù)值仿真，分析磨料顆粒在后混合式磨料水射流噴嘴內(nèi)的基本運(yùn)動規(guī)律。

1 物理模型和基本假設(shè)

后混合式磨料水射流噴嘴一般由水噴嘴、混合腔、聚焦管三部分結(jié)構(gòu)組成，如圖1所示。該試驗(yàn)只研究聚焦管內(nèi)磨料顆粒的運(yùn)動情況，采用了廣泛使用的錐形噴嘴模型。

圖1 后混合式磨料水射流噴嘴結(jié)構(gòu)示意圖

為了研究磨料顆粒的運(yùn)動情況，作以下兩個(gè)假設(shè):(1)磨料顆粒為相同質(zhì)量、體積的剛性等直徑小球，顆粒與流體的相對運(yùn)動遵循Stokes阻力定律; (2)固液兩相之間不存在質(zhì)量交換。對于流體作以下假設(shè):(1)水為連續(xù)介質(zhì);(2)水為不可壓縮流體。同時(shí)，對于固液兩相流做以下假設(shè):(1)兩相流和外界無熱量交換且溫度不變;(2)固液兩相流是定常的紊流［3］。

2 磨料顆粒在噴嘴內(nèi)的加速過程理論分析

在固液兩相流中，聚焦管內(nèi)軸線附近的單個(gè)磨料顆粒的受力情況，主要包括兩部分:(1)與流體－顆粒相對運(yùn)動無關(guān)的力;(2)與流體－顆粒相對運(yùn)動直接相關(guān)的力。當(dāng)磨料顆粒進(jìn)入噴嘴后，以牛頓第二定律為基礎(chǔ)的歐拉型兩相流運(yùn)動微分方程可表示為:

式中:下標(biāo)f代表流體;下標(biāo)p表顆粒;D為顆粒直徑;μ為流體黏性系數(shù);up為顆粒速度;uf為顆粒周圍流體速度;Fp為外應(yīng)力。

式 (1)中存在較多復(fù)雜的未知變量，需要進(jìn)行大量假設(shè)，給運(yùn)算帶來了極大的不便，且不適合針對具體模型的數(shù)學(xué)分析。研究表明，在噴嘴內(nèi)磨料顆粒所受到的慣性力、重力、壓差力和Stokes力對磨料顆粒的運(yùn)動有著重要的影響，必須予以考慮。Basset力、虛擬質(zhì)量力和Stokes力處于同一量級不可忽略。而流體的黏性極大地限制了磨料顆粒的旋轉(zhuǎn)，所以相對于磨料顆粒所受到的Stokes力和Mangus力可以忽略。同時(shí)，由于流速梯度沿徑向變化較緩慢，磨料顆粒所受到的Saffman力相對于Stokes力也可忽略不計(jì)［4－5］?；谌缟系募僭O(shè)分析，且不考慮因存在聚焦管薄壁而產(chǎn)生不均勻流場的影響，可以得到如下磨料顆粒運(yùn)動的控制方程:

考慮到壁面附近存在不均勻流場，對上式中平行于壁面方向的Stokes阻力項(xiàng)做以下修正:

計(jì)算表明，當(dāng)Y/D=2時(shí)，阻力為無界時(shí)阻力的1.3倍，將磨料顆粒在聚焦管內(nèi)所受到的定常阻力設(shè)為無界條件下的1.3倍。因而，Stokes阻力可表達(dá)為:

式中:CD為磨料顆粒阻力系數(shù)，與磨料顆粒的雷諾數(shù)Re直接相關(guān)，Re可通過下式計(jì)算:

由式 (6)可得Re＜＜1，可以忽略N－S方程中的慣性項(xiàng)，推導(dǎo)出低Re下的CD的計(jì)算公式:

對于十分稀疏的固液兩相流，流體速度梯度和壓強(qiáng)存在如下關(guān)系:

綜上所述，最終優(yōu)化的顆粒運(yùn)動控制方程可表達(dá)為:

其中等式右邊第一項(xiàng)為Stokes阻力，第二項(xiàng)為壓差力，第三項(xiàng)為附加質(zhì)量力，第四項(xiàng)為Basset力，第五項(xiàng)為外應(yīng)力。

由于難以得到相對精確的時(shí)間參數(shù)τ，導(dǎo)致求解精確的Basset力難度增加，但Basset力和Stokes阻力大小處于同一量綱，通過其比值可估算出相對準(zhǔn)確的Basset力［6－7］。

后混合式磨料水射流系統(tǒng)中，流體處于高速運(yùn)動狀態(tài)，而系統(tǒng)時(shí)間和相對速度弛豫時(shí)間不可知，由聲音在水中的傳播速度為1 500 m/s，可以估算出t/τ約為2。由圖2可知，該情況下 Basset力大小約為Stokes阻力的0.2倍左右，不妨設(shè) Basset力大小為Stokes阻力大小的0.2倍，則化簡后式 (9)可表示為

在高壓射流中，外力的作用的影響很小 (主要是重力)，因此可以忽略掉Fp項(xiàng)?？傻?/p>

將上式兩邊同除以μf，得到兩相之間的速度相對于時(shí)間參量的變化關(guān)系:

由式 (14)可以看出，如果高壓水和磨料顆粒進(jìn)入聚焦管時(shí)的初始速度已知，就能得到磨料顆粒速度的精確解。聚焦管內(nèi)固相和液相的速度比與時(shí)間參量滿足一定的指數(shù)函數(shù)關(guān)系，用MATLAB繪制其函數(shù)圖象(圖3所示)。

圖3 兩相速度與時(shí)間參量關(guān)系

圖2 FB/FS與t/τ的關(guān)系

由圖可知，聚焦管軸線附近磨料顆粒的速度將以指數(shù)形式逼近高壓水的速度，最終保持一致。但在實(shí)際情況下，磨料顆粒的速度并不會無限接近高壓水的速度，兩者之間存在一定速度滑移，在顆粒數(shù)學(xué)模型計(jì)算中沒有考慮沿程能量損失，是造成上述差異的主要原因。

3 磨料顆粒與流體的固液兩相流動仿真分析

后混合式磨料水射流系統(tǒng)主要依靠磨料自身重力和高壓水射流在混合腔中產(chǎn)生的負(fù)壓，將磨料卷吸進(jìn)入混合腔，并通過混合腔的混合和聚焦管的集中，最終從噴嘴口射出高速磨料射流［8］。與前混合式磨料水射流相比，后混合式磨料水射流的混合全過程都發(fā)生在混合腔內(nèi)，磨料所具有的初速度只來源于自身的重力和混合腔中高壓水產(chǎn)生負(fù)壓的吸力，其速度遠(yuǎn)小于高壓水的速度，導(dǎo)致了大多數(shù)磨料顆粒需要經(jīng)過多次碰撞，才能獲得進(jìn)入射流所需的動能［9－10］。

運(yùn)用FLUENT模塊對噴嘴進(jìn)行三維建模，如圖4所示。模型尺寸為:水噴嘴直徑φ0.1 mm，混合腔直徑φ6 mm，混合腔長度8 mm，收斂角30°，聚焦管直徑φ0.3 mm，聚焦管長度40 mm。采用MIXTURE模型仿真磨料在混合腔中的混合過程，基本參數(shù)見表1，仿真結(jié)果如圖5—8所示。

圖4 噴嘴三維模型圖

表1 基本參數(shù)設(shè)置

由圖5—6可知，磨料顆粒在混合腔中混合較為均勻，這是由于在水噴嘴入口處 (遠(yuǎn)離磨料入口附近)形成了較強(qiáng)的漩渦并產(chǎn)生負(fù)壓 (圖7—8所示)，促使磨料和高壓水混合。但混合腔腔壁附近磨料分布密度相對較大，這主要是因?yàn)榛旌锨惠S線處流體徑向速度較大，高壓水射流中心剛度較大，顆粒受到的阻力也相對較大，阻礙磨料顆粒進(jìn)入射流中心端。由于MIXTURE模型采用顆粒云方式運(yùn)算顆粒運(yùn)動，不計(jì)算單個(gè)顆粒之間相互的碰撞，導(dǎo)致磨料顆粒和高壓水混合的真實(shí)情況與仿真結(jié)果存在偏差，需要進(jìn)一步的試驗(yàn)驗(yàn)證。

圖5 磨料和水混合的體積分布

圖6 磨料顆粒密度分布

圖7 噴嘴內(nèi)壓力分布圖

圖8 速度矢量示意圖

圖9 軸線方向顆粒速度變化

在MIXTURE模型中計(jì)算離散相，跟蹤對象為混合腔軸線方向上的顆粒云，速度變化如圖9所示。磨料顆粒進(jìn)入混合腔后，在高壓水的攜帶作用下速度激增，但由于混合腔內(nèi)部結(jié)構(gòu)的突變性、收斂段錐面、顆粒之間的相互作用等因素的影響，導(dǎo)致顆粒始終處于非穩(wěn)態(tài)，顆粒紊動。直到顆粒進(jìn)入收斂段才逐漸趨于穩(wěn)態(tài)，進(jìn)入聚焦管后，速度基本保持不變。另外，顆粒大部分速度是在混合腔中通過高壓水的攜帶過程獲得，這并不表示顆粒在聚焦管中沒有加速過程，進(jìn)入聚焦管的顆粒，相互間的碰撞甚小，雖然其運(yùn)動狀態(tài)趨向穩(wěn)定，但仍然會以指數(shù)形式逼近液相速度。

4 結(jié)論

通過理論、數(shù)學(xué)模型及建模仿真分析，得到后混合式磨料水射流系統(tǒng)中磨料顆粒的運(yùn)動情況:

(1)磨料顆粒通過自重和高壓水產(chǎn)生的負(fù)壓作用被吸入混合腔，進(jìn)入混合腔后，磨料通過磨料和磨料、磨料和高壓水射流以及磨料和壁面之間的相互碰撞，獲得進(jìn)入高壓水射流所需的動能，完成磨料和高壓水的混合，最終形成固液兩相流。但由于射流中心的剛度較大，導(dǎo)致混合腔腔壁附近磨料分布密度較大。

(2)當(dāng)混合腔中的磨料進(jìn)入高壓水射流后，在高速水的攜帶作用下磨料速度急劇上升，但由于混合腔內(nèi)部突變性、錐面結(jié)構(gòu)等影響，磨料顆粒的速度會發(fā)生較大的紊動，直到顆粒進(jìn)入收斂段逐漸進(jìn)入穩(wěn)態(tài)。

(3)磨料顆粒進(jìn)入聚焦管后將以指數(shù)形式逼近高速水的速度，由于沿程能量損失等影響，磨料顆粒與高速水射流存在一定的速度滑移，即磨料顆粒的速度始終小于高壓水的速度。

［1］李志榮，封志明，郭宗環(huán).微磨料水射流和多晶硅切割的試驗(yàn)研究和模型分析［J］.礦山機(jī)械，2012，40(7):127－130.

［2］沈忠厚.水射流理論和技術(shù)［M］.北京:石油大學(xué)出版社.

［3］王明波，王瑞和.磨料水射流中磨料顆粒的受力分析［J］.中國石油大學(xué)學(xué)報(bào):自然科學(xué)版，2006，30(4):47－50.

［4］陸國盛，龔烈航，王強(qiáng)，等.前混合式磨料水射流磨料顆粒加速機(jī)理分析［J］.解放軍理工大學(xué)學(xué)報(bào):自然科學(xué)版，2006，7(3):275－280.

［5］萬繼偉，牛爭鳴，廖偉麗，等.后混式射流中多相流特性及加速分區(qū)理論研究［J］.水動力學(xué)研究與進(jìn)展，2012，27(5):546－550.

［6］王軍，于超，耿鵬飛.基于量綱分析的一種磨料水射流打孔模型［J］.機(jī)床與液壓，2012，40(9):49－51.

［7］熊佳，雷玉勇，楊志峰，等.基于FLUENT磨料水射流噴嘴內(nèi)流場可視化研究［J］.潤滑與密封，2008，33(6):51－54.

［8］陳琳，雷玉勇，郭宗環(huán)，等.基于FLUENT后混合式磨料水射流噴嘴內(nèi)流場數(shù)值模擬［J］.潤滑與密封，2012，37 (4):66－69.

［9］NANDURI Madhusarathi，TAGGART David G，KIM Thomas J.The Effects of System and Geometric Parameters on Abrasive Water Jet Nozzle Wear［J］.International Journal of Machine Tools＆Manufacture，2002，42(5):615－623.

［10］王榮娟，雷玉勇，蔣代君，等.基于FLUENT的水射流聚焦管集合參數(shù)優(yōu)化［J］.機(jī)床與液壓，2011，39(11):24－26.

Research on Motion of Abrasive in Post-mixed Abrasive Water Jet

PAN Zhengzheng，WAN Qingfeng，LEI Yuyong，JIANG Daijun，CHEN Zhongmin
(School of Mechanical Engineering＆Automation，Xihua University，Chengdu Sichuan 6100039，China)

In order to improve the cutting performance of the post-mixed abrasive water jet，the mathematical model of the abrasive motion in the focusing tube was established and numerically simulated on the basis of theory of solid-fluid two phase flow by using FLUENT software module.The traditional control model of particle phase was optimized by a researching object of single particle in the focusing tube，and the basic movement condition of abrasive particle in the nozzle was obtained.Research results show that the abrasive particles join the high-speed jet by the collision between each phase after they enter into the mixing chamber.Therefore two phase jet flow in the end is formed.Under the action of high speed water jet carrying，the speed of the particles will rise rapidly on initial forming stage of the two-phase flow.Due to causes of sudden change of internal structure in mixing chamber，phase velocity of particle has a big difference and is unsteady.After the particles enter into the focus tube through convergence section，the phase velocity of particle will exponentially approach the phase velocity of liquid jet and it will gradually be in steady.But because of causes of linear energy loss along the path，phase velocity of particle has always slide away relatively to that of the liquid jet.The velocity of abrasive particles is slightly slower than liquid jet speed after all.

Post-mixed abrasive water jet;Abrasive particle;Mathematical simulation;FLUENT

TP601

A

1001－3881(2014)9－109－4

10.3969/j.issn.1001－3881.2014.09.030

2013－04－25

四川省科技廳項(xiàng)目資助 (2011JYZ017);西華大學(xué)研究生創(chuàng)新基金資助

潘崢正 (1988—)，男，碩士研究生，主要研究方向?yàn)樗淞魈胤N精密加工技術(shù)。E－mail:276010304@qq.com。