摘 要： 濃縮接受釜作為一種常用的化工設(shè)備，其內(nèi)部固液懸浮均勻性對下游生產(chǎn)穩(wěn)定性至關(guān)重要。應(yīng)用多重參考系法和混合多相流模型對濃縮接受釜進行計算流體力學(xué)（CFD）模擬，分析了現(xiàn)有攪拌條件下固液懸浮狀態(tài)，并探究了攪拌轉(zhuǎn)速對固液懸浮的影響，最后提出了槳葉優(yōu)化設(shè)計方案。果表明：現(xiàn)有3層斜葉槳軸流推動能力弱，無法在釜內(nèi)形成主體循環(huán)，導(dǎo)致固液懸浮均勻性差，且通過增加攪拌轉(zhuǎn)速無法滿足生產(chǎn)要求。優(yōu)選的雙層軸流槳軸流推動能力強，可形成中心向下、邊壁向上的主體循環(huán)流動，有效提高了釜內(nèi)固液懸浮均勻性，保證了下游生產(chǎn)運行的穩(wěn)定。

關(guān) 鍵 詞：固液懸浮；混合；CFD；多相流；多層槳

中圖分類號：TQ027.2 文獻標志碼： A 文章編號： 1004-0935（2024）09-1465-05

固液懸浮屬于典型的非均相攪拌混合，其目的是借助攪拌器的作用，使固體顆粒懸浮于液相中，形成固液混合物或懸浮液^[1]，廣泛應(yīng)用于石化、制藥、冶金等領(lǐng)域^[2]。在結(jié)晶過程中，良好的固液懸浮狀態(tài)是保證結(jié)晶產(chǎn)品質(zhì)量的關(guān)鍵^[3]；在冶金工業(yè)中，維持固液懸浮是調(diào)漿浮選過程的重點，可提升浮選藥劑的利用率^[4]；在化學(xué)反應(yīng)過程中，顆粒懸浮可增大相間接觸面積、強化傳質(zhì)，以提升反應(yīng)效率^[5]。濃縮接受釜作為常見的化工設(shè)備，在接收固液漿料的同時，需將釜內(nèi)漿料均勻輸送至下游生產(chǎn)單元。釜內(nèi)通過攪拌槳維持固液懸浮狀態(tài)，但攪拌槳設(shè)計往往依賴經(jīng)驗，運行過程中易出現(xiàn)出料固含率不穩(wěn)定，因此深入研究濃縮接受釜內(nèi)固液懸浮過程對提高下游單元運行穩(wěn)定具有重要意義。

隨著仿真技術(shù)的快速發(fā)展，計算流體力學(xué)模擬（CFD）已廣泛應(yīng)用于固液懸浮過程的研究^[6]。湯俊杰等^[3]采用CFD模擬的方法對結(jié)晶器內(nèi)攪拌流場進行了研究，優(yōu)化了槳葉類型、槳葉層數(shù)以及安裝高度。茍炳中等^[4]采用CFD模擬了強力攪拌機內(nèi)的固液混合過程，并進一步優(yōu)選了攪拌轉(zhuǎn)速等操作參數(shù)。胡越等^[7]則基于流體力學(xué)模擬對比了釜底形狀對固液懸浮過程的影響。JIANG等^[8]采用CFD模擬的手段，研究了槳葉安裝高度對固液混合過程的影響，同時優(yōu)選了槳葉直徑等參數(shù)。

以工業(yè)濃縮接受釜為模型，采用CFD對現(xiàn)有釜內(nèi)固液懸浮狀態(tài)進行模擬，并探究操作參數(shù)對固液懸浮狀態(tài)的影響，最后基于模擬結(jié)果對現(xiàn)有攪拌型式進行優(yōu)化設(shè)計。

1 計算模型和模擬方法

1.1 計算模型

1.1.1 濃縮接受釜模型

工業(yè)上運行的濃縮接受釜三維模型如圖1所示，釜體直徑為3.2 m，運行液位高度為3.6 m，上部為進料口，釜底為出料口，物料連續(xù)進出體積流量均為40 m³·h^-1。入口處漿料固含率為1.7%，顆粒平均粒徑為0.25 mm，顆粒密度為2200 kg·m^-3，液相密度為1300 kg·m^-3，黏度為1 mPa·s。

釜內(nèi)攪拌器為3層槳式攪拌器，其中上2層斜葉上翻、最底層斜葉下翻。槳葉直徑1.2 m，槳葉寬為0.2 m，斜葉角度為45°，相鄰2層槳間距為0.8 m，底層槳離底高度為0.35 m。攪拌器順時針轉(zhuǎn)速，設(shè)計轉(zhuǎn)速為75 r·min^-1，攪拌電機裝機功率為11 kW。

1.1.2 攪拌區(qū)域設(shè)置及網(wǎng)格劃分

為模擬濃縮接受釜內(nèi)攪拌槳葉轉(zhuǎn)動對固液懸浮的影響，需對接受釜進行區(qū)域劃分，攪拌槳葉附近區(qū)域為轉(zhuǎn)動區(qū)域，其余為靜止區(qū)域，如圖2（a）所示。在轉(zhuǎn)動區(qū)域和靜止區(qū)域設(shè)置交界面，交界面邊界類型為Interface，從而實現(xiàn)轉(zhuǎn)動區(qū)域和靜止區(qū)域計算結(jié)果的數(shù)值交換^[9]。采用非結(jié)構(gòu)化網(wǎng)格對濃縮接受釜三維模型進行網(wǎng)格劃分，為了提高計算精度，對攪拌槳葉所在轉(zhuǎn)動區(qū)域進行網(wǎng)格加密^[10]，如圖2（b）所示，整體模型網(wǎng)格數(shù)約為170萬。

1.2 模擬方法

利用Fluent19.0模擬濃縮接受釜內(nèi)固液懸浮。采用多層參考系法（MRF）定義轉(zhuǎn)動區(qū)域，采用Mixture模型對固液兩相進行計算，液固間曳力模型采用Wen-Yu模型，湍動模型選用考慮旋轉(zhuǎn)修正的RNG模型^[11]。求解過程采用PISO壓力-速度耦合方法，離散方法選用二階迎風(fēng)格式，殘差設(shè)置為0.0001，起始計算步長為0.0001s，計算步長根據(jù)收斂情況動態(tài)調(diào)整。為了加快模擬進度，在初始化階段，將濃縮接受釜內(nèi)固含率設(shè)置為1.7%。

2 模擬結(jié)果分析

2.1 固液懸浮過程分析

基于上述CFD模型以及求解參數(shù)設(shè)置，模擬濃縮接受釜在連續(xù)進出料過程中固液懸浮過程的變化，攪拌轉(zhuǎn)速為75 r·min^-1。模擬過程監(jiān)測攪拌槳攪拌功率以及進出口固體質(zhì)量流量，結(jié)果如圖3所示。

由圖3（a）可知，在模擬監(jiān)測時間范圍內(nèi)攪拌功率趨于穩(wěn)定，最終攪拌功率在5.3 kW左右。根據(jù)前期調(diào)研，濃縮接受釜實際運行功率在5.6 kW左右，即采用CFD模擬得到的攪拌功率與實際運行功率相接近，由此證明CFD模擬結(jié)果的可靠性。由圖3（b）可知，隨著攪拌的進行出口固體量高于進口固體量，表明顆粒在釜底附近出現(xiàn)堆積，造成出口附近固含率升高，進而影響下游生產(chǎn)的穩(wěn)定性。

為分析攪拌過程固液懸浮過程的變化，提取不同時間段顆粒濃度分布，結(jié)果如圖4所示。由圖4可知，隨著攪拌的進行，濃縮接受釜內(nèi)顆粒相向下沉積，釜底封頭區(qū)顆粒量不斷增多，最終形成圖4（c）所示的顆粒濃度分布。圖4（c）所示的固液懸浮均勻性較差，底部顆粒濃度高、頂部顆粒濃度低，攪拌槳葉周圍區(qū)域顆粒濃度小。

為進一步剖析固液懸浮效果差的原因，提取=117s時垂直剖面上速度分布及速度矢量分布，結(jié)果如圖5所示。由圖5（a）可知，槳葉轉(zhuǎn)動過程中，槳葉葉端附近形成高速區(qū)；由圖5（b）可知，槳葉對流體的推動范圍較小，僅在槳葉附近出現(xiàn)明顯的速度矢量，且上2層槳葉起到斜向上推動的作用，底層槳葉起到斜向下推動的作用?，F(xiàn)有攪拌器在75r·min^-1條件下，無法在濃縮接受釜形成整體環(huán)流，進而無法將沉積至釜底的顆粒均勻分散。

2.2 攪拌轉(zhuǎn)速對固液懸浮的影響

為改善釜底顆粒堆積區(qū)域顆粒分散性，擬提高攪拌轉(zhuǎn)速。采用CFD模擬探究110、150 r·min^-1條件下的固液懸浮特性。圖6為不同轉(zhuǎn)速下濃縮接受釜垂直剖面上顆粒濃度分布。由圖6可知，隨著轉(zhuǎn)速的增加，釜底沉積的顆粒量變少，釜內(nèi)顆粒分布均勻性有所增加。

提取不同轉(zhuǎn)速下垂直剖面速度及速度矢量分布，結(jié)果如圖7所示。由圖7（a）可知，隨著攪拌轉(zhuǎn)速的增加釜內(nèi)流場速度整體增加，槳葉周圍高速區(qū)面積增加明顯。由圖7（b）可知，攪拌轉(zhuǎn)速增加后槳葉附近的速度矢量更密度且數(shù)值增大，槳葉推動范圍更廣，由此改善了釜底顆粒濃度分布^[12]。

統(tǒng)計不同轉(zhuǎn)速下接受釜出口固體質(zhì)量流量以及相應(yīng)攪拌功率，結(jié)果如圖8所示。由圖8可知，隨著攪拌轉(zhuǎn)速的增加，出口固體質(zhì)量流量逐漸減小，同時攪拌功率大幅增加。當(dāng)轉(zhuǎn)速為110 r·min^-1時，出口質(zhì)量流量為0.36 kg·s^-1，攪拌功率為17 kW，此時出口固體質(zhì)量流量與進口固體質(zhì)量流量相接近，但攪拌功率已經(jīng)超過裝機功率，且此時釜底仍有較高的固體沉積。當(dāng)轉(zhuǎn)速達到150 r·min^-1時，出口質(zhì)量流量為0.16 kg·s^-1，攪拌功率為48 kW。在轉(zhuǎn)速為150 r·min^-1條件下，釜底沉積的固體量減少，但此時出口固體質(zhì)量嚴重偏低，且此時攪拌功率遠遠超過電機功率。由上述分析可知，現(xiàn)有攪拌條件下，無法通過提高轉(zhuǎn)速改善濃縮接受釜運行穩(wěn)定性。現(xiàn)有運行攪拌器為3層斜葉槳，斜槳葉以剪切力為主，軸流推動能力較弱；與此同時，上2層槳葉斜葉上翻、底層槳葉斜葉下翻，軸流推動方向不一致，使得現(xiàn)有攪拌器無法在釜內(nèi)形成中心向下邊壁向上的大循環(huán)，進而導(dǎo)致固液懸浮效果較差。

2.3 槳葉類型對固液懸浮的影響

為改善濃縮接受釜內(nèi)固液懸浮效果，提出圖9所示推進式攪拌器。該攪拌器采用雙層槳葉，每層均為下翻式三葉槳，層間高度1 m，槳徑1.2 m。槳葉離底距離0.75 m，槳葉根部寬0.24 m，安裝折角為45°；槳葉末端寬0.17 m，折角為20°。

采用CFD模擬圖9所示三維模型固液懸浮過程，攪拌轉(zhuǎn)速為75 r·min^-1，得到圖10所示模擬結(jié)果。由圖10可知，ttZ/incG478yctUpon6b88UIzhTtZqUfMLHT7sIvKXQ=槳葉類型更換為兩層推進式槳葉后，固液懸浮均勻性大幅改善，僅有少量顆粒在釜底出料口附近出現(xiàn)沉積。攪拌槳葉功耗為4 kW，出口處固體質(zhì)量流量為0.42 kg·s^-1。由此可知，變更槳葉類型后一方面能夠降低攪拌功耗，另一方面可以保持均勻的固液懸浮狀態(tài)，提升后續(xù)反應(yīng)單元運行穩(wěn)定性。

進一步提取穩(wěn)定后垂直剖面速度及速度矢量分布，得到圖11。圖11（a）所示的速度分布表明，在兩層軸流槳的作用下，釜內(nèi)速度場分布基本均勻，僅在槳葉末端出現(xiàn)少量高速區(qū)；圖11（b）為速度矢量分布圖，此時全釜內(nèi)均有速度矢量出現(xiàn)，并基本形成中心向下邊壁向上的大循環(huán)，從而避免顆粒在釜底大量沉積，提升了固液懸浮的均勻性^[13-14]。

3 結(jié) 論

采用CFD模擬研究了現(xiàn)有濃縮接受釜內(nèi)固液懸浮過程，并通過攪拌功率模擬值和實際值的對比驗證了模擬結(jié)果的可靠性，同時考察了攪拌轉(zhuǎn)速對固液懸浮過程的影響，最后結(jié)合模擬結(jié)果提出了優(yōu)化的槳型設(shè)計方案。

1）現(xiàn)有3層斜葉槳軸流推動力較弱，且上2層推動方向與底層推動方向相反，無法在釜內(nèi)形成整體環(huán)流，從而造成釜底顆粒沉底嚴重，導(dǎo)致進出口固體質(zhì)量流量偏差較大，影響下游生產(chǎn)運行穩(wěn)定。

2）針對現(xiàn)有3層斜葉槳，增大攪拌轉(zhuǎn)速可改善固液懸浮效果，但攪拌功耗遠超現(xiàn)有裝機功率，且出口處固體質(zhì)量流量小于進口，同樣會引起下游生產(chǎn)單元運行不穩(wěn)定。

3）基于前期模擬分析，將槳葉類型優(yōu)化為雙層軸流槳，在降低攪拌功耗的同時，有效提升了槳葉軸流推動能力，在釜內(nèi)基本形成了中心向下邊壁向上的大循環(huán)，大幅提升了固液懸浮均勻性，實現(xiàn)了進出料固體質(zhì)量流量的平衡。

參考文獻：

[1] 周勇軍，盧源，陳明濠，等. 改進型INTER-MIG攪拌槽內(nèi)固液懸浮特性的數(shù)值模擬[J]. 過程工程學(xué)報，2014，14（5）：744-750.

[2] 楊鋒苓，周慎杰. 攪拌固液懸浮研究進展[J].化工學(xué)報，2017，68（6）： 2233-2248.

［[3］]湯俊杰，易爭明，郭璐.基于CFD對DTB結(jié)晶器內(nèi)攪拌槳構(gòu)型的優(yōu)選[J].現(xiàn)代化工，2019，39（1）：209-212.

［[4］]茍炳中，孫春寶，寇玨.強力攪拌機攪拌調(diào)漿固液混合數(shù)值模擬[J].礦業(yè)研究與開發(fā)，2022，42（3）：134-139.

[5] 許雯婧，王璐，甄衛(wèi)軍，等. 聚苯硫醚聚合攪拌釜結(jié)構(gòu)優(yōu)化及流場分析[J]. 化學(xué)反應(yīng)工程與工藝，2020，36（3）：203-212.

[6] 李良超，徐斌，楊軍. 基于計算流體力學(xué)模擬的下沉與上浮顆粒在攪拌槽內(nèi)的固液懸浮特性[J]. 機械工程學(xué)報，2014，50（12）：185-191.

［[7］]胡越，張一敏，薛楠楠，等.攪拌槽形對含釩頁巖固液混合特性影響的數(shù)值模擬[J].有色金屬（冶煉部分），2022（11）：68-77.

[8] JIANG H W， YUAN S， LIU H， et al. Numerical analysis and optimization of key parts in the stirred tank based on solid-liquid flow field[J]. ， 2022， 105（1）： 003685042110672.

［[9］]袁慎峰，金玲貞，陳志榮，等.基于CFD的攪拌萃取塔結(jié)構(gòu)優(yōu)化與流場分析[J].過程工程學(xué)報，2023，23（5）：681-690.

[10] 劉麗艷，李晨辰，朱國瑞，等. 粉煤灰-氫氧化鈉漿液槽放大準則數(shù)值模擬研究[J]. 現(xiàn)代化工，2014，34（2）：148-151.

[11] 李青云. 基于CFD的小型反應(yīng)釜中不同湍流模型數(shù)值模擬比較[J]. 當(dāng)代化工，2020，49（7）：1483-1487.

[12] 曹小爽，鄭海燕，王琦，等. 濕法提釩攪拌反應(yīng)釜的流場分布[J]. 鋼鐵釩鈦，2021，42（4）：6-11.

[13] 肖志鵬，陶蘭蘭，陳帥，等. 帶穩(wěn)定翼四斜葉-Rushton組合槳釜內(nèi)流場實驗研究[J]. 化工機械，2022，49（5）：756-763.

［[14］]趙珣，肖志鵬，林偉振.基于計算流體動力學(xué)的改進型三斜 Rushton組合槳攪拌釜內(nèi)流場研究[J].石油化工設(shè)備，2022，51（4）：15-20.

CFDSimulation ofSolid-LiquidSuspension in

ConcentrateReceivedStirredTank

LIANGXiaofeng，YANG Wenqiong， LIU Xiaobang， ZI Can， LI Jianchang

（Zhejiang Institute of Tianjin University（Shaoxing）， Shaoxing Zhejiang 312300， China）

Abstract：Solid-liquid suspension uniformity in concentrate received stirred tank is the key factor to determine the stability of downstream operation. In the present work， the solid-liquid suspension process and the effect of stirred speed were simulated by computational fluid dynamics using multiple reference frame method and mixture multiphase flow model. The simulated results showed that the axial flow capacity of current three-layer agitator was not enough， leading to solid deposition. Besides， the non-uniform suspension condition was not be solved by increasing the stirred speed. Furthermore， an optimized two-layer axial agitator with strong axial flow capacity was proposed， which could satisfy efficient bulk circulation in the tank and relatively uniform distribution of solids.

Key words：Solid-liquid suspension; Mixing; CFD; Multiphase flow; Multilayer agitator