半徑比對氣液噴射器吸氣量的影響研究

， ，，永軍

(1.青島科技大學 機電工程學院，山東 青島 266061；2.甘肅藍科石化高新裝備股份有限公司，甘肅 蘭州 730070；3.長慶油田分公司 第一采氣廠，陜西 靖邊 718500)

氣液噴射器是石油化工等行業(yè)中廣泛使用的混合、傳質(zhì)和反應設(shè)備，其工作原理是利用高速流動的工作流體去卷吸引射流體，實現(xiàn)流體之間的能量交換和混合傳質(zhì)，具有無轉(zhuǎn)動部件、傳質(zhì)效率高、密封性好、結(jié)構(gòu)簡單、維護費用低、連續(xù)性好以及安全性高等優(yōu)點[1]。

近年來，眾多學者運用實驗和數(shù)值方法對氣液噴射器進行了相關(guān)研究。耿利紅等[2]利用數(shù)值模擬的方法研究了噴射器幾何結(jié)構(gòu)對壓縮/噴射制冷循環(huán)性能的影響，并與傳統(tǒng)壓縮制冷循環(huán)的性能進行比較，發(fā)現(xiàn)噴射器存在一個最優(yōu)面積比。Kandakure等[3]的研究表明，混合段長徑比為0、混合段與噴嘴面積比為4時，噴射器能夠取得最大的吸氣量。Li等[4]以吸氣量為評價指標對噴射器混合段長徑比進行模擬研究，得出氣氣噴射器最優(yōu)長徑比為5～7，氣液噴射器最優(yōu)長徑比為1～2。Zhang等[5]的模擬研究顯示，噴嘴位置位于距離混合段入口6 mm時，噴射器能夠得到最大的吸氣量和最高的效率。Sharma等[6]利用簡化的噴射器三維數(shù)值模擬模型，研究了噴射器結(jié)構(gòu)參數(shù)對其吸氣量的影響，其研究結(jié)果為吸氣室結(jié)構(gòu)尺寸及噴嘴與混合段間距對吸氣量影響較小，擴散角為5°時，取得最大吸氣量。秦敬軒等[7]利用數(shù)值模擬的方法對不同出入口條件下氣液噴射器的噴射性能進行了研究，并得到了最優(yōu)工況。姚云等[8]以水和CO2為物系，利用平面激光誘導熒光技術(shù)，分析了操作條件對氣液噴射器內(nèi)流型的影響。范軍領(lǐng)等[9]對應用在回路反應器中的氣液噴射器進行了總結(jié)介紹，并對噴射器未來的應用做了展望。李韶璞[10]依照噴射器的工作原理，提出了氣液噴射器吸氣量的估算公式，并探究了噴射器的結(jié)構(gòu)參數(shù)對其混合性能的影響。從這些研究中可以看出，因應用領(lǐng)域和應用物系等方面的不同，得到的最優(yōu)幾何尺寸往往有較大的差別，無法對實際化工生產(chǎn)中氣液噴射器的制造提供直接指導。而且現(xiàn)有的氣液噴射器研究多集中于氣液輸送和冷卻循環(huán)系統(tǒng)領(lǐng)域，在化工生產(chǎn)領(lǐng)域的研究還不充分[11-16]。

文中內(nèi)容為用于三氟丙基甲基二氯硅烷生產(chǎn)工藝的噴射回路反應器開發(fā)項目中的一部分，為了確定最優(yōu)氣液噴射器的結(jié)構(gòu)參數(shù)，先建立了氣液噴射器的三維模型，基于該模型并利用FLUENT軟件對混合室與噴嘴出口的半徑比對噴射器吸氣量的影響規(guī)律進行了研究，探討了噴嘴出口流速對最優(yōu)半徑比值的影響，旨在對今后氣液噴射器的工程應用和設(shè)計優(yōu)化提供理論指導。

1 氣液噴射器三維模型建立

1.1 數(shù)學模型

1.1.1控制方程

質(zhì)量守恒方程：

(1)

動量守恒方程：

(2)

式(1)～式(2)中，ρ為密度，kg/m3；u、v、w分別為x、y、z方向的速度，m/s；ui、uj分別為i方向和j方向速度分量，m/s；xi、xj分別代表i、j方向的坐標；μ為空氣動力黏度，Pa·s；k為湍動能；ε為耗散率；Cμ為經(jīng)驗常數(shù)；Si為動量源項，N/m3；p為壓力，MPa；δij為克羅內(nèi)克爾張量。

1.1.2湍流模型

采用標準k-ε方程模擬三維流場，k-ε方程由湍動能k方程和耗散率ε方程組成。

k方程：

(3)

ε方程：

(4)

(5)

(6)

(7)

(8)

(9)

式(3)～式(9)中，Gk為平均速度梯度引起的湍動能的產(chǎn)生項，Gb為由于浮力引起的湍動能的生成項；Pri為湍動普朗特準數(shù)；μt為湍流黏度，Pa·s；β為熱膨脹系數(shù)，K-1；gi為重力加速度在i方向的分量，m/s2；T為溫度，K；YM為可壓縮流中脈動擴張項；a為導溫系數(shù)，m2/s；C1ε=1.44、C2ε=1.92、C3ε=0.09、σε=1.0、σk=1.3、Pri=0.85。

1.1.3半徑比m

半徑比m指氣液噴射器中混合室半徑與噴嘴出口半徑的比值，是影響噴射器性能的重要結(jié)構(gòu)參數(shù)之一。

(10)

式中，R為氣液噴射器中混合室的半徑，r為噴嘴出口半徑，mm。

1.2 物理模型

氣液噴射器的主要結(jié)構(gòu)包括噴嘴、吸氣室、混合室和擴散室，其結(jié)構(gòu)示意圖見圖1。

圖1 氣液噴射器結(jié)構(gòu)示圖

為了更好地對氣液噴射器吸氣量影響因素進行研究，對氣液噴射器進行了簡化處理，所用初始模型具體結(jié)構(gòu)尺寸見圖2，然后通過Solidworks軟件進行三維建模。

圖2 氣液噴射器初始結(jié)構(gòu)簡圖

1.3 網(wǎng)格劃分和邊界條件

1.3.1網(wǎng)格劃分

利用ICEM CFD軟件對氣液噴射器的三維模型進行網(wǎng)格劃分，網(wǎng)格類型為六面體結(jié)構(gòu)網(wǎng)格，具體見圖3。

圖3 氣液噴射器網(wǎng)格劃分

然后通過不斷增加網(wǎng)格數(shù)量對氣液噴射器進行網(wǎng)格加密，分別建立了4組網(wǎng)格，其網(wǎng)格數(shù)量分別為536 222、1 007 580、1 665 254和2 400 168。選取氣相入口體積流量進行網(wǎng)格無關(guān)性驗證，1 665 254網(wǎng)格與2 400 168網(wǎng)格的氣相入口體積流量相差僅為0.94%。

綜合考慮模擬結(jié)果的準確性以及計算的時間成本，最終選擇1 665 254作為模擬氣液噴射器模型的計算網(wǎng)格數(shù)量。

1.3.2邊界條件

液相入口為速度入口邊界。氣相入口為壓力入口邊界，壓力為0.98 MPa?；旌狭黧w出口為壓力出口邊界，壓力為1 MPa。壁面條件為光滑無滑移邊界條件。

1.4 物性參數(shù)

本文的模擬是基于甲基二氯硅烷(CH4Cl2Si)與三氟丙烯(C3H3F3)反應生產(chǎn)三氟丙基甲基二氯硅烷的生產(chǎn)工藝，工作流體為液相的甲基二氯硅烷，引射流體為氣相的三氟丙烯，反應溫度為100 ℃，反應壓力1 MPa，此時氣液兩相流體的部分物性參數(shù)見表1。

表1 兩種流體部分物性參數(shù)

2 氣液噴射器模擬過程

2.1 求解方法及收斂精度設(shè)置

定義氣液噴射器的模型為穩(wěn)定、壓力基以及隱式求解器。

采用基于壓力的離散求解法，利用SIMPLE算法，選擇Mixture模型且無滑移及標準k-ε湍流方程，壓力項求解采用PRESTO格式。初始化時，動量項、湍動能項和耗散率項求解采用二階迎風格式，收斂精度設(shè)為10-5。

2.2 模擬方案

為研究混合室與噴嘴出口的半徑比對氣液噴射器吸氣量的影響規(guī)律，在初始氣液噴射器模型基礎(chǔ)上，保持氣液噴射器的其它結(jié)構(gòu)尺寸不變，僅通過改變混合室半徑實現(xiàn)半徑比的變化，而且在不同液相入口流速下對噴嘴出口半徑為8.0 mm、10.0 mm、11.5 mm的噴射器進行了模擬，具體模擬計算的參數(shù)條件見表2。

表2 氣液噴射器模擬計算參數(shù)條件

3 氣液噴射器模擬結(jié)果與討論

3.1 初步模擬結(jié)果

氣液噴射器壓力分布云圖、速度分布云圖和液相體積分率分布云圖分別見圖4～圖6。

由圖4～圖6所示的分布云圖可以知道，氣液噴射器的壓力和速度變化趨勢總體上同步， 工作流體流經(jīng)噴嘴，壓力能轉(zhuǎn)化為動能，在較短位置上產(chǎn)生較大的速度梯度，在噴嘴出口處形成高速低壓區(qū)，為卷吸氣相提供條件。氣相在吸氣室進、出口壓差的作用下，被壓入混合室。在混合室內(nèi)，氣液進行混合，純氣相與純液相的體積逐漸減小，而氣液混合相的體積逐漸增大，出現(xiàn)壓力最小值和速度最大值。最后混合流體經(jīng)過擴散室，動能轉(zhuǎn)化為壓力能，混合流體壓力升高，流速緩慢降低，完成氣液兩相混合及噴射過程，得到的初步模擬結(jié)果符合噴射器的工作機理。

圖4 氣液噴射器壓力分布云圖

圖5 氣液噴射器速度分布云圖

圖6 氣液噴射器液相體積分率分布云圖

3.2 半徑比對噴射器吸氣量影響規(guī)律

對利用FLUENT軟件初步模擬所得的各項數(shù)據(jù)進行歸納整理，得到了噴嘴出口半徑r為8.0 mm、10.0 mm和11.5 mm時，在不同噴嘴出口流速v下，氣液噴射器吸氣量隨半徑比的變化曲線，分別見圖7～圖9。

圖7 r=8.0 mm 時氣液噴射器吸氣量隨半徑比變化曲線

圖8 r=10.0 mm時氣液噴射器吸氣量隨半徑比變化曲線

圖9 r=11.5 mm時氣液噴射器吸氣量隨半徑比變化曲線

由圖7～圖9所示的變化曲線可以知道，隨著半徑比m的增加，氣液噴射器的吸氣量呈先增大后減小趨勢，而且對于同一氣液噴射器結(jié)構(gòu)，當噴嘴出口流速不同時，得到最大吸氣量時的半徑比也不相同。這是由于半徑比增加不僅增大了混合室的直徑，還間接增大了吸氣室出口高度。吸氣室出口高度的增大會增大氣相進入混合室的流通面積，進而增大氣相的體積流量，有利于吸氣量的增加。而混合室直徑的增大，會增大混合室內(nèi)流體的流通面積，進而降低混合室內(nèi)的流體流速，使得混合室內(nèi)最低壓力得到提高?；旌鲜覂?nèi)最低壓力的升高會減小氣相通過吸氣室進入混合室的動力，不利于吸氣量的增加。在半徑比增加的初始階段，吸氣室出口高度增大對吸氣的促進作用占主導，吸氣量逐漸增大。而當吸氣量逐漸達到最大值以后，繼續(xù)增加半徑比比值，混合室半徑增大對吸氣的不利作用占主動，吸氣量開始降低。

噴嘴出口流速v為10 m/s、13 m/s和20 m/s時，不同噴嘴出口半徑噴射器的吸氣量隨半徑比m的變化曲線分別見圖10～圖12。

圖10 v=10 m/s時氣液噴射器吸氣量隨半徑比變化曲線

圖11 v=13 m/s時氣液噴射器吸氣量隨半徑比變化曲線

圖12 v=20 m/s時氣液噴射器吸氣量隨半徑比變化曲線

由圖10～圖12可知，對出口流速相同而噴嘴出口半徑不同的噴射器，取得最大吸氣量時的半徑比相近，且該比值隨著噴嘴出口流速的增加而逐漸增大。噴嘴出口流速為10 m/s、13 m/s和20 m/s時的最優(yōu)半徑比值分別為1.5、1.6～1.8和1.9～2.0。這是由于噴嘴出口流速越小，混合室內(nèi)的流體越容易受到混合室半徑增大的影響，進而使得吸氣量隨半徑比的變化越早地出現(xiàn)降低。

4 結(jié)語

利用CFD軟件對多組不同混合室與噴嘴出口半徑比的氣液噴射器進行了數(shù)值模擬，分析了半徑比對吸氣量的影響規(guī)律，結(jié)果表明：①隨著半徑比m的增加，氣液噴射器的吸氣量呈先增大后減小趨勢。對同一氣液噴射器結(jié)構(gòu)，當噴嘴出口流速不同時，得到最大吸氣量時的半徑比也不相同。②噴嘴出口處液相流速影響取得最大吸氣量時的半徑比，噴嘴出口流速越高，該比值越大，噴嘴出口流速為10 m/s、13 m/s、20 m/s時的最優(yōu)半徑比值為1.5、1.6～1.8、1.9～2.0。