李介普， 趙 斐

（北京航天試驗(yàn)技術(shù)研究所， 北京100047）

0 引言

在石油行業(yè)中，石油在開采、煉制等過程中會(huì)有大量的油泥產(chǎn)生，油泥[1]的組成充分很難明確定義，是一種多相的混合物，主要的成分為原油、污水、油砂等。油泥既是生產(chǎn)中的廢物，又是可以循環(huán)利用的二次能源，需要一套完整的油泥清洗裝置，其中油泥清洗器的作用是將油水固三相充分混合，通過攪拌器的攪拌，使流體得到動(dòng)勢(shì)能，在攪拌槽形成合適預(yù)期效果的流場(chǎng)[2-4]。本文采用CFD 方法，研究了得到油泥清洗器攪拌槽內(nèi)的局部流動(dòng)狀態(tài)和混合狀況，對(duì)于油泥清洗器中攪拌器的設(shè)計(jì)提供了一定的參考。

1 數(shù)學(xué)模型

1.1 理論模型

本論文采用基于歐拉-歐拉模型下的歐拉兩相流模型，根據(jù)實(shí)際工藝參數(shù)，罐內(nèi)實(shí)際為油水固三相，由于含油量較少，將水和泥沙作為液-固兩相。 兩相初始狀態(tài)簡(jiǎn)化為固體泥沙沉于清洗罐底，泥沙沉底高度為1200mm。兩相動(dòng)量交換系數(shù)選用Gidaspow 模型，湍流模型采用標(biāo)準(zhǔn)k-ε模型[5]，模型的湍動(dòng)能k 和湍動(dòng)耗散率ε 方程如下式：

式中：Gk—由速度梯度變化而產(chǎn)生的湍流動(dòng)能；Gb—受浮力的作用產(chǎn)生的湍流動(dòng)能；YM—可壓縮流體中的脈動(dòng)膨脹對(duì)耗散率的作用；μl—湍流粘性系數(shù)。在Fluent 里，作為默認(rèn)值，C1ε=1.44，C2ε=1.92，C3ε=0.09；湍流普朗克數(shù)分別為σε=1.3，σk=1.0。

液固兩相流的模擬求解方法采用基于非穩(wěn)態(tài)的隱性非耦合算法求解，壓力與速度之間的耦合采用Phase Coupled SIMPLE 算法，動(dòng)量方程、湍動(dòng)能方程等方程全部按一階迎風(fēng)的格式進(jìn)行離散[6]。

1.2 計(jì)算模型

設(shè)備采用平底、平蓋攪拌罐，為常壓容器，體積75m3，容器直徑4000mm，容器高度6000mm，攪拌器為三段槳葉式攪拌器，兩瓣槳葉，槳葉長(zhǎng)度3000mm，槳葉高度和槳葉寬度的尺寸分別為0.1 和0.01 倍槳葉長(zhǎng)度的尺寸， 最下一層槳葉離罐底200mm，每一層相距1400mm，攪拌轉(zhuǎn)速18r/min。

本文采用滑移網(wǎng)格模型進(jìn)行網(wǎng)格劃分?；凭W(wǎng)格模型在計(jì)算時(shí)，將計(jì)算流域分為兩大部分：一個(gè)區(qū)域包含運(yùn)動(dòng)的槳葉和下半部分?jǐn)嚢栎S，即轉(zhuǎn)子區(qū)；另一區(qū)域包靜止的容器和上半部分?jǐn)嚢栎S，即定子區(qū)；兩個(gè)區(qū)域之間有網(wǎng)格分界面。 模型采用非結(jié)構(gòu)化網(wǎng)格進(jìn)行劃分，采用TGrid 類型的三面體網(wǎng)格，同時(shí)為了保證網(wǎng)格質(zhì)量，在動(dòng)區(qū)域的槳葉的網(wǎng)格先進(jìn)行面網(wǎng)格劃分，后采用結(jié)構(gòu)化的四邊形網(wǎng)格，攪拌罐的網(wǎng)格劃分圖如圖1 所示。

圖1 攪拌罐的網(wǎng)格劃分圖Fig.1 Geometric diagram of stirred tank

1.3 物料參數(shù)

工藝參數(shù)為含固10%進(jìn)料，油水固三相進(jìn)料，物料停留時(shí)間3 小時(shí)，含油量較少可忽略。 各物料的質(zhì)量流量分別為：水19662kg/h，固體（泥沙）2280kg/h。 固體的相體積分?jǐn)?shù)為90%。取泥沙顆粒粒徑約為111μm，密度2500kg/m3，粘度0.001003Pa·s。

2 模擬結(jié)果分析

2.1 攪拌器垂直截面的速度云圖和速度矢量圖

圖2 所示為三層槳葉式攪拌器垂直截面的速度云圖和速度矢量圖。由速度云圖可知，攪拌器葉片周圍的流體速度最大，槳葉端部的速度比槳葉中心靠近攪拌軸的地方大。兩個(gè)槳葉間的區(qū)域速度較小，最上層的槳葉上部有集中的高亮區(qū)域，此處的吸入速度較大，葉片噴射的液流進(jìn)入低速的流體中，卷吸周圍的流體，隨后沿著徑向擴(kuò)散。由速度矢量圖可知，在槳葉的上方和下方均形成了一個(gè)渦旋流動(dòng)，三層槳葉附近均出現(xiàn)雙渦旋環(huán)流場(chǎng)結(jié)構(gòu)，槳葉與槳葉間有相互作用， 流體在葉輪端部的出口處有很強(qiáng)烈的徑向運(yùn)動(dòng)，徑向速度較大。

圖2 垂直截面速度云圖和速度矢量圖Fig.2 Velocity contours and velocity vectors of vertical section

在葉輪區(qū)域，葉輪旋轉(zhuǎn)直接推動(dòng)流體流動(dòng)，形成主流，流動(dòng)特征具有射流特征，射流出去的流體撞得攪拌槽壁，分成兩個(gè)分支，分別向上、下流動(dòng)，在遠(yuǎn)離葉輪的區(qū)域，射流流體卷帶周圍的流體形成次流。 在葉輪區(qū)域的三個(gè)速度分量中，主要是徑向速度和軸向速度為主。槳葉下方靠近攪拌軸中心的區(qū)域，流體出現(xiàn)反向的流動(dòng)，這個(gè)區(qū)域叫誘導(dǎo)錐形區(qū)，該區(qū)域附近的流體速度很小，容易出現(xiàn)攪拌混合不均勻現(xiàn)象。

2.2 攪拌器的軸向速度云圖和徑向速度云圖

圖3 所示為軸向速度云圖和底層槳葉的徑向速度云圖。由軸向速度云圖可知，槳葉端部上下兩側(cè)的軸向速度較大，上下兩側(cè)的軸向速度方向相反。由于頂層槳葉上方無槳葉區(qū)域的環(huán)流渦旋現(xiàn)象較明顯，最頂層槳葉上端區(qū)域的軸向速度最大，液固兩相混合較充分。由底層槳葉的徑向速度云圖可知，沿著攪拌器的旋轉(zhuǎn)方向，葉輪的端部的徑向速度最大，槳葉兩側(cè)的流體徑向速度分布較為均勻，整體形成一個(gè)大環(huán)流， 攪拌器的周圍比攪拌槽邊壁處的徑向速度大，有利于液固兩相流混合。

圖3 垂直截面的軸向速度云圖和底層槳葉徑向速度云圖Fig.3 Axial velocity contours of vertical section and radial velocity of bottom blade

3 結(jié)論

油泥清洗器攪拌漿葉周圍的流體速度最大， 葉片附近的流動(dòng)具有射流特征，槳葉上下方出現(xiàn)雙渦旋環(huán)流場(chǎng)結(jié)構(gòu)，有利于兩相流混合。 在槳葉下方靠近攪拌軸中心的誘導(dǎo)錐形區(qū)，流體速度很小，容易出現(xiàn)攪拌混合不均勻現(xiàn)象。

受環(huán)流渦旋的影響， 槳葉端部上下兩側(cè)的軸向速度較大，葉輪的端部的徑向速度最大，槳葉兩側(cè)的流體徑向速度分布較為均勻，有利于液固兩相流混合。