王澤林，劉 芳，耿文廣，高 玲，張大鵬，李子淳，張 瀟

(1.山東建筑大學(xué) 熱能工程學(xué)院，山東 濟(jì)南 250101；2.齊魯工業(yè)大學(xué)山東省科學(xué)院能源研究所，山東 濟(jì)南 250353；3.齊魯工業(yè)大學(xué) 能源與動(dòng)力工程學(xué)院， 山東 濟(jì)南 250353；4.中華全國(guó)供銷合作總社濟(jì)南果品研究院，山東 濟(jì)南 250220)

1 概述

由于降低了物料的水分活性和微生物活性，并最大限度地減少了儲(chǔ)存期間的物理和化學(xué)變化，脫水干燥成為延長(zhǎng)果蔬儲(chǔ)存期最常用的方法之一[1]。對(duì)流干燥是最常見(jiàn)的干燥方式，在干燥過(guò)程中，熱量由干燥介質(zhì)傳遞至物料，物料表面溫度率先升高，表面水分蒸發(fā)形成內(nèi)外濕度梯度，最終使內(nèi)部水分傳遞到物料表面蒸發(fā)[2-3]。

預(yù)測(cè)干燥過(guò)程中物料內(nèi)部的濕度和溫度分布，對(duì)了解干燥機(jī)理、提高干燥效率具有重要的意義。然而，由于測(cè)溫設(shè)備易破壞物料的質(zhì)構(gòu)，因此通過(guò)實(shí)驗(yàn)測(cè)量干燥過(guò)程中物料內(nèi)部溫度的瞬態(tài)分布十分困難。數(shù)值模擬方法基于對(duì)流干燥過(guò)程中發(fā)生的多物理場(chǎng)耦合數(shù)學(xué)模型，能準(zhǔn)確預(yù)測(cè)干燥過(guò)程中物料內(nèi)部溫度場(chǎng)、濕度場(chǎng)的變化[4-6]，對(duì)深入認(rèn)識(shí)物料的對(duì)流干燥過(guò)程機(jī)制、優(yōu)化對(duì)流干燥工藝具有重要的理論研究和實(shí)際指導(dǎo)意義。

目前，基于宏觀水平對(duì)果蔬干燥進(jìn)行研究的方法主要基于擴(kuò)散模型[7]，擴(kuò)散模型假設(shè)導(dǎo)熱為能量傳遞方式，擴(kuò)散為水分傳遞方式[8]，許多學(xué)者在此基礎(chǔ)上對(duì)不同果蔬(如哈密瓜、芒果、獼猴桃、洋蔥等)進(jìn)行了大量研究[9-12]。從驗(yàn)證結(jié)果可以發(fā)現(xiàn)，模擬結(jié)果與實(shí)驗(yàn)結(jié)果具有良好的匹配性，但并不能對(duì)多孔介質(zhì)內(nèi)流體流動(dòng)的驅(qū)動(dòng)力如毛細(xì)力等其他輸運(yùn)機(jī)理做出解釋。因此，不能將所有的水輸運(yùn)過(guò)程都?xì)w為擴(kuò)散[13]。

Datta[14]最早提出多相多孔介質(zhì)模型(MPMM)，即以果蔬的多孔介質(zhì)特性為基礎(chǔ)，基于連續(xù)性假設(shè)，將多孔介質(zhì)的所有變量和參數(shù)在一個(gè)具有代表性的基本體積(即REV尺度)上平均，同時(shí)考慮多孔介質(zhì)中的傳質(zhì)以及水分蒸發(fā)過(guò)程。多相多孔介質(zhì)模型以果蔬的多孔介質(zhì)特性為基礎(chǔ)，基于宏觀的REV尺度，同時(shí)考慮多孔介質(zhì)中壓力驅(qū)動(dòng)與毛細(xì)擴(kuò)散[15]導(dǎo)致的水分與水蒸氣輸運(yùn)過(guò)程及氣液相轉(zhuǎn)換。Purli[16]將多相多孔介質(zhì)模型應(yīng)用于食品干燥過(guò)程中，并推導(dǎo)出適用于吸濕多孔材料的低強(qiáng)度、低溫度對(duì)流干燥的傳熱和傳質(zhì)的多相多孔介質(zhì)模型。Chen等人[17]將多相多孔介質(zhì)模型加以改進(jìn)，結(jié)合果蔬的多孔介質(zhì)屬性，將其用于棗片干燥的模擬中，通過(guò)設(shè)置對(duì)流邊界條件，建立棗片內(nèi)部的傳質(zhì)、傳熱方程，并進(jìn)行數(shù)值模擬。

目前，將多相多孔介質(zhì)模型應(yīng)用于果蔬切片干燥研究時(shí)，多未考慮空氣流動(dòng)對(duì)干燥的影響，導(dǎo)致模擬結(jié)果產(chǎn)生一定偏差。本文以蘋(píng)果切片為研究對(duì)象，模擬研究對(duì)流干燥過(guò)程中，蘋(píng)果切片平均含水率的影響因素以及蘋(píng)果切片內(nèi)部溫度、含水率的變化。

2 物理模型

將圓柱體蘋(píng)果切片作為研究對(duì)象，將其放在流道的底面中心，三維物理模型見(jiàn)圖1。流道長(zhǎng)×寬×高為150 mm×130 mm×50 mm。蘋(píng)果切片底面半徑為20 mm，高為5 mm。

圖1 三維物理模型

3 數(shù)學(xué)模型

3.1 基本設(shè)定

對(duì)流干燥過(guò)程中，熱濕傳遞過(guò)程較為復(fù)雜，影響因素多。為提高模擬計(jì)算效率，進(jìn)行以下設(shè)定：蘋(píng)果切片內(nèi)部孔隙均勻，各向同性，等效孔隙率為0.92[18]。蘋(píng)果切片內(nèi)部固、氣、液三相的分布視為連續(xù)，內(nèi)部水分僅以水蒸氣的方式離開(kāi)切片。忽略干燥過(guò)程中切片因體積收縮對(duì)內(nèi)部熱質(zhì)傳遞產(chǎn)生的影響。所有氣體均視為理想氣體，滿足理想氣體狀態(tài)方程。

3.2 控制方程

① 流體區(qū)域

a.流動(dòng)方程

空氣流道內(nèi)干燥空氣流動(dòng)為層流，連續(xù)性方程、動(dòng)量方程分別為[19]：

式中ρ——干燥空氣密度，kg/m3

t——時(shí)間，s

div——散度

u——干燥空氣在x、y、z軸方向上的速度矢量，m/s

ui——干燥空氣在x、y、z軸方向上的速度，m/s

p——干燥空氣壓力，Pa

μ——干燥空氣動(dòng)力黏度，Pa·s

F——體積力項(xiàng)

b.能量方程

空氣流道內(nèi)干燥空氣的熱量傳遞符合傅里葉定律，非穩(wěn)態(tài)熱量傳遞方程為[20]：

式中cp,a——干燥空氣氣相比定壓熱容，J/(kg·K)

T——干燥空氣氣相溫度，K

λa——干燥空氣氣相熱導(dǎo)率，W/(m·K)

r——水汽化潛熱，J/kg

meva——水蒸發(fā)率，kg/(m3·s)

② 蘋(píng)果切片多孔介質(zhì)區(qū)域

a.傳質(zhì)方程

蘋(píng)果切片內(nèi)部及干燥空氣中的水分傳遞符合斐克第二擴(kuò)散定律，非穩(wěn)態(tài)質(zhì)量傳遞方程為[20]：

液態(tài)水：

水蒸氣：

式中cL——蘋(píng)果切片內(nèi)部液態(tài)水濃度，mol/m3

ρL——液態(tài)水密度，kg/m3

kL——水的固有滲透率，m2

kr,L——水的相對(duì)滲透率，m2

μL——液態(tài)水動(dòng)力黏度，Pa·s

Dc——水分?jǐn)U散系數(shù)，m2/s

cv——蘋(píng)果切片內(nèi)部水蒸氣濃度，mol/m3

Dva——水蒸氣在空氣中二元擴(kuò)散系數(shù)，m2/s

ux——干燥空氣在x軸上的速度，m/s

b.能量方程

蘋(píng)果切片內(nèi)部的熱量傳遞符合傅里葉定律，非穩(wěn)態(tài)熱量傳遞方程為[20]：

式中ρs——蘋(píng)果切片多孔介質(zhì)區(qū)域有效密度，kg/m3

cp——蘋(píng)果切片多孔介質(zhì)區(qū)域有效比定壓熱容，J/(kg·K)

λ——蘋(píng)果切片多孔介質(zhì)區(qū)域有效熱導(dǎo)率，W/(m·K)

3.3 求解條件設(shè)置

邊界條件、初始條件見(jiàn)表1。進(jìn)風(fēng)口為速度入口，出風(fēng)口為壓力出口。除進(jìn)出風(fēng)口外，流道其他面均為絕熱無(wú)滑移邊界。

表1 模型邊界條件和初始條件

3.4 網(wǎng)格劃分

采用非結(jié)構(gòu)化網(wǎng)格對(duì)模型進(jìn)行網(wǎng)格劃分[21]，由蘋(píng)果切片邊界向氣流通道邊界增長(zhǎng)因子為1.2，為確保模擬準(zhǔn)確性，對(duì)蘋(píng)果切片表面進(jìn)行網(wǎng)格加密。將蘋(píng)果切片表面溫度變化率作為評(píng)價(jià)目標(biāo)，對(duì)比網(wǎng)格數(shù)35×104、51×104、69×104、80×104后發(fā)現(xiàn)，當(dāng)網(wǎng)格數(shù)超過(guò)69×104時(shí)，蘋(píng)果切片表面溫度變化率小于0.93%。綜合考慮精度及運(yùn)算速率，模型網(wǎng)格數(shù)確定為69×104。

3.5 模擬軟件與相關(guān)設(shè)定

基于COMSOL Multiphysics軟件的有限元法進(jìn)行模擬，相關(guān)設(shè)定如下。

① 物理場(chǎng)接口：a.添加層流、流體傳熱接口，用于模擬流道中干燥空氣的流動(dòng)以及與蘋(píng)果切片多孔介質(zhì)間的傳熱。軟件自動(dòng)添加非等溫流動(dòng)多物理場(chǎng)耦合接口，用于計(jì)算流動(dòng)傳熱的多物理場(chǎng)耦合過(guò)程。b.分別添加液態(tài)水、水蒸氣的稀物質(zhì)傳遞接口，用于模擬液態(tài)水、水蒸氣在蘋(píng)果切片多孔介質(zhì)中的擴(kuò)散，以及水蒸氣在蘋(píng)果切片多孔介質(zhì)表面與干燥空氣間的傳遞。

② 求解設(shè)置：同時(shí)求解動(dòng)量方程、能量方程、傳質(zhì)方程。在求解器配置中，選擇終止技術(shù)為容差，將相對(duì)容差設(shè)置為1×10-4。將最大迭代次數(shù)設(shè)為100次。

3.6 模型實(shí)驗(yàn)驗(yàn)證

將新鮮蘋(píng)果洗凈、削皮、切樣，得到底面半徑為20 mm、高為5 mm的蘋(píng)果切片。干燥空氣進(jìn)風(fēng)溫度為60 ℃，流速為1 m/s，進(jìn)風(fēng)相對(duì)濕度為20%。

將蘋(píng)果切片平均含水率隨時(shí)間變化的實(shí)驗(yàn)結(jié)果與模擬結(jié)果進(jìn)行比較。結(jié)果表明，蘋(píng)果切片平均含水率的實(shí)驗(yàn)結(jié)果與模擬結(jié)果隨時(shí)間的變化趨勢(shì)一致，最大相對(duì)誤差為4%，說(shuō)明模型的準(zhǔn)確性比較理想。

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

4.1 對(duì)平均含水率的影響

① 進(jìn)風(fēng)溫度

果蔬的干燥溫度一般不宜過(guò)高，高溫易導(dǎo)致褐變反應(yīng)以及破壞物料內(nèi)部組織結(jié)構(gòu)[22]。干燥空氣流速1 m/s、進(jìn)風(fēng)相對(duì)濕度20%時(shí)，不同進(jìn)風(fēng)溫度蘋(píng)果切片平均含水率隨時(shí)間的變化見(jiàn)圖2。由圖2可知，干燥空氣流速、進(jìn)風(fēng)相對(duì)濕度一定時(shí)，不同進(jìn)風(fēng)溫度蘋(píng)果切片平均含水率隨時(shí)間的變化趨勢(shì)基本一致，均為先下降然后趨于穩(wěn)定。進(jìn)風(fēng)溫度為60 ℃時(shí)，蘋(píng)果切片平均含水率率先達(dá)到穩(wěn)定，進(jìn)風(fēng)溫度為40 ℃時(shí)最后達(dá)到穩(wěn)定。增大進(jìn)風(fēng)溫度可有效提高干燥速率。

② 干燥空氣流速

進(jìn)風(fēng)溫度60 ℃、進(jìn)風(fēng)相對(duì)濕度20%時(shí)，不同干燥空氣流速蘋(píng)果切片平均含水率隨時(shí)間的變化見(jiàn)圖3。由圖3可知，進(jìn)風(fēng)溫度、進(jìn)風(fēng)相對(duì)濕度一定時(shí)，不同干燥空氣流速蘋(píng)果切片平均含水率隨時(shí)間的變化趨勢(shì)基本一致，均為先下降然后趨于穩(wěn)定。干燥空氣流速為2.0 m/s時(shí)，蘋(píng)果切片平均含水率率先達(dá)到穩(wěn)定，干燥空氣流速為1.0 m/s時(shí)最后達(dá)到穩(wěn)定。增大干燥空氣流速，對(duì)促進(jìn)蘋(píng)果切片平均含水率快速達(dá)到穩(wěn)定的作用并不明顯。

圖2 干燥空氣流速1 m/s、進(jìn)風(fēng)相對(duì)濕度20%時(shí)不同進(jìn)風(fēng)溫度蘋(píng)果切片平均含水率隨時(shí)間的變化

圖3 進(jìn)風(fēng)溫度60 ℃、進(jìn)風(fēng)相對(duì)濕度20%時(shí)不同干燥空氣流速蘋(píng)果切片平均含水率隨時(shí)間的變化

③ 進(jìn)風(fēng)相對(duì)濕度

進(jìn)風(fēng)溫度60 ℃、干燥空氣流速1.0 m/s時(shí)，不同進(jìn)風(fēng)相對(duì)濕度蘋(píng)果切片平均含水率隨時(shí)間的變化見(jiàn)圖4。由圖4可知，進(jìn)風(fēng)溫度、干燥空氣流速一定時(shí)，不同進(jìn)風(fēng)相對(duì)濕度蘋(píng)果切片平均含水率隨時(shí)間的變化趨勢(shì)基本一致，均為先下降然后趨于穩(wěn)定。進(jìn)風(fēng)相對(duì)濕度為20%時(shí)，蘋(píng)果切片平均含水率率先達(dá)到穩(wěn)定，進(jìn)風(fēng)相對(duì)濕度為40%時(shí)最后達(dá)到穩(wěn)定。減小進(jìn)風(fēng)相對(duì)濕度可有效提高干燥速率。

圖4 進(jìn)風(fēng)溫度60 ℃、干燥空氣流速1.0 m/s時(shí)不同進(jìn)風(fēng)相對(duì)濕度蘋(píng)果切片平均含水率隨時(shí)間的變化

4.2 內(nèi)部物理場(chǎng)變化

由以上分析可知，進(jìn)風(fēng)溫度、進(jìn)風(fēng)相對(duì)濕度對(duì)干燥速率影響比較大，而干燥空氣流速影響比較小。因此，考慮減少風(fēng)機(jī)能耗，選取在送風(fēng)溫度60 ℃、干燥空氣流速1 m/s、進(jìn)風(fēng)相對(duì)濕度20%條件下，研究蘋(píng)果切片干燥過(guò)程中內(nèi)部物理場(chǎng)的變化。

沿干燥空氣流向，在蘋(píng)果切片50%高度的徑向軸線上均勻設(shè)置觀測(cè)點(diǎn)1～5，測(cè)點(diǎn)1靠近進(jìn)風(fēng)口，測(cè)點(diǎn)5靠近出風(fēng)口。蘋(píng)果切片觀測(cè)點(diǎn)布置見(jiàn)圖5。

圖5 蘋(píng)果切片觀測(cè)點(diǎn)布置1～5.觀測(cè)點(diǎn)

① 蘋(píng)果切片內(nèi)部溫度

觀測(cè)點(diǎn)1～5溫度隨時(shí)間的變化見(jiàn)圖6。由圖6可知，蘋(píng)果切片內(nèi)部溫度經(jīng)歷快速上升、緩慢上升、趨于平穩(wěn)的過(guò)程。由于干燥初期蘋(píng)果切片與干燥空氣溫差比較大，各觀測(cè)點(diǎn)溫度上升明顯。隨著干燥進(jìn)行，蘋(píng)果切片與干燥空氣溫差減小，各觀測(cè)點(diǎn)溫度上升速率變慢并逐漸趨于穩(wěn)定。迎風(fēng)面及邊緣溫度高于中間位置。

圖6 觀測(cè)點(diǎn)1～5溫度隨時(shí)間的變化

② 蘋(píng)果切片內(nèi)部含水率

觀測(cè)點(diǎn)1～5含水率隨時(shí)間的變化見(jiàn)圖7。由圖7可知，蘋(píng)果切片內(nèi)部含水率先下降然后趨于穩(wěn)定，迎風(fēng)面及邊緣含水率下降速率高于中間位置。這種情況易導(dǎo)致蘋(píng)果切片外表面干燥過(guò)快，發(fā)生硬化甚至產(chǎn)生裂縫。

圖7 觀測(cè)點(diǎn)1～5含水率隨時(shí)間的變化

5 結(jié)論

① 干燥空氣流速、進(jìn)風(fēng)相對(duì)濕度一定時(shí)，增大進(jìn)風(fēng)溫度可有效提高干燥速率。進(jìn)風(fēng)溫度、進(jìn)風(fēng)相對(duì)濕度一定時(shí)，增大干燥空氣流速，對(duì)促進(jìn)蘋(píng)果切片平均含水率快速達(dá)到穩(wěn)定的作用并不明顯。進(jìn)風(fēng)溫度、干燥空氣流速一定時(shí)，減小進(jìn)風(fēng)相對(duì)濕度可有效提高干燥速率。

② 在對(duì)流干燥過(guò)程中，蘋(píng)果切片內(nèi)部溫度隨時(shí)間延長(zhǎng)逐漸升高直至穩(wěn)定，表面溫度先于內(nèi)部達(dá)到穩(wěn)定。

③ 干燥是由外向內(nèi)逐漸緩慢進(jìn)行的，模擬中蘋(píng)果切片外表面與中心存在溫差，而蒸發(fā)與溫度密切相關(guān)，因此切片中心與切片外表面的含水率不同。