種翠娟,朱文學,劉云宏,吳建業(yè),樊金玲,羅 磊
胡蘿卜薄層干燥動力學模型研究
種翠娟,朱文學*,劉云宏,吳建業(yè),樊金玲,羅 磊
(河南科技大學食品與生物工程學院,河南 洛陽 471023)
為探索胡蘿卜熱風干燥過程中水分的變化規(guī)律,本研究以胡蘿卜為干燥對象,進行薄層干燥特性及模型研究,探討不同溫度、風速及物料厚度條件下胡蘿卜水分比與干燥時間的關系,建立動力學模型;以Fick擴散定律為依據(jù),確定胡蘿卜一維傳熱傳質的有效水分擴散系數(shù)并建立其數(shù)學模型。結果表明:胡蘿卜薄層干燥動力學模型可用Page方程來描述,并通過回歸分析確定方程系數(shù)m、k,通過多元線性回歸方法得到有效水分擴散系數(shù)(Deff)與溫度、風速和厚度的表達式,實驗得到的Deff值在0.84×10-9~6.69×10-9m2/s范圍內隨著干燥溫度、風速和物料厚度的升高而增大。
胡蘿卜;干燥;動力學模型;水分擴散系數(shù)
胡蘿卜含有豐富的營養(yǎng)成分、無機鹽類及氨基酸等[1-2],具有很高的營養(yǎng)價值和多種保健功效,被譽為“土人參”[3-4],深受消費者喜愛。胡蘿卜的干燥方式多種多樣,但熱風干燥仍是目前普遍使用的干燥技術[4-6]。熱風干燥作為傳統(tǒng)的干燥方法,是用加熱后的空氣作媒介,將物料進行加熱以促進水分的蒸發(fā),具有熱效率高,干燥速率快,設備結構簡單、投資少、適應性強、操作和控制簡單,衛(wèi)生條件較好等優(yōu)勢。但由于物料長時間與高溫空氣接觸,容易導致干制品色澤褐變及營養(yǎng)成分損失[7-8]。為有效提高干燥產(chǎn)品的品質,優(yōu)化干燥工藝,有必要使用數(shù)學模型來模擬不同干燥條件下的干燥過程。目前,常用半經(jīng)驗、經(jīng)驗模型描述干燥過程物料中平均含水率的變化情況[9-13]。要進一步分析和體現(xiàn)物料內部水分的擴散過程,需確定物料內部水分的有效擴散系數(shù)[14-16],它反映物料在一定干燥條件下的脫水能力,也是干燥工藝優(yōu)化的重要指標之一。
本研究通過對數(shù)轉化與線性回歸,求得Page模型[17-20]的參數(shù)及相關系數(shù),并以Fick定律[20-21]為依據(jù),得到胡蘿卜沿軸向的水分擴散系數(shù)值。實驗目標:1)分析風速、溫度、切片厚度等參數(shù)對干燥特性的影響;2)利用回歸分析得到動力學模型;3)利用DPS軟件通過多元線性回歸分析,獲得胡蘿卜有效水分擴散系數(shù)模型。
1.1 材料與儀器
胡蘿卜,購于當?shù)爻?,要求新鮮無損傷。
MP4000B型電子秤(精度0.01 g) 上海第二天平廠;BK150型游標卡尺 上海申工量具有限公司;GIS-Ⅱ型熱風干燥機 黑龍江省農(nóng)業(yè)儀器設備制造廠;202型恒溫干燥箱 北京永光明醫(yī)療儀器廠;QDF-30M/S熱球式風速儀 北京市檢測儀器廠。
1.2 方法
1.2.1 實驗方案
實驗中將胡蘿卜洗凈去皮后切片,切片厚度為2~10 mm,側面用聚四氟乙烯密封帶完全密封,在物料盤上均勻平鋪一層,確保物料僅在上下表面進行傳質。在干燥過程中,由于物料厚度遠小于其直徑,所以物料可被假定為大平板。忽略溫度梯度導致的水分擴散,物料中的水分擴散可被假設為沿物料內部厚度方向的一維擴散過程。具體實驗方案如表1所示。干燥過程中,每隔20 min將物料迅速取出并稱質量,記錄不同時刻胡蘿卜的質量后將物料放回繼續(xù)干燥,直至物料含水率小于0.05,即認為已達到平衡含水率,停止干燥。計算不同時刻的水分比(MR),作圖并將MR線性轉換后用Origin軟件進行線性回歸,以確定方程中有關參數(shù)。
表1 胡蘿卜薄層干燥方案Table 1 Scheme of thin layer drying for carrots
1.2.2 熱風干燥
實驗所用干燥機為GIS-Ⅱ型熱風干燥機主要由風機、風量調節(jié)板、風筒、加熱器、控制器、干燥筒及操作臺等組成。在干燥過程中,氣流經(jīng)過風筒加熱后進入干燥筒,穿過干燥筒內的物料層,將物料加熱干燥。風量大小由調節(jié)板調節(jié)。
1.2.3 含水率的測定
物料初始含水率采用常壓干燥法計算。設定干燥箱溫度為105℃,將鮮胡蘿卜切成薄片,稱取一定量的樣品,放入干燥箱除濕干燥,待其質量不再發(fā)生變化后干燥結束,根據(jù)干燥前后樣品的質量變化,計算得到胡蘿卜初始含水率為8.73(以干基計)。本研究中物料含水率均以干基表示。按照式(1)計算干基含水率。
式中:M為物料干基含水率;m為物料質量/kg;md為絕干物料的質量/kg。
1.2.4 動力學模型
目前干燥研究大多認為擴散是水分遷移的主要途徑。物料內部水分擴散是一個復雜的傳質過程,涉及分子擴散、毛細流動、水力學流動和表面擴散等現(xiàn)象[21]。在大多數(shù)情形下,F(xiàn)ick第二定律被廣泛用于描述水分擴散過程[18,22],同時適用于本實驗的計算。本實驗的假設條件如下:1)物料內部的水分運動是在水分梯度作用下的液態(tài)擴散;2)忽略干燥過程中物料內部溫度梯度對水分擴散的影響;3)忽略干燥過程中胡蘿卜的收縮變形;4)干燥開始時,物料內部水分和溫度均勻分布[9-10]。
Fick第二定律表達式為:
式中:Deff為有效水分擴散系數(shù)/(m2/s);x為厚度/m;t為干燥時間/s。
Fick定律的數(shù)值解為:
式中:MR為水分比;H為物料厚度/m。
其簡化形式為:
該數(shù)值解可進一步簡化為:
式中:A為常數(shù);z為干燥常數(shù)。
此式即為常用的Henderson和Pabis經(jīng)驗模型。
Page模型也廣泛應用于薄層干燥模型[23-26]。其模型方程為:
式中:k為干燥常數(shù);m為常數(shù)。
由于平衡水分含量Me相對于M和M0很小,可忽略不計,所以(M-Me)/(M0-Me)可近似認為M/M0。
1.2.5 有效水分擴散系數(shù)計算
有效水分擴散系數(shù)是表示物料中水分擴散情況的重要指標,反映物料在一定干燥條件下的脫水能力。作出-lnMR和t的關系圖,并進行線性回歸,其斜率為F,即可得到干燥過程的有效水分擴散系數(shù)的計算式為:
2.1 工藝參數(shù)對干燥特性的影響
圖1 不同切片厚度下的MR--t關系曲線Fig.1 Relationship betweebn MR and t at different sample thickness
由圖1可知,在溫度60℃、風速0.6 m/s條件下,隨著胡蘿卜切片厚度的減少,胡蘿卜水分下降速率加快,說明胡蘿卜片越薄,單位體積的表面積增大,與熱空氣的對流換熱面積增大,傳熱速率增加。此外,物料厚度降低,會導致內部的傳遞阻力減小,水分擴散速率增大,最終提高干燥速率。胡蘿卜片厚度為10 mm時,干燥時間顯著增加,而在切片厚度為2 mm時,會發(fā)生嚴重變形,品質嚴重下降。因此,物料切片厚度過厚或過薄均不利于干燥的最終結果,在研究中宜選用4~8 mm的胡蘿卜片進行干燥。
圖2 不同干燥溫度下的MR--t關系曲線Fig.2 Relationship between MR and t at different drying temperatures
由圖2可知,對于4 mm的胡蘿卜切片,在0.6 m/s的風速條件下,40 ℃需要干燥11.5 h,而70 ℃僅需4.5 h,說明干燥溫度越高,傳熱速率越快,水分下降越快,水分擴散速率越大。但溫度過高,胡蘿卜干燥品質明顯變差,而溫度過低,干燥時間過長,降低干燥效率,所以干燥過程中采用50~70 ℃的干燥溫度較好。從曲線的斜率可以看出,中前期干燥速率快,后期逐漸減少,這是由于隨著干燥的進行,物料的水分不斷減少,水分梯度的作用越來越小,干燥界面內移,干燥過程由表面?zhèn)髻|控制轉化為內部擴散控制,干燥曲線的斜率越來越平坦。
由圖3可知,在溫度60 ℃、切片厚度4 mm條件下,隨著干燥風速的增大,胡蘿卜片水分比下降速率增大,說明風速越大,傳熱傳質邊界層越薄,物料與干燥介質的對流傳質越快,氣流帶走的水分越多,有效水分擴散系數(shù)也越大。但當表面水分的汽化量與內部滲出的水分量達到平衡后,干燥速率不再隨著風速的增大而加快,所以增加干燥介質的風速,對物料的前期干燥速率影響極顯著,但后期影響不顯著。
圖3 不同風速下的MR--t關系曲線Fig.3 Relationship between MR and t at different air flow rates
2.2 動力學模型比較
由圖1~3所示的MR-t圖可以看出,其變化呈曲線關系,將式(5)、(6)進行線性轉換,可得到以下方程:
將實驗數(shù)據(jù)進行分析和處理,分別作出lnMR-t和ln(-lnMR)-lnt的關系曲線圖如圖4~6所示。
圖4 不同切片厚度下lnMR--t和ln(-lnMR)--lnt關系曲線Fig.4 lnMR vs. t and ln(-lnMR)vs. lnt at different sample thickness
圖5 不同干燥溫度下lnMR-t和lnn(-lnMR)--lnt關系曲線Fig.5 lnMR vs. t and ln(-lnMR)vs. lnt at different drying temperatures
圖6 不同風速下lnMR-t和ln(-lnMR)--lnt關系曲線Fig.6 lnMR vs. t and ln(-lnMR)vs. lnt at different air flow rates
表2 回歸方程及相關系數(shù)Table 2 regression equations and correlation coefficients
將圖4~6所示的ln(-lnMR)-lnt和-lnMR-t關系曲線進行線性回歸,得到回歸方程及相關系數(shù)如表2所示。胡蘿卜薄層干燥過程中,實驗數(shù)據(jù)在lnMR-t坐標系內相關系數(shù)在0.984 7~0.999 1之間變化,在ln(-lnMR)-lnt坐標系內相關系數(shù)在0.996 8~0.999 6之間變化。由相關系數(shù)可知,實驗數(shù)據(jù)在lnMR-t和ln(-lnMR)-lnt坐標系內均具有較好的線性關系,但對比分析可得,實驗數(shù)據(jù)在ln(-lnMR)-lnt坐標系內線性擬合精度更高,對ln(-lnMR)-lnt進行指數(shù)轉換即得Page方程,說明Page方程更適合描述本實驗胡蘿卜薄層干燥過程。為了更精確的選擇干燥模型,對幾種不同實驗條件下得到的ln(-lnMR)-lnt實驗數(shù)據(jù)用Origin軟件進行線性回歸,得到模型參數(shù)k、m及決定系數(shù)R2值,結果如表3所示。
表3 不同參數(shù)下的模型參數(shù)值及有效水分擴散系數(shù)值Table 3 Model parameters and Deeffffat different drying parameters
由表3不同參數(shù)下的模型參數(shù)值及有效水分擴散系數(shù)值可知,決定系數(shù)R2值在0.982 8~0.999 0之間變化,說明線性回歸得到的結果較好,即Page模型可以很好地描述胡蘿卜薄層干燥過程。由表3可知,模型參數(shù)k和m的值均隨實驗條件的改變而變化,即在不同的干燥溫度(T)、風速(v)及厚度(H)條件下,它們的值也隨之發(fā)生變化,因此,干燥常數(shù)k和m是T、v及H的函數(shù)。根據(jù)實驗數(shù)據(jù),用DPS軟件對k和m進行回歸分析,得到參數(shù)的估計值及回歸方程的顯著性分析如表4、5所示。
表4 Page方程參數(shù)k和m回歸參數(shù)估計Table 4 Regression parameter estimates (k and m) of Page equation
由參數(shù)估計量可以獲得k和m與T、v和H的關系模型:
表5 方差分析表Table 5 Analysis of variance
由表5可知,Page模型參數(shù)k和m均有F>F0.01,方程回歸效果顯著。為了驗證模型的精確性,將50 ℃、0.3 m/s、5 mm及60 ℃、0.5 m/s、4 mm條件下的實驗數(shù)據(jù)與模擬值進行比較,其結果如圖7所示。在整個干燥過程中,Page模型預測值與實驗值擬合較好,任意時刻MR的預測值與實驗值相對偏差(相對偏差=|實驗值-預測值|/實驗值)均小于6%,說明胡蘿卜薄層干燥特性曲線為典型的指數(shù)曲線,Page模型能較準確地描述胡蘿卜在薄層干燥過程中的水分變化規(guī)律。
圖7 不同條件下MR計算值與實驗值的比較Fig.7 Comparison between the experimental and calculated MR under different conditions
2.3 有效水分擴散系數(shù)模型
研究有效水分擴散系數(shù)對深入分析物料內部水分擴散過程及優(yōu)化干燥工藝具有重要意義,主要包括液相擴散、水蒸氣擴散以及其他可能的質量傳遞現(xiàn)象,與物料成分、結構、溫度和含水率等因素有關。實驗研究了Deff和溫度T、風速v和物料厚度H的關系,結果如圖8~10所示。Deff和溫度T、風速v和物料厚度H呈較好的線性關系。用Origin軟件分別對曲線進行線性回歸,得到不同切片厚度、溫度及風速條件下的回歸方程及相關檢驗值如式(12)~(14)所示。
通過式(12)~(14)可知,有效水分擴散系數(shù)與實驗參數(shù)具有較好的線性關系,不同實驗參數(shù)對有效水分擴散系數(shù)的影響程度不同,其中切片厚度對有效水分擴散系數(shù)的影響最顯著。
圖8 不同厚度下的有效水分擴散系數(shù)Fig.8 Deffat different sample thickness
圖9 不同溫度下的有效水分擴散系數(shù)Fig.9 Deffat different drying temperatures
圖10 不同風速下的有效水分擴散系數(shù)Fig.10 Deffat different air flow rates
將Deff與T、v和H用DPS軟件進行多元線性回歸分析,得到Deff的回歸方程為:
由回歸方程可知,方程的相關系數(shù)較高,可用于簡化的傳質擴散系數(shù)方程。Deff隨著物料厚度、風速、溫度的增大而增大。一定溫度及風速下,物料越厚,水分梯度越大,水分擴散系數(shù)越快;一定物料厚度及溫度下,風速越大,物料表面的空氣流動速度越快,蒸發(fā)面上的空氣層越薄,即傳熱傳質邊界層越薄,物料與干燥介質的對流傳質加強,有效水分擴散系數(shù)越大;一定物料厚度及風速下,溫度越高,水分子運動越劇烈,水分遷移越快,有效水分擴散系數(shù)越大。實驗計算得到的胡蘿卜有效水分擴散系數(shù)在0.84×10-9~6.69×10-9m2/s之間變化,與大多數(shù)文獻中得到的胡蘿卜有效水分擴散系數(shù)結果類似[13,17-19]。
3.1 物料的干燥速率與干燥介質的溫度、速度及物料的切片厚度密切相關。提高干燥介質溫度、減小物料厚度、增大干燥介質風速,均有利于加快物料中的水分擴散,提高干燥速率。在實際生產(chǎn)干燥過程中,應根據(jù)實際情況采取合適的措施來提高干燥速率,強化有效水分擴散系數(shù)。
3.2 胡蘿卜薄層干燥曲線為典型的指數(shù)曲線。在實驗條件范圍內,Page模型預測值與實驗值擬合精度較高,能準確預測胡蘿卜薄層干燥過程中的水分比變化規(guī)律。對不同溫度、風速及物料厚度條件下的水分比數(shù)據(jù)進行回歸分析,建立了描述干燥規(guī)律的Page數(shù)學模型如下:MR=exp(-ktm),其中,m= 0.525 8+0.158 9v+0.007 6T+0.007 2H;k=0.019 7+0.002v+0.000 026T-0.002 15H。
3.3 根據(jù)物料的干燥特性曲線,確定其有效水分擴散系數(shù)如式(15)所示,實驗得到的胡蘿卜有效水分擴散系數(shù)在0.84×10-9~6.69×10-9m2/s之間隨著干燥溫度、風速和物料厚度的升高而增大。
[1] SUMNU G, TURABI E, OZTOP M. Drying of carrots in microwave and halogen lamp-microwave combination ovens[J]. LWT-Food Science and Technology, 2005, 38(5): 549-553.
[2] LI Ruijie, ZHANG Min. Drying methods affect the quality and hygroscopic capacity of carrot chips[J]. 食品與生物技術學報, 2010, 29(3): 342-349.
[3] 郭寬, 胡云峰, 李喜宏, 等. 不同品種胡蘿卜切分加工物理性能研究[J]. 食品科學, 2004, 25(11): 184-187.
[4] 李偉榮, 任愛清, 陳國寶. 胡蘿卜干燥技術研究進展[J]. 現(xiàn)代農(nóng)業(yè)科技, 2010(13): 388-390.
[5] 徐小東, 崔政偉. 農(nóng)產(chǎn)品和食品干燥技術及設備的現(xiàn)狀和發(fā)展[J].農(nóng)業(yè)機械學報, 2005, 36(12): 171-174.
[6] 李瑞杰, 張慜. 不同干燥方式對胡蘿卜片吸濕性及品質的影響[J].食品與生物技術學報, 2010, 29(3): 342-349.
[7] MASKAN M. Drying, shrinkage and rehydration characteristics of kiwifruits during hot air and microwave drying[J]. Journal of Food Engineering, 2001, 48(2): 177-182.
[8] O’NEILL M B, RAHMAN M S, PERERA C O, et al. Colour and density of apple cubes dried in air and modified atmosphere[J]. International Journal of Food Properties, 1998, 1(3): 197-205.
[9] 肖美添, 朱艷, 劉華信. 薄層介質干燥過程內部水分擴散[J]. 華僑大學學報: 自然科學版, 2003, 24(2): 184-188.
[10] 劉榴, 路倩倩, 劉相東. 基于干燥過程的生物材料有效水分擴散系數(shù)的模型與測定[C]//第十三屆全國干燥會議論文集. 葫蘆島: 中國化工學會, 2011: 215-221.
[11] SABAREZ H T, PRICE W E. A diffusion model for prune dehydration[J]. Journal of Food Engineering, 1999, 42(3): 167-172.
[12] AZZOUZ S, GUIZANI A, JOMAA W, et al. Moisture diffusivity and drying kinetic equation of convective drying of grapes[J]. Journal Food Engineering, 2002, 55(4): 323-330.
[13] DOYMAZ I. Convective air drying characteristics of thin layer carrots[J]. Journal of Food Engineering, 2004, 61(3): 359-364.
[14] TONG C H, LUND D B. Effective moisture diffusivity in porous materials as a function of temperature and moisture content[J]. Biotechnology Progress, 1990, 6(1): 67-75.
[15] 曾令彬, 趙思明, 熊善柏, 等. 風干白鰱的熱風干燥模型及內部水分擴散特性[J]. 農(nóng)業(yè)工程學報, 2008, 24(7): 280-283.
[16] 關志強, 王秀芝, 李敏, 等. 荔枝果肉熱風干燥薄層模型[J]. 農(nóng)業(yè)機械學報, 2012, 43(2): 151-158.
[17] RUIZ-CABRERA M A, SALGADO-CERVANTES M A, WALISZEWSKI-KUBIAK K N, et al. Effect of path diffusion on the effective moisture diffusivity in carrot slabs[J]. Drying Technology, 1997, 15(1): 169-181.
[18] ZIELINSKA M, MARKOWSKI M. Air drying characteristics and moisture diffusivity of carrots[J]. Chemical Engineering and Processing: Process Intensification, 2010, 49(2): 212-218.
[19] MULET A. Drying modeling and water diffusivity in carrots and potatoes[J]. Journal of Food Engineering, 1994, 22(1): 329-348.
[20] 吉永奇. 金銀花干燥特性及工藝參數(shù)研究[D]. 洛陽: 河南科技大學, 2008.
[21] 李陽春, 陳光明, 王劍鋒, 等. 含濕多孔生物材料干燥過程的傳熱傳質分析[J]. 農(nóng)業(yè)機械學報, 2003, 34(4): 82-85.
[22] 楊彬彬. 多孔介質干燥分形孔道網(wǎng)絡模擬及實驗研究[D]. 北京: 中國農(nóng)業(yè)大學, 2006.
[23] 趙義. 熱及濕敏性物料吸附式除濕干燥特性研究[D]. 廣州: 華南理工大學, 2005.
[24] MADAMBA P S, DRISCOLL R H, BUCKLE K A. The thin-layer drying characteristics of garlic slices[J]. Journal of Food Engineering, 1996, 29(1): 75-97.
[25] DIAMANTE L M, MUNRO P A. Mathematical modelling of the thin layer solar drying of sweet potato slices[J]. Solar Energy, 1993, 51(4): 271-276.
[26] DOYMAZ I, PALA M. The effects of dipping pretreatments on airdrying rates of the seedless grapes[J]. Journal of Food Engineering, 2002, 52(4): 413-417.
Dynamic Modeling of Thin Layer Drying of Carrots
CHONG Cui-juan, ZHU Wen-xue*, LIU Yun-hong, WU Jian-ye, FAN Jin-ling, LUO Lei
(College of Food and Bioengineering, Henan University of Science and Technology, Luoyang 471023, China)
The thin layer drying characteristics of carrot slices were investigated and moisture content (relative to dry matter) was modeled as a function of drying time under different conditions of temperature, air flow rate and sample thickness. Meanwhile, based on Fick’s diffusion law, the effective diffusion coefficients of water (Deff) in carrots were calculated and modeled as well. Results showed that the thin layer drying kinetics of carrots could be described by the Page equation, the coefficients of equation were obtained through regression analysis, and the equations of Deffvs. air temperature, air flow rate and sample thickness were obtained by multiple linear regression. The experimentally observed Deffvalues increased from 0.84 × 10-9to 6.69 × 10-9m2/s with increasing drying temperature, air flow rate and sample thickness.
carrot; drying; dynamic modeling; moisture diffusion coefficient
TS205.1
A
1002-6630(2014)09-0024-06
10.7506/spkx1002-6630-201409006
2013-04-23
國家自然科學基金面上項目(31171723)
種翠娟(1986—),女,碩士研究生,研究方向為農(nóng)產(chǎn)品加工。E-mail:905497479@qq.com
*通信作者:朱文學(1967—),男,教授,博士,研究方向為農(nóng)產(chǎn)品加工。E-mail:zwx@mail.haust.edu.cn