陳年金

（溫州市工業(yè)科學(xué)研究院)

0 引言

攪拌是一種廣泛應(yīng)用的單元操作，普遍應(yīng)用于化工、食品、石油等行業(yè)中。攪拌過程就是通過槳葉的旋轉(zhuǎn)向攪拌槽內(nèi)輸入機(jī)械能，使流體獲得適宜的流場形態(tài)，從而在流場中進(jìn)行動量、熱量、質(zhì)量傳遞及化學(xué)反應(yīng)的過程。因此流場形態(tài)問題和攪拌能量問題一直是攪拌過程所研究的主要課題。目前，對于各種型式的攪拌器能提供什么樣的流動場、需供給多大的能量，并沒有很好的計算方法，大多是依據(jù)經(jīng)驗(yàn)公式的計算結(jié)果，并參考已投入生產(chǎn)的設(shè)備來決定攪拌器功率。對于為滿足工藝要求而開發(fā)的新型攪拌器，則更是沒有參考依據(jù)。

生產(chǎn)實(shí)踐中攪拌器功率不足的問題易于被查覺，而攪拌器功率過大造成浪費(fèi)的問題則容易被忽視，最理想的狀態(tài)當(dāng)然是攪拌器功率正好等于攪拌作業(yè)功率[1]。因此攪拌器功率的精確計算是至關(guān)重要的。

CFD技術(shù)的出現(xiàn)極大地促進(jìn)了攪拌過程的研究，它可以模擬任意型式攪拌器下的流場。本文將以廣泛應(yīng)用的六平直葉圓盤渦輪攪拌器為研究對象，利用CFD技術(shù)計算攪拌器功率，并對比經(jīng)驗(yàn)公式計算結(jié)果，為各種攪拌器的功率計算提供較為精確的方法。

1 模型建立

1.1 攪拌槽結(jié)構(gòu)

本文采用六平直葉圓盤渦輪，攪拌槽為圓柱形筒體，槽體邊緣均布四塊擋板，擋板寬度W為150 mm，攪拌槽直徑D為1500 mm，高度H為1700 mm，攪拌槳離槽底高度c為500 mm，dj/D=0.33， dj∶l∶b=20∶5 ∶4， 攪拌轉(zhuǎn)速 n 為 120 r/min，其中dj為攪拌槳直徑；l為槳葉長度；b為槳葉寬度。工作介質(zhì)為單一不可壓性的水。

1.2 網(wǎng)格劃分

在對攪拌槽進(jìn)行模擬時，解決運(yùn)動的槳葉和靜止的擋板之間的相互作用有不同的模型方法：黑箱模型法，動量源法，內(nèi)外迭代法，多重參考系法[2]和滑移網(wǎng)格法。本文采用的是多重參考系法。這種方法是將計算域分為兩個區(qū)域，內(nèi)部區(qū)域包括旋轉(zhuǎn)的槳葉，外部區(qū)域包括靜止的擋板。它與內(nèi)外迭代法不同之處是沒有重疊的部分，而是通過交界面將內(nèi)外流體區(qū)域連接在一起，流體參數(shù)的匹配直接通過交界面參數(shù)的轉(zhuǎn)換實(shí)現(xiàn)。

根據(jù)流動的對稱性，計算域選取了槽體的一半。采用ANSYS軟件劃分網(wǎng)格，將外部計算域進(jìn)行分區(qū)，采用六面體結(jié)構(gòu)化網(wǎng)格劃分模型，如圖1、圖2所示；內(nèi)部旋轉(zhuǎn)域采用非結(jié)構(gòu)網(wǎng)格劃分，如圖3所示。單元數(shù)分別為413 816和405 326，節(jié)點(diǎn)數(shù)分別為437 076和75 956。

圖1 靜止域網(wǎng)格頂視圖

1.3 邊界條件

合并后的網(wǎng)格如圖4所示。本文應(yīng)用CFX軟件進(jìn)行模擬計算[3]。槽內(nèi)自由液面定義為對稱邊界條件 （symmetric）；1/2槽體的軸向縱截面定義為周期性邊界條件 （rotational periodicity）；罐體、 攪拌軸、攪拌槳及擋板表面均采用無滑移壁面；內(nèi)外流體區(qū)域的分界面定義為交界面 （interface），以保證計算過程中內(nèi)外計算域相互耦合。

2 數(shù)值模擬

2.1 計算方法

圖2 靜止域網(wǎng)格

圖3 旋轉(zhuǎn)域網(wǎng)格

圖4 合并后的網(wǎng)格

CFX軟件采用全隱式耦合求解的技術(shù)，避免了傳統(tǒng)算法需要 “假設(shè)壓力項(xiàng)——求解——修正壓力項(xiàng)”的反復(fù)迭代過程，CFX的計算速度和穩(wěn)定性比傳統(tǒng)方法提高了許多。攪拌器與槽體間的相對運(yùn)動采用穩(wěn)態(tài)的多重參考系法進(jìn)行模擬，對流項(xiàng)的離散采用一階迎風(fēng)差分格式，湍流模型選用標(biāo)準(zhǔn)κ-ε模型，使用Auto Timescale，所有變量的收斂殘差均小于 1×10-4。

2.2 流場分析

圖5、圖6分別為六平直葉圓盤渦輪攪拌器的軸向速度矢量圖和橫截面速度矢量圖。其中軸截面選在擋板正中間的軸向平面，橫截面取位于葉片中部的平面。

圖5 軸向速度矢量圖

圖6 橫截面速度矢量圖

由圖5可以觀察到，六平直葉圓盤渦輪是典型的徑向流槳，葉片將流體沿徑向甩出，流體到達(dá)壁面形成上下兩個循環(huán)。葉輪區(qū)形成高速徑向射流，流體的速度在槳葉尖端處最大。從圖6可以看到,整個動區(qū)域流場的速度都很高，但攪拌槳內(nèi)部速度很小，在攪拌槳葉以外流體的速度逐漸降低。

2.3 攪拌器功率

攪拌器在流場中受壓力和黏性力作用，因而攪拌器的功率可以通過流體作用在槳葉上的扭矩進(jìn)行計算。由數(shù)值模擬得到攪拌器在x方向的扭矩值，如表1所示。所以攪拌器所受總扭矩為：

表1 1/2攪拌器的扭矩

對于旋轉(zhuǎn)的系統(tǒng)，功率為：

式中M——扭矩；

ω——角速度；

n——轉(zhuǎn)速。

通過數(shù)值模擬，可得到攪拌器的軸功率為：

3 公式計算法

對于渦輪式攪拌器，基于Rushton算圖[4-5]的公式法是比較精確的。計算攪拌功率前必須計算雷諾數(shù)：

式中ρ——介質(zhì)密度；

μ——介質(zhì)動力黏度，查化學(xué)化工物性數(shù)據(jù)手冊[6]，常溫下ρ=0.997×103kg/m3，

μ=0.893 7×10-3Pa·s；

dj——攪拌器直徑，dj=0.5 m；

n——攪拌器轉(zhuǎn)速，n=120 r/min=2 r/s。

代入數(shù)值，得到雷諾數(shù)為Re=2×0.52×0.997×103/(0.893 7×10-3)=5.58×105。 根 據(jù) 雷 諾 數(shù) ， 查Rushton算圖[1]，可得到相應(yīng)的功率準(zhǔn)數(shù)NP=6。因此，通過公式法計算得到的攪拌軸功率為：

將上述計算得到的軸功率與通過CFX軟件模擬得到的軸功率相比，其結(jié)果基本一致，模擬結(jié)果略小于公式計算的結(jié)果。因此，以數(shù)值模擬方法計算攪拌功率是合適的。

4 結(jié)論

本文針對六平直葉圓盤渦輪攪拌器，利用CFX軟件，模擬了攪拌槽內(nèi)的流場分布情況，并得出了攪拌軸功率，再應(yīng)用公式法進(jìn)行了計算驗(yàn)證，得到以下結(jié)論：

（1）用CFX軟件模擬出來的流場流向和規(guī)律跟實(shí)際情況相符，說明該模擬結(jié)果是合適的。

（2）利用模擬得到的數(shù)據(jù)計算的功率跟利用Rushton算圖計算的功率相差不大，說明模擬結(jié)果是比較可靠的，這對于那些不適用公式法計算的攪拌器提供了一個很好的設(shè)計參考。

[1]王凱，虞軍.化工設(shè)備設(shè)計全書：攪拌設(shè)備 [M].北京：化學(xué)工業(yè)出版社，2003.

[2]劉敏珊，張麗娜，董其伍.圓盤渦輪式攪拌槽的數(shù)值模擬 [J].鄭州大學(xué)學(xué)報，2007,28(1):122-124.

[3]馬青山，聶毅強(qiáng).攪拌槽內(nèi)三維流場的數(shù)值模擬 [J].化工學(xué)報，2005,54(5)：612-618.

[4]陳志平，章序文，林興華，等.攪拌與混合設(shè)備設(shè)計選用手冊 [M].北京：化學(xué)工業(yè)出版社，2002.

[5]李英慧.反應(yīng)釜攪拌功率的常規(guī)計算和軟件設(shè)計的比較[J].石油和化工設(shè)備，2008（6）:14-17.

[6]劉光啟，馬連湘，劉杰.化學(xué)化工物性數(shù)據(jù)手冊 [M].北京：化學(xué)工業(yè)出版社，2002.