張立棟，劉若云，王 擎，李少華，秦 宏

?

不同偏心距滾筒中大顆粒的平均速度脈動特性

張立棟，劉若云，王 擎，李少華，秦 宏

(東北電力大學(xué) 能源與動力工程學(xué)院，吉林 132012)

采用離散元數(shù)值模擬法對多元顆粒在不同偏心距工況下進行數(shù)值模擬，研究偏心滾筒對全部大顆粒的運動初始階段、穩(wěn)定階段的平均速度波動的影響。結(jié)果表明：在不同偏心距的偏心滾筒中，同粒徑顆粒的平均速度波動不同；偏心滾筒中顆粒平均速度曲線在每個運動周期都會出現(xiàn)波峰和波谷現(xiàn)象，且偏心距越大，該現(xiàn)象越明顯；在顆粒運動的初始階段，大粒徑顆粒平均速度波動范圍大，偏心距為20 mm的滾筒中所有大顆粒平均速度波動范圍最大；在顆粒運動的穩(wěn)定階段，隨著時間的延長，粒徑對于顆粒平均速度的影響不明顯。

偏心距；大顆粒；顆粒運動；平均速度

顆粒物質(zhì)由有大量離散單體固體顆粒構(gòu)成的復(fù)雜系統(tǒng)。顆粒物質(zhì)被廣泛應(yīng)用于工農(nóng)業(yè)生產(chǎn)和日常生活中，如金屬冶煉、礦物加工、藥品生產(chǎn)、食品加工 等[1?2]。滾筒是運輸和處理顆粒物質(zhì)的最常見且最重要的設(shè)備，廣泛應(yīng)用在干燥冷卻、熱解干餾、顆?；旌系裙I(yè)生產(chǎn)過程中，并且顆粒物質(zhì)體系在回轉(zhuǎn)筒內(nèi)運動混合的狀態(tài)直接影響生產(chǎn)的效率和工藝產(chǎn)品的質(zhì)量[3]。因此，研究和掌握滾筒內(nèi)顆粒物質(zhì)體系的運動混合狀態(tài)規(guī)律具有重要的理論意義和工程應(yīng)用價值。

正因為滾筒內(nèi)顆粒物質(zhì)的運動規(guī)律和混合規(guī)律在工農(nóng)業(yè)上的廣泛應(yīng)用，國內(nèi)外的科研人員對此進行了大量的研究。陳輝等[4]針對物料的散體顆粒性質(zhì)，采用離散單元法建立物料的運動模式，探究一元散體顆粒物料在回轉(zhuǎn)窯截面的運動與混合，得出回轉(zhuǎn)窯轉(zhuǎn)速越高，物料動能的耗散性越高。GUI等[5?7]將滾筒壁摩擦力對二元顆?；旌嫌绊戇M行數(shù)值模擬，發(fā)現(xiàn)筒壁摩擦力對于顆粒的運動產(chǎn)生阻力，并且抑制顆粒的混合程度；滾筒速度決定顆粒的主動驅(qū)動的摩擦力進而加快顆粒的運動和混合；除此之外，還發(fā)現(xiàn)顆粒間的大小比例對于顆粒流整體運動及混合的影響是最小的。還有一些科研人員探究筒形對于顆粒物質(zhì)混合的影響，例如葛良等[8]設(shè)計波型邊界的旋轉(zhuǎn)滾筒，探究顆粒的混合和導(dǎo)熱分布形態(tài)特性，相對于圓形滾筒，波型滾筒在較低速度下就可以達到圓型滾筒較高速度下的混合效果，波型滾筒較圓型滾筒有更好的混合效率。本文作者前期曾采用實驗與理論分析方法，以Lacey混合指數(shù)作為混合評價指標，研究了轉(zhuǎn)速、填充率對橢圓型混合器內(nèi)二元顆粒的徑向混合，為橢圓型混合器的生產(chǎn)應(yīng)用提供了重要的理論依據(jù)[9]。除針對圓型滾筒和圓顆粒的運動混合探究，白玨明等[10]研究通過離散元法模擬直觀獲得滾筒顆粒運動狀態(tài)，分析顆粒運動的特點[9]。張立棟等[11]采用DEM方法探究變轉(zhuǎn)速下滾筒內(nèi)二元顆粒的混合特性，認為最佳混合工況下的顆粒運動都處于拋落的運動模式，且無論是勻速還是變轉(zhuǎn)速滾動，長軸筒混合效果始終優(yōu)于短軸筒混合效果。GENG等[12]在滾筒干燥器中對桿狀顆粒特性進行三維數(shù)值模擬，模擬結(jié)果得到顆粒的運動周期越小，越有利于顆粒床內(nèi)部的擴散和混合，并且主動層顆粒最大速率隨轉(zhuǎn)速、填充率、滾筒直徑的增大而增大。王會等[13]則進行了圓錐料倉顆粒周期脈動特征研究，以明晰顆粒脈動特征，整個卸料過程中顆粒速度始終處于波動變換中，而且料倉的半錐角越小，卸料速度越穩(wěn)定。

滾筒中顆粒運動混合方面已經(jīng)有了豐富的研究成果可借鑒，本文作者將重點探究多元顆粒在滾筒中運動狀態(tài)和速度波動特性。

1 模型的建立

本實驗中以油頁巖固體熱載體干餾為應(yīng)用背景，顆粒間的接觸基于軟球模型簡化，顆粒的接觸模型采用經(jīng)典的無滑移接觸模型對滾筒內(nèi)的顆粒運動進行數(shù)值模擬，將顆粒與滾筒壁面以及顆粒與顆粒之間的接觸過程簡化為類似于彈簧振子的阻尼振動。

1.1 幾何模型

本文作者采用DEM(離散單元法)對滾筒內(nèi)5種粒徑顆粒運動過程進行研究，滾筒內(nèi)徑為84 mm、厚度為20 mm，研究滾筒幾何中心偏離滾筒旋轉(zhuǎn)中心的距離(簡稱偏心距)對4種粒徑顆粒速度脈動的影響。其中滾筒內(nèi)徑為84 mm、厚度為20 mm；3種偏心距分別為0=0 mm，1=10 mm，2=20 mm；滾筒以15 r/min進行定速轉(zhuǎn)動；轉(zhuǎn)動前，轉(zhuǎn)軸在顆粒體系的垂直正上方，黑點代表滾筒轉(zhuǎn)軸在徑向切面的位置[14]，如圖1所示。

圖1 滾筒內(nèi)顆粒初始位置及偏心距示意圖

1.2 DEM模擬條件

模擬中滾筒的密度設(shè)置為7800 kg/m3，彈性模量是1.82×109N/m2，泊松比為0.3[15?16]，小粒徑顆粒在下大粒徑顆粒在上的原則向滾筒內(nèi)填充，顆粒的相關(guān)數(shù)據(jù)見表1。

1.3 速度脈動波動分析方法

關(guān)于滾筒內(nèi)顆粒流的研究中，通常選取隨機的顆粒進行顆粒速度的探究[17?18]，但為了研究整體大顆粒在滾筒中的運動狀態(tài)，本實驗中提取所有大顆粒(包括直徑為2、3、4和5 mm的顆粒)的速度數(shù)據(jù)，處理之后，進而研究大顆粒速度脈動特性。

若以u表示滾筒內(nèi)部第個顆粒的速度，為顆粒數(shù)，則該時刻滾筒內(nèi)部所有顆粒的平均速度為

表1 顆粒的物理參數(shù)及力學(xué)性質(zhì)

2 模擬結(jié)果與分析

滾落模式是工程實際中最常見的顆粒運動形態(tài)，在本文中的模擬工況均為滾落模式[19?20]。由于顆粒在對滾筒運動的過程中從起始到結(jié)束會有不同的運動特性，根據(jù)顆粒的運動特性，起始時顆粒的平均速度波動頻率和波動范圍都具有不穩(wěn)定性，而之后階段的速度波動頻率和波動范圍都逐漸穩(wěn)定。本實驗中將顆粒的整個運動過程分為兩個階段單獨研究，即初始階段(1~17 s，即平均速度的前四個周期)和穩(wěn)定階段(18 s~結(jié)束，即第五周期至結(jié)束)，這樣有利于研究整體顆粒在運動過程中運動特性及平均速度的運動周期。

2.1 工況的初始階段速度分析

2.1.1 同粒徑顆粒平均速度在不同偏心工況下的比較分析

顆粒在滾筒中運動，活動層的顆粒在重力的作用下隨著滾筒的轉(zhuǎn)動滾落至平流層的底部[21]，顆粒的平均速度就會在每個周期內(nèi)出現(xiàn)速度波峰與速度波谷。

圖2所示為同粒徑顆粒初始階段平均速度比較。由圖2可以看出，同一粒徑下，偏心距越大其平均速度的波動范圍越大，平均速度周期幾乎相同，達到波峰與波谷的頻率幾乎相同；偏心距越大，顆粒距旋轉(zhuǎn)中心的距離越大，則其整體顆粒的平均速度波動范圍就越大。

在圖2所示的相同粒徑顆粒初始階段平均速度中，同粒徑不同偏心距的顆粒平均速度波峰的增長趨勢幾乎相同。但是隨著偏心距的增大，相鄰波谷會出現(xiàn)明顯的速度差，且速度差隨著偏心距的增大而增加。這是因為，隨著滾筒轉(zhuǎn)動，大部分相同粒徑顆粒在相鄰兩次平流層運動過程中，顆粒體系整體位置與轉(zhuǎn)軸距離不同，相應(yīng)的線速度差異較大，以至于顆粒的平均速度出現(xiàn)了兩個波谷的現(xiàn)象。以圖2(c)為例，其第三個周期中兩個波谷中較大的速度為9.71 s時的0.06 m/s，較小速度為11.34 s時的0.045 m/s。速度偏高的波谷，本實驗中稱為大波谷，速度偏低的波谷，稱為小波谷。大波谷速度隨偏心距增大而增高，依舊如圖2(c)中的第三周期，偏心距為20 mm的大波谷速度明顯高于偏心距為10 mm的大波谷速度，偏心距為0 mm的大波谷速度最低。

2.1.2 零偏心距滾筒中大顆粒平均速度分析

在偏心距為=0 mm的工況下，滾筒的旋轉(zhuǎn)半徑就是滾筒的半徑。因為其旋轉(zhuǎn)半徑穩(wěn)定，顆粒的運動周期也穩(wěn)定，平均速度沒有大小波谷的現(xiàn)象存在，見圖3(a)。

在初始階段，2 mm顆粒的平均速度較其他粒徑的顆粒先達到波峰，其次是3 mm顆粒，4 mm顆粒，最后是5 mm顆粒。這是因為啟動前的顆粒位置排布順序從上至下依次是5、4、3和2 mm，如圖1所示。由于重力的原因，粒徑大的顆粒，即質(zhì)量更大的顆粒在重力的作用下會獲得更大的動能，到達滾筒的底部；而粒徑較小的顆粒則無法到達滾筒的最底部，然后粒徑較小的顆粒提前到達該周期的速度峰值，粒徑較大顆粒提前進入平流層隨滾筒做圓周運動。

在5 mm顆粒平均速度上升的過程中，其速度幾乎呈直線上升，速度的波動較少，直接達到波峰速度0.17 m/s。隨著時間的延長，5 mm顆粒的平均速度波動逐漸穩(wěn)定，在6 s后其波谷速度穩(wěn)定在0.55~0.65 m/s之間，在11 s后其波峰速度穩(wěn)定在0.95~1.05 m/s之間。其原因是在起始時刻，5 mm顆粒沒有和其他粒徑的顆?；旌?，其他粒徑粒子對其幾乎沒有影響，所以其運動自由速度上升快。隨著5 mm顆粒與其他粒徑顆粒的混合，運動受其他粒徑顆粒的影響，速度波動變小，并與整體的顆粒運動周期逐漸一致。

圖2 同粒徑顆粒初始階段的平均速度比較

圖3 不同粒徑顆粒初始階段平均速度比較

同樣在初始階段，4 mm顆粒的平均速度波動范圍較5 mm 的小，其波峰速度在4 s后就穩(wěn)定在0.95~1.05 m/s之間，其波谷速度在4 s后穩(wěn)定在0.55~0.65 m/s之間。

3 mm顆粒平均速度的波動范圍更小，從啟動開始，其波峰速度穩(wěn)定在0.85~0.95 m/s，自11 s后3 mm顆粒的波峰速度穩(wěn)定在0.90~0.95 m/s，其波谷速度在初始階段緩慢穩(wěn)步增加，直至11 s后，其值穩(wěn)定在0.60~0.65 m/s。

與以上3種粒徑顆粒的平均速度相比，2 mm顆粒的波峰速度在初始階段是緩慢增加的，直至到9 s后，其波峰速度穩(wěn)定在0.90~0.95 m/s之間，其波谷速度也是在初始階段緩慢增加，直至9 s后，其波谷速度穩(wěn)定在0.60~0.65 m/s之間，速度波動更加平穩(wěn)。

綜上所述，顆粒的粒徑越大，在初始階段其速度最大峰值越大，隨著時間的延長，所有顆粒的平均速度逐漸平穩(wěn)，15 s后所有顆粒的波峰速度穩(wěn)定在0.90~0.95 m/s；啟動后，隨時間的增加，所有顆粒的波谷速度緩慢增加，逐漸穩(wěn)定，15 s后穩(wěn)定在0.60~0.65 m/s之間。

2.1.3 10 mm偏心距滾筒中大顆粒平均速度分析

在10 mm偏心距工況下，隨著滾筒繞旋轉(zhuǎn)中心旋轉(zhuǎn)，滾筒的運動半徑是變化的(見圖3(b))。觀察圖3(b)所有顆粒的平均速度波動范圍是趨于穩(wěn)定的，5 mm顆粒的波峰速度逐漸下降，從啟動開始的0.12 m/s到10 s后穩(wěn)定在0.95~1.05 m/s的范圍，在5 s后出現(xiàn)明顯的大小波谷，大波谷速度在7 s后穩(wěn)定在0.65~0.70 m/s之間，其小波谷速度在初始階段便穩(wěn)定在0.55~ 0.60 m/s之間，7 s之后小波谷速度穩(wěn)定在0.50~ 0.55 m/s間。

4 mm顆粒的波峰速度0.15 m/s從初始階段開始小幅度下降，其數(shù)值在10 s后穩(wěn)定在1.00~1.05 m/s之間，其波谷速度在5 s后便出現(xiàn)了清晰的大小之分，自10 s后大波谷速度穩(wěn)定在0.07 m/s左右，小波谷速度穩(wěn)定在0.50~0.55 m/s。

3 mm顆粒的波峰速度穩(wěn)定在0.10~0.11 m/s之間，5 s后同樣出現(xiàn)了清晰的大小波谷，大波谷速度穩(wěn)定0.07~0.75 m/s之間，小波谷速度則一直穩(wěn)定在0.50~0.55 m/s之間。

2 mm顆粒則與在偏心距為0 mm情況下的速度趨勢相似，從啟動開始其波峰速度呈緩慢增長趨勢，在9 s后穩(wěn)定在0.90 m/s左右，波谷速度在剛啟動的階段比較不穩(wěn)定，在5 s之后出現(xiàn)明顯穩(wěn)定的大小波谷之分，大波谷速度在0.65~0.70 m/s之間，小波谷速度在0.50~0.55 m/s之間。

2.1.4 20 mm偏心距滾筒中大顆粒平均速度分析

在偏心距是20 mm的工況下(見圖3(c))，5 mm顆粒平均速度的波動最大，波峰速度在初始階段波動范圍也大，這是因為在物料最初的放置5 mm顆粒在最上層，滾筒開始運動5 mm顆粒便處于活動層，由于重力的作用在滑落過程中獲得動能。由圖3可以看出，5 mm顆?；涞木嚯x最大，相應(yīng)地也獲得更多的動能，其平均速度最大，其次是4 mm顆粒，3 mm顆粒，最后是2 mm顆粒。

在顆粒隨滾筒運動的過程中，由于粒徑大的顆粒質(zhì)量大，慣性也大，大粒徑的顆粒較小粒徑顆粒相比，會遠離渦心，貼近筒壁，獲得更大的速度差(波峰速度與波谷速度之差)。如表2初始階段前4個周期內(nèi)速度差圖像，第一個周期顆粒粒徑越大其速度差越大，即平均速度波動隨顆粒粒徑的增大而增大，反之亦然；第二周期之后，不同粒徑顆粒平均速度差之間的差值變小，這是因為隨著運動的進行，不同粒徑顆粒開始混合，其運動狀態(tài)就會受到影響，不同粒徑顆粒之間的速度差就會越來越小。

5 mm顆粒在4 s之后速度波谷出現(xiàn)明顯的大小波谷現(xiàn)象(其他粒徑的顆粒也是在4 s之后出現(xiàn)這一現(xiàn)象)，這是因為在4 s后由于顆粒整體隨滾筒的旋轉(zhuǎn)進入運動軌跡，5 mm顆粒的絕大部分進入了兩個轉(zhuǎn)軸距離不同的平流層，便出現(xiàn)了大小波谷的現(xiàn)象。大波谷速度9 s后穩(wěn)定在0.75~0.85 m/s之間，小波谷速度穩(wěn)定在0.40~0.50 m/s之間。而其波峰速度變化不大，逐漸平穩(wěn)于12 s穩(wěn)定在0.95~1.05 m/s之間。

4 mm顆粒平均速度較5 mm顆粒平均速度波動的范圍小且平穩(wěn)，其波峰速度沒有較大的波動，在9 s后便穩(wěn)定在1.00~1.05m/s之間，大波谷速度在9 s后穩(wěn)定在0.80~0.85 m/s之間，小波谷速度穩(wěn)定在0.45~0.50 m/s之間。3 mm顆粒與4 mm、5 mm顆粒相比，其平均速度波動更加平穩(wěn)，波動的范圍更小，在啟動開始，其波峰速度就穩(wěn)定在0.95~0.10 m/s之間，大波谷速度穩(wěn)定在0.40~0.50 m/s之間，小波谷速度則在5 s以后穩(wěn)定在0.80~0.85 m/s之間。

2 mm顆粒的平均速度波動與3 mm、4 mm、5 mm顆粒的整體趨勢不同，在啟動開始到穩(wěn)定階段，其平均速度呈緩慢上升趨勢，波峰速度和波谷速度都在緩慢上升；波峰速度在8 s后逐漸穩(wěn)定在0.9~1.0 m/s之間，大波谷速度則在8 s后穩(wěn)定在0.8 m/s左右，小波谷速度則一直緩慢增加。這是因為隨滾筒開始轉(zhuǎn)動，大部分2 mm顆粒處于渦心的位置，速度變化較慢；隨著滾筒的轉(zhuǎn)動，2 mm顆粒逐漸分散至滾筒的各個位置，其顆粒的整體速度便慢慢增大，逐漸穩(wěn)定。因為在偏心距10 mm工況下，不同粒徑顆粒平均速度的變化趨勢與在偏心距20 mm工況下的情形幾乎相同，只是速度變化沒有那么明顯，所以本文主要分析在偏心距20 mm工況下的速度變化趨勢。

2.2 穩(wěn)定階段平均速度分析

2.2.1 不同偏心距下的同粒徑顆粒平均速度分析

該部分分析穩(wěn)定階段平均速度的波動特性，由圖4(a)可知，偏心距為0 mm時，同粒徑的顆粒平均速度波動頻率是不一致的，速度波峰和速度波谷出現(xiàn)位置的偏差較大；隨時間的延長則速度波峰和速度波谷位置逐漸一致，到50 s左右，速度波峰和速度波谷位置幾乎相同。當偏心距為10和20 mm時，同樣出現(xiàn)該現(xiàn)象，只是偏心距為10 mm時在40 s左右其速度波峰和速度波谷的位置達到一致，而偏心距為20 mm時則在30 s左右其速度波峰和速度波谷位置達到一致。3種不同工況下的同粒徑顆粒在到達波動頻率一致后，其平均速度的波動范圍幾乎相同。這是因為顆粒在隨滾筒運動的過程中，大顆粒之間會有較大的接觸空隙，小顆粒由于尺寸較小則滲透到大顆粒與大顆粒之間，并且隨著滾筒的運動，大顆粒將逐漸被小顆粒包圍，這時大小顆粒的運動會被彼此束縛，致使所有顆粒的速度波動頻率逐漸一致，最后達到幾乎相同。

由圖4(a)、(b)、(c)、(d)對比可以明顯看出，不同工況下的同粒徑顆粒偏心距越大，其速度的波動范圍就越大，即?v0 mm≤?v10 mm≤?v20 mm。

2.2.2 同偏心距下的不同粒徑顆粒平均速度分析

由圖5(a)可以看出，在同偏心距下，不同粒徑顆粒的平均速度周期是比較混亂，不同粒徑顆粒都有自己的運動周期，彼此之間的影響較小。平均速度的波動在0.60~0.90 m/s范圍內(nèi)時，這種現(xiàn)象一直持續(xù)。

圖4 相同粒徑顆粒穩(wěn)定階段平均速度比較

圖5 不同粒徑穩(wěn)定階段平均速度比較

如圖5(b)所示，在偏心距為10 mm的工況下，不同粒徑顆粒的平均速度運動周期開始逐漸一致，隨著時間的延長，顆粒的平均速度周期明顯一致，并且大致穩(wěn)定在0.60~0.90 m/s之間。隨著偏心距的增加，這種現(xiàn)象越明顯。

在偏心距為20 mm的工況下，顆粒的平均速度周期明顯比其他兩種工況下的要統(tǒng)一得多，并且其平均速度的波動范圍為0.05~0.10 m/s，明顯大于前兩種工況時的。這是因為隨偏心距的增大，顆粒受滾筒的影響也增加，所有顆粒在偏心滾筒中隨滾筒的轉(zhuǎn)動而運動，所以整體顆粒的平均速度波動周期逐漸一致。

圖6所示為所有顆粒在整個運動階段平均速度的標準差。由圖6可以發(fā)現(xiàn)，隨著顆粒粒徑的增大，其平均速度標準差也是呈增加的趨勢。比較相同粒徑顆粒在不同偏心距滾筒中的平均速度標準差可知，偏心距為20 mm的大于偏心距為0 mm的大于偏心距為10 mm的。

圖6 不同偏心距下不同粒徑平均速度的標準差

3 結(jié)論

1) 在不同偏心距工況下，相同粒徑顆粒的平均速度波動不同，其平均速度曲線在每個運動周期內(nèi)會出現(xiàn)波峰和波谷現(xiàn)象，且偏心距越大，平均速度的波動范圍就越大。

2) 在顆粒運動的初始階段，粒徑越大，顆粒平均速度曲線波動范圍越大；偏心距越大，顆粒平均速度波動范圍越大。

3) 在顆粒運動的穩(wěn)定階段，隨時間的延長，粒徑對于顆粒平均速度影響不大；偏心距為20 mm時顆粒平均速度的波動范圍最大。

[1] ISABEL F, VARGAS W L, MCCARTAHY J J. Mixing and heat conduction in rotating tumblers[J]. Chemical Engineering Science, 2010, 65(2): 1045?1054.

[2] TADA E F R, GRAJALES L M, LEMOS Y P, THOMEO J C. Mixture and motion of sugar cane bagasse in a rotating drum[J]. Powder Technology, 2017, 317: 301?309.

[3] 李少華, 張立棟, 張 軒, 余侃勝, 王 擎, 郝萬軍. 回轉(zhuǎn)式干餾爐內(nèi)影響顆?；旌线\動因素的數(shù)值分析[J]. 中國電機工程學(xué)報, 2011, 31(2): 32?38. LI Shao-hua, ZHANG Li-dong, ZHANG Xuan, YU Kan-sheng, WANG Qing, HAO Wan-jun. Numerical analysis of particle mixing and movement in rotary retorting[J]. Proceedings of the CSEE, 2011, 31(2): 32?38.

[4] 陳 輝, 劉義倫, 趙先瓊, 肖友剛, 劉 穎. 一元散體物料在回轉(zhuǎn)窯截面上的運動和混合[J]. 中國有色金屬學(xué)報, 2015, 25(5): 2575?2581. CHEN Hui, LIU Yi-lun, ZHAO Xian-qiong, XIAO You-gang, LIU Ying. Motion and mixing of mono-disperse granular material in cross section of rotary kiln[J]. The Chinese Journal of Nonferrous Metals, 2015, 25(5): 2575?2581.

[5] ZHANG Z W, GUI N, GE L, LI Z L. Numerical study of mixing of binary-sized particles in rotating tumblers on the effects of end-walls and size ratios[J]. Powder Technology, 2016, 314: 164?174.

[6] AISSA A A, DUCHESNE C, RODRIGUE D. Effect of friction coefficient and density on mixing particles in the rolling regime[J]. Powder Technology, 2010, 212: 340?347.

[7] GENG F, LI Y M, YUAN L J, LIU M, WANG X Y, YUAN Z L,YAN Y-M, LUO D S. Experimental study on the space time of flexible filamentous particles in a rotary dryer[J]. Experimental Thermal and Fluid Science, 2011, 44: 708?715.

[8] 葛 良, 桂 南, 徐文凱, 閆 潔. 波形滾筒內(nèi)顆?；旌虾蛯?dǎo)熱分布形態(tài)特性的研究[J]. 中國科學(xué), 2014, 44(1): 62?70. GE Liang, GUI Nan, XU Wen-kai, YAN Jie. Numerical study of the distribution characteristics on particle mixing and thermal conduction in wavy drum[J]. Science China, 2014, 44(1): 62?70.

[9] 張立棟, 李連好, 王 擎, 秦 宏, 李少華. 橢圓形混合型內(nèi)二元顆粒徑向混合[J]. 中國有色金屬學(xué)報, 2017, 27(4): 825?832. ZHANG Li-dong, LI Lian-hao, WANG Qing, QIN Hong, LI Shao-hua. Transverse mixing of binary granular in elliptical mixer[J]. The Chinese Journal of Nonferrous Metals, 2017, 27(4): 825?832.

[10] 白玨明, 劉柏謙, 白玨瑩, 譚培來. 滾筒冷凝器內(nèi)顆粒床運動狀態(tài)的離散單元法模擬[J]. 熱力發(fā)電, 2016, 45(2): 75?85. BAI Jue-yu, LIU Bo-qian, BAI Jue-ying, TAN Pei-lai. DEM simulations of solid motion in rotary drum ash cooler[J]. Thermal Power Generation, 2016, 45(2): 75?85.

[11] 張立棟, 程 碩, 李少華, 王 擎. 變轉(zhuǎn)速回轉(zhuǎn)筒內(nèi)二元顆粒混合特性數(shù)值分析[J]. 東北電力大學(xué)學(xué)報, 2018, 38(3): 39?45. ZHANG Li-dong, CHENG Suo, LI Shao-hua, WANG Qing. Numerical simulation study on the mixing characteristics of the binary particle in a variable speed rotating device[J]. Journal of Northeast Dianli University, 2018, 38(3): 39?45.

[12] GENG F, GANG L Y, WANG Y C, LI Y M, YUAN Z L. Numerical investigation on particle mixing in a ball mill[J]. Powder Technology, 2015, 292: 64?73.

[13] 王 會, 賈富國, 韓燕龍, 張亞雄, 曹 斌. 圓錐料倉內(nèi)顆粒周期性脈動特征研究[J]. 物理學(xué)報, 2016, 66(1): 192?202. WANG Hui, JIA Fu-guo, HAN Yan-long, ZHANG Ya-xiong, CAO Bin. Cyclical pulsation properties of particles in cone sile[J]. Acta Physica Sinica, 2016, 66(1): 192?202.

[14] 張立棟, 李連好, 秦 宏, 王 擎. 多粒徑顆粒在圓形偏心滾筒內(nèi)的運動混合[J]. 化工進展, 2017, 36(2): 451?456. ZHANG Li-dong, LI Lian-hao, QIN Hong, WANG Qing. Research on the motion and mixing of polydisperse particles in a circular eccentric drum[J]. Chemical Industry and Engineering Progress, 2017, 36(2): 451?456.

[15] BASINSKAS G, SAKA M. Numerical study of the mixing efficiency of a batch mixer using the discrete element method[J]. Powder Technology, 2016, 301: 815?829.

[16] SANTOS D A, BARROZO M A S, DUARTE C R, WEIGLER F, MELLMANN J. Investigation of particle dynamics in a rotary drum by means of experiments and numerical simulations using DEM[J]. Advanced Powder Technology, 2016, 27: 692?703.

[17] HASSANPOUR A, TAN H, BAYLY A, GOPALKRIDHNAN P, NG B, ALI M G. Analysis of particle motion in a paddle mixer using discrete element method (DEM)[J]. Powder Technology, 2010, 206: 189?194.

[18] GUI N, YANG X T, TU J Y, JIANG S Y. Effect of roundness on the discharge flow of granular particles[J]. Powder Technology, 2016, 314: 140?147.

[19] HALIDAN M, CHANDRATILLEKE G R, CHAN S L I, YU A B, Bridgwater J. Prediction of the mixing behaviour of binary mixtures of particles in a bladed mixer[J]. Chemical Engineering Science, 2014, 120: 38?48.

[20] PASHA M, HASSANPOUR A, AHMADIAN H, TANC H S, BAYLY A, GHADIRI M. A comparative analysis of particle tracking in a mixer by discrete element method and positron emission particle tracking[J]. Powder Technology, 2015, 270: 569?574.

[21] YANG R Y, ZOU R P, YU A B. Microdynamic analysis of particle flow in a horizontal rotating drum[J]. Powder Technology, 2003, 130: 138?146.

Pulsation characteristics of larger particle average velocity in different eccentric rollers

ZHANG Li-dong, LIU Ruo-yun, WANG Qing, LI Shao-hua, QIN Hong

(School of Energy and Power Engineering, Northeast Electric Power University, Jilin 132012, China)

A numerical study on multiple particles under different eccentric conditions by DEM(Discrete element method) was presented. The aim is to research the influence of average speed fluctuation of larger particles by the eccentric rollers in initial and stable stages. The results show that the particles with the same size have different fluctuations of average velocity in different eccentric rollers; the average speed curves of particles in eccentric cylinder appear the peaks and troughs in a motion period, and greater eccentric rollers arise more obvious phenomenon. At the initial stage, bigger size particles have a wide range of average velocity fluctuation, the average speeds of all particles have the widest range in the 20 mm eccentricity drum. At the stable stage, the sizes of particles have no obvious impact on the average speed of particles when the particle movement gradually become stable.

eccentricity; largeparticle; particle motion; average velocity

Project(IRT17R19) supported by Program of Changjiang Scholars and Innovative Research Team in University, China; Project(20150101033JC) supported by the Natural Science Foundation of Jilin Province, China; Project(201464044) supported by the Science and Technology Development Projects of Jilin City, China

2017-09-26;

2018-05-08

ZHANG Li-dong; Tel: +86-13843225181; E-mail: nedu1015@aliyun.com

教育部長江學(xué)者和創(chuàng)新團隊發(fā)展計劃資助項目(IRT_17R19)；吉林省自然科學(xué)基金資助項目(20150101033JC)；吉林市科技計劃項目(201464044)

2017-09-26；

2018-05-08

張立棟，副教授，博士；電話：13843225181；E-mail: nedu1015@aliyun.com

10.19476/j.ysxb.1004.0609.2018.10.25

1004-0609(2018)-10-2164-08

A

(編輯 龍懷中)