磁場對濕顆粒流化床系統(tǒng)中介尺度結構影響機制研究

唐天琪，何玉榮

（1 哈爾濱工業(yè)大學能源科學與工程學院，黑龍江 哈爾濱 150001； 2 黑龍江省新型儲能材料與儲能過程研究重點實驗室，黑龍江 哈爾濱 150001）

引 言

流態(tài)化技術由于具有傳熱傳質(zhì)效率高、反應器內(nèi)溫度均勻等優(yōu)點，在化工行業(yè)、太陽能轉(zhuǎn)化與利用以及清潔能源轉(zhuǎn)化等各個方面有著廣泛的應用[1-2]。其中，反應器內(nèi)部氣固兩相流動在氣固兩相之間非線性相間曳力和固相應力、不同反應器結構和操作條件的影響下，系統(tǒng)內(nèi)很容易產(chǎn)生非均勻的介尺度流動結構，如顆粒聚團、氣泡等[3-4]。這些介尺度結構的存在對于反應器內(nèi)的流動及傳熱傳質(zhì)特性有著明顯的影響。因此，研究者們針對介尺度結構對氣固兩相流動特性[5]、傳熱特性[6]影響展開了相關研究。其中，對介尺度結構演化過程規(guī)律的探索、優(yōu)化和調(diào)控成為了研究的重點。

含濕氣固兩相系統(tǒng)作為實際化工領域常見的現(xiàn)象之一，液體的存在導致了氣固兩相系統(tǒng)中流動結構發(fā)生明顯的變化，從而導致濕顆粒系統(tǒng)中的介尺度結構與非含濕氣固兩相系統(tǒng)有著明顯差異。例如，Boyce 等[7]利用CFD-DEM（computational fluid dynamics-discrete element method）方法研究了濕顆粒在流化床中的流動特性，并通過特征數(shù)表征了不同階段顆粒聚團流動行為。Zhao 等[8]應用CFDDEM 方法探究了埋管分布、黏結力以及流化氣速等因素對流化床中濕顆粒系統(tǒng)流動特性的影響，并指明了流化床中容易產(chǎn)生氣體通道和顆粒聚團的區(qū)域。Song 等[9]應用CFD-DEM 方法預測了不同碰撞速度、液體含量、液體黏度和表面張力等因素影響下的顆粒碰撞和黏結行為，研究結果表明當液體黏結作用進一步提高時，濕顆粒的內(nèi)循環(huán)受到抑制，氣體通道變得更加穩(wěn)定。Liu 等[10]、Balakin 等[11]采用DEM 方法分別分析了單顆粒碰撞、顆粒聚團碰撞特性和顆粒聚團流動行為。Cheng 等[12]結合粒子圖像測速技術并基于“形成-破碎”競爭機制，提出了濕顆粒聚團演化模型來預測不同條件下顆粒聚團的演化特征。Wang 等[13]通過實驗測量的方法探究了濕顆粒流化床中的氣泡動力學，結果表明液橋力和曳力的相互影響機制是改變濕顆粒系統(tǒng)最小流化速度的主要原因。由此可見，液體的存在導致氣固兩相系統(tǒng)形成顆粒聚團和氣泡等介尺度結構、宏觀氣固兩相流流體動力學特性發(fā)生改變。因此，對濕顆粒系統(tǒng)中介尺度結構演化規(guī)律進行分析和研究具有重要的研究意義。

目前，國內(nèi)外研究者通過添加外加物理場，如電場、聲場、磁場等，對流化床內(nèi)的氣固兩相流動特性進行調(diào)控和研究。如van Willigen 等[14-15]分析了顆粒對電場的定性響應特性，探究了電場減小流化床中氣泡尺寸的原因和機理。Zhu 等[16]研究了聲場對納米顆粒聚團的流化狀態(tài)、最小流化速度、床層壓降和床層膨脹等流化特性的影響。此外，由于磁流化床在物料篩分等實際工程領域的廣泛應用，磁場對流化床中氣固兩相流動特性的影響同樣得到了研究者的關注。Zhu 等[17-18]對Geldart-B 可磁化顆粒及可磁化顆粒和不可磁化顆粒的雙組分系統(tǒng)在磁流化床中的流體動力學行為進行了研究，分析了磁場改善流化質(zhì)量的原因，研究結果表明隨著磁場強度的增加氣泡尺寸減小，但氣泡生成頻率增加，流化質(zhì)量逐步提高。李響[19]、楊慧等[20]分別通過數(shù)值模擬的方法研究了磁場自身特性對宏觀流動參數(shù)的影響，以及磁場強度對氣固流化床內(nèi)氣含率、壓力波動均方根等參數(shù)的影響規(guī)律。此外，磁場的特性對流態(tài)化特性也有著明顯影響。例如，Espin 等[21]研究了磁場強度、磁場方向以及流化氣速等因素對流態(tài)化的影響，研究指出當磁場方向與氣流方向一致時，流態(tài)化的穩(wěn)定性達到最佳狀態(tài)。毛志等[22]針對二維潰壩問題進行DEM 數(shù)值模擬，分別研究磁場和液體對顆粒行為的影響。由此可見，磁場對于氣固兩相流動特性以及其內(nèi)部的聚團、氣泡等介尺度結構有著明顯的影響，而目前對于磁場對濕顆粒系統(tǒng)中介尺度結構影響機制仍有待進一步研究。

本文采用DEM 方法對鼓泡流化床中典型的介尺度結構演化過程展開研究，并分析外加勻強磁場對含濕顆粒系統(tǒng)中介尺度結構演化過程的影響。首先，在沒有外部磁場的情況下，對干、濕顆粒系統(tǒng)中氣泡的演化過程進行分析。在此基礎上，引入外加勻強磁場，分析濕顆粒系統(tǒng)中介尺度結構在磁場作用下的演化機制。通過分析氣泡受力、氣泡邊緣處顆粒接觸力、磁場力以及曳力，獲得不同磁場對濕顆粒系統(tǒng)中氣泡結構演化過程的主要影響因素，為后續(xù)介尺度結構的演化和調(diào)控提供理論基礎。

1 數(shù)學模型

1.1 氣相控制方程

本文中采用CFD-DEM 數(shù)值模擬方法研究磁場對濕顆粒流化床內(nèi)介尺度影響機制。在CFD-DEM方法當中，流體相的流動過程通過流體動力學方法進行計算，計算過程中視為連續(xù)相，質(zhì)量守恒方程和動量守恒方程基于體積平均Navier-Stokes 偏微分方程組進行數(shù)值求解。氣相運動控制方程為[23]：

如果氣體為牛頓流體，則黏性應力張量τg為：

單位張量I為：

1.2 顆?？刂品匠?/h3>

DEM 方法中，顆粒的運動通過牛頓第二運動定律進行求解，無外加物理場情況下，顆粒受到顆粒之間或顆粒與壁面之間的碰撞作用、氣體對顆粒的作用以及顆粒自身受到重力的作用。在外加磁場下，顆粒受到磁場力的作用；當系統(tǒng)中存在液體時，顆粒間受到液橋力作用，如式（5）所示[24]：

顆粒在轉(zhuǎn)矩的驅(qū)動下進行旋轉(zhuǎn)，轉(zhuǎn)矩Tp包含由接觸力和液橋力在切向方向上的分量以及滾動摩擦產(chǎn)生的轉(zhuǎn)矩，顆粒的旋轉(zhuǎn)速度ωp為[24]：

在本文中，采用的軟球碰撞模型，允許顆粒碰撞瞬間產(chǎn)生輕微形變，而顆粒發(fā)生形變導致顆粒旋轉(zhuǎn)過程中受到阻力，顆粒的旋轉(zhuǎn)運動發(fā)生改變。該部分通過滾動摩擦轉(zhuǎn)矩Tr表示[25]：

1.2.1 接觸力 如式（5）所示，顆粒合外力項包括氣體對顆粒的作用、顆粒碰撞產(chǎn)生的接觸力、液體對顆粒的液橋力以及磁場力等作用。在本文中，顆粒間和顆粒與壁面間碰撞力通過軟球碰撞模型進行求解，該模型定量描述了顆粒碰撞過程的動態(tài)運動過程和能量損耗，并分解為切向和法向兩個方向[26]：

1.2.2 液橋力 對于濕顆粒系統(tǒng)，需要考慮液體對顆粒運動特性的影響。目前通過建立液橋力模型描述該過程。液橋力包括毛細力和黏性力兩部分，根據(jù)相對液體量的大小，毛細力和黏性力對顆粒流動特性主導作用有所不同：當液體含量較小時，毛細力作用更為明顯，當液體含量較高或液體黏度較大時，黏性力對顆粒的流動行為占主導作用。本文中，假設液體均勻包覆在顆粒表面，顆粒碰撞過程產(chǎn)生液橋結構，液體對顆粒運動的影響通過液橋力進行反饋。

其中，毛細力是有關于液體物理性質(zhì)的函數(shù)，顆粒與壁面間法向靜態(tài)毛細力如式（10）所示[27]：

與濕顆粒/壁面間相互作用相似，濕顆粒間的法向毛細力Fcp,n如式（11）所示：

基于Lambert 等[28]和Liu 等[29]的工作，黏性力同樣從切向和法向兩個方向進行求解[30-31]，其中切向方向上黏性力表達式為：

法向黏性力[32]表達式為：

1.2.3 磁場力 通常情況下，磁場力包括磁場梯度力和粒子間磁場作用力兩部分。其原理如圖1所示。

當顆粒被磁化后，顆粒未發(fā)生直接碰撞時，顆粒間受到磁場力的作用，并分為徑向（Fr）和切向（Fθ）兩個方向，如圖1（a）所示；當顆粒間未發(fā)生碰撞但小于液橋臨界斷裂距離時，顆粒受到磁場力、液橋力以及曳力的作用，如圖1（b）所示；當顆粒間發(fā)生直接碰撞產(chǎn)生形變時，顆粒受到接觸力、液橋力、磁場力以及曳力的作用，如圖1（c）所示。

圖1 磁場作用下含濕磁性顆粒受力情況示意圖Fig.1 Schematic diagram of forces imposing on wet particles under effect of magnetic field

磁場力模型的應用基于以下假設：（1）顆粒為球形；（2）顆粒之間間距不大于兩倍顆粒直徑之和；（3）磁性顆粒適用于軟磁性材料的顆粒；（4）顆粒為理想偶極子；（5）忽略顆粒之間的靜電效應。磁化后的鐵磁顆粒會產(chǎn)生磁感應場，其不僅受到外磁場力的作用，還受到臨近顆粒的作用力。如果相鄰磁性顆粒相對位置與磁場方向一致時，顆粒間磁感應力表現(xiàn)為引力；如果相鄰磁性顆粒相對位置與磁場方向垂直時則表現(xiàn)為排斥力[19]。顆粒間磁場力及磁場梯度力計算公式為[33]：

1.3 氣固之間動量交換

氣相運動通過體積平均N-S 方程進行求解，氣體動量方程無法直接閉合，氣固之間動量交換需要通過建立相應的動量交換方程進行求解，相間動量交換速率是與氣固兩相之間速度差、空隙率等因素相關的函數(shù)，F(xiàn)gp如式（24）所示：

其中，動量交換系數(shù)β對于準確描述流化床內(nèi)氣固兩相之間動量交換過程是十分重要的。但直接從理論推導得到其數(shù)值存在一定難度，目前研究者們通過實驗以及直接數(shù)值模擬等方法建立了不同曳力系數(shù)半經(jīng)驗公式。本文中選用Beetstra等[34]提出的曳力模型，目前已廣泛應用于流化床內(nèi)氣固兩相流動特性的計算中[35-36]。求解方程如下所示：

在所有曳力模型中，顆粒Reynolds數(shù)的定義為：

2 邊界條件及模擬工況

本文中選用的磁場力模型、液橋力模型均在前期的研究中得到了驗證[37-38]。因此，本文應用建立的數(shù)學模型對Collier 等[39]研究中的鼓泡流化床中氣固兩相流動特性進行研究，具體結構及尺寸如表1和圖2 所示。底部為布風板，作為流化氣體入口。對于氣相，頂部和壁面分別采用壓力出口和無滑移邊界條件。對于固相，模擬中顆?？倲?shù)為16000個，顆粒平均直徑為3 mm，假設顆粒含有軟磁性物質(zhì)并可以被均勻磁化[40]，是理想的可被磁化的球形顆粒，顆粒磁化率、密度相同。模擬過程中，所有工況模擬時間持續(xù)15 s，模擬參數(shù)如表1 所示。本文中僅考慮勻強磁場對含濕磁性顆粒系統(tǒng)中介尺度結構演化過程的影響，磁場強度分別為0.005、0.010 和0.020 T，濕顆粒系統(tǒng)中液體含量為0.01%。

表1 模擬參數(shù)設置Table 1 Parameters used in the simulation

圖2 流化床結構示意圖Fig.2 Structure of the bubbling fluidized bed

3 結果與討論

3.1 無磁場作用下濕顆粒系統(tǒng)介尺度結構演化特性研究

圖3中給出了無磁場作用下一個氣泡生長周期內(nèi)干顆粒系統(tǒng)和濕顆粒系統(tǒng)中顆?？臻g分布情況。從圖中可以發(fā)現(xiàn)，在干顆粒系統(tǒng)中流化床中形成明顯的氣泡結構，氣泡結構邊界較為清晰（氣泡1），氣泡內(nèi)部攜帶顆粒向上運動。隨著氣泡向上運動，氣泡到達床層表面破裂，同時新的氣泡（氣泡2）在床層底部生成，具體演化結構如圖3（a）所示。與干顆粒相比，濕顆粒系統(tǒng)中形成的氣泡結構內(nèi)部存在明顯的聚團結構，氣泡結構內(nèi)部聚團與氣體進行動量交換，氣泡從床層底部產(chǎn)生，隨著氣泡向上運動，氣泡結構長徑比不斷增大，氣泡邊緣處的聚團結構更不規(guī)則，氣泡逐漸向上運動到達床層表面后推動床層表面形成的顆粒聚團向兩側(cè)運動，如圖3（b）所示。

圖3 無磁場作用下一個氣泡生長周期內(nèi)顆?？臻g分布情況Fig.3 Instantaneous particle distribution and bubble forming process without magnetic field in dry and wet particle systems (a cycle)

圖4中給出了無磁場作用下干顆粒系統(tǒng)和濕顆粒系統(tǒng)中壓降脈動情況以及對應氣泡結構空間分布情況。其中，圖4（b）中氣泡結構通過空隙率進行識別。根據(jù)不同時刻空隙率云圖分布情況，當空隙率位于為[0.8，1.0]區(qū)間時，認為該區(qū)域為氣泡結構。根據(jù)該方法，將不同瞬時氣泡結構進行圖像疊加繪制得到不同工況下不同時刻氣泡結構，直觀追蹤氣泡生長過程。由圖4（a）中的壓降分布情況可以發(fā)現(xiàn)，與干顆粒相比濕顆粒系統(tǒng)中床層壓降脈動更為劇烈，壓降脈動幅度更大。這一現(xiàn)象說明，針對本文中的濕顆粒工況，形成的濕顆粒聚團結構增加了床層的阻力，但床層內(nèi)形成了完整的氣泡結構，因此床層的脈動與干顆粒系統(tǒng)相比較為劇烈，氣固兩相之間的動量交換更為劇烈。同時，氣泡內(nèi)部的聚團結構增大了床層壓降的不穩(wěn)定性。如圖4（b）所示，在干顆粒系統(tǒng)中，氣泡寬度隨著氣泡上升過程逐漸增大，當氣泡運動到床層表面時，氣泡寬度接近床層寬度，氣泡結構充分發(fā)展，在床層表面破裂。而在濕顆粒系統(tǒng)中氣泡結構長徑比明顯大于干顆粒系統(tǒng)中氣泡長徑比，并沿著床層中線向上運動。同時在濕顆粒系統(tǒng)中由于受到液橋力的限制，氣泡兩側(cè)受到擠壓，氣泡寬度無法生長到床層寬度，因此氣泡以更高的長徑比結構突破床層阻力向上運動。

圖4 無磁場作用下床層壓降脈動及一個氣泡生長周期內(nèi)氣泡生長情況Fig.4 Distribution of pressure drop and bubble evolution process without magnetic field (a cycle)

圖5中給出了無磁場作用下干顆粒與濕顆粒系統(tǒng)中氣泡顆粒溫度隨著顆粒體積分數(shù)變化情況，其中氣泡顆粒溫度與氣泡速度脈動的二階矩相關[41]：

圖5 干、濕顆粒系統(tǒng)中氣泡顆粒溫度分布情況Fig.5 Distribution of bubble granular temperature in dry and wet particle systems

從圖中可以發(fā)現(xiàn)，濕顆粒系統(tǒng)氣泡顆粒溫度峰值略高于干顆粒系統(tǒng)中氣泡顆粒溫度峰值，位于顆粒體積分數(shù)[0.25，0.35]對應區(qū)域，而該區(qū)域?qū)獨馀葸吘壩恢?，這說明在氣泡邊緣位置由于受到液體的限制，氣泡周圍受到阻力增強，導致該區(qū)域的脈動增強，同時也導致氣泡結構與干顆粒系統(tǒng)中氣泡結構有明顯不同，氣泡長寬比更大。而在空隙率[0.8，1.0]區(qū)域，即對應氣泡內(nèi)部區(qū)域，氣泡顆粒溫度與干顆粒系統(tǒng)較為接近，但分布更不規(guī)則。這是由于氣泡內(nèi)部存在更多分散的“小聚團”結構，導致氣泡內(nèi)部結構的速度脈動更不規(guī)則。

3.2 勻強磁場作用下濕顆粒系統(tǒng)介尺度結構演化特性研究

圖6中給出了不同磁場強度下一個氣泡生長周期內(nèi)顆粒空間分布情況。從圖中可以發(fā)現(xiàn)，當磁場強度為0.005 T 時，如圖6（a）所示，氣泡從床層底部產(chǎn)生，氣泡邊界不規(guī)則且氣泡內(nèi)部存在明顯的顆粒聚團，隨著氣泡不斷向床層上部運動，氣泡尺寸逐漸增大，并在床層表面破碎推動塊狀顆粒聚團向兩側(cè)運動。隨著磁場強度增大，顆粒聚團逐漸由塊狀轉(zhuǎn)變?yōu)檠刂鸥芯€方向的鏈條狀。如圖6（b）所示，當磁場強度為0.010 T 時，氣泡尺寸明顯減小，氣泡邊界更不清晰。當磁場強度增強到0.020 T 時，如圖6（c）所示，氣泡結構明顯減小，由于磁場作用增強，氣泡和顆粒運動被限制，少量顆粒以鏈狀顆粒聚團沿著磁感線向上運動，流動結構發(fā)生明顯變化。

圖7中給出了不同磁場強度下時均氣泡顆粒溫度分布情況以及氣泡顆粒溫度隨著顆粒體積分數(shù)變化情況。從圖7（a）中可以發(fā)現(xiàn)，當磁場強度為0.005 T 時，氣泡顆粒溫度存在一個位于床層上部氣泡核心，說明在床層上部氣泡脈動相對較為劇烈。隨著磁場強度增強到0.010 T，氣泡顆粒溫度明顯減小，氣泡顆粒溫度核心幾乎消失，說明氣泡脈動現(xiàn)象減弱。當磁場強度進一步增強到0.020 T時，氣泡脈動幾乎消失，這一現(xiàn)象說明隨著磁場強度增強，氣泡脈動現(xiàn)象受到明顯限制。由圖7（b）可知，隨著磁場強度增強，氣泡顆粒溫度明顯減小，當磁場強度為0.005 T 和0.010 T 時，氣泡顆粒溫度在顆粒體積分數(shù)[0.25，0.35]范圍內(nèi)存在峰值，該區(qū)域?qū)獨馀葸吔缣幬恢?，說明氣泡邊界處氣固之間動量交換劇烈，脈動較為強烈。隨著顆粒體積分數(shù)增大，由于顆粒同時受到液橋力、磁場力的限制，氣泡脈動明顯被削弱。當磁場強度為0.020 T時，在較低顆粒分數(shù)處存在峰值，結合圖6中顆?？臻g分布發(fā)現(xiàn)，該區(qū)域?qū)矊由喜孔杂煽臻g區(qū)域。

圖6 不同磁場強度下一個氣泡生長周期內(nèi)顆?？臻g分布情況Fig.6 Instantaneous particle distribution and bubble forming process under different magnetic field intensities in wet particle system

圖8中給出了不同磁場強度下壓降脈動情況以及對應氣泡空間分布情況。由圖8（a）中的壓降分布情況可以發(fā)現(xiàn)，隨著磁場強度增強，壓降脈動周期逐漸減小，這一現(xiàn)象說明由于受到顆粒間液橋力以及磁場力的限制，氣泡生長運動周期明顯減小。此外，隨著磁場強度的增強，壓降脈動振幅逐漸減小，這一現(xiàn)象說明流化床內(nèi)氣泡形成和演化過程明顯受到削弱。圖8（b）給出了氣泡生長過程中的氣泡結構演化過程，從圖中可以發(fā)現(xiàn)，隨著磁場強度的增強，氣泡尺寸逐漸減小，但氣泡生成頻率增加，在流化床內(nèi)運動空間逐漸受到限制，無法達到床層頂部，氣泡邊界更不規(guī)則。這一現(xiàn)象說明，引入勻強磁場限制了氣泡結構的生成和演化過程。

圖8 不同磁場強度下床層壓降脈動及一個氣泡生長周期內(nèi)氣泡生長情況Fig.8 Distribution of pressure drop and bubble evolution process under different magnetic field intensities

從圖8可知，引入勻強磁場后，流化床內(nèi)氣泡尺寸減小但生成頻率提高。為了進一步說明磁場濕顆粒系統(tǒng)中氣泡演化過程作用機制，圖9 在勻強磁場B=0.005 T 工況基礎上，引入梯度磁場（k=0.005 T/m）。從圖中可以發(fā)現(xiàn)，引入與重力方向相反的梯度磁場后氣泡長徑比增大，通過提高顆粒向上的作用，促進了氣泡的生長。因此，通過施加合適的磁場強度或者磁場梯度可以有效控制濕顆粒系統(tǒng)氣泡演化過程，突破顆粒之間液橋力對顆粒的運動阻力。

圖9 梯度磁場對一個氣泡生長周期內(nèi)氣泡生長情況的影響Fig.9 Effects of magnetic field gradient on bubble evolution process

3.3 勻強磁場作用下濕顆粒系統(tǒng)中氣泡邊界受力競爭機制研究

由于氣泡周圍受到顆粒間接觸力、磁場力、氣固之間交互作用等因素共同的影響，通過將顆粒受到的作用力與重力進行對比，并對顆粒受力進行無量綱化處理，對比了氣泡區(qū)域垂直方向上碰撞力與重力比值（Fcon/Fmg）、外加磁場力與重力比值（Fmag/Fmg）、液橋力與重力比值（Fliq/Fmg）以及曳力（氣固之間作用力）與重力比值（Fdrag/Fmg）分布情況，如圖10所示。當磁場強度為0.005 T時，流化床內(nèi)存在一個正在上升中的氣泡和底部剛剛產(chǎn)生的氣泡，底部產(chǎn)生的氣泡初始階段氣泡尺寸較小，氣泡周圍氣固之間的動量交換明顯高于顆粒碰撞、顆粒間液橋力和磁場力產(chǎn)生的運動阻力，此時液橋力與磁場力大小幾乎相同；當氣泡上升到床層中間高度時，氣泡的寬度和高度明顯增大，此時氣泡受到的氣固之間的阻力值略有減小，顆粒之間碰撞產(chǎn)生的阻力、顆粒間液橋力略有增大，此外位于較高位置處的氣泡發(fā)生破裂；隨著氣泡繼續(xù)向上運動，氣泡長徑比進一步增大，氣泡邊緣處顆粒碰撞過程產(chǎn)生的阻力進一步增強，氣泡整體受到的運動阻力隨著氣泡逐漸運動到床層表面處逐漸減小，在該工況下可以看出磁場力未能完全克服顆粒之間液橋形成的阻力。當磁場強度為0.010 T時，床層內(nèi)同樣同時存在兩個氣泡結構，且運動周期更短，這與圖8中壓降脈動分布一致。與磁場強度0.005 T相比，此時磁場對顆粒運動阻力增強，因此進一步阻礙了氣泡的運動，氣固之間的動量交換逐漸減小。由于氣泡內(nèi)部存在明顯的顆粒聚團結構，氣泡內(nèi)部受到的阻力波動增強，此時液橋力明顯增大，這也說明在氣泡內(nèi)部有明顯的聚團結構，液橋力與磁場力共同控制著氣泡的運動。當磁場強度為0.020 T時，一個氣泡演化周期床層內(nèi)僅有一個氣泡存在，根據(jù)顆粒的受力分布情況可以發(fā)現(xiàn)，勻強磁場對顆粒運動阻力增強，從而對氣泡的運動阻力也進一步增大，明顯限制了氣泡的運動，結合顆?？臻g分布情況氣泡內(nèi)部存在更多不規(guī)則的顆粒聚團結構，氣泡內(nèi)部受到的阻力增強，磁場力明顯高于顆粒之間的液橋力，這也說明磁場和液橋力進一步限制了氣泡的演化過程，此時磁場作用對氣泡的限制作用更為明顯。

4 結 論

本文采用CFD-DEM 方法對鼓泡流化床中典型的介尺度結構演化過程展開研究，并分析外加勻強磁場對含濕顆粒系統(tǒng)中介尺度結構演化過程的影響。分析了沒有外部磁場的情況下，干、濕顆粒系統(tǒng)中氣泡的演化過程。在此基礎上，引入了外加勻強磁場，研究了濕顆粒系統(tǒng)中介尺度結構在磁場作用下的演化機制。并通過分析氣泡受力、氣泡邊緣處顆粒接觸力、磁場力以及曳力，獲得了不同磁場強度對濕顆粒系統(tǒng)中氣泡結構演化過程的主要影響因素，獲得了如下主要結論。

（1）隨著磁場強度增大，顆粒聚團逐漸由塊狀轉(zhuǎn)變?yōu)檠刂鸥芯€方向的鏈條狀，氣泡尺寸明顯減小，氣泡邊界更不清晰，氣泡演化過程受到抑制。同時，隨著磁場強度增強，壓降脈動周期逐漸減小，說明由于受到顆粒間液橋力以及磁場力的限制，氣泡生長運動周期明顯減小。

（2）氣泡結構生成初期，氣泡尺寸較小，受到氣泡上部床層重力以及阻力的影響，氣泡受到的阻礙作用較大，隨著氣泡結構不斷生長并沿床層上升，氣泡邊緣處受到的阻力逐漸減小，由于氣泡內(nèi)部在該工況下存在不規(guī)則聚團，氣泡中心區(qū)域受到的阻力高于氣泡邊緣位置。

（3）對比了氣泡區(qū)域接觸力與重力比值（Fcon/Fmg）、外加磁場力與重力比值（Fmag/Fmg）、液橋力與重力比值（Fliq/Fmg）以及曳力（氣固之間作用力）與重力比值（Fdrag/Fmg）分布情況。隨著磁場強度增強，勻強磁場對顆粒運動阻力增強，氣泡內(nèi)部存在更多不規(guī)則的顆粒聚團結構，氣泡內(nèi)部受到的阻力增強，限制了氣泡的演化過程。

本文通過添加勻強磁場，限制了濕顆粒流化床內(nèi)氣泡結構的生成，減小了氣泡的運動周期，加快了氣泡生成頻率，添加適當?shù)奶荻却艌隹梢源龠M氣泡的生成。后續(xù)將進一步開展梯度磁場促進或抑制氣泡結構生長及演化機制研究，為通過磁場調(diào)節(jié)濕顆粒系統(tǒng)傳熱傳質(zhì)特性提供一定理論基礎。

符 號 說 明

B——磁感應強度，T

d——浸沒高度，m

dp——顆粒直徑，m

Fc——接觸力，N

Fc,t——切向接觸力，N

Fgp——氣相與顆粒相間的動量交換率，kg/(m2·s2)

Flb——液橋力，N

Flb,t——切向液橋力，N

Fmag——磁場力，N

Fme——磁場梯度力，N

g——重力加速度，m/s2

H——顆粒/壁面間的分離距離，m；

Ip——轉(zhuǎn)動慣量，kg·m2

k——彈簧剛度，N/m

Mf——用于計算時間平均值的幀數(shù)

M——磁場強度，A/m

?M——磁場強度梯度，A/m2

m——磁矩，A/m2

mp——顆粒質(zhì)量，kg

Np——模擬單元內(nèi)的顆粒數(shù)

p——氣體壓力，Pa

r——顆粒之間相對距離，m

r——顆粒中心位置，m

rp——顆粒半徑，m

Tp——轉(zhuǎn)矩，N·m

Tr——滾動摩擦轉(zhuǎn)矩，N·m

t——時間，s

ug——氣體速度，m/s

Vcell——模擬單元體積，m3

Vp——顆粒體積，m3

vp——顆粒速度，m/s

vr——相對碰撞速度，m/s

β——相間動量交換系數(shù)，kg/(m3·s）

εg——每個計算網(wǎng)格內(nèi)的氣相空隙率

δ——彈性變形量，m

η——阻尼系數(shù)，N·s/m

Θ——氣泡顆粒溫度，m2/s2

θ——接觸角/粒子連線與磁感線間的夾角，rad

μf——滑動摩擦系數(shù)

μg——氣體剪切黏度，Pa·s

μr——滾動摩擦系數(shù)

μ0——真空磁導率，N/A2

ρg——氣體密度，kg/m3

σ——表面張力系數(shù)，N/m

φ——半填充角，rad

χe——顆粒的有效磁化率

ωp——顆粒旋轉(zhuǎn)速度，rad/s

下角標

n——法向

t——切向