段潔利，呂恩利，邱 漢，張馨予，趙文鋒，劉 偉，陸華忠,2

(1.華南農業(yè)大學 a.工程基礎教學與訓練中心；b.工程學院；c.電子工程學院，廣州 510642；2.廣東省農業(yè)科學院，廣州 510640)

0 引言

空氣濕度控制是倉儲生產中的重要環(huán)節(jié)，尤其在南方多濕條件下，倉儲除濕可有效保證物資貯藏質量[1]。谷物及種子對貯藏環(huán)境濕度有較為嚴格的要求，種子的含水量直接受環(huán)境中相對濕度的影響，控制環(huán)境的相對濕度、降低種子的含水量可適當保持種子活性，延長種子貯藏壽命[2-4]。種子貯藏環(huán)境的降濕比降溫相對較易, 造價也低。目前，科研人員大多采用CFD軟件模擬庫體內溫度場和速度場的分布，而有關除濕過程中的濕度場分布情況及貯藏物水分影響因素的文獻較少。COMSOL Multiphysics以計算性能高、能進行多物理場耦合等優(yōu)點近年來得到較多的關注[5-6]。Huang Zhi等人[7]采用COMSOL模擬大豆粉的干燥過程。許明杰等[8]采用COMSOL建立一個耦合CO變換反應傳質、傳熱、動量傳遞和化學反應的二維軸對稱擬均相固定床反應器模型。因此，基于COMSOL Multiphysics可以進行流場數值模擬，了解除濕過程中濕度的變化情況和分布規(guī)律。

本文以基于壓差原理的貯藏箱體為研究對象，以雜交水稻種子為試驗物料，基于COMSOL Multiphysics 5.0，通過k-ε湍流模型、多孔介質模型和稀物質傳遞模型等對除濕過程濕度場和不同包裝種子在不同濕度場中不同時間的含水率進行數值模擬，得出不同包裝種子在不同濕度下貯藏7天后的濕度場云圖。

1 物理模型

除濕系統試驗平臺結構如圖1所示。箱體尺寸為1.60m×1.10m×1.50m(長×寬×高)。在風機作用下，箱內空氣由頂部回風道進入制冷室，通過蒸發(fā)器，將降溫氣體通過開孔隔板進入貯藏室，實現降溫循環(huán)；同時，空氣在除濕機作用下也沿此循環(huán)。抽濕機(DY-618LV/A德業(yè)除濕機，廣州昊德儀器有限公司)靠近開孔隔板，通過箱體外控制電箱調節(jié)貯運室中的相對濕度；試驗平臺主要由執(zhí)行機構、可編程控制器、溫度傳感器及相對濕度傳感器等組成。通過傳感器監(jiān)測箱體相應的溫度和相對濕度，可編程控制器控制執(zhí)行機構的開啟與關閉，從而根據試驗要求調節(jié)箱內的環(huán)境，實現智能化。其中，溫度傳感器精度±0.3℃，測量范圍-20℃～80℃；相對濕度傳感器精度±3%，測量范圍0～100%。

試驗物料為雜交水稻種子，種子外形完好，無病蟲害，色彩光澤正常。設定箱體貯藏溫度為15℃，設置5個濕度水平(W1為<35%±3%,W2為50%±3%,W3為65%±3%，W4為80%±3%，W5為>95%±3%)，即5個相同的試驗平臺同時運行，采用不銹鋼盤(散裝)、牛皮紙袋、編織袋、塑料袋及鋁箔袋5種包裝形式。每個試驗平臺內有1個0.83m×0.64m×0.21m(高×長×寬)的5層貨架，每層放裝1kg的2袋種子，每層包裝不同。種子袋布置圖如圖2所示。每天同一時間在不同濕度的平臺內從不同包裝的袋內各取樣一次，測量其含水率作為試驗值。根據邵學良文獻[9]，種子含水率選用整粒法來測量，國外文獻報道已普遍使用此方法進行測試[10]。種子初始含水率為11.28%，試驗中對各濕度環(huán)境的能耗進行了統計。

1.回風道 2.壓差式箱體 3.進氣閥 4.繼電器盒 5.可編程控制器 6.無紙記錄儀 7.冷凝器 8.壓縮機 9.排水管 10.補水箱 11.進水管 12.積水槽 排氣閥 13.加濕器 14.開孔隔板 15.檢測溫度傳感器 16.汽化盤管 17.蒸發(fā)器 18.風機安裝版 19.風機 20.傳感器盒 21.抽濕機 22.排氣閥 23.氣流導軌 24.抽濕機控制箱

1.貯藏室 2.貨架 3.不銹鋼盤裝 4.牛皮紙袋包裝 5.編織袋包裝 6.塑料袋包裝 7.鋁箔袋包裝

選取箱體所在中縱截面，通過AutoCAD進行二維建模，運用COMSOL Multiphysics5.0軟件對模型結構進行網格劃分，采用三角形網格和自由網格劃分方法離散求解域，如圖3所示。

圖3 貯藏箱體網格模型圖Fig.3 Grid model of storage and transportation chamber

2 數學模型

為簡化模型和方便計算，對模型進行了如下假設[11-14]：

1)假設氣體為牛頓流體；

2)箱體內部氣體在流動過程中不可壓縮，符合Boussinesq假設；

3)假設該箱體結構密閉絕熱。

通過k-ε模型、多孔介質模型和稀物質傳遞模型，對各不同濕度箱體的流場變化進行仿真計算?？刂品匠滩捎糜邢摅w積法的通用控制方程，稀物質傳遞模型擴散系數為e-9m2/s，其不同包裝多孔介質模型擴散系數如表1所示，k-ε模型經驗常數如表2所示[15]。

表1 多孔介質模型流體擴散系數Table 1 The fluid diffusion coefficient of porous media model m2/s

表2 k-ε模型中采用的經驗常數Table 2 Coefficients of model k-ε

3 邊界條件設置及求解

根據求解公式[15]對貯藏箱體進行雷諾數計算，得到雷諾數Re均達到103以上，該模型適用雷諾數紊流模型[16]。

Re=UL/ηk

(1)

其中，U為截面平均速度(m/s)；L為特征長度(m)；ηk為流體運動學黏性系數(m2/s)。

3.1 入口邊界條件

箱體內除濕機工作時，從開孔隔板中部出來進入箱體貯藏室的是干空氣，定義此口為入口。入口速度為7 m/s，溫度15℃。輸入紊流參數為湍流強度I和水力直徑DH。湍流強度I計算公式為[17]

I=0.16(ReDH)-1/8

(2)

式中ReDH—以水力直徑求出的雷諾數。

干空氣入口條件采用質量流量，通過試驗測得15℃時的除濕量約為0.014g/s。參照文獻[13]方法，將相對濕度轉化為水蒸氣的摩爾濃度輸入。

3.2 出口條件

濕空氣順著頂部回風道流出貯藏室，定義回風道口條件為壓力出口。

3.3 貨物單元

試驗物料為不同包裝的水稻種子，為簡化運算，文章將種子和包裝袋看作是一個貨物單元，并令該單元為多孔介質模型，貨物內部為層流模式。物性參數值如表3所示。

表3 物性參數設置Table 3 Setting of material parameters

3.4 求解設置

運用COMSOL Multiphysics瞬態(tài)求解器進行求解，設置環(huán)境重力加速度為9.8m/s2，設置時間步長1s。結合試驗，設回風口風速為7m/s。在貯藏室前部開孔隔板相應尺寸處設置濕度監(jiān)測點，當相對濕度達到35%時停止迭代[18]。種子含水率模擬設置環(huán)境重力加速度為9.8m/s2，設置時間步長0.2d，貯藏時間到達7d時停止迭代計算。

4 模擬結果分析與試驗驗證

4.1 箱體內不同包裝不同相對濕度種子含水率與濕度場云圖分布

4.1.1不同相對濕度環(huán)境對相同包裝水稻種子含水率的影響

設定貯運環(huán)境溫度為15℃，所有種子都放在箱體內貨架上的不銹鋼盤中，依次改變相對濕度環(huán)境，通過comsol 迭代計算，到達7d時停止迭代計算。不銹鋼盤裝的種子在不同相對濕度下種子含水率同濕度場云圖如圖4所示。

由圖4(a)、(b)、(c)可以明顯看出：環(huán)境相對濕度越大，包圍在種子周圍的水氣濃度越大，與盤內的低水分含量的種子形成一定的水勢差；隨著時間的變化，環(huán)境中的水分就會在此水勢差作用下進入到種子內部，導致種子含水量升高。從圖4(d)、(e)可看出：隨著相對濕度的降低，種子周圍的水氣含量逐漸降低；當低于種子內部的水分含量時，種子會在水勢差的作用下向周圍環(huán)境失水，從而使種子的含水率逐漸下降；50%以下相對濕度環(huán)境下，不銹鋼盤裝的種子能保持自身的含水率而變化較小。當然，低濕干燥環(huán)境利于種子的保存。試驗表明：環(huán)境相對濕度越低，抽濕的能耗越高。綜合考慮， 50%的相對濕度環(huán)境綜合效果較好。

4.1.295%高濕環(huán)境對不同包裝水稻種子含水率的影響

為了模擬南方的高濕環(huán)境(95%以上)，采用上述同樣的方法迭代計算，得其對應的濕度云圖,如圖5所示。由圖5可知：高濕環(huán)境下，盡管鋁箔袋和塑料袋周圍有很高的水汽，但包裝內的種子仍能保持較低的含水率。這是因為兩種包裝對水分有較好的阻隔作用，在此阻擋下包裝內外的水勢差難以形成通路，包裝內種子與環(huán)境間的水分交換少，導致種子含水率變化較小。從模擬數據可知，平均吸濕速率僅有0.2×10-2mol/d和0.36×10-2mol/d。不銹鋼盤裝的種子周圍的水汽含量較大，由于敞開著口對水分的阻礙小，種子會在水勢差的作用下與環(huán)境間的水分進行較多的交換，導致盤內種子的含水率升高，其平均吸濕速率可達1.79×10-2mol/d。在盤界面處出現的低濕是由于模擬是在靜止狀態(tài)下進行，種子吸收界面的水分而使界面處相對濕度較小還沒來得及和外界進行交換造成的。因此，高濕環(huán)境下，鋁箔袋和塑料袋包裝可以較好地阻隔濕環(huán)境中的水分進入種子內部。

4.1.335%超干環(huán)境對不同包裝水稻種子含水率的影響

對35%超干燥環(huán)境進行迭代計算，其對應的濕度云圖如圖6所示。

由圖6可知：不銹鋼盤、牛皮紙袋包裝的種子顏色與環(huán)境顏色沒有形成明顯的對比，種子內部的水分較多地與周圍環(huán)境進行了交換進入到周圍環(huán)境中；界面處之所以水分含量較高，是因為在貯運過程中種子內部釋放的水分到達界面處由于接觸界面面積較小暫時沒來得及與周圍環(huán)境進行交換造成的；種子的含水率總體表現在逐漸下降，其兩種包裝的平均除濕速率分別是0.68×10-2mol/d、0.45×10-2mol/d。鋁箔袋和塑料袋裝的種子含水率基本保持原來狀態(tài)，水分被隔絕在包裝袋內，沒有和外界環(huán)境進行交換，其阻隔率達97.06%。但無論何種包裝，種子的含水率都仍在12%安全貯藏范圍內。因此，超干環(huán)境下包裝對種子含水率的影響較小。貯運7天后發(fā)現：此箱的除濕機耗電達91.75kW·h，而50%相對濕度耗電僅為22.04kW·h，前者耗能大概是后者的4倍多。

4.2 試驗驗證

為驗證模擬結果是否準確，通過調節(jié)箱體貯運室內溫度為15 ℃，濕度為對應的5種濕度，采用5種不同的包裝，進行了不同包裝、不同濕度環(huán)境種子含水量的試驗。每天測量1次各種環(huán)境下種子的含水率，用無紙記錄儀監(jiān)控記錄貯運過程中縱截面上各點的濕度變化。從模擬結果中取出測試點對應的種子含水率值，與試驗值進行對比，結果如圖7所示。

由圖7(a)可知：貯藏過程中，散裝種子在不同濕度環(huán)境下模擬值與試驗值隨時間變化的規(guī)律基本吻合，模擬值與試驗值種子含水率最大偏差值僅為1%；在50%以上的相對濕度環(huán)境中，從第1天到第5天，種子含水率模擬值與試驗值吻合度較好，最大偏差值小于0.4%。環(huán)境濕度越大，種子的吸水速率會越大，50%相對濕度環(huán)境下種子的含水率處于一個相對平衡狀態(tài)。

由圖7(b)可知：高濕環(huán)境下，包裝對種子含水率的影響較大，不同包裝種子含水率測試點的模擬值與試驗值隨時間變化規(guī)律基本一致，模擬值與試驗值的最大偏差值僅為0.5%；不銹鋼盤裝的種子平均吸濕速率能達到1.79×10-2mol/d，鋁箔袋包裝的種子平均吸濕速率僅有0.2×10-2mol/d，兩者相差約9倍。

由圖7(c)可知：超干環(huán)境下，包裝對種子含水率的影響沒有高濕環(huán)境下明顯，測試點種子含水率模擬值與試驗值出現了一定的偏差，但最大偏差也僅為1.2%。由于塑料袋和鋁箔袋的阻隔作用，種子含水率降低速率明顯低于其余3種包裝。塑料袋和鋁箔袋包裝種子含水率試驗值稍偏高的原因是：除濕機的不停運轉使箱體內的空氣溫度升高，導致這兩種包裝的袋子內微小環(huán)境溫度升高，種子內的水氣不能及時排出袋子外面，致使種子含水率略有些偏高，而其余3種包裝直接與干燥空氣進行了交換出現含水率偏低的情況。

5 結論

1)高濕環(huán)境下包裝對種子含水率的影響較大。95%高濕環(huán)境下散裝的種子平均吸濕速率可達1.79×10-2mol/d，鋁箔袋包裝平均吸濕速率僅有0.2×10-2mol/d，是散裝包裝的1/9；高濕環(huán)境下宜采用水分阻隔率較大的鋁箔袋或塑料袋包裝，利于種子的貯藏。

2)35%超干燥環(huán)境下，散裝的種子平均解吸速率可達0.68×10-2mol/d，貯運7天后含水量降到9%左右；而50%相對濕度環(huán)境水稻種子的含水率基本穩(wěn)定在11.5%，前者能耗是后者能耗的4倍多。因此，50%相對濕度環(huán)境適于種子貯藏。

3)從模擬值與試驗值對比來看，不同濕度環(huán)境、不同包裝種子含水率模擬值與試驗值隨時間變化的規(guī)律基本吻合，模擬值與試驗值最大偏差僅為1.2%。由此表明，該模型是有效的。