譚 歡，鄭維師，王正金

（1.中國石油集團川慶鉆探工程有限公司鉆采工程技術(shù)研究院，西安 710018；2.低滲透油氣田勘探開發(fā)國家工程實驗室，西安 710018）

0 引言

在石油化工等領(lǐng)域，經(jīng)常需要槳式攪拌器進行固液兩相攪拌，將固相顆粒分散到液體中，形成固液混合物，增強固液相間的傳質(zhì)傳熱[1]。槳式攪拌器作為一種混合機械，可用于黏性或有凝聚性的粉粒體的混合，廣泛應(yīng)用于農(nóng)藥、食品、化學(xué)品、生物、養(yǎng)殖業(yè)等固-固、固-液的混合[2-4]。由于流體中存在密度較大的固相顆粒，攪拌過程中容易出現(xiàn)沉積，從而影響攪拌效果。因此防止槳式攪拌器的固相顆粒沉積，使其在液相中達(dá)到懸浮狀態(tài)，是槳式攪拌器的攪拌特性研究的重要內(nèi)容[5]。前人在雙層攪拌器流動特性分析了做了許多研究，陳逸范等[6]通過實驗的方法,在改變多種雙層平斜槳式攪拌器的直徑和轉(zhuǎn)速情況下，得出攪拌器的攪拌效果以D/T=0.5雙層平槳為最佳；陳凱等[7]通過實驗的方法系統(tǒng)研究和比較了徑流槳和徑流槳組合、徑流槳和斜葉槳組合以及斜葉槳和斜葉槳組合3 類不同的雙層攪，在氣液分散攪拌過程中的優(yōu)劣。拌器組合；莫瑾等[8]使用Fluent對雙層攪拌器進行數(shù)值模擬，分析雙層攪拌器安裝間距在900 mm時攪拌器內(nèi)流動狀態(tài)最好；鄧華軍等[9]利用計算流體動力學(xué)（CFD）的方法,對攪拌器混合油漆混合液體的過程進行數(shù)值模擬,分析了這種自制的攪拌器層間距的改變對混合時間的影響；張慧敏等[10]利用Fluent 軟件模擬了單、雙層開啟渦輪式攪拌器在釜內(nèi)的流場，研究表明選用單層渦輪式攪拌器時流體流動范圍較窄，選用雙層渦輪式攪拌器時流體流動范圍較廣。

雖然對雙層攪拌器內(nèi)流動特性有了許多研究，但對雙層攪拌器固液混合效果方面的研究較少，因此本文針對雙層攪拌槳的攪拌器，利用Fluent軟件對攪拌器內(nèi)的固液混合流場特性進行數(shù)值模擬，研究了該攪拌器內(nèi)的速度流場情況，以及在不同轉(zhuǎn)速下固體顆粒的懸浮狀態(tài)。

1 建立模型

1.1 幾何模型

研究對象為弧形底圓柱式攪拌器，用三維建模軟件Solid-Works 建立的模型如圖1 所示，建模中心為攪拌器頂部中心。對于完全離底懸浮，只需一層葉輪，而對于均勻懸浮，則必須使用多層葉輪，因此攪拌器由兩層葉輪組成，同軸布置。模型參數(shù)：攪拌器直徑D=900 mm，高度H=900 mm，攪拌軸直徑D1=65 mm，上層葉輪長度L1=255 mm、寬度W1=40 mm，下層葉輪長度L2=200 mm、寬度W2=40 mm。

圖1 幾何模型

1.2 網(wǎng)格模型

將模型導(dǎo)入到Mesh 模塊中進行網(wǎng)格劃分，數(shù)值模擬涉及到旋轉(zhuǎn)計算因此將模型劃分為動區(qū)域和靜區(qū)域，其中動區(qū)域?qū)嚢铇ㄔ诶锩?，其余區(qū)域為靜區(qū)域，總體采用四面體非結(jié)構(gòu)化網(wǎng)格劃分，并且對內(nèi)層動區(qū)域進行網(wǎng)格加密，單元數(shù)為311 198，網(wǎng)格質(zhì)量為0.8，結(jié)果如圖2所示。

圖2 網(wǎng)格模型

2 數(shù)值模擬

2.1 控制方程

Mixture 模型可用于模擬各相之間存在速度差異的多相流動，該模型假定了在短空間尺度上局部的平衡[11]。由于固相顆粒與流體之間的相互作用力是未知的，且固相顆粒需要廣泛的分布于液相之中，對于這種情況可以Mixture 模型進行模擬。并且Mixture 將固液兩相都視為連續(xù)相，對于不涉及化學(xué)反應(yīng)的固液混合模擬，也需要考慮質(zhì)量和動量守恒方程。

連續(xù)性方程：

質(zhì)量守恒方程：

動量守恒方程：

2.2 邊界條件與求解方法

Fluent 處理攪拌問題時，主要提供了3 種解決辦法：多重參考系模型（MRF）、滑移網(wǎng)格模型SMM）和動網(wǎng)格模型[12]。本文模擬采用的MRF 模型，對不同的區(qū)域可以設(shè)置不同的轉(zhuǎn)速，而且每個運動區(qū)域網(wǎng)格中使用運動參考系方程求解流動場。本文中葉輪和攪拌器邊界均設(shè)置為壁面邊界條件，葉輪半徑范圍內(nèi)的流體區(qū)域定義為旋轉(zhuǎn)參考系，設(shè)置該區(qū)域的轉(zhuǎn)速，旋轉(zhuǎn)方向定位Y軸正方向。采用壓力基穩(wěn)態(tài)求解，壓力-速度耦合方法選用SIMPLE 算法，空間離散化梯度設(shè)置為Green-Gauss Node Nased，采用一階迎風(fēng)差分格式。設(shè)置底部300 mm 深度區(qū)域的固體顆粒作為初始化條件，固相初始分布如圖3 所示。

圖3 初始狀態(tài)固液相分布

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

本文中雙層攪拌槳的目的是讓固體顆粒均勻的懸浮在液態(tài)溶液中，雙層攪拌槳旋轉(zhuǎn)過程中將自身的速度和動能不斷地傳遞到流體域中來實現(xiàn)固體顆粒的均勻懸浮，因此需要分析轉(zhuǎn)速對固液均勻懸浮的影響，以及分析轉(zhuǎn)速一定時不同位置的速度場和顆粒的體積分布情況，進而為相關(guān)攪拌器的設(shè)計研究提供一定的科學(xué)依據(jù)。

3.1 攪拌器軸向速度場分析

在轉(zhuǎn)速為300 r/min 時攪拌器軸截面的速度云圖如圖4 所示。從圖中可以看出，攪拌器的總體流場特征呈現(xiàn)為對稱分布，在攪拌槳端部處有明顯的徑向流動且速度較高，故攪拌流場可以認(rèn)為是一種典型的徑向流分布。在攪拌器的頂部和底部位置，流體的流動基本沒有影響，要想固液混合均勻，因此要增大轉(zhuǎn)速以此改變流場運動。

圖4 軸向截面速度云圖

3.2 攪拌槳處截面速度分析

本文中攪拌器為雙層攪拌槳，其所在截面的速度分布云圖如圖5 所示。從圖中可以看出，上層攪拌器周圍的流動速度比下層攪拌器周圍的流動速度較大，并且在高速旋轉(zhuǎn)下都存在打旋的現(xiàn)象。

圖5 攪拌器截面速度云圖

3.3 轉(zhuǎn)速對固相顆粒分布的影響

分別在轉(zhuǎn)速為50、100、150、200、250、300、350 r/min下分析攪拌器內(nèi)固液兩項分布情況，結(jié)果如圖6 所示。從圖中可以看出，在轉(zhuǎn)速較低的時候固相幾乎都沉積在底部沒有發(fā)生流動，當(dāng)轉(zhuǎn)速增大到150 r/min 時靠近攪拌器避免的固體顆粒向著中心發(fā)生流動，轉(zhuǎn)速繼續(xù)增大時流動現(xiàn)象更加明顯，當(dāng)轉(zhuǎn)速到達(dá)350 r/min 時可以看出底部沒有明顯的顆粒沉積了。

圖6 不同轉(zhuǎn)速下固相分布云圖

4 結(jié)束語

（1）通過Fluent 對雙層攪拌槳的攪拌器進行了流場分析，結(jié)果表明攪拌器為徑向流攪拌器，并且在高速選裝下，攪拌器流體存在明顯的打旋現(xiàn)象。

（2）分析了轉(zhuǎn)速在50～350 r/min 之間，攪拌器內(nèi)固相顆粒的分布狀態(tài)，結(jié)果表明轉(zhuǎn)速較低時攪拌器內(nèi)的固體幾乎都沉積在底部，隨著轉(zhuǎn)速的增加固體顆粒才開始向攪拌器上部流動，當(dāng)轉(zhuǎn)速達(dá)到350 r/min 時，攪拌器底部已無明顯的顆粒沉積。