磷酸反應槽內(nèi)臨界離底懸浮轉(zhuǎn)速的CFD模擬

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

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

磷酸反應槽內(nèi)臨界離底懸浮轉(zhuǎn)速的CFD模擬

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

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

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

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

0 引 言

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

1 計算體系與模擬方法

1.1 計算體系

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

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

表1 攪拌槳參數(shù)表Table.1 parameters of the PBT airfoil impeller

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

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

圖2 網(wǎng)格劃分示意圖Fig.2 Meshing for the stirred tank

1.2 模擬方法

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

2 結(jié)果與討論

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

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

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

2.1 均勻懸浮狀態(tài)

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

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

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

2.2 不同轉(zhuǎn)速下的顆粒濃度分布

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

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

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

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

2.3 臨界懸浮轉(zhuǎn)速

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

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

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

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

3 結(jié) 語

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

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

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

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

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

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

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

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

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

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

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

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

［12］ 包麗云.固液攪槽中液相速度測量技術開發(fā)及攪拌裝置在水處理過程中的應用研究［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），男，湖北安陸人，碩士研究生.研究方向：過程裝備模擬與仿真.

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