張 巖 任萬龍,? 張旭輝,?,2) 魯曉兵,?

* (中國科學(xué)院力學(xué)研究所流固耦合系統(tǒng)力學(xué)重點實驗室,北京 100190)

? (中國科學(xué)院大學(xué)工程科學(xué)學(xué)院,北京 100049)

引言

全球工業(yè)高速發(fā)展,各國對能源的需求也越來越大,深海能源的開發(fā)與利用逐漸成為世界各國的研究重點.深海礦產(chǎn)資源是深海能源的重要組成部分,我國“十四五”規(guī)劃和2035 年遠(yuǎn)景目標(biāo)綱要也將深海礦產(chǎn)資源的開發(fā)列為重要發(fā)展方向[1].集礦機(jī)與管道結(jié)合的水力輸送系統(tǒng)是目前最有應(yīng)用前景的深海采礦系統(tǒng),它具有回采率高,產(chǎn)能高,操作方便等優(yōu)點[2].垂直管水力輸送是其中的關(guān)鍵環(huán)節(jié),輸送過程含“粗顆?！薄ⅰ皩捈壟洹焙汀案邼舛取钡湫吞卣鞯亩嘞嗔鲃右?guī)律是亟需解決的核心科學(xué)問題之一,是提升輸送效率與保障輸送安全的理論基礎(chǔ).

與細(xì)顆粒漿體相比,含粗顆粒的固液混合物所具有的一個顯著特征就是由固液速度滑移導(dǎo)致的顆粒滯留現(xiàn)象.顆粒滯留會使管道中的局部顆粒濃度顯著增加,根據(jù)相關(guān)文獻(xiàn)報道,垂直管水力輸送中顆粒局部濃度相比于輸送入口濃度增幅可達(dá)150%～250%[2-3],很可能造成管道堵塞.長距離管道輸送過程一旦發(fā)生堵塞,必然造成系統(tǒng)癱瘓.因此,研究輸送過程的顆粒滯留規(guī)律至關(guān)重要.一方面可在稠密粗顆粒固液兩相流動規(guī)律這一基礎(chǔ)科學(xué)問題上有所拓展,另一方面可為深海采礦輸送設(shè)備的設(shè)計選用以及系統(tǒng)的穩(wěn)定運行提供科學(xué)依據(jù).

由于顆粒滯留產(chǎn)生于固-液相間速度滑移,因此早期顆粒滯留的研究主要關(guān)注固-液相間速度滑移的變化規(guī)律.如夏建新等[4]提出顆粒滯留效應(yīng)對垂直管顆粒輸送過程影響較大,當(dāng)入口輸送速度過小時,很大概率造成管道堵塞.他們采用流態(tài)化原理研究了粗顆粒體積濃度對固液滑移速度的影響規(guī)律,并給出垂直管中顆粒局部濃度的預(yù)測關(guān)系式.Van Wijk 等[5]在實驗中發(fā)現(xiàn)粗顆粒滯留會導(dǎo)致垂直管內(nèi)出現(xiàn)顆粒栓塞,他們將顆粒栓塞的產(chǎn)生歸因于密度波的傳播,并通過實驗測量了密度波的波速.宋躍文等[6]著重研究粗顆粒在垂直管道中的滑移速度,基于高速攝影技術(shù)和粒子圖像測速技術(shù),采用實際深海礦石包括錳結(jié)核,富鈷結(jié)殼,多金屬硫化物等進(jìn)行室內(nèi)實驗,發(fā)現(xiàn)相同材料和粒徑時,滑移速度隨入口輸送速度增加而減小;相同粒徑和入口輸送速度時,滑移速度隨顆粒密度增大而增大.

近年來,在深海采礦試采中發(fā)現(xiàn)垂直管粗細(xì)顆?；旌陷斔涂赡軙龃蟪霈F(xiàn)栓塞,堵塞的可能性[7],一些研究人員在此方面開展了相關(guān)研究.如Li 等[8]基于歐拉-矩方法研究混合顆粒運動過程,發(fā)現(xiàn)粗細(xì)顆粒群出現(xiàn)明顯分層.Dai 等[9]在數(shù)值模擬中發(fā)現(xiàn)粗細(xì)顆粒混合輸送(粗細(xì)顆粒粒徑比為6)時顆粒局部濃度相比于均勻粒徑輸送時更大.而劉磊[10]則提出,在相同輸送濃度時,相比于均勻粒徑條件,有級配顆粒群輸送時管道內(nèi)局部顆粒濃度降低.Cú?ez 等[11-12]單獨研究兩種粒徑顆粒(粒徑 6 mm 和4.8 mm)混合輸送時顆粒分離現(xiàn)象,發(fā)現(xiàn)粗細(xì)顆?；旌想A段顆粒間碰撞力和顆粒濃度顯著增加.Ren 等[13]采用類似[11]的研究方案,從顆粒受力角度解釋了粗細(xì)顆粒流動分離的機(jī)理.Zhang 等[2]在前人的基礎(chǔ)上進(jìn)一步分析粗細(xì)顆?；旌陷斔蛯︻w粒滯留的影響規(guī)律,發(fā)現(xiàn)粗細(xì)顆?；旌陷斔蜁r由于粗細(xì)顆粒滑移速度不同,會出現(xiàn)細(xì)顆粒繞過粗顆粒向上運動導(dǎo)致粗細(xì)顆粒流動分離的現(xiàn)象,此過程會增強(qiáng)顆粒滯留效應(yīng).主要是因為細(xì)顆粒在粗顆粒群空隙中運動時會碰撞擠壓,一方面增大能量損耗,使顆粒群運動速度降低.另一方面,顆粒群局部體積濃度增大,壁面的限制效應(yīng)增強(qiáng),從而增加了顆粒在壁面處的碰撞能量損耗.

以上研究可以發(fā)現(xiàn),目前針對混合粒徑條件下的顆粒滯留機(jī)理的研究仍未完善,現(xiàn)有研究出現(xiàn) 了一些相互矛盾的結(jié)論.此外,目前的研究大多著眼于宏觀的顆粒濃度變化規(guī)律,統(tǒng)計局部顆粒濃度時均值受輸送參數(shù)如輸送速度,輸送濃度,顆粒粒徑等的影響規(guī)律,未從顆粒運動特性及固-液兩相相互作用角度解釋粗細(xì)顆粒分離的力學(xué)機(jī)制.因此,本文擬采用計算流體力學(xué)-離散元方法(CFD-DEM)開展相關(guān)研究.CFD-DEM 方法[14]兼顧了顆粒流直接數(shù)值模擬[15]和歐拉-歐拉方法[16]的優(yōu)點,可以追蹤單個顆粒的運動軌跡,求解顆粒之間的碰撞,并且所需計算資源適中.基于數(shù)值模擬,著重探究雙尺寸顆粒群輸送過程顆粒體積濃度,顆粒群運動特性,兩相相互作用等的變化,解釋雙尺寸顆粒群混合及分離過程的力學(xué)機(jī)制,為實際深海采礦水力輸送的施工設(shè)計提供參考.

1 方法

由于本文研究的問題存在“寬級配”和“高濃度”特征,以實際深海采礦試采為例,顆粒體積濃度范圍約為0.05～0.2,顆粒尺寸范圍約為1～30 mm[1,2,17-18].如此高固相含量條件下,顆粒運動對流場反作用不可忽略[19].此外,由于顆粒尺寸大小不一,在CFDDEM 數(shù)值模擬中單個顆?？赡苷紦?jù)多個流體網(wǎng)格,導(dǎo)致局部顆粒濃度的計算出現(xiàn)非物理情況,如顆粒體積濃度大于1.0.目前最常用的顆粒體積濃度計算方法為顆粒中心法(particle centroid method,PCM)[20],即將單個流體網(wǎng)格內(nèi)的顆粒總體積除以網(wǎng)格體積即為此網(wǎng)格內(nèi)的顆粒體積濃度.此方法需要保證顆粒尺寸與網(wǎng)格尺寸比小于0.2,否則計算誤差較大[21].本文作者在前期研究中拓展了顆粒中心方法的適用范圍,使其可以應(yīng)用于粗顆粒固液兩相流的數(shù)值模擬中[22].考慮到本文應(yīng)用場景中粗細(xì)顆粒尺寸差異較大,作者對前期數(shù)值方法進(jìn)一步優(yōu)化,保證同時存在粗細(xì)顆粒時顆粒體積濃度計算的準(zhǔn)確性.下面介紹本文數(shù)值模型.

1.1 控制方程

連流體相采用歐拉方法描述,通過求解流體相的質(zhì)量守恒方程和動量守恒方程獲得速度和壓力分布場.在管道流中,兩相流是低速流動,因此流體相可看作不可壓縮流體.具體控制方程為[23]

其中ρf是流體密度,αf是流體體積分?jǐn)?shù),u為流體速度,Fpf是表征固液相間耦合的動量交換源項,τf是流體黏性應(yīng)力張量,g為重力加速度.

固液動量交換源項為

其中Vc是網(wǎng)格體積,fd,i是顆粒所受拖曳力,flift,i是升力,fvm是虛擬質(zhì)量力.

流體黏性應(yīng)力張量為

其中 νf是流體黏度,I為單位張量.

顆粒相采用拉格朗日方法描述,以單個顆粒為例,顆粒在流體中的運動主要包括平動和轉(zhuǎn)動.顆粒運動的驅(qū)動力來自于流體相的拖曳力,壓力梯度力,虛擬質(zhì)量力以及顆粒之間的碰撞力等.在數(shù)值模擬中,求解顆粒運動方程的時間步長通常為 10-4～10-6s量級,在如此小的時間步長內(nèi),可以只考慮顆粒與其附近的流體及顆粒的相互作用[24].因此,顆粒運動的控制方程為

其中vi是顆粒的平動速度,ωi是顆粒的轉(zhuǎn)動速度,fpf,i是固-液相互作用力,fc,ij是顆粒之間的碰撞力,kc為與顆粒i同時碰撞的顆粒數(shù)量,Mij為顆粒之間碰撞力產(chǎn)生的力矩.

固液相互作用力主要包括流體拖曳力、壓力梯度力、虛擬質(zhì)量力、Basset 力、Magnus 力和Saffman 力等.通過量級比較,垂直管水力輸送中,流體拖曳力、壓力梯度力和虛擬質(zhì)量力起主要作用,因此本文只考慮上述3 個力的作用.流體拖曳力模型采用ErgunWenYu 模型[25]

其中kd是拖曳力作用系數(shù),Cd是阻力系數(shù),Vp是顆粒體積,Re為雷諾數(shù).

壓力梯度力和虛擬質(zhì)量力的計算式分別為

其中D/Dt是物質(zhì)導(dǎo)數(shù).

1.2 顆粒體積濃度計算

前期研究中,作者采用虛擬密度分布函數(shù)方法求解粗顆粒的體積濃度場.具體實現(xiàn)方式為將顆粒的密度進(jìn)行重分布,使密度分布為空間中任意點與顆粒中心點距離的連續(xù)函數(shù)

其中ρp是顆粒密度,r是顆粒半徑,β是空間任意一點與顆粒中心點的距離.隨后對重分布的密度函數(shù)在流體域內(nèi)積分獲得顆粒的體積濃度場[22].在程序?qū)崿F(xiàn)時,求解以PCM 方法計算得到的顆粒體積濃度場為初值的擴(kuò)散方程,即可得到粗顆粒的體積濃度分布

其中αp是顆粒體積濃度,αp,PCM是基于PCM 方法計算得到的顆粒體積濃度,τ為擴(kuò)散方程的時間參數(shù).上述顆粒體積濃度場的計算效果見圖1.

圖1 虛擬密度分布函數(shù)求解粗顆粒體積濃度場示意圖Fig.1 Diagram of the particle concentration field based on the virtual mass distribution method

本文中涉及的研究問題同時存在兩種尺寸顆粒,并且前期研究中發(fā)現(xiàn),若顆粒尺寸較小,虛擬密度分布函數(shù)方法不適用,主要是因為密度重分布的方式需要顆粒尺寸與網(wǎng)格尺寸相當(dāng).因此,本文將式(14)修改為

其中Γ是擴(kuò)散系數(shù),為空間的函數(shù).若顆粒尺寸與網(wǎng)格尺寸之比小于 0.2,Γ設(shè)為較小的值,本文設(shè)為0.001,否則Γ值為 1.0.通過控制Γ的取值,實現(xiàn)了PCM 與虛擬密度分布函數(shù)方法的切換.

1.3 歐拉場到拉格朗日場映射

在歐拉-拉格朗日方法中,流體相采用有限體積方法求解,流體相場量如速度,壓力等存儲于網(wǎng)格單元中心點(內(nèi)部區(qū)域)或網(wǎng)格面中心點(邊界區(qū)域).但是,顆粒在流場中運動時顆粒中心點與網(wǎng)格單元中心點或網(wǎng)格面中心點一般不重合,需要插值獲得顆粒中心點的流體相場量以計算顆粒受力.這主要是因為顆粒受力計算模型中需要獲得顆粒中心位置處的流體相速度、壓力、黏度等(式(7)～式(12)).顆粒尺寸不同,使得顆粒所占據(jù)的網(wǎng)格數(shù)量也會不同,而顆粒中心位置處的流體相場量一般與顆粒的影響范圍有關(guān).因此,本文提出一種基于顆粒尺寸的歐拉場到拉格朗日場的權(quán)函數(shù)插值方法.具體實現(xiàn)思路如下: 首先根據(jù)顆粒中心位置及顆粒尺寸獲取顆粒所占據(jù)網(wǎng)格的編號,隨后計算顆粒中心位置與所占據(jù)網(wǎng)格中心點的距離,接著計算權(quán)函數(shù),權(quán)重為顆粒中心位置與網(wǎng)格中心點之間距離的倒數(shù)

其中wi為權(quán)重,xp為顆粒中心點位置矢量,xi為網(wǎng)格中心點位置矢量,Nc為顆粒所占據(jù)的網(wǎng)格數(shù)量.顆粒所占據(jù)網(wǎng)格數(shù)量的判斷標(biāo)準(zhǔn)為網(wǎng)格中心點位于顆粒尺寸范圍之內(nèi),如圖2 所示,顆粒共占據(jù)7 個網(wǎng)格.獲得各網(wǎng)格編號后,顆粒中心位置處的流體相場量為(以速度為例)

圖2 權(quán)函數(shù)插值法中顆粒所占據(jù)網(wǎng)格示意圖Fig.2 Diagram of the weighted function method

在某些特殊條件下,如顆粒尺寸較小,按圖可能無法獲得對應(yīng)的占據(jù)網(wǎng)格.此時,令wi=0,即u(xp)=u(xi).綜合上述兩種情況即可獲得不同顆粒尺寸條件下,顆粒中心位處的流體相場量.

本文采用C++程序語言,基于CFDEM 開源平臺進(jìn)行前述顆粒體積濃度計算模型和歐拉場到拉格朗日場映射模型的程序植入.CFDEM 為流體計算開源軟件OpenFOAM 與顆粒計算開源軟件LIGGGHTS的耦合平臺.

1.4 模型驗證

本文研究的是粗細(xì)顆粒群混合輸送過程中的分離問題,其中最重要的影響因素之一就是顆粒尺寸.并且數(shù)值模型中也針對顆粒尺寸的影響統(tǒng)一了粗細(xì)顆粒的濃度計算方法以及歐拉場到拉格朗日場的映射法.因此,本文選用流化床案例進(jìn)行模型驗證,分別考慮不同顆粒尺寸對最小流化速度以及床層壓降的影響.圖3 為驗證案例的幾何尺寸及顆粒初始狀態(tài).幾何模型為一垂直管,直徑 27.6 mm,高度55 mm.初始時刻在管道底部自由堆積一定數(shù)量的顆粒,顆粒直徑取值分別為 1 mm,2 mm,3 mm 和4 mm,初始顆粒堆積高度為 15 mm.從圓管底部通入水,密度 1000 kg/m3,黏度 0.001 Pa·s,入口表觀速度從 0 m/s 逐漸增大,統(tǒng)計顆粒床層壓降隨時間的變化.圖3 中圓管底部壓力梯度為零,圓管頂部壓力為固定值 0 Pa,管壁為無滑移邊界.當(dāng)?shù)撞咳肟诹魉佥^小時,流體拖曳力小于顆粒浮重力,流體不足以托起顆粒床層,因此顆粒群高度保持不變,顆粒床層壓降隨著入口流速增加而增加.當(dāng)流體拖曳力剛好與顆粒浮重力相等時,顆粒床層到達(dá)最小流化狀態(tài),此時的入口流速稱為最小流化速度,并且顆粒床層壓降到達(dá)最大值.此狀態(tài)后若進(jìn)一步增大入口速度,只會使顆粒床層高度增加但床層壓降保持不變.因此,可通過顆粒床層壓降隨入口速度的變化曲線來確定最小流化速度.

圖3 垂直管流化床算例示意圖Fig.3 Geometry representation of the simulated fluidized bed

Ergun 等[26]基于實驗研究了顆粒床層壓降與入口流速的關(guān)系,并結(jié)合理論推導(dǎo)給出壓降與入口流速的計算式

其中uz為入口流速.此外,最小流化速度可通過下式計算[13]

其中Remf為最小流化速度雷諾數(shù).基于式 (18)～式(20) 即可得到顆粒床層壓降隨入口流速變化的解析解曲線.圖4 為數(shù)值模擬結(jié)果與解析解的對比,其中圖4(a) 為顆粒尺寸為 1 mm 時顆粒床層的壓降曲線,根據(jù)壓降曲線的拐點可以確定數(shù)值模擬的最小流化速度值.圖4(b) 為不同顆粒尺寸時,根據(jù)壓降曲線拐點獲得的最小流化速度與解析解的對比.根據(jù)對比結(jié)果可以發(fā)現(xiàn)本文數(shù)值模型可以準(zhǔn)確預(yù)測不同尺寸顆粒床的壓降和最小流化速度.

圖4 顆粒床層壓降與最小流化速度對比Fig.4 Comparison of particle pressure drop and minimum fluidization velocity

2 數(shù)值模擬工況設(shè)計

圖5 為本文數(shù)值模擬工況示意圖,由于本文著重關(guān)注的點為顆粒尺寸效應(yīng)引起的顆粒分離.為了簡化研究過程,只考慮兩種顆粒尺寸.初始狀態(tài)時,細(xì)顆粒群在下粗顆粒群在上,粗細(xì)顆粒的尺寸分別為d1和d2,初始濃度0.3,初始高度0.1 m,初始顆粒群速度為零.水平方向為重力加速度方向.與驗證算例設(shè)置類似(圖3),管道左側(cè)為入口,入口流體速度恒定,管道右側(cè)為出口.管道直徑 100 mm,長度3.6 m,顆粒密度2000 kg/m3,流體密度 1000 kg/m3,流體黏度0.001 Pa·s.粗細(xì)顆粒尺寸不定.

圖5 垂直管顆粒分離算例示意圖Fig.5 Geometry representation of the simulated particle segregation case

3 結(jié)果討論

3.1 粗細(xì)顆粒群分離過程

圖6 為管道內(nèi)粗細(xì)顆粒群在不同時刻的位置圖.細(xì)顆粒直徑 5 mm,粗顆粒直徑 10 mm,管道入口流速1 m/s.粗細(xì)顆粒群在水流攜帶作用下向前運動,同時細(xì)顆粒運動速度大于粗顆粒,因此位于后方的細(xì)顆粒群逐漸與前方粗顆粒群混合,使得粗顆粒群位置處顆粒濃度增大.由圖6(b)和圖6(c)可以發(fā)現(xiàn),細(xì)顆粒群并未整體越過粗顆粒群,而是前部與粗顆粒群混合,后部逐漸與前部分離.這是因為前部混合后,混合區(qū)局部濃度增大,流體在管道內(nèi)過流面積減小,進(jìn)而增大局部流體速度以及流體拖曳力.因此混合區(qū)域的顆粒群運動速度增大,與后方細(xì)顆粒群分離.管道內(nèi)出現(xiàn)兩個區(qū)域: 細(xì)顆粒群區(qū)域和粗細(xì)顆粒群混合區(qū)域.隨著顆粒群繼續(xù)向前方運動,粗細(xì)顆粒群混合區(qū)域雖然整體運動速度增大,但是其內(nèi)部兩種粒徑顆粒由于滑移速度不同,細(xì)顆粒群同樣會逐漸超越粗顆粒群最終運動到粗顆粒群前方,出現(xiàn)粗細(xì)顆粒群徹底分離的現(xiàn)象(圖6(l)).前方混合區(qū)域的顆粒逐漸分離后,后方細(xì)顆粒群運動速度大于前方粗顆粒群,又會逐漸出現(xiàn)新的混合區(qū).新混合區(qū)同樣會逐漸分離,最終細(xì)顆粒群完全越過粗顆粒群,這就是粗細(xì)顆粒群的分離過程.

圖6 垂直管內(nèi)粗細(xì)顆粒群在不同時刻的位置圖Fig.6 Distribution of the binary particles at different time

為了更清晰地顯示粗細(xì)顆粒混合及分離過程的濃度變化,圖7 給出了垂直管軸線上不同時刻顆粒群濃度曲線.需要注意的是,雖然粗細(xì)顆粒群初始時刻濃度都為0.3,但是在輸送過程中顆粒群存在軸向擴(kuò)散,使得顆粒群寬度增加,局部濃度相應(yīng)減少.因此顆粒群濃度在輸送過程中一些區(qū)域小于0.3.從圖7(b)可以清晰看到軸向高度z=0.75 m 左右濃度曲線出現(xiàn)尖峰,最大濃度值約為0.37.這就是因為細(xì)顆粒群前部與粗顆粒群后部混合.此外,由于顆粒群擴(kuò)散,細(xì)顆粒群后方濃度降低,前方混合區(qū)濃度高而后方細(xì)顆粒群濃度低,導(dǎo)致前方混合區(qū)運動速度增加,引起混合區(qū)與細(xì)顆粒群分離 (圖7(d)和圖7(f)).從圖7(j)可以看到,前方高濃度區(qū)域極值降低,寬度增加,并且濃度尖峰處分為兩段.這意味著此時細(xì)顆粒群越過粗顆粒群,并且開始分離.最終粗細(xì)顆粒群完全分離,前方濃度曲線由單峰完全變?yōu)殡p峰.并且同樣由于顆粒群擴(kuò)散作用,粗細(xì)顆粒群分離后,濃度約為 0.1.

圖7 垂直管軸線處不同時刻顆粒體積群濃度曲線Fig.7 Particle concentration at different times at the vertical pipe axis

從圖7(f)到圖7(k)可以發(fā)現(xiàn)濃度曲線的后方尖峰極值在逐漸增大,正如前文所述,這是由于后方的細(xì)顆粒群與前方分離后的粗顆粒群相交,重復(fù)最初粗細(xì)顆粒群相混到完全分離的過程.

3.2 粗細(xì)顆粒群分離機(jī)理

稠密固液兩相流動中,有兩種物理效應(yīng)需要著重關(guān)注,分別是固液相互作用與顆粒碰撞作用[27].本節(jié)將從上述兩個物理效應(yīng)角度探究粗細(xì)顆粒群分離過程的顆粒運動特性,并進(jìn)一步研究粗細(xì)顆粒群的分離機(jī)理.首先定義兩個無量綱數(shù): 流固相互作用應(yīng)力和顆粒碰撞應(yīng)力,分別表征流固相互作用強(qiáng)度以及顆粒碰撞作用強(qiáng)度

其中Np為統(tǒng)計顆粒數(shù)量,Vpn為顆粒n的體積.在顆粒雷諾數(shù)遠(yuǎn)大于1 時,兩相流中顆粒弛豫時間為[28-29]

其中s=ρp/ρf為顆粒與流體密度比.因此,后文以顆粒弛豫時間為基準(zhǔn)將時間無量綱化.圖8 給出了流固相互作用強(qiáng)度和顆粒碰撞強(qiáng)度隨時間的變化.由于統(tǒng)計了從零時刻開始的變化,而初始狀態(tài)顆粒速度為零且濃度較高,因此可忽略初始時刻附近的統(tǒng)計數(shù)據(jù)(t/tp＜ 2).當(dāng)t/tp=6 時,顆粒碰撞應(yīng)力突然增大,這是因為粗細(xì)顆粒群相交,顆粒群濃度增加,顆粒之間的碰撞頻率也更大.顆粒碰撞應(yīng)力增大過程持續(xù)到t/tp=13,隨后碰撞應(yīng)力逐漸減小并趨于穩(wěn)定值.顆粒碰撞應(yīng)力曲線可以清晰顯示顆粒群混合以及分離的過程.圖將粗細(xì)顆粒群的流固相互作用應(yīng)力分開統(tǒng)計,以清楚顯示混合輸送過程中粗細(xì)顆粒群所受流體作用的區(qū)別.在粗細(xì)顆粒群相交時,細(xì)顆粒群流固相互作用應(yīng)力顯著增加,而粗顆粒群流固相互作用應(yīng)力輕微降低,這是因為粗細(xì)顆粒群混合時粗顆粒群后方由于顆粒擴(kuò)散導(dǎo)致粗顆粒群寬度增加,粗顆粒群后方濃度降低,使此區(qū)間粗顆粒所受流體拖曳力降低,進(jìn)而使整體粗顆粒群的流固相互作用應(yīng)力輕微減小.在整個輸送階段,細(xì)顆粒群的流固相互作用都大于粗顆粒群,使得細(xì)顆粒群運動速度更快,從而出現(xiàn)細(xì)顆粒群越過粗顆粒群分離的現(xiàn)象.此外,可以發(fā)現(xiàn),從粗細(xì)顆粒群混合到完全分離的過程中,粗細(xì)顆粒群的流固相互作用強(qiáng)度之差明顯大于完全分離后的流固相互作用強(qiáng)度之差.這是因為混合區(qū)域流速增大,粗細(xì)顆粒群受到的流體拖曳力都增加,但細(xì)顆粒群的速度提升更大,即粗細(xì)顆粒混合階段細(xì)顆粒群從流體中獲得的能量更大.因此,粗細(xì)顆?；旌虾髸铀賰深w粒群的分離過程.當(dāng)粗細(xì)顆粒群完全分離后,細(xì)顆粒群的流固相互作用強(qiáng)度只略微大于粗顆粒群,從而保證細(xì)顆粒群的速度永遠(yuǎn)略大于粗顆粒群.最后,在本文的數(shù)據(jù)統(tǒng)計范圍內(nèi),顆粒碰撞應(yīng)力顯著大于流固相互作用應(yīng)力,這與常理相悖.因為顆粒群不斷向上運動的驅(qū)動力來源于流體攜帶作用,因此流固相互作用理應(yīng)為主導(dǎo)效應(yīng).但是實際統(tǒng)計數(shù)據(jù)恰恰相反,顆粒輸送過程由碰撞效應(yīng)占主導(dǎo).這是因為算例中顆粒濃度較大,顆粒之間的碰撞頻率與顆粒濃度呈正相關(guān),使得碰撞效應(yīng)超過了流固相互作用效應(yīng).

圖8 流固相互作用強(qiáng)度和顆粒碰撞強(qiáng)度隨時間的變化Fig.8 Changes of the fluid-particle interaction stress and the particle collision stress with time

3.3 影響因素分析

3.3.1 入口流速

入口輸送速度是工程中最受關(guān)注的物理量之一,因為它與輸送能耗息息相關(guān).本小節(jié)討論入口速度對粗細(xì)顆粒群分離過程的影響.圖9 為輸送速度分別為1.0 m/s,1.5 m/s 以及 2.0 m/s 時顆粒群碰撞應(yīng)力隨時間的變化.與上文相同,顆粒碰撞應(yīng)力增加的持續(xù)時間即為粗細(xì)顆粒群混合分離過程的持續(xù)時間.因此,圖9 用線條標(biāo)識出顆粒碰撞應(yīng)力增大部分,可以看出隨著輸送速度的增加,粗細(xì)顆粒群分離過程持續(xù)時間增加,并且輸送速度從1.5 m/s 增大到2.0 m/s 時,分離過程持續(xù)時間延長了幾乎兩倍.這是因為粗細(xì)顆粒分離產(chǎn)生于顆粒尺寸差異引起的滑移速度差異,而滑移速度不僅僅受顆粒尺寸影響,流體輸送速度同樣是重要的影響因素.Zhang 等[2]發(fā)現(xiàn),當(dāng)輸送速度較大時,顆粒在管道中的積聚明顯減少,就是因為顆粒滑移速度降低,流體的攜帶能力增強(qiáng).圖10 給出了不同輸送速度時,顆粒群軸向速度與輸送速度之比(wp/u0)的時間平均值的變化.隨著入口速度的增加,顆粒群時均速度越來越接近入口速度,即顆粒與流體之間的速度差越來越小.粗細(xì)顆粒分離與粗細(xì)顆粒之間的速度差有關(guān),而本文算例中粗細(xì)顆?；旌想A段碰撞頻率增加,顆粒速度的瞬時值波動較大,單獨統(tǒng)計粗細(xì)顆粒的速度時均值進(jìn)行對比發(fā)現(xiàn)差異不明顯.但是基于前人研究結(jié)論[6],在其他條件不變時,顆?；扑俣入S著顆粒尺寸的增大而增大,因此細(xì)顆粒與粗顆粒之間的速度差在輸送速度增加時減小,從而導(dǎo)致分離過程持續(xù)時間增大.

圖9 不同輸送入口速度時,顆粒群碰撞應(yīng)力隨時間的變化Fig.9 Change of the particle collision stress over time at different inlet velocities

圖10 顆粒群軸向時均速度與入口流速之比隨著入口速度的變化Fig.10 The ratio of the vertical time-averaged velocity of the particles to the inlet velocity varies with the inlet velocity

此外,對比發(fā)現(xiàn),隨著入口速度增大,顆粒之間的碰撞應(yīng)力也逐漸增大,即輸送過程會產(chǎn)生更多的能量損耗.這在工程中是需要避免的,因此,實際水力輸送方案設(shè)計中,從經(jīng)濟(jì)角度考慮,不宜采用過大的輸送速度.

3.3.2 粗細(xì)顆粒尺寸

本文著重研究的就是粗細(xì)顆粒的分離過程,因此粗細(xì)顆粒的尺寸同樣是重要影響因素.為了對比明顯,算例設(shè)置中令粗顆粒尺寸為 10 mm,細(xì)顆粒尺寸取值分別為2 mm,3 mm,4 mm 和 5 mm,其他參數(shù)保持不變.圖11 為細(xì)顆粒取不同尺寸時,顆粒群碰撞應(yīng)力隨時間的變化.首先,隨著細(xì)顆粒尺寸的減少,由于顆粒群初始濃度保持不變,因此顆粒數(shù)量增加.如細(xì)顆粒尺寸為 2 mm 時,細(xì)顆粒數(shù)量為 56000,而細(xì)顆粒尺寸增大到 5 mm 時,顆粒數(shù)量為3500.顆粒數(shù)量的增加意味著顆粒之間的碰撞頻率更高,因此按照式(22)計算得到的顆粒碰撞應(yīng)力增加.隨著顆粒尺寸從5 mm 降低到 2 mm,顆粒碰撞應(yīng)力從2.5 左右增大到 34.5.此外,與上節(jié)分析方法相同,顆粒碰撞應(yīng)力的增大部分即為粗細(xì)顆粒群混合分離的過程.可以看出,細(xì)顆粒尺寸變大后,粗細(xì)顆粒群混合分離過程持續(xù)時間更長.這是因為管道中顆粒與流體的滑移速度隨著顆粒尺寸的增大而增大.Richardson 等[30]曾基于實驗提出,兩相流中顆粒與流體速度差滿足如下關(guān)系

圖11 不同細(xì)顆粒尺寸時,顆粒群碰撞應(yīng)力隨時間的變化Fig.11 Change of the particle collision stress over time at different particle sizes

其中wf為流體軸向速度,wt為單顆粒在靜水中沉降速度,n為與顆粒雷諾數(shù)Rep相關(guān)的函數(shù)[31].從式 (24)也可以顯示出顆粒尺寸對兩相滑移速度的影響規(guī)律.因此,當(dāng)粗細(xì)顆粒尺寸差異逐漸增大時,粗細(xì)顆粒之間的軸向運動速度差異也同樣增大,即細(xì)顆粒群相對于粗顆粒群運動更快,從而分離過程持續(xù)時間也就更短.

圖12 為不同細(xì)顆粒尺寸時,細(xì)顆粒群固液相互作用應(yīng)力隨時間的變化圖,可以看出固液相互作用應(yīng)力受細(xì)顆粒尺寸的影響較小.這是因為固液相互作用驅(qū)動著顆粒群持續(xù)向前運動,顆粒運動呈現(xiàn)先加速后勻速運動的趨勢,因此一段時間后顆粒垂向受力會逐漸趨于平衡,即軸向流體拖曳力與顆粒浮重力相等.

圖12 不同細(xì)顆粒尺寸時,細(xì)顆粒群固液相互作用應(yīng)力隨時間的變化Fig.12 Change of the fluid-particle interaction stress of the fine particles over time at different particle sizes

4 結(jié)論

本文著眼于深海采礦中的關(guān)鍵環(huán)節(jié)——垂直管水力輸送,此環(huán)節(jié)的典型特征是顆粒尺寸大,顆粒級配寬,并且顆粒體積濃度高.這些特征會使輸送過程中管道內(nèi)顆粒局部體積濃度顯著大于輸送入口濃度,很可能造成輸送堵塞,影響輸送安全.近年來,相關(guān)研究發(fā)現(xiàn)寬級配特征下粗細(xì)顆粒存在流動混合分離現(xiàn)象,可能會增大出現(xiàn)栓塞,堵塞的可能性.因此,本文對垂直管中雙尺寸顆粒群的混合及分離過程展開研究,探究雙尺寸顆粒群混合及分離時的濃度變化特征,顆粒運動特征,顆粒群受力特性變化等,以解釋顆粒群的分離機(jī)理.

研究發(fā)現(xiàn),雙尺寸顆粒群混合后會使顆粒群前后出現(xiàn)明顯分層.這是因為顆?；旌虾鬂舛仍黾?流體拖曳力增大,混合區(qū)域的運動速度也隨之增大,逐漸與后方顆粒群脫離.而混合區(qū)域內(nèi)兩種尺寸的顆粒運動速度存在差異,細(xì)顆粒運動速度比粗顆?？?會逐漸越過粗顆粒運動到前方,隨后粗細(xì)顆粒徹底分離.雙尺寸顆粒群從初始混合到徹底分離的過程中顆粒之間的碰撞頻率會增加,從而使顆粒群整體碰撞應(yīng)力增加,因此可通過碰撞應(yīng)力增加及降低點來確定顆粒群混合及分離的時刻.對比發(fā)現(xiàn),入口輸送速度增大時,雙尺寸顆粒群從初始混合到徹底分離的時間更長,這是因為入口輸送速度增大后,流體攜帶能力增強(qiáng),顆粒速度越來越接近流體速度,兩種尺寸之間的顆粒速度差異也更小,從而使顆粒分離變慢.此外,研究發(fā)現(xiàn)兩種尺寸顆粒之間尺寸差異越大,顆粒群越容易分離.主要是因為顆粒尺寸對滑移速度影響較大,尺寸差異越大,細(xì)顆粒相比于粗顆粒運動就越快.

最后,綜合研究結(jié)果,本文提出流固相互作用的差異是引起雙尺寸顆粒群分離的根本因素,在混合分離階段,細(xì)顆粒的流固相互作用應(yīng)力遠(yuǎn)大于粗顆粒,使得細(xì)顆粒群能夠快速越過粗顆粒群.當(dāng)完全分離后,細(xì)顆粒流固相互作用應(yīng)力略大于粗顆粒,因此粗細(xì)顆粒群之間的距離會越來越遠(yuǎn).