|
張立棟，韋慶文，秦宏，王擎
（東北電力大學油頁巖綜合利用教育部工程研究中心，吉林吉林132012）
摘要：采用離散單元法DEM（discrete element method）對圓柱形生物質顆粒和鋼球顆粒在滾筒中的混合進行了數值模擬，分析了滾筒轉速和顆粒數量比對混合質量的影響。結果表明：在本文設定的工況下，顆粒的混合模式為階梯模式，并且顆粒在混合時可以分成3個區域，即左面的單層鋼球顆粒區、中間的鋼球顆粒和生物質顆?；旌蠀^、右面的生物質顆粒堆積區。左右兩邊的顆?；旌闲Ч^差，中間的顆?；旌闲Ч^好。當滾筒轉速相同時，鋼球顆粒和生物質顆粒數量比為3000∶200時的顆?；旌闲Ч蠕撉蝾w粒和生物質顆粒數量比為3000∶100時的好，即當鋼球顆粒數量遠大于生物質顆粒數量時，增加生物質顆粒的數量可以提高混合效果。在鋼球顆粒和生物質顆粒數量比相同的情況下，當滾筒轉速在5～25r/min的范圍內，滾筒轉速越高，顆粒的混合質量越好，并且顆粒混合達到穩定的時間就越短。
生物質主要包括植物廢棄物、禽畜糞便、城市垃圾等。生物質的利用多種多樣，如直接燃燒[1-3]、沼氣發酵、生物質熱解[4-5]等。顆粒物質是由大量相互作用的顆粒組成的復雜體系[6]。顆粒物質在自然界很常見，如積雪、泥石流、土壤、沙漠等，在日常生活中也很常見，如堆積的糧食、食鹽等。
目前，國內外不少研究學者對顆粒在滾筒中的混合進行了研究，主要研究了顆粒在滾筒中的運動狀態以及混合效果，分析了滾筒大小、滾筒轉速、物料填充率等因素對運動狀態以及混合效果的影響[7-12]。在DEM仿真模擬中，仿真參數對于DEM仿真十分敏感，ALIZADEH等[13]建立了以楊氏模量和摩擦系數為重點的DEM分析模型，推導出量綱歸一化運動方程和相應的量綱歸一化特征數，以此來研究仿真參數對顆粒運動的影響。
非球形顆粒在滾筒中的混合會表現出與球形顆?；旌喜煌奶匦裕跞鸱嫉?sup>[14]利用EDEM軟件對水平轉筒內大豆顆粒的運動進行了模擬。朱立平等[15]運用離散單元法建立了絲狀顆粒傳熱傳質數學模型。DUBE等[16]則運用粒子示蹤法研究非球形顆粒在滾筒中的運動。HOHNER等[17]通過比較實驗研究和數值模擬來分析顆粒形狀對顆?；旌系挠绊?。LU等[18]研究了非球形顆粒在水平滾筒內橫截面中的流動特性。陶賀等[19]運用球形顆粒拼接的方法建立了異徑玉米形顆粒模型、異徑橢球形顆粒模型以及異徑生物質顆粒模型，研究了它們在移動床中的運動情況。
顆粒在滾筒中混合的評價方法多種多樣，李少華等[20]對變異系數、接觸數以及Lacey指數這3種常用的混合度評價方法進行了分析，變異系數適合評價軸向混合，接觸數法主要適用于計算機數值模擬中涉及的混合，Lacey指數算法適合評價徑向混合。嚴建華等[21]用圖像法測量技術來獲得混合指標從而評價顆粒的混合程度。LIAO等[22]將混合過程用相機拍攝下來，并通過計算像素的方法計算黑白顆粒所占的比例，以此來得到顆粒的混合指數。呂春旺等[23]則研究了顆粒在滾筒冷渣機中的徑向擴散運動，通過引入擴散系數的概念來探討顆粒物料的擴散規律以及徑向擴散對傳熱的影響。
本文通過EDEM軟件模擬了圓柱形生物質顆粒和鋼球顆粒在滾筒中的混合，用接觸數指數作為混合指標比較了不同工況下顆粒的混合效果，分析了滾筒轉速和填充率的變化對顆粒混合效果的影響，并對滾筒劃分區域，研究了不同區域顆粒的混合效果。
1EDEM仿真
在EDEM仿真中生物質顆粒為圓柱形。本文通過球形拼接的方法建立了圓柱形顆粒模型，并改變生物質顆粒和鋼球顆粒數量比以及滾筒轉速來得到6種不同的工況，從而對其進行模擬。
1.1仿真參數
模擬中的滾筒采用圓形滾筒，其半徑為93mm，深度為57mm，滾筒材料為鋼材，其泊松比為0.25，剪切模量為7.5×1010Pa，密度為7800kg/m3，其與生物質顆粒的靜摩擦系數和動摩擦系數分別為0.3和0.01，與鋼球顆粒的靜摩擦系數和動摩擦系數分別為0.2和0.01。在仿真中，鋼球顆粒直徑為3mm，生物質顆粒的底面直徑為4mm、高度為6mm。生物質顆粒采用小麥秸稈的性質，仿真中用到的參數如表1和表2所示。
鋼球顆粒和生物質顆粒的數量比分別選取3000∶100和3000∶200，滾筒轉速分別取5r/min、15r/min和25r/min。一共形成6種工況，如表3所示。
1.2生物質顆粒模型的建立
由于EDEM軟件中顆粒模型的基本單元是球形顆粒，所以生物質顆粒要用球形顆粒拼接來合成。本文建立了45球元模型（球形顆粒直徑為2mm），中間1個小球外層8個小球，外層小球的表面與中間小球的球心相接觸，中間小球的表面也與外層小球的球心相接觸。一共5組這樣的顆粒（9×5=45），45球元模型如圖1所示。
1.3實驗對比
用離散單元法模擬顆粒的混合可以得到與實驗接近的顆粒分布規律，在前期工作中，對于顆粒在干餾爐以及帶抄板干餾爐內的混合，模擬與實驗結果基本吻合[10,24]。本實驗中的滾筒采用半徑為93mm、深度為47mm的中碳鋼滾筒，以直徑為3mm的鋼球和底面直徑為4mm、長度為6mm的圓柱形木條為填料，以與模擬相同的填充率和滾筒轉速，對6種工況做對比實驗。圖2對比了20s時實驗與模擬過程中各工況下兩種顆粒在混合時的顆粒分布。從圖2中可以看出，在各個工況下，實驗和模擬過程中的顆粒分布都可以分成3個區域：左側的單層鋼球顆粒區，在這個區域只存在鋼球顆粒，并且只有一層；中間的鋼球顆粒和生物質顆?；旌蠀^，在這個區域鋼球顆粒處于中心，生物質顆粒分布在顆粒自由層表面以及壁面處，即生物質顆粒分布在鋼球顆粒的外圍；右邊的生物質顆粒堆積區，這個區域只有生物質顆粒。模擬與實驗結果吻合良好，即本文所建立的模型是合理的。
2仿真結果分析
2.1顆粒在滾筒內的運動描述
顆粒在圓形滾筒內的運動共有6種模式，隨著滾筒轉速增加，分別經歷滑移、階梯、滾動、泄落、拋落、離心6種運動模式。通過對混合運動過程的觀察，可以看出在各工況下顆粒的運動模式都是階梯模式。整個顆粒群隨著滾筒的轉動交替地上升、下落。圖3是20s時的仿真顆?；旌线\動圖，可以明顯地看到顆粒在混合時可以分成3個區域，左面的單層鋼球顆粒區、中間的鋼球顆粒和生物質顆?；旌蠀^、右面的生物質顆粒堆積區。之所以會形成左面的單層鋼球顆粒區，是因為本文中的滾筒轉速和填充率較低。從圖3中可以看出，隨著滾筒轉速的提高，單層鋼球顆粒區域的長度在減小。所以低轉速和低填充率的情況下會形成單層鋼球顆粒區。在顆粒填充時鋼球顆粒在下，生物質顆粒在上，當滾筒轉動時，下方處于平流層中的鋼球顆粒隨著滾筒一起向上運動，而在上方的生物質顆粒處于活動層中，隨著滾筒的轉動向下運動，便堆積在右下方，而顆粒的混合時間不是足夠長，生物質顆粒不能充分混合在鋼球顆粒中，由此便形成了右邊的生物質顆粒堆積區。
2.2整體混合效果的比較
本文采用接觸數指數作為混合質量對顆粒的混合程度進行分析。接觸數指數用M來表示，見式(1)[20]。
圖4為滾筒轉速相同時接觸數指數隨時間的變化曲線。從圖4(a)中可以看出，滾筒轉速為5r/min時，鋼球顆粒和生物質顆粒的數量比為3000∶200時的M值一直比鋼球顆粒和生物質顆粒的數量比為3000∶100時的大，到20s時，顆粒數量比為2000∶200時的M值已經達到0.050，而20s時顆粒數量比為3000∶100時的M值只有0.027，因此在本文設定的工況下鋼球顆粒和生物質顆粒數量比為3000∶200時的顆?；旌闲Ч谩Ｒ驗樯镔|顆粒的數量比鋼球顆粒的數量小很多，在這種情況下，生物質顆粒所占的比例越大，接觸數指數越大，顆粒的混合效果越好。從圖4(a)中可以看出，鋼球顆粒和生物質顆粒的數量比為3000∶100時接觸數指數處于下降狀態，說明此時顆粒混合效果很差，并且到20s時仍然沒有上升的趨勢。鋼球顆粒和生物質顆粒的數量比為3000∶200時，接觸數指數的數值隨著時間的增加先下降后增加，以11s為界。11s之前，鋼球顆粒和生物質顆粒的混合質量在下降；11s之后，鋼球顆粒和生物質顆粒的混合質量在上升；到達20s時接觸數指數仍然有上升的趨勢，說明即使到了20s顆粒的混合仍然沒有穩定。
從圖4(b)和(c)中同樣可以看出，在滾筒轉速相同的情況下，鋼球顆粒和生物質顆粒的數量比為3000∶200時的接觸數指數數值比鋼球顆粒和生物質顆粒的數量比為3000∶100時的接觸數指數數值要大，因此混合效果也就更好。從圖4中可以看出，在滾筒轉速為15r/min時不管鋼球顆粒和生物質顆粒的數量比是3000∶200還是3000∶100，接觸數指數的數值都隨著時間的增加而上升，到20s時接觸數指數分別達到了0.074和0.043，但是仍都有上升的趨勢，說明在滾筒轉速為15r/min時，兩種顆粒數量比下的混合在20s時都沒有穩定。當滾筒轉速為25r/min時，當鋼球顆粒和生物質顆粒的數量比為3000∶100時，M值在10s之前呈現出上升趨勢，在10s之后趨于平穩，顆粒的接觸數指數穩定在0.045上下，顆粒的混合已經穩定下來。當鋼球顆粒和生物質顆粒的數量比為3000∶200時，接觸數指數在14s之前呈現出上升趨勢，在14s之后趨于平穩，顆粒的M值穩定在0.083上下，說明兩種顆粒在14s之后的混合已經穩定。所以在滾筒轉速相同的情況下，盡管鋼球顆粒和生物質顆粒的數量比為3000∶200時的接觸數指數數值比鋼球顆粒和生物質顆粒的數量比為3000∶100時的大，混合效果更好，但是鋼球顆粒和生物質顆粒數量比為3000∶200時兩種顆?；旌戏€定所需要的時間也比鋼球顆粒和生物質顆粒的數量比為3000:100時的要長。
鋼球顆粒和生物質顆粒的數量比相同而滾筒轉速不同的情況下接觸數指數隨時間的變化曲線如圖5所示。從圖5中可以看出，在本文設定的工況下，當鋼球顆粒和生物質顆粒數量比相同時，滾筒轉速越高，接觸數指數越高。同時滾筒轉速越高，顆?；旌线_到穩定的時間就越短。在鋼球顆粒和生物質顆粒的數量比為3000∶100的情況下，當滾筒轉速為25r/min時，顆粒在10s時已經混合穩定。而滾筒轉速為15r/min時，接觸數指數數值一直上升，顆粒的混合尚未穩定。滾筒轉速為5r/mim時，接觸數指數的數值一直呈現出下降趨勢，說明離顆?；旌戏€定還需要更長的時間。在鋼球顆粒和生物質顆粒的數量比為3000∶200的情況下，當滾筒轉速為25r/min時，顆粒的接觸數指數數值在14s后已經漸趨平緩，說明顆粒在14s時已經混合穩定。滾筒轉速為15r/min時，接觸數指數數值一直在上升。當滾筒轉速為5r/min時，接觸數指數的數值在11s之前呈現出下降趨勢，而在11s之后呈現出上升趨勢，但離顆?；旌戏€定仍然需要一段時間。
2.3各區域混合效果的比較
將滾筒進行4×1×4的網格劃分，將混合過程中有顆粒的部分劃分成4個區域，如圖6所示。并選取工況6對4個區域內顆粒的混合情況進行分析比較。
3結論
在本文設定的工況下，對滾筒轉速和顆粒數量比對生物質顆粒和鋼球顆粒在滾筒中混合特性的影響進行了研究，得出以下結論。
?。?）6種工況下，顆粒的混合模式為階梯模式，并且顆粒在混合時可以分成3個區域：左面的單層鋼球顆粒區，中間的鋼球顆粒和生物質顆?；旌蠀^，右面的生物質顆粒堆積區。
（2）當滾筒轉速相同時，鋼球顆粒和生物質顆粒數量比為3000∶200時的接觸數指數一直比鋼球顆粒和生物質顆粒數量比為3000∶100時的接觸數指數高。到20s時，當滾筒轉速為5r/min時，顆粒數量比為3000∶200和3000∶100時的接觸數分別達到0.050和0.027，在20s，滾筒轉速為15r/min時，顆粒數量比為3000∶200和3000∶100時的M值分別達到0.074和0.043，在滾筒轉速為25r/min時，在顆粒數量比為3000∶200和3000∶100這兩種情況下，顆?；旌戏€定后的M值分別穩定在0.083上下和0.045上下。即滾筒轉速相同時，顆粒數量比為3000∶200時的顆?；旌闲Ч阮w粒數量比為3000∶100時的顆粒混合效果好，因此當鋼球顆粒數量遠大于生物質顆粒數量時，增加生物質顆粒的數量可以提高混合效果。
?。?）在鋼球顆粒和生物質顆粒數量比相同的情況下，當滾筒轉速在5～25r/min的范圍內，滾筒轉速越高，接觸數指數越高，并且顆粒混合達到穩定的時間也各不相同，滾筒轉速為5r/min時M值的上升并不明顯，甚至還會出現下降的趨勢，轉速為15r/min時，兩種顆粒數量比下的M值一直呈上升趨勢，而滾筒轉速為25r/min，顆粒數量比為3000∶200和3000∶100時顆粒分別在14s和10s時混合穩定，即轉速越高，顆?；旌线_到穩定的時間就越短，顆粒的混合效果也越好。
?。?）生物質顆粒和鋼球顆粒在滾筒中混合時，區域1到區域4的均衡鋼球顆粒比例的穩定值分別為1、0.83、0.67和0.21，區域3的Pb穩定值最接近0.5，其混合效果最好，即左右兩邊的顆?；旌闲Ч^差，中間的顆粒混合效果較好。
參 考 文 獻
[1]壽恩廣，李詩媛，任強強，等.生物質與煤富氧混合燃燒特性研究[J].可再生能源，2014，32（10）：1551-1557.
[2]張永亮，趙立欣，姚宗路，等.生物質固體成型燃料燃燒顆粒物的數量和質量分布特性[J].農業工程學報，2013，29 （19）：185-191.
[3]王月喬，田宜水，侯書林，等.生物質顆粒燃燒器燃燒適應性試驗[J].農業工程學報，2014，30（7）：197-203.
[4]王建飛，趙建濤，李風海，等.煙煤與生物質快速共熱解產物特性分析[J].燃料化學學報，2015，43（6）：642-647.
[5]李帥，張居兵，樸桂林，等.木質生物質定溫熱解炭化特性研究[J].動力工程學報，2015，35（6）：495-500.
[6]孫其誠，王光謙.顆粒物質力學導論[M].北京：科學出版社，2009：1-10.
[7]金輝霞，栢娜，楊格蘭.滾筒內顆粒混合過程的DEM仿真研究[J].控制工程，2013，20（3）：566-570.
[8]張立棟，王麗偉，朱明亮，等.回轉干餾爐內油頁巖與固體熱載體顆粒軸向混合特性實驗[J].東北電力大學學報，2012，32（6）：63-66.
[9]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.
[10]張立棟，李少華，朱明亮，等.回轉干餾爐內抄板形式與雙組元顆?；旌线^程冷模數值研究 [J].中國電機工程學報，2012，32（11）：72-78.
[11]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.
[12]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.
[13]ALIZADEH E，BERTRAND F，CHAOUKI J.Comparison of DEM results and Lagrangian experimental data for the flow and mixing of granules in a rotating drum[J].AIChE Journal，2014，60（1）：60-75.
[14]王瑞芳，李占勇，竇如彪，等.水平轉筒內大豆顆粒隨機運動與混合特性模擬[J].農業機械學報，2013，44（6）：93-99.
[15]朱立平，秦霞，袁竹林，等.絲狀顆粒在滾筒橫向截面中的傳熱傳質特性[J].東南大學學報，2014，44（4）：756-763.
[16]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.
[17]HOHNER 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.
[18]LU G，JAMES R T，CHRISTOPH R M.Effect of wall rougheners on cross-sectional flow characteristic for non-spherical particle in a Horizontal rotating cylinder[J].Particuology，2014（12）：44-53.
[19]陶賀，鐘文琪，金保昇，等.異徑混合非鋼球顆粒在移動床中流動特性的數值模擬[J].東南大學學報，2012，42（6）：1122-1126.
[20]李少華，朱明亮，張立棟，等.回轉裝置內三組元顆粒徑向混合評價方法分析[J].化工進展，2013，32（6）：1224-1229.
[21]嚴建華，朱建新，馬增益，等.圖像法用于流化床顆?；旌咸匦缘难芯縖J].高?；瘜W工程學報，2006，20（5）：745-751.
[22]LIAO C C，HSIAU S S，KIWING T.Granular dynamics of a slurry in a rotating drum [J].Physical Review E，2010，82（1）：1-4.
[23]呂春旺，譚培來，劉柏謙，等.滾筒冷渣機中的顆粒徑向擴散運動研究[J].熱能動力工程，2014，29（5）：532-538.
[24]李少華，張立棟，張軒，等.回轉干餾爐內影響顆?；旌线\動因素的數值分析[J].中國電機工程學報，2011，31（2）：32-38. |
