孫棟，陸海峰，曹嘉琨，吳雨婷，郭曉鐳，龔欣

（華東理工大學(xué)上海煤氣化工程技術(shù)研究中心，上海200237）

引 言

作為儲(chǔ)運(yùn)粉體的重要單元操作設(shè)備，料倉(cāng)簡(jiǎn)單的結(jié)構(gòu)獲得了工業(yè)界的青睞[1-2]，在其中復(fù)雜的顆粒流動(dòng)現(xiàn)象也受到科學(xué)研究者的廣泛關(guān)注[3-5]。近50年間，圍繞顆粒靜態(tài)堆積、動(dòng)態(tài)流動(dòng)和塑形流變等行為的研究大量開(kāi)展[6]，但由于受到顆粒間非線性和不可逆的受力-變形行為[7]，力鏈結(jié)構(gòu)強(qiáng)度迥異、互相交接[8]等復(fù)雜因素的影響，顆粒的流動(dòng)機(jī)理仍沒(méi)有得到充分掌握[9]。受限于復(fù)雜的顆粒特性、迥異的操作條件和變化的結(jié)構(gòu)參數(shù)，顆粒在料倉(cāng)中展現(xiàn)出不同的流動(dòng)行為，例如質(zhì)量流、中心流，在某些極端的情況下，會(huì)出現(xiàn)結(jié)拱、鼠洞、死區(qū)等異常現(xiàn)象[10-11]。

顆粒流動(dòng)順暢、連續(xù)以及料倉(cāng)結(jié)構(gòu)穩(wěn)定、耐久[12]是評(píng)價(jià)下料過(guò)程的重要指標(biāo)。通氣[2,13]、振動(dòng)[14-15]和改流體[16-17]是調(diào)節(jié)料倉(cāng)中顆粒流動(dòng)形態(tài)，改善不正常操作的常用手段。Lu 等[18]在中試裝置中研究了通氣對(duì)粉煤下料的影響，依據(jù)下料過(guò)程的團(tuán)聚模型成功預(yù)測(cè)了下料流率。Matchett[19]建立了拱穩(wěn)定模型，考慮了主動(dòng)態(tài)和被動(dòng)態(tài)不同的作用規(guī)律，指導(dǎo)了振動(dòng)料倉(cāng)的下料過(guò)程。作為一種簡(jiǎn)便的料倉(cāng)流動(dòng)調(diào)節(jié)手段，改流體對(duì)流型的影響得到了廣泛關(guān)注[20-21]，但其對(duì)流率的促進(jìn)作用往往被忽視，缺乏改流體料倉(cāng)下料流率的預(yù)測(cè)方法。

在平底圓柱料倉(cāng)中，Beverloo 方程[22]被廣泛用來(lái)預(yù)測(cè)下料流率。為了提高預(yù)測(cè)精度，擴(kuò)展適用范圍，Beverloo 方程被反復(fù)修正。Tierrie 等[23]采用攝像捕捉的方法，對(duì)顆粒的形狀進(jìn)行了描述，給出了修正Beverloo 方程中形狀參數(shù)的指導(dǎo)方法。Lu 等[4]從顆粒間作用力的角度出發(fā)，考察了表面粗糙度對(duì)下料流率的影響，基于Bond 數(shù)劃分了下料的流型，預(yù)測(cè)了下料流率。在錐形料倉(cāng)中，顆粒沿壁面滑動(dòng)甚至滾落，必將對(duì)下料流率產(chǎn)生影響，Khashayar 等[24]認(rèn)為錐形料倉(cāng)不能視為平底料倉(cāng)的特例。在半錐角α 的料倉(cāng)中，Rose 等[25]引入校正因子f = (tanα tanφd)-0.35來(lái)反映料倉(cāng)半錐角的作用效果?；诹ζ胶獾脑恚珻arleton[26]得到了錐形料倉(cāng)出口處顆粒的平均速度，對(duì)下料流率進(jìn)行了預(yù)測(cè)。作為調(diào)節(jié)流動(dòng)的常用手段，掌握并預(yù)測(cè)改流體對(duì)粉體下料流率的影響規(guī)律十分必要。然而目前對(duì)料倉(cāng)改流體的研究大多是關(guān)于調(diào)節(jié)流型的報(bào)道，缺乏對(duì)粉體下料流率變化的定量分析，更沒(méi)有成熟的理論模型能夠預(yù)測(cè)加裝改流體后的復(fù)雜流道結(jié)構(gòu)料倉(cāng)的粉體下料流率。

本文借助實(shí)驗(yàn)室有機(jī)玻璃料倉(cāng)，考察了傳統(tǒng)料倉(cāng)(無(wú)改流體)和加裝改流體料倉(cāng)(封閉改流體和開(kāi)放改流體)的粉體下料特性，研究了不同結(jié)構(gòu)的復(fù)雜流道對(duì)玻璃微珠、煤粉和聚氯乙烯三種粉體下料流率的作用規(guī)律；然后基于不同流道下粉體特有的流動(dòng)現(xiàn)象，并綜合考慮物料性質(zhì)以及料倉(cāng)結(jié)構(gòu)的影響，在對(duì)最小能量理論方程進(jìn)行修正的基礎(chǔ)上，建立包含粉體物性、下料流型和流道結(jié)構(gòu)影響的粉體下料流率預(yù)測(cè)模型。本研究對(duì)下料過(guò)程的強(qiáng)化與料倉(cāng)結(jié)構(gòu)優(yōu)化設(shè)計(jì)具有一定的指導(dǎo)意義。

1 實(shí)驗(yàn)物料及方法

1.1 實(shí)驗(yàn)物料

本文選取不同粒徑的玻璃微珠(glass bead-a、glass bead-b、glass bead-c)、褐煤(lignite)和聚氯乙烯顆粒(pvc)作為實(shí)驗(yàn)物料，通過(guò)氣動(dòng)篩分后(Retsch AS200 Jet Air)，得到Span 指數(shù)較小的窄粒徑分布的樣品，規(guī)避寬的、不規(guī)則的粒度分布對(duì)下料流率的影響。借助馬爾文激光粒度儀(Malvern 2000)對(duì)樣品的粒徑大小與粒度分布進(jìn)行測(cè)試，結(jié)果如圖1 所示，其中g(shù)lass bead-c、lignite 和pvc 的粒徑分布近似相同。表1給出了顆粒的表面積平均粒徑(d32)、體積平均粒徑(d43)和分布寬度Span指數(shù)。

圖1 粒徑累積分布Fig.1 Accumulative particle size distribution

表1 物性數(shù)據(jù)Table 1 Physical properties of experimental materials

借助FT4 粉體流變儀(Freeman technology)和PT-X 型粉體流動(dòng)性-噴流性測(cè)量?jī)x(Hosokawa Micron Corporation)對(duì)實(shí)驗(yàn)物料的流動(dòng)性質(zhì)和堆積性質(zhì)進(jìn)行了表征，得到了如表1 所示物料的堆積密度(ρb)、休止角(AOR)和內(nèi)摩擦角(AIF)等流動(dòng)性參數(shù)。以休止角作為表征物料流動(dòng)性的參數(shù)，可知，glass bead-b 的流動(dòng)性最好，屬于自由流動(dòng)粉體；lignite 的流動(dòng)性最差，屬于黏附性粉體。實(shí)驗(yàn)前在105℃下對(duì)物料進(jìn)行烘干處理，烘干后使用紅外水分測(cè)試儀(Sartorius MA150)測(cè) 量glass bead-a、glass bead-b、glass bead-c、lignite 的水分含量分別為0.02%、0.08%、0.08%和0.21%，可以忽略濕含量對(duì)顆粒流動(dòng)的影響。

1.2 實(shí)驗(yàn)方法

下料實(shí)驗(yàn)在有機(jī)玻璃料倉(cāng)中進(jìn)行，實(shí)驗(yàn)裝置如圖2 所示，打開(kāi)滑動(dòng)閥后，料倉(cāng)內(nèi)的物料從出口流出，沿斜槽滑落到接料盤(pán)中。整個(gè)下料過(guò)程的質(zhì)量變化由連接到計(jì)算機(jī)上的電子秤實(shí)時(shí)采集，精度為0.1 g，采樣頻率為8 Hz。從離散點(diǎn)繪制的下料曲線中可以得到物料的質(zhì)量流率。

圖2 下料實(shí)驗(yàn)裝置Fig.2 Diagram of discharge device

為了考察不同流道結(jié)構(gòu)料倉(cāng)內(nèi)的粉體流動(dòng)規(guī)律，在30°半錐角的有機(jī)玻璃料倉(cāng)中內(nèi)置了半錐角為15°的改流體，并且改流體分為封閉式和開(kāi)放式兩種。圖3 給出了料倉(cāng)、改流體的結(jié)構(gòu)參數(shù)及其裝配方式。改流體位于料倉(cāng)出口上端30 mm 處，如圖中虛線所示，兩者在料倉(cāng)出口處具有相同的口徑。

圖3 設(shè)備尺寸參數(shù)Fig.3 Device parameters

根據(jù)料倉(cāng)和改流體的裝配方式，本文研究了圖4 所示的三種流動(dòng)通道下粉體的下料規(guī)律，分別為無(wú)改流體料倉(cāng)(No-In)、封閉式改流體料倉(cāng)(Con-In)和開(kāi)放式改流體料倉(cāng)(Ucon-In)。無(wú)改流體料倉(cāng)是常用的錐形料倉(cāng)，整個(gè)錐體作為粉體的流動(dòng)通道；所謂封閉式改流體料倉(cāng)，意味著粉體無(wú)法從改流體內(nèi)部流過(guò)，只能在改流體和料倉(cāng)構(gòu)成的夾層間流動(dòng)；所謂開(kāi)放式改流體料倉(cāng)，表明粉體可以同時(shí)在改流體內(nèi)部、改流體和料倉(cāng)夾層間流動(dòng)。

圖4 下料操作模式(陰影部分為非流動(dòng)區(qū))Fig.4 Discharge operation modes(hatched domain is notavailable for flow)

2 結(jié)果與討論

2.1 復(fù)雜流道下料

圖5所示為粉體通過(guò)不同流道結(jié)構(gòu)的流動(dòng)示意圖。在No-In 的作用下，整個(gè)料倉(cāng)作為流動(dòng)通道，物料以中心流動(dòng)快、邊壁流動(dòng)慢的方式從料倉(cāng)出口卸出；引入改流體后，Con-In 作用下的物料流動(dòng)，在徑向方向上速度均勻性得到增強(qiáng)，在圖5 中表現(xiàn)為代表大速度梯度的深色陰影區(qū)面積減?。辉赨con-In的作用下，改流體內(nèi)部的物料先于夾層間的物料流出，存在流動(dòng)競(jìng)爭(zhēng)的特殊流動(dòng)現(xiàn)象。分析攝像拍攝的下料過(guò)程可知，兩個(gè)流動(dòng)通道內(nèi)的過(guò)程獨(dú)立進(jìn)行，不存在物料的相互混合。

Janssen 效應(yīng)的存在解釋了不隨時(shí)間變化的下料流率。本實(shí)驗(yàn)體系下，在No-In 和Con-In 的作用下，料倉(cāng)內(nèi)物料的剩余質(zhì)量隨時(shí)間均勻減少，不隨時(shí)間改變的下料流率，可通過(guò)線性擬合下料曲線獲得。在Ucon-In 的作用下，下料曲線存在轉(zhuǎn)折點(diǎn)，同樣采用線性擬合的方式，平均化處理兩段過(guò)程，得到下料流率。

表2給出了不同流道結(jié)構(gòu)下各種物料的下料流率。在所有的流道結(jié)構(gòu)下，三種glass bead 的下料流率接近，遠(yuǎn)大于lignite和pvc的下料流率。值得注意的是，盡管glass bead-c、lignite 和pvc 粒徑分布相當(dāng)，但三者的下料流率差異較大。glass bead-c 作為理想介質(zhì)，大的真密度、小的黏附性，重力驅(qū)動(dòng)下的流動(dòng)效果最好。對(duì)于lignite，一方面表面粗糙，孔隙結(jié)構(gòu)發(fā)達(dá)，顆粒間作用力較為復(fù)雜，黏聚性較強(qiáng)；另一方面，lignite 的大顆粒含量更多，大小顆粒鑲嵌流動(dòng)，不利于順暢下料。對(duì)于pvc，真密度小，下料過(guò)程重力驅(qū)動(dòng)的作用不顯著；同時(shí)，易摩擦帶電的性質(zhì)增加了顆粒間的靜電相互作用，不利于顆粒的流動(dòng)。

圖5 下料流動(dòng)順序Fig.5 Diagram of sequential flow order during discharge

表2 下料流率實(shí)驗(yàn)值Table 2 Experimental value of discharge rate

引入改流體調(diào)整流道結(jié)構(gòu)后，在Con-In 和Ucon-In 的作用下，物料的下料流率得到了增大。Con-In 作用下，物料的下料流率最大，促進(jìn)流動(dòng)的效果最為明顯。下料流率的提升程度如圖6 所示，性質(zhì)不同的物料，下料流率的提升程度也不盡相同，但Con-In 和Ucon-In 具有相同的影響趨勢(shì)。在Con-In 的作用下，glass bead-b 下料流率的提升幅度為最小的16%，lignite 的提升幅度為最大的58%。結(jié)合表1 可知，對(duì)休止角大、流動(dòng)性差的物料，流道結(jié)構(gòu)改變對(duì)下料流率的提升幅度更加明顯。

2.2 復(fù)雜流道下料流率預(yù)測(cè)

圖6 流道結(jié)構(gòu)對(duì)下料流率提升幅度的影響Fig.6 Relationship between increase of discharge rate and flow channel structure

盡管研究者們對(duì)粉體下料流率的建模研究已經(jīng)開(kāi)展了大量工作，但針對(duì)改流體料倉(cāng)，目前還沒(méi)有這方面的相關(guān)報(bào)道。本文從傳統(tǒng)料倉(cāng)的下料模型出發(fā)，綜合考慮粉體物性、下料流型和流道結(jié)構(gòu)的影響，最終建立具有普適性的復(fù)雜流道結(jié)構(gòu)的料倉(cāng)下料流量預(yù)測(cè)模型。

2.2.1 無(wú)改流體料倉(cāng)（No-In）下料模型修正 基于“free-fall arch”和“minimum energy theory”的原理，Brown 等[27]在圖7 所示的坐標(biāo)系統(tǒng)下，提出以下假設(shè)：①固體顆粒連續(xù)性方程中“free-fall arch”附近的物質(zhì)膨脹可以忽略；②在“free-fall arch”下方固體顆粒應(yīng)力σ的變化相對(duì)于其他部分對(duì)總能量的影響是可以忽略不計(jì)的；③假設(shè)顆粒徑向流動(dòng)；④“freefall arch”理想表面位置是假定在料倉(cāng)出料口邊界徑向坐標(biāo)r=r0處。

圖7 坐標(biāo)軸系統(tǒng)結(jié)構(gòu)Fig.7 Schematic diagram of coordinate system

基于這些假設(shè)，得到了最小能量理論方程式(1)，可用于預(yù)測(cè)理想物料以理想流動(dòng)狀態(tài)流出料倉(cāng)時(shí)的下料流率。

式中，Ws是下料流率，kg·s-1；α 是料倉(cāng)半錐角，（°）；ρb是物料的堆積密度，kg·m-3；g 是重力加速度，m·s-2；Do和dp分別是料倉(cāng)出口直徑和顆粒直徑，m；k是無(wú)量綱的Beverloo常數(shù)。

然而，一方面，受到顆粒形狀、粒徑分布和表面材質(zhì)等因素的影響，粉體的性質(zhì)會(huì)偏離理想物料。另一方面，細(xì)顆粒在料倉(cāng)中心附近膨脹快，邊壁附近膨脹慢，引發(fā)了料倉(cāng)內(nèi)的速度梯度。所以，“freefall arch”附近的物質(zhì)膨脹不再可以忽略，物料不再服從理想流動(dòng)。因此，如表3所示，式(1)對(duì)實(shí)際過(guò)程下料流率的預(yù)測(cè)產(chǎn)生了偏差，預(yù)測(cè)值偏大，對(duì)lignite的預(yù)測(cè)偏差達(dá)到了39.7%，不能正確地反映實(shí)驗(yàn)規(guī)律。下面依據(jù)實(shí)驗(yàn)中的流動(dòng)現(xiàn)象修正理論模型，將其拓展至實(shí)際下料過(guò)程。

Freyssingeas 等[6]描述了黏性顆粒下料過(guò)程中自由表面的波動(dòng)情況，圖8 所示為時(shí)間疊合下的物料主要流動(dòng)區(qū)的變化。下料過(guò)程中，物料并非以理想質(zhì)量流的方式流向出口，而是分層分級(jí)向出口流去。

圖8 下料過(guò)程物料自由表面的變化[6]Fig.8 Free surface during powder discharge process[6]

如圖9所示，本實(shí)驗(yàn)體系下，出現(xiàn)了同樣的流動(dòng)現(xiàn)象。glass bead-c 的自由表面存在凹陷，料倉(cāng)中心處的流動(dòng)速度快于壁面附近，物料并非以理想質(zhì)量流的方式卸出。因此，將料倉(cāng)中的流動(dòng)區(qū)域分為快速流動(dòng)區(qū)(fast flow zone)和剪切摩擦區(qū)(shear flow zone)。所謂的快速流動(dòng)區(qū)分布在料倉(cāng)中心處，顆粒運(yùn)動(dòng)速度快且分布均勻，粉體以理想質(zhì)量流的方式流向出口；所謂的剪切摩擦區(qū)則分布在料倉(cāng)邊壁附近，存在速度梯度。

圖9 玻璃微珠的下料過(guò)程Fig.9 Diagram of discharge of glass bead

摩擦特性的存在是實(shí)際物料偏離理想物料的的主要原因，物料性質(zhì)的差異造成了實(shí)際流動(dòng)與理想流動(dòng)的不一致。假定剪切摩擦區(qū)對(duì)下料流率沒(méi)有貢獻(xiàn)，借助實(shí)驗(yàn)值與式(1)所得理論值的比值來(lái)反映料倉(cāng)中剪切摩擦區(qū)的面積大小，lignite 剪切摩擦區(qū)的面積高達(dá)40%。式(1)對(duì)下料流率的預(yù)測(cè)偏大，可歸結(jié)為剪切摩擦區(qū)的影響。

以此為依據(jù)，對(duì)下料流率方程進(jìn)行修正。基于Rose 等[25]死區(qū)與水平方向夾角Φd的理念，以反映顆粒間的摩擦特性和抗剪特性的內(nèi)摩擦角φ近似作為快速流動(dòng)區(qū)與剪切摩擦區(qū)的劃分邊界，對(duì)式(1)進(jìn)行修正，引入如式(2)所示的校正因子F

表3 No-In下料模式下實(shí)驗(yàn)值與預(yù)測(cè)值的偏差Table 3 Deviation between experimental value and theoretical value under No-In mode

對(duì)于實(shí)際物料，受到剪切流動(dòng)的影響，下料流率低于理想流動(dòng)，所以校正因子F＜1。隨著內(nèi)摩擦角的增大，物料受剪切區(qū)的影響增大，下料流率減小，校正因子F也相應(yīng)減小。

將校正因子F 代入式(1)中，便得到下料流率預(yù)測(cè)公式

式中，Ws'是下料流率，kg·s-1;φ是物料的內(nèi)摩擦角，(°); α 是料倉(cāng)半錐角，(°); ρb是物料的堆積密度，kg·m-3；g 是重力加速度，m·s-2；Do和dp分別是料倉(cāng)出口直徑和顆粒直徑，m; k 是無(wú)量綱的Beverloo常數(shù)。

表3 給出了所有物料下料流率的實(shí)驗(yàn)值，式(1)計(jì)算所得的理論值，式(3)計(jì)算所得的修正值以及相對(duì)應(yīng)的偏差?？紤]了剪切摩擦區(qū)影響的式(3)提高了下料流率的預(yù)測(cè)精度。

如圖10所示，使用式(3)對(duì)文獻(xiàn)[28-30]中的下料流率值進(jìn)行預(yù)測(cè)。文獻(xiàn)涉及煤粉、玻璃微珠等多種物料，10°、20°等不同的錐角以及22.7、32.8 mm 不同的開(kāi)口直徑。下料流率覆蓋范圍從100 g·s-1左右到800 g·s-1左右。修正后計(jì)算公式的精度提高，并將理想的料倉(cāng)下料模型拓展至實(shí)際下料過(guò)程。

2.2.2 封閉式改流體料倉(cāng)（Con-In）下料模型修正

在料倉(cāng)中引入改流體，形成復(fù)雜的流道結(jié)構(gòu)，強(qiáng)化了流動(dòng)過(guò)程，提升了下料流率，具有潛在的應(yīng)用價(jià)值。準(zhǔn)確預(yù)測(cè)下料流率可以對(duì)強(qiáng)化下料過(guò)程與優(yōu)化料倉(cāng)結(jié)構(gòu)做出指導(dǎo)。

如圖11 所示，流道結(jié)構(gòu)從No-In 變?yōu)镃on-In后，流動(dòng)區(qū)域的發(fā)展衍變規(guī)律發(fā)生了改變。對(duì)于快速流動(dòng)區(qū)：①由料倉(cāng)中心向邊壁位置轉(zhuǎn)移，顆粒的徑向速度獲得了增大；②不對(duì)稱(chēng)壁面的剪切作用更加充分。這使得Con-In 流道內(nèi)的物料具有了更大的能量密度，快速流動(dòng)區(qū)得到更加充分的發(fā)展，有利于下料流率的提升。對(duì)于剪切摩擦區(qū)：狹小的體積使得剪切摩擦區(qū)的發(fā)展空間受到限制，無(wú)法充分衍變，相較之下，快速流動(dòng)區(qū)所占比重增大，物料的下料流率獲得提升。

圖10 實(shí)驗(yàn)值與預(yù)測(cè)值的比較Fig.10 Comparison of experiment and model values of discharge rate

圖11 流動(dòng)區(qū)域轉(zhuǎn)變示意圖Fig.11 Schematic diagram of flow zones transition

確定料倉(cāng)中剪切摩擦區(qū)的影響以及流道結(jié)構(gòu)的特征參數(shù)對(duì)提高預(yù)測(cè)公式的準(zhǔn)確性具有較大的意義。

由上文可知，式(2)中校正因子F的引入，可以較好地反映剪切摩擦區(qū)的影響。所以，預(yù)測(cè)Con-In 流道結(jié)構(gòu)下的下料流率時(shí)，仍然借用F 因子反映剪切摩擦區(qū)的影響。

流道結(jié)構(gòu)改變后，錐角是對(duì)結(jié)構(gòu)特征的直接反映，特征錐角的確定將尤為關(guān)鍵。圖12 分析了Con-In 作用下粉體的受力情況，微元粉體受到來(lái)自不對(duì)稱(chēng)壁面的支撐力和摩擦力以及重力的作用，水平方向的力平衡關(guān)系如式(4)所示

圖12 Con-In作用下粉體受力示意圖Fig.12 Schematic diagram of force analysis of powders under Con-In

式中，F(xiàn)i、Fh為壁面對(duì)粉體的摩擦力，N；Ni、Nh為壁面對(duì)粉體的支撐力，N；αi為改流體半錐角，αh為料倉(cāng)半錐角，(°)。

豎直方向的力平衡關(guān)系如式(5)所示

根據(jù)Coulomb-Mohr 定律[25]，不考慮物料的內(nèi)聚力，則

式中，μ是物料的摩擦系數(shù)。

將式(6)代入式(4)、式(5)，聯(lián)立求解，則

其中

由此可知，料倉(cāng)半錐角(αh)與改流體半錐角(αi)的差值決定了粉體在復(fù)雜流道中受到的支撐力與摩擦力的大小，是影響粉體運(yùn)動(dòng)的關(guān)鍵參數(shù)。借用此差值對(duì)原有模型中的錐角項(xiàng)進(jìn)行修正。

綜上，考慮剪切摩擦區(qū)以及料倉(cāng)結(jié)構(gòu)的影響，可借助式(10)對(duì)Con-In 流道結(jié)構(gòu)的下料流率進(jìn)行預(yù)測(cè)。

式中，特征錐角α'修正為料倉(cāng)半錐角(αh)與改流體半錐角(αi)之差。

表4給出了不同物料下料流率的實(shí)驗(yàn)值以及式(10)計(jì)算所得的修正值。結(jié)果表明，考慮了剪切摩擦區(qū)以及復(fù)雜流道結(jié)構(gòu)影響后的式(10)可以較好地預(yù)測(cè)下料流率。

2.2.3 開(kāi)放式改流體料倉(cāng)（Ucon-In）下料模型修正

相比于No-In，Ucon-In 可以提高物料的下料流率，Con-In 的作用效果更加明顯。對(duì)Ucon-In 而言，物料流動(dòng)通道的增加，并沒(méi)有帶來(lái)下料流率的增長(zhǎng)。由圖5 可知，在Ucon-In 的作用下，存在競(jìng)爭(zhēng)流動(dòng)導(dǎo)致的分階段下料現(xiàn)象。圖13給出了Ucon-In作用下glass bead-c 的瞬時(shí)下料流率，存在流率變化的轉(zhuǎn)折點(diǎn)，拐點(diǎn)前后下料流率平均值不同，第二階段的下料流率值與Con-In 作用下的下料流率數(shù)值吻合。低速下料的第一階段的存在，使得Ucon-In 流道下的下料流率小于Con-In 作用下的下料流率。因此，進(jìn)行流率預(yù)測(cè)時(shí)，需要充分考慮這一特殊的流動(dòng)現(xiàn)象。

表4 Con-In下料模式下實(shí)驗(yàn)值與預(yù)測(cè)值的偏差Table 4 Deviation between experimental value and theoretical value under Con-In mode

圖13 Ucon-In作用下的瞬時(shí)流率Fig.13 Diagram of variation of sample instantaneous discharge rate with time in silo with unconfined insert

從改流體內(nèi)部流出的物料受到出口直徑變化的影響，下料流率小于從夾層間流出的物料。理想情況下，物料在料倉(cāng)中的堆積密度分布均勻，計(jì)算得出質(zhì)量占比為m 的物料受到出口直徑變化的影響而導(dǎo)致流率減小。因此，可借助式(11)對(duì)Ucon-In流道結(jié)構(gòu)的下料流率進(jìn)行預(yù)測(cè)。

如圖14 所示，對(duì)不同物料，兩個(gè)階段下料流率的比值近似滿(mǎn)足線性關(guān)系，因此式(11)中f 因子可視為一個(gè)不受物性影響，與料倉(cāng)結(jié)構(gòu)有關(guān)的參數(shù)。

圖14 兩個(gè)階段下料流率的關(guān)系Fig.14 Relationship of discharge rate between two stages

由式(1)可知，下料流率與出口直徑呈2.5 次方，所以考慮了出口直徑影響的校正因子f 以式(12)的形式給出。

因此，便可獲得如式(13)所示的，Ucon-In 流道結(jié)構(gòu)作用下的下料流率預(yù)測(cè)公式。

式中，W?s是下料流率，kg·s-1；m 是改流體內(nèi)物料的質(zhì)量占比；Di、Dh分別為改流體和有機(jī)玻璃料倉(cāng)的出口直徑，m。

2.2.4 下料流率綜合預(yù)測(cè) 綜合上述三種不同料倉(cāng)結(jié)構(gòu)下的下料流動(dòng)現(xiàn)象和流率預(yù)測(cè)，對(duì)于傳統(tǒng)料倉(cāng)和復(fù)雜流道結(jié)構(gòu)料倉(cāng)中的顆粒流動(dòng)，可以使用式(14)進(jìn)行流率預(yù)測(cè)。

圖15 給出了復(fù)雜流動(dòng)通道作用下實(shí)驗(yàn)值和預(yù)測(cè)值的比較結(jié)果，預(yù)測(cè)公式符合實(shí)驗(yàn)規(guī)律，預(yù)測(cè)誤差在±10%以?xún)?nèi)。

3 結(jié) 論

本文研究了不同粉體在No-In、Con-In和Ucon-In 流道結(jié)構(gòu)作用下的流動(dòng)特性和流率規(guī)律，并分別建立模型對(duì)其下料流率進(jìn)行了預(yù)測(cè)，得到的主要結(jié)論如下。

圖15 實(shí)驗(yàn)值與預(yù)測(cè)值的比較Fig.15 Comparison of experiment and model values of discharge rate

(1)改流體具有顯著提升粉體下料流率的作用，且相比于Ucon-In，Con-In 促進(jìn)流動(dòng)的效果更佳。Ucon-In 流道結(jié)構(gòu)料倉(cāng)由于存在流動(dòng)競(jìng)爭(zhēng)現(xiàn)象，盡管實(shí)際流動(dòng)通道較大，但粉體的下料流率低于Con-In 流道結(jié)構(gòu)料倉(cāng)。對(duì)于黏附性煤粉，Con-In 對(duì)下料流率提升幅度可高達(dá)58%，表明加裝封閉式改流體是改善黏附性粉體下料的有效手段之一。

(2)物料的摩擦性質(zhì)導(dǎo)致了料倉(cāng)實(shí)際下料過(guò)程偏離理想流動(dòng)，使得用于預(yù)測(cè)粉體下料流率的傳統(tǒng)模型出現(xiàn)偏差。本文定量計(jì)算的剪切摩擦區(qū)面積解釋了偏差存在的原因，并進(jìn)一步提出流率校正因子，綜合考慮物料性質(zhì)和區(qū)域流動(dòng)的作用，對(duì)最小能量理論方程進(jìn)行了修正，將理想的料倉(cāng)下料模型拓展至實(shí)際下料過(guò)程。

(3)針對(duì)加裝改流體后的復(fù)雜流道結(jié)構(gòu)料倉(cāng)，分別建模分析了Con-In和Ucon-In流道結(jié)構(gòu)下粉體的下料流率。對(duì)于Con-In，受到不對(duì)稱(chēng)壁面的剪切作用，料倉(cāng)半錐角無(wú)法反映倉(cāng)壁的影響，引入了特征參數(shù)α'對(duì)下料模型進(jìn)行修正；對(duì)于Ucon-In，基于粉體流動(dòng)競(jìng)爭(zhēng)機(jī)制和內(nèi)外流道內(nèi)的流率關(guān)系對(duì)下料模型進(jìn)行修正。在上述基礎(chǔ)上，建立了綜合考慮了粉體物性、下料流型和流道結(jié)構(gòu)影響的下料流率預(yù)測(cè)模型，預(yù)測(cè)偏差＜10%。