楊鵬 郭帥帥 賀素艷 趙有信 馬永志 尹義金

摘要： ?針對低溫冷庫進行開門存取貨物時，存在庫門處與外界環(huán)境發(fā)生強烈的熱濕交換，并且隨著時間的增加造成庫內(nèi)溫度較大波動而影響存儲貨物質(zhì)量的問題。本文利用ANSYS Fluent有限元分析軟件，對初始庫溫為243 K冷庫存取貨物開門過程進行溫度場分析，同時對存取結(jié)束后庫內(nèi)降溫過程進行溫度場分析。研究結(jié)果表明，在冷庫庫門開啟30 s，冷庫與庫外環(huán)境熱交換比較強烈，庫內(nèi)溫度波動較大，且大部分在260 K以上。冷庫冷風機開始工作前10 s，冷庫上部降溫速度要比下部快，在10 s之后冷庫上部和下部降溫速度放緩。冷風機工作30 s時，庫內(nèi)溫度不能降低到開門以前的溫度，為了消除開門對冷庫溫度波動的影響，冷風機工作時間必須在30 s以上。冷庫開門時外界環(huán)境的熱濕空氣會因密度差而聚集在冷庫上部，造成冷庫上部溫度高，下部溫度低的現(xiàn)象，為了消除溫度波動帶來的影響，庫內(nèi)冷風機需要工作較長時間。該研究對有效控制冷庫溫度的波動具有重要意義。

關(guān)鍵詞： ?低溫冷庫; 溫度波動; 數(shù)值模擬

中圖分類號： TB661 ?文獻標識碼： A

隨著人們生活水平的提高以及醫(yī)療水平的不斷發(fā)展，冷庫憑借可以大規(guī)模儲存冷凍食品和生物樣本的優(yōu)勢成為低溫存儲的最佳選擇[1]。在冷庫使用過程中，由于冷庫開關(guān)門、庫體保溫強度、庫溫控制方式和制冷系統(tǒng)匹配度等諸多因素都會造成庫內(nèi)溫度的波動，工作人員頻繁進出也是造成冷庫溫度波動的主要因素[24]。冷庫溫度不均勻性加大，會對溫度變化敏感的食品或藥品帶來不可逆的損壞。濕空氣進入到低溫冷庫，庫內(nèi)溫度在空氣露點溫度以下時還會造成濕空氣在蒸發(fā)器表面結(jié)霜，降低換熱器的傳熱效率[56]，結(jié)霜嚴重時還會影響庫內(nèi)空氣的流通。為研究影響冷庫溫度波動的因素，蒲亮[7]以實驗數(shù)據(jù)為基礎(chǔ)，建立了冷庫溫度場和流場數(shù)學模型，討論了冷風機送風速度、送風溫度、冷風機位置、貨物位置及貨物堆積方式對冷庫溫度均勻性的影響。而冷庫開關(guān)門對庫內(nèi)溫度波動影響的研究也是近年來研究的重點，顧成華等人[8]利用VALSALA溫濕度變送器對疫苗冷庫在開關(guān)門過程的冷庫溫度場進行測量，分析結(jié)果顯示，對于小型疫苗冷庫，開關(guān)門造成的溫度波動最大值小于5 K，而通過改善冷風機布置方式，可以達到疫苗存放的要求;張川等人[9]研究了冷庫開關(guān)門頻率對溫度波動的影響，結(jié)果表明每小時開關(guān)門每減少一次，系統(tǒng)整體能耗就會降低20%左右。除了利用實驗方法對冷庫溫度場進行測量，運用計算機模擬的方法對冷庫進行溫度場分析也越來越普遍，謝晶等人[1011]利用Fluent有限元軟件，對冷庫開門過程進行穩(wěn)態(tài)和非穩(wěn)定狀態(tài)下的模擬分析，分析結(jié)果表明，冷庫在開關(guān)門的過程中，庫內(nèi)溫度場發(fā)生強烈波動，在庫門處和庫門上部波動更加強烈，當冷庫關(guān)門之后，庫內(nèi)溫度需要很長時間才能恢復(fù)，但該文并沒有對關(guān)門之后庫內(nèi)的溫度場變化進行分析?；谀壳按蟛糠盅芯考性陂_關(guān)門對冷庫溫度的影響上，對開門過程溫度的非穩(wěn)態(tài)變化過程研究較少，同時對冷庫關(guān)門之后非穩(wěn)態(tài)降溫過程的研究也是空缺。因此，本文基于模擬軟件Fluent，對低溫冷庫在開門過程中內(nèi)部與外界環(huán)境的非穩(wěn)態(tài)溫度變化過程進行模擬研究，同時對存取貨物結(jié)束之后，冷風機工作庫內(nèi)溫度的變化進行分析，研究冷庫開關(guān)門后溫度變化規(guī)律。該研究為有效控制冷庫溫度波動提供了理論依據(jù)。

1 模型構(gòu)建

1.1 物理模型

本研究以海爾中試冷庫機組為實物進行模型繪制。保溫庫板選擇150 mm厚的聚氨酯泡沫夾心隔熱板，并在冷庫底部鋪設(shè)槽鋼，起防凍隔熱作用。實際上，庫體材料導熱系數(shù)很小，因此忽略當冷庫門關(guān)閉時與外界環(huán)境進行的熱量交換[1213]。本文重點研究開關(guān)門對冷庫溫度波動的影響，而忽略冷庫內(nèi)部結(jié)構(gòu)對氣流分布的影響。模擬所用冷庫實物圖如圖1所示。冷庫模型尺寸設(shè)置如表1所示。

為了有效表達冷庫與外界環(huán)境進行的熱量交換，在冷庫外部搭建一個長寬高為3 000 mm×1 500 mm×2 700 mm的環(huán)境區(qū)域。在Gambit軟件中進行模型的構(gòu)建及網(wǎng)格的劃分，冷庫模型如圖2所示，網(wǎng)格劃分圖如圖3所示。將劃分網(wǎng)格之后的模型導入Fluent中進行仿真計算。

1.3 邊界條件

設(shè)置外界環(huán)境空氣溫度為300 K，壁面采用第一類邊界條件，且屬于氣體無滑移狀態(tài)。設(shè)置冷庫未開門時庫內(nèi)溫度為243 K，壁面條件采用第一類邊界條件，并默認壁面無滑移[19]。當冷庫門打開時，將冷庫門邊界類型設(shè)為interior;當取貨完成后，冷庫門設(shè)為wall邊界。在進行冷庫開門溫度模擬時，其屬于封閉腔內(nèi)的自然對流，符合Boussinesq假設(shè)條件，以便于計算溫差引起的浮升力項[2021]。風機統(tǒng)一采用速度進口邊界類型，風速為5 m/s，初始溫度為243 K。湍流強度為10%，水力直徑設(shè)定為0.5 m，回風口邊界類型設(shè)置為outflow，計算過程考慮重力因素，所以勾選重力條件[22]。

2 模擬結(jié)果及分析

2.1 庫門開啟時間對冷庫溫度波動的影響

三維模擬結(jié)果很難清晰得到冷庫內(nèi)部溫度變化規(guī)律，因此為準確對結(jié)果進行分析，選取冷庫中的一個代表性截面（在庫門中間處取一個縱向）作為參考面，并分別將庫門開啟10，20，25，30 s時，對冷庫附近溫度場進行模擬和分析。庫門開啟時冷庫與外界溫度等值線如圖4所示。

由圖4可以看出，當冷庫門開啟時，冷庫與外界高溫環(huán)境發(fā)生激烈熱量交換，庫內(nèi)溫度發(fā)生較劇烈的波動。由圖4a可以看出，當庫門打開10 s時，外界環(huán)境的高溫高濕空氣進入冷庫后，因為密度較小往冷庫上部進行流動。庫門處的溫度交換最劇烈，庫門上部的溫度最高達到外界環(huán)境的300 K，庫頂最高溫度達到275 K，同時在開啟10 s時，庫內(nèi)溫度大部分還維持在較低水平;由圖4b可以看出，當庫門開啟20 s時，庫內(nèi)原始低溫區(qū)域面積減少一半，庫頂最高溫度達到285 K，庫內(nèi)溫度大部分在260～270 K之間。開門時間繼續(xù)延長，當達到25 s時，庫內(nèi)溫度基本維持在265～285 K之間，原始低溫區(qū)域面積繼續(xù)減小。此時，對于庫頂位置，因為平均溫度升高約37 K，會對存放的貨物造成嚴重損害。當庫門開啟30 s時，庫內(nèi)的平均溫度繼續(xù)升高，原始低溫區(qū)域進一步減小，庫底的溫度也隨之升高。由此可以看出，當冷庫開啟30 s及更長時間后，庫內(nèi)溫度已不能滿足存儲貨物的要求，為了減小庫門開啟對庫內(nèi)溫度波動的影響，應(yīng)將冷門開啟時間控制在20 s以內(nèi)。

2.2 風機開啟時間對冷庫溫度波動的影響

在冷庫開啟期間，庫內(nèi)的溫度不斷升高，極大損害冷庫內(nèi)存儲貨物的使用價值，所以當庫門開啟30 s之后關(guān)閉庫門，此時將冷庫門的邊界條件類型由interior改為wall，并對冷風機開啟時對庫內(nèi)溫度的影響進行模擬研究。同時，為了更好地觀察庫內(nèi)溫度變化規(guī)律，在Z=1 000 mm處進行橫截面，Z=1 000 mm截面溫度等值線如圖5所示。

由圖5a可以看出，當冷庫門開啟，冷風機未開始工作時，冷庫Z=1 000 mm橫截面，即冷庫下部區(qū)域溫度都在260 K以上。同時，庫外溫度比庫內(nèi)溫度低，說明庫內(nèi)外氣體還要進行進一步的熱量交換，此時關(guān)閉冷庫門，將庫內(nèi)冷風機打開;由圖5b可以看出，當冷風機開啟10 s之后，在冷庫Z=1 000 mm截面，溫度變化不是特別明顯，這是因為冷風機位于冷庫上部，開啟時由于冷風機的直吹作用，首先對冷庫上部的氣體進行降溫，并且在冷庫中形成溫度逐漸增高的同心圓等溫線。當冷風機開啟20 s之后，Z=1 000 mm截面冷庫內(nèi)的溫度基本降到260 K以下，當時間達到30 s時，溫度降至254～256 K之間，一方面是由冷庫冷風機的降溫作用，使低溫空氣由于密度大而流向冷庫底部;另一方面，因為冷庫底部的溫度原本就比較低，所以冷庫底部的降溫效果比較明顯。

同樣，為了更好地觀察庫內(nèi)溫度變化規(guī)律，在Z=2 500 mm處進行橫截面，Z=2 500 mm截面溫度等值線圖如圖6所示。

由圖6a可以看出，截面平均溫度比Z=1 000 mm截面的平均溫度高。當風機開啟10 s時，溫度基本維持在260～265 K之間，與圖6b相比，降溫效果比較明顯，這是因為Z=2 500 mm截面正好位于冷風機處，受冷風機直接降溫的影響比較大;由圖6c和圖6d可以看出，當冷風機開啟20 s和30 s之后，其降溫效果不明顯，溫度同樣在255～265 K之間，與圖6b不同的是，溫度在260 K以下的區(qū)域面積增大。

由圖5和圖6可以看出，風機開啟時，庫底溫度降低較明顯，庫頂降溫比較緩慢，當經(jīng)過冷風機工作30 s之后，庫底溫度可以降到256 K以下，而庫頂?shù)臏囟葏s只能降到260 K左右，所以要想達到庫內(nèi)溫度均勻，需要將冷風機工作時間延長至30 s以后。

3 結(jié)束語

本文主要利用有限元分析軟件Fluent，對冷庫開門過程庫內(nèi)溫度的變化及關(guān)門之后冷庫冷風機工作對庫內(nèi)溫度的影響進行模擬仿真。仿真結(jié)果表明，冷庫開門初期庫內(nèi)冷空氣與外界熱濕空氣在門口處發(fā)生強烈的熱量和質(zhì)量交換，外界溫度較高的氣體會大量進入冷庫，并且因為密度差的原因，氣體聚集在冷庫上部，使冷庫上部溫度升溫劇烈，對位于冷庫上部的存儲貨物影響較大。在冷庫使用過程中，可以降低冷庫上部的貨物堆積密度，將貨物集中在貨架下部存放，同時將開門時間控制在20 s以內(nèi)，最大限度降低開門對冷庫溫度波動的影響。冷風機工作初始時刻，冷庫上部降溫比較快，下部因為初始溫度比較低，降溫速度比較慢。對于冷庫后期降溫緩慢的現(xiàn)象，可以適當增加冷風機轉(zhuǎn)速，增強庫內(nèi)氣體流動。

參考文獻：

[1] 李藝哲， 謝晶. 小型冷庫內(nèi)溫度場的數(shù)值模擬[J]. 山東農(nóng)業(yè)大學學報： 自然科學版， 2017， 48（4）： 525527.

[2] 賴偉彬， 劉文利， 陳偉群. 小型風冷型裝配式冷庫熱回收改裝[J]. 制冷， 2014， 33（2）： 7982.

[3] 楊磊. 微型保鮮庫氣體流場的數(shù)值模擬與試驗研究[D]. 江蘇： 南京農(nóng)業(yè)大學， 2008.

[4] 吳天. 冷庫內(nèi)氣體流場分析及優(yōu)化設(shè)計[D]. 上海： 上海水產(chǎn)大學， 2006.

[5] 趙文鋒， 楊洲. 微型冷庫貨物降溫特性的CFD數(shù)值模擬[J]. 中國科技論文， 2015， 10（5）： 546551.

[6] 繆晨， 謝晶. 冷庫空氣幕流場的非穩(wěn)態(tài)數(shù)值模擬及驗證[J]. 農(nóng)業(yè)工程學報， 2013， 29（7）： 246253.

[7] 蒲亮. 冷庫溫度場的實驗研究及數(shù)值模擬[D]. 陜西： 西安交通大學， 2003.

[8] 顧呈華， 張皓月， 張忠斌， 等. 小型疫苗冷庫內(nèi)部溫度場試驗研究[J]. 中國設(shè)備工程， 2020（19）： 7779.

[9] 張川， 司春強. 開門頻率對餐廳用室內(nèi)裝配式冷庫溫度及能耗影響[J]. 冷藏技術(shù)， 2019， 42（3）： 2833.

[10] 謝晶， 吳天. 小型冷庫開關(guān)門過程溫度場的數(shù)值模擬[J]. 上海水產(chǎn)大學學報， 2006， 15（3）： 333339.

[11] 翟玉玲， 鐘曉暉， 勾昱君， 等. 小型冷庫開關(guān)門過程的數(shù)值模擬[J]. 河北理工大學學報： 自然科學版， 2011， 33（2）： 144149.

[12] 芮文琴. 小型裝配式冷庫庫內(nèi)氣流組織與數(shù)值模擬[D]. 河北： 華北電力大學， 2012.

[13] 湯毅. 計算流體力學（CFD）技術(shù)用于冷庫氣流優(yōu)化的研究[D]. 上海： 上海海洋大學， 2013.

[14] 王福軍. 計算流體動力學分析一CFD軟件原理與應(yīng)用[M]. 北京： 清華大學出版社， 2004.

[15] 陶文銓. 數(shù)值傳熱學[M]. 2版. 西安： 西安交通大學出版社， 2002.

[16] 楊強生， 浦保榮. 高等傳熱學 [M]. 2版. 上海： 上海交通大學出版社， 2001.

[17] Hu Z H， Sun D W. CFD simulation of heat and moisture transfer for predicting cooling rate and weight loss of cooked ham during air-blast chilling process[J]. Journal of Food Engineering， 2000， 46（3）： 189197.

[18] Mirade P S， Daudin J D. A numercial study of the air flow patterns in a sausage dryer[J]. Drying Technology， 2000， 18（1/2）： 8197.

[19] 張銘， 艾宇， 張少杰， 等. 冷庫不同布置方案對冷庫溫度均勻性的影響[J]. 江蘇農(nóng)業(yè)科學， 2020， 48（2）： 210221.

[20] 張俊杰， 安毓輝. 船舶伙食冷庫內(nèi)氣流組織模擬與分析[J]. 船舶工程， 2019， 41（11）： 5661.

[21] 田甜， 李福良， 張雯， 等. 冷風機不同布置位置對冷庫氣流組織影響的模擬及實驗研究[J]. 冷藏技術(shù)， 2019， 42（1）： 1823.

[22] 王金鋒， 李文俊， 謝晶， 等. 速凍機寬度方向上流場的實驗分析[J]. 熱能動力工程， 2018， 33（4）： 110114.