孔令菲，陳延佩，王維

（中國科學院過程工程研究所，多相復雜系統(tǒng)國家重點實驗室，北京 100190）

引 言

由于顆粒-顆粒非彈性碰撞和顆粒-氣體之間黏性耗散的存在，氣固流態(tài)化系統(tǒng)是一種典型的非線性非平衡系統(tǒng)。特別是在鼓泡流化、快速流化等流域中，大量的顆粒聚集形成非均勻結構[1]。這些非均勻結構具有較大的時間和空間跨度，其尺寸可以小到顆粒微觀尺寸，大到反應器尺寸，因此也稱為介尺度結構(the mesoscale structure)[1-2]。介尺度結構的出現顯著影響氣固流態(tài)化過程中反應速率、動量傳遞、熱量傳遞以及質量傳遞等宏觀特征[3]。例如，在制藥、礦石焙燒和造粒等工業(yè)應用中，顆粒團聚調控可影響特定粒徑分布的化工產品的生產[4]；在催化裂化、煤氣化、礦石還原以及生物質燃燒等工業(yè)應用中，顆粒團聚產生的介尺度結構若不能有效控制或消除，會導致床層去流化，造成反應器停車和經濟損失，產生不利影響[5]。因此，為了應用顆粒團聚機理實現目標產品可控制備，抑制其不良影響，研究人員需要進一步了解顆粒流態(tài)化介尺度結構形成與演化機理，以實現對團聚過程的預測和調控。然而，由于介尺度結構的非平衡特性，經典的顆粒動力學理論無法準確預測介尺度結構的形成和演化過程，因此，介尺度結構研究成為工業(yè)與學術界的前沿研究問題之一，尤其是本文將關注細顆粒與超細顆粒形成的介尺度結塊結構，總結細顆粒介尺度結構的研究現狀和發(fā)展動態(tài)，以及形成的相應介尺度結構的時間與空間依賴特性特征。重點從結塊的力學成因出發(fā)，分析顆粒之間范德華力對結塊形成的影響，即細顆粒力平衡理論模型，探究顆粒之間相互作用力對介尺度結構動力學影響。

1 非均勻介尺度結構

顆粒團聚物介尺度結構在氣固流態(tài)化過程中主要有兩種存在形式：一類是由流體不穩(wěn)定性或顆粒間非彈性碰撞引起顆粒能量耗散[6]，進而形成的顆粒短暫接觸、局部區(qū)域濃度高的松散顆粒集合，亦稱作團聚(cluster)[7]。團聚多存在于鼓泡流化[2]、快速流化[1]等流域，相互間不斷發(fā)生碰撞、碰撞后反彈、破碎和再次團聚等動態(tài)演化，并不穩(wěn)定。另一類是由于黏性細顆粒或超細納米顆粒(nanometer particles,NPs)間相互作用導致的持久接觸而形成的黏性結塊(agglomerate)。從微觀角度看，結塊是多孔道、分形的穩(wěn)定結構[8-9]。結塊可以自然形成，也會在流態(tài)化過程中由于碰撞黏附形成。結塊一旦形成，會在流化床中穩(wěn)定存在，不會因為流化過程的停止而立即消失。

另外，針對介尺度結構的命名問題，部分文獻[10]將黏性細顆?；蚣{米顆粒形成的穩(wěn)定結塊稱為聚合體(aggregate)，或者是硬結塊(hard agglomerate)[11-12]。近年來，研究人員主要使用結塊[13-14]描述黏性細顆粒與納米顆粒形成的穩(wěn)定結構，而聚合體在文獻[15-16]中是指由強力（如共價鍵、離子鍵、金屬鍵或其他燒結鍵）剛性結合在一起的初級顆粒組合。本綜述主要關注顆粒物質在冷態(tài)氣固流態(tài)化中的介尺度結構成因與動態(tài)演化特征，不涉及化學反應或高溫過程。為了避免歧義，本文在論述中不出現聚合體這一名詞概念，使用團聚和結塊來描述兩類非均勻介尺度結構。

在氣固流態(tài)化系統(tǒng)[17]中，介尺度結構不僅在同一空間位置隨時間交替出現，而且在不同空間位置會同時出現稀密分布，具有時間和空間依賴性[18]。介尺度結構的形成被認為是氣固系統(tǒng)動力學不穩(wěn)定性與非線性演化的根源[19]。目前由于實驗測量與分析能力的限制，研究人員很難捕捉單顆粒水平的運動規(guī)律，因此對介尺度結構的成因及運動規(guī)律的物理本質沒有深入的了解?，F有理論大多描述宏觀現象，基于顆粒層次相互作用力、能量耗散機制和動態(tài)演化特征的機理研究尚未成熟。因此，從時間和空間尺度探究顆粒層次的微觀運動，可以幫助研究人員理解介尺度結構成因與動態(tài)演化的影響，進而預測宏尺度行為。

2 介尺度結構成因

介尺度結構的形成是顆粒間相互作用（黏附、碰撞摩擦等耗散作用）[15,20]與氣相流體力學不穩(wěn)定性[21]共同作用的結果。這些相互作用過程發(fā)生在短時間間隔、很小的空間內。在這種微小的時空范圍內，顆粒間、顆粒與氣相間的相互作用和各種物理因素（比如，操作條件[21]、顆粒表面性質[22-23]、是否引入外力[24-25]等）有關，所以對介尺度結構成因與動態(tài)演化分析是一個十分復雜的課題。本文將其分成流體-顆粒、顆粒-顆粒之間相互作用分別展開說明。

2.1 流體-顆粒相互作用：流體力學不穩(wěn)定性

氣固兩相流中，流體的黏性耗散引起氣固兩相間相對運動[26]。相對運動會釋放出低壓渦流，吸引周圍顆粒遷移，導致局部氣速和氣壓波動，由此產生的黏滯阻尼效應會使顆粒在流體中漲落衰減。在湍流渦旋下，顆粒并非完全跟隨氣體的流線發(fā)展，它們可能螺旋離開或穿過氣體渦流，在流化床形成傾向性分布[27]。當顆粒與流化氣體產生相對運動時，顆粒會受到Saffman 力、Magnus力、Basset力和曳力等作用力[15]。不同尺寸、表面性質、堆積排布方式的顆粒在氣固系統(tǒng)中受到的主要作用力不同，介尺度結構、動態(tài)演化過程以及動力學特征也存在差異。在稀顆粒濃度下(如氣力輸送提升管[6,28]中)，顆粒碰撞概率很低，此時，流體力學不穩(wěn)定性成為顆粒團聚的主要因素。Stokes 數(St)表示為顆粒在流場中運動的弛豫時間與流體特征時間之比，用于衡量顆粒所受流體力學作用，表達式為

式中，U和L分別代表流場的特征速度與特征尺度。對于黏附性強的細顆粒（小St），顆粒慣性較弱，基本跟隨氣相流線運動；對于中等St的情況[St=Ο(1)]，顆粒相對流體運動有較明顯的漂移；對粗顆粒流(通常dp> 1 mm)，一般St?1，此時流體力學不穩(wěn)定性相對顆粒慣性較弱，通?？梢院雎?。在氣固流化床中，相較于其他作用力，氣固相間曳力對顆粒的夾帶和輸送過程影響較大。關于氣固兩相間的曳力計算問題，國內外學者已做了大量的重要的研究工作[29-32]，這里不再贅述。本文將主要關注顆粒-顆粒之間相互作用。

2.2 顆粒-顆粒（結塊）：顆粒間作用力

顆粒間相互作用力一般指由分子間范德華力[15]、靜電力[33]以及液橋力[34]引起的黏附力(cohesion or adhesion)。顆粒間相互作用力不僅是單分散顆粒形成團聚的驅動力，也是團聚物動態(tài)演化[35]的基礎。其中，靜電力是由顆粒間或顆粒與壁柱摩擦促進表面間的電子轉移產生的相互作用。當凈電荷在顆粒上積累足夠大時，靜電力導致顆粒團聚[36]，影響顆粒的流動行為。通常，研究人員采用模擬的方法研究靜電力，例如Pei 等[37]采用DEM-CFD 方法模擬了流態(tài)化過程中接觸帶電和靜電相互作用，分析比較顆粒的電荷累積和動力學行為。研究發(fā)現，當考慮靜電相互作用時，雙電荷顆粒具有不同功函數，形成團聚。當氣體速度增大到足以促使顆粒間發(fā)生摩擦運動時，電荷的積累速度加快，此時靜電作用不能忽略。實驗方面，Royer 等[35]對比導電顆粒與非導電顆粒團聚現象，證明在100 μm玻璃球真空自由落體體系下，靜電力對團聚現象沒有影響?？梢姡o電力與顆粒尺寸密切相關。

而范德華力[38-39]是兩個相距極近顆粒之間偶極子相互作用的近程力，由顆粒表面的分子和原子之間的電力學力組成，包括定向(Keesom)電位、誘導(Debye)電位和色散(London)電位。范德華力的大小與顆粒作用距離、凈勢能，以及顆粒的介電常數等性質有關。值得注意的是，Royer 等[35]類比液體的Plateau-Rayleigh 不穩(wěn)定作用[40-41]，觀察到即使在排除流體力學不穩(wěn)定性影響的真空條件下，自由落體的粒徑為100 μm的干燥顆粒也會發(fā)生團聚，破裂形成簇狀團聚物。為解釋這一現象，Royer等[35]使用高速視頻成像技術和原子力顯微鏡分析表面性質對團聚現象的影響，發(fā)現在該實驗研究體系下，測得的顆粒間黏附力比液體Plateau-Rayleigh 不穩(wěn)定作用的表面張力小約5 個數量級。即便如此，微小的范德華力和毛細力也會導致團聚的形成。

Lamarche 等[42]將范德華力除以單個顆粒重力，得到的無量綱的范德華力，也被稱為黏結數(Bond number,Bo)[43]，用于評價顆粒間范德華力主導的黏附力對顆粒運動行為的影響。從簡單分析角度來看，顆粒的重力與其粒徑的立方成正比；Hamaker 理論解釋說，對于完全光滑球顆粒，范德華力與顆粒粒徑成正比[38]。Bo可以表示為

從式（2）可以分析出，Bo與顆粒粒徑的平方成反比，隨著粒徑減小，Bo以粒徑的平方的數量級增大，范德華力相對于顆粒重力作用效果更加顯著；隨著粒徑增大，Bo急劇下降，此時顆粒運動受力更多來自于氣固曳力、碰撞耗散以及表觀重力。盡管目前研究人員沒有對顆粒粒徑、Bo與受力特征有明顯的定量歸類，仍然可以從大量氣固體系下顆粒物質運動動力學分析實驗[42,44-47]中歸納流化區(qū)發(fā)生轉變與相應的Bo變化有關，如圖1所示。

圖1 無量綱數Bond數及介尺度結構形成主導作用劃分[17]Fig.1 Dimensionless Bond number and the dominant forces in mesoscale structure formation

從Bo出發(fā)，可以將介尺度結構形成分為兩大類：對于單分散非黏性粗顆粒[48]，粒徑較大，此時Bo< 1，范德華力[49]相對于顆粒重力作用效果不顯著；此外，范德華力屬于近程力，但是粗顆粒間相互作用距離遠大于范德華力作用范圍(0.3～0.4 nm)[38]，碰撞耗散效應[28]在這一類顆粒動態(tài)演化中占主導作用。對于黏性細顆粒或納米顆粒，Bo≥1，范德華力主導的黏附力作用顯著增強，此時黏附性細顆粒和納米顆粒呈現不同于粗顆粒流態(tài)化特征，形成穩(wěn)定結塊。隨著流化過程進行，結塊發(fā)生碰撞耗散，形成類似于粗顆粒的流化團聚。Molerus[50]采用黏結數Bo來區(qū)分Geldart A 類顆粒和C 類顆粒[5]?？梢岳斫鉃?，顆粒之間的范德華力產生黏附作用，這一黏附作用對于粗顆粒可以忽略，但對于細顆?；蚣{米顆粒的黏附接觸起著主導作用。因此研究者有必要更深入探究不同Bo下，顆粒的動力學及其介尺度結構演化行為。且目前對顆粒流體兩相流系統(tǒng)相圖的研究，研究人員一般只考慮其重力與曳力因素，尚未考慮顆粒間范德華力，而依賴于Bo的流化相圖也有待進一步展開研究。

2.3 顆粒-顆粒：碰撞耗散作用

當顆粒之間相互作用力可以忽略時，團聚形成[51]物理過程通常描述為，顆粒漲落使其局域數密度增加，縮短了這一局域內顆粒運動的平均自由程，顆粒碰撞概率增加；同時，由于顆粒間非彈性碰撞耗散使其相對運動速度降低。此時，顆粒碰撞的衰減速度大于周圍局域數密度較低區(qū)域內衰減速度，顆粒溫度隨之降低[6,28]；固相壓力與其溫度成正比[52]，由此產生的壓力梯度使周圍高壓區(qū)流場的自由顆粒趨向于流向低壓區(qū)，形成團聚。在Geldart A、B、D 三類顆粒[5,19]流態(tài)化過程中，均有出現團聚。團聚的形成機理以及演化影響因素復雜，因此發(fā)展實驗和計算機模擬研究手段對該類介尺度行為進行深入定量研究具有重要意義。

根據對被測氣固流場有無擾動，現有的測量方式主要分為兩大類：接觸式測量[53]和非接觸式測量[54-55]。接觸式測量使用直接接觸流化介質的探頭來采集信號（例如光纖探頭[56]、壓力探頭[57]等），然后分析處理采集信號。Liu等[56]使用PV6D反射式光纖探頭測量循環(huán)流化床內的局部固相體積分數εc，該參數影響氣固相間動量傳遞；這一研究提出的基于識別因子n的團聚識別準則可用于識別流化床中的介尺度結構?？紤]到團聚具有時間依賴性，而且形狀不規(guī)則（流線形、拱橋形、馬蹄形等），接觸式有限點測量和低頻率采集無法完全表征團聚幾何特征和動態(tài)演化特性；同時，探頭干擾流場的自然流動，給測試結果帶來不可避免的誤差。

近年來，隨著圖像采集技術和圖像處理算法的快速發(fā)展，團聚分析更多采用非接觸式可視化粒子測速技術[58]。該技術對流場流動無干擾，簡單易操作，可實現全場瞬態(tài)測量。隨著研究角度的不斷深入，許多學者使用高速攝像機攝取流化床內顆?；驓馀葸\動圖像序列，對其運動特性進行深入研究，以獲得流化床內更為詳實的氣固兩相流動信息。Yang 等[59]使用高速攝像機記錄了二維矩形CFB 提升管內的氣固兩相流動，以固體含率直方圖中的類間方差最大值作為判斷團聚的閾值。McMillan等[21,28]采用PTV技術，統(tǒng)計分析A類FCC催化劑顆粒與B類玻璃球顆粒在流化床中顆粒運動速度和顆粒濃度分布，討論了不同工況下團聚形成機制，結果表明，顆粒間范德華力在鼓泡床中對團聚起主要作用，而顆粒碰撞耗散與氣固曳力在湍動床中占主導地位；Mondal 等[60-61]重點分析區(qū)別不同粒徑玻璃球顆粒（粒徑范圍256～440 μm）團聚與單顆粒的灰度閾值分割方法，定量分析團聚在空間尺度從毫米級單顆粒到床層直徑的數量級變化，刻畫介尺度團聚的空間依賴特征。

以上實驗和模擬都是將團聚作為一個整體來研究，缺乏從單顆粒層次定量統(tǒng)計團聚演化的數據。目前采用時間分辨的分析方法(time-resolved analysis)[62]，使用高速攝像機準確檢測擬二維流化床中D 類顆粒的單顆粒的位置坐標和速度，以及湍動流化形成的介尺度團聚。本研究使用PTV，獲得單顆粒級別的物理量隨時間演化規(guī)律，例如顆粒位置、速度、動量流率、平均加速度和平均動能（Voronoi 單元面積平均），進一步量化團聚的動力學特征，發(fā)現能量損失模型(?E∝t32)，從而獲得了單顆粒層次上團聚動態(tài)行為數據。這一實驗研究將有助于后續(xù)對團聚隨時間演化開展更深入模擬建模工作。

除實驗外，借助計算機模擬手段，研究人員的關注點從團聚的穩(wěn)態(tài)逐漸轉變到動態(tài)演化[17,62-63]上。Liu 等[64]對A 類顆粒團聚演化進行CFD-DEM 模擬，重現重力驅動的A 類顆粒團聚過程。模擬結果表明，團聚局域數密度增加會提高顆粒運動相對速度，顆粒間碰撞概率增大，碰撞后能量耗散，團聚會破碎分裂成更小的團聚或分散顆粒，團聚程度降低。Liu 等[64]將這一過程總結為團聚誘導解聚(cluster-induced deagglomeration)效應，并提出“團聚飽和(saturation of agglomeration)”的概念，該理論可用于發(fā)展團聚演化的動力學建模。然而該模型分析氣固間曳力仍建立在平均曳力基礎上，細顆粒循環(huán)流態(tài)化模擬與實際偏差較大，捕捉介尺度團聚仍有難度，難以實現對真實的流化床進行模擬。Li等[65]構建了考慮顆粒團聚影響的非均勻動態(tài)多尺度氣固曳力模型(dynamic cluster structure dependent drag model, DCSD drag model)。該模型考慮顆粒間相互作用力的時空演化特性，構造動量守恒方程，模擬描述了提升管內顆粒團聚的生成、破碎，捕捉到了提升管內環(huán)-核流動結構特征，較好地預測了提升管內顆粒的非均勻分布。Jiang 等[66]繼續(xù)對DCSD 模型進行擴展研究。模型考慮顆粒間非彈性碰撞和流體黏性引起的能量耗散，通過雙流體模型耦合顆粒動理論，依據最小能量耗散率識別團聚現象，討論了提升管入口脈動氣速頻率和幅值對團聚的尺寸與弛豫時間的調節(jié)和控制作用。完善的顆粒動力學理論對于顆粒流體系統(tǒng)的多尺度建模起著至關重要的作用。顆粒間相互作用力使顆粒物質在氣固系統(tǒng)中產生豐富的動力學行為，以下將重點討論范德華力主導的介尺度結構。

3 范德華力主導的介尺度結構

黏性細顆粒（又稱Geldart C 類顆粒[5]）與超細納米顆粒(NPs)應用的迅速發(fā)展，例如，Ni/Al2O3細顆粒作催化材料合成天然氣[67]，金屬-有機骨架(MOF)作為新型微孔材料用作固體吸附劑[68]，MgO顆粒形成高表面積支撐結構促進碳納米管生長[69]等，氣固流態(tài)化的研究熱點從粗顆粒輻射到更細小尺寸顆粒。

由于黏性細顆粒與納米顆粒難以流態(tài)化，目前應用大多采用[70-72]振動或微射流等輔助方法來改善其流態(tài)化特性，以實現高效運輸與混合。但是，細顆粒與納米顆粒氣固體系內復雜作用力以及介尺度結構的動態(tài)演化動力學及過程調控尚未深入研究。

3.1 黏性細顆粒流態(tài)化

與單分散無黏性粗顆粒不同，黏性細顆粒隨著粒徑減?。叫∮?0 μm[5]），在氣固流態(tài)化過程中所受范德華力愈加顯著，同時黏性細顆粒表面能大，接觸剛度小，碰撞時形變程度大，需要更多的能量打破顆粒間的接觸。在氣流作用下，細顆粒極易相互黏附，形成具有分形結構的多孔道穩(wěn)定結塊。經典Hamaker 理論[73]假設表面無形變兩光滑剛性球顆粒，接觸距離在范德華力作用范圍內（0.2～0.4 nm）時的范德華力對應兩顆粒組成分子之間的范德華力的積分累積[74]，且只計入范德華力引力部分（Lennard-Jones 關系式），通過積分疊加計算兩剛性球顆粒的相互作用力。

式中，A是Hamaker 常數，與材料的介電常數和折射率有關；R為球形顆粒粒徑；D是兩球之間最近點距離，即作用距離。Hamaker 理論未考慮表面粗糙度等表面性質對兩顆粒接觸狀態(tài)下發(fā)生接觸形變的范德華力的影響。發(fā)生接觸形變的范德華力的大小不是簡單分子間范德華力的積分[75]，在計算時需要引入黏附接觸模型（如JKR[75]、DMT[76]模型等），使用壓縮表面的表面能確定黏附力。顆粒的團聚和破碎主要通過動態(tài)碰撞過程完成，黏附接觸理論描述了靜態(tài)或準靜態(tài)的物理作用。同時，范德華力主導的顆粒間黏附力的存在為系統(tǒng)引入新的能量耗散機制。系統(tǒng)考慮條件不同，黏附模型選擇也有相應的差異。在實驗中，難以確定黏附接觸過程中考慮的極性相互作用和接觸形變的關鍵物理量，黏附接觸模型適用范圍較窄。

在實際生活與工業(yè)應用中，顆粒通常不是理想光滑的，而是具有一定粗糙度，粗糙顆粒表面范德華力顯著減弱[40]。因此，研究人員開始關注顆粒表面粗糙度對范德華力的影響[77]。Lamarche 等[47]考慮A 類顆粒表面粗糙度，對Hamaker 范德華力預測公式進行修正，解釋兩個相交粗糙表面顆粒的相互作用，并提出預測新理論，如圖2所示。該理論預測不規(guī)則形狀或多粗糙度接觸的顆粒間黏附力，使用原子力顯微鏡以及快速傅里葉變換濾波得到相應大、小尺度表面粗糙度的定量數據，即粗糙度波長λ（相鄰的峰-峰之間的距離）和粗糙度的均方根高度RMS。對于不同粗糙度的兩個表面，范德華力主導的顆粒間黏附力的計算公式如式(4)所示。

圖2 粗糙顆粒接觸作用力新理論[47]Fig.2 A new theory of adhesive force for rough particles[47]

式中，A代表Hamaker 常數；D是指兩個粗糙體分離距離，通常取兩顆粒間范德華力作用距離，即D=D0= 0.3 nm；λ代表原始模型中定義的峰間距離；r指顆粒粗糙度半徑；y指粗糙度高度，下角標L、S分別代表濾波后大尺度和小尺度。

式（4）等號右邊第1 項指兩個顆粒光滑部分的黏附力；第2、3 項指顆粒i與顆粒j的大小凸起之間的黏附力；第4、5 項指顆粒j與顆粒i的大小凸起之間的黏附力；第6、7 項指顆粒i的大凸起與顆粒j大小凸起之間的黏附力；第8、9項指顆粒i的小凸起與顆粒j大小凸起之間的黏附力。

考慮到顆粒間相互作用黏附力，細顆粒團聚或結塊碰撞動力學比粗顆粒的碰撞動力學更復雜。Durhuus 等[16]提出基于Langevin 動力學算法[78]的磁性納米顆粒相互作用模型，該模型用以模擬膠體懸浮液的時間演化和聚集。模擬結果表明，促使形成顆粒團聚的范德華力的取向性會影響團聚過程。例如，各向同性顆粒間范德華力最大限度地增加了近距離表面的數量，有利于結塊緊密堆積，而各向異性顆粒間范德華力則促進團聚向線性結構發(fā)展。Liu 等[79]將考慮表面粗糙度預測的黏附力結果代入CFD-DEM 模擬中，與去流化實驗結果進行對比，引入了表面粗糙度的范德華力預測模型，修正后的Hamaker 理論可以更好地復現顆粒的宏觀流動行為。本研究構建了一個從單顆粒微觀尺度的受力分析到顆粒宏觀流動行為預測的橋梁，使得實驗與相應模擬工作直接驗證成為可能。然而這一實驗與模擬體系仍然基于穩(wěn)態(tài)下的宏觀流態(tài)化物理量的對比，反映介尺度結構時間依賴性的動態(tài)演化的分析模型仍然缺乏。

3.2 超細納米顆粒流態(tài)化

超細納米顆粒的表面結構和晶格結構發(fā)生變化，具有豐富的物理化學性質。在液固體系中，超細納米顆?？梢酝ㄟ^液固流態(tài)化，由電解質和溶解聚合物的空間位阻形成的雙層結構有效削弱范德華力的作用，阻止NPs 聚集[80]。與液固體系相比，這種分離機制在氣固體系中并不明顯，同時決定范德華力大小的Hamaker常數通常比液固系統(tǒng)大，因此，NPs在儲存或流態(tài)化過程中會出現大而致密的結塊結構，難以均勻流態(tài)化。

納米顆粒通過兩種不同黏聚機理，經三步最終形成復雜結塊結構，也稱為多級團聚體(multi-stage agglomerates,MSA)結構。Wang 等[14]研究了6 種高孔隙率、尺寸介于7～16 nm 的SiO2粉體的流態(tài)化過程，詳細分析多級團聚體結構對流態(tài)化特性的影響。一級MSA結構是粒徑小于20 nm的非多孔球體聚集形成的具有極低堆積密度的三維網狀骨架鏈狀團聚(chain clusters)。二級MSA 結構是三維網狀結構合并形成更大尺寸球形或橢球形聚集物(conglomeration)， 又 稱 簡 單 結 塊 (simple agglomerates)，尺寸分布在1～100 μm。簡單結塊體積小，表面密度低，表面粗糙，易黏附；連接的簡單結塊之間存在大量的空隙，這是納米粉體堆積密度低的另一個重要原因。流態(tài)化前后的簡單結塊無明顯差異，其形成與氣固作用無關，僅與顆粒表面性質與顆粒間作用力有關。流態(tài)化開始，在氣固相間曳力作用與簡單結塊之間碰撞與黏附力的作用下，簡單結塊進一步結合形成三級MSA 結構，又稱復雜結塊(complex agglomerates)。在流態(tài)化過程中，復雜結塊是簡單結塊的合并與破碎的動態(tài)平衡狀態(tài)的結果，因此結塊流態(tài)化過程可以看作是簡單結塊的自我重排過程。上述的三維骨架聚集物、簡單結塊，發(fā)展至復雜結塊，即為納米顆粒流態(tài)化過程中的三級MSA 結構。Durhuus 等[16]考慮磁性納米膠體顆粒的團聚動力學，將80個磁性納米顆粒分散在流體中，使用Langevin 動力學方法描述偶極相互作用、范德華力和布朗運動的競爭效應對聚集動力學以及結構特征的影響。結果表明，主導作用力不同，相應形成的團聚體結構特征（線性鏈狀、環(huán)狀、三角形或四面體形等）也不同：在以范德華力占據優(yōu)勢的體系中，顆粒具有較高的團聚矩，形成的長鏈結構進一步團聚成結塊；而更大的團聚物主要通過布朗運動生長，這是一個比鏈生長更漫長的過程。在受力復雜的細顆粒與納米顆粒流態(tài)化體系中，分析作用力對結塊形成、動態(tài)演化以及系統(tǒng)宏觀流動行為影響對深入探究介尺度結構有重要意義。

3.3 力平衡模型研究進展

在細顆?；蚣{米顆粒氣固流態(tài)化系統(tǒng)中，結塊尺寸是反映流態(tài)化性質的重要參數之一。不同類型的力作用于單顆粒及結塊，促使結塊不斷碰撞、合并、破碎，呈現一個典型的時間依賴與空間依賴的結塊動態(tài)演化過程。上述作用力可以分為兩大類：黏附作用力(adhesion)和分離作用力(separation)。通常認為黏附作用力主要來源于顆粒間強相互作用的范德華力與顆粒表觀重力，顆粒黏附聚集，形成結塊；分離作用力使結塊從外部分離破裂，主要包括氣固曳力與碰撞力。如果黏附作用力大于分離作用力，結塊仍能黏附更多離散顆粒，不斷發(fā)展長大；反之，結塊發(fā)生碰撞后破碎成更小的碎片。當分離作用力和黏附作用力處于動態(tài)平衡時，此時流化床內結塊尺寸將會穩(wěn)定。

對于穩(wěn)態(tài)結塊，研究人員提出不同的力平衡模型來計算結塊尺寸。Chaouki 等[81]將結塊受到的黏附力簡化為原生顆粒單點接觸的范德華力。事實上，在實際流化過程中，結塊間是多點接觸，且結塊尺寸比原生顆粒尺寸約大5 個數量級，這一模型低估了結塊間黏附力，僅適用于低密度材料。為了估算細顆粒氣固流態(tài)化中結塊的尺寸，Zhou 等[82]對不同尺寸結塊間的碰撞引入力平衡模型。實驗結果表明，結塊的尺寸對計算作用力的未知參數非常敏感。通過兩個碰撞結塊的作用力平衡分析，獲得黏性細顆粒結塊的碰撞動力學，并提出一個團聚準則，這一準則定量解釋了高表面氣速、高流體密度、低顆粒黏性以及結塊間高頻率碰撞有利于黏性顆粒的結塊流態(tài)化。van Ommen 等[9]進一步完善力平衡模型，考慮內部孔隙結構結塊的滲透率對曳力的影響，以及氣固體系中其他化學物質（如異丙醇）對范德華力和納米顆粒結塊流化行為的影響。實驗結果表明，隨著流化氣速的增加，曳力和碰撞力增大，結塊破碎，粒徑變?。划惐嫉拇嬖谑笻amaker常數減小，從而范德華力減小，黏附作用力減小導致結塊尺寸減小。使用該模型計算的平衡結塊尺寸與實驗數據吻合較好。Tamadondar 等[83]認為分離力包括氣固間曳力、碰撞力和表觀重力，范德華力則是促使結塊間結合的內部黏附力的主要來源，并利用該理論預測結塊分離作用力與黏附作用力動態(tài)平衡時流化床內結塊的尺寸。該模型計算曳力時考慮結塊滲透率以及結塊周圍流場對結塊慣性效應的影響，對結塊進行了較全面的受力分析，預測的結塊尺寸和實驗數據之間的平均誤差在20%以內。但上述模型均缺乏結塊形成機理的分析，研究人員引入的校正因子存在大量假設，因此仍需不斷完善氣固流態(tài)化中力平衡模型以提高適用性。

細顆粒體系下的結塊對流動性能影響較大，但對原始顆粒的物理化學特性與控制結塊宏觀行為的定量關系研究較少，顆粒間鍵的性質與結塊的結構之間的關系還沒有確定。目前大多數力平衡模型只是分析結塊尺寸的平均值。事實上，流化床內結塊尺寸分布較寬。各種因素，包括顆粒表面特性、幾何形狀、毛細作用，是否影響或如何影響顆粒的團聚結塊行為，還有待進一步展開研究。

4 展 望

流化床中介尺度結構的研究是解決重油催化裂化、選礦冶礦、鍋爐燃燒等領域中多尺度放大和優(yōu)化工業(yè)過程設計的關鍵點和難點。隨著測量技術、流體力學、固體力學、數學、統(tǒng)計物理學以及計算機科學等相關學科和技術手段的發(fā)展，介尺度相關實驗和模擬方法都取得了重大進展，研究人員進一步了解了流化床中介尺度結構形成原因和動力學行為。然而，對單顆粒層次的離散動力學研究仍然處于初步階段。對于顆粒間相互作用力，研究人員主要考慮顆粒非接觸時的范德華力產生的黏附力，即采用簡單的解析式進行分析，未計入接觸形變的影響，忽略了顆粒接觸狀態(tài)下顆粒法向運動或滾動摩擦過程。其次，大部分介尺度結構研究都是對穩(wěn)態(tài)的介尺度結構進行受力分析，仍缺乏動態(tài)演化過程的受力分析。因此，在未來，研究人員有必要研究微觀受力與介尺度動態(tài)演化行為內在聯系，發(fā)展時間分辨的帶有顆粒力相互作用的介尺度結構動力學研究，為復雜的工業(yè)應用提供有效可控的介尺度結構調控技術。

符 號 說 明

A——Hamaker常數

AFM——原子力顯微鏡

Bo——黏結數

D——范德華力作用距離

dp——固相顆粒粒徑，m

L——流場特征尺度

PIV——粒子圖像測速法

PTV——粒子跟蹤測速法

R——光滑球形顆粒粒徑

RMS——顆粒粗糙度均方根高度

r——顆粒粗糙度半徑

St——Stokes數

U——流場特征速度

V——結塊相對碰撞速度

λ——表面粗糙度波長（峰間距離）

μ——流體黏度，Pa·s

ρp——固相顆粒密度,kg/m3

下角標

g——氣相

i——顆粒i

j——顆粒j

L——濾波后大尺度

p——固（顆粒）相

S——濾波后小尺度