夏國正　宋紅攀　郭明超

摘 要：通過對比攪拌容器的不同直徑、不同的槳葉結(jié)構以及攪拌轉(zhuǎn)速等參數(shù)，利用多相流模型進行攪拌器內(nèi)部流場仿真分析，以獲得并優(yōu)化攪拌混合器的攪拌效果，為設備設計及操作提供參考。結(jié)果表明，桶徑選取900 mm、攪拌轉(zhuǎn)速140～150 r/min;槳葉選用S型混合結(jié)構，且直槳段及S槳段選用適宜的傾角可兼顧排出性能、剪切性能、混合性能，并有較低的最高液面高度及較優(yōu)的攪拌效果;隨著下料速度的增加，水泥的混合性能會逐漸變差。

關鍵詞：桶外攪拌混合器;多相流;攪拌槳;攪拌流場仿真

中圖分類號：TL941? 文獻標志碼：A? 文章編號：1671-0797（2022）11-0032-05

DOI：10.19514/j.cnki.cn32-1628/tm.2022.11.009

0? ? 引言

桶外攪拌混合器是核設施中利用桶外水泥固化技術處理低、中放廢液的核心設備，其工作原理是將水泥及干粉添加劑與料液在攪拌器內(nèi)充分攪拌均勻后再將水泥漿排出，裝入固化桶內(nèi)進行固化及養(yǎng)護處理[1]。

本文對某工程中使用的400 L攪拌混合器的攪拌過程進行了數(shù)值模擬，分析了攪拌容器直徑、攪拌槳葉結(jié)構及尺寸、攪拌轉(zhuǎn)速等影響攪拌器內(nèi)部流場的因素，從而進行攪拌混合器結(jié)構與參數(shù)的優(yōu)化。

1? ? 流動控制物理模型

水泥干粉與料液攪拌的使用工況涉及兩相和三相計算（忽略攪拌過程中水泥與水的化學反應對攪拌流場的影響），在工程中采用的多相流模型主要包括VoF（Volume of Fluid）模型、混合（Mixture）模型、歐拉（Eulerian）模型等[2]。VoF模型主要用來模擬自由液面問題;混合模型是一種簡化的兩（多）相流模型，它使用單流體方法來模擬各相有不同速度的兩（多）相流;歐拉模型將連續(xù)相與分散相視為連續(xù)的一體，對每一相都建立動量方程和連續(xù)性方程，通過壓力和相間交換系數(shù)的耦合來計算求解[3-4]。

三相計算的水泥顆粒體積分數(shù)較大，計算量較大，選取混合模型作為多相流模型，既能保證對攪拌器流場的準確分辨，又能提高計算效率。

2? ? 數(shù)值模擬

攪拌流場涉及空氣—料液—水泥顆粒組成的三相，對于不同槳型特定工況下的流場運動特性對比，可簡化為兩相流處理，以提高效率;對于瞬態(tài)工況、轉(zhuǎn)速優(yōu)化分析等部分，研究目標為固體顆粒的相關特征，無法忽略氣液兩相對固態(tài)的影響，此時采用混合模型計算氣體、液體和固體顆粒的物理特征。

2.1? ? 幾何建模

選取不同槳葉模型進行對比計算，如圖1所示，簡化設備局部結(jié)構，可以得到計算域的完整模型。

2.2? ? 網(wǎng)格劃分

多重參考系（Multiple Reference Frame，簡稱MRF）法適用于穩(wěn)態(tài)計算，對非穩(wěn)態(tài)流場也能得到良好的模擬結(jié)果。攪拌器內(nèi)部流場結(jié)構復雜，可使用MRF法進行處理。將攪拌器內(nèi)部流場分為兩個區(qū)域：（1）攪拌槳及其周圍區(qū)域為動區(qū)域，包含高速運動的攪拌器及旋轉(zhuǎn)坐標系，設定動區(qū)域內(nèi)的流體速度和攪拌槳轉(zhuǎn)速相同;（2）靜區(qū)域部分為攪拌器內(nèi)部計算域減動區(qū)域，處于靜止坐標系。

靜止區(qū)域采用結(jié)構網(wǎng)格，旋轉(zhuǎn)區(qū)域采用帶邊界層的非結(jié)構網(wǎng)格，以提高計算精度。網(wǎng)格模型如圖2所示。計算域主體采用六面體網(wǎng)格、四面體網(wǎng)格、邊界層網(wǎng)格等，充分考慮交界面網(wǎng)格密度匹配、槳葉附近區(qū)域網(wǎng)格處理以及網(wǎng)格尺度連續(xù)性。

2.3? ? 評價標準

攪拌效果的性能評價主要圍繞攪拌的排出性能、剪切性能和混合性能三方面，具體包括：液相固含率、混合時間、攪拌功率、相同位置的徑向及軸向速度分布、剪切性能等。其中，液相固含率、混合時間是攪拌均勻性最直觀的體現(xiàn)，是相對更重要的評價指標，其次流場的速度分布也能從流場方面體現(xiàn)攪拌器的攪拌均勻性。

3? ? 結(jié)果與討論

3.1? ? 不同槳型

以相同介質(zhì)的穩(wěn)定工況，對不同槳型（2.1節(jié)所述）的攪拌流場進行分析。三種槳型在X=0截面的速度矢量分布如圖3所示。

分析攪拌桶同一高度截面的軸向、徑向、切向速度，得到以下結(jié)論：

（1）相同轉(zhuǎn)速下，直槳型軸向排出性能與徑向排出性能均較好，但剪切性能較弱，且液面過高，速度梯度較大，攪拌過程易引起飛濺。

（2）S型槳徑向排出性能較強，軸向排出性能弱，且軸心位置料分布少，容易出現(xiàn)混合不均勻現(xiàn)象，影響混合性能。

（3）S混合型槳軸向排出性能較強，徑向排出性能較弱，但剪切性能較強。

3.2? ? 不同桶徑瞬態(tài)工況

選取S型混合槳，采用氣液固三相模型（混合模型），在特定攪拌轉(zhuǎn)速下對攪拌桶直徑分別為850 mm、900 mm、950 mm時的攪拌性能參數(shù)和流場分布情況進行分析。

對比不同桶徑攪拌器模型在相同的混合時間內(nèi)達到的均勻程度，判斷桶徑的選取。選取Y=0截面，高度分別為300 mm和600 mm，距離軸心距離分別為0、0.5r、0.9r的6個監(jiān)測點，如圖4所示，監(jiān)測各點的樹脂體積分數(shù)隨時間的變化，樹脂占料液的體積分數(shù)變化如圖5所示。

隨著攪拌時間的推進，樹脂在三種攪拌器中的變化趨勢相近;由于攪拌導致空氣的混入，樹脂的體積分數(shù)有一定減小。從圖5中樹脂在料液中的體積分數(shù)變化可知，6個監(jiān)測點樹脂占料液的體積分數(shù)變化上限和下限的差別為0.02%左右，可見樹脂在料液中的均勻性沒有隨著料液加入水泥后整體物性的變化而受到過多影響，不同位置的體積分數(shù)相差不大，樹脂還是均勻分布在料液之中。隨著攪拌桶直徑的增大，模型的自由液面逐漸降低。5D75C0B2-D069-4725-8631-2297753DD30B

3.3? ? 攪拌速度分析

選取S型混合槳，針對不同轉(zhuǎn)速下水泥下料后到攪拌均勻狀態(tài)，開展空氣—水泥—料液三相數(shù)值模擬，分析判斷轉(zhuǎn)速對攪混均勻性的影響。以水泥加入后在液面上的落點區(qū)域作為初始分布。

選取Y=0截面，高度分別為200 mm、500 mm、800 mm，距離軸心距離分別為0、0.5r、0.95r的9個監(jiān)測點，如圖6所示，監(jiān)測水泥體積分數(shù)隨時間的變化。

選取攪拌75 s后的監(jiān)測數(shù)據(jù)，對比不同轉(zhuǎn)速下攪拌流場內(nèi)不同監(jiān)測點的水泥體積分數(shù)差別。Y=0截面水泥體積分數(shù)分布如圖7所示。

結(jié)果表明，轉(zhuǎn)速越大，攪拌越能快速趨近均勻狀態(tài);攪拌75 s后，9個監(jiān)測點的水泥體積分數(shù)趨近;轉(zhuǎn)速增加有助于攪混效率的提高，但轉(zhuǎn)速較高時（170 r/min、200 r/min），液面位置較高（分別達到1.15 m和1.19 m），容易造成飛濺。綜合考慮攪拌混合均勻性、攪拌液面高度等因素，攪拌槳轉(zhuǎn)速在140～150 r/min可較好地匹配下料與攪拌，且不會造成攪拌過程的飛濺。

3.4? ? 槳型角度參數(shù)的影響

在S型混合槳的基礎上，分析直槳段、S槳段不同角度值下的攪拌性能參數(shù)和流場分布情況，槳葉兩部分傾角示意圖如圖8所示，各槳在已有模型的基礎上對直槳段及S槳段型傾角進行調(diào)整，如表1所示。

圖9為各槳葉在X=0截面處速度矢量圖，從軸向、徑向、切向速度的變化趨勢分析兩段槳葉的傾斜角度對攪拌流場的影響：

（1）軸向速度：直槳段下壓之后，直槳段的斜率增大，攪拌過程中對泥漿軸向的沖擊增加;S段斜度可較好地對徑向及軸向速度進行分配，斜度過大或過小均會降低攪拌槳排出性能;攪拌槳軸向排出性能受S槳段影響較大，攪拌槳軸向排出性能與攪拌槳角度之間存在最佳角度使軸向排出性能最好。

（2）徑向速度：直槳段角度減小，攪拌過程中徑向作用力增加，從而使徑向速度增加;S段角度增加使攪拌槳與泥漿徑向接觸面積增大，攪拌過程中徑向作用力增加，從而使徑向速度增加。S槳段上抬+直槳段下壓，最大限度地增加了攪拌槳與泥漿的徑向接觸面積，從而使攪拌過程中徑向速度最大。綜合分析，攪拌槳徑向排出性能受S槳段角度影響較大，攪拌槳徑向排出性能隨著攪拌槳角度增大而呈逐漸增大趨勢。

（3）切向速度：直槳上抬后角度增大，攪拌過程易出現(xiàn)渦流，直槳段下壓后剪切角變小，剪切力下降;S槳段傾角過大會造成渦流，剪切角減小則剪切力下降。

4? ? 結(jié)論

此文應用FLUENT軟件，針對桶徑參數(shù)、槳葉結(jié)構及選型、攪拌轉(zhuǎn)速進行下料過程仿真研究，主要結(jié)論如下：

（1）關于桶徑：相同轉(zhuǎn)速下，攪拌桶直徑越大，自由液面越低;各桶徑下，固相在料液中的體積分數(shù)占比相差不大。

（2）關于轉(zhuǎn)速：隨著轉(zhuǎn)速增加，各點水泥體積分數(shù)與平均值間的差距逐漸縮小，轉(zhuǎn)速越高，攪拌器內(nèi)流場越能快速趨近攪混均勻;轉(zhuǎn)速增加有助于攪混效率的提高，但轉(zhuǎn)速過高，液面位置較高，容易造成攪拌過程中的飛濺。

綜合混合性能、液面高度等因素，桶徑選取900 mm，攪拌轉(zhuǎn)速為140～150 r/min。

（3）關于槳型：直槳槳型的排出性能較好，但剪切性能差;S型槳徑向排出性能較強，軸向排出性能弱，且軸心位置料分布少，容易出現(xiàn)混合不均勻現(xiàn)象，影響混合性能;S型混合槳軸向排出性能較強，徑向排出性能較弱，剪切性能較強。

（4）關于槳型傾角：攪拌槳軸向排出性能受S槳段角度影響較大，軸向排出性能隨攪拌槳角度增加，先增加后減小;攪拌槳徑向排出性能受S槳段角度影響較大，徑向排出性能隨著攪拌槳角度增大而逐漸增大;攪拌槳剪切性能受直槳段角度影響較大，攪拌槳剪切性能與攪拌槳角度存在最佳角度使攪拌槳剪切性能最好。

綜合以上，丙C型槳葉可綜合兼顧排出性能、剪切性能、混合性能及液面高度，為較優(yōu)的槳型選擇。

[參考文獻]

[1] 劉國彪，劉天斌，石兵.新型桶外混合器的理論與實驗研究[J].核標準計量與質(zhì)量，2012（3）：46-54.

[2] 韓占忠，王敬，蘭小平.FLUENT流體工程仿真計算實例與應用[M].北京：北京理工大學出版社，2010.

[3] 密路祥，徐立軍，胡兵.基于Fluent的水電解槽電場數(shù)值模擬[J].科學技術與工程，2021，21（21）：8883-8888.

[4] 周俊杰，徐國權，張華俊.FLUENT工程技術與實例分析[M].北京：中國水利水電出版社，2010.

收稿日期：2022-03-21

作者簡介：夏國正（1985—），男，江西上饒人，碩士研究生，高級工程師，研究方向：核化工非標機械設備。5D75C0B2-D069-4725-8631-2297753DD30B