孟慶元 王長(zhǎng)軍 羅立臣

（中海油能源發(fā)展裝備技術(shù)有限公司 天津300452）

0 引 言

攪拌器在石油、化工等行業(yè)應(yīng)用廣泛，一般分為機(jī)械式與流體式。根據(jù)不同的應(yīng)用場(chǎng)景和工況，會(huì)選取不同類型的攪拌器。機(jī)械攪拌器的工作原理是利用電機(jī)帶動(dòng)攪拌軸旋轉(zhuǎn)向釜內(nèi)介質(zhì)輸入機(jī)械能從而達(dá)到攪拌目的。海洋石油平臺(tái)絕大部分采用機(jī)械式攪拌器，其優(yōu)點(diǎn)是具有較大的攪拌能量、流體微團(tuán)分散較好、溶解效果好，適用于中低黏度的流體混合；但是海洋石油平臺(tái)空間小，設(shè)備的安裝受到危險(xiǎn)區(qū)域的影響，安全可靠性要求較高。

氣體攪拌是流體攪拌中最常見的一種形式，通過把氣體通入攪拌液內(nèi)形成氣液混合物，由于氣體密度小且形成的氣泡向上運(yùn)動(dòng)，從而對(duì)液體進(jìn)行攪拌。與機(jī)械攪拌相比，氣體攪拌具有運(yùn)行穩(wěn)定、結(jié)構(gòu)簡(jiǎn)單和能耗低等優(yōu)點(diǎn)，但有時(shí)過高的氣速使傳質(zhì)系數(shù)下降而不利于傳質(zhì)。有些行業(yè)為了緩解單一形式攪拌器帶來的缺陷，采取使用氣體分布器通入氣體，并在機(jī)械攪拌的作用下加強(qiáng)氣液傳質(zhì)效果。

近期，美國(guó)DERRICK 公司針對(duì)海洋泥漿混合系統(tǒng)研發(fā)了1 種化學(xué)混合罐，該混合罐中的攪拌器就是利用氣體推動(dòng)攪拌器旋轉(zhuǎn)，達(dá)到機(jī)械攪拌和流體攪拌雙重效果。

本文探索新的攪拌形式，設(shè)計(jì)并提出一種以壓力氣體為動(dòng)力的具有雙重特性的攪拌器，不僅氣體直接作用于介質(zhì)進(jìn)行氣體攪拌，而且氣體噴出后的反作用力也可推動(dòng)槳葉轉(zhuǎn)動(dòng)進(jìn)行機(jī)械攪拌。為了進(jìn)一步了解新型攪拌器的性能，文中利用CFD數(shù)值模擬和激光多普勒測(cè)速法（ Laser Doppler Velocimeter，LDV）對(duì)攪拌流場(chǎng)進(jìn)行了測(cè)試驗(yàn)證。

1 新型槳葉的結(jié)構(gòu)、工作原理及優(yōu)點(diǎn)

圖1所示是新型氣動(dòng)攪拌器。其攪拌軸為空心通氣管，空心管上端固定于釜體上，不跟隨槳葉旋轉(zhuǎn)。空心管下端開有若干孔供氣體輸出，管外部套有一個(gè)“C”型圓環(huán)（內(nèi)部結(jié)構(gòu)見圖2）。

圖1 新型氣動(dòng)攪拌器模型圖

圖2 槳葉內(nèi)部剖面圖

“C”型圓環(huán)上連接著空心槳葉，葉片上設(shè)置有氣體分布孔。在“C”型圓環(huán)上下設(shè)置有限位板（見圖1），防止運(yùn)行過程中上下移動(dòng)，為減少限位板與圓環(huán)間的氣體泄漏和摩擦阻力，之間粘有帶自潤(rùn)滑作用的墊片。

如圖1 所示，空氣壓縮機(jī)產(chǎn)生的氣體首先經(jīng)由空心通氣管 (攪拌軸)輸入，然后經(jīng)過通氣管下端的方孔進(jìn)入“C”型圓環(huán)，最后進(jìn)入與“C”型圓環(huán)相連接的槳葉，由槳葉上的氣孔排出。氣體排出后直接作用于介質(zhì)進(jìn)行攪拌，而產(chǎn)生的反推力推動(dòng)槳葉轉(zhuǎn)動(dòng)，最終達(dá)到氣體攪拌和機(jī)械攪拌的雙重作用。

新型攪拌器因?yàn)槭菤怏w作為驅(qū)動(dòng)力，所以只需要1 根氣管線與之連接，精簡(jiǎn)了以往電機(jī)驅(qū)動(dòng)攪拌器所需要的電纜、接線盒等附屬設(shè)備，更無需因?yàn)榭紤]危險(xiǎn)區(qū)因素而提高設(shè)備的防爆等級(jí)。因?yàn)橥夤芎附釉谕搀w上，采用靜密封結(jié)構(gòu)取代以往攪拌軸的動(dòng)密封結(jié)構(gòu)，不僅有效精簡(jiǎn)了設(shè)備結(jié)構(gòu)、降低制造成本，而且氣體帶動(dòng)攪拌器工作可以減少電機(jī)驅(qū)動(dòng)所產(chǎn)生的熱量以及潤(rùn)滑油泄漏和污染。

2 攪拌器流場(chǎng)的數(shù)值模擬

以攪拌釜中點(diǎn)為坐標(biāo)原點(diǎn)建模，豎直方向?yàn)檩S。模擬所采用的攪拌釜為平底圓柱筒，釜體內(nèi)直徑為300 mm、液面高度300 mm，攪拌槳中心面位于釜體的幾何中心，模擬介質(zhì)為水。

2.1 網(wǎng)格劃分

本文利用Gambit軟件進(jìn)行建模。由于實(shí)驗(yàn)?zāi)Ｐ徒Y(jié)構(gòu)不規(guī)則不利于劃分網(wǎng)格，建模時(shí)在不影響效果的前提下對(duì)葉片進(jìn)行了簡(jiǎn)化，采用多重考系法（Multi-Reference Frame，MRF）解決了靜止部件和轉(zhuǎn)動(dòng)的攪拌槳之間的相對(duì)運(yùn)動(dòng)問題，整個(gè)計(jì)算域分為動(dòng)區(qū)域和靜區(qū)域。槳葉及其附近流體區(qū)域設(shè)置為動(dòng)區(qū)域，其他流體區(qū)域設(shè)置為靜區(qū)域。通過交界面（Interface）兩部分區(qū)域?qū)崿F(xiàn)動(dòng)量、能量的交換。網(wǎng)格劃分采用適應(yīng)性較強(qiáng)的四面體非結(jié)構(gòu)化網(wǎng)格，對(duì)動(dòng)區(qū)域網(wǎng)格的劃分進(jìn)行加密，分別生成約70萬個(gè)動(dòng)網(wǎng)格和90萬個(gè)靜網(wǎng)格，見圖3。

圖3 網(wǎng)格劃分圖

2.2 模擬方法

計(jì)算采用Fluent軟件。本文主要研究了氣體和液體混合的二相流，且離散相為氣相容比較低的泡狀流，計(jì)算域分布較廣，所以選擇計(jì)算模型為Mixture模型。湍流模型選用標(biāo)準(zhǔn)模型，主相為水、第二相為空氣，流動(dòng)狀態(tài)設(shè)置為非定常流動(dòng)，壓力-速度耦合采用SIMPLE算法，對(duì)流項(xiàng)的離散使用二階迎風(fēng)差分格式，迭代時(shí)間步長(zhǎng)為1 ms，設(shè)置殘差迭代到10時(shí)為收斂。設(shè)置邊界條件，基于前期試驗(yàn)運(yùn)行的基礎(chǔ)，設(shè)置動(dòng)區(qū)域以60 r/min轉(zhuǎn)動(dòng)，氣體出口設(shè)置為相對(duì)于槳葉面流速為15 m/s的速度出口，頂部液面為自由流動(dòng)液面。

2.3 數(shù)值模擬結(jié)果分析

圖4為=0面的軸向速度圖，圖5為=0面放大后的軸向矢量速度圖。

圖4 軸向速度圖

圖5 軸向速度矢量圖

圖4中：槳葉以上區(qū)域除了靠近壁面小部分呈現(xiàn)黃色和紅色，其他大部分呈現(xiàn)濃綠色，槳葉斜上方不遠(yuǎn)處有小部分的藍(lán)色；槳葉以下區(qū)域除了靠近壁面小部分和正對(duì)軸的下方呈現(xiàn)藍(lán)色，其他大部分呈現(xiàn)黃綠色，2個(gè)肋部黃色比較明顯。由圖譜對(duì)比可得出結(jié)論：槳葉以上的區(qū)域，靠近壁面區(qū)域的軸向速度高于其他區(qū)域的軸向速度；槳葉下方區(qū)域，靠近壁面和軸正下方的區(qū)域軸向速度不明顯，而2個(gè)肋部軸向速度最明顯。整體來看，槳葉下方的軸向速度優(yōu)于槳葉上方的軸向速度，分析其原因是槳葉上方區(qū)域由于氣泡的螺旋上升和離散對(duì)流場(chǎng)的循環(huán)流動(dòng)起到不定向的擾動(dòng)。

觀察圖5整體的流場(chǎng)分布，槳葉所在平面主要以徑向流動(dòng)為主，流動(dòng)到壁面處分為向上和向下的2股流體，將顏色對(duì)照?qǐng)D譜可知向上的流體速度高于向下流動(dòng)的流體；槳葉下方呈現(xiàn)2個(gè)類似“C”型的流場(chǎng)，“C”型流場(chǎng)在靠近壁面處與徑向流向下的支流相遇、相互作用，速度相互抵消形成了低速的亂流，呈現(xiàn)出藍(lán)色。槳葉上方的流場(chǎng)主要是氣泡螺旋上升和離散所造成的湍流，方向沒有規(guī)律且速度較慢；而靠近壁面區(qū)域的淺黃色主要是徑向流向上的1股分流形成的，速度稍高于其他區(qū)域。

圖6是模擬條件完全相同的普通機(jī)械式攪拌器（45°折葉槳）=0面軸向速度模擬圖，圖7是與圖6對(duì)應(yīng)的速度矢量圖。

圖6 軸向速度圖

對(duì)比圖6和圖4后發(fā)現(xiàn)，部分區(qū)域顏色有很大差別。圖6中的槳葉上方區(qū)域藍(lán)色部分較多，軸與槳葉的連接處藍(lán)色部分最明顯，軸向速度比圖4相同區(qū)域要差，但是靠近壁面區(qū)域黃色和紅色面積較大，軸向速度優(yōu)于圖4相應(yīng)區(qū)域。槳葉下方區(qū)域與圖4有很大的反差，大面積連續(xù)黃色和紅色區(qū)域聚集于軸下方，除靠近壁面小部分區(qū)域外，其他區(qū)域軸向速度良好。

圖5和圖7進(jìn)行流場(chǎng)對(duì)比后發(fā)現(xiàn)，無氣孔的攪拌器呈現(xiàn)4個(gè)明顯的“O”型循環(huán)流場(chǎng),高速區(qū)域集中在軸正下方，流場(chǎng)的連續(xù)性和規(guī)律性明顯。

圖7 軸向速度矢量圖

綜上觀察分析得出：從顏色上看，新型攪拌器的整體軸向速度分布沒有45°折葉攪拌器好，但是分別對(duì)比左側(cè)的圖譜可以發(fā)現(xiàn)：新型攪拌器軸向速度較高，并且新型攪拌器槳葉以上區(qū)域的速度分布比較均勻，在氣泡的攪拌作用下低速區(qū)較少，形成穩(wěn)定的湍流，靠近壁面的高速區(qū)域小，防止高速流體集中在壁面區(qū)域影響攪拌效率；而在槳葉以下區(qū)域，新型攪拌器注重2個(gè)肋部的攪拌，攪拌作用范圍較廣，在氣泡的作用下流場(chǎng)略微不規(guī)律。

3 LDV 實(shí)驗(yàn)測(cè)試

激光多普勒測(cè)速法（LDV）是一種針對(duì)流體速度的非接觸測(cè)量手段，具有不干擾流體流動(dòng)、空間分辨率高、動(dòng)態(tài)響應(yīng)快和精度高等優(yōu)點(diǎn)，使LDV系統(tǒng)在流體力學(xué)領(lǐng)域得到了廣泛應(yīng)用。本實(shí)驗(yàn)主要選取中心立面中的若干點(diǎn)，對(duì)新型攪拌裝置的軸向速度進(jìn)行測(cè)量，并得出了相應(yīng)區(qū)域流場(chǎng)分布和流動(dòng)速度等實(shí)驗(yàn)結(jié)果。

3.1 實(shí)驗(yàn)過程

空氣壓縮機(jī)規(guī)格：公稱容積流量0.36 m/min，額定排氣壓力0.8 MPa，功率3 kW。將各個(gè)實(shí)驗(yàn)裝置連接，空氣壓縮機(jī)運(yùn)轉(zhuǎn)后將氣閥逐漸打開，直到攪拌器達(dá)到最高轉(zhuǎn)速的穩(wěn)定狀態(tài)。此時(shí)，空壓機(jī)閥門全開，攪拌器速度為60 r/min，計(jì)算得到氣孔噴氣速度為15 m/s。

試驗(yàn)主要研究槳葉以下區(qū)域的流動(dòng)狀況，槳葉以上由于氣泡原因產(chǎn)生劇烈的湍流而未進(jìn)行測(cè)試。實(shí)驗(yàn)取點(diǎn)：水平方向從槳葉底部所在平面壁面開始每隔5 mm 取1 個(gè)點(diǎn)，總共取23 個(gè)點(diǎn)；豎直方向從槳葉底部所在平面開始每隔10 mm 取1 個(gè)面；總共23×10=230 個(gè)點(diǎn)，每個(gè)點(diǎn)測(cè)試40 s，取平均值并得出數(shù)據(jù)。

3.2 實(shí)驗(yàn)數(shù)據(jù)處理

根據(jù)所得數(shù)據(jù)，本文選3 個(gè)面與相應(yīng)的數(shù)值模擬導(dǎo)出數(shù)據(jù)對(duì)比，每個(gè)面從23 個(gè)數(shù)據(jù)中篩選10 個(gè)數(shù)據(jù)作圖。圖8 為槳葉下部10 mm 平面的對(duì)比圖，圖9 為槳葉下部50 mm 平面的對(duì)比圖，圖10 為槳葉下部90 mm 面的對(duì)比圖。圖中橫坐標(biāo)為所測(cè)點(diǎn)距桶壁的距離，縱坐標(biāo)為速度，正負(fù)代表方向，向上為正。紅色曲線代表實(shí)驗(yàn)數(shù)據(jù)，黑色曲線代表數(shù)值模擬數(shù)據(jù)。

圖8 槳葉下部10 mm 液面的軸向速度

圖9 槳葉下部50 mm 液面的軸向速度

圖10 槳葉下部90 mm 液面的軸向速度

由圖8—圖10 可知：實(shí)驗(yàn)數(shù)據(jù)得到的速度明顯要小于模擬速度，靠近壁面處實(shí)驗(yàn)所得數(shù)據(jù)要比數(shù)值模擬得到的數(shù)據(jù)波動(dòng)大。圖8 中，在距桶壁9 cm 處的點(diǎn)實(shí)驗(yàn)所得數(shù)據(jù)對(duì)比模擬數(shù)據(jù)要先變?yōu)樨?fù)方向；圖9 中，在距桶壁為3 cm 處的數(shù)據(jù)比模擬數(shù)據(jù)要先變?yōu)檎较?；圖10 中，在距桶壁3 cm處的點(diǎn)同樣也是實(shí)驗(yàn)數(shù)據(jù)首先變正，并且在距桶壁6 cm 處首先達(dá)到波峰。究其原因，主要是：

（1）數(shù)值模擬是在理想狀態(tài)下得到的數(shù)據(jù)。試驗(yàn)中，攪拌器轉(zhuǎn)動(dòng)連接處會(huì)有一定量的氣體漏出，影響氣孔出氣量進(jìn)而使速度降低；試驗(yàn)中，經(jīng)過壓縮的氣體進(jìn)入水中后，由于溫度降低、氣泡收縮，也會(huì)使速度降低。

（2）由于光射入不同介質(zhì)以及有機(jī)玻璃曲率折射影響，使實(shí)驗(yàn)的測(cè)試點(diǎn)發(fā)生偏移，造成數(shù)據(jù)偏差。

（3）由于試驗(yàn)中的桶壁不光滑等因素，導(dǎo)致黏滯阻力較大，對(duì)靠近桶壁的流體流動(dòng)有一定的干擾。

4 結(jié) 論

本研究利用 Fluent 軟件對(duì)新型攪拌器的攪拌特性進(jìn)行數(shù)值模擬，獲得了桶內(nèi)的流體流動(dòng)和速度場(chǎng)分布，并與45°折葉槳進(jìn)行模擬對(duì)比，最后利用LDV 實(shí)驗(yàn)得出的數(shù)據(jù)對(duì)數(shù)值模擬結(jié)果進(jìn)行了驗(yàn)證，得出以下結(jié)論：

（1）氣體驅(qū)動(dòng)攪拌器與傳統(tǒng)45°折葉槳對(duì)比，槳葉以上部分速度分布比較均勻，在氣體的攪拌作用下低速區(qū)較少，并且靠近壁面的高速區(qū)域小，能更好地防止因高速流體集中在壁面區(qū)域而影響攪拌效率。

（2）與45°折葉攪拌槳相比，氣體驅(qū)動(dòng)攪拌器更注重2 個(gè)肋部的攪拌,攪拌范圍更大。

（3）由于一些不可避免因素的影響，導(dǎo)致實(shí)驗(yàn)數(shù)據(jù)和數(shù)值模擬數(shù)據(jù)有一定誤差，但是流場(chǎng)趨勢(shì)基本相同。

綜上所述，本文對(duì)新型氣體驅(qū)動(dòng)攪拌器進(jìn)行初步分析研究、模擬和試驗(yàn)，對(duì)攪拌器設(shè)備的更新?lián)Q代具有一定的參考價(jià)值，然而在實(shí)際應(yīng)用中可能還存在一些問題，將在今后的研究中不斷完善。