數值仿真在自吸式浮選機開發上的應用

郭天宇， 姚 心， 李 鵬， 周鳳娟， 張 穎， 劉 誠

(1.中國恩菲工程技術有限公司， 北京 100038； 2.中國科學院過程工程研究所， 北京 100190)

0 引言

19世紀末，泡沫浮選正式成為一種工業選礦方法而得到應用[1]，經過100多年的發展，浮選技術和浮選設備得到了快速發展而日趨完善[2]。但是隨著近年來礦產資源的日趨匱乏和惡化，不斷對浮選機的設計和性能提出了新的要求[3-6]。

浮選機大型化過程中通過浮選機的設計原則、放大方法、設備結構參數等研究可以設計出大型浮選機，但是必須要通過試驗來驗證，這個過程周期長、費用大。為了縮短開發周期，降低研究費用，減少工業試驗工作量，需借助計算流體力學技術，建立浮選機流體動力學仿真模型，對大型浮選機的結構及不同操作條件下的流體動力學特性進行模擬計算[7-9]。本文利用CFD技術針對優化設計后的浮選機在設計工況下的運行性能進行驗證，通過礦漿流態、吸附壓力分布、結構設計以及功耗計算等一系列關鍵數據對浮選機內結構及運行參數進行優化和驗證。

1 計算模型

1.1 幾何模型

某項目浮選機由電機、傳動裝置、葉輪、定子、主軸、槽體、支架和電動推桿等組成。從CFD計算角度，浮選機可分為攪拌混合區和泡沫區2個區域。本文所有針對浮選機內流場數值模擬研究都是對2個區域范圍內的氣- 液兩相流場的考察。本文浮選機幾何模型在ANSYS Workbench中建立，選用ANSYS軟件中ICEM- MESHING模塊軟件進行網格劃分。

圖1 浮選機結構圖

圖2 浮選機模型

為了最大限度地減小浮選機整體模型的交互面的數量，并根據槽體內流體流動特征不同，將槽體內部分為空氣域、旋轉域和槽內區域，并分別進行網格劃分，同時進行數值模擬。浮選機主體結構分為旋轉區和固定區，其中葉輪- 定子旋轉區采用四面體非結構網格，其余部分劃分為混合型網格。

1.2 數學模型

浮選機模擬計算的質量守恒方程的表達式為：

(1)

對應的動量守恒方程為：

(2)

在Standardk-εModel模型中，關于k-ε的輸送方程如下：

(3)

(4)

上式中，Yk和Yw分別代表湍動生成的k和w，Sk和Sw為自定義源項。Gk代表的是由于平均速度產生湍流動能；計算如一般k-ε方程湍動生成一樣。Gb是由湍動形成的湍流動能。σk與σε分別是k和ε的湍動普朗特數。

2 邊界條件

對于浮選機模擬而言，邊界條件的正確設置也是仿真計算準確的關鍵。本文模擬計算的浮選機仿真邊界條件設置如表1所示。

表1 仿真邊界條件設置

3 結果與討論

3.1 浮選機全流體運行軌跡驗證

圖3所示為浮選機內全流體運行軌跡模擬結果，包括浮選機內礦漿粒子的運動情況，以及不同位置處的速度大小和分布。模擬結果可以直觀反映設備結構、尺寸和渦輪設計對礦粒在浮選機內運動的影響。礦漿被吸入浮選機后，隨轉子轉動速度達到最大值，隨后甩出而成徑向射流最后發展直至槽壁。礦漿粒子運行軌跡遍布整個浮選機槽體，這驗證了浮選機結構設計的合理性。

3.2 浮選機內速度場驗證

圖4所示為浮選機內中心剖面的液相速度云圖，作為浮選機槽內流體流動的主體，混合區內較強的液相速度有利于礦漿和氣泡的混合，而運輸區內液相速度較弱可使礦化氣泡穩定上升而不致因擾動過大而脫落。液相離開定子葉片后成為徑向射流，發展直至槽壁。圖4右側所示為葉輪- 定子區內橫截面速度云圖，葉輪葉片的旋轉帶動流體流動，故葉片迎風面均產生了很高的流體速度，定子葉片處流速有所降低，說明定子的穩流和導流設計有效。整體流場分布有一定的對稱性。

3.3 浮選機負壓形成及分布驗證

浮選機由電動機三角代傳動帶動葉輪旋轉，產生離心作用形成負壓，礦漿從而被吸入浮選機內。圖5所示為浮選機設備的壓力分布云圖，仿真結果顯示浮選機中心筒內和葉輪葉片上存在明顯的負壓，這對于浮選設備內上、下循環的產生、吸氣、吸漿作用等都是有益的。通過對浮選機內負壓的形成和分布的仿真模擬，從而對浮選機在設計工況下的運行性能進行了驗證。

圖4 浮選機內速度流場分布

圖5 浮選機內負壓形成及分布

模擬結果中紅色向藍色的變化表示壓力由高到低的變化過程。圖6中在葉輪的迎風面處存在高壓分布區，分別位于葉輪葉片轉角的前緣和上下部接近葉輪頂端的部位，這與工業現場葉輪磨損情況相吻合。同樣的，在定子表面的高壓區域是與葉輪高壓區相對的位于定子上部的迎風區域，這里承受的流場壓力要比其他區域大，磨損也更多，與實際中測試的情況也一致。

3.4 浮選機內氣泡分布及氣泡形成驗證

浮選機內氣、液兩相的分布情況對礦物浮選是十分重要的，也是評價設備性能和浮選操作優劣的重要指標。圖7為浮選機內液相和氣相的兩相分布情況。隨著葉輪的轉動，浮選機旋轉區形成的負壓不斷吸取空氣，進入旋轉區的空氣又不斷被旋轉葉片打碎分散。仿真內氣液兩相分布結果表明，氣泡能夠比較均勻的存在于整個浮選機設備內，并在液相中不斷渦旋上升，這說明礦漿與氣泡混合比較充分，有助于礦粒與氣泡的接觸、碰撞和黏附，有利于氣泡的礦化。通過分析計算浮選氣泡的產生過程和氣泡分布，從而對浮選機設計工況下的浮選性能進行了驗證。

圖6 浮選機轉子的受壓圖

圖7 浮選機內氣泡分布

3.5 浮選機擋板優化驗證

設計人員對浮選機內部結構及擋板進行了優化，通過仿真模擬對比了有無擋板兩種浮選機設計的氣液兩相分布及速度場分布。圖8a氣液兩相分布仿真結果表明，優化擋板后，浮選機吸氣能力和破碎生成氣泡量都明顯增加，從而有助于對礦粒的吸附和收集。圖8b速度場仿真對比結果表明，優化擋板后浮選機內速度分布更均勻，有利于混合區流場循環的穩定。通過對浮選機內部結構及擋板設計的優化、驗證，最大化提高了浮選機對礦漿的吸附和收集。

圖8 浮選機有無擋板模擬結果對比

3.6 浮選機渦輪理論功率驗證

通過浮選機模擬過程中，對渦輪剪切力、軸向受力等的數值監控和分析，進一步計算得浮選機在不同工況下的理論能耗值。圖9計算了礦漿工況下浮選機渦輪在轉動、吸液、吸氣液三種運行狀態下的理論能耗值，仿真結果表明浮選機運行初始時渦輪功率達到最大值，吸漿、吸氣運行平穩后渦輪功率值依次降低。渦輪能耗仿真結果與經驗設計值非常接近。同時我們對比了不同轉速下浮選機的渦輪理論能耗。圖9仿真結果表明，渦輪轉速越快能耗值越高。通過對浮選機渦輪的理論能耗計算，最終得到不同工況下浮選機的最優運行功率和能耗。

圖9 不同轉速下的理論功率值

4 結論

本文通過模擬計算浮選機內全流體運行軌跡，直觀驗證設備結構、尺寸和渦輪設計對礦粒在浮選機內運動軌跡的影響；通過浮選機工作過程中的氣液兩相流動、氣泡產生和均布、吸附壓力分布分析對浮選機在設計工況下的運行性能進行驗證；對浮選機內部結構及擋板進行了優化，提升對礦粒的吸附和收集；模擬計算浮選機在鎳礦礦漿中運行時的理論功率值，從而獲得浮選機不同工況下的最優運行功率和能耗。