物料-干燥箱一體化的香菇熱風干燥模擬方法及應用

胡如響，孫東亮,*，朱躍強，盧星宇，王 鵬，宇 波

（1.北京石油化工學院機械工程學院，北京 102617；2.西安交通大學能源與動力工程學院，陜西 西安 710049）

果蔬類農產品是人類生產生活所必需的食物，其新鮮采摘后的含水率大部分高達85%以上，高含水率給長期貯存和長距離運輸帶來了巨大困難，農產品干燥技術應運而生。干燥作為一種延長水果、蔬菜等高含水率食品貯存時間的方法，主要通過降低含水率以降低水的活性和抑制微生物的生長繁殖，從而達到延長貯存時間的目的[1]。目前最常用的干燥技術為熱風干燥，但干燥過程中存在產品品質降低、耗時等缺陷，且干燥過程是一個復雜的傳熱傳質過程，影響因素眾多，一直備受學者們的關注[2-6]。

在現有的農產品熱風干燥研究中，主要分為試驗研究和模擬研究。試驗研究中大部分聚焦于通過大量試驗探究提升產品品質和降低干燥能耗的方法[7-11]，存在周期長、耗材量大等缺陷。而在熱風干燥的模擬研究中，大部分的研究聚焦于物料層面，通過描述單個或幾個物料的熱風干燥過程，探究物料內部的傳熱傳質過程以及物料的變形、干燥條件對產品品質的影響等[12-15]。而在實際工程應用中干燥條件多變且分布不均勻，所以物料層面的研究無法考慮到干燥箱內流速、溫度和濕度不均勻性的影響，導致結果與實際情況不相符。基于此，部分學者對干燥箱層面的熱風干燥模擬展開了研究。Smolka[16]和Kj?r[17]等通過模型模擬研究了干燥箱內干燥參數的分布情況，并采用不同方法降低了干燥參數的分布不均勻性，但是他們采用的模型中沒有考慮物料對流場、溫度場和濕度場的影響，與實際情況偏差較大。這是因為在干燥過程中物料內的液態水吸熱蒸發成水蒸氣，會加劇溫度場和濕度場的分布不均勻性。為了解決這一問題，Wang Peng等[18]嘗試將物料和干燥箱進行耦合，并將人參堆積區域假定為多孔介質，考慮了人參對流場、溫度場和濕度場的影響。但在這項研究中人參干燥過程中的水分蒸發速率被設置為常數，沒有考慮干燥速率在時間和空間上的變化，與實際的干燥過程不符。

綜上，目前干燥箱規模的熱風干燥模擬研究重點主要集中在干燥箱方面，缺乏對物料本身干燥過程的分析，沒有真正將物料與干燥箱耦合起來。為解決這一問題，本研究構建了物料-干燥箱一體化的香菇熱風干燥模擬方法，并通過模擬結果分析了香菇干燥過程中的干燥不均勻性。

1 材料與設備

1.1 材料

新鮮香菇（Lentinus edodes）購自北京某生鮮超市，初始濕基含水率約為88.56%，選取質量和體積相似的香菇裝入密封袋后貯存在醫用冷藏箱中（（4±1）℃）。

1.2 儀器與設備

多功能干燥機（圖1） 常州一步干燥設備有限公司：干燥箱（1.42 m×1.20 m×1.65 m）、YX3-905-2循環風機（額定功率為1.5 kW） 江蘇運中電機股份有限公司；LDR0.016-0.7蒸汽加熱器（額定功率為12 kW）張家港威孚熱能股份有限公司；XKY-1H熱泵（額定功率為1.05 kW） 中山愛美泰電器有限公司；YC-100醫用冷藏箱 澳柯瑪股份有限公司；BSA220.4電子分析天平上海卓精電子科技有限公司；XY-100MW鹵素水分分析儀常州市幸運電子設備有限公司。

2 物料-干燥箱一體化的香菇熱風干燥模擬方法

2.1 定工況條件下香菇干燥動力學模型

2.1.1 定工況條件下香菇熱風干燥試驗

待試驗前取出香菇晾至室溫并去掉菇柄，保留約1 cm蒂根，單層平鋪于干燥架隔板上，每一層平板平鋪9 列13 行，總共12 層（圖2）。在實際的工業生產過程中，香菇熱風干燥的溫度多為50～60 ℃，風速多為2～4 m/s，故本試驗的空氣溫度設置5 個水平，分別為45、50、55、60 ℃和65 ℃，風速設置5 個水平，分別為1、2、3、4 m/s和5 m/s，相對濕度設置了5 個水平，分別為0.25、0.30、0.35、0.40和0.45，進行3因素5水平的正交試驗，如表1所示。試驗過程中，選取第1層上第1列的香菇，利用電子分析天平每0.5 h測定其質量變化，重復測定3 次取平均值。干燥至GB 7096—2014《食品安全國家標準 食用菌及其制品》規定的香菇安全貯存標準以下，即濕基含水率不高于13%[19]，隨后將樣品密封、常溫保存。

表1 熱風干燥試驗工況Table 1 Experimental conditions of hot air drying

圖2 香菇在干燥箱中的排列方式Fig.2 Arrangement of Lentinus edodes in the drying chamber

2.1.2 試驗參數的獲取

干基含水率Mt的計算如公式（1）所示。

式中：mt為干燥至t時刻香菇的質量/kg；md為香菇的干基質量/kg；Mt為干燥到t時刻香菇的干基含水率/（kg/kg）。

水分比（moisture rate，MR）的計算如公式（2）所示。

式中：Me為干燥結束時香菇的干基含水率/（kg/kg）；M0為初始時刻香菇的干基含水率/（kg/kg）。

農產品干燥過程中水分比隨時間的變化通常用干燥動力學模型描述，目前常見的干燥動力學模型主要有Newton、Page等半理論模型，如表2所示。

表2 常見的干燥動力學模型Table 2 Common drying kinetic models

根據表2中常見的干燥動力學模型，可以將其寫成統一的形式，如式（3）所示。

式中：t為干燥時間/h；A、B、C均為擬合參數。本研究將采用式（3）擬合試驗數據。

為了評價數據擬合后的效果，通常使用決定系數R2、離差平方和χ2、均方根誤差（root mean square error，RMSE）和平均相對誤差e來說明擬合結果的準確性，R2的值越接近1，χ2、RMSE和e的值越接近于0，說明擬合結果的準確度越高，效果越好[26]。上述4 個評價指標分別采用式（4）～（7）進行計算[27]。

式中：MRexp為水分比的試驗值；MRpre為水分比的擬合預測值；n為試驗數據的總數；m為模型中參數的個數。

2.1.3 定工況條件下干燥動力學模型

通過上述試驗已經得到了不同工況下香菇的質量隨干燥時間的變化情況，利用Origin軟件對試驗結果進行處理，采用式（3）進行干燥動力學模型擬合，可以得到不同試驗工況下香菇的干燥動力學模型和對應模型的統計分析結果，如表3所示。

表3 不同工況下的模型擬合方程及統計分析結果Table 3 Fitting model equations under different conditions and their figures of merit

通過表3中的統計分析結果可以看出，R2的值均大于0.99，且χ2和RMSE的值均接近于0，以干燥試驗數據為對照，對回歸模型進行驗證，試驗數據和擬合值的平均相對誤差e小于1%，說明常見干燥動力學模型的通式（3）對于香菇在定工況條件下的熱風干燥試驗數據具有較高的預測精度。

分別建立式（3）中擬合參數A、B、C關于空氣溫度、空氣相對濕度和風速的回歸關系，如式（8）～（10）所示。

式中：T為空氣溫度/℃；RH為空氣的相對濕度；U為風速/（m/s）。

以空氣溫度50 ℃、相對濕度0.35和風速3 m/s的干燥試驗數據為對照，對模型通式及各擬合參數的回歸關系進行驗證，驗證結果的R2為0.993 7，表明模型的預測效果較好。說明模型通式（3）及其參數的多元回歸模型描述香菇在定工況條件下熱風干燥的干燥動力學是可行的。

2.2 變工況條件下香菇干燥動力學導數模型

上述香菇的干燥動力學模型雖然能夠預測在定工況條件下香菇干燥過程的水分比變化，但是在實際工程應用中，由于物料較多，吸熱蒸發產生的水蒸氣會影響干燥箱內的溫度場和濕度場，致使干燥箱中的工況處于動態變化中，所以上述定工況下的香菇干燥動力學模型不再適用，而且難以應用于干燥過程的計算流體力學（computational fluid dynamics，CFD）數值模擬之中。為了更好地解決上述問題，本試驗推導了適用于變工況下的香菇干燥動力學導數模型，并基于該導數模型得到了香菇內部液態水的蒸發速率。

以下基于定工況條件下的干燥動力學模型（3），給出變工況條件下香菇干燥動力學導數模型的推導過程。首先將通式（3）對時間求導，得到求導后的方程，如式（11）所示。

再將通式（3）代入到式（11）中，得到式（12）。

最后將擬合參數與干燥工況參數之間的多元回歸關系式（8）～（10）代入式（12）中，即得到能夠適用于變工況參數條件下的香菇熱風干燥動力學導數模型。由式（12）可知，該導數模型中不含有時間項，只含有水分比和工況參數，所以能夠保證在實際求解應用過程中水分比的連續性，可以用來預測變工況條件下香菇熱風干燥過程中的水分比變化。

在上述變工況條件下香菇熱風干燥動力學導數模型的基礎上，本研究提出了一種能夠適用于變工況條件下香菇內部液態水蒸發速率的計算方法。首先，單位體積內香菇液態水的蒸發速率可以用式（13）進行表示[28]。

然后將變工況條件下的導數模型（12）帶入到式（13），最終可以得到香菇內部液態水的蒸發速率，如式（14）所示。

由于擬合參數A、B、C是空氣溫度、空氣相對濕度和風速的多元回歸關系，因此能夠將香菇內的水分蒸發速率與干燥箱內部的工況變化過程聯系在一起，充分考慮時間和空間上的變化，進而能夠用于準確計算在變工況條件下香菇熱風干燥過程中的水分蒸發速率。

2.3 物料-干燥箱一體化的模型構建

2.3.1 物理模型及網格劃分

針對該試驗所用的干燥箱尺寸和香菇的排列方式（圖2），對物料和干燥箱建立一體化的物理模型，如圖3所示。主要將干燥箱內部劃分為兩個區域，分別為空氣自由流動區域和香菇平鋪區域。

圖3 物料-干燥箱一體化的物理模型Fig.3 Physical model of material-drying oven integration

由于計算區域相對規整，所以選用結構化網格劃分計算區域。另外考慮到香菇平鋪區域同時存在流動、傳熱傳質和香菇失水形變等復雜的物理現象，所以對此區域進行了加密處理，網格劃分結果如圖4所示。

圖4 物料-干燥箱網格劃分Fig.4 Mesh generation of material-drying oven integration

2.3.2 數學模型

在熱風干燥過程中，熱空氣從香菇表面流過，通過對流換熱作用使香菇的溫度上升，內部的水分不斷蒸發，產生的水蒸氣在壓差和擴散的作用下轉移到香菇的表面并被熱空氣不斷帶走。所以干燥過程是一個流動、傳熱和傳質的耦合過程。

控制方程基于如下假設建立：在香菇平鋪區域，香菇之間存在充滿空氣的孔隙，所以將此區域假設為多孔介質，由熱空氣和香菇組成，其中香菇被假設為固體；熱空氣為不可壓縮理想氣體；香菇平鋪區域保持局部非熱平衡，即熱空氣與香菇的溫度不同；由于空氣流速較高，以強迫對流為主，故不考慮重力的影響。

2.3.2.1 流動方程

在干燥過程中，干燥箱入口風速在1 m/s及以上，雷諾數大于等于5.6×104，所以干燥箱內部為湍流流動[29]。連續性方程和動量方程分別如式（15）和式（16）所示。

式中：u為矢量風速/（m/s）；ρg為熱空氣的密度/（kg/m3）；p為干燥箱內的相對壓強/Pa；νg為熱空氣的運動黏度/（m2/s）；νT為熱空氣的湍流黏度/（m2/s）；F為動量源項。

在自由流動區域，熱空氣沒有額外的阻力，所以該區域的動量源項F為零。在香菇平鋪區域，由于被設定為多孔介質，所以熱空氣在該區域的流動受多孔介質阻力的影響，根據Darcy-Forchheimer定律[30]，多孔介質內的流動阻力分為與流速成正比的黏性阻力和與流速的平方成正比的慣性阻力兩部分，采用Ergun公式（式（17））計算動量源項F。

式中：φ為孔隙率；dep為當量直徑/m。

2.3.2.2 能量方程

在自由流動區域，只有熱空氣存在，根據能量守恒定律，該區域的能量方程如式（18）所示。

式中：Cp,g為熱空氣的比熱容/（J/（kg·℃））；λg為熱空氣的導熱系數/（W/（m·℃））。

在香菇平鋪區域，存在氣相的空氣和固相的香菇，由于兩者的溫度相差較大，所以采用局部非熱平衡模型描述該區域的傳熱過程。根據局部非熱平衡假設[31]，多孔介質內各相溫度不相等，所以要為熱空氣和香菇分別建立能量方程，建立的能量方程分別如式（19）和式（20）所示。

式中：hAlg（T1－Tg）表示多孔介質內熱空氣與物料的對流換熱量/J；h為空氣與物料之間的對流換熱系數；Alg為界面面積密度/m－1；T1為香菇的溫度/℃；Tg為熱空氣的溫度/℃；ρ1為香菇的初始密度/（kg/m3）；Cp,1為香菇的比熱容/（J/（kg·℃））；λ1為香菇的導熱系數/（W/（m·℃））；表示香菇內部液態水蒸發消耗的蒸發潛熱/（J/kg）。

2.3.2.3 傳質方程

熱風干燥過程中，自由流動區域和香菇平鋪區域都存在水蒸氣的傳質過程，這主要是由于壓差形成的對流和濃度擴散造成的[32]。因為香菇內部的水分蒸發后會被周圍的熱空氣帶走，所以只有香菇平鋪區域的傳質方程存在質量源項。自由流動區域和香菇平鋪區域的傳質方程分別如式（21）和式（22）所示。

式中：ω為空氣中水蒸汽的質量分數/%；Dva為擴散系數/（m2/s）；為質量源項，即香菇內部液態水的蒸發速率/（kg/（m3·s））。

本研究中，每次干燥箱中香菇的數量都達到上千個，而在實際的工程應用中，干燥的香菇數量會更多，因此在這種情況下對每個香菇進行建模計算其內部的水分蒸發速率是不現實的，所以這里采用前文推導得到的式（14）對香菇內部液態水的蒸發速率進行計算。由式（14）可知，香菇內部液態水的蒸發速率的計算與擬合參數相關，而擬合參數又是空氣溫度、空氣相對濕度和風速的回歸關系，因此實現了香菇內部液態水的蒸發速率與傳熱傳質方程的耦合，充分考慮了香菇吸熱蒸發過程和干燥箱內工況環境變化的相互影響，實現了物料與干燥箱的真正耦合，最終構建了符合實際干燥過程的物料-干燥箱一體化香菇熱風干燥模型。

2.3.3 邊界條件

入口邊界和出口邊界的位置如圖3所示，其中入口邊界的空氣溫度、相對濕度和風速為固定值，湍動能和湍流耗散率分別由式（23）和（24）確定，出口邊界為自由流動出口。干燥箱其余壁面具有良好的保溫特性，為絕熱壁面，速度和水蒸氣質量分數分別為無滑移和零法向梯度條件；壁面湍流采用壁面函數，湍動能的法向梯度為零，湍流耗散率由式（25）確定。

式中：ωwall為壁面的湍流耗散率；ωVis為黏性底層（壁面的法向距離的無量綱量y＋＜5）的湍流耗散率；ωlog為對數律層（y＋＞30）的湍流耗散率。

2.3.4 模型參數

2.3.4.1 孔隙率

香菇在干燥過程中會發生收縮，且收縮率與水分比之間呈如式（26）的函數關系[33]，這使得香菇的體積會隨著干燥過程不斷減小，對應的孔隙率會不斷增大，所以孔隙率由式（27）計算。

式中：SR為香菇的收縮率；V1為香菇的體積/m3；Vg為香菇平鋪區域空氣的體積/m3；V1,0為香菇初始時刻的體積/m3。

2.3.4.2 等效直徑

對于非球形粒子填充的多孔介質，用式（28）、（29）求解等效直徑可使多孔介質阻力的計算誤差最低[33]。

式中：φ為球形度；A1為單個香菇物料的表面積/m2；As1為與單個香菇體積相同的球的表面積/m2。

在干燥過程中，香菇會發生收縮，對應的等效直徑也會發生變化。香菇的菌蓋直徑約為0.07 m，高度約為0.02 m，形狀和球冠類似，所以本試驗采用球冠的相關公式計算香菇的體積和表面積。假設香菇各向同性，在各個方向上均勻收縮，則等效直徑可進一步表示為式（30）。

式中：dep,0為初始時刻香菇的等效直徑/m。

2.3.4.3 界面面積密度

界面面積密度是固體與流體交界面的面積和多孔介質總體積的比值，即本研究中香菇的表面積和香菇平鋪區域體積的比值。香菇在干燥過程中的收縮會導致表面積減小，對應的界面面積密度也會減小，變化后的界面面積密度由式（31）計算。

式中：Alg,0為初始時刻的界面面積密度/m－1。

2.3.4.4 香菇的熱物性參數

香菇主要由固體基質、水和氣體組成，隨著干燥過程的進行，各組分的體積分數和質量分數也會隨之變化，對應香菇的密度、導熱系數和比熱容分別由式（32）、（33）和（34）計算。

式中：ρs為香菇固體基質的密度/（kg/m3）；ρw為香菇中液態水的密度/（kg/m3）；ρg為香菇孔隙中空氣的密度/（kg/m3）；φs,0為初始時刻香菇中固體基質的體積分數/%；φw,0為初始時刻香菇中液態水的體積分數/%；λs為香菇固體基質的導熱系數/（W/（m·℃））；λw為香菇中液態水的導熱系數/（W/（m·℃））；λg為香菇孔隙中空氣的導熱系數/（W/（m·℃））；Cp,s為香菇固體基質的比熱容/（J/（kg·℃））；Cp,w為香菇中液態水的比熱容/（J/（kg·℃））；Cp,g為香菇孔隙中空氣的比熱容/（J/（kg·℃））。物性參數的具體值如表4所示[34]。

表4 物料-干燥箱一體化模型的物性參數Table 4 Physical parameters adopted in the model of material-drying oven integration

2.3.5 模型求解

本試驗采用商業軟件FLUENT對所建模型進行瞬態模擬求解。非穩態項、對流項分別采用一階隱式格式、二階迎風格式進行離散，采用SIMPLE算法進行求解。壓力和速度的亞松弛系數分別為0.3和0.7，湍動能和比耗散率的亞松弛系數為0.8，溫度的亞松弛系數為1.0。總模擬的干燥時間為14 h，時間步長設置為6 s。

3 結果與分析

3.1 模擬結果與試驗數據對比

首先，進行一組進口空氣溫度為55 ℃、相對濕度為0.35、風速為2 m/s條件下的驗證實驗，將數值模擬結果與干燥箱內第一層上不同位置香菇的水分比實際變化情況進行對比，結果如圖5所示。可以看出，試驗數據和模擬結吻合度較高，相對偏差很小。為了更加準確地描述二者之間的微小偏差，計算了模擬結果與試驗數據的R2，第1、5、9列的R2分別為0.995 6、0.998 3和0.998 6，3 個不同位置香菇水分比的R2都大于0.99，說明本研究所提出的模擬方法能夠用于準確模擬香菇的熱風干燥過程。

圖5 第一層上不同位置香菇的水分比模擬值與實際值的對比結果Fig.5 Comparison between simulated and actual values of moisture ratios at different locations on the first layer

3.2 模擬結果分析

本試驗對進口空氣溫度為55 ℃、相對濕度為0.35、風速2 m/s下的模擬結果展開研究，主要分析香菇干燥過程中的干燥不均勻性，包括干燥過程中不同層、不同列、不同行之間的水分比分布情況。考慮到物理模型在Z軸方向上的對稱性，所以只對干燥箱下半部分的第1、3、6層的模擬結果進行分析。

3.2.1 整體的干燥不均勻性

圖6給出了不同干燥時期香菇的水分比分布情況。根據水分比的大小分布，將干燥過程分為干燥前期、干燥中期、干燥后期和干燥末期4 個階段。從圖6可以看出，干燥前期的干燥不均勻性最為顯著，隨著干燥過程的進行，干燥不均勻性逐漸降低。

圖6 干燥過程中香菇水分比分布的變化Fig.6 Variation in moisture ratio distribution of Lentinus edodes during drying

在X軸方向即空氣流動方向上，干燥前期和中期香菇的干燥不均勻性比較顯著，水分比呈現逐漸增大的規律。這主要是因為隨著對流換熱過程的進行，空氣的熱量不斷向香菇傳遞，使得空氣的溫度下降，同時香菇干燥后產生的水蒸氣使得空氣的相對濕度逐漸升高，所以在空氣流動方向上干燥速率逐漸降低，水分比分布呈現逐漸增大的現象。而在干燥后期和末期，香菇的干燥不均勻性逐漸減小，且水分比呈現先增大后減小的分布規律。水分比先增大的原因與干燥前、中期水分比逐漸增大的原因一致，而后逐漸減小是由于在干燥后期和末期，香菇的體積大幅縮小，香菇間隙變大，結合出口處自由流動區域的體積突然變大，引起尾部香菇間的流速增大，導致尾部香菇干燥速率變快。此外，通過比較不同層香菇的水分比分布可以發現，越靠近干燥箱的頂部或底部，空氣流動方向上的干燥不均勻性越大。

從圖6中還可以看出，在Y軸方向即橫向方向上，兩側的水分比明顯高于中間位置。這是由于在兩側的位置靠近干燥箱壁面，干燥架與兩側壁面之間空氣的自由流動區域體積較大，導致大部分熱空氣流過自由流動區域，使得流過兩側香菇的熱空氣流速較小，同時，緩慢的空氣流動使得這一區域的空氣溫度更低，相對濕度更高，所以兩側香菇的干燥速率降低，水分比高于中間位置。隨著干燥過程的進行，橫向方向上空氣的溫度和相對濕度的不均勻性逐漸減小，但流速的不均勻性依然存在，所以即使在干燥后期和末期，橫向方向上的干燥不均勻性依然較大。此外，通過比較不同層香菇的水分比分布可以發現，越靠近干燥箱中間層，橫向方向上的干燥不均勻性越大。

綜上，干燥過程中香菇在空氣流動方向上和橫向方向上都存在明顯的干燥不均勻性。在空氣流動方向上，越靠近干燥箱頂部或底部干燥不均勻性越明顯；在橫向方向上，越靠近干燥箱中間層干燥不均勻性越明顯。

3.2.2 空氣流動方向上的干燥不均勻性

為了更加準確地描述香菇在空氣流動方向上的干燥不均勻性，選取第1、3、6層上的第1、3、5、7、9列香菇，并計算出不同時刻每列13 個香菇的水分比平均值，得到變化結果如圖7所示。每一層上不同列的香菇具有明顯的干燥不均勻性，其中第1層上不同列的香菇水分比差異最大，最大水分比偏差為0.16。另外，根據GB 7096—2014[19]的規定，香菇的水分比為0.026時達到貯存要求，視為干燥過程完成。根據此標準，由圖7可以得到第1、3、6層上不同列香菇達到干燥要求的最長時間差值分別為1.42、0.89、0.57 h。

圖7 干燥過程中空氣流動方向上不同位置香菇水分比的變化情況Fig.7 Variation in moisture ratios of Lentinus edodes in the direction of air flow at different positions during drying

圖8、9、10分別為不同層上不同列風速、空氣相對濕度和空氣溫度隨干燥過程的變化情況。可以看出，靠近干燥箱出口的第9列區域空氣相對濕度較高（圖9）、溫度較低（圖10），但由于風速較大（圖8），第9列香菇的干燥速率并沒有明顯降低，和第7列相差不多（圖7）。這一現象在第3層和第6層更加明顯，所以空氣流動方向上這兩層上香菇的干燥結果比第1層更為均勻。結合圖7可以看出，在干燥末期階段，最晚達到干燥要求的是第3列，而不是更加靠近出口的第7或第9列，這主要是因為干燥末期空氣的相對濕度和溫度趨于一致，而第3列的風速最小，導致其干燥速率降低，最晚達到干燥要求。

圖8 干燥過程中空氣流動方向上不同位置風速的變化情況Fig.8 Variation in air flow rate in the direction of air flow at different positions during drying

圖9 干燥過程中空氣流動方向上不同位置空氣相對濕度的變化情況Fig.9 Variation in air relative humidity in the direction of air flow at different positions during drying

圖10 干燥過程中空氣流動方向上不同位置空氣溫度的變化情況Fig.10 Variation in air temperature in the direction of air flow at different positions during drying

3.2.3 橫向方向上的干燥不均勻性

在橫向方向上，同樣對第1、3、6層上的香菇進行分析，由于在橫向方向呈中間對稱，所以選取每一層上的第1、3、5、7行進行研究。首先計算出不同層上每一行9 個香菇的水分比取平均值，得到水分比變化結果，如圖11所示。在相同時刻，同一層上越靠近邊側香菇的水分比越高，其中第1行香菇的水分比明顯高于其余幾行，這與3.2.1節中的分析結果一致。對比不同層上的水分比變化可以發現，越靠近干燥箱中間層，不同行的水分比差異越大，最大水分比偏差出現在第6層上，差值為0.08。另外，由圖11可以得到，第1、3、6層上不同行香菇達到干燥要求的最大時間差值分別為1.16、1.24、1.35 h。

圖11 干燥過程中橫向方向上不同位置香菇水分比的變化情況Fig.11 Variation in moisture ratios in the horizontal direction at different positions during drying

圖12、13、14分別給出了不同層上不同行的風速、空氣相對濕度和空氣溫度隨干燥過程的變化情況。可以看出，在同一層上越靠近中間位置，空氣的流速和溫度越高，空氣的相對濕度越低，所以同一層上越靠近中間位置，香菇的水分比越低。從圖中還可以發現，越靠近干燥箱中間層，在橫向方向上流速偏差越大，最大偏差達到0.28 m/s，導致越靠近干燥箱中間層香菇的干燥不均勻性越大。

圖12 干燥過程中橫向方向上不同位置風速的變化情況Fig.12 Variation in air flow rate in the horizontal direction at different positions during drying

圖13 干燥過程中橫向方向上不同位置空氣相對濕度的變化情況Fig.13 Variation in air relative humidity in the horizontal direction at different positions during drying

圖14 干燥過程中橫向方向上不同位置空氣溫度的變化情況Fig.14 Variation in air temperature in the horizontal direction at different positions during drying

4 結 論

本實驗以香菇為研究對象，開發了物料與干燥箱一體化的熱風干燥數值模擬方法，并通過實驗驗證了模擬方法的精確度，能夠用于準確模擬大批量的香菇熱風干燥過程。通過分析該方法的模擬結果得出以下結論：1）干燥過程中，整體存在明顯的干燥不均勻性，在干燥前期干燥不均勻性最為顯著，隨著干燥的進行，不均勻程度逐漸降低；2）在空氣流動方向上，越靠近干燥箱頂部或底部干燥不均勻性越明顯，同一層上達到干燥要求的最大時間偏差為1.42 h；3）在橫向方向上，越靠近干燥箱中部干燥不均勻性越明顯，達到干燥要求時的最大干燥時間偏差出現在第一層，為1.16 h。本研究所開發的香菇熱風干燥模擬方法可用于干燥工藝優化、干燥箱結構設計和產品質量提升等方面，可為香菇熱風干燥技術的進一步發展提供支撐。