許冰洋，朱文魁，潘廣樂(lè)，申玉軍，李 斌，于川芳

中國(guó)煙草總公司鄭州煙草研究院，鄭州高新技術(shù)產(chǎn)業(yè)開發(fā)區(qū)楓楊街2號(hào) 450001

基于收縮特性分析的葉絲快速對(duì)流干燥動(dòng)力學(xué)模型

許冰洋，朱文魁，潘廣樂(lè)，申玉軍，李 斌，于川芳*

中國(guó)煙草總公司鄭州煙草研究院，鄭州高新技術(shù)產(chǎn)業(yè)開發(fā)區(qū)楓楊街2號(hào) 450001

為準(zhǔn)確表征葉絲脫水過(guò)程中的干燥動(dòng)力學(xué)特性，分析了烤煙和白肋煙兩種葉絲在下行床快速對(duì)流干燥中的孔隙結(jié)構(gòu)變化特征，對(duì)葉絲干燥收縮過(guò)程進(jìn)行了數(shù)學(xué)模型擬合，建立了考慮收縮形變的葉絲干燥過(guò)程水分?jǐn)U散模型，并對(duì)模型進(jìn)行了驗(yàn)證。結(jié)果表明：①隨著干燥過(guò)程中兩種葉絲內(nèi)孔容積的減小，葉絲中孔徑大于0．5μm的孔容比例呈降低趨勢(shì)，而孔徑小于0．05μm的孔容比例呈升高趨勢(shì)；在不同干燥階段兩種葉絲孔徑分布的分形維數(shù)介于2．45～2．71之間，表明其孔隙結(jié)構(gòu)具有分形特征；②線性疊加式收縮模型能夠較好地描述葉絲干燥過(guò)程中的收縮現(xiàn)象，烤煙葉絲和白肋煙葉絲體積比V/V0和含水率比X/X0的線性相關(guān)系數(shù)均大于0．99，采用該收縮模型對(duì)基于Fick第二定律的水分?jǐn)U散模型進(jìn)行修正，實(shí)驗(yàn)數(shù)據(jù)的擬合精度從修正前的0．906 9提升到修正后的0．958 0；③采用修正的水分?jǐn)U散模型描述葉絲快速干燥動(dòng)力學(xué)發(fā)現(xiàn)，考慮了干燥過(guò)程中的葉絲收縮現(xiàn)象后，得到的葉絲水分有效擴(kuò)散系數(shù)降低，表明葉絲干燥過(guò)程中的體積收縮不利于傳質(zhì)過(guò)程。

葉絲；對(duì)流干燥；收縮；動(dòng)力學(xué)模型

干燥脫水是煙草物料加工中的關(guān)鍵熱濕處理環(huán)節(jié)，貫穿了煙葉原料從采后初烤、打葉復(fù)烤到制絲的整個(gè)加工流程。與已報(bào)道的果蔬、種子、木材等植物基多孔材料一致［1］，片煙、葉絲等煙草物料在干燥加工過(guò)程中，由于溫濕度變化引起的應(yīng)力作用，均存在物理結(jié)構(gòu)上的收縮形變。對(duì)于葉絲干燥過(guò)程，這一物理結(jié)構(gòu)上的變化一方面可能影響最終煙支的燃燒特性［2］；另一方面，由于干燥過(guò)程中水分主要通過(guò)物料內(nèi)部的孔隙進(jìn)行擴(kuò)散和遷移［3］，物理結(jié)構(gòu)上的變化也會(huì)顯著影響其干燥動(dòng)力學(xué)。關(guān)于物理結(jié)構(gòu)收縮對(duì)多孔介質(zhì)干燥過(guò)程傳質(zhì)、傳熱的影響，在食品、果蔬等領(lǐng)域已有研究報(bào)道。Brasiello等［4］研究了收縮對(duì)茄子干燥過(guò)程中熱質(zhì)傳遞過(guò)程的影響，通過(guò)收縮系數(shù)對(duì)干燥動(dòng)力學(xué)模型進(jìn)行修正，能更真實(shí)地描述茄子干燥過(guò)程。Dissa等［5］以芒果片為材料，考慮了體積收縮對(duì)其干燥動(dòng)力學(xué)的影響，發(fā)現(xiàn)材料的體積收縮不利于材料干燥過(guò)程中的傳質(zhì)。目前，對(duì)于煙草原料對(duì)流干燥動(dòng)力學(xué)的研究，相關(guān)文獻(xiàn)［6-8］較多采用基于Fick第二定律的擴(kuò)散模型來(lái)描述煙草干燥動(dòng)力學(xué)，所建立的動(dòng)力學(xué)模型中均假定被干燥物料脫水過(guò)程中保持剛性，即忽略了葉絲脫水干燥過(guò)程中的體積收縮形變，也未考慮物料物理結(jié)構(gòu)的變化對(duì)干燥脫水過(guò)程的影響，難以準(zhǔn)確反映葉絲在實(shí)際脫水過(guò)程中的干燥動(dòng)力學(xué)特性?；诖?，在探討兩種葉絲快速對(duì)流干燥過(guò)程收縮特性的基礎(chǔ)上，建立了考慮物料物理結(jié)構(gòu)收縮的葉絲干燥動(dòng)力學(xué)模型，旨在為煙草物料干燥動(dòng)力學(xué)特性的分析和評(píng)價(jià)提供更為準(zhǔn)確的模型和方法依據(jù)。

1 材料與方法

1.1 材料和儀器

以2012年河南平頂山烤煙C3F和2011年重慶白肋煙C3F為原料。片煙切絲后，在恒溫恒濕間將葉絲含水率調(diào)節(jié)至30%待用。

ALPHA 1-4 LD plus型真空冷凍干燥機(jī)（德國(guó)Christ公司）；DHG-9623型烘箱（上海精宏試驗(yàn)設(shè)備有限公司）；BSA124S-CW電子天平（感量：0.000 1 g，德國(guó)Sartorius公司）；UltraPYC 1200e真密度儀（美國(guó)Quantachrome公司）；AutoPore IV 9500型壓汞儀（美國(guó)Micromeritics公司）。

1.2 測(cè)試方法

采用鄭州煙草研究院自行設(shè)計(jì)的下行床干燥試驗(yàn)裝置進(jìn)行葉絲干燥實(shí)驗(yàn)，裝置具體結(jié)構(gòu)及試驗(yàn)方法見(jiàn)文獻(xiàn)［6］。不同干燥階段煙絲樣品含水率采用烘箱法測(cè)定［7］。將干燥后得到的葉絲樣品分為兩部分，一部分根據(jù)行業(yè)標(biāo)準(zhǔn)方法YC/T 473—2013［8］進(jìn)行表觀密度和真密度測(cè)試；另一部分經(jīng)過(guò)真空冷凍干燥后，用于壓汞實(shí)驗(yàn)以得到樣品的孔徑分布數(shù)據(jù)。壓汞測(cè)試壓力最高為414 MPa，浸潤(rùn)角取130°，采用儀器自帶軟件進(jìn)行孔徑分布數(shù)據(jù)處理。

2 結(jié)果與討論

2.1 干燥過(guò)程中葉絲孔隙結(jié)構(gòu)的變化特征分析

2.1.1 葉絲內(nèi)孔容積的變化

葉絲屬膠體毛細(xì)管多孔介質(zhì)，干燥過(guò)程中極易發(fā)生收縮變形導(dǎo)致葉絲內(nèi)部孔隙結(jié)構(gòu)發(fā)生改變。如表1所示，在下行床快速干燥過(guò)程中，烤煙和白肋煙葉絲的內(nèi)孔容積總體上均隨含水率的降低呈減小趨勢(shì)。與白肋煙葉絲相比，烤煙葉絲的毛細(xì)管收縮應(yīng)力作用所導(dǎo)致的內(nèi)孔容積減小更為顯著。

表1 下行床干燥過(guò)程中葉絲內(nèi)孔容積的變化

根據(jù)壓汞實(shí)驗(yàn)結(jié)果，進(jìn)一步分析了葉絲干燥過(guò)程中孔容在不同孔徑區(qū)間分布的變化。對(duì)初始樣品及干燥至含水率25%、15%和10%的3個(gè)不同脫水階段的葉絲，統(tǒng)計(jì)分析了＞20μm、0.5～20μm、0.05～0.50μm以及＜0.05μm等4個(gè)孔徑范圍內(nèi)的孔容比例，結(jié)果見(jiàn)圖1。

圖1 干燥過(guò)程中葉絲各孔徑區(qū)間內(nèi)對(duì)應(yīng)的孔容比例

由圖1可以看出：①對(duì)烤煙葉絲，干燥至含水率10%時(shí)孔徑大于20μm的大孔容積比例明顯降低，0.05～0.50μm的孔容比例也有所降低，而孔徑小于0.5μm的細(xì)微孔容積比例則明顯增加。對(duì)白肋煙葉絲，干燥過(guò)程中隨著含水率的不斷降低，孔徑大于20μm的大孔容積比例呈先升高后降低趨勢(shì)，而孔徑小于0.5μm的細(xì)微孔容積比例增加明顯。葉絲干燥過(guò)程中一方面由于毛細(xì)管的收縮應(yīng)力使葉絲內(nèi)部部分大孔及中孔孔隙會(huì)收縮成為小孔，另一方面由于葉絲含水率降低，脆性加大，葉絲內(nèi)大孔及中孔孔隙出現(xiàn)坍塌、閉合現(xiàn)象，導(dǎo)致干燥后期烤煙葉絲和白肋煙葉絲在孔徑小于0.5μm范圍內(nèi)的孔容比例有所增大。②與烤煙相比，白肋煙中未見(jiàn)孔徑小于0.05μm的細(xì)微孔，說(shuō)明白肋煙葉絲孔徑分布范圍較窄，這與前期工作［9］的分析結(jié)果是一致的。③烤煙及白肋煙葉絲孔隙中均以大孔所占孔容比例較高，干燥過(guò)程中大孔容積的變化決定了葉絲整體內(nèi)孔容積的變化趨勢(shì)，其中白肋煙葉絲中孔徑大于20μm的大孔容積比例的變化趨勢(shì)也與前期研究工作［11］中白肋煙葉絲干燥過(guò)程中內(nèi)孔容積出現(xiàn)的階段性上升現(xiàn)象相吻合。

2.1.2 葉絲孔隙結(jié)構(gòu)的分形表征

為了進(jìn)一步定性表征葉絲干燥過(guò)程中孔隙結(jié)構(gòu)復(fù)雜程度的變化，計(jì)算了兩種葉絲不同干燥階段的分形維數(shù)。對(duì)于孔隙具有分形特征的多孔物料，孔隙結(jié)構(gòu)的分形維數(shù)D一般介于2～3之間，分形維數(shù)越大表明物料孔隙表面在空間的形貌特性偏離光滑表面的程度也越遠(yuǎn)，孔隙結(jié)構(gòu)復(fù)雜程度越高［10］。根據(jù)文獻(xiàn)［11］所述方法，計(jì)算初始樣品及25%、15%和10%等4個(gè)不同含水率階段葉絲分形維數(shù)，結(jié)果見(jiàn)表2。由表2可知，兩種葉絲的孔隙分形維數(shù)處于2～3之間，與文獻(xiàn)［12］報(bào)道馬鈴薯、凍干牛肉等生物基多孔材料類似，具有明顯的分形特征。對(duì)多孔介質(zhì)材料，分形維數(shù)能夠反映孔隙結(jié)構(gòu)的復(fù)雜程度，包括材料的孔徑分布、孔隙表面的粗糙度及孔隙形狀規(guī)則程度等方面。根據(jù)本實(shí)驗(yàn)結(jié)果可知，烤煙和白肋煙兩種不同材料葉絲的分形維數(shù)存在一定差異，但同種材料在干燥過(guò)程中分形維數(shù)變化不明顯。這是由于干燥中材料內(nèi)部的孔徑分布雖然發(fā)生一定程度的變化，但與材料本身特性相關(guān)的孔隙表面的粗糙度及孔形狀規(guī)則程度等因素，在干燥過(guò)程中變化不大，導(dǎo)致反映孔隙結(jié)構(gòu)復(fù)雜度的分形維數(shù)變化并不明顯。Macías-García等［13］在研究不同材料制成的活性炭的分形特征時(shí)也發(fā)現(xiàn)，不同材料制成的活性炭分形維數(shù)差異較大，而同種材料制成的活性炭分形維數(shù)波動(dòng)范圍較小，這也說(shuō)明材料本身的結(jié)構(gòu)特性是決定多孔介質(zhì)分形維數(shù)的關(guān)鍵因素。

表2 葉絲孔隙結(jié)構(gòu)的分形維數(shù)

2.2 葉絲干燥動(dòng)力學(xué)模型分析

2.2.1 葉絲干燥收縮模型表征

干燥過(guò)程中葉絲內(nèi)部在溫濕度變化所引起的毛細(xì)管收縮應(yīng)力等作用下，內(nèi)部孔隙結(jié)構(gòu)會(huì)發(fā)生改變并影響干燥過(guò)程。由于多孔介質(zhì)內(nèi)水分一般在孔隙內(nèi)傳輸，葉絲干燥過(guò)程中的孔隙結(jié)構(gòu)變化會(huì)導(dǎo)致水分在葉絲內(nèi)部的遷移路徑發(fā)生改變，對(duì)葉絲干燥過(guò)程中水分?jǐn)U散過(guò)程產(chǎn)生影響。由于含濕多孔介質(zhì)的體積主要由固體骨架體積和孔隙中水的體積組成，則干燥特定階段的葉絲體積比可由式（1）確定［5］：

式中：ms—干基質(zhì)量，g；ρs—葉絲干燥過(guò)程固體真密度，g/cm3；ρs0—葉絲干燥前固體真密度，g/cm3；ρw—水的密度，g/cm3；V0—葉絲初始體積，cm3；V—干燥過(guò)程中葉絲體積，cm3；X0—葉絲初始含水率，%；X—葉絲干燥過(guò)程中的含水率，%。

根據(jù)不同干燥階段的葉絲真密度、含水率以及水的密度，由式（1）即可確定該階段的葉絲體積比的試驗(yàn)值。葉絲的體積比V/V0即反映了葉絲干燥過(guò)程中的收縮程度。其中，式（1）中水的密度ρw取定值為1.000 g/cm3。表3為不同干燥階段葉絲含水率及真密度測(cè)試結(jié)果。

表3 下行床干燥過(guò)程中葉絲含水率及真密度的變化

植物性多孔介質(zhì)干燥過(guò)程中的收縮模型已有較多研究。Mayor等［14］總結(jié)了不同種類多孔介質(zhì)的收縮模型，分為線性、非線性的經(jīng)驗(yàn)?zāi)Ｐ秃途€性、非線性的理論模型。Vacarezza［15］在研究收縮對(duì)甜菜根干燥過(guò)程的影響時(shí)選用了薄片狀的甜菜根為材料，材料的物理特征與葉絲相近，所以本實(shí)驗(yàn)中使用文獻(xiàn)［16］中提出的疊加式線性理論模型，如式（2）所示，對(duì)干燥過(guò)程中葉絲的體積收縮現(xiàn)象進(jìn)行描述：

式中：p1，p2—常數(shù)。

干燥過(guò)程中，在不同含水率階段，兩種葉絲體積比V/V0的試驗(yàn)結(jié)果見(jiàn)圖2所示。同時(shí)，根據(jù)式（2）所述的線性收縮模型，在圖2中對(duì)不同干燥階段的葉絲體積比V/V0和含水率比X/X0進(jìn)行了線性擬合。可以看出，兩種葉絲的體積比在干燥過(guò)程中均呈下降趨勢(shì)，即葉絲體積收縮程度隨含水率的降低而增大。不同干燥階段的葉絲體積比V/V0和含水率比X/X0的線性相關(guān)系數(shù)均大于0.99，表明此線性模型能夠較好地描述葉絲干燥過(guò)程中的收縮現(xiàn)象。在下行床干燥過(guò)程中烤煙葉絲的干燥溫度較白肋煙葉絲高，而烤煙葉絲的體積收縮程度仍低于白肋煙，說(shuō)明葉絲本身的物性對(duì)其收縮程度影響較大。

2.2.2 葉絲水分?jǐn)U散模型的修正

圖2 葉絲體積收縮程度隨含水率的變化情況

目前，對(duì)于滿足薄層干燥條件的片煙、葉絲等煙草物料對(duì)流干燥過(guò)程，通常采用基于水分?jǐn)U散理論的Fick第二定律來(lái)分析干燥動(dòng)力學(xué)［16-17］，其數(shù)學(xué)描述如式（3）所示。

式中：M—擴(kuò)散介質(zhì)的濃度，kg/m3；t—干燥時(shí)間，s；x—擴(kuò)散距離，m；De—有效擴(kuò)散系數(shù)，m2/s。

采用該理論模型對(duì)薄層干燥動(dòng)力學(xué)分析過(guò)程中，普遍將被干燥物料視為剛性體，即不考慮物料的收縮形變，進(jìn)而通過(guò)式（4）所示解析解［18］的形式，求解得到物料的有效擴(kuò)散系數(shù)。擴(kuò)散系數(shù)通常用于反映物料在一定干燥條件下的脫水能力。

式中：Xeq—平衡含水率，%；L—葉絲寬度的1/2，m；n—級(jí)數(shù)項(xiàng)數(shù)。

基于上述剛性物料的假設(shè)，目前在式（4）解析解方程中均將物料尺寸參數(shù)L視為常量。顯然，這一假定與實(shí)際干燥過(guò)程中存在的收縮現(xiàn)象有較大差異。由于干燥過(guò)程中葉絲產(chǎn)生的體積收縮會(huì)導(dǎo)致葉絲內(nèi)部孔隙的形態(tài)、尺寸發(fā)生改變，必然會(huì)影響干燥過(guò)程中水分的傳遞過(guò)程。因此，考慮葉絲體積收縮對(duì)傳質(zhì)過(guò)程的影響，假設(shè)水分在葉絲內(nèi)的傳遞為一維擴(kuò)散，在與葉絲切面垂直方向上進(jìn)行單向傳質(zhì)，切絲寬度為1 mm，則L可取葉絲寬度的一半，即5×10-4m。假設(shè)干燥過(guò)程中葉絲的體積收縮為一維收縮，且收縮方向與水分傳遞方向保持一致，則可得到傳質(zhì)路徑L與含水率X的關(guān)系式：

由上式可知，干燥過(guò)程中傳質(zhì)路徑L為含水率X的函數(shù)，含水實(shí)驗(yàn)中干燥介質(zhì)為絕干空氣，相對(duì)濕度近似為零，該條件下葉絲最終的平衡含水率很小，可以忽略，則式（4）可整理為：

因此，由式（5）和（6）結(jié)合，即可在考慮物料實(shí)際收縮的條件下對(duì)煙草干燥過(guò)程中的水分變化特征進(jìn)行描述，并獲取相應(yīng)的有效擴(kuò)散系數(shù)。由實(shí)驗(yàn)室前期研究［16］可知，當(dāng)n≥3時(shí)，有效擴(kuò)散系數(shù)計(jì)算值的相對(duì)誤差趨于穩(wěn)定，故取n=3。

2.3 葉絲干燥動(dòng)力學(xué)模型的驗(yàn)證

由式（6）可得到修正后葉絲干燥過(guò)程中的有效擴(kuò)散系數(shù)?？紤]葉絲體積收縮對(duì)傳質(zhì)過(guò)程的影響前后對(duì)比結(jié)果見(jiàn)圖3。如圖3所示，考慮體積收縮對(duì)葉絲干燥過(guò)程的影響并對(duì)擴(kuò)散方程進(jìn)行修正后，得到的理論值與實(shí)測(cè)值的相關(guān)系數(shù)從修正前的0.906 9提高到0.958 0。

圖3 考慮葉絲體積收縮對(duì)傳質(zhì)過(guò)程的影響前后對(duì)比

葉絲體積收縮對(duì)傳質(zhì)過(guò)程中水分?jǐn)U散系數(shù)的影響見(jiàn)表4。由表4可知，考慮干燥過(guò)程中葉絲體積收縮對(duì)傳質(zhì)過(guò)程的影響，烤煙葉絲的水分有效擴(kuò)散系數(shù)由2.97×10-8m2/s降至2.12×10-8m2/s，對(duì)白肋煙葉絲，水分有效擴(kuò)散系數(shù)由3.99×10-8m2/s降至2.94×10-8m2/s?？梢园l(fā)現(xiàn)，修正后葉絲的有效擴(kuò)散系數(shù)明顯降低，說(shuō)明干燥過(guò)程中葉絲的體積收縮行為不利于傳質(zhì)過(guò)程，與已報(bào)道的收縮對(duì)芒果片干燥過(guò)程影響的研究結(jié)果一致［5］。

表4 葉絲體積收縮對(duì)傳質(zhì)過(guò)程中水分?jǐn)U散系數(shù)的影響

葉絲內(nèi)部孔隙結(jié)構(gòu)復(fù)雜，干燥過(guò)程中水分從葉絲內(nèi)部向表面擴(kuò)散遷移，但并不能在任意方向上擴(kuò)散，孔隙結(jié)構(gòu)是決定水分傳遞過(guò)程的關(guān)鍵因素。葉絲干燥過(guò)程中發(fā)生的體積收縮現(xiàn)象，使內(nèi)部水分傳質(zhì)路徑在收縮方向上發(fā)生變化，但葉絲的收縮導(dǎo)致內(nèi)部孔隙結(jié)構(gòu)發(fā)生改變，使無(wú)序排列的孔隙更加復(fù)雜，內(nèi)部水分在向葉絲表面遷移擴(kuò)散過(guò)程中受到更多阻滯，導(dǎo)致葉絲內(nèi)水分的有效擴(kuò)散系數(shù)降低，從而也說(shuō)明干燥過(guò)程中葉絲發(fā)生的體積收縮現(xiàn)象不利于傳質(zhì)過(guò)程。

3 結(jié)論

（1）干燥過(guò)程中兩種葉絲的內(nèi)孔容積呈降低趨勢(shì)，葉絲在大于0.5μm孔徑范圍內(nèi)的孔容比例呈降低趨勢(shì)，小于0.05μm孔徑范圍內(nèi)的孔容比例呈升高趨勢(shì)，葉絲的體積收縮行為對(duì)葉絲的孔隙結(jié)構(gòu)分布具有顯著影響。在不同干燥階段兩種葉絲孔徑分布的分形維數(shù)介于2.45～2.71之間，表明其孔隙結(jié)構(gòu)具有分形特征。

（2）Vacarezza提出的線性收縮模型能夠較好描述葉絲干燥過(guò)程中收縮現(xiàn)象，葉絲體積比V/V0和含水率比X/X0的線性相關(guān)系數(shù)均大于0.99。采用該收縮模型對(duì)基于Fick第二定律的水分?jǐn)U散模型進(jìn)行修正，實(shí)驗(yàn)數(shù)據(jù)的擬合精度從修正前的0.906 9提升到修正后的0.958 0，得到了能夠更準(zhǔn)確描述葉絲實(shí)際干燥過(guò)程的動(dòng)力學(xué)模型。

（3）采用修正的水分?jǐn)U散模型描述葉絲快速干燥動(dòng)力學(xué)發(fā)現(xiàn)，考慮實(shí)際收縮條件下得到的葉絲水分有效擴(kuò)散系數(shù)較未考慮收縮時(shí)低26%～28%。這也表明葉絲干燥過(guò)程中的體積收縮不利于傳質(zhì)過(guò)程。

［1］LIU Guohong，CHEN Junruo，LIU Meihong，etal. Shrinkage，porosity and density behaviour during convective drying of bio-porous material［J］.Procedia Engineering，2012，31：634-640.

［2］李巧靈，陳曉東，李軍，等.卷煙燃燒模型研究進(jìn)展［J］.中國(guó)煙草學(xué)報(bào)，2013，19（2）：115-122.

［3］袁越錦，楊彬彬，焦陽(yáng)，等.多孔介質(zhì)干燥過(guò)程分形孔道網(wǎng)絡(luò)模型與模擬：Ⅰ.模型建立［J］.中國(guó)農(nóng)業(yè)大學(xué)學(xué)報(bào)，2007，12（3）：65-69.

［4］Brasiello A，Adiletta G，Russo P，etal.Mathematical modeling of eggplant d rying：Shrinkage effect［J］. Journal of Food Engineering，2013，114（1）：99-105.

［5］Dissa A O，Desmorieux H，Bathiebo J，etal.Convective drying characteristics of Amelie mango（Mangifera indicaL.cv.‘Amelie’）w ith correction for shrinkage［J］. Journal of Food Engineering，2008，88（4）：429-437.

［6］王鵬飛，劉朝賢，朱文魁，等.在氣固并流下行床中煙絲的快速干燥動(dòng)力學(xué)分析［J］.煙草科技，2014（4）：5-9.

［7］YC/T 31—1996煙草及煙草制品試樣的制備和水分的測(cè)定烘箱法［S］.

［8］YC/T 473—2013煙絲表觀密度、真密度和內(nèi)孔容積的測(cè)定［S］.

［9］許冰洋，朱文魁，申玉軍，等.不同干燥方式下葉絲孔隙結(jié)構(gòu)的變化特征分析［J］.煙草科技，2015，48（4）：60-65.

［10］Mandelbrot B B.Fractals：Form，chance and dimension［M］.San Francisco：Freeman，1977.

［11］ZHANG Baoquan，LI Shaofen.Determ ination of the surface fractal dimension for porous media by mercury porosimetry［J］.Industrial&Engineering Chem istry Research，1995，34（4）：1383-1386.

［12］劉永忠，陳三強(qiáng)，孫皓，等.凍干物料孔隙特性表征的分形模型與分形維數(shù)［J］.農(nóng)業(yè)工程學(xué)報(bào)，2004，20（6）：41-45.

［13］Macías-García A，Díaz-Díez MA，Cuerda-Correa E M，etal.Study of the pore size distribution and fractal dimension of HNO3-treated activated carbons［J］.Applied Surface Science，2006，252（17）：5972-5975.

［14］Mayor L，Sereno A M.Modelling shrinkage during convective drying of food materials：A review［J］. Journal of Food Engineering，2004，61（3）：373-386.

［15］Vacarezza L.Cinetica y mecanismo de transporte de agua durante la deshidrataci on de la remolachaazucarera［D］.Buenos Aires：Universidad de Buenos Aires，1975.

［16］張合明，閆亞明，劉朝賢，等.壓力對(duì)煙絲干燥有效擴(kuò)散系數(shù)的影響［J］.煙草科技，2011（6）：5-8.

［17］劉澤，李斌，于川芳.復(fù)合傳熱對(duì)薄層煙絲干燥強(qiáng)化傳質(zhì)作用的動(dòng)力學(xué)分析［J］.煙草科技，2009（11）：5-10.

［18］Crank J.The mathematics of diffusion［M］.London：Oxford University Press，1975.

責(zé)任編輯 周雅寧

Kinetic Model for Rapid Convective Drying of Cut Strips Based on Shrinkage Property Analysis

XU Bingyang,ZHU Wenkui,PAN Guangle,SHEN Yujun,LI Bin,and YU Chuanfang*
Key Laboratory of Tobacco Processing Technology of CNTC,Zhengzhou Tobacco Research Institute of CNTC, Zhengzhou 450001,China

In order to accurately characterize the kinetic characteristics of cut strips inthe process of dehydration,the variation of pore structure in cut strips of flue-cured and burley tobaccos during rapid convective drying in a downflowbed was analyzed.The process o f cut tobacco drying and shrinking was fitted by a mathematic model,and the model for moisture diffusion during drying taking into account the shrinking deformation was developed and verified.The results showed that：1)W ith the decrease of pore volume of cut strips of flue-cured and burley tobaccos during drying,the proportion of volume of pores of＞0.5μm in size decreased,while that of＜0.05μm increased.The fractal dimension of the two tobaccos ranged from 2.45 to 2.71 at different drying stages,which indicated that their pore structure possessed fractal characteristics.2)The linear superposition model described the shrinkage of cut strips during drying pretty well,the linear correlation coefficients for the volume ratio(V/V0)and moisture content ratio(X/X0)between cut strips of the two tobaccos were higher than 0.99.The shrinkage model was used to revise the moisture diffusion model based on Fick’s second law,so that the fitting accuracy ofexperimental data increased from 0.906 9 to 0.958 0.3)From the revised moisture diffusion model,it was found that the effective diffusivity of moisture was overestimated when the shrinkage of cut strips during drying was neglected,which indicated that the shrinkage of cut strips impacted a negative influence on mass transfer p rocess during d rying.

Cut strip;Convective drying;Shrinkage;Kinetic model

TS412

A

1002-0861（2015）09-0069-06

10.16135/j.issn1002-0861.20150912

2014-11-07

2015-06-11

國(guó)家自然科學(xué)基金項(xiàng)目“植物多孔介質(zhì)快速干燥過(guò)程微細(xì)結(jié)構(gòu)變化的分形特征及對(duì)熱質(zhì)傳遞影響機(jī)制”（51306213）。

許冰洋（1990—），在讀碩士研究生，研究方向：食品科學(xué)/煙草工藝。E-mail：hedaxby@163.com；*

于川芳，E-mail：yuchuanfang0208@126.com

許冰洋，朱文魁，潘廣樂(lè)，等.基于收縮特性分析的葉絲快速對(duì)流干燥動(dòng)力學(xué)模型［J］.煙草科技，2015，48（7）：69-74.

XU Bingyang,ZHU Wenkui,PAN Guangle,etal.Kinetic model for rapid convective drying of cut strips based on shrinkage property analysis［J］.Tobacco Science&Technology，2015，48（9）：69-74.