魏化中，陶保林，舒安慶，金 戀

(1.武漢工程大學(xué)機(jī)電工程學(xué)院，湖北 武漢 430205； 2.武漢市壓力容器壓力管道安全技術(shù)研究中心，湖北 武漢 430074)

0 引 言

磷酸是重要的基本化工原料之一，又是磷化工中重要的中間產(chǎn)品.磷酸采用濕法生產(chǎn)的占多數(shù).用硫酸、硝酸或鹽酸分解磷礦制得的磷酸統(tǒng)稱為濕法磷酸，而用硫酸分解磷礦制得磷酸的方法是濕法磷酸生產(chǎn)中最基本的方法[1-3].濕法磷酸生產(chǎn)技術(shù)的關(guān)鍵在于獲得粗大、均勻、穩(wěn)定易于過濾的硫酸鈣水合結(jié)晶，以利于磷石膏與磷酸的分離[4]；有資料表明攪拌器的結(jié)構(gòu)型式和混合特性很大程度上決定了磷酸陳化槽內(nèi)固液物質(zhì)的混合效果[5]，深入研究磷酸陳化槽攪拌器的結(jié)構(gòu)對(duì)進(jìn)一步改善磷酸陳化槽固-液混合效果、提高后續(xù)的過濾效果及減低功率消耗，避免磷石膏在槽底沉積結(jié)垢都具有重要的現(xiàn)實(shí)意義.

然而，實(shí)際情況是攪拌槽內(nèi)物質(zhì)的流動(dòng)通常都是非常復(fù)雜的.對(duì)攪拌裝置的設(shè)計(jì)和放大方法還停留在傳統(tǒng)的經(jīng)驗(yàn)設(shè)計(jì)階段，這種方法的短板是不能觸及攪拌槽內(nèi)復(fù)雜流體流動(dòng)過程的本質(zhì)，所以這種方法始終不能圓滿解決攪拌裝置的設(shè)計(jì)和放大問題.另外，用于測(cè)量攪拌流場(chǎng)的實(shí)驗(yàn)儀器價(jià)格一般都很昂貴，并且流場(chǎng)測(cè)量的工作是一項(xiàng)非常費(fèi)時(shí)的工作；有時(shí)還會(huì)因?yàn)閿嚢璨蹆?nèi)介質(zhì)特性的原因，從而造成對(duì)該物質(zhì)的攪拌混合實(shí)驗(yàn)無法進(jìn)行.

近年來，計(jì)算流體力學(xué)(CFD)技術(shù)發(fā)展迅速，運(yùn)用CFD的方法獲取攪拌槽中的信息成為可能.本文以某公司的磷酸陳化槽裝置為原型，利用商業(yè)CFD軟件ANSYS13.0里的FLUENT為工具，對(duì)磷酸陳化槽內(nèi)固體顆粒與高粘度溶液的攪拌混合過程進(jìn)行了數(shù)值模擬、解析，以期為磷酸陳化槽攪拌器的后續(xù)設(shè)計(jì)、優(yōu)化提供參考.

1 模型與數(shù)值方法

1.1 磷酸陳化槽物理模型

模擬所采用的磷酸陳化槽為圓筒形，其結(jié)構(gòu)如圖1所示.槽體直徑為T=1 000 mm，液面高度為H=1 000 mm，沿圓周方向均布四塊擋板，擋板寬度為100 mm.所采用的液固物系為磷酸溶液和磷石膏顆粒，固體相的體積分?jǐn)?shù)為5.0%，固體顆粒形狀設(shè)為圓形；由于介質(zhì)中磷石膏顆粒直徑大于0.045 mm的占70%，故取固相顆粒直徑為0.045 mm，密度為2 320 kg/m3.磷酸溶液的密度為1 285 kg/m3，動(dòng)力粘度為0.015 Pa·s.

圖1 磷酸陳化槽幾何模型Fig.1 Geometric model of phosphate agin tank

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

1.2 攪拌轉(zhuǎn)速的確定

在攪拌設(shè)備的工業(yè)應(yīng)用中，固-液兩相混合的使用目的主要有以下兩個(gè)：a.使固體顆粒物質(zhì)在攪拌槽中均勻懸浮起來，簡(jiǎn)稱均勻懸浮.b.使固體顆粒物質(zhì)在攪拌槽中完全離底懸浮起來，簡(jiǎn)稱為完全離底懸浮.

完全使固體顆粒物質(zhì)達(dá)到離底懸浮狀態(tài)時(shí)的最低攪拌轉(zhuǎn)速叫做臨界轉(zhuǎn)速Njs.自從在50年前引入臨界轉(zhuǎn)速Njs概念后，它一直都是人們?cè)谔幚砉?液兩相流問題時(shí)所采用的最主要設(shè)計(jì)參數(shù)之一.在攪拌裝置設(shè)計(jì)過程中，對(duì)于固-液兩相流體系，有不少人發(fā)表了關(guān)于臨界轉(zhuǎn)速的關(guān)聯(lián)式的文章.在實(shí)際的生產(chǎn)過程中發(fā)現(xiàn)Zwietering的關(guān)聯(lián)式最為精確.

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

Zwietering的完全離底關(guān)聯(lián)式可寫為以下形式：

其中μ1是液體粘度(Pa·s)；ρ1是液體密度(kg/m3)；Δρ是固體和液體的密度差；dp是懸浮粒子直徑(m)；X是固體顆粒濃度；s是無因次幾何因子，從手冊(cè)中查到PBT的s值約為10.可以計(jì)算出磷石膏顆粒在陳化槽內(nèi)臨界轉(zhuǎn)速為1.6 r/s，即90 r/min，工程上通常要乘以安全系數(shù)1.5，結(jié)果為135 r/min，取為150 r/min.

1.3 固-液懸浮判據(jù)

1.4 模擬策略

1.4.1 網(wǎng)格劃分 采用fluent MixSim對(duì)模型進(jìn)行網(wǎng)格生成，生成的網(wǎng)格是非結(jié)構(gòu)四面體網(wǎng)格，然后通過FLUENT軟件進(jìn)行處理最后轉(zhuǎn)化為六面體網(wǎng)格，這樣處理就能使網(wǎng)格的數(shù)量大大減少，從而提高了軟件的計(jì)算速度.由于攪拌槳槳葉結(jié)構(gòu)比較復(fù)雜，為了提高近槳區(qū)域的流場(chǎng)的計(jì)算精度，所以對(duì)此區(qū)域內(nèi)的網(wǎng)格進(jìn)行加密處理，最后生成的總網(wǎng)格數(shù)量為170 302.

1.4.2 模擬方法 使用ANSYS 13.0里的Fluent軟件進(jìn)行計(jì)算，使用有限體積法求解離散方程，選用Eulerian-Eulerian模型對(duì)陳化槽的固液體系進(jìn)行數(shù)值模擬，液相湍流模型選用 standardk-ε湍流模型，近壁區(qū)域流動(dòng)的處理方法采用標(biāo)準(zhǔn)壁面函數(shù)法，由于計(jì)算體系中斯托克斯數(shù)要比1小很多，因此多相k-ε湍流模型選用Dispersed湍流模型[7-8].顆粒最大堆積體積取0.6，相間阻力系數(shù)使用Gidaspow模型，采用多重參考系法(MRF)對(duì)陳化槽內(nèi)固-液兩相流流場(chǎng)進(jìn)行數(shù)值模擬，槳葉及其附近的區(qū)域定義為內(nèi)計(jì)算域，其它區(qū)域定義為外計(jì)算域，內(nèi)計(jì)算域采用旋轉(zhuǎn)坐標(biāo)系，外計(jì)算域采用靜止坐標(biāo)系.控制方式的傳送項(xiàng)采用SIMPLE算法，差分格式采用一階迎風(fēng)格式，時(shí)間步長(zhǎng)設(shè)為0.01 s，流體流動(dòng)為定常流動(dòng)，所有項(xiàng)的殘差收斂標(biāo)準(zhǔn)均采用10[9-10].

2 模擬結(jié)果及分析

2.1 流場(chǎng)分布

圖3為磷酸陳化槽CFD模擬所得出的攪拌槽x=0軸截面的速度矢量圖，其所顯示的為陳化槽的流場(chǎng)，其中矢量的大小代表液相速度的大小.從圖3中可以得出以下幾點(diǎn)結(jié)論，①葉輪尖端附近的液流速度最大，自由液面處的速度最小.②葉輪兩側(cè)各有一個(gè)大的漩渦，說明在葉輪的高速轉(zhuǎn)動(dòng)下，磷酸陳化槽內(nèi)形成了一個(gè)大的循環(huán)區(qū)域.③葉輪附近的液體在電機(jī)所提供的外力作用下經(jīng)過葉輪的作用沖向槽底，最后沿著陳化槽的壁面流向液面，在此循環(huán)過程中磷石膏顆粒在葉輪所提供的沖擊力和溶液的速度帶動(dòng)在攪拌槽中流動(dòng)，以致混合或懸浮，最終達(dá)到攪拌混合的目的.④在葉輪的下方形成了一個(gè)倒錐體區(qū)域，在這個(gè)區(qū)域中液相的流動(dòng)方向與周圍其他區(qū)域內(nèi)的液流方向相反，為自下向上，同時(shí)這個(gè)區(qū)域的流速與周圍區(qū)域的流速相比比較小，因此在攪拌轉(zhuǎn)速不夠的情況下這個(gè)區(qū)域相對(duì)于其它區(qū)域來說更容易發(fā)生磷石膏固體顆粒堆積的現(xiàn)象.同時(shí)還可以看出在槽底靠近壁面區(qū)域的液流速度與周圍速度相比相差兩個(gè)數(shù)量級(jí)幾乎為0，此區(qū)域可以說是死區(qū)，固體顆粒在此區(qū)域極易發(fā)生堆積現(xiàn)象.

圖3 軸截面速度矢量分布Fig.3 Velocity vector distribution on axis cross-section

為了更詳細(xì)地來描述磷酸陳化槽的攪拌流場(chǎng)，現(xiàn)截取磷酸陳化槽橫截面的速度矢量圖，圖4所示平面為距離陳化槽底部垂直方向上方z=10 mm處的水平截面速度矢量圖.

圖4 水平截面速度矢量分布Fig.4 Velocity vector distribution on horizontal cross-section

從圖4中可以看出，水平平面的速度矢量方向是以攪拌軸為中心向四周發(fā)散的，這與圖3中所觀察到的情形是一致的，即液流匯聚成體積流后，在剪切力的作用下沿葉輪推進(jìn)方向運(yùn)動(dòng)，最終使槽內(nèi)介質(zhì)達(dá)到攪拌混合的目的.

2.2 濃度分布

圖5是磷酸陳化槽內(nèi)x=0軸截面上磷石膏固體顆粒體積分?jǐn)?shù)分布的相圖，圖6是磷酸陳化槽內(nèi)z=10mm水平截面上的固體顆粒體積分?jǐn)?shù)分布的相圖.從圖5和6圖中可以看出，在擋板附近以及槽底的角落位置和中心位置的磷石膏固體顆粒體積分?jǐn)?shù)較大，整個(gè)大循環(huán)漩渦區(qū)域以及攪拌槳尖端附近的磷石膏固體顆粒體積分?jǐn)?shù)相對(duì)來說較低，液面附近的磷石膏固體顆粒體積分?jǐn)?shù)幾乎為0.

圖5 軸截面固體顆粒體積分?jǐn)?shù)分布Fig.5 Axis cross-section volume fraction contours of the solid

圖6 水平截面固體顆粒體積分?jǐn)?shù)分布Fig.6 Horizontal cross-section volume fraction contours of the solid

通過上述分析可以認(rèn)為，在攪拌過程中，磷石膏固體顆粒受到兩個(gè)力的作用：由于地球引力作用所產(chǎn)生的重力和由于粘性作用而使流體在流動(dòng)過程中所產(chǎn)生的剪切力.重力使固體顆粒在運(yùn)動(dòng)過程中具有向下運(yùn)動(dòng)的重力加速度，剪切力使固體顆粒能夠隨液體流動(dòng).磷石膏固體顆粒在攪拌槽中最終所沉淀的位置是由這兩個(gè)力疊加作用的結(jié)果.在葉輪正下方由于存在一個(gè)倒錐體流型，此流型形成液流絕緣區(qū)，故此區(qū)域固體顆粒體積分?jǐn)?shù)較高；槽底壁面是整個(gè)陳化槽內(nèi)循環(huán)的轉(zhuǎn)向地點(diǎn)，又由于槽底與壁面接觸處為直角，故此處也會(huì)沉積固體顆粒；在液面上，由于離槳葉處較遠(yuǎn)，又由于磷酸為高粘度流體，因此液面附近流體流動(dòng)的速率較小，故而在液面附近區(qū)域重力起主導(dǎo)作用，又由于重力方向是指向陳化槽底面的，磷石膏固體顆粒在重力的作用下會(huì)向下運(yùn)動(dòng)，所以液面處固體顆粒體積分?jǐn)?shù)很低，在攪拌轉(zhuǎn)速不大的情況下此區(qū)域固體顆粒體積分?jǐn)?shù)幾乎為0.

3 結(jié) 語(yǔ)

a. 從陳化槽槽底固體顆粒的體積分?jǐn)?shù)的最大值為0.072 2，小于0.52，說明陳化槽底部沒有發(fā)生磷石膏固體顆粒堆積現(xiàn)象，從模擬結(jié)果所得出的固體顆粒體積分?jǐn)?shù)分布情況來看槽底的顆粒是處于完全懸浮狀態(tài)的，從而說明利用Zwietering的完全離底關(guān)聯(lián)式在以磷酸溶液和磷石膏固體顆粒為介質(zhì)的情形下適用于磷酸陳化槽攪拌設(shè)備臨界轉(zhuǎn)速的估算.

b. 從整個(gè)磷酸陳化槽內(nèi)液流的流型來看，磷酸陳化槽內(nèi)流動(dòng)以軸向流為主，液體從液面經(jīng)葉輪流下，在葉輪外力作用下沖擊陳化槽的槽底，然后從陳化槽的壁面流向液面，在整個(gè)循環(huán)流動(dòng)過程中，陳化槽利用葉輪所產(chǎn)生的沖擊力和流體的速度帶動(dòng)磷石膏固體顆粒流動(dòng)，從而達(dá)到攪拌混合的目的.

參考文獻(xiàn)：

[1] 張?jiān)氏? 磷肥及復(fù)合肥料工藝學(xué)[M]. 北京：化學(xué)工業(yè)出版社，2008.

[2] 王樹林，黃志良，羅欣，等.宜昌磷礦工藝礦相學(xué)研究[J].武漢工程大學(xué)學(xué)報(bào),2012，34(2)：55-59.

[3] 田興，石和彬，池汝安，等.湖北某中低品位硅質(zhì)磷礦工藝礦物學(xué)研究[J].武漢工程大學(xué)學(xué)報(bào),2012，34(3)：45-50.

[4] 歐陽(yáng)峰．新型穿流式攪拌槳在濕法磷酸反應(yīng)器中的應(yīng)用[J]．西南交通大學(xué)學(xué)報(bào)，2000,35(2):145-150.

[5] 陳志平,章序文,林興華,等．?dāng)嚢枧c混合設(shè)備設(shè)計(jì)選用手冊(cè)[M]．北京：化學(xué)工業(yè)出版社，2004.

[6] 鐘 麗,黃雄斌,賈志剛,等．固-液攪拌槽內(nèi)顆粒離底懸浮臨界轉(zhuǎn)速的CFD模擬[J]．北京化工大學(xué)學(xué)報(bào)，2003,30(6):18-22.

[7] 李鵬飛,徐敏義,王飛飛,等．精通CFD工程仿真與案例實(shí)戰(zhàn)[M]．北京：人民郵電出版社，2011.

[8] 舒安慶，劉凱，魏化中，等.生物反應(yīng)器內(nèi)攪拌器的優(yōu)化設(shè)計(jì)[J].武漢工程大學(xué)學(xué)報(bào)，2011，33(11)：66-69.

[9] He X， Yan S， Hu Y. Face recognition using Lapla-cianfaces[J]. IEEE Transaction on Patterm Analysis and Machine Intelligence,2005，27(3):328-340.

[10] 趙洪亮，張延安，張超，等．改進(jìn)Interming槳種分槽攪拌性能的數(shù)值模擬[J]．過程工程學(xué)報(bào)，2011,11(1):15-19.