攪拌槽內固液兩相流的數值模擬及功率計算

董厚生 魏化中 舒安慶 劉 凱

（1.武漢工程大學機電工程學院 2.武漢市壓力容器壓力管道安全技術研究中心)

攪拌槽內固液兩相流的數值模擬及功率計算

董厚生*1魏化中1，2舒安慶1，2劉 凱2

（1.武漢工程大學機電工程學院 2.武漢市壓力容器壓力管道安全技術研究中心)

使用計算流體動力學的方法對攪拌槽中的流場進行模擬，得到攪拌槽中液體的流動狀況和體積分數分布。對流場分布規律、固體顆粒體積分數特點加以分析，進而利用模擬出的數據計算攪拌軸的功率，為攪拌器的設計提供參考。

攪拌槽 功率 顆粒 計算流體動力學 數值模擬

0 前言

攪拌設備廣泛應用于石油化工等行業中。攪拌設備其內部流動非常復雜，但對其研究目前還未形成完善的理論體系，設計者往往是依靠經驗和實驗數據來進行設計的。實踐證明，按照這種傳統方法設計出來的攪拌器大多不是處于最佳狀態，可靠性低。如何準確地描述攪拌槽中的流動狀況和混合過程，這是許多設計者所關心的問題。近年來，計算流體動力學 （CFD）技術發展迅速，運用數值模擬的方法獲得攪拌槽中的信息成為可能。本文利用計算流體動力學的方法模擬出攪拌槽中液體的流動狀況和固體顆粒的體積分數分布情況，并利用模擬出來的數據計算攪拌軸的功率，以期望對攪拌器的設計研究提供參考。

1 數值模擬

1.1 計算體系

計算采用的攪拌槽槽體為圓筒形，攪拌槽直徑D=15.2 m，液位高H=12.5 m，槽內均布四塊擋板。介質為磷酸和固體小顆粒，顆粒的密度ρ=1 800 kg/m3，顆粒的平均直徑為0.045 mm，物料中顆粒體積百分數為3.7%。攪拌槳采用PBT槳，槳葉直徑d=2.92 m,槳葉離底高度h=1.6 m。攪拌器轉速n=25 r/min。

1.2 網格劃分

計算中采用的網格是非結構化網格。將攪拌槳附近區域設為動區域 （如圖1所示），攪拌槽內動區域以外區域設為靜區域（如圖2所示）。其中，動區域網格數為188 134，靜區域網格數為665 408。

圖1 攪拌槳附近的動區域

圖2 攪拌槽內的靜區域

1.3 模擬方法

使用Fluent 6.3軟件進行計算，將液相設為連續相，固體顆粒設為分散相。計算中假設固液兩相間無質量交換，由于固相體積分數較低，故固液兩相間動量交換系數使用Wen-Yu模型。液相湍流模型采用κ-ε模型，使用多重參考系法，流動狀況設為定常流動，壓力速度耦合使用SIMPLE算法，差分格式采用一階迎風格式。

2 數值模擬結果與分析

2.1 流場分布

圖3為通過攪拌軸垂直截面的速度矢量圖。由圖3可以看出，液流從液面經葉輪流出，沖擊槽底，再沿槽壁流向液面，在攪拌槽中形成一個大的循環。正是利用這種沖擊力及液流的夾帶作用才使得固體顆粒在攪拌槽中混合或懸浮，以達到攪拌的目的。

為了更清楚地描述攪拌槽內流場的分布情況，選取與z軸平行并在zx平面上的4條直線，其距離z軸分別為 200 mm、1 460 mm、2 000 mm、6 000 mm。流場中速率分布如圖4所示。

圖3 通過攪拌軸垂直截面的速度矢量圖

圖4 流場中速率分布

從圖4可以看出，攪拌槽內的高速區主要集中在攪拌槳附近，而靠近液面位置的流速非常小，由此可以說明攪拌槳輸出的能量大多耗散在槳葉附近，而遠離槳葉的位置沒有得到很好的混合。

2.2 濃度分布

圖5所截取的平面為經過攪拌軸并與水平面垂直的平面，圖6顯示的4條直線分別為距離z軸200 mm、1 460 mm、2 000 mm、6 000 mm。從圖中固體顆粒的體積分數分布可以看出，液面的濃度分布最小，攪拌槽下方的濃度分布均勻，且固體顆粒體積分數最大。

圖5 垂直面內固體顆粒體積分數

圖6 不同距離的軸向濃度分布

經分析認為，固體顆粒在攪拌槽中存在兩種運動：顆粒本身的隨機運動和液流對其的夾帶運動。固體顆粒的運動主要取決于局部液流速度、顆粒的沉降速度和局部湍流強度。固體顆粒所處位置是其隨機運動和液流對其的夾帶運動綜合作用的結果。在液面流速較小，重力起主導作用，因而液面固體顆粒體積分數很低；在槽的下方位置，液流速度起主要作用，因此，固體顆粒在槽的下方位置能得到充分混合。

3 功率計算

對于旋轉的系統，功率為P=Mω。式中M為扭矩，ω為角速度。

攪拌槳在流場中受壓力和黏性力作用，因而作用在攪拌槳上的扭矩也由兩部分組成，M=M1+M2，其中M是攪拌槳所受的總的扭矩，M1是由壓力產生的扭矩，M2是由切應力產生的扭矩。由模擬得到攪拌槳的扭矩值，如表1所示。

表1 攪拌槳的扭矩

從表1可以看出，由壓力產生的扭矩M1占總扭矩的大部分，由切應力產生的扭矩M2占總扭矩的22.26%，故M2不可忽略，計算功率時應予以考慮。攪拌軸功率為

P=17.599 78×2π×25/60=46 kW

4 結論

利用Fluent軟件研究了攪拌槽中流場分布及濃度分布，并利用模擬得出的數值計算攪拌軸功率。

（1）用CFD的方法模擬出來的流場流向及規律跟以往研究者的研究結果相同，說明該模擬結果跟設計要求基本吻合。

（2）由固體顆粒體積分數的分布可以看出，固體顆粒在整個槽中基本成均勻分布，但在液面位置，固體顆粒體積分數很低。因此可以考慮通過增大槳葉的離底距離或者使用雙層槳葉來改善攪拌效果。

（3）利用模擬得到的數據計算出的功率跟利用《攪拌與混合設備設計選用手冊》計算出的功率相差不大，因此，模擬結果可為攪拌器的設計提供參考。

[1]陳志平,章序文,林興華，等．攪拌與混合設備設計選用手冊 [M]．北京：化學工業出版社，2004.

[2]侯栓弟．攪拌槽內三維流場的實驗研究與數值模擬．[D].北京：北京化工大學，1997.

[3]黃雄斌，施力田，等.固液攪拌槽中液相速度的分布[J].化工學報，2002（7）.

[4]馬青山，聶毅強.攪拌槽內三維流場的數值模擬 [J].化工學報，2005，54 （5）：612-618.

[5]張林進，葉旭初．攪拌器內湍流場的CFD模擬研究[J]．南京工業大學學報，2005，27(2)：59-63.

Numerical Simulation and Power Calculation of Solid-liquid Two-phase Flow in Stirred Tank

Dong Housheng Wei Huazhong Shu Anqing Liu Kai

Using computational fluid dynamics to simulate the flow field in stirred tank,the flow situation and volume fraction distribution of the liquid in stirred tank were obtained.The flow field distribution laws of liquid and the volume fraction characteristics of solid particles were analysed,and then using the date of simulation to calculate the power of agitator shaft,so as to provide reference for the design of mixer.

Stirred tank;Power;Particle;Computational fluid dynamics;Numerical simulation

TQ 051.7

*董厚生，男，1985年生，碩士研究生。武漢市，430205。

2011-09-05）