含塵氣泡形變及相互作用對(duì)顆粒脫除率的影響

梁 纖， 潘偉童， 陳雪莉

（ 華東理工大學(xué)國(guó)家能源煤氣化技術(shù)研發(fā)中心，上海煤氣化工程技術(shù)研究中心，上海 200237）

煤氣化作為現(xiàn)代煤化工的龍頭，氣化界區(qū)的合成氣中細(xì)灰顆粒含量需嚴(yán)格控制在1 mg/Nm3以下[1]，因此煤氣化粗合成氣中細(xì)灰的高效脫除尤為重要。在已工業(yè)化應(yīng)用的煤氣化技術(shù)中，合成氣中細(xì)灰的脫除多采用濕洗分離方法，而鼓泡接觸型濕式分離又是其中的首選[2]。含塵氣泡是鼓泡接觸型濕式分離體系中的單元集塵體，對(duì)含塵氣泡的運(yùn)行特性進(jìn)行深入研究，將為細(xì)顆粒的高效脫除和鼓泡洗滌器的優(yōu)化設(shè)計(jì)提供重要的理論支撐。

在前人的研究中，多將集塵氣泡分為兩部分討論：氣泡的單獨(dú)上升運(yùn)動(dòng)與理想化流場(chǎng)下氣泡攜帶顆粒的上升運(yùn)動(dòng)。1976 年，Grace 等[3]給出了單個(gè)不含塵氣泡上升的形狀譜圖，Bhaga 等[4]對(duì)氣泡圖譜進(jìn)行了完善。之后又有許多學(xué)者對(duì)單個(gè)氣泡進(jìn)行了研究，包括上升速度、界面變化、添加表面活性劑等[5-8]。對(duì)于同步上升的兩氣泡，其受力和運(yùn)動(dòng)情況與單氣泡不同，其中研究最多的是并行雙氣泡的上升過(guò)程[9-11]，對(duì)其他組合方式雙氣泡的共同上升過(guò)程也有研究[12-14]。

在系統(tǒng)研究鼓泡體系中氣液兩相相互作用機(jī)制的基礎(chǔ)上，學(xué)者們開始關(guān)注氣液固三相間相互作用機(jī)制。顆粒在氣泡內(nèi)部的運(yùn)動(dòng)傳輸主要受不同的傳輸機(jī)制控制，如慣性擴(kuò)散、布朗擴(kuò)散、對(duì)流等機(jī)制[15-17]。1964 年Fuchs 等[18]通過(guò)實(shí)驗(yàn)研究首先預(yù)測(cè)了氣泡在上升過(guò)程中會(huì)出現(xiàn)內(nèi)部循環(huán)，該循環(huán)將大大增強(qiáng)粒子的擴(kuò)散和慣性作用，并據(jù)此提出了球形上升氣泡內(nèi)顆粒脫除的數(shù)學(xué)模型。Li 等[19]運(yùn)用CFD 計(jì)算模擬軟件對(duì)氣液固三相流態(tài)化系統(tǒng)建立了模型，并模擬了液固兩相床中單氣泡上升過(guò)程的顆粒夾帶和尾跡脫落現(xiàn)象，結(jié)果表明其所建立的模型與文獻(xiàn)報(bào)道結(jié)果吻合良好。Xu 等[20]運(yùn)用電聲技術(shù)考察了不同粒徑和潤(rùn)濕性的石英及玻璃微珠顆粒在振蕩氣泡中的分離過(guò)程，結(jié)果發(fā)現(xiàn)氣泡-顆粒聚集體在高黏度液體中的穩(wěn)定性與顆粒的動(dòng)態(tài)接觸角有關(guān)。Pan 等[21]通過(guò)流體體積-離散相模型（VOF-DPM）耦合法模擬研究了含塵氣泡在靜止液體中上升時(shí)的形變特征，分析了氣泡形變對(duì)顆粒脫除的影響，結(jié)果表明氣泡形變能夠增強(qiáng)微細(xì)顆粒的脫除率。

綜上分析可知，關(guān)于含塵氣泡上升過(guò)程的動(dòng)態(tài)變化與顆粒脫除機(jī)制間的內(nèi)在關(guān)系尚不明晰，且大部分模擬研究?jī)H是基于二維物理模型，無(wú)法反映含塵氣泡上升過(guò)程動(dòng)態(tài)變化的全貌。基于此，本文將通過(guò)建立三維物理模型模擬研究靜止液體中不同含塵氣泡上升過(guò)程對(duì)顆粒脫除的影響。

1 數(shù)值模擬

1.1 物理模型

基于二維物理模型的局限性，本文將采用三維物理模型對(duì)含塵氣泡上升過(guò)程進(jìn)行模擬研究。三維物理模型如圖1 所示，圖中顯示的為單氣泡，雙氣泡時(shí)的計(jì)算域模型與此相同。相；α=0 時(shí)，表示該區(qū)域沒(méi)有該相；0＜α＜1 時(shí)為混合物[23]。不可壓縮流體的流動(dòng)控制方程如下：

連續(xù)性方程：

動(dòng)量方程：

其中u,ρ 分別為計(jì)算相的速度和密度；Fsf是表面張力；Fpg是來(lái)自顆粒的作用于流體相的動(dòng)量交換力，作用于含有顆粒的晶胞。在VOF 模型中通過(guò)求解整個(gè)區(qū)域內(nèi)的單一的動(dòng)量方程，得到由各相共享的速度場(chǎng)，而動(dòng)量方程取決于通過(guò)所有相的容積比率。

體積分?jǐn)?shù)方程：

計(jì)算域中充滿了水(溫度25 ℃，密度9.982×102kg/m3，黏度10-3Pa·s，表面張力系數(shù)0.072 8 N/m)。初始時(shí)在計(jì)算區(qū)域底部中心放置一個(gè)含塵氣泡(氣體密度1.225 kg/m3，氣體黏度1.789×10-5Pa·s；含塵顆粒粒徑分別為1 μm 和2.5 μm，密度2×103kg/m3)。模擬計(jì)算時(shí)根據(jù)氣泡尺寸大小調(diào)整計(jì)算域的寬度，計(jì)算域的高度始終保持在60 mm，網(wǎng)格劃分使每個(gè)氣泡內(nèi)至少包含10 000 個(gè)單元以確保計(jì)算的精確性。

模擬計(jì)算基于以下幾點(diǎn)假定：

(1) 氣體被認(rèn)為是不可壓縮的。

(2) 顆粒的脫除通過(guò)水膜擦洗方式實(shí)現(xiàn)，當(dāng)顆粒到達(dá)液體體積分?jǐn)?shù)0.5 的區(qū)域即認(rèn)為被脫除(在已有文獻(xiàn)報(bào)道中認(rèn)為空泡體積分?jǐn)?shù)大于0.5 的部分是氣泡[22])。因此本文忽略了液體-顆粒之間的相互作用和顆粒通過(guò)液膜時(shí)的界面力。

(3) DPM 模型將微米級(jí)顆粒視為運(yùn)動(dòng)質(zhì)點(diǎn)，忽略顆粒間的相互作用及其對(duì)氣相體積的影響。

1.2 數(shù)值模型

1.2.1 流體相模型 模擬中氣相與液相均視為連續(xù)相，通過(guò)歐拉法計(jì)算其連續(xù)性方程與動(dòng)量方程，數(shù)值模型采用VOF 模型。VOF 模型是一種表面跟蹤技術(shù)，用于處理沒(méi)有相互穿插的多相流問(wèn)題，通過(guò)求解一組動(dòng)量方程并在整個(gè)計(jì)算域內(nèi)跟蹤每一相的體積分?jǐn)?shù)(α)來(lái)區(qū)分兩相，當(dāng)α=1 時(shí)，表示該區(qū)域充滿該

其中，下標(biāo)l 表示每個(gè)控制體積中的液相，由于所有相的體積分?jǐn)?shù)之和應(yīng)為1，因此可以根據(jù)式（4）計(jì)算氣相的體積分?jǐn)?shù)( αg)：

表面張力是一種作用在表面上的力，它平衡了內(nèi)部放射狀的分子引力和跨過(guò)表面的放射狀的外部壓力梯度。表面張力對(duì)處于水中的氣泡的變形影響顯著。另外，F(xiàn)luent 中表面張力模型采用的是Brackbill等[24]提出的連續(xù)表面力模型。

1.2.2 固體相模型 由于固體顆粒的體積分?jǐn)?shù)較低，將其視為離散相，采用拉格朗日方法，選用DPM 模型，并假定不考慮顆粒-顆粒之間的相互作用以及顆粒體積分?jǐn)?shù)對(duì)連續(xù)相的影響。Fluent 中通過(guò)積分拉氏坐標(biāo)系下的顆粒作用力微分方程來(lái)求解離散相顆粒的軌道。顆粒的作用力平衡方程(顆粒慣性＝作用在顆粒上的各種力)在笛卡爾坐標(biāo)系下的形式為（x方向）：

其中μl為液體動(dòng)力黏度，dp為顆粒粒徑，Cc為Stokes'曳力公式的Cunningham 修正(考慮稀薄氣體力學(xué)的顆粒壁面速度滑移的修正)，其計(jì)算公式為：

其中λ為氣體分子平均自由程：

式中T為溫度；dm為空氣氣體分子的直徑，dm=3.5×10-10m；P是氣相壓力。由式（8）計(jì)算得平均自由程λ=6.9×10-8m。

在整個(gè)計(jì)算過(guò)程中，顆粒脫除率η可定義為：

式中Np是該時(shí)刻氣泡內(nèi)剩余顆粒數(shù)，Np0是初始時(shí)刻氣泡內(nèi)總顆粒數(shù)。

1.2.3 氣泡上升形態(tài)驗(yàn)證 圖2 所示為單氣泡上升形態(tài)實(shí)驗(yàn)裝置。透明有機(jī)玻璃水箱尺寸為200 mm×200 mm×600 mm，有機(jī)玻璃板厚度為1 cm，箱內(nèi)裝有適量的純水。實(shí)驗(yàn)裝置中的氣泡發(fā)生部分主要由空氣注射器、測(cè)量管道、可旋轉(zhuǎn)調(diào)節(jié)的杯罩和單向閥組成。在實(shí)驗(yàn)裝置中分別進(jìn)行了尺寸為6、8、10、12 mm氣泡的上升實(shí)驗(yàn)。

圖2 單氣泡上升實(shí)驗(yàn)裝置Fig. 2 Device of single bubble rise experiment

圖3(a)所示為實(shí)驗(yàn)拍攝的氣泡形態(tài)變化，氣泡尺寸從左至右依次為6、8、10、12 mm；圖3(b)所示為模擬條件下含塵單氣泡在上升過(guò)程中的形態(tài)變化結(jié)果，對(duì)比圖3(a)可知，與實(shí)驗(yàn)結(jié)果基本一致。

圖3 氣泡形變的實(shí)驗(yàn)結(jié)果（a）和模擬結(jié)果（b）比較Fig. 3 Comparison of bubble deformation experimental results (a)with simulation results (b)

2 結(jié)果與討論

2.1 含塵氣泡上升過(guò)程

2.1.1 含塵單氣泡上升過(guò)程 模擬計(jì)算了不同尺寸(6、8、10、12、14 mm)含塵氣泡的上升過(guò)程，氣泡中顆粒粒徑均一，分別為1.0 μm 和2.5 μm，結(jié)果如圖4所示。圖4 中顆粒粒徑為1.0 μm，紅色部分為顆粒(非顆粒實(shí)際大小，為了便于觀察將顆粒放大處理)。

由圖4 可以看出，對(duì)于給定的氣泡尺寸范圍，氣泡在上升過(guò)程中均發(fā)生形變，且形變總是從氣泡下表面向內(nèi)凹陷開始，最終都發(fā)展成扁橢球的形狀，但隨著氣泡尺寸增大，氣泡的形狀變化越來(lái)越穩(wěn)定。尺寸較小(6～8 mm)的氣泡呈現(xiàn)出復(fù)雜的變形和不規(guī)則的振蕩，形變強(qiáng)度較高；當(dāng)氣泡尺寸大于12 mm時(shí)，形狀變化被抑制，當(dāng)下部凹陷變?yōu)椤邦^盔”狀后，氣泡由于倒流渦旋而使其尾部變平，隨后成為帶有凹入底部的球帽并上升。

圖4 不同尺寸含塵單氣泡上升過(guò)程Fig. 4 Rising process of dust-containing single bubbles of different sizes

2.1.2 并行含塵雙氣泡上升過(guò)程 單氣泡及其上升過(guò)程動(dòng)力學(xué)的研究已較多，而在雙氣泡共同上升過(guò)程中，由于氣泡間的相互作用，其形變必然不同于單氣泡，從而對(duì)氣泡內(nèi)顆粒的脫除產(chǎn)生影響，因此有必要對(duì)雙氣泡的共同上升過(guò)程進(jìn)行模擬。雙氣泡共同上升的模式較多，其中等尺寸并行上升氣泡是較典型且有代表性的一種。本文模擬計(jì)算了不同尺寸(6、8、10、12、14 mm)含有1.0 μm 顆粒的并行雙氣泡上升過(guò)程，單個(gè)氣泡中顆粒數(shù)量與前文相同，雙氣泡間初始間距均為5 mm。圖5 所示為不同尺寸含塵雙氣泡上升過(guò)程示意圖，圖中紅色部分為顆粒(非顆粒實(shí)際大小，為了便于觀察將顆粒放大處理)。

圖5 不同尺寸并行含塵雙氣泡上升過(guò)程Fig. 5 Rising process of parallel dust-containing double bubbles of different sizes

計(jì)算域高度一定，在固定初始間距為5 mm 情況下，并行雙氣泡在上升過(guò)程中發(fā)生了形變但均未發(fā)生聚并；尺寸較小的氣泡形變與單氣泡上升時(shí)形變相似，氣泡首先變化為扁橢球形，呈搖擺狀上升，氣泡對(duì)之間呈現(xiàn)出靠近-分離-靠近-分離的趨勢(shì)。分析認(rèn)為，氣泡一旦開始上升，周圍流體將產(chǎn)生渦旋，伴隨渦旋脫落將形成尾流區(qū)，尾流內(nèi)流體上升速度較大，隨著尾流的增長(zhǎng)，氣泡的搖擺運(yùn)動(dòng)使得兩個(gè)氣泡可以正向運(yùn)動(dòng)也可以反向運(yùn)動(dòng)。并行上升的大尺寸氣泡因氣泡間相互作用較強(qiáng)，使其氣泡形變與單氣泡上升時(shí)的形變顯著不同。上升初期，兩氣泡間交匯處形成剪切變稀流體區(qū)域，該黏度特征推動(dòng)氣泡在水平方向相互靠近，而靠近后的大尺寸氣泡對(duì)間的渦旋擠壓作用增強(qiáng)，氣液界面向兩側(cè)凹陷，氣泡對(duì)被推動(dòng)分離，呈現(xiàn)V 型對(duì)稱式分布；隨間距的增加，內(nèi)側(cè)渦旋作用減弱，氣泡對(duì)再次靠近。

2.2 不同含塵氣泡上升過(guò)程中顆粒的脫除率

2.2.1 含塵單氣泡上升過(guò)程中顆粒的脫除率 圖6 所示為含塵單氣泡上升過(guò)程中顆粒脫除率隨氣泡尺寸的變化，氣泡中顆粒粒徑分別為1.0 μm 和2.5 μm。

由圖6 可知，當(dāng)顆粒粒徑為1.0 μm 時(shí)，隨著氣泡尺寸增加，顆粒的脫除率單調(diào)下降；當(dāng)氣泡中顆粒粒徑為2.5 μm 時(shí)，不同氣泡上升過(guò)程中顆粒的脫除率幾乎都高于粒徑1.0 μm 顆粒的脫除率，且隨氣泡尺寸增大顆粒脫除率曲線不再呈單調(diào)變化，而是在氣泡尺寸為8 mm 左右時(shí)出現(xiàn)峰值(97.5%)。綜上說(shuō)明，小尺寸氣泡中的顆粒更易脫除；對(duì)不同粒徑的顆粒，存在使其高效脫除的適宜的氣泡尺寸范圍，超過(guò)該范圍，則顆粒脫除率顯著下降。

圖6 不同尺寸含塵單氣泡中顆粒的脫除率Fig. 6 Removal efficiency of particles in different sizes of dustcontaining single bubbles

圖7 所示為含不同粒徑顆粒時(shí)單氣泡上升過(guò)程的形態(tài)變化。由圖可以看出，氣泡中顆粒粒徑增大對(duì)其上升過(guò)程形態(tài)變化幾乎沒(méi)有影響，顆粒脫除率不同的原因在于顆粒本身，粒徑較大的顆粒具有更大的慣性，更容易被脫除。

圖7 含不同粒徑顆粒的單氣泡上升過(guò)程形態(tài)變化Fig. 7 Morphological change of dust-containing single bubbles of different sizes during the rising process

圖8 所示為含兩種粒徑顆粒的不同尺寸單氣泡內(nèi)顆粒脫除過(guò)程。由圖8 可以看出，氣泡尺寸越小，其內(nèi)顆粒在上升初期脫除速率越快，在較短時(shí)間內(nèi)即能達(dá)到較高的脫除率；當(dāng)顆粒粒徑由1.0 μm 增大至2.5 μm 時(shí)，6、8、10 mm 3 種尺寸氣泡呈現(xiàn)出不同的脫除趨勢(shì)，8 mm 和10 mm 氣泡中的顆粒脫除更快，這與圖6 中所表現(xiàn)的顆粒脫除率趨勢(shì)相符合。從圖8 還可以看出，對(duì)于所有尺寸的氣泡，顆粒的脫除速率并未隨時(shí)間的增加而單調(diào)減小，在某些時(shí)刻會(huì)增大，該現(xiàn)象可在圖9 中直觀體現(xiàn)。

圖8 不同尺寸單氣泡內(nèi)顆粒脫除過(guò)程Fig. 8 Removal of particles in single bubbles of different sizes

圖9 所示為含兩種粒徑顆粒時(shí)，不同尺寸氣泡在上升過(guò)程中的顆粒脫除速率。由圖8 和圖9 可知，在氣泡上升過(guò)程中，顆粒的脫除速率并不單調(diào)，不同尺寸氣泡存在快速脫除時(shí)段且各不相同，總體表現(xiàn)為6～8 mm 氣泡中顆粒脫除速率變化相似，10～12 mm氣泡中顆粒脫除速率變化相似，這與氣泡上升過(guò)程中的變形規(guī)律相符。

圖9 顆粒的脫除速率隨時(shí)間變化Fig. 9 Removal rate of particles over time

為探究快速脫除階段的機(jī)理，取各尺寸氣泡快速脫除時(shí)段內(nèi)顆粒濃度分布并作圖(氣泡中顆粒粒徑為1.0 μm)，結(jié)果如圖10 所示。結(jié)果表明，快速脫除階段為氣泡由下凹陷球體變化為扁橢球的階段，這可能是由于該階段氣泡形變最顯著，在變形階段，由于氣泡渦旋的變化，氣相對(duì)顆粒的作用力方向(式5 中Fg)改變，從而加快顆粒與氣液界面發(fā)生碰撞，與已報(bào)道文獻(xiàn)[21]中提出的變形對(duì)顆粒脫除有增強(qiáng)作用相符。

圖10 不同上升時(shí)間單氣泡內(nèi)顆粒濃度分布變化Fig. 10 Changes of particle concentration distribution in single bubbles with different rise time

2.2.2 并行含塵雙氣泡上升過(guò)程中顆粒脫除率 圖11所示為并行含塵雙氣泡上升時(shí)粒徑為1.0 μm 顆粒的脫除過(guò)程。由圖可知，8 mm 與10 mm 并行含塵雙氣泡的脫除過(guò)程基本相同，這與前文含塵單氣泡中顆粒的脫除過(guò)程不同(8 mm 含塵單氣泡具有顯著優(yōu)勢(shì))。該現(xiàn)象表明在一定氣泡尺寸范圍內(nèi)，兩氣泡的共同作用可以增強(qiáng)脫除率；但隨著氣泡尺寸進(jìn)一步增大，該增強(qiáng)作用減弱。

圖11 不同尺寸并行含塵雙氣泡內(nèi)顆粒脫除過(guò)程(dp=1.0 μm)Fig. 11 Removal of particles in parallel dust-containing double bubbles of different sizes(dp=1.0 μm)

圖12 所示為不同尺寸含塵單氣泡與并行含塵雙氣泡上升過(guò)程中粒徑為1.0 μm 顆粒的脫除率。由圖可知，不同尺寸并行含塵雙氣泡中，粒徑為1.0 μm顆粒的脫除率并未全部高于或低于含塵單氣泡，而是交錯(cuò)變化。對(duì)含塵單氣泡和并行含塵雙氣泡而言，小尺寸(6 mm)氣泡均表現(xiàn)出較高的脫除率；但從8 mm 左右的氣泡開始，并行含塵雙氣泡中1.0 μm 顆粒的脫除率增加，氣泡尺寸在10 mm 左右時(shí)1.0 μm 顆粒的脫除率達(dá)到最大，隨后隨氣泡尺寸增大，1.0 μm顆粒的脫除率快速降低，甚至低于含塵單氣泡中顆粒的脫除率。

圖12 不同尺寸含塵單氣泡與并行含塵雙氣泡顆粒脫除率(dp=1.0 μm)Fig. 12 Removal efficiency of particles with different sizes of dust-containing single-bubble and parallel dustcontaining double bubbles(dp=1.0 μm)

圖13 所示為顆粒粒徑為1.0 μm 時(shí)，并行含塵雙氣泡顆粒脫除速率曲線。由圖可知，各尺寸氣泡均有快速脫除階段，6 mm 氣泡在上升初期即快速脫除并達(dá)到高脫除率；8 mm 氣泡與10 mm 氣泡表現(xiàn)出相似的脫除速率變化趨勢(shì)，12 mm 氣泡與14 mm 氣泡的脫除速率規(guī)律相似。由前文可知，快速脫除的原因一方面是由氣泡發(fā)生劇烈變形所造成，而另一方面必然與雙氣泡間的相互作用相關(guān)。

圖13 并行含塵雙氣泡顆粒脫除速率隨時(shí)間變化(dp=1.0 μm)Fig. 13 Removal rate of parallel dust-containing double bubbles varies with time(dp=1.0 μm)

圖14 所示為8 mm 含塵雙氣泡在t為0.020～0.045 s(圖13 中脫除速率增加階段)與t為0.045～0.065 s(圖13 中脫除速率減小階段)的氣泡形態(tài)與相對(duì)運(yùn)動(dòng)趨勢(shì)。由圖可知，當(dāng)t為0.020～0.045 s 時(shí)，兩氣泡相互靠近，呈現(xiàn)出倒八字形；當(dāng)t為0.045～0.065 s時(shí)，在前時(shí)段內(nèi)兩氣泡已非常接近，但受中心射流作用無(wú)法聚并，二者各自的尾流使氣泡發(fā)生遠(yuǎn)離，呈現(xiàn)出類八字形。取靠近過(guò)程的流場(chǎng)矢量圖分析并行含塵雙氣泡上升過(guò)程中的流場(chǎng)變化，結(jié)果如圖15 所示。由圖15 可知，當(dāng)兩氣泡靠近時(shí)，靠近部分流場(chǎng)的擠壓作用可能更有利于顆粒向氣液界面運(yùn)動(dòng)；相比于含塵單氣泡，并行含塵雙氣泡周圍的多個(gè)渦旋強(qiáng)化了顆粒的脫除。

圖14 不同時(shí)段內(nèi)并行含塵雙氣泡運(yùn)動(dòng)趨勢(shì)變化(db=8 mm)Fig. 14 Trend change of parallel dust-containing double bubbles movement during different periods(db=8 mm)

圖15 并行含塵雙氣泡上升過(guò)程中流場(chǎng)變化(db=8 mm)Fig. 15 Flow field changes during the rising process of parallel dust-containing double bubbles(db=8 mm)

12 mm 與14 mm 的并行含塵雙氣泡也經(jīng)歷了靠近過(guò)程，但總體顆粒脫除率卻下降，該現(xiàn)象可從兩氣泡相互作用導(dǎo)致雙氣泡的形變與含塵單氣泡不同進(jìn)行解釋。由圖5 可知，當(dāng)氣泡尺寸較大(12 mm與14 mm)時(shí)，兩氣泡靠近會(huì)使靠近側(cè)被擠壓變形，與含塵單氣泡的對(duì)稱變形及小尺寸并行含塵雙氣泡的對(duì)稱變形不同，內(nèi)側(cè)所含空泡體積減小。由前文分析可知，脫除強(qiáng)化的機(jī)制在于靠近側(cè)的流場(chǎng)擠壓作用，但當(dāng)受擠壓側(cè)變形空泡體積減小后，顆粒會(huì)更集中到兩氣泡的外側(cè)區(qū)域，內(nèi)側(cè)顆粒濃度降低，如圖16 所示，使得大尺寸并行含塵雙氣泡中顆粒的脫除率下降。

圖16 t=0.045 s 與0.065 s 時(shí)含塵雙氣泡內(nèi)顆粒濃度分布(db=12 mm)Fig. 16 Particle concentration distribution in dust-containing double bubbles at t=0.045 s and 0.065 s(db=12 mm)

3 結(jié) 論

本文通過(guò)DPM-VOF 耦合方法模擬研究了不同尺寸含塵單氣泡和并行含塵雙氣泡在初始靜止液體中的上升過(guò)程，以及對(duì)不同粒徑顆粒脫除率的影響，得出結(jié)論如下：

(1) 當(dāng)單氣泡尺寸為6～14 mm 時(shí)，上升氣泡的變形對(duì)顆粒的脫除具有增強(qiáng)作用，形變最復(fù)雜的6 mm氣泡中顆粒的脫除率最高。

(2) 不同粒徑的顆粒中存在能夠使顆粒被高效脫除的適宜的氣泡尺寸范圍，1.0 μm 顆粒在6 mm 單氣泡中表現(xiàn)出最佳的脫除率，2.5 μm 顆粒在8 mm 單氣泡中脫除率最佳。

(3) 并行含塵雙氣泡的初始間距一定時(shí)，氣泡在初始靜止液體中的上升過(guò)程呈現(xiàn)“靠近-分離-靠近”的循環(huán)；受氣泡間相互作用影響，并行含塵大氣泡上升過(guò)程的形變與含塵單氣泡不同。

(4) 不同尺寸并行雙氣泡上升過(guò)程中的相互作用對(duì)顆粒脫除的影響存在差異，8～10 mm 雙氣泡的相互作用強(qiáng)化了粒徑為1.0 μm 顆粒的脫除，而12～14 mm 雙氣泡相互作用則弱化了粒徑為1.0 μm 顆粒的脫除。