戚 振

（泰安市特種設(shè)備檢驗研究院，山東泰安 271000）

0 引言

攪拌槳的功率準數(shù)是攪拌裝置設(shè)計所用到的最基本的參數(shù)[1]，它決定著設(shè)備的投入以及運行成本[2]。槳式攪拌器是最簡單的攪拌槳形式，在工業(yè)中的應用也比較普遍。組合槳的氣液攪拌可采用計算流體動力學（Computational Fluid Dynamics，CFD）方法進行研究。CFD是一種研究流體運動的輔助模擬計算軟件，可以清晰地呈現(xiàn)流體運動云圖，也可以分析熱傳導等相關(guān)特性，可較為真實地反映現(xiàn)場工況，根據(jù)現(xiàn)場工況設(shè)置邊界條件得出接近實際的計算結(jié)果。其具體操作步驟如下：在時域和空間領(lǐng)域，將非線性的連續(xù)場離散成為一系列有限點，這些有限點上的計算數(shù)值組合起來可近似為連續(xù)解，這樣就將連續(xù)體的無限解轉(zhuǎn)換為有限點上的近似解，離散點越多計算數(shù)值越精確，最終求出連續(xù)場的數(shù)值。

1 數(shù)值模擬方案設(shè)計

1.1 CFD理論

流體動力學涉及范圍廣，問題解決棘手，目前常用的計算方法有3種，分別為CFD數(shù)值計算、理論分析以及實驗模擬，此3種方法各有其適用的環(huán)境。理論分析方法有其優(yōu)缺點。其優(yōu)點是：計算結(jié)果的適用性，可以揭示流體流動的固有規(guī)律，是指導實驗測試分析和展示新的CFD數(shù)值模擬方法的理論依據(jù)。而其缺點是：為了解決非線性問題，分析結(jié)果往往只能在幾個流程問題中得到；理論分析要求研究人員必須掌握更高的數(shù)學能力和理論素質(zhì)。

CFD研究方法需要較少的時間和費用投入，可以很好地解決理論上的不足、試驗上的缺陷，計算結(jié)果可靠性高。所以，在流體動力學研究中，CFD方法正在發(fā)揮極為重要的作用[3]。例如，波音777飛機是一種“Paper Design”飛機，其機翼的繞流，大量的風洞測試和整機應力實驗，這些都是使用計算機完成CFD模擬和有限元等數(shù)值分析的，結(jié)果可以直觀地顯示在屏幕上，主要計算結(jié)果有各位置流場壓力的變化和分布，渦流分布和變化等。通過CFD方法模擬流體的流場，以FLUENT為例，解決方案通常包括如圖1所示步驟。

圖1 FLUENT數(shù)值模擬的通用步驟

1.2 槳葉區(qū)的處理

有3種方法來模擬槳葉區(qū)的攪動，分別為多重參考坐標系 法 （Multiple Reference Frame），滑 移 網(wǎng) 格 法 （Sliding Mesh）以及內(nèi)外迭代法（Inner-Outer Method）。多重參考坐標系法：首先設(shè)定兩個不相同的參考坐標系，1個動坐標系和1個靜坐標系，在旋轉(zhuǎn)動坐標系中計算攪拌槳存在的區(qū)域，在靜坐標系中計算不包括攪拌槳的區(qū)域，分別計算兩個坐標系中的數(shù)值[4]。轉(zhuǎn)動區(qū)域和靜止區(qū)域不重疊，交界面上的數(shù)據(jù)轉(zhuǎn)換可實現(xiàn)兩個區(qū)域的數(shù)據(jù)傳遞。設(shè)動區(qū)域為攪拌槳所在的區(qū)域，不包含攪拌槳的區(qū)域設(shè)為靜區(qū)域，則動區(qū)域選擇旋轉(zhuǎn)動坐標系作為坐標參考系，靜區(qū)域選擇靜坐標系為參考系。

多重參考坐標系法雖然是種穩(wěn)態(tài)流場的近似計算方法，但是在許多非穩(wěn)態(tài)流場計算中也得到了很好的應用[5]。滑移網(wǎng)格法與多重參考坐標系法的相同之處是都將計算流體區(qū)域分成兩塊，不同點是滑移網(wǎng)格法設(shè)置的兩個坐標系都是靜坐標系，包含攪拌槳的區(qū)域隨攪拌槳運動，其余靜止。攪拌槳區(qū)域與無攪拌槳區(qū)域重合部分會發(fā)生相對滑動，可利用插值函數(shù)進行質(zhì)量和動量的轉(zhuǎn)換匹配?；凭W(wǎng)格法多用于瞬態(tài)模擬，它的不足之處是計算量大，耗費時間長，需要高配置的虛擬環(huán)境進行運算。

內(nèi)外迭代法是Brucato在1994年提出的，這種方法是Brucato在深入研究“黑箱”模型之后提出迭代計算方法。該模型計算前需要設(shè)置2個可以互相重疊的區(qū)域，一般統(tǒng)稱為外區(qū)域和內(nèi)區(qū)域。外區(qū)域包括模型的內(nèi)壁、擋板、氣體分布器等結(jié)構(gòu)，內(nèi)區(qū)域主要包括攪拌葉槳。內(nèi)外迭代法的優(yōu)點是可直接進行迭代計算，不需要輸入實驗數(shù)據(jù)。所以這種方法有其明顯的缺點，收斂條件不容易控制，需要試差迭代[6]。

在以上3種方法中，多重參考坐標系法和滑移網(wǎng)格法模擬攪拌反應器內(nèi)轉(zhuǎn)動區(qū)域常用的兩種方法[7]。而由于本文在模擬通氣攪拌反應器內(nèi)流場的流體動力學問題時需要加入群體平衡模型（PBM）考慮氣泡的聚并與破碎對攪拌反應器內(nèi)部流場的影響，固此處用多重參考坐標系法（MRF）進行數(shù)值計算[8]。

2 仿真分析與結(jié)論

2.1 模型構(gòu)建

使用三維軟件Pro-E將建成的流體模型另存為x-t格式文件，然后將x-t格式文件導入Gambit進行網(wǎng)格劃分以及設(shè)置邊界條件。本課題中共設(shè)有3種攪拌槳類型，所以設(shè)定3種不同的模型進行分析，用Gambit軟件對3個模型進行網(wǎng)格劃分。實驗中攪拌反應器高度為1 000 mm，內(nèi)徑為380 mm，從底層到頂層共有3層攪拌槳，4塊擋板均勻分布在內(nèi)壁周圍，反應器的最低端是氣體分布器，直徑為160 mm，其上均勻分布16個孔徑為3 mm的氣孔。為提高網(wǎng)格劃分的質(zhì)量，并且縮短計算時間，在此選用四面體和六面體混合網(wǎng)格，并對攪拌反應器進行分塊處理。對靜區(qū)域進行切分處理，最大限度地使分割后的流體結(jié)構(gòu)規(guī)整以方便六面體網(wǎng)格的劃分，用四面體網(wǎng)格對不規(guī)則和動區(qū)域進行網(wǎng)格劃分，包保證網(wǎng)格的劃分質(zhì)量和計算精度，對攪拌槳做加密處理。圖2所示為網(wǎng)格質(zhì)量檢查截面圖，圖2（a）為低位組合槳網(wǎng)格軸截面圖，網(wǎng)格總數(shù)為762 786；圖2（b）為中位組合槳網(wǎng)格軸截面圖，網(wǎng)格總數(shù)為763 236；圖2（c）為高位組合槳網(wǎng)格軸截面圖，網(wǎng)格總數(shù)為755 241。由于最差網(wǎng)格質(zhì)量在動區(qū)域處，這3個模型的槳型完全一致，畫網(wǎng)格所用的尺寸也一致，因此，3個模型的最差網(wǎng)格質(zhì)量均為0.957。

圖2 網(wǎng)格軸向截面圖

2.2 模擬結(jié)果

層槳葉輪區(qū)的流體在槳葉的作用下形成徑向流動，流向攪拌反應器內(nèi)壁，在這個過程中混合了大量液體，在碰到反應器內(nèi)壁后分成兩部分，一部分沿內(nèi)壁向上，另一部分沿內(nèi)壁向下，最終會回流到槳葉區(qū)，底部的徑向槳葉圓盤將兩個循環(huán)分隔，形成“雙循環(huán)”運動。如圖3所示，處于中部位置槳葉的作用是控制渦心位置，保持旋轉(zhuǎn)穩(wěn)定進行，改變槳葉位置對下部流體的循環(huán)特性影響較小，即下循環(huán)不受槳葉位置的影響，對于上循環(huán)來說，中層槳葉位置升高，循環(huán)范圍開始擴大，在這3種情況下，底層槳的上循環(huán)作用也隨著上循環(huán)中心上移。中層槳為混合流式槳葉，上層槳為四寬頁軸流式槳，由于混合流式中層槳的存在，上層槳的軸流作用被放大，使得循環(huán)圈變大。如圖3（a）所示，當中層槳葉位置較低時，其與上層槳葉的距離較遠，上層槳葉受到的控制作用較小，釜體上層槳至液面處形成了一個較大的循環(huán)，這和單層軸向流型槳的流動有相似之處。當中層槳在中位時，如圖3（b）所示，中層槳與上層槳之間大循環(huán)的循環(huán)中心開始上移，上層槳的葉輪作用范圍要減小。當中層槳處于高位時，如圖3（c）所示，上層槳與中層槳連接地更加密切，因為上層4個寬頁攪拌槳的軸流作用與中層槳的混合流作用結(jié)合效果較強，反應器內(nèi)中上層流域軸向流效果比較明顯，在兩槳間形成一個速度比較大的軸向流動循環(huán)，這對混合比較有利。中層槳與上層槳的循環(huán)中心均位于對應攪拌槳的斜上方，盡管上層槳位置已經(jīng)固定，但由于循環(huán)中心隨著中層槳位置的上移，循環(huán)中心作用范圍減小，而且循環(huán)中心的位置隨著中層槳的上升而升高。

3 結(jié)束語

多層攪拌槳的攪拌過程中，上下槳之間距離的變化會引起流場循環(huán)變化。通過實驗驗證以及數(shù)值模擬對比，3層組合槳反應器內(nèi)均形成了循環(huán)流場，改變中層槳位置，查看中層槳位置變化對反應器內(nèi)流場的影響。中層槳位置變化引起3層組合槳在反應器內(nèi)產(chǎn)生了比較復雜的流場，當中層槳處于高位時，反應器內(nèi)中、上層槳之間的區(qū)域產(chǎn)生的流場受影響大，中層槳與上層槳位置較近，減弱了上層槳的軸向流作用，當中層槳位置下移時，中層槳對上層槳的抑制作用削弱，而且對流場產(chǎn)生了傳遞作用，將流場傳遞到下層攪拌槳，這使得形成的循環(huán)流場變大，底層槳下循環(huán)形態(tài)特征基本沒有變化。

圖3 速度場