柱狀生物質(zhì)顆粒在回轉(zhuǎn)干餾爐中的運動及導熱特性

張立棟，韋慶文，秦宏，王擎

（東北電力大學油頁巖綜合利用教育部工程研究中心，吉林 吉林 132012）

柱狀生物質(zhì)顆粒在回轉(zhuǎn)干餾爐中的運動及導熱特性

張立棟，韋慶文，秦宏，王擎

（東北電力大學油頁巖綜合利用教育部工程研究中心，吉林 吉林 132012）

為探究柱狀生物質(zhì)顆粒在固體熱載體（鋼球顆粒）回轉(zhuǎn)干餾爐（回轉(zhuǎn)爐）內(nèi)的運動及導熱特性，采用離散單元法模擬了柱狀生物質(zhì)顆粒和鋼球顆粒在回轉(zhuǎn)爐中的運動和導熱；分析了轉(zhuǎn)速和生物質(zhì)顆粒數(shù)量（鋼球顆粒數(shù)量固定）對顆粒運動和導熱的影響；探討了顆粒在回轉(zhuǎn)爐中的分布規(guī)律；用生物質(zhì)顆粒的平均溫度和溫度標準差來評價鋼球顆粒對生物質(zhì)顆粒的加熱速率和生物質(zhì)顆粒溫度分布的均勻程度。結(jié)果表明：在本文設(shè)定的工況下，顆粒在回轉(zhuǎn)爐中可以分成3個區(qū)域，即左側(cè)的單層鋼球顆粒區(qū)、中間的顆?；旌蠀^(qū)和右側(cè)的生物質(zhì)顆粒堆積區(qū)；單層鋼球顆粒區(qū)的不穩(wěn)定因素（不穩(wěn)定區(qū)、不穩(wěn)定線、不穩(wěn)定空穴）誘發(fā)了此區(qū)的周期性塌落；顆粒的運動模式為階梯模式；隨著轉(zhuǎn)速的升高，鋼球顆粒的加熱速率升高。在前20s，相較5r/min和25r/min，轉(zhuǎn)速為15r/min時的生物質(zhì)顆粒溫度分布均勻性最差。

柱狀生物質(zhì)顆粒；離散單元法；運動；導熱

生物質(zhì)主要包括植物廢棄物、禽畜糞便、城市垃圾等，生物質(zhì)的利用多種多樣，如直接燃燒、發(fā)酵、熱解等。

顆粒物質(zhì)是由大量相互作用的顆粒組成的復雜體系，顆粒物質(zhì)在化工過程中很常見，如冶金、食品工業(yè)、制藥工業(yè)、農(nóng)業(yè)等。目前，國內(nèi)外不少研究學者對顆粒在回轉(zhuǎn)裝置中的混合進行了研究，主要研究了顆粒在回轉(zhuǎn)裝置中的運動狀態(tài)以及混合效果，分析了回轉(zhuǎn)裝置的尺寸、轉(zhuǎn)速、物料填充率等因素對運動狀態(tài)以及混合效果的影響[1-6]。

非球形顆粒在回轉(zhuǎn)裝置中的混合會表現(xiàn)出與球形顆粒混合不同的特性，王瑞芳等[7]利用離散單元仿真軟件對水平轉(zhuǎn)筒內(nèi)大豆顆粒的運動進行了模擬。朱立平等[8]運用離散單元法（discrete element method，DEM）建立了絲狀顆粒傳熱傳質(zhì)數(shù)學模型。DUBE等[9]則運用粒子示蹤法研究非球形顆粒在滾筒中的運動。H?HNER等[10]通過實驗研究和數(shù)值模擬的比較來分析顆粒形狀對顆?；旌系挠绊?。

可以通過直接或間接方法對回轉(zhuǎn)裝置中的顆粒進行加熱，KOMOSSA等[11]通過實驗與模擬研究了單一粒徑顆粒在直接加熱轉(zhuǎn)鼓中的傳熱。GUI等[12]研究了顆粒在波形轉(zhuǎn)鼓內(nèi)的導熱分布形態(tài)特性。FIGUEROA等[13]通過改變滾筒的截面形狀、轉(zhuǎn)速和填充率來研究顆粒在滾筒中的混合速率和傳熱速率。

本文基于吉林省汪清縣龍騰能源開發(fā)有限公司新能源的綜合利用，以生物質(zhì)顆粒在固體熱載體（鋼球顆粒）回轉(zhuǎn)干餾爐（以下簡稱回轉(zhuǎn)爐）內(nèi)干餾為背景，以生物質(zhì)顆粒分解前在回轉(zhuǎn)爐內(nèi)與固體熱載體的混合運動和導熱為重點。采用離散單元法模擬了柱狀生物質(zhì)顆粒和固體熱載體鋼球顆粒在回轉(zhuǎn)爐中的運動和導熱。在工程實際中，因顆粒在回轉(zhuǎn)爐內(nèi)運動的復雜性，追蹤單個物料顆粒較為困難，難以獲得精確的顆粒運動及傳熱規(guī)律。而采用離散單元法（DEM）進行數(shù)值模擬可以方便地追蹤每個顆粒，進而可以對單個顆?；蛘哳w粒群的運動及傳熱進行研究。近年來離散單元法已經(jīng)成為研究顆粒體系的重要數(shù)值試驗工具[14]。因此本文借助離散單元法進行數(shù)值模擬，對柱狀生物質(zhì)顆粒在回轉(zhuǎn)爐中的運動規(guī)律和導熱特性進行初步探索。分析了顆粒在回轉(zhuǎn)爐中的分布規(guī)律；分別用生物質(zhì)顆粒的平均溫度和生物質(zhì)顆粒的溫度標準差來評價鋼球顆粒對生物質(zhì)顆粒的加熱速率和生物質(zhì)顆粒溫度分布的均勻程度。探討了回轉(zhuǎn)爐轉(zhuǎn)速和生物質(zhì)顆粒數(shù)量對顆粒運動及導熱的影響。為生物質(zhì)顆粒在回轉(zhuǎn)爐內(nèi)的混合運動特點提供理論依據(jù)，為確定合理的回轉(zhuǎn)爐轉(zhuǎn)速及生物質(zhì)顆粒填充量來干餾生物質(zhì)提供參考。

1 DEM模型

離散單元法由CUNDALL和STRACK提出[15]，它建立在基本的牛頓運動定律上，用于研究非連續(xù)性顆粒物質(zhì)的結(jié)構(gòu)和運動規(guī)律。顆粒在回轉(zhuǎn)爐中的運動屬于密集顆粒流，因此采用基于軟球模型的離散單元法。軟球模型將顆粒間的法向力簡化為彈簧和阻尼器，切向力簡化為彈簧、阻尼器和滑動器。不考慮顆粒表面變形，依據(jù)顆粒間的法向重疊量和切向位移計算接觸力。單個顆粒的運動控制方程為式(1)。

式中，mi、a、Ii、β分別為顆粒i的質(zhì)量、加速度、轉(zhuǎn)動慣量和角加速度；g為重力加速度；為顆粒i受到的合力（法向分量加切向分量）；為顆粒i受到的合力矩。

1.1 接觸力和傳熱模型

不考慮顆粒間的粘連作用，顆粒間的接觸力模型采用Hertz-Mindlin（No Slip）模型。在Hertz-Mindlin(No Slip)模型中，法向力基于Hertz接觸理論，切向力基于Mindlin-Deresiewicz理論[15]。

對于在回轉(zhuǎn)爐中用固體熱載體加熱生物質(zhì)顆粒，傳熱機制有：顆粒與空氣之間的對流傳熱、顆粒與空氣之間的輻射傳熱、顆粒與顆粒之間的導熱、顆粒之間的熱對流、顆粒內(nèi)部的導熱、顆粒與筒壁之間的導熱、顆粒與筒壁之間的輻射傳熱、顆粒與顆粒之間的輻射傳熱。根據(jù)DING等[13,16]的研究，在轉(zhuǎn)鼓中，通過空氣的傳熱只占很小的一部分，因此不考慮空氣的傳熱。熱對流是指由于流體的宏觀運動而引起的流體各部分之間發(fā)生相對位移，冷熱流體相互摻混所導致的熱量傳遞過程。熱對流僅能發(fā)生在流體中。單個顆粒為固體，當顆粒物質(zhì)運動時，顆粒整體表現(xiàn)出流體的性質(zhì)，因此顆粒之間的熱對流是顆粒之間的熱傳導和混合運動的綜合效果。在加入復雜的輻射模型之前，可首先探究主要的傳熱方式，忽略輻射傳熱[17]。本文主要考察顆粒之間的熱傳導和顆粒之間的熱對流，不考慮輻射傳熱，并且本文不研究顆粒內(nèi)部的溫度分布，所以視顆粒內(nèi)部的溫度分布均勻，忽略顆粒內(nèi)部的導熱。熱傳導模型為Hertz-Mindlin with Heat Conduction模型。

1.2 柱狀生物質(zhì)顆粒模型

由于離散單元法軟件（EDEM）中顆粒模型的基本單元是球形顆粒，所以柱狀生物質(zhì)顆粒要用球形顆粒拼接來合成[18]。本文建立了45球元模型（球形顆粒直徑為2mm），中間1個小球外層8個小球，外層小球的表面與中間小球的球心相接觸，中間小球的表面也與外層小球的球心相接觸。一共5組這樣的顆粒（9×5=45），45球元模型如圖1所示。

圖1 柱狀生物質(zhì)顆粒模型

1.3 仿真參數(shù)

考慮到合理的計算機計算量以及實驗室規(guī)模下回轉(zhuǎn)爐模型為小尺度回轉(zhuǎn)爐[19]，因此本文對小尺度物理模型，少物料填充量下，柱狀生物質(zhì)顆粒在固體熱載體法回轉(zhuǎn)干餾爐中的運動及導熱進行了數(shù)值模擬。作為一種研究方法，一些學者[1,7,14]做了相似的工作，對小尺度物理模型下的顆粒運動進行了研究。本文也采用這種方法對柱狀生物質(zhì)顆粒在回轉(zhuǎn)爐中的運動規(guī)律以及導熱特性進行了初步探討。對于實際回轉(zhuǎn)爐尺寸下的物料運動規(guī)律和導熱特性有待進一步的深入探討和研究。模擬中的回轉(zhuǎn)爐采用圓形回轉(zhuǎn)爐，其半徑為93mm，深度為57mm，回轉(zhuǎn)爐材料為鋼材，其泊松比為0.25，剪切模量為7.5×1010Pa，密度為7800kg/m3，其與生物質(zhì)顆粒的靜摩擦系數(shù)和動摩擦系數(shù)分別為0.3和0.01，與鋼球顆粒（固體熱載體）的靜摩擦系數(shù)和動摩擦系數(shù)分別為0.2和0.01。在仿真中，鋼球顆粒直徑為3mm，生物質(zhì)顆粒的底面直徑為4mm、高度為6mm。生物質(zhì)顆粒采用小麥秸稈的性質(zhì)，仿真中用到的參數(shù)如表1和表2所示。

固體熱載體鋼球顆粒數(shù)量為3000，鋼球顆粒的初始溫度為1073.15K，生物質(zhì)顆粒的數(shù)量分別選取100和200，生物質(zhì)顆粒的初始溫度為373.15K，回轉(zhuǎn)爐轉(zhuǎn)速分別取5r/min、15r/min和25r/min。一共分形成6種工況，如表3所示。前1.25s為顆粒生成過程，先生成鋼球顆粒，再生成生物質(zhì)顆粒?；剞D(zhuǎn)爐從1.25s開始轉(zhuǎn)動。

表1 顆粒的物理性質(zhì)

表2 顆粒的力學性質(zhì)

表3 各工況的對比

2 實驗驗證

冷態(tài)實驗中采用半徑為93mm，深度為47mm的中碳鋼回轉(zhuǎn)爐，以直徑為3mm的鋼球和底面直徑為4mm、長度為6mm的圓柱形木條為填料，以與模擬相同的顆粒數(shù)量和回轉(zhuǎn)爐轉(zhuǎn)速，對6種工況做對比實驗。圖2對比了工況1在20s時實驗與模擬過程中兩種顆粒在混合時的顆粒分布。從圖中可以看出，實驗和模擬過程中的顆粒分布都可以分成3個區(qū)域：左側(cè)的單層鋼球顆粒區(qū)，在這個區(qū)域只存在鋼球顆粒，并且只有一層；中間的鋼球顆粒和生物質(zhì)顆?；旌蠀^(qū)，在這個區(qū)域鋼球顆粒處于中心，生物質(zhì)顆粒分布在顆粒自由層表面以及壁面處，即生物質(zhì)顆粒分布在鋼球顆粒的外圍；右邊的生物質(zhì)顆粒堆積區(qū)，這個區(qū)域只有生物質(zhì)顆粒。通過對其他5個工況的比較，可以得到相同的結(jié)論。

圖2 顆粒在回轉(zhuǎn)爐中的分布

圖3為隨機選取的工況（工況5）下回轉(zhuǎn)爐中單層鋼球顆粒區(qū)的右視圖，圖中(a)為模擬，(b)為實驗，圖中有各種形態(tài)的顆粒分布，在單層鋼球顆粒區(qū)中，有顆粒緊密排布的穩(wěn)定區(qū)（如圖中矩形所示），顆粒無規(guī)則排布的不穩(wěn)定區(qū)（如圖中正六邊形所示），還有不穩(wěn)定線（如圖中橢圓所示）和不穩(wěn)定空穴（如圖中圓形所示）。

實驗與模擬結(jié)果吻合良好，驗證了模型的合理性。

圖3 單層鋼球顆粒區(qū)的顆粒分布

3 結(jié)果與分析

3.1 顆粒的運動

隨著回轉(zhuǎn)爐轉(zhuǎn)速的提高，顆粒在回轉(zhuǎn)爐中主要經(jīng)歷6種運動模式，依次為：滑移、階梯、滾動、泄落、拋落和離心[20]。在本文中，6種工況下顆粒的運動皆為階梯模式，即當顆粒整體受到的筒壁摩擦力和支持力的合力大于重力時，顆粒整體隨筒壁上升，當顆粒整體的重力大于筒壁摩擦力和支持力的合力時，顆粒整體下滑。顆粒在回轉(zhuǎn)爐中運動受到混合和偏析的共同作用?；旌蠙C理主要有對流混合、擴散混合和剪切混合[21]，偏析主要由體積差和密度差引起，偏析作用使體積大的顆粒在上、體積小的顆粒在下，密度小的顆粒在上、密度大的顆粒在下[22]。在本文中，在混合和偏析的共同作用下，顆粒在回轉(zhuǎn)爐中逐漸形成3個區(qū)域，左側(cè)的單層鋼球顆粒區(qū)、中間的鋼球顆粒和生物質(zhì)顆粒混合區(qū)和右側(cè)的生物質(zhì)顆粒堆積區(qū)。

3.1.1 單層鋼球顆粒區(qū)的塌落

在回轉(zhuǎn)爐運動過程中，當單層鋼球顆粒區(qū)受到擾動時，不穩(wěn)定區(qū)、不穩(wěn)定線、不穩(wěn)定空穴處的顆粒（如圖3所示）容易發(fā)生塌落。擾動的來源主要有兩種：在階梯運動模式下，顆粒從上升變成下滑的瞬間，顆粒的運動方向發(fā)生變化而產(chǎn)生擾動；中間顆?；旌蠀^(qū)的顆粒運動對單層鋼球顆粒區(qū)底部產(chǎn)生擾動。擾動促使了不穩(wěn)定因素（不穩(wěn)定區(qū)、不穩(wěn)定線和不穩(wěn)定空穴）上方的顆粒發(fā)生塌落，若不穩(wěn)定因素集中在單層鋼球顆粒區(qū)的底部或者較多不穩(wěn)定因素均勻分布在整個單層鋼球顆粒區(qū)時，則會導致單層鋼球顆粒區(qū)整體的塌落。當顆粒塌落后又有新的單層鋼球顆粒區(qū)生成，單層鋼球顆粒區(qū)一直重復著塌落-生成-塌落-生成這樣的過程。

3.1.2 生物質(zhì)顆粒在x軸方向上的分布概率

圖4為回轉(zhuǎn)爐的位置坐標。

圖4 回轉(zhuǎn)爐的位置坐標

定義生物質(zhì)顆粒在x軸方向上的分布概率Fx為式(2)。

式中，表示生物質(zhì)顆粒在某一時刻在x方向某一區(qū)間內(nèi)出現(xiàn)的個數(shù)（以0.01s為一個單位時刻，以10mm為一個單位區(qū)間，即表示生物質(zhì)顆粒在1.25s這一時刻在x方向上0～10mm區(qū)間內(nèi)出現(xiàn)的數(shù)量，表示生物質(zhì)顆粒在1.26s這一時刻在x方向上10～20mm區(qū)間內(nèi)出現(xiàn)的數(shù)量，依次類推）；Nt為生物質(zhì)顆粒的總數(shù)量；n為取樣時刻總數(shù)；1.25表示回轉(zhuǎn)爐從1.25s開始轉(zhuǎn)動，所以從1.25s開始取樣，每隔0.01s取一次樣。F10表示生物質(zhì)顆粒在1.25s到20s這段時間內(nèi)，在0～10mm區(qū)間內(nèi)平均出現(xiàn)的概率，F(xiàn)20表示生物質(zhì)顆粒在1.25s到20s這段時間內(nèi)，在10～20mm區(qū)間內(nèi)平均出現(xiàn)的概率，依次類推。

圖5為Fx隨x的變化曲線。從圖5中可以看出在0～110mm區(qū)間內(nèi)，生物質(zhì)顆粒出現(xiàn)的概率緩慢上升，在110～130mm區(qū)間內(nèi)急劇升高，各工況下的概率峰值都出現(xiàn)在110～130mm區(qū)間內(nèi)，即生物質(zhì)顆粒堆積區(qū)主要出現(xiàn)在110～130mm區(qū)間內(nèi)。F10不等于0，是因為在初始時刻，生物質(zhì)顆粒均勻分布在鋼球顆粒上方，隨著回轉(zhuǎn)爐的轉(zhuǎn)動，少量生物質(zhì)顆粒被帶到0～10mm區(qū)間內(nèi)。之后隨著回轉(zhuǎn)爐的轉(zhuǎn)動，顆粒在回轉(zhuǎn)爐的帶動下，整體上從左上方向右下方傾斜，生物質(zhì)顆粒逐漸向右下方堆積。而顆粒之間的混合作用又使部分生物質(zhì)顆粒與左側(cè)的鋼球顆?；旌希诨剞D(zhuǎn)爐的帶動下部分生物質(zhì)顆粒向左移動，但生物質(zhì)顆粒因體積比鋼球顆粒大，相比鋼球顆粒，生物質(zhì)顆粒的重心離筒壁的距離比鋼球顆粒遠，所以在生物質(zhì)顆粒剛要到達單層鋼球顆粒區(qū)時便滾落下來。生物質(zhì)顆粒幾乎不能再到達0～10mm的區(qū)間，所以F10很小。生物質(zhì)顆粒的體積比鋼球顆粒大，密度比鋼球顆粒小，所以在回轉(zhuǎn)爐轉(zhuǎn)動時，偏析作用使生物質(zhì)顆粒趨向于分布在鋼球顆粒上方，而在上方的生物質(zhì)顆粒隨著顆粒隨著顆粒整體的運動以及自身的滾落，又回到右下方。不斷有生物質(zhì)顆粒從右下方左移混入鋼球顆粒以及鋼球顆粒上方的生物質(zhì)顆粒運動到右下方，于是便形成了中間的生物質(zhì)顆粒和鋼球顆?；旌蠀^(qū)。而回轉(zhuǎn)爐的帶動僅能使少部分生物質(zhì)顆粒左移混入鋼球顆粒，仍有大量的生物質(zhì)顆粒堆積在右下方，所以在110～130mm區(qū)間內(nèi)Fx急劇上升，F(xiàn)x的峰值也出現(xiàn)在此區(qū)間內(nèi)。當轉(zhuǎn)速相同時，生物質(zhì)顆粒數(shù)量為100時的峰值比數(shù)量為200時的高，即顆粒的分布更加集中。當生物質(zhì)顆粒數(shù)量為100時，隨著轉(zhuǎn)速的升高，生物質(zhì)顆粒在110～120mm區(qū)間出現(xiàn)的概率減小，在130～140mm區(qū)間出現(xiàn)的概率增大，在回轉(zhuǎn)爐轉(zhuǎn)速為15r/min和25r/min時的峰值分別為0.288和0.286，比轉(zhuǎn)速為5r/min時的0.246高。因此隨著轉(zhuǎn)速的升高，有更多的顆粒集中在120～140mm區(qū)間內(nèi)。當生物質(zhì)顆粒數(shù)量為200時，隨著回轉(zhuǎn)爐轉(zhuǎn)速的升高，概率峰值在減小，并且轉(zhuǎn)速為5r/min和15r/min時，生物質(zhì)顆粒在120～130mm區(qū)間內(nèi)出現(xiàn)的概率最高，而回轉(zhuǎn)爐轉(zhuǎn)速為25r/min時，生物質(zhì)顆粒在110～120mm區(qū)間內(nèi)出現(xiàn)的概率最高，即隨著回轉(zhuǎn)爐轉(zhuǎn)速的升高，峰值左移。6種工況下，在130～140mm區(qū)間內(nèi)，生物質(zhì)顆粒出現(xiàn)的概率急劇下降，140 mm之后已經(jīng)基本上沒有顆粒了，即在x方向上，顆粒的運動范圍限制在140mm之前。

3.2 鋼球顆粒對生物質(zhì)顆粒的加熱效果

圖5 不同回轉(zhuǎn)爐轉(zhuǎn)速及生物質(zhì)顆粒數(shù)量在x方向上的分布概率

本文用生物質(zhì)顆粒的平均溫度來評價鋼球顆粒對生物質(zhì)顆粒的加熱速率，用生物質(zhì)顆粒的溫度標準差來評價生物質(zhì)顆粒溫度分布的均勻程度。

3.2.1 生物質(zhì)顆粒的溫升速率

圖6為生物質(zhì)顆粒的平均溫度隨時間的變化曲線。生物質(zhì)顆粒數(shù)量為100，回轉(zhuǎn)爐轉(zhuǎn)速為5r/min、15r/min和25r/min時曲線的平均斜率分別為0.120、0.189和0.177；生物質(zhì)顆粒數(shù)量為200，回轉(zhuǎn)爐轉(zhuǎn)速為5r/min、15r/min和25r/min時曲線的平均斜率分別為0.112、0.158和0.157。曲線的平均斜率越高，生物質(zhì)顆粒的溫升速率越快，表明鋼球顆粒對生物質(zhì)顆粒的加熱速率越快。顯然，鋼球顆粒數(shù)量固定，生物質(zhì)顆粒越少其溫升速率越快。當生物質(zhì)顆粒數(shù)量相同時，回轉(zhuǎn)爐轉(zhuǎn)速為15r/min和25r/min時的平均溫度明顯比轉(zhuǎn)速為5r/min時上升的快，即生物質(zhì)顆粒的溫升速率快。但15r/min和25r/min這兩種轉(zhuǎn)速下，生物質(zhì)顆粒的溫升速率相當。因為隨著回轉(zhuǎn)爐轉(zhuǎn)速增加，顆粒的運動更加劇烈，生物質(zhì)顆粒與鋼球顆粒接觸的概率增加，但生物質(zhì)顆粒與鋼球顆粒的接觸時間縮短，所以增加回轉(zhuǎn)爐轉(zhuǎn)速并不能一直提高生物質(zhì)顆粒的溫升速率。

圖6 生物質(zhì)顆粒平均溫度隨時間的變化曲線

3.2.2 生物質(zhì)顆粒溫度分布的均勻程度

標準差的公式為式(3)。

式中，為生物質(zhì)顆粒的平均溫度；Ti為第i個生物質(zhì)顆粒的溫度；N生物質(zhì)顆粒總數(shù)。

圖7為生物質(zhì)顆粒的溫度標準差隨時間的變化曲線，從圖7中可以看出當回轉(zhuǎn)爐轉(zhuǎn)速相同時，當生物質(zhì)顆粒的數(shù)量為200時的溫度標準差比生物質(zhì)顆粒數(shù)量為100時的小，即生物質(zhì)顆粒數(shù)量為200時，顆粒的溫度分布更加均勻。這與前面生物質(zhì)顆粒在x方向上分布概率所得到的結(jié)論相符，當回轉(zhuǎn)爐轉(zhuǎn)速相同時，生物質(zhì)顆粒數(shù)量為200時概率峰值比生物質(zhì)顆粒數(shù)量為100時小，即生物質(zhì)顆粒的分布較均勻。當生物質(zhì)顆粒的數(shù)量為100時，在16.5s之前，回轉(zhuǎn)爐轉(zhuǎn)速為5r/min時的標準差最小，16.5s之后，轉(zhuǎn)速為25r/min時的標準差最小。當生物質(zhì)顆粒數(shù)量為200時，在16.8s之前，回轉(zhuǎn)爐轉(zhuǎn)速為5r/min時的標準差最小，16.8s之后，轉(zhuǎn)速為25r/min時的標準差最小。在前20s任意時刻，當回轉(zhuǎn)爐轉(zhuǎn)速為15r/min時，生物質(zhì)顆粒的溫度標準差均不能達到最小，即顆粒溫度分布的均勻性不能達到最好。所以雖然當回轉(zhuǎn)爐轉(zhuǎn)速為5r/min時生物質(zhì)顆粒溫升較轉(zhuǎn)速為15r/min時慢，顆粒整體溫度較低，但溫度分布較均勻。而轉(zhuǎn)速為15r/min和25r/min時顆粒溫升速率相當，但轉(zhuǎn)速為25r/min時顆粒運動更加劇烈，有更多的生物質(zhì)顆粒有機會接觸到鋼球顆粒，所以轉(zhuǎn)速為25r/min時顆粒的溫度分布更加均勻。

圖7 生物質(zhì)顆粒溫度標準差隨時間的變化曲線

4 結(jié)論

在本文設(shè)定的工況下，通過對柱狀生物質(zhì)顆粒和鋼球顆粒（固體熱載體，數(shù)量固定為3000）混合運動和導熱的分析，可得到如下結(jié)論。

（1）顆粒在回轉(zhuǎn)爐中的運動模式為階梯模式。在混合和偏析共同作用下，顆粒在回轉(zhuǎn)爐中的分布分成3個區(qū)域，即左側(cè)的單層鋼球顆粒區(qū)、中間的鋼球顆粒和生物質(zhì)顆粒混合區(qū)、右側(cè)的生物質(zhì)顆粒堆積區(qū)。單層鋼球顆粒區(qū)的不穩(wěn)定因素（不穩(wěn)定區(qū)、不穩(wěn)定線、不穩(wěn)定空穴）誘發(fā)了此區(qū)的塌落，單層鋼球顆粒區(qū)重復著塌落-生成-塌落-生成這樣的過程。

（2）生物質(zhì)顆粒在x方向上集中分布在110～130mm區(qū)間內(nèi)。轉(zhuǎn)速相同時，生物質(zhì)顆粒數(shù)量為100時的顆粒分布更加集中。當生物質(zhì)顆粒數(shù)量為200時，隨著轉(zhuǎn)速的提高，生物質(zhì)在x方向的分布概率峰值左移。

（3）回轉(zhuǎn)爐轉(zhuǎn)速相同，生物質(zhì)顆粒越少，溫升速率越快。生物質(zhì)顆粒數(shù)量相同時，轉(zhuǎn)速越高，溫升速率越快，但轉(zhuǎn)速為15r/min和25r/min時生物質(zhì)顆粒平均溫度的變化趨勢相差不大，即提高轉(zhuǎn)速并不能一直提高生物質(zhì)顆粒的溫升速率；在前20s，轉(zhuǎn)速為15r/min時的生物質(zhì)顆粒溫度分布均勻性一直最差。

[1] 金輝霞，栢娜，楊格蘭.滾筒內(nèi)顆粒混合過程的DEM仿真研究[J].控制工程，2013，20（3）：566-570.JIN H X，BAI N，YANG G L. The DEM simulation study of the particle mixing process in the rotating cylinder[J]. Control Engineering of China，2013，20（3）：566-570.

[2] 張立棟，王麗偉，朱明亮，等.回轉(zhuǎn)干餾爐內(nèi)油頁巖與固體熱載體顆粒軸向混合特性實驗[J].東北電力大學學報，2012，32（6）：63-66.ZHANG L D，WANG L W，ZHU M L，et al. Experiment on axis mixing degree of oil shale with solid heat carrier particles in rotary retorting [J]. Journal of Northeast Dianli University，2012，32（6）：63-66.

[3] SUNKARA K A，HERZ F，SPECHT E，et al. Influence of flight design on the particle distribution of a flight rotating drum[J].Chemical Engineering Science，2013，90：101-109.

[4] 張立棟，李少華，朱明亮，等.回轉(zhuǎn)干餾爐內(nèi)抄板形式與雙組元顆?；旌线^程冷模數(shù)值研究[J].中國電機工程學報，2012，32（11）：72-78.ZHANG L D，LI S H，ZHU M L，et al. Cold mode numerical analysis of flights forms and two component particles mixing in rotary retorting[J]. Proccedings of the CSEE，2012，32（11）：72-78.

[5] HILL K M，GIOIA G，AMARAVADI D，et al. Moon patterns，sun patterns，and wave breaking in rotating granular mixtures [J].Complexity，2005，10（4）：79-86.

[6] PUYVELDE D R V，YOUNG B R，WILSON M A，et al. Experimental determination of transverse mixing kinetics in a rolling drum by image analysis[J]. Powder Technology，1999，106（3）：183-191.

[7] 王瑞芳，李占勇，竇如彪，等.水平轉(zhuǎn)筒內(nèi)大豆顆粒隨機運動與混合特性模擬[J].農(nóng)業(yè)機械學報，2013，44（6）：93-99.WANG R F，LI Z Y，DOU R B，et al. Simulation on random motion and mixing characteristic for soybean in rotary drum[J]. Transactions of the Chinese Society for Agricultural Machinery，2013，44（6）：93-99.

[8] 朱立平，秦霞，袁竹林，等.絲狀顆粒在滾筒橫向截面中的傳熱傳質(zhì)特性[J].東南大學學報.2014，44（4）：756-763.ZHU L P，QIN X，YUAN Z L，et al. Heat and mass transfer characteristics of filamentous particles in transverse section of rotary dryer [J]. Journal of Southeast University （Nature Science Edition），2014，44（4）：756-763.

[9] DUBE O，ALIZADEH E，CHAOUKI J，et al. Dynamics of non-spherical particles in a rotating drum[J]. Chemical Engineering Science，2013，101：486-502.

[10] H?HNER D，WIRTZ S，SCHERER V. A study on the influence of particle shape and shape approximation on particle mechanics in a rotating drum using the discrete element method[J]. Powder Technology，2014，253：256-265.

[11] KOMOSSA H，WIRTZ S，SCHERER V，et al. Heat transfer in indirect heated rotary drums filled with monodisperse spheres:Comparison of experiments with DEM simulations[J]. Powder Technology，2015，286：722-731.

[12] GUI N，F(xiàn)AN J R. Numerical study of heat conduction of granular particles in rotating wavy drums[J]. International Journal of Heat and Mass Transfer，2015，84：740-751.

[13] FIGUEROA I，VARGAS W L，MCCARTHY J J. Mixing and heat conduction in rotating tumblers[J]. Chemical Engineering Science，2010，65：1045-1054.

[14] 陳輝，劉義倫，趙先瓊，等.一元散體顆粒物料在回轉(zhuǎn)窯截面上的運動與混合[J].中國有色金屬學報，2015，25（5）：2575-2581.CHEN H，LIU Y L，ZHAO X Q，et al. Motion and mixing of mono-disperse granular material in cross section of rotary kiln[J]. The Chinese Journal of Nonferrous Metal，2015，25（5）：2575-2581.

[15] CUNDALL P A，STRACK O L. A discrete numerical model for granular assemblies [J]. Geotechnique，1979，29（1）：47-65.

[16] DING Y L，F(xiàn)ORSTER R N，SEVILLE J P K，et al. Some aspects of heat transfer in rolling mode rotating drums operated at low to medium temperature[J]. Powder Technology，2001，121：168-181.

[17] CHAUDHURI B，MUZZIO F J，TOMASSONE M S. Modeling of heat transfer in granular flow in rotating vessels[J]. Chemical Engineering Science，2006，61：6348-6360.

[18] 陶賀，鐘文琪，金保昇.采用多元顆粒模型對圓柱形顆粒在移動床中流動的離散單元法直接數(shù)值模擬[J].中國電機工程學報，2012，32（17）：13-19.TAO H，ZHONG W Q，JIN B S. Discrete element modeling of cylindrical particle flowing in the moving bed by multi-element particle model [J]. Proceedings of the CSEE，2012，32（17）：13-19.

[19] 王擎，肖冠華，孔祥釗，等.固體熱載體干餾樺甸油頁巖試驗研究[J].東北電力大學學報，2013，33（5）：15-21.WANG Q，XIAO G H，KONG X Z，et al. Experimental investigation of solid heat carrier retorting of Huadian oil shale [J]. Journal of Northeast Dianli University，2013，33（5）：15-21.

[20] HENEIN H，BRIMACOMBE J K，WATKINSON A P.Experimental study of transverse bed motion in rotary kilns[J]. Metallurgical Transactions B，1983，14（2）：191-205.

[21] OTTINO J M，KHAKHAR D V. Mixing and segregation of granular materials [J]. Annual Review of Fluid Mechanics，2000，32：55-91.

[22] 陳輝，肖友剛，趙先瓊，等.回轉(zhuǎn)窯內(nèi)二元顆粒物料的徑向混合[J].工程科學學報，2016，38（2）：194-199.CHEN H，XIAO Y G，ZHAO X Q，et al. Transverse mixing of binary solid material in a rotating kiln[J]. Chinese Journal of Engineering，2016，38（2）：194-199.

Motion and heat conduction in rotary retorting filled with cylindrical biomass particles

ZHANG Lidong，WEI Qingwen，QIN Hong，WANG Qing
（Engineering Research Center of Ministry of Education for Comprehensive Utilization of Oil Shale，Northeast Electric Power University，Jilin 132012，Jilin，China）

In order to understand the characteristics of motion and heat conduction of biomass particles and solid thermal carriers（steel particles） in rotary retorting，the motion and heat conduction of biomass particles and steel particles in rotary retorting were simulated using discrete element methods.The effect of the rotational speed and quantity of biomass particles on particle motion and heat conduction was analyzed. The distribution of particles in the rotary retorting was investigated. The mean temperature of biomass particles and the temperature standard deviation of biomass particles were used to evaluate the heating rate of steel particles and the uniformity of biomass particles＇temperature. The result showed that particles in the rotary retorting was divided into three regions: the monolayer steel particles region in the left of the particles region，the mixing region of biomass particles and steel particles in the middle of the particles region，and the accumulating region of biomass particles in the right of the particles region. The unstable factors（unstable regions，unstable lines and unstable holes）led to the periodic collapse of the monolayer steel particles region. The motion mode was slumping in all six situations. With the increasing of rotational speed，the heating rate of steel particles was increased.When the rotational speed was 15r/min，the uniformity of biomass particles＇temperature was worst in the first 20s comparing to the rotational speed of 5r/min and 25r/min.

cylindrical biomass particles；discrete element method；motion；heat conduction

TQ 051

A

1000–6613（2017）11–3993–07

10.16085/j.issn.1000-6613.2017-0333

2017-03-01；修改稿日期2017-03-20。

教育部長江學者和創(chuàng)新團隊發(fā)展計劃(IRT13052）、吉林省自然科學基金（20150101033JC）、吉林市科技計劃（201464044）及吉教科合字（2015-237）項目。

及聯(lián)系人：張立棟（1980—），男，博士，副教授，主要研究方向為油頁巖綜合利用及回轉(zhuǎn)裝置混合與分離。E-mail：nedu1015@aliyun.com。