離心浮選機(jī)結(jié)構(gòu)對(duì)流場(chǎng)影響的數(shù)值模擬

魏月友，孫姣，陳文義

（1.河北工業(yè)大學(xué)化工學(xué)院，天津 300130；2.河北工業(yè)大學(xué)工程流動(dòng)與過(guò)程強(qiáng)化研究中心，天津 300130）

離心浮選機(jī)結(jié)構(gòu)對(duì)流場(chǎng)影響的數(shù)值模擬

魏月友1,2，孫姣1,2，陳文義1,2

（1.河北工業(yè)大學(xué)化工學(xué)院，天津 300130；2.河北工業(yè)大學(xué)工程流動(dòng)與過(guò)程強(qiáng)化研究中心，天津 300130）

為探討浮選機(jī)結(jié)構(gòu)對(duì)流場(chǎng)的影響，對(duì)整體浮選機(jī)流場(chǎng)的分布特性進(jìn)行了數(shù)值模擬．以自吸式離心浮選機(jī)為模型，使用計(jì)算流體動(dòng)力學(xué)軟件STAR CCM+，運(yùn)用Realizable k-湍流模型和VOF多相流模型，采用SIMPLE算法對(duì)速度和壓力進(jìn)行耦合求解，得到了速度、壓力、湍動(dòng)能以及氣液相含率的分布規(guī)律．模擬結(jié)果表明，吸氣量、噴嘴出口動(dòng)能和湍動(dòng)能值隨噴嘴直徑增大而減小，當(dāng)噴嘴直徑為19 mm時(shí)效果較好．外筒體流場(chǎng)存在較大的速度梯度，速度值從7～8m/s降至1m/s以下．礦化氣體在外筒體的彌散程度較低，入料初始階段氣含率均一度僅為20%～30%．

STAR CCM+；離心浮選機(jī)；噴嘴；全流場(chǎng)；數(shù)值模擬

0 引言

自吸式離心浮選機(jī)是礦物加工領(lǐng)域中分選細(xì)粒物料的重要設(shè)備之一，其綜合利用了重力場(chǎng)以及離心力場(chǎng)的優(yōu)點(diǎn)，具有單位處理量大、占地面積小、能耗小、浮選時(shí)間短、分選精度高和分選粒度下限低等特點(diǎn)，克服了傳統(tǒng)浮選中較長(zhǎng)浮選時(shí)間與較大浮選空間的限制，從而提高了生產(chǎn)能力[1-3]．常規(guī)浮選主要以碰撞礦化為主導(dǎo)，離心浮選則多種礦化方式并存：湍流碰撞礦化、微泡析出礦化和微泡兼并礦化等[4-8]．自吸氣式離心浮選機(jī)較好地滿足“動(dòng)態(tài)礦化、靜態(tài)浮選”的理念，但工廠生產(chǎn)運(yùn)行中存在著泡沫層不穩(wěn)定、顆粒跑粗嚴(yán)重致使噴嘴阻塞現(xiàn)象，影響分選效率．

由于離心浮選機(jī)內(nèi)部流場(chǎng)十分復(fù)雜、仍未摸清分布狀況，現(xiàn)有的研究建立普遍性理論和計(jì)算方法比較困難，并且周期較長(zhǎng)、成本較高．為了進(jìn)一步研究、改進(jìn)和優(yōu)化浮選機(jī)的結(jié)構(gòu)，提高自吸氣式離心浮選機(jī)的浮選效率，結(jié)合浮選機(jī)在選煤廠運(yùn)行狀況，本文對(duì)浮選機(jī)整機(jī)流場(chǎng)及速度、壓力、湍動(dòng)能等分布特性做了較為系統(tǒng)的模擬研究．

1 計(jì)算模型與工作原理

自吸式離心浮選機(jī)和氣泡發(fā)生器結(jié)構(gòu)如圖1所示，其工作原理是有壓礦漿自入料管進(jìn)入環(huán)形給料槽，然后分配給氣泡發(fā)生器，由噴嘴噴出，在混合室產(chǎn)生負(fù)壓，同時(shí)由吸氣口吸入大量空氣，空氣和礦漿充分混合，空氣被破碎成大量微細(xì)氣泡，研究表明[9-10]，氣泡發(fā)生器產(chǎn)生的氣泡尺寸小于0.384mm的約為95%，適用于微泡浮選、分選細(xì)粒物料．氣泡發(fā)生器內(nèi)極高的湍流強(qiáng)度促進(jìn)了氣泡和礦粒的充分混合接觸，氣泡和礦粒粘附，從而實(shí)現(xiàn)礦化．礦化后的礦漿切向進(jìn)入浮選槽外筒體，同時(shí)由于壓力降低使溶解在礦漿中的氣泡再次析出[11]，礦化氣泡在浮力與離心力場(chǎng)的綜合作用下，在浮選槽外筒體中做低速旋轉(zhuǎn)，升浮至內(nèi)筒體從精礦出口排出．親水矸石則不與氣泡發(fā)生作用，保留在礦漿中，沿浮選槽外筒體錐段迅速螺旋下降，由尾礦出口排出，實(shí)現(xiàn)了浮選槽中精礦和尾礦的分離．

本計(jì)算以汾西礦務(wù)局水峪煤礦自吸式離心浮選機(jī)為研究對(duì)象：結(jié)構(gòu)為柱錐形，外筒直徑2 000mm，高2 000mm，內(nèi)筒直徑500mm，柱段高550mm，錐段高550mm，由8個(gè)氣泡發(fā)生器礦化給料，切向給料點(diǎn)距外筒體柱段下邊緣為50 mm．浮選機(jī)網(wǎng)格模型示意如圖2所示，采用多面體網(wǎng)格，體網(wǎng)格數(shù)為2 262 723，坐標(biāo)原點(diǎn)位于柱錐段連接處Z=0平面．

圖1 離心浮選機(jī)結(jié)構(gòu)簡(jiǎn)圖及氣泡發(fā)生器局部圖Fig.1 Structurediagram of centrifugal flotation and bubblegenerator

圖2 浮選機(jī)網(wǎng)格模型示意圖Fig.2 Meshmodelof centrifugal flotation

2 控制方程及求解方法

2.2 求解方法與邊界條件

采用有限體積法計(jì)算連續(xù)介質(zhì)模型，選用VOF多相流模型捕捉相界面，由Realizable k-湍流模型封閉方程組．選擇分離求解器，基于SIMPLE算法求解流場(chǎng)中的壓力和速度耦合，對(duì)流項(xiàng)采用二階迎風(fēng)差分格式．

浮選機(jī)實(shí)際運(yùn)行時(shí)，存在著礦粒、氣泡、水、浮選藥劑等多相的流動(dòng)，各相之間又相互作用，還涉及到礦粒與氣泡間的碰撞、粘附、升浮或脫落，受計(jì)算流體動(dòng)力學(xué)軟件限制，目前不能提供一種多相流模型來(lái)完整描述浮選過(guò)程的粘附與脫落現(xiàn)象，因此必須通過(guò)簡(jiǎn)化來(lái)模擬，用擬均相方法處理氣液固三相流動(dòng)，將粘附顆粒的礦化氣泡視為連續(xù)相氣體，假設(shè)每個(gè)氣泡上粘附5個(gè)精礦顆粒，求得=59.1kg/m3、=5.84×104Pa s，通過(guò)研究氣液運(yùn)動(dòng)和分布狀況，考察流動(dòng)特性對(duì)浮選動(dòng)力學(xué)的影響．

礦漿入口為質(zhì)量流量進(jìn)口，可根據(jù)礦漿處理量計(jì)算，介質(zhì)設(shè)為水；氣體入口為停滯入口，初始狀態(tài)為大氣壓，通過(guò)形成負(fù)壓完成自吸；尾礦出口為壓力出口，根據(jù)初始液面高度計(jì)算，P=gh；精礦出口為大氣壓．其它邊界設(shè)為壁面，壁面按無(wú)滑移條件處理，并用壁面函數(shù)法求解．邊界條件設(shè)置見(jiàn)表1．

表1 邊界條件參數(shù)設(shè)置Tab.1 Parametersetting ofboundary conditions

3 模擬結(jié)果分析與討論

3.1 環(huán)形給料槽流場(chǎng)分析

環(huán)形給料槽作為礦漿分配器，為礦化裝置提供有壓礦漿．環(huán)形給料槽與氣泡發(fā)生器的8個(gè)連接出口，由入料管逆時(shí)針旋轉(zhuǎn)依次編號(hào)為1～8．表2為環(huán)形給料槽不同出口的質(zhì)量流量和壓力值，根據(jù)表中壓力報(bào)告值顯示，環(huán)形給料槽出口壓力基本保持均一狀態(tài)，礦漿壓力及流量的均一有利于液面和泡沫層的穩(wěn)定．

表2 不同出口的質(zhì)量流量與壓力值Tab.2 Mass flow and pressureatdifferentoutlets

3.2 氣泡發(fā)生器流場(chǎng)分析

3.2.1 氣泡發(fā)生器工作狀況與流場(chǎng)模擬

氣泡發(fā)生器是離心浮選機(jī)的關(guān)鍵部件，自吸氣、破碎成氣泡、氣泡礦化同時(shí)完成．該礦化裝置充分利用了射流泵原理，有壓礦漿自噴嘴噴出，根據(jù)伯努利方程，壓能向動(dòng)能轉(zhuǎn)化，混合室內(nèi)形成負(fù)壓，空氣經(jīng)吸氣管進(jìn)入混合室，在喉管內(nèi)的強(qiáng)湍流作用下，空氣與礦漿相互剪切、破碎，充分礦化[13-15]．根據(jù)速度矢量圖3可以看出，由于氣泡發(fā)生器自身結(jié)構(gòu)以及吸氣管與混合室的連接方式，混合室內(nèi)出現(xiàn)多處漩渦以及氣體對(duì)噴嘴的直接沖擊，對(duì)能量的損耗以及吸氣量有一定的影響．氣泡發(fā)生器的噴嘴直徑、收縮角、面積比和喉嘴距等參數(shù)都會(huì)影響礦化效果．根據(jù)選煤廠生產(chǎn)運(yùn)行中反饋，由于噴嘴直徑較小，顆粒跑粗致使噴嘴阻塞嚴(yán)重，因此增大噴嘴直徑的同時(shí)保證吸氣量、湍動(dòng)強(qiáng)度以及給入旋流的壓力及速度等流體動(dòng)力學(xué)性能尤為關(guān)鍵，設(shè)

計(jì)合理的氣泡發(fā)生器結(jié)構(gòu)有助于提高礦化效果，并減小使用中故障的頻率．

3.2.2 不同噴嘴直徑下氣泡發(fā)生器性能研究

選煤廠氣泡發(fā)生器的噴嘴直徑為18 mm，針對(duì)工廠反饋由于噴嘴直徑小引起的易堵塞現(xiàn)象，分別對(duì)（18mm、19 mm、20 mm、21mm、22mm）5組氣泡發(fā)生器充氣性能進(jìn)行模擬研究，吸氣量和噴嘴出口速度隨噴嘴直徑變化曲線如圖4所示．從圖4中可以看出隨著噴嘴直徑的增大，吸氣量與噴嘴出口速度呈現(xiàn)明顯的下降趨勢(shì)，吸氣量是表征氣泡發(fā)生器性能的重要參數(shù)，氣含率越大，浮選機(jī)內(nèi)氣體含量增大，氣體的比表面積越大，所能提供顆粒粘附的載體越多，浮選能力越強(qiáng)．噴嘴出口速度代表著氣泡發(fā)生器吸氣動(dòng)力，動(dòng)能與壓能的相互轉(zhuǎn)化，隨著噴嘴直徑加大，動(dòng)能減小，隨之形成的負(fù)壓能越小，吸氣動(dòng)力減弱．噴嘴直徑大于19 mm后，噴嘴出口速度下降較為明顯，噴嘴出口直徑對(duì)物料在噴嘴處的出口速度影響較為敏感．

圖3 氣泡發(fā)生器內(nèi)速度矢量場(chǎng)及局部放大圖Fig.3 Velocity vectordistribution of bubblegenerator and partially enlarged

圖5為不同噴嘴出口直徑下沿氣泡發(fā)生器縱向坐標(biāo)的湍動(dòng)能變化值，此處坐標(biāo)軸以噴嘴進(jìn)口為坐標(biāo)原點(diǎn)，縱軸向下為負(fù)值．從圖5中可以看出，從礦漿進(jìn)入氣泡發(fā)生器至噴嘴出口處，湍動(dòng)能逐漸達(dá)到最大值，氣液相在此處發(fā)生激烈的碰撞，進(jìn)入喉管后湍動(dòng)能逐漸降低，影響了三相混合及礦化．隨著噴嘴直徑加大，氣泡發(fā)生器內(nèi)的湍動(dòng)程度呈降低趨勢(shì)，大于19mm后下降顯著．湍動(dòng)程度是影響氣泡礦化中慣性碰撞機(jī)理的重要因素，細(xì)粒浮選中，脫落功較小，紊流度的增加有利于微細(xì)粒礦物的浮選．

圖4 吸氣量與噴嘴出口速度隨噴嘴直徑變化圖Fig.4 Inspiratory capacity and nozzleoutervelocity atdifferent nozzleouter diameters

圖5 不同噴嘴出口直徑下沿氣泡發(fā)生器縱向坐標(biāo)的湍動(dòng)能變化值Fig.5 Turbulentkineticenergy in z axisofbubblegeneratoratdifferent nozzle outer diameters

3.3 浮選機(jī)筒體流場(chǎng)分析

氣泡礦化后的礦漿切向入料至浮選槽，某切向入料管至浮選槽所建立探針監(jiān)測(cè)的速度和壓力變化值見(jiàn)圖6，由于直徑突然擴(kuò)大以及三相梯度發(fā)生突變，速度及壓力迅速下降、變化趨勢(shì)一致，到達(dá)外筒壁（徑向坐標(biāo)為1.0m）時(shí)，速度由7m/s降至1m/s以下．氣泡發(fā)生器形成負(fù)壓時(shí)所吸入的空氣，由于受到壓力作用，一部分溶解在礦漿中，當(dāng)進(jìn)入分選槽后，動(dòng)壓能向靜壓能轉(zhuǎn)變，壓力降低．礦漿中的過(guò)飽和空氣以微泡形式優(yōu)先于疏水礦粒上析出，從而使未礦化顆粒再次礦化．微泡析出浮選相比于傳統(tǒng)機(jī)械攪拌式充氣浮選，具有選擇性好、浮選速度快等優(yōu)點(diǎn)．

浮選槽縱截面（X=0）的礦化氣體含率如圖7所示，礦化氣泡進(jìn)入浮選槽外筒體后，在離心力及浮力的綜合作用下，迅速螺旋上升形成穩(wěn)定的泡沫層，而泡沫層液面呈內(nèi)凹錐形，因此無(wú)需刮板，礦化泡沫即可內(nèi)旋轉(zhuǎn)入內(nèi)筒體精礦槽，由精礦出口排出，通過(guò)調(diào)節(jié)精礦槽溢流堰高度即可調(diào)節(jié)精礦品味，而親水矸石則由尾礦出口排出．該結(jié)構(gòu)取消了常規(guī)浮選機(jī)中的刮泡裝置，結(jié)構(gòu)簡(jiǎn)化，動(dòng)力消耗減小．

圖6 切向入料口處速度和壓力值Fig.6 Velocity and pressurevaluesin tangential inletofouter flotation cell

圖7 縱截面X=0氣含率云圖Fig.7 Contour diagram ofairholdup at X=0

浮選槽外筒橫截面Z=0.05m、0.25m、0.5m速度云圖見(jiàn)圖8，顯示流場(chǎng)中存在較大的速度梯度，切線入料段進(jìn)入外筒體壁面處達(dá)到7～8m/s，靠近精礦收集槽內(nèi)筒體處降至1m/s左右，可能引起嚴(yán)重的偏流問(wèn)題以及氣泡兼并和液面翻花現(xiàn)象．另外可以看出，隨Z軸向上增大，速度梯度逐漸減?。V化氣泡由切向入料段進(jìn)入浮選槽后，與礦將混合速度有所減緩，隨液面上升，礦化氣泡慣性減小，速度梯度減?。?/p>

沿浮選槽Z軸方向氣含率均一度的變化趨勢(shì)如圖9所示，外筒的氣體彌散程度在入料初始階段并不大，氣含率僅為20%左右，隨著Z軸縱向高度的增大，外筒體中氣含率均一度呈逐漸上升的趨勢(shì)，礦漿中氣體分散度增大．空氣分散度越高，表明礦化氣泡在浮選機(jī)中分散越均勻．團(tuán)聚現(xiàn)象使礦化氣泡分布不均勻，易造成氣泡的兼并和細(xì)泥的機(jī)械夾帶，減小精礦的回收率．

圖8 Z=0.05m、Z=0.25m、Z=0.5m截面速度云圖Fig.8 Velocity contour diagram atthree crosssections (Z=0.05m,Z=0.25m,Z=0.5m)

圖9 氣含率均一度沿Z軸方向的變化Fig.9 Uniformity distribution of airholdup along the Z axis

4 結(jié)論

1）環(huán)形給料槽與氣泡發(fā)生器各連接出口流量和壓力基本處于均一狀態(tài)，利于液面及泡沫層的穩(wěn)定．2）通過(guò)對(duì)氣泡發(fā)生器不同噴嘴出口直徑下的流場(chǎng)模擬，可知噴嘴直徑增大，吸氣量、噴嘴出口速度和湍動(dòng)能值隨之減小，19mm時(shí)效果較好．同時(shí)氣泡發(fā)生器流場(chǎng)中由于漩渦的存在，造成了較大的能耗．

3）由于切向進(jìn)料，浮選槽外筒體空間，流場(chǎng)存在較大的速度梯度，速度值從外筒壁面處的7～8m/s到內(nèi)筒壁面降至1m/s以下．礦化氣體在外筒體的彌散程度較低，入料初始階段氣含率均一度僅為20%～30%．可適當(dāng)添加擋板或?qū)Я餮b置，避免較大的速度梯度及偏流等問(wèn)題．

4）浮選機(jī)流場(chǎng)的數(shù)值模擬能夠?yàn)榻Y(jié)構(gòu)優(yōu)化提供借鑒意義，并對(duì)指導(dǎo)工業(yè)生產(chǎn)具有一定的參考價(jià)值．

[1]王靜波，陳文義．離心浮選機(jī)的研究進(jìn)展[J]．選煤技術(shù)，2008（2）：288-318．

[2]方坤，陳文義，鄔領(lǐng)弟．離心力場(chǎng)浮選機(jī)械的研究[J]．煤礦機(jī)械，2005（12）：72-73．

[3]ChenWY．Aspirationspiral-flowtypecentrifugalflotationmachine[J]．JournalofCoalScience&Engineering，2002，8（1）：110-112．

[4]KohP，SchwarzM．CFDmodelofaself-aeratingflotationcell[J]．InternationalJournalofMineralProcessing，2007，85（1-3）：16-24．

[5]KohP，SmithL．Theeffectofstirringspeedandinductiontimeonflotation[J]．MineralsEngineering，2011，24（5）：442-448．

[6]呂玉庭，丁繼昌，孫庶．圓形離心浮選機(jī)的研究與應(yīng)用[J]．煤礦機(jī)械，2002（3）：16-18．

[7]PeterM，GraemeJ．Collisionofarisingbubble-particleaggregatewithagas-liquidinterface[J]．InternationalJournalofMineralProcessing，2014，130：1-7．

[8]YaG，GeoffreyM，EricaJ，etal．DEMsimulationofsinglebubbleflotation：Implicationsforthehydrophobicforceinparticle-bubbleinteractions[J]．AdvancedPowderTechnology，2014，25（4）：1177-1184．

[9]劉炯天，王永田．自吸式微泡發(fā)生器充氣性能研究[J]．中國(guó)礦業(yè)大學(xué)學(xué)報(bào)，1998，27（1）：29-33．

[10]LiL，LiuJT，WangLJ．Numericalsimulationofaself-absorbingmicrobubblegeneratorforacyclonic-staticmicrobubbleflotationcolumn[J]．MiningScienceandTechnology，2010，20（1）：88-92．

[11]胡為柏．浮選[M]．北京：冶金工業(yè)出版社，1983．

[12]王福軍．計(jì)算流體動(dòng)力學(xué)分析-CFD軟件原理與應(yīng)用[M]．北京：清華大學(xué)出版社，2004．

[13]閆小康，劉炯天，周長(zhǎng)春．旋流-靜態(tài)浮選柱管流段的兩相流數(shù)值模擬[J]．煤炭學(xué)報(bào)，2012，37（3）：506-510．

[14]桂夏輝，劉炯天，曹亦俊，等．氣泡發(fā)生器結(jié)構(gòu)性能的研究與進(jìn)展[J]．選煤技術(shù)，2009（2）：66-70．

[15]陸宏圻．射流泵技術(shù)的理論及應(yīng)用[M]．北京：水利電力出版社，1989．

[責(zé)任編輯 田豐]

Numerical simulation on theeffectof centrifugal flotationmachine structureon flow field

WEIYue-you1,2，SUN Jiao1,2，CHENWen-yi1,2

(1.School of Chemical Engineering,HebeiUniversity of Technology,Tianjin 300130,China;2.Research Center of Engineering Fluid and Process Enhancement,HebeiUniversity of Technology,Tianjin 300130,China)

In order to discuss the effectof the structureof Centrifugal flotationmachine on the flow field,the flow field distribution inw hole flotationmachinewas investigated using STAR CCM+.A Realizable k-turbulentand VOFmultiphasemodelw ereem ployed to simulate the presentmodel,and the SIMPLE algorithm wasused to acquire thepressurevelocity coupling solution.The simulation results indicate that inspiratory capacity,nozzleouter velocity and Turbulent Kinetic Energy decreased by increasingnozzleouter diameter,and thebetternozzle diameterwas19mm.High velocity gradientexisted inouter flotation cell,from 7～8m/s down to 1m/sor less.The dispersion ofmineral laden bubbles in outer cellwas relatively low,and at the initialstages itwasonly 20%～30%.

STARCCM+；centrifugal floatingmachine；nozzle；whole flow field；numericalsimulation

TD456

A

1007-2373(2015)01-0050-06

10.14081/j.cnki.hgdxb.2015.01.010

2014-10-13

河北省自然科學(xué)基金（E2012202109）

魏月友（1988-），男（漢族），碩士生．通訊作者：陳文義（1963-），男（漢族），博士，教授，E-mail：cwy63@126.com．

數(shù)字出版日期：2015-01-24數(shù)字出版網(wǎng)址：http://www.cnki.net/kcms/detail/13.1208.T.20150124.0953.006.htm l