基于Mixture模型的磷酸萃取攪拌器固液兩相流研究

張正陽，張國鵬

（中水北方勘測設計研究有限責任公司，天津 300000）

0 引言

在濕法制磷酸的過程中，強酸與磷礦石的反應發(fā)生在磷酸萃取槽內，而萃取槽中的攪拌器對于兩種物質是否可以充分接觸并完成反應起到至關重要的作用，攪拌器的優(yōu)劣直接與磷礦石的萃取率與固相物的二次攪拌率息息相關，這關系到磷酸最終的生產(chǎn)效率，合理的攪拌器類型選擇對于磷酸萃取效率有著決定性作用，對于整個磷酸制造過程與企業(yè)經(jīng)濟效益的實現(xiàn)都有著重要意義[1-2]。

在萃取磷酸的過程中，核心的參數(shù)除了契合的槽體大小之外，最重要的就是萃取槽攪拌器本身的攪拌效率。如果攪拌器運轉不暢或者效率較低，萃取產(chǎn)生的副產(chǎn)品磷石膏，原材料磷礦石及磷酸和其他礦物酸等無法充分接觸，會在萃取槽下部產(chǎn)生一定量的沉積，這會進一步影響整個制酸過程的效率，導致電動機效率和功率因數(shù)大大下降，直至設備無法正常運行。因此對于萃取槽不同類型攪拌器的分析及功率轉速的選擇都有著重要的意義[3]。

攪拌槽內的液體通過循環(huán)流動使物料充分混合，而攪拌器的形狀和運轉狀況是決定攪拌槽中液體流動狀態(tài)最基本的因素，攪拌槽內的動態(tài)機理為液體攜帶物料通過循環(huán)流動使其充分混合。而攪拌器的形狀和運轉狀況是決定攪拌槽中液體流動狀態(tài)最基本的因素。本文將通過數(shù)值模擬的方法對磷酸萃取過程中不同類型攪拌器（圓盤渦輪式、折葉型、雙折葉型）的攪拌性能進行分析。

1 基本理論

1.1 控制方程

任何流體運動都遵循質量守恒定律、動量守恒定律和能量守恒定律。一般不將能量守恒方程作為控制方程之一，從而通過減少方程個數(shù)來提高穩(wěn)定性、收斂性及節(jié)約時間成本。對不可壓縮流體，其質量和動量守恒方程為：

其中，動量守恒使用的是Navier-Stokes方程，同時湍流流場需要遵循湍流控制方程。

1.2 Mixture模型選擇

在FLUENT中提供解決多相流問題的模型主要有VOF、Mixture及Eulerian三種模型。其中VOF模型主要探討氣液多相流中自由表面流的問題，本文不做考慮。下面對Eulerian模型和Mixture模型進行探討。

Eulerian多相流模型在兩相中都對連續(xù)性、動量和能量方程進行了求解。每相都設置了體積分數(shù)表征方程式。Eulerian模型適用于從稀疏到密集的混合密度，也適用于從小到大的漿料顆粒。強耦合使得這個模型比Mixture模型使用起來更加困難。與Mixture模型相比，其計算精度方面雖然有些優(yōu)勢，但是其計算穩(wěn)定性差。

對于Mixture模型，將混合過程視作連續(xù)相和分散相的組合，模型的控制方程與均質模型類似，分散相對連續(xù)相的相位滑移通過平衡從密度差所產(chǎn)生的阻力來計算。模型的基本假定是基于一個局部平衡，其中的分散粒子相對于連續(xù)相總是以其終端速度進行移動。Mixture模型可以應用的范圍十分廣泛，只要內部力矩達到平衡，Mixture模型便可應用到不同速度、不同粒徑、不同密度的場合。它最適用于在液體中的小顆?；驓馀?。而磷酸萃取過程中磷酸鈣顆粒在濃磷酸中的運動便與這一點十分契合。因此本次計算選用Mixture模型[4-5]。

2 模型的建立及網(wǎng)格的劃分

本文將通過數(shù)值模擬的方法對磷酸萃取過程中圓盤渦輪式、折葉型、雙折葉型3種器的攪拌性能進行分析。3種攪拌器形狀如圖1～圖3所示。

圖1 圓盤渦輪式攪拌器

圖2 折頁型攪拌器

圖3 雙層折頁型攪拌器

先用CFD前處理軟件UG進行模型的建立，然后使用ICEM進行網(wǎng)格的劃分。之后輸出網(wǎng)格并在FLUENT 14.0軟件里進行計算。依據(jù)多重參考系的思路，將整個模型劃分為兩個部分。分別是攪拌器葉輪域及葉輪外域。對于兩個不同區(qū)域采取不同大小的思路來劃分網(wǎng)格，對于攪拌區(qū)域的網(wǎng)格進行加密以提高計算的精確性。經(jīng)過調試比較，以折葉型攪拌器為例，整體攪拌模型的網(wǎng)格數(shù)量達到180萬左右。劃分網(wǎng)格時對攪拌葉輪區(qū)域和葉輪外區(qū)域采用不同大小網(wǎng)格劃分，并對葉輪進行網(wǎng)格局部加密。網(wǎng)格劃分示意如圖4所示。

圖4 攪拌槽與折葉型攪拌器的網(wǎng)格劃分示意圖

在磷酸萃取槽中，主要的攪拌混合物為固相形式出現(xiàn)的磷石膏（其主要成分是硫酸鈣），以及以液相形式出現(xiàn)的濃磷酸。其中液相濃磷酸的密度取1650 kg/m3，黏度取44 mPa·s；一般濃硫酸與磷礦石反應過程中生成的磷石膏，80%部分的直徑在0.04 mm左右，因此取固相顆粒的直徑為0.04 mm，密度取2000 kg/m3。

為了直觀地了解對比兩種不同葉型的攪拌器對于磷酸萃取過程中的攪拌懸浮固相物的能力，本文在葉輪上方，攪拌槽水體以內，設置了一個體積分數(shù)約為1%的區(qū)域，這個區(qū)域中磷酸鈣的含量高達95%，整體反應物系中，固相顆粒的體積占總體積的0.95%。

使用FLUENT 14.0并用滑移網(wǎng)格技術（SG）完成動靜區(qū)域之間的數(shù)據(jù)傳遞匹配。其中旋轉葉輪區(qū)為動區(qū)域，剩余水體部分為靜止區(qū)域。使用Mixture模型來實現(xiàn)磷酸萃取過程中的兩相流模擬。應用PISO算法以及標準k-ε湍流模型。兩種葉型設置相同的初始條件：轉速選取60 r/min，時間步長選取0.004 s，為了全面地了解整個攪拌過程，選取了30個周期作為整體的攪拌時間。

3 結果分析

1）徑流式攪拌器代表圓盤渦輪式攪拌器與軸流式攪拌器代表折葉型攪拌器對比。由圖5可以看出，t=10 s時，折葉型攪拌器攪拌槽內的固相物由重力的作用向下擴散，而圓盤渦輪式攪拌器攪拌槽中的固相物已經(jīng)蔓延至槽底；t=14 s時，折葉型攪拌器攪拌槽中的固相物由攪拌器的混合作用有向槽頂向上流動的趨勢時，圓盤渦輪式攪拌器的固相物已經(jīng)大部分聚集于槽底，主要集中于葉輪下方。t=18～22 s的時間段，折葉型攪拌器攪拌槽內的固相物漸漸散開擴散于槽內而圓盤渦輪式攪拌器還有相當部分的固相物沉淀于葉輪下方，槽內其他地方的固相物也沒有完全打開；t=30 s時，折葉型攪拌器固體顆粒已經(jīng)相對均勻地分布于槽內，而圓盤渦輪式攪拌槽內的硫酸鈣顆粒依然集中于局限的區(qū)域。由此看出，在磷酸萃取過程中，折葉型攪拌槳對于固相物顆粒（硫酸鈣為主）的擴散混合作用要優(yōu)于圓盤渦輪式攪拌器。

2）折葉型攪拌器與雙層折葉攪拌器對比。由圖5可以觀察到，在t=10 s時，雙層折葉攪拌器與普通折葉攪拌器的固相物顆粒分布情況基本相似，均受到重力的作用，密度較大的硫酸鈣顆粒向槽底移動；當t=14 s時，雙層折葉攪拌器攪拌槽中的固相物有著明顯的從槽底向上運動的趨勢，而普通折葉攪拌器固體顆粒向上運動的趨勢并不明顯；t=18 s時，雙層折葉的攪拌槽中已形成了一個濃度相對集中的區(qū)域，且完全處于懸浮狀態(tài)，而折葉型攪拌器攪拌槽中還有相當一部分固相物十分接近槽底，但已有向上運動的趨勢；t=22 s時，雙層折葉的攪拌槽中固相物已分布較廣，固相物被分散于初始固相物的另一端，而單折葉型攪拌器攪拌槽中的固相物只有約為10%的部分到達槽的另一端；t=30 s時，雙層折葉攪拌器的攪拌槽內固相物分布已十分均勻，固相物顆粒基本平均分散于截面的每個角落，而單層折葉攪拌器槽內的硫酸鈣顆粒分布相對較差。

圖5 0～30 s內不同攪拌器萃取槽中固液分布對比圖

以轉速為n=60 r/min為例，普通折葉型攪拌器整體功率為7.89 W，而雙層折葉式攪拌器的運轉功率為9.15 W，對比普通折葉型攪拌器功率上升不到16%，從整體攪拌效果的改善來看，功率的增加值在可接受范圍之內。

4 結論

1）通過非定常的形式，利用Mixture模型，以隨時間推移的眼光對圓盤渦輪式攪拌器及折葉型攪拌器在介質為濃磷酸和磷石膏的混合物中進行了固液兩相流數(shù)值模擬，得出結論：磷酸萃取過程的相同時間內，折葉型攪拌器的攪拌效果優(yōu)于圓盤渦輪式攪拌器。

2）對比模擬計算了雙層折葉型攪拌器與折葉型攪拌器在磷酸萃取過程中的攪拌情況，發(fā)現(xiàn)雙折葉攪拌器的攪拌性能明顯優(yōu)于折葉型攪拌器。

（編輯 馬忠臣）

作者簡介：張正陽（1989—），男，碩士，工程師，主要從事水力機械設計工作。

收稿日期：2021-11-29