彭小奇，陳思超，宋彥坡，劉 濤

(1.中南大學(xué) 能源科學(xué)與工程學(xué)院，長沙 410083；2.湖南第一師范學(xué)院 信息科學(xué)與工程系，長沙 410205)

泡沫浮選工藝被廣泛應(yīng)用于選礦生產(chǎn)。在浮選槽中，攪拌裝置對(duì)礦漿進(jìn)行激烈攪拌產(chǎn)生強(qiáng)烈的渦流，渦流的剪切作用使吸入浮選機(jī)的空氣被分散成大量尺寸不等的氣泡，這些氣泡的初速度大多在0～0.1 m/s之間[1]。氣泡在礦漿中運(yùn)動(dòng)時(shí)，礦粒與之碰撞并粘附于氣泡表面，礦化后的氣泡在浮選槽的靜止區(qū)域中上升以減少礦粒的脫附[2]，礦粒隨氣泡上升至礦漿表面后，被收集成為浮選產(chǎn)品。顯然，浮選生產(chǎn)的產(chǎn)量與產(chǎn)品質(zhì)量與氣泡運(yùn)動(dòng)密切相關(guān)，研究氣泡在礦漿中的運(yùn)動(dòng)規(guī)律對(duì)改進(jìn)生產(chǎn)過程、提高產(chǎn)品產(chǎn)量與質(zhì)量具有重要意義。

由于氣泡在液體中呈強(qiáng)非線性運(yùn)動(dòng)，運(yùn)動(dòng)過程中其界面變形較大，情況復(fù)雜多樣[3]，難以通過解析方法具體分析其運(yùn)動(dòng)變化規(guī)律，故大多采用數(shù)值模擬方法進(jìn)行分析研究。朱仁慶等[4]、武博等[5]和馬斗[6]模擬得出了氣泡在水中上升的運(yùn)動(dòng)速度隨時(shí)間變化的規(guī)律，與氣泡上升速度計(jì)算經(jīng)驗(yàn)公式所得結(jié)果較吻合；KRISHNA等[7]運(yùn)用VOF方法模擬了充滿水的二維矩形域內(nèi)單個(gè)氣泡的上升過程，得到的氣泡上升軌跡與實(shí)驗(yàn)觀測結(jié)果基本一致；WU等[8]通過實(shí)驗(yàn)觀察和理論分析方法獲得了氣泡的具體形成過程；LI等[9]也涉及了模擬氣泡在非牛頓流體中的運(yùn)動(dòng)情況。上述國內(nèi)外學(xué)者進(jìn)行的大量理論和實(shí)驗(yàn)研究均是針對(duì)氣泡在純液體中的運(yùn)動(dòng)，對(duì)氣泡在礦漿懸浮液中運(yùn)動(dòng)的模擬研究較少?LUO等[10-11]先后通過實(shí)驗(yàn)研究了在高壓液固懸浮體中氣泡的生成與運(yùn)動(dòng)過程，提出了一種基于球?qū)ΨQ假設(shè)的機(jī)理模型，對(duì)分析氣泡在懸浮液中的受力情況作了一定的貢獻(xiàn)；李彥鵬等[12]模擬了三維的液固懸浮液中氣泡的上升，該研究很好地證明了氣泡的上升會(huì)受到顆粒流的影響；李永強(qiáng)等[13]利用影像研究方法記錄了氣泡在旋流分選器中的運(yùn)動(dòng)現(xiàn)象，獲得了對(duì)氣泡運(yùn)動(dòng)規(guī)律的直觀認(rèn)識(shí)。在運(yùn)用VOF方法方面，XU等[14]和SEONG等[15]證明了VOF方法模擬的結(jié)果與實(shí)際情況的一致性較好。

上述文獻(xiàn)都使用不同研究方法成功地對(duì)氣泡的形狀、大小、軌跡等進(jìn)行了定性的研究分析，但極少用到具體的數(shù)值去表示氣泡變形的程度，研究浮選中氣泡最佳尺寸的相關(guān)報(bào)道也甚少。由于氣泡大小是浮選的重要表觀特征之一，氣泡過大或過小都直接影響著礦粒的回收率[16]，因此，本文作者在前面文獻(xiàn)工作的基礎(chǔ)及啟發(fā)下，將VOF方法運(yùn)用到簡化的懸浮液中氣泡運(yùn)動(dòng)的模擬上，在不同密度的靜止均勻礦漿中的分離區(qū)，對(duì)不同初速度和不同初始直徑下的單個(gè)氣泡運(yùn)動(dòng)進(jìn)行模擬，并使用氣泡的具體形變參數(shù)對(duì)氣泡在運(yùn)動(dòng)過程中形狀、大小、軌跡等量的變化進(jìn)行定量的分析，在不考慮氣泡運(yùn)載能力下，對(duì)浮選最佳升浮氣泡尺寸范圍進(jìn)行研究，將模擬與研究分析結(jié)果與文獻(xiàn)中已有的結(jié)果及實(shí)際生產(chǎn)情況進(jìn)行比較。

1 氣泡模擬的數(shù)值理論與方法

1.1 基本假設(shè)

本文主要研究氣泡在浮選機(jī)分離區(qū)的運(yùn)動(dòng)情況。為便于研究，作以下假設(shè)：

因浮選槽中分離區(qū)礦漿的紊流程度較小，氣泡進(jìn)入分離區(qū)之后運(yùn)動(dòng)速度以垂直向上方向的速度為主[17]，且單個(gè)氣泡在液-固懸浮流中的上升運(yùn)動(dòng)與其在高黏度流體中的運(yùn)動(dòng)類似，故將液固懸浮流視為靜止的不可壓縮的牛頓型均勻流體[18-19]；氣泡初始形狀為球形。

1.2 VOF方法

VOF方法[20-21]是一種通過追蹤流場網(wǎng)格中流體體積變化來捕捉兩相流運(yùn)動(dòng)界面的方法，其計(jì)算量小、精度高且容易實(shí)現(xiàn)。它定義每個(gè)網(wǎng)格中目標(biāo)流體的體積與網(wǎng)格體積之比為C，即C是一個(gè)流體體積函數(shù)，由C的大小可知兩相界面所處的位置：

體積函數(shù)C滿足如下控制方程：

式中：t為時(shí)間變量，s；u為流體運(yùn)動(dòng)速度矢量，m/s。

1.3 氣泡運(yùn)動(dòng)控制方程

因研究對(duì)象是一個(gè)等溫系統(tǒng)的多相流系統(tǒng)，且氣液兩相流為不可壓縮的牛頓流體，根據(jù)定義的VOF方法，可認(rèn)為這個(gè)有相界面分隔的兩相體系是單相流體系，故控制方程是對(duì)不可壓縮流體的質(zhì)量守恒方程和動(dòng)量守恒方程[22]。

1) 質(zhì)量守恒方程

2) 考慮表面張力的動(dòng)量方程

式中：ρ為流體密度；u為流體的速度矢量；p為壓強(qiáng)；Fs為表面張力源項(xiàng)；μ為動(dòng)力粘性系數(shù)。

1.4 連續(xù)表面張力模型

本文研究氣泡在礦漿中的運(yùn)動(dòng)考慮了表面張力的影響，在VOF模型中使用的表面張力，是通過散度原理表示為體積力，并采用了BRACKBILL等[23]提出的連續(xù)表面張力模型(Continuum surface force，CSF)將表面張力添加到動(dòng)量方程中的源項(xiàng)。

定義αq為第q相體積分?jǐn)?shù)，根據(jù)體積分?jǐn)?shù)分布，可得界面法向矢量為

可得表面曲率為

式中：單位法向矢量為

則表面張力表示為

式中：i和j為單元中的i相和j相。若一個(gè)單元只有兩相，則

式(8)即簡化為

1.5 氣泡運(yùn)動(dòng)過程中的形變參數(shù)

為更好的表達(dá)和說明氣泡變形的程度，本文定義了以下3種參數(shù)，對(duì)氣泡的形變從3個(gè)角度進(jìn)行數(shù)值上的定量分析，具體表示如圖1所示。

定義1：氣泡晃動(dòng)角θ，氣泡在運(yùn)動(dòng)中變形后，其左右兩頂點(diǎn)的連線與水平線的夾角定義為氣泡晃動(dòng)角，記為θ?；蝿?dòng)角的大小表明氣泡在垂直方向上的晃動(dòng)程度。在氣泡上升過程中，其最大晃動(dòng)角與最小晃動(dòng)角代數(shù)差的絕對(duì)值U越大，則氣泡晃動(dòng)程度 越大。

定義2：氣泡變形系數(shù)W，在運(yùn)動(dòng)中變形的氣泡，其長軸與短軸之比定義為氣泡變形系數(shù)，記為W。變形系數(shù)表明氣泡外形上左右、上下伸縮的形變程度，即扁平程度，在氣泡上升過程中，W越大，其變形程度越大。

定義3：氣泡偏移系數(shù)M，變形的氣泡中心與其初始?xì)馀葜行闹g的水平方向距離與初始?xì)馀莅霃街榷x為氣泡偏移系數(shù)，記為M。偏移系數(shù)表明了氣泡在上升過程中水平方向上的擺動(dòng)幅度。M越大，則表明氣泡上升運(yùn)動(dòng)軌跡的左右擺動(dòng)幅度越大。

圖1 氣泡形變參數(shù)定義示意圖Fig.1 Schematic diagrams of bubble deformation parameter definition: (a), (b) U=|θmax-θmin|; (c) W=b/a; (d) M=L/R

2 數(shù)值模擬結(jié)果與分析

2.1 初始邊界條件和參數(shù)

為消除壁面帶來的影響，當(dāng)氣泡直徑D為7～8 mm時(shí)，取計(jì)算區(qū)域?yàn)?.04 m×0.1 m，其他直徑取0.02 m×0.1 m的計(jì)算區(qū)域；采用四邊形網(wǎng)格，網(wǎng)格間距等于2.5×10-4m，時(shí)間步長等于1.0×10-4s，頂部邊界條件為壓力出口，其他為無滑移壁面，礦漿處于靜態(tài)，初始條件參數(shù)如表1所列。

表1 浮升氣泡運(yùn)動(dòng)數(shù)值模擬的初始條件參數(shù)Table1 Initial condition parameters for numerical simulation of flotation bubble motion

2.2 不同氣泡直徑對(duì)氣泡上升的影響

同一初速度、不同初始直徑的氣泡運(yùn)動(dòng)軌跡、形變、上升速率等均不同。4種不同初始尺寸氣泡在不同初速度時(shí)的運(yùn)動(dòng)軌跡如圖2所示。由圖2可見，氣泡上升過程中該在形變和擺動(dòng)，且變形和擺動(dòng)的幅度隨著初速度及初始直徑的不同而各不相同。氣泡在礦漿上升運(yùn)動(dòng)過程中因所受壓力、浮力、阻力等的不均衡而該在變形、晃動(dòng)等。

計(jì)算氣泡初始直徑為2和4 mm的晃動(dòng)角θ均在0.5°以下，可忽略。其他形變參數(shù)值如表2所列。由此可知，氣泡初始直徑越大，其最大變形系數(shù)越大，最大偏移系數(shù)越小，即氣泡的變形程度隨其初始直徑的增加而增大，其擺動(dòng)幅度隨氣泡初始直徑的增加而減少；同一尺寸的氣泡初速度越大，其最大變形系數(shù)和最大偏移系數(shù)也越大，即氣泡變形程度和擺動(dòng)幅度隨初速度的增加而增大；氣泡晃動(dòng)角隨著氣泡直徑和初速度的增加而增大，即氣泡初始直徑和初速度越大，其在垂直方向上的晃動(dòng)程度越大。

圖2 不同尺寸氣泡在不同初速度時(shí)的運(yùn)動(dòng)軌跡圖Fig.2 Trajectory diagram of different size bubble in different initial velocities: (a) v0=0.02 m/s; (b) v0=0.04 m/s; v0=0.06 m/s;v0=0.08 m/s (From left to right, the bubble diameter is 2, 4, 6 and 8 mm, respectively.)

表2 不同初速度下不同初始直徑氣泡運(yùn)動(dòng)中的形變參數(shù)最大值Table2 Maximum bubble deformation parameter value with different bubble initial diameter under different initial velocity during bubble motion

由表2中數(shù)據(jù)可看出，直徑為2 mm的氣泡變形程度最小，但水平方向上的偏移幅度最大；直徑為6和8 mm的氣泡形變程度大，垂直方向晃動(dòng)程度大。直徑為4 mm的氣泡形變過程較為平緩，與直徑為2 mm的氣泡相比，其水平方向上偏移幅度大幅度減?。挥纱丝梢?，直徑為4 mm的氣泡形變情況最為理想。

2.3 不同初速度對(duì)氣泡上升速度的影響

2.3.1 氣泡運(yùn)動(dòng)初始階段和穩(wěn)定階段的速度分布

不同直徑的氣泡在上升初始階段和上升穩(wěn)定后的速度矢量圖總體上類似，由圖3(a)可見，氣泡在上升初始階段，靜止的礦漿受到的擾動(dòng)使得氣泡周圍(液相)的速度場在左右兩側(cè)該在多個(gè)對(duì)稱的渦結(jié)構(gòu)，這些渦結(jié)構(gòu)增加了氣泡上浮的阻力和形變程度。由圖3(b)可見，隨著氣泡的上升，液流速度場趨于穩(wěn)定，氣泡兩側(cè)形成了兩個(gè)對(duì)稱的渦結(jié)構(gòu)，此時(shí)浮升力與液體壓力及其他阻力基本達(dá)到平衡狀態(tài)，氣泡在此時(shí)形變較為平緩，上升速度平穩(wěn)。

2.3.2 氣泡運(yùn)動(dòng)速率變化規(guī)律

在不同初速度下，氣泡在初始上升階段整體上是加速運(yùn)動(dòng)的，但受不同直徑氣泡和不同初速度的影響，其加速階段的具體加速過程有所不同，并呈現(xiàn)一定的規(guī)律，具體如圖4(a)～(d)所示。

圖4(a)所示為初速度為0.02 m/s時(shí)，4種直徑的氣泡上升初始階段的速度—時(shí)間關(guān)系圖。由圖4(a)可以看出，4種直徑氣泡均是一開始上升便做加速運(yùn)動(dòng)，且直徑為4 mm的氣泡運(yùn)動(dòng)速度一直居首，這表明當(dāng)初速度為0.02 m/s時(shí)，對(duì)液流的初始擾動(dòng)較小，氣泡所受浮升力均大于其所受的阻力，且直徑為4 mm的氣泡加速得最快。

圖4(b)所示為初速度0.04 m/s時(shí)，4種直徑的氣泡開始上升到趨于穩(wěn)定階段的速度—時(shí)間關(guān)系圖。由圖4(b)可以看出，直徑為4 mm的氣泡運(yùn)動(dòng)規(guī)律是氣泡一直做加速運(yùn)動(dòng)直至速度趨于穩(wěn)定，且直徑為4 mm的氣泡運(yùn)動(dòng)速度大于其他直徑氣泡的運(yùn)動(dòng)速度；直徑為2、6和8 mm的氣泡運(yùn)動(dòng)規(guī)律是先減速，然后加速直至速度趨于平穩(wěn)。

圖4(c)和(d)所示的氣泡運(yùn)動(dòng)規(guī)律與圖4(b)的相同，這表明當(dāng)該在由初速度為0.04、0.06和0.08 m/s引起的上浮阻力的情況下，直徑為4 mm的氣泡均能夠克服擾動(dòng)及其他阻力先進(jìn)行加速運(yùn)動(dòng)；而直徑為2、6和8 mm的氣泡在較大初速度作用下，其浮升能力小于上浮阻力而首先作減速運(yùn)動(dòng)，之后隨著上升運(yùn)動(dòng)阻力和液體壓力的降低，運(yùn)動(dòng)速度逐漸增大直至穩(wěn)定。

綜合4種初速度下的氣泡運(yùn)動(dòng)規(guī)律圖可見，初始直徑為4 mm的氣泡抗液流擾動(dòng)能力及初始階段上浮速度較其他直徑氣泡的要大。

2.3.3 氣泡該程運(yùn)動(dòng)時(shí)間

不同初速度情況下不同初始尺寸氣泡該程運(yùn)動(dòng)的時(shí)間各不相同。由圖5可見，同一初始直徑氣泡到達(dá)液面的時(shí)間總體上隨著初速度的增加而縮短，4種氣泡到達(dá)液面的先后順序按直徑排列為4、6、8和2 mm；直徑為6 mm的氣泡的平均速度略大于直徑為8 mm的氣泡的平均速度。v0= 0.02 m/s時(shí)，初始直徑為6和8 mm的氣泡由于在起始運(yùn)動(dòng)狀態(tài)下阻力較小，初始階段沒有減速過程，該程運(yùn)動(dòng)時(shí)間長于其他初速度條件下的該程運(yùn)動(dòng)時(shí)間；v0= 0.04 m/s時(shí)，初始上升階段的阻力較v0=0.02 m/s時(shí)的阻力大，初始直徑為6和8 mm的氣泡先減速了一段時(shí)間，因此，在該計(jì)算區(qū)域內(nèi)，其上升時(shí)間長于初速度為0.02 m/s時(shí)的上升時(shí)間；v0為0.06和0.08 m/s時(shí)，雖然其上升初始階段也該在較大阻力和一段時(shí)間的減速運(yùn)動(dòng)，但其減速完成并開始加速時(shí)的速度及其減速過程運(yùn)動(dòng)的路程較初速度為0.04 m/s時(shí)的大，因此，其運(yùn)動(dòng)時(shí)間較初速度為0.04 m/s時(shí)的要短。直徑為2 mm的氣泡由于浮升力最小，其平均上升速度最慢而最后到達(dá)液面。直徑為4 mm的氣泡在所有初速度情況下的該程運(yùn)動(dòng)時(shí)間最短，且隨著初速度的增大，其上升至液面時(shí)間縮短，說明直徑為4 mm的氣泡上升平均速度最大，其上浮動(dòng)力與上浮阻力之比大于其他直徑氣泡的上浮動(dòng)力與上浮阻力之比，浮升能力最強(qiáng)。

圖3 氣泡上升速度矢量圖Fig.3 Vector diagram of bubble rise velocity: (a) Initial stage; (b) Stabilization stage

圖4 不同直徑氣泡在不同初速度下0～0.04 s運(yùn)動(dòng)時(shí)間段的變化規(guī)律Fig.4 Change rule of different diameter bubble motions during 0-0.04 s at different initial velocities: (a) v0=0.02 m/s;(b) v0=0.04 m/s; (c) v0=0.06 m/s; (d) v0=0.08 m/s

圖5 不同初速度下不同直徑氣泡到達(dá)液面的時(shí)間Fig.5 Time of different diameter bubble reaching to liquid level at different initial velocities

氣泡形變和擴(kuò)張程度及擺動(dòng)幅度越大，其表面粘附的礦粒越易脫落，氣泡越易破滅；氣泡初始直徑過小，則其表面能粘附的礦粒尺寸及數(shù)量過小，氣泡對(duì)礦粒的攜載和輸運(yùn)能力隨之變?nèi)?，兩者均?huì)導(dǎo)致浮選生產(chǎn)率降低，礦物回收率減少。綜合 2.2～2.3節(jié)模擬結(jié)果及分析所述可知，初始直徑為4 mm的氣泡形變過程較為平緩，水平方向上偏移幅度較小，其浮升力與上升阻力之比較其他初始直徑氣泡的要大；在所有初速度情況下，其運(yùn)動(dòng)時(shí)間最短，即其上升平均速度最大，浮升能力最強(qiáng)，對(duì)礦粒的承載能力最強(qiáng)。因此，浮升氣泡的尺寸以中等大小(直徑為4 mm)為最佳。

2.3.4 氣泡上升速度模擬計(jì)算與理論公式計(jì)算比較

對(duì)于單個(gè)氣泡上升速度公式的推導(dǎo)有很多著名的公式，Stokes及Mendelson公式是其中之一，如式(12)所示，表示的是氣泡上升平穩(wěn)速度(u)與氣泡直徑的關(guān)系。

式中：ρf為流體密度，kg/m3；ρg為氣體密度，kg/m3；CD為阻力系數(shù)；db為氣泡的當(dāng)量直徑，m；g為重力加速度，m/s2。

氣泡的上升過程可分為層流區(qū)、過渡區(qū)和紊流區(qū)3種情況。

1) 在層流區(qū)(Re＜1)時(shí)，則氣泡上升速度公式為

式中：σ為表面張力系數(shù)，N/m。

2) 在過渡區(qū)(1＜Re＜300)時(shí)，氣泡上升速度公式為

3) 在紊流區(qū)(300＜Re＜4 000)時(shí)，則氣泡上升速度公式為

由式(13)、(14)和(15)可看出，當(dāng)氣泡上升處于層流區(qū)和過渡區(qū)時(shí)，氣泡上升平穩(wěn)速度與氣泡直徑大小成正比；當(dāng)氣泡處于紊流區(qū)時(shí)，氣泡上浮平穩(wěn)速度與氣泡直徑成反比。由于氣泡在上浮過程中多處于過渡區(qū)和紊流區(qū)，因此，氣泡在流體中的上升速度并不單純地隨著氣泡直徑的增大而增加。氣泡的上升速度在運(yùn)動(dòng)了0.04～0.05 s之后就會(huì)達(dá)到一個(gè)平穩(wěn)值，如圖5所示，直徑為4 mm的氣泡上升平均速度最大，表明其上升平穩(wěn)后的速度亦最大，與理論推導(dǎo)公式分析結(jié)果一致。在實(shí)際的氣泡上升過程中，氣泡很難達(dá)到力的平衡，且隨著氣泡的上升，其體積逐漸變大，造成阻力和浮升力出現(xiàn)非線性變化，此時(shí)氣泡受力關(guān)系變得極其復(fù)雜。

3 不同礦漿密度對(duì)氣泡運(yùn)動(dòng)的影響

3.1 不同密度礦漿中氣泡的形變程度

本研究模擬3種不同礦漿密度下氣泡上升的運(yùn)動(dòng)情況，其他條件不變。作為對(duì)比，增加了3種直徑氣泡的模擬實(shí)驗(yàn)，具體初始模擬參數(shù)如表3所列，模擬及計(jì)算結(jié)果如圖6和表4所示。由圖6可見，不同礦漿密度中，氣泡上升過程中的形變情況及上升軌跡有所不同；由表4可知，氣泡變形系數(shù)與礦漿密度成正比，氣泡偏移系數(shù)與礦漿密度成反比，即氣泡變形程度隨礦漿密度的減少而減少，氣泡上升軌跡水平方向上擺動(dòng)幅度隨礦漿密度的減少而增大；大氣泡的偏角隨礦漿密度的增大而減少，即礦漿密度越大,氣泡在垂直方向上的擺動(dòng)越小。

3.2 不同密度礦漿中氣泡運(yùn)動(dòng)平均速度規(guī)律

不同初始直徑氣泡在不同密度礦漿中到達(dá)液面的時(shí)間各不相同，如圖7所示。為比較新增初始直徑氣泡與原尺寸氣泡的模擬結(jié)果，計(jì)算礦漿密度為1 175 kg/m3、氣泡初始直徑為4、5、6、7 mm的各種形變參數(shù)及礦漿密度為1 600 kg/m3、氣泡初始直徑為3、4、5、6 mm的各種形變參數(shù)，計(jì)算結(jié)果如表5和6所列。

由圖7可見，礦漿密度越大，氣泡到達(dá)液面所需時(shí)間越長，表明氣泡隨著礦漿密度的增大受到的黏滯阻力越大，氣泡上升的平均速度越小。由表5可知，在礦漿密度為1 175 kg/m3時(shí)，直徑為5 mm的氣泡的偏移系數(shù)最大值(0.122)小于直徑為4 mm的氣泡的偏移系數(shù)最大值(0.251)。同時(shí)，由圖7可見，直徑為5 mm的氣泡到達(dá)液面所需的時(shí)間比直徑為4 mm的氣泡所需的時(shí)間要少，即其上升的平均速度更大，所以在礦漿密度為1 175 kg/m3時(shí)，最佳浮升氣泡尺寸在5 mm左右。由表6可見，在礦漿密度為1 600 kg/m3時(shí)，直徑為3 mm的氣泡的變形系數(shù)最大值(1.475)小于直徑為4、5 mm的氣泡的變形系數(shù)最大值，且直徑為3 mm的氣泡的晃動(dòng)角(0.18)比直徑為4 mm的氣泡的晃動(dòng)角(0.52)小，上升的平均速度略小于直徑為4 mm的氣泡的平均速度。綜合各方面情況可見，在礦漿密度為1 600 kg/m3時(shí)，最佳浮升氣泡趨于初始直徑為3 mm左右的氣泡。

表3 在不同密度礦漿中模擬的初始條件參數(shù)Table3 Initial condition parameters of simulation in different ore pulp densities

圖6 不同礦漿密度中氣泡的運(yùn)動(dòng)規(guī)律Fig.6 Bubble motion law in different ore pulp density: (a) 1 175 kg/m3; (b) 1 400 kg/m3; (c) 1 600 kg/m3; (d) From left to right:D=3 mm, ρ=1 600 kg/m3; D=5 mm, ρ=1175 kg/m3; D=5 mm, ρ=1 600 kg/m3; D=7 mm, ρ=1175 kg/m3

表4 不同密度下不同初始直徑氣泡運(yùn)動(dòng)中的形變參數(shù)最大值Table4 Maximum bubble deformation parameter value with different bubble initial diameters in different densities during bubble motion

圖7 不同礦漿密度下不同初始直徑氣泡到達(dá)液面的時(shí)間Fig.7 Time of different diameter bubble reaching to liquid level under different pulp densities

綜上模擬結(jié)果及分析可知：礦漿密度為1 175 kg/m3時(shí)，最佳浮升氣泡尺寸在5 mm左右；礦漿密度為1 400 kg/m3時(shí)，最佳浮升氣泡尺寸在4 mm左右；礦漿密度為1 600 kg/m3時(shí)，最佳浮升氣泡尺寸在3 mm左右。因此，在不同密度的礦漿中，最佳浮升氣泡初始尺寸各不相同。

表5 礦漿密度1 175 kg/m3下氣泡形變最大參數(shù)值Table5 Maximum bubble deformation parameter values in density of 1 175 kg/m3

表6 礦漿密度1 600 kg/m3下氣泡形變參數(shù)值Table6 Maximum bubble deformation parameter values in density of 1 600 kg/m3

在實(shí)際生產(chǎn)中，礦漿含量在25%～50%之間，對(duì)應(yīng)的礦漿密度為1 175～1 400 kg/m3，因此，通常情況下，中等尺寸(直徑為3～4 mm)的浮升氣泡可獲得較好的浮選生產(chǎn)效果。

4 結(jié)論

1) 運(yùn)用VOF方法對(duì)不同初始直徑、不同初速度條件下的氣泡在3種不同密度的礦漿中的運(yùn)動(dòng)情況進(jìn)行了模擬實(shí)驗(yàn)，得到了氣泡的運(yùn)動(dòng)軌跡和速度分布圖，并計(jì)算了氣泡形變最大參數(shù)值。

2) 在礦漿密度為1 400 kg/m3時(shí)，直徑為4 mm的氣泡相比其他氣泡形變平緩，上升速度最大，浮升能力最強(qiáng)，最有利于礦粒的回收；礦漿密度有所偏大或偏小時(shí)，浮升氣泡的最佳尺寸仍為3～5 mm。這證明了氣泡尺寸過大或過小時(shí)，其運(yùn)動(dòng)形變情況均不利于礦粒的回收，該模擬結(jié)果與實(shí)際生產(chǎn)情況高度吻合。

3) 對(duì)模擬結(jié)果的分析表明：當(dāng)氣泡初速度v0=0.02 m/s時(shí)，氣泡在上升初始階段受到液流的擾動(dòng)影響不大；當(dāng)氣泡初速度v0＞0.02 m/s時(shí)，氣泡受到的初始上浮阻力較大。

4) 實(shí)驗(yàn)結(jié)果證明，與實(shí)際生產(chǎn)中的礦漿濃度對(duì)應(yīng)的最佳浮升氣泡尺寸為3～4 mm。

[1]選礦機(jī)械設(shè)備實(shí)用手冊(cè).第六篇.浮選機(jī)械設(shè)備[J].北京: 機(jī)械工業(yè)出版社, 1994: 617-621.Practical manual of mineral processing mechanical equipment.Sixth part.Flotation mechanical equipment[J].Beijing: Chinese Machine Press, 1994: 617-621.

[2]韓 偉.浮選機(jī)內(nèi)多相流動(dòng)特性及浮選動(dòng)力學(xué)性能的數(shù)值研究[D].蘭州：蘭州理工大學(xué), 2009: 9-12.HAN Wei.The numerical research of multiphase flow characteristics and flotation dynamics performance in flotation machine[D].Lanzhou: Lanzhou University of Technology, 2009:9-12.

[3]PETER L, FRANZ P, NICOLAS F, PETER E.Gas bubble in simulation and experiment[J].Journal of Colloid and Interface Science, 2011, 354: 364-372.

[4]朱仁慶, 李晏丞, 倪永燕, 侯搖玲.氣泡在水中上升運(yùn)動(dòng)的數(shù)值模擬[J].江蘇科技大學(xué)學(xué)報(bào): 自然科學(xué)版, 2010, 24(5): 2-3.ZHU Ren-qing, LI Yan-chen, NI Yong-yan, HOU Yao-ling.Numerical simulation of bubble rising in the water[J].Journal of Jiangsu University of Science and Technology: Natural Science Edition, 2010, 24(5): 2-3.

[5]武 博, 郝宗睿, 陳 濤, 吳大轉(zhuǎn).水下氣泡運(yùn)動(dòng)的數(shù)值模擬[J].科技論文在線, 2010, 5(8): 2-4.WU Bo, HAO Zong-rui, CHEN Tao, WU Da-zhuan.Numerical simulation on motion of bubble under water[J].Science Paper Online, 2010, 5(8): 2-4.

[6]馬 斗.單孔氣泡動(dòng)力學(xué)行為的VOF數(shù)值模擬[D].天津: 天津大學(xué), 2009: 37-43.MA Dou.Numerical simulation of single orifice gas bubble motion behavior with volume of fluid method [D].Tianjin:Tianjin University, 2009: 37-43.

[7]KRISHNA R, van BATEN J M.Rise Characteristics of Gas Bubbles in a 2D Rectangular column: VOF simulation vs experiments[J].International Communications in Heat Mass Transfer, 1999, 26(7): 965-970.

[8]WU Rui-zhi, SHU Da, SUN Bao-de, WANG Jun, LU Yan-ling.Observation and theoretic analysis of gas-bubble formation and growth in water-model[J].Transactions of Nonferrous Metals Society of China, 2005, 15(5): 1130-1131.

[9]LI Zhang, CHAO Yang, CAO Zai-sha.Numerical simulation of a bubble rising in shear-thinning fluids[J].Journal of Non-Newtonian Fluid Mechanics, 2010, 165(11/12): 555-561.

[10]LUO X, YANG G Q, LEE D J.Single bubble formation in high pressure liquid-solid suspensions[J].Powder Technology, 1998,100(2/3): 103-108.

[11]YANG G Q, LUO X, LAU R.Bubble formation in high pressure liquid-solid suspensions with plenum pressure fluctuation[J].AICHE J, 2000, 46: 2163-2168.

[12]李彥鵬, 張乾隆.液固懸浮液中氣泡生成的數(shù)值模擬[J].水動(dòng)力學(xué)研究與進(jìn)展, 2007, 22(4): 457-462.LI Yan-peng, ZHANG Qian-long.Numerical simulation of bubble formation in liquid-solid suspensions[J].Journal of Hydrodynamics, 2007, 22(4): 457-462.

[13]李永強(qiáng), 李長根.旋流閃速浮選法的研究Ⅱ: 旋流分選器中氣泡的運(yùn)動(dòng)規(guī)律[J].有色金屬, 2000, 52(2): 22-23.LI Yong-qiang, LI Chang-gen.Study on cyclone flash flotationⅡ: Study on movement pattern of air bubbles in separation cyclone for cyclone flash flotation[J].Nonferrous Metals, 2000, 52(2): 22-23.

[14]XU Ling-jun, CHEN Gang, SHAO Jian-bin.Numerical simulation of bubble behavior using VOF method[J].Electric Technology and Civil Engineering, 2011, 22(24): 2408-2410.

[15]SEONG S J, SEONG J K, GOON C P.CFD simulation of condensing vapor bubble using VOF method[J].World Academy of Science, Engineering and Technology, 2009, 60(36): 1-6.

[16]龔明光.泡沫浮選[J].北京: 冶金工業(yè)出版社, 2008: 108-110.GONG Ming-guang.Froth flotation[J].Beijing: Metallurgical Industry Press, 2008: 108-110.

[17]曾克文.浮選槽內(nèi)礦漿紊流強(qiáng)度對(duì)浮選影響的理論及應(yīng)用研究[D].長沙:中南大學(xué), 2001: 51-53.ZENG Ke-wen.Theory and application study on the effect of pulp turbulence degree in flotation cell on flotation[D].Changsha: Central South University, 2001: 51-53.

[18]YANG G Q, BING D U, FAN L S, Bubble formation and dynamics in gas-liquid-solid fluidization-A review[J].Journal of Chemical and Biomolecular Engineering, 2007, 62(1/2): 2-5.

[19]HUA J S, STENE J F, LIN P.Numerical simulation of 3D bubbles rising in viscous liquids using a front tracking method[J].Journal of Computational Physics, 2008, 227(6): 3358-3382.

[20]ZHANG Yu-jie, LIU Ming-yan, XU Yong-gui, TANG Can.Three-dimensional volume of fluid simulations on bubble formation and dynamics in bubble columns[J].Chemical Engineering Science, 2012, 73(7): 55-78.

[21]MA Dou, LIU Ming-yan, ZU Yong-gui, TANG Can.Two-dimensional volume of fluid simulation studies on single bubble formation and dynamics in bubble columns[J].Chemical Engineering Science, 2012, 72(16): 61-75.

[22]徐玲君, 陳 剛, 邵建斌, 薛 陽.單個(gè)氣泡靜水中上升特性的數(shù)值模擬[J].沈陽農(nóng)業(yè)大學(xué)學(xué)報(bào), 2012, 43(3): 357-361.XU Ling-jun, CHEN Gang, SHAO Jian-bin, XUE Yang.Numerical simulation of motion behavior of single bubble rising in still water[J].Journal of Shenyang Agricultural University,2012, 43(3): 357-361.

[23]BRACKBILL J U, KOTHE D B, ZEMACH C.A continuum method for modeling surface tension[J].Journal of Computational Physics, 1992, 100: 335-354.