磷酸反應槽內臨界離底懸浮轉速的CFD模擬

（1.武漢工程大學機電工程學院，湖北武漢430205；

2.武漢市壓力容器壓力管道安全技術研究中心，湖北武漢430074）

磷酸反應槽內臨界離底懸浮轉速的CFD模擬

（1.武漢工程大學機電工程學院，湖北武漢430205；

2.武漢市壓力容器壓力管道安全技術研究中心，湖北武漢430074）

使用計算流體軟件Fluent對磷酸反應槽內顆粒的臨界離底懸浮轉速進行了數值模擬.攪拌槽直徑.T＝0.5 m，四塊擋板均布，攪拌槳采用45°斜葉槳.兩相物系為磷石膏—硫酸，固體體積百分比濃度φ＝4.65%.使用濃度判據得到顆粒離底臨界懸浮轉速NJS，模擬計算結果的誤差在工業允許的范圍內.模擬得出攪拌槽中液體的流動狀況和固體體積分數的分布；同時對6個不同攪拌轉速下的固體顆粒懸浮狀況進行比較，得出均勻懸浮臨界轉速.

計算流體力學（CFD）；臨界離底懸浮轉速；攪拌槽

0 引 言

如何以最小的能耗獲得所需的攪拌效果一直以來是固液懸浮技術研究的重點.因此，找出固液體系的臨界離底轉速顯得尤為重要［1］.當顆粒全部處于運動時，且顆粒在槽底停留時間不超過1～2 s，即認為達到了完全離底懸浮，而達到此狀態所需的最低攪拌轉速稱為臨界轉速NJS.在固液懸浮操作中，對完全離底懸浮的研究較多，也發表了不少有關臨界轉速的關聯式.其中，Zwietering提出的完全離底懸浮臨界轉速關聯式應用的最為廣泛.目前得出的NJS值大多是在實驗室尺度的槽中通過冷模實驗估算獲得的，再放大到工業反應器中.而對于工業上應用的高溫、高壓等復雜條件下的大尺度反應器，是無法通過實驗獲得所需的詳盡信息的.因而需要尋找一種新的途徑來對工業中應用的攪拌槽內的臨界懸浮轉速NJS進行預測.計算流體力學（CFD）技術的發展讓這種預測成為可能.關于完全離底懸浮臨界轉速NJS方面的報道比較少，本文的目的是預測攪拌槽內固液懸浮的臨界懸浮轉速，將模擬結果和文獻數據進行對比，驗證Fluent軟件對攪拌槽內臨界離底懸浮轉速進行預測的合理性.

1 計算體系與模擬方法

1.1 計算體系

1.1.1 攪拌槽結構 模擬采用的攪拌槽槽體為圓柱形，四塊擋板均布，攪拌槽直徑T＝0.5 m，液面高H＝0.46 m，擋板寬W＝，擋板離底距離為0.01 m.漿型采用45度斜葉槳（攪拌槳示意圖如圖1所示，參數如表1所示），槳葉離底距離C＝.

圖1 PBT攪拌槳示意圖Fig.1 Configuration of PBT airfoil impeller

表1 攪拌槳參數表Table.1 parameters of the PBT airfoil impeller

1.1.2 計算物系 介質為磷石膏和硫酸，固體相的體積分數為4.65%，顆粒的平均直接為72 μm，密度為2 320 kg／m3.流體操作密度為1 500kg／m3，操作粘度為5 MPa·s.

1.1.3 網格劃分 根據流動的對稱性，選取槽體的四分之一為計算域.由于模擬所采用的攪拌槳的結構不規則，網格劃分采用非結構化的網格，共1 014 262個網格節點，如圖2所示.

圖2 網格劃分示意圖Fig.2 Meshing for the stirred tank

1.2 模擬方法

計算使用的軟件是Fluent 6.3，利用Euleian模型模擬該固—液兩相流［2－3］.Euleian模型分別解每一相動量方程，并通過連續相對分散相的曳力與動量交換項來模擬兩相間相互作用.并采用標準的壁面函數法來處理近壁區域流動的計算.而液相湍流模型采用標準k－ε方程模型.計算中假設固—液兩相間無質量傳遞，由于固相濃度較低，所以相應的阻力系數采用wen－yu模型［4－5］.流動場的計算使用多重參考系法，并設定流動性質為定常流動.相應的壓力速度耦合選擇Phase Coupled SIMPLE法，并采用一階迎風的差分格式［6－7］.

2 結果與討論

均勻懸浮臨界轉速的確定最常用的方法是通過測釜內各點的固相濃度，根據釜內固相濃度分布（固體顆粒的體積分數分布情況反映了固體顆粒的的濃度分布狀況）的均勻度來判斷.

衡量攪拌釜內固體顆粒濃度分布均勻性的判據很多.廣泛采用的是濃度分布的標準偏差σ.

本文研究選取數值模擬的攪拌轉速為20 r／min、30 r／min、40 r／min、50 r／min、60 r／min，70 r／min，其對應的雷諾數分別為0.83×105，1.25× 105，1.67×105，2.08×105，2.5×105，2.92×105.

2.1 均勻懸浮狀態

均勻懸浮狀態下的速度、濃度分布如圖3所示.

圖3 固液懸浮流場分布Fig.3 Flow field for dilute suspension

圖3為60 r／min轉速下的Fluent數值模擬后的混合流的液相速度矢量圖與固相體積分數分布圖.圖3（a）為經過攪拌軸垂直面內的流場，其中矢量長度代表速度大小，攪拌槳區域的速度值較大，而自由液面區的速度值較小.圖3（b）為固相體積分數分布圖，表示了槽內的固體分布情況.

2.2 不同轉速下的顆粒濃度分布

圖4列舉了20～70 r／min轉速下的濃度分布圖.

從圖4中固體顆粒的體積分數分布可以看出，20 r／min轉速時，液面處的濃度大小跟其他轉速下的濃度大小不在一個數量級上.這是由于在液面處，速度過小時，重力起主導作用造成的.

圖4 不同轉速下垂直截面固相體積分數分布Fig.4 Vertical direction volume fraction contours of the solids with different rotation speeds

另外，隨著轉速的增加，槽內顆粒濃度更加趨于均勻狀態，在攪拌槽的自由液面區顆粒濃度也隨之增大.隨著轉速的不斷增加，槽內顆粒濃度分布的標準偏差趨于定值時，即認為達到了均勻懸浮.此時所需的最小轉速即為均勻懸浮臨界轉速［9］.

2.3 臨界懸浮轉速

歸納以上幾種轉速下槽內的顆粒濃度分布狀況，繪制沿圖5中測試線L中的幾個測試點隨著轉速的增加其顆粒濃度的標準偏差，結果如圖6所示.

圖5 測試點示意圖Fig.5 The distribution of testpoints

圖6 顆粒濃度偏差圖Fig.6 the concentration variance of solids

從圖6中可以得出轉速在20～70 r／min內，槽內上層顆粒濃度偏差均較大，往下則偏差越來越小，當轉速達到40 r／min時，槽內的濃度分布標準偏差減小并趨于穩定，即認為達到了均勻懸浮狀態，由此可推斷得到本體系的的顆粒懸浮轉速NJS＝40 r／min，根據文獻經驗公式計算得出NJS為38～45 r／min，經計算兩者誤差為5%～12.5%［10－12］.

3 結 語

以上是使用計算流體力學的方法對攪拌槽內的低濃度固液懸浮狀況進行了模擬，通過模擬計算，得出不同攪拌轉速下的固體顆粒濃度分布圖.利用濃度判據的方法推斷出臨界離底懸浮轉速NJS.模擬結果與文獻結果比較吻合，不過，此結果仍需用實驗來驗證，初步驗證了CFD技術對NJS預測的可能性，模擬方法與結果可對設備的設計和放大起到一定的指導作用.

［1］ 陳志平.攪拌與混合設備設計選用手冊［M］.北京：化學工業出版社，2004.

［2］ 馬青山，包云雨，聶毅駿.攪拌槽內的三維流場的數值模擬［J］.化工學報，2003，54（5）：612－618.

［3］ 朱榮生，李維斌，黃道見，等.攪拌器攪拌流場的三維數值模擬［J］.農機化研究，2003，10（4）：75－77.

［4］ 王振松，黃雄斌，李良超，等.固—液攪拌槽內槽底流場的CFD模擬［J］.北京化工大學學報，2005，32（4）：5－9.

［5］ 王軍，吳鐸，楊志春，等.熱聲諧振管壓比影響因素的數值模擬［J］.武漢工程大學學報，2010，32（2）：80－83.

［6］ 孫亞忠，陳作炳，董新菅.立磨內部流場的數值模擬［J］.武漢工程大學學報，2011，33（4）：89－93.

［7］ 常玉鋒，陳鋒.排煙方式對火場人員疏散影響的數值模擬［J］.武漢工程大學學報，2010，32（9）：46－49.

［8］ 黃雄斌，閆憲斌，施力田，等.固—液攪拌槽內液相速度的分布［J］.化工學報，2002，53（7）：717－722.

［9］ 李良超，黃雄斌.固液攪拌槽內近壁區液相速度研究［J］.北京化工大學學報，2005，32（1）：33－38.

［10］ 金中麗，黃雄武，賈志剛.固液攪拌槽內顆粒離底懸浮轉速的CFD模擬［J］.北京化工大學學報，2003， 30（60）：18－22.

［11］ 祝鈴鈺.固液懸浮攪拌器中的局部速度測量［D］.北京：北京化工大學，1996.

［12］ 包麗云.固液攪槽中液相速度測量技術開發及攪拌裝置在水處理過程中的應用研究［D］.北京：北京化工大學，1995.

CFD simulation of critical suspension lmpeller speed in a phosphoric acid reaction tank

DONG Hou－sheng1，WEI Hua－zhong1，2，SHU An－qing1，2，LIU Kai2
（1.School of Mechanical and Electrical Engineering，Wuhan Institute of Technology，Wuhan 430205，China；
2.Wuhan research center of pressure vessel and pipeline safety engineering，Wuhan 430074，China）

Based on the function provided by Fluent，the critical suspension impeller speed，Njs，in a phosphoric acid reaction tank was simulated.The investigation was carried out in a fully baffled，flat bottom，cylindrical vessel with 500 mm diameter，equipped with PBT impeller.The phosphogypsum of 4.65%was chosen as the dispersed phase.The original criterion was used to determine Njsfrom the numerical results，and the error of simulation is acceptable to the industry standard.Solid concentration distribution and the velocity distribution of solid and liquid were predicted.The performance of solid suspension under the different agitation speeds was also studied，the critical suspension impeller speed was obtained.

computational fluid dynamics（CFD）；critical suspension speed of impeller；stirred tank

陳小平

TQ027.1

A

10.3969／j.issn.1674－2869.2011.10.020

16742869（2011）10008604

20110721

董厚生（1985），男，湖北安陸人，碩士研究生.研究方向：過程裝備模擬與仿真.

指導老師：魏化中，男，副教授，碩士研究生導師.研究方向：壓力容器壓力管道、流體機械和攪拌設備的開發、新型高效過程裝備及仿真技術.＊通信聯系人