程 航，喻九陽，汪 威，孟觀林，王家全

武漢工程大學(xué) 湖北省綠色化工裝備工程技術(shù)研究中心，湖北 武漢 430205

氣浮技術(shù)是一種利用高度分散的微氣泡作為載體粘附污水中的污染顆粒形成混合絮凝體，使得絮凝體的浮力大于其自身重量和上浮阻力，從而上浮至污水表面，達(dá)到分離浮選效果的分離方法［1-2］。氣浮領(lǐng)域中，由于加壓溶氣技術(shù)在加壓條件下空氣溶解度大，并且能產(chǎn)生足量的微氣泡，能夠滿足不同要求的多相混合介質(zhì)的分離，因此應(yīng)用最為廣泛［3］。然而傳統(tǒng)的加壓溶氣氣浮技術(shù)采用接觸式附著方式，該技術(shù)對(duì)絮體顆粒與氣泡直徑要求高，絮體之間容易發(fā)生托附現(xiàn)象，極細(xì)的氣泡很難保持，隨著壓力的降低，氣泡迅速膨脹，難以附著在污染物顆粒上［4］。實(shí)際生產(chǎn)過程中加壓溶氣氣浮技術(shù)一般經(jīng)過加壓溶氣系統(tǒng)、氣泡釋放系統(tǒng)和懸浮物去除系統(tǒng)等系統(tǒng)步驟，操作復(fù)雜過程繁瑣［5-6］。

為了提升氣浮系統(tǒng)的浮選性能，本文基于氣體能量混合管理（gas energy mixing management，GEM）溶氣氣浮技術(shù)設(shè)計(jì)了1種三相混合器裝置［7］。三相混合器通過其特殊結(jié)構(gòu)引入弱旋流場(chǎng)，將氣浮旋流一體化技術(shù)結(jié)合應(yīng)用于污水處理領(lǐng)域［8-9］。該裝置能夠?qū)σ后w、懸浮物、溶解空氣進(jìn)行充分混合，并完成高壓空氣溶解、藥劑分子拉伸、混凝絮凝攪拌、絮凝形成、氣泡晶核成長(zhǎng)和超輕中空絮體形成等一系列步驟。其原理是通過三相混合器，直接在廢水中完成化學(xué)藥劑的投加和加壓溶氣過程，隨著壓力的降低，飽和溶氣在絮體內(nèi)部和周邊生長(zhǎng)成大量微氣泡，微氣泡的生長(zhǎng)擠壓絮體中的水分，氣體和固體絮體形成多孔中空形態(tài)，含水率降低并且自身比重越來越輕，在不借助外力的作用下自行上?。?0-11］。較于常規(guī)溶氣氣浮技術(shù)，三相混合器氣浮技術(shù)產(chǎn)生的氣泡直徑小數(shù)量多，且分布范圍集中，粘附效果更佳，形成的絮體結(jié)實(shí)且含水率低，因此氣浮分離效果更佳［12］。目前，工業(yè)應(yīng)用中只能應(yīng)用經(jīng)驗(yàn)公式對(duì)三相混合器內(nèi)部流場(chǎng)進(jìn)行分析，由于混合器內(nèi)部流場(chǎng)為氣液固三相介質(zhì)混合的弱旋流流態(tài)，難以運(yùn)用精確的數(shù)學(xué)公式對(duì)其復(fù)雜流場(chǎng)進(jìn)行描述。本文將利用計(jì)算流體動(dòng)力學(xué)（computational fluid dynamics，CFD）數(shù)值模擬技術(shù)，對(duì)三相混合器充氣前后的內(nèi)部流場(chǎng)進(jìn)行分析，并引入群體平衡模型對(duì)混合器內(nèi)部氣泡分布特征進(jìn)行研究，探究三相混合器的工作性能［13-15］。

1 數(shù)學(xué)模擬與計(jì)算方法

1.1 幾何模型及網(wǎng)格劃分

根據(jù)三相混合器的結(jié)構(gòu)建立如圖1（a）所示的混合器幾何模型結(jié)構(gòu)示意圖。三相混合器主體結(jié)構(gòu)由內(nèi)筒、外筒、污水入口、充氣入口和出口5 個(gè)部分組成，其中內(nèi)筒上對(duì)稱分布有6 組與內(nèi)筒內(nèi)壁相切的渦流通道，流體通過渦流通道后切向進(jìn)入混合器后在腔壁的作用下形成旋流場(chǎng)。為了直觀地對(duì)混合器內(nèi)部的流場(chǎng)進(jìn)行分析，選取不同截面處的流場(chǎng)進(jìn)行研究，取混合器出口Z0=0 截面為參考截面，分別選取Z1～Z5截面為混合器研究截面，其截面選取示意圖如圖1（b）所示。應(yīng)用數(shù)值計(jì)算前處理軟件ICEM 對(duì)模型計(jì)算域進(jìn)行網(wǎng)格劃分，根據(jù)混合器模型的結(jié)構(gòu)特點(diǎn)，對(duì)混合器模型進(jìn)行了非結(jié)構(gòu)網(wǎng)格劃分，劃分后的混合器網(wǎng)格示意圖如圖1（c）和圖1（d）所示，其中圖1（c）為混合器整體網(wǎng)格示意圖，圖1（d）為局部放大網(wǎng)格示意圖。網(wǎng)格共劃分了190 069 個(gè)節(jié)點(diǎn)和557 251 個(gè)網(wǎng)格單元，為提高網(wǎng)格質(zhì)量，對(duì)網(wǎng)格進(jìn)行了光順處理，并且為了提高網(wǎng)格的計(jì)算精度，對(duì)混合器以及模型入口、腔壁等湍流強(qiáng)度較大等流場(chǎng)特征區(qū)域進(jìn)行網(wǎng)格加密處理。

圖1 混合器幾何模型示意圖：（a）結(jié)構(gòu)，（b）截面選取，（c）整體網(wǎng)格，（d）局部放大網(wǎng)格Fig.1 Schematic diagrams of mixer geometry model：（a）structure，（b）cross section，（c）overall grid，（d）partially enlarged grid

1.2 邊界條件設(shè)置及計(jì)算方法

根據(jù)混合器模型的模擬條件，污水入口采用速度入口邊界，定義入口速度為3 m/s；出口采用壓力出口邊界條件，定義出口相對(duì)壓強(qiáng)為0 MPa；充氣入口采用速度入口邊界條件，設(shè)置充氣入口速度為20 m/s。充氣入口氣體體積分?jǐn)?shù)定義為100%；混合器壁面采用無滑移邊界條件，壁面粗糙度采用默認(rèn)值0.5，同時(shí)采用標(biāo)準(zhǔn)壁面函數(shù)處理邊界湍流。

在數(shù)值模擬計(jì)算過程中，選用雙精度求解器。第一階段對(duì)單相液體進(jìn)行數(shù)值模擬，將充氣入口定義為壁面wall，湍流模型選用RNG k-ε模型，計(jì)算模型采用歐拉模型，利用非穩(wěn)態(tài)求解器進(jìn)行求解，求解計(jì)算混合器流動(dòng)時(shí)間15 s 后停止計(jì)算。第二階段利用單相混合器的計(jì)算結(jié)果作為初始狀態(tài)，激活群體平衡模型，設(shè)置液相為第一相，氣體為第二相，設(shè)置充氣入口為速度入口，進(jìn)行混合器氣液兩相非穩(wěn)態(tài)求解。

1.3 網(wǎng)格無關(guān)性驗(yàn)證

為了驗(yàn)證計(jì)算模型的可靠性，排除網(wǎng)格數(shù)量對(duì)模擬結(jié)果的影響，針對(duì)混合器模型進(jìn)行了網(wǎng)格無關(guān)性驗(yàn)證。如圖2 所示，其中圖2（a）為不同網(wǎng)格數(shù)量模型整體切向速度分布數(shù)值圖，圖2（b）為切向速度數(shù)值圖虛線框所示部位放大示意圖。從圖中可以觀察到，當(dāng)模型網(wǎng)格數(shù)量達(dá)到557 251 個(gè)時(shí)，流場(chǎng)中的切向速度變化甚微，進(jìn)一步增加網(wǎng)格并未對(duì)計(jì)算結(jié)果造成影響，兼顧計(jì)算精度以及計(jì)算效率，選取網(wǎng)格單元數(shù)為557 251 個(gè)時(shí)的模型進(jìn)行計(jì)算。

2 計(jì)算結(jié)果與分析

氣浮旋流一體化技術(shù)中，微氣泡與粒子碰撞粘附過程是影響分離效率的主要因素，對(duì)碰撞粘附過程的研究有助于揭示氣浮技術(shù)機(jī)理。在實(shí)際應(yīng)用中，混合器流場(chǎng)分布特征以及流場(chǎng)中氣泡的分布特性是制約微氣泡與粒子碰撞過程的關(guān)鍵因素，因此本文將針對(duì)未充氣與充氣條件下混合器流場(chǎng)特征以及充氣條件下氣泡分布規(guī)律進(jìn)行研究。

2.1 未充氣與充氣條件下混合器內(nèi)流場(chǎng)分析

2.1.1 混合器內(nèi)部流場(chǎng)速度矢量圖 圖3 為充氣條件下混合器內(nèi)部流場(chǎng)橫截面與縱截面速度矢量圖。從圖3 中可以觀察到，污水從入口進(jìn)入混合器后運(yùn)動(dòng)方向發(fā)生改變且運(yùn)動(dòng)速度降低，這是因?yàn)槲鬯M(jìn)入混合器后在圓柱面作用下產(chǎn)生分流，并且在腔壁處受到了較大的阻力所致。由于混合器的特殊結(jié)構(gòu)，污水通過引流槽和渦流通道切向進(jìn)入內(nèi)筒并在內(nèi)筒中形成旋流。充氣條件下，空氣高速射流進(jìn)入混合器后在旋流場(chǎng)的作用下在中心處形成空氣柱，高速運(yùn)動(dòng)的空氣柱帶動(dòng)周圍流體在軸心處形成高速區(qū)。由于離心場(chǎng)的存在，高速區(qū)中的氣體不斷與旋流場(chǎng)中的流體發(fā)生物質(zhì)之間的交換，并且由于污水通過渦流通道射流進(jìn)入內(nèi)筒中，在運(yùn)動(dòng)過程中與高速區(qū)中的流體發(fā)生碰撞，因此極大地增大了污水中的粒子與氣泡的碰撞幾率。

2.1.2 切向速度 旋流場(chǎng)流動(dòng)中，由切向速度引起的離心場(chǎng)一方面導(dǎo)致微氣泡產(chǎn)生自轉(zhuǎn)，使氣泡間產(chǎn)生聚并和破碎，另一方面旋轉(zhuǎn)離心場(chǎng)加速了各相之間的運(yùn)移過程，加快混合進(jìn)程，因此切向速度是旋轉(zhuǎn)混合分離技術(shù)中重要參數(shù)。圖4 為未充氣與充氣條件下混合器內(nèi)筒各截面的切向速度分布數(shù)值圖。從圖4（a）和圖4（b）可以發(fā)現(xiàn)由于近壁面的粘度作用，靠近近壁面的切向速度躍遷至較大值，隨著半徑的減小，切向速度呈現(xiàn)遞減的趨勢(shì)，至軸心處切向速度為零，由此表明混合器內(nèi)流場(chǎng)中的切向速度符合強(qiáng)制渦的運(yùn)動(dòng)特征。

圖2 網(wǎng)格無關(guān)性驗(yàn)證：（a）不同網(wǎng)格數(shù)量模型切向速度，（b）局部放大示意圖Fig.2 Verification of grid-independent：（a）tangential velocity of different mesh quantity models，（b）schematic diagram of partial magnification

圖3 混合器內(nèi)部流場(chǎng)速度矢量圖：（a）橫截面速度，（b）縱截面速度Fig.3 Velocity vectographs of mixer internal flow field：（a）cross-sectional velocity，（b）longitudinal sectional velocity

圖4 混合器各截面切向速度分布圖：（a）未充氣條件下，（b）充氣條件下Fig.4 Tangential velocity profiles of each cross section of mixer：（a）non-aerated conditions，（b）aerated conditions

充氣狀態(tài)下，截面Z3=110 mm、Z4=160 mm、Z5=200 mm 處切向速度略有減小，而截面Z1=10 mm、Z2=60 mm 處切向速度未發(fā)生明顯變化。由于氣體軸向進(jìn)入混合器后在旋流場(chǎng)作用下形成中心氣柱，氣柱與流體之間的氣液耦合包括氣液摩擦，氣液混合以及氣液運(yùn)移與擴(kuò)散導(dǎo)致了混合器上半段處的能量消耗，因此充氣后混合器上半段切向速度略有減小。對(duì)比充氣與未充氣條件下流場(chǎng)特征，未充氣與充氣條件下切向速度在流場(chǎng)中的分布規(guī)律基本相同，表明充氣條件下并未對(duì)切向速度場(chǎng)產(chǎn)生劇烈擾動(dòng)。穩(wěn)定的切向速度場(chǎng)提供了穩(wěn)定的離心場(chǎng)，有利于氣泡的遷移擴(kuò)散和氣液間的混合，并且穩(wěn)定的離心場(chǎng)避免了更多的微氣泡發(fā)生剪切破碎，有助于氣泡與粒子間的碰撞粘附。

圖5 混合器各截面軸向速度分布圖：（a）未充氣條件下，（b）充氣條件下Fig.5 Axial velocity profile of each cross section of mixer：（a）non-aerated conditions，（b）aerated conditions

2.1.3 軸向速度 圖5 為未充氣和充氣條件下混合器內(nèi)部流場(chǎng)中各截面軸向速度分布數(shù)值圖。從圖5（a）中可以發(fā)現(xiàn)，近壁面處軸向速度為零，隨著半徑的減小軸向速度呈現(xiàn)遞增的趨勢(shì)，在中心軸線處軸向速度達(dá)到最大值。軸向方向上，越靠近出口處軸向速度越大，在Z5=200 mm 截面處軸向速度接近于零，然而在靠近出口處的Z1截面處軸向速度最大值可達(dá)到5.6 m·s-1。軸向速度的速度梯度分布以及不同截面處軸向速度分布跨度大，導(dǎo)致了下層的流體總是快于上層的流體流出混合器，因此軸向分布上不同截面處的粒子間難以發(fā)生有效的碰撞粘附。

充氣條件下軸向速度分布如圖5（b）所示，由于氣體通過從上至下射流的方式注入混合器，氣體進(jìn)入混合器后在離心場(chǎng)的作用下在中心軸線附近形成中心氣柱，并在混合器中心氣柱處形成一個(gè)高速區(qū)，高速區(qū)內(nèi)的流體軸向速度遞增至較大值并帶動(dòng)周圍流體產(chǎn)生加速。充氣后流場(chǎng)中的軸向速度整體分布上有明顯提升，徑向方向上軸向速度依舊呈現(xiàn)逐漸遞增的趨勢(shì)，但是靠近充氣入口處的增長(zhǎng)幅度明顯大于下半段，靠近中心氣柱處Z4與Z5截面處的軸向速度反而大于Z1和Z2截面處的軸向速度，在離心場(chǎng)的共同作用下，上層流體中的粒子更快地進(jìn)入下層流體中并與粒子發(fā)生接觸，因此提供了更多的碰撞幾率。而混合器中軸向速度場(chǎng)的強(qiáng)化使流場(chǎng)中產(chǎn)生了更強(qiáng)的湍流強(qiáng)度，導(dǎo)致流場(chǎng)中湍流彌散作用加強(qiáng)，彌散現(xiàn)象的加強(qiáng)促進(jìn)了粒子間的運(yùn)移擴(kuò)散，有利于提高微氣泡與粒子之間的碰撞聚并粘附效率，因此在對(duì)軸向速度場(chǎng)的研究表明充氣條件下提升了混合器氣浮選性能。

2.1.4 徑向速度 圖6 為混合器流場(chǎng)中各截面徑向速度分布圖，如圖6（a）所示，徑向速度在3 個(gè)速度分量中速度最小，同于切向速度與軸向速度分布特征，徑向速度數(shù)值上也呈對(duì)稱分布與速度梯度分布，由于混合器內(nèi)部流場(chǎng)符合強(qiáng)制渦的流場(chǎng)特征，因此靠近中心軸線處徑向速度大小趨近于零。圖6（b）為充氣條件下徑向速度分布圖，在充氣條件下，徑向速度整體大小上未發(fā)生明顯的變化，但是由于中心氣柱的形成以及中心軸線處軸向速度和湍流強(qiáng)度的突然增加帶動(dòng)周邊流體在中心處產(chǎn)生方向相反位置相對(duì)的徑向流動(dòng)。徑向流動(dòng)的產(chǎn)生表明中心氣柱處存在著氣泡的運(yùn)移與擴(kuò)散，中心氣柱處的氣泡不斷地向流體中擴(kuò)散，在流場(chǎng)中的離心場(chǎng)和軸向速度場(chǎng)等作用下與流體中的粒子發(fā)生碰撞粘附，因此充氣條件下極大地提高了微氣泡與粒子的碰撞幾率。

圖6 混合器各截面徑向速度分布圖：（a）未充氣條件下，（b）充氣條件下Fig.6 Radial velocity profiles of each cross section of mixer：（a）non-aerated conditions，（b）aerated conditions

通過對(duì)充氣前后混合器內(nèi)流場(chǎng)分布特征進(jìn)行分析，可以發(fā)現(xiàn)充氣前后混合器內(nèi)三相速度均呈現(xiàn)速度梯度分布，在旋流場(chǎng)中由速度梯度引起的速度差能夠?yàn)椴煌Ｗ犹峁┎煌那邢蚣袅Α㈦x心力、徑向力，因此的激增和湍流強(qiáng)度的突變，以及在中心氣柱處產(chǎn)生的徑向流動(dòng)促進(jìn)了粒子間的遷移擴(kuò)散，有利于氣泡與顆粒間的碰撞粘附。并且軸向速度的增大，減少了流體在混合器內(nèi)的水力停留時(shí)間，加速混合進(jìn)程。

2.2 混合器內(nèi)流場(chǎng)氣泡分布特征

溶氣氣浮技術(shù)中，氣泡分布、氣泡大小及氣泡的數(shù)量等條件都是制約混合分離效果的重要因素，而混合器流場(chǎng)對(duì)氣泡的作用則影響著氣泡的特性，因此本文將針對(duì)充氣條件下混合器流場(chǎng)中的微氣泡進(jìn)行研究。圖7 為釋壓后混合器內(nèi)流場(chǎng)微氣泡粒徑分布柱狀圖，通過觀察可以發(fā)現(xiàn)流場(chǎng)內(nèi)產(chǎn)生的微氣泡粒徑D 主要分布在100～200 μm之間，并且該粒徑范圍下氣泡數(shù)量占?xì)馀菘偭康?8%左右，與傳統(tǒng)的氣浮方法相比，該粒徑下的微氣泡直徑小，分布范圍窄，粘附效果更佳，更適用于溶氣氣浮技術(shù)。

圖7 氣泡粒徑分布百分比柱狀圖Fig.7 Histogram of bubble particle size distribution

利用數(shù)值模擬軟件對(duì)混合器內(nèi)部氣泡數(shù)量進(jìn)行追蹤，通過統(tǒng)計(jì)得到混合器內(nèi)部溶解氣體含量為0.002 2 kg/m3，流體的含量為1.786 4 kg/m3，利用公式計(jì)算出流場(chǎng)內(nèi)溶解氣含率達(dá)到50.46%，并且對(duì)出口處的氣泡含量進(jìn)行追蹤統(tǒng)計(jì)，得到出口處氣泡攜帶率為19.77%，表明通過混合器流場(chǎng)的作用，氣體的溶解度高，產(chǎn)生氣泡的數(shù)量多，并且混合完成后氣泡攜帶率高。通過對(duì)流場(chǎng)中氣泡的追蹤研究，驗(yàn)證了在混合器作用下用下流場(chǎng)產(chǎn)生的微氣泡直徑小、分布范圍窄、氣泡數(shù)量多，并且氣泡攜帶率高，因此在混合粘附過程中效果更佳。

3 結(jié) 論

借助數(shù)值模擬技術(shù)對(duì)未充氣與充氣條件下混合器內(nèi)部流場(chǎng)進(jìn)行模擬，通過對(duì)內(nèi)流場(chǎng)的對(duì)比研究分析，可得出以下結(jié)論：

1）充氣前后，混合器模型內(nèi)部流場(chǎng)的分布特征表明了流場(chǎng)在流動(dòng)過程中粒子之間能夠很好的發(fā)生碰撞粘附等過程。充氣條件下對(duì)流場(chǎng)的強(qiáng)化作用，表明了充氣后的流場(chǎng)能夠極大地促進(jìn)粒子的遷移與氣泡的擴(kuò)散，有利于提高微氣泡與粒子之間的碰撞聚并粘附效率。

2）在充氣條件下，混合器內(nèi)產(chǎn)生的微氣泡粒徑主要分布在100～200 μm 之間，該粒徑范圍下氣泡數(shù)量占?xì)馀菘偭康?8%，通過在混合器旋流場(chǎng)的作用下，流體中氣溶率高達(dá)50.46%，并且出口處的氣泡攜帶率達(dá)到19.77%。結(jié)果表明通過混合器作用下流場(chǎng)產(chǎn)生的微氣泡直徑小，分布范圍窄，并且氣泡數(shù)量多，因此在實(shí)際的工業(yè)應(yīng)用中氣泡的粘附攜帶性能更好，氣浮選的效果更佳。

3）從數(shù)值模擬的角度驗(yàn)證了三相混合器內(nèi)部流場(chǎng)分布特征與工作特性，為下一步三相混合器混合性能的實(shí)驗(yàn)研究以及結(jié)構(gòu)優(yōu)化設(shè)計(jì)提供理論基礎(chǔ)。