王文軒，高亞男，洪 坤

(淮陰工學(xué)院化工學(xué)院，江蘇 淮安 223003)

能源是人類生存和社會發(fā)展的基石，其變革歷程大致經(jīng)歷了薪柴、煤炭、石油、天然氣等階段。我國具有豐富的生物質(zhì)資源，生物質(zhì)種類眾多、分布廣泛、儲量巨大，包括能源植物、農(nóng)林廢棄物、城市和工業(yè)廢物、畜禽糞便等。當(dāng)前，對大量可再生的生物質(zhì)資源進行高效清潔利用是我國比較理想的新能源發(fā)展途徑之一。因此，生物質(zhì)能作為一種極具開發(fā)和利用潛能的可再生清潔能源，對促進能源發(fā)展戰(zhàn)略、環(huán)境保護和農(nóng)村發(fā)展都具有重要意義。

1 生物質(zhì)能轉(zhuǎn)化技術(shù)

近年來，生物質(zhì)轉(zhuǎn)化技術(shù)正成為國內(nèi)外研究與開發(fā)的熱點。不同種類的生物質(zhì)在物理化學(xué)特性方面存在差異，其轉(zhuǎn)化利用的途徑也不盡相同。歸納而言，生物質(zhì)能源轉(zhuǎn)化利用的主要方法可分為三類[1]，如圖1所示，即燃燒法、生物化學(xué)轉(zhuǎn)化法和熱化學(xué)轉(zhuǎn)化法。

圖1 生物質(zhì)轉(zhuǎn)化技術(shù)分類Fig.1 Classification of biomass conversion technology

因為本文主要聚焦生物質(zhì)熱解過程的流動和混合特性，故此以下章節(jié)重點對生物質(zhì)熱解的流動特性模擬等方面進行系統(tǒng)綜述。

2 生物質(zhì)熱解模擬進展

生物質(zhì)熱解技術(shù)是一種清潔高效的生物質(zhì)利用新方法，其熱解產(chǎn)物能源利用率高，被認為是一種極具發(fā)展?jié)摿Φ纳镔|(zhì)能轉(zhuǎn)化技術(shù)。而生物質(zhì)熱解的核心場所是熱解反應(yīng)器裝置，它決定了熱解效率、熱解產(chǎn)物組成及其含量分布等。常見的熱解裝置類型包括固定床、流化床、氣流床、旋轉(zhuǎn)錐反應(yīng)器、旋轉(zhuǎn)爐等。其中，流化床反應(yīng)器因具有高效的混合與傳熱性能，滿足快速熱解對溫度及升溫速率的要求而被廣泛使用，在生物質(zhì)熱解領(lǐng)域也受到更多的關(guān)注與重視。隨著快速熱解工藝的日益發(fā)展，亟需對熱解流化床進行放大優(yōu)化以提高其生物質(zhì)的處理量，從而實現(xiàn)生物質(zhì)快速熱解的工業(yè)化應(yīng)用。在新型流化床反應(yīng)器開發(fā)過程中，實驗研究和理論分析都是不可或缺的重要手段。隨著計算能力的提升及計算理論模型的發(fā)展，計算機模擬已逐步成為熱解流化床研究的重要輔助手段。朱玉琴[2]、熊勤鋼[3]等都對熱解流化床的計算流體力學(xué)(Computational Fluid Dynamics，CFD)模擬進行了全面的綜述，系統(tǒng)總結(jié)分析了CFD模擬中常用的理論模型，包括氣-固相間曳力模型、熱解動力學(xué)模型、顆粒內(nèi)傳熱模型和顆粒頸縮模型。在熱解流化床模擬中，氣-固相間曳力顯得尤其重要，不僅影響生物質(zhì)顆粒流動行為，而且還會影響熱解產(chǎn)物分布?，F(xiàn)有熱解模擬文獻一般采用基于均勻分布的曳力模型進行模擬預(yù)測。然而，流化床內(nèi)具有顯著的非均勻結(jié)構(gòu)，如氣泡、顆粒聚團等結(jié)構(gòu)，這些結(jié)構(gòu)又與流動、熱量傳遞、質(zhì)量傳遞及化學(xué)反應(yīng)等“三傳一反”過程存在復(fù)雜的相互作用，導(dǎo)致這類問題的模擬放大一直是學(xué)術(shù)界和工業(yè)界極具挑戰(zhàn)性的難題。

2.1 流化床流動結(jié)構(gòu)特征

作為生物質(zhì)快速熱解過程的重要反應(yīng)場所，流化床反應(yīng)器內(nèi)部通常呈現(xiàn)出時空多尺度結(jié)構(gòu)特征，并伴有復(fù)雜的動態(tài)演變。除了流化氣體外，床內(nèi)一般會有一種或多種生物質(zhì)顆粒，有時也會加入沙子等惰性顆粒進行強化傳熱，因此流化床式生物質(zhì)熱解過程屬于典型的氣-固多相流系統(tǒng)。具體而言，如圖2所示，單個生物質(zhì)或惰性顆粒代表系統(tǒng)的微尺度，反應(yīng)器代表系統(tǒng)的宏尺度，介于兩者之間的局部結(jié)構(gòu)(顆粒聚團、氣泡等)則是典型的介尺度結(jié)構(gòu)。在通常的操作條件(流化氣速和生物質(zhì)顆粒入口流量等)下，惰性顆粒一般易處于床底部，呈現(xiàn)鼓泡流化狀態(tài)，而生物質(zhì)顆粒則更易處于床上部，呈現(xiàn)快速流化狀態(tài)。這種動態(tài)多變的結(jié)構(gòu)又會對“三傳一反”等過程產(chǎn)生重要影響，如影響生物質(zhì)和惰性顆粒的分離與混合、生物質(zhì)顆粒的停留時間、熱解產(chǎn)率及產(chǎn)品分布等等。近年來相關(guān)研究[4-5]逐漸將此類難題歸結(jié)為時空多尺度問題，其突破的關(guān)鍵在于認知和調(diào)控這種多尺度結(jié)構(gòu)及其耦合影響。

圖2 熱解流化床的多尺度結(jié)構(gòu)[3]Fig.2 The multi-scale structures of pyrolysis fluidized bed

2.2 熱解流化床模擬方法

在早期有關(guān)流化床研究的文獻中，報道了一些經(jīng)典的無量綱模型和兩相模型。例如，Davidson提出的氣泡模型[6]用來分析單個氣泡上升的流動行為。受限于當(dāng)時的計算條件，不能大規(guī)模數(shù)值求解，因而這些模型建模時都引入了一些合理的簡化近似。隨著計算機軟硬件及數(shù)值計算方法的快速發(fā)展，近年來借助計算機模擬可以輔助解決一些多相流工程問題，它是傳統(tǒng)逐級試驗放大的有力補充。與傳統(tǒng)試驗放大相比，使用計算機模擬進行工程放大及優(yōu)化，不僅可以顯著地縮短開發(fā)周期，而且可以節(jié)約大量人力、物力和財力。

根據(jù)兩相的描述方式以及網(wǎng)格尺度大小來劃分，如圖3所示，多相流體力學(xué)模型主要分為三種[7]：直接數(shù)值模擬(Direct Numerical Simulation，DNS)、顆粒軌道模型(Discrete Particle Modelling，DPM)和歐拉多流體模型(Eulerian Multi-Fluid Model，MFM)。其中，DNS屬于最微觀的模擬方法，氣相控制方程采用單相流的Navier-Stokes方程求解，固相跟蹤所有顆粒的運動邊界，從而揭示微元空間內(nèi)所有的流場細節(jié)。由于DNS所需計算量極大，一般僅用它來研究局部平衡狀態(tài)下的微元流動(微元內(nèi)顆粒數(shù)量級約為103)。而在顆粒軌道模型中，氣相采用歐拉坐標(biāo)系求解，固相采用拉格朗日坐標(biāo)系求解，可以跟蹤所有顆粒的運動軌跡(顆粒數(shù)量級約為106)。同樣受限于計算量，顆粒軌道模型仍難以用于顆粒達數(shù)十億級(顆粒數(shù)量級為>109)的實際工程問題計算。在歐拉多流體模型中，顆粒相被視為“擬流體”的連續(xù)介質(zhì)，這樣處理后就可以在歐拉-歐拉坐標(biāo)系下求解氣相和固相控制方程，它能夠捕捉到宏觀和局部的流動行為，但無法獲得每個顆粒的運動細節(jié)。由于歐拉多流體模型的計算量大大減少，故而被廣泛用于大尺度工程問題的研究。

圖3 三種不同尺度的多相流體力學(xué)模型[7]Fig.3 Three multi-phase hydrodynamic models for different geometries

在顆粒軌道模型中，應(yīng)用最為廣泛的是離散元法(Discrete Element Modelling，DEM)，DEM賦予了顆粒很多特性，包括顆粒碰撞、顆粒旋轉(zhuǎn)、顆粒形狀和粘結(jié)力等。將DEM模型與CFD結(jié)合便形成了CFD-DEM方法，可用于流-固兩相的數(shù)值仿真?；跉W拉-拉格朗日坐標(biāo)的CFD-DEM耦合法同樣也受計算量的限制，它多用于小規(guī)模實驗系統(tǒng)的研究，在工業(yè)應(yīng)用領(lǐng)域仍受限于計算能力和模型真實度等因素。在工業(yè)流化床模擬中取得成果并被廣泛使用的是歐拉-歐拉坐標(biāo)系下的雙流體模型(Two-Fluid Model，TFM)，然而，顆?！皵M流體”處理后，雙流體模型需要封閉顆粒粘度和壓力項，以及氣-固相間的作用力。近年來，在傳統(tǒng)CFD-DEM基礎(chǔ)上，新發(fā)展了粗?；嘞嗔Ｗ訂卧▉頊p少計算量，它采用虛擬粒子來打包表征若干具有相同動力學(xué)特性的顆粒，例如多相流質(zhì)點網(wǎng)格模型(Multi-Phase Particle-In-Cell，MP-PIC)[8]以及Coarse-grained CFD-DEM方法[9]。與傳統(tǒng)顆粒軌道模型相比，粗?；椒梢燥@著減少系統(tǒng)的計算量，但其模型應(yīng)用于工程實際問題仍處于發(fā)展階段。

下面簡要介紹DEM、MP-PIC和TFM這三種方法的特點：

(1)離散元法(DEM)

離散元法追蹤每個固體顆粒的運動軌跡，并在空間尺度上將氣相建模為至少比粒子直徑大一個數(shù)量級的連續(xù)體，如圖4(a)所示。DEM可以提供關(guān)于顆粒流場非常詳細的細節(jié)信息，例如顆粒軌跡、固體與壁面相互作用，這對于表征流化床中復(fù)雜的流體動力學(xué)是非常有用的。然而，其工業(yè)應(yīng)用還仍受計算能力和模型真實程度等方面限制。

圖4 兩種典型多相流建模示意圖[7]Fig.4 Schematic diagram of two typical multi-phase modelling

(2)多相流質(zhì)點網(wǎng)格模型(MP-PIC)

與傳統(tǒng)離散顆粒模型不同的是，MP-PIC采用顆粒相壓力梯度力取代顆粒碰撞來表征顆粒間相互作用力，并采用計算粒子來表征具有相同動力學(xué)特性的顆粒群。粗?；幚砗螅琈P-PIC能大大地降低粒子模擬的計算量，提高了工程問題的計算速度，但其適用性和準確性還仍有待發(fā)展和驗證，已有用于模擬半工業(yè)規(guī)模粒子系統(tǒng)的驗證研究[10-11]。

(3)雙流體模型(TFM)

在雙流體模型中，固體顆粒被視為以“擬流體”形式存在的連續(xù)相，流體和顆粒兩相是能夠充滿整個流場的連續(xù)介質(zhì)，且在流場中兩相能夠相互滲透、相互作用，如圖4(b)所示。與DEM和MP-PIC相比，TFM可以節(jié)省大量的計算資源，非常適合用于工程問題的數(shù)學(xué)建模，也是流化床模擬放大的首選方法[12-13]。

前已特別提及，流化床熱解過程涉及到復(fù)雜的多尺度問題。此類反應(yīng)器放大時，TFM模擬的計算精度主要取決于TFM守恒方程中本構(gòu)關(guān)系式的適用性或準確度。若本構(gòu)關(guān)系式缺失結(jié)構(gòu)特征的合理描述，如采用傳統(tǒng)的平均化方法處理，就會湮滅了結(jié)構(gòu)及其所帶來的影響，那么將會導(dǎo)致預(yù)測結(jié)果的可靠性大大降低，甚至與實際情況相悖[14]。因此，構(gòu)建與本征物理結(jié)構(gòu)相一致的數(shù)學(xué)模型，是認識和調(diào)控這類多尺度耦合問題的重要手段。

2.3 多尺度建模方法

多尺度方法是一種解決不同尺度間關(guān)聯(lián)的有效方法，粗略可劃分為三類[15]：描述型、關(guān)聯(lián)型和變分型。其中，描述型方法僅僅用于描述不同尺度上的各種結(jié)構(gòu)特征，不關(guān)注結(jié)構(gòu)的形成機理；關(guān)聯(lián)型方法通過微尺度相互作用分析進行大尺度上的粗粒化處理，屬于單向關(guān)聯(lián)；而變分型方法通過控制機理(如穩(wěn)定性條件)反映不同尺度之間的雙向關(guān)聯(lián)，這其中就涉及到介尺度結(jié)構(gòu)問題。

熱解流化床中，常見的介尺度結(jié)構(gòu)主要包括氣泡和顆粒聚團等。針對此類問題建模時，Glasser等[16]曾指出這些非均勻結(jié)構(gòu)具有相同的形成與傳播機理，屬于同一類非線性方程的求解范疇。認識到結(jié)構(gòu)及其影響的重要性后，近年來，越來越多的研究者開始聚焦介尺度結(jié)構(gòu)建模。國際上得到廣泛關(guān)注和應(yīng)用的當(dāng)屬變分型的能量最小多尺度(EMMS)方法[15]，它通過穩(wěn)定性條件分析不同尺度間的聯(lián)系，揭示產(chǎn)生結(jié)構(gòu)的主要機理，得到含有結(jié)構(gòu)參數(shù)的本構(gòu)關(guān)系式。縱觀EMMS模型30多年的發(fā)展歷程，從循環(huán)流化床中的顆粒聚團，到經(jīng)典鼓泡床中的氣泡，介尺度結(jié)構(gòu)及其耦合效應(yīng)始終是貫穿其模型研究的核心內(nèi)容。當(dāng)前EMMS研究大多集中于以顆粒聚團為特征的氣-固兩相流，然而應(yīng)用于以氣泡為特征的氣固多相流研究報道相對較少。

對于生物質(zhì)快速熱解流化床，情況也更為復(fù)雜，生物質(zhì)顆粒(一般為Geldart A或B顆粒)和惰性顆粒(一般為Geldart B或D顆粒)因物性差別大，兩者常處于不同的流化狀態(tài)，即惰性顆粒常見鼓泡流化，而生物質(zhì)顆粒常見快速流化。Xue等[17]利用歐拉多流體模型模擬了實驗室規(guī)模的快速熱解流化床并考察了生物質(zhì)組成的影響。也有研究者模擬考察了操作溫度、流化氣速和生物質(zhì)粒徑等對動力學(xué)和熱解產(chǎn)率的影響[18]。Lee等[19]采用模擬考察了反應(yīng)器幾何結(jié)構(gòu)對熱解過程中氣泡大小和數(shù)量的影響。然而，上述所有研究仍采用平均化方法或一些經(jīng)驗關(guān)聯(lián)式封閉本構(gòu)關(guān)系，以致無法準確表達內(nèi)秉性的非均勻結(jié)構(gòu)及其傳遞現(xiàn)象，其放大規(guī)律的可靠性也大大降低。本文將重點圍繞熱解流化床內(nèi)的結(jié)構(gòu)進行數(shù)學(xué)建模，以期為熱解流化床的放大及優(yōu)化積累理論基礎(chǔ)。

3 結(jié) 語

生物質(zhì)快速熱解是一種很有發(fā)展?jié)摿Φ纳镔|(zhì)能利用新技術(shù)。然而，生物質(zhì)熱解過程異常復(fù)雜，受到很多的因素影響，也存在著各種需要解決的問題。計算模擬方法對熱解流化床內(nèi)的流動特性研究具有非常重要的指導(dǎo)意義。本文系統(tǒng)綜述總結(jié)了生物質(zhì)熱解過程的特征結(jié)構(gòu)、各種熱解流化床模擬方法的優(yōu)缺點，以及近年發(fā)展的多尺度建模方法在其中的應(yīng)用前景。