裝設(shè)均勻送風(fēng)管道對冷藏庫氣流流場特性的改善

劉曉菲，南曉紅（西安建筑科技大學(xué)環(huán)境與市政工程學(xué)院，西安710055）

?

裝設(shè)均勻送風(fēng)管道對冷藏庫氣流流場特性的改善

劉曉菲，南曉紅※
（西安建筑科技大學(xué)環(huán)境與市政工程學(xué)院，西安710055）

摘要：果蔬冷藏環(huán)境對庫內(nèi)溫濕度及速度場的均勻性有較高要求，采用冷風(fēng)機(jī)加均勻送風(fēng)管道的形式可以更好地滿足果蔬貯藏的條件。該文針對某200 t蘋果冷藏庫，研究了庫內(nèi)裝設(shè)與不裝設(shè)均勻送風(fēng)管道時庫內(nèi)的溫度場、速度場及相對濕度場的分布情況。文中建立了冷風(fēng)機(jī)（加均勻送風(fēng)管道）-貨物-庫內(nèi)空氣環(huán)境的氣體流動、傳熱與傳質(zhì)的三維耦合數(shù)學(xué)求解模型。在冷風(fēng)機(jī)送風(fēng)速度為5.24 m/s，送風(fēng)溫度為271.15 K，送風(fēng)相對濕度為90%的條件下，數(shù)值模擬研究了冷藏庫內(nèi)速度場、溫度場及相對濕度場，并通過計(jì)算不均勻系數(shù)及空氣分布特性指標(biāo)（air diffusion performance index, ADPI）對流場特性進(jìn)行了評價。結(jié)果顯示，當(dāng)庫內(nèi)裝設(shè)均勻送風(fēng)管道時貨物區(qū)速度不均勻系數(shù)比不設(shè)送風(fēng)管道時降低22.65%，溫度不均勻系數(shù)降低20%，相對濕度不均勻系數(shù)降低22.73%，空氣分布特性指標(biāo)增大11.13%。因此，裝設(shè)均勻送風(fēng)管道時庫內(nèi)速度場、溫度場及相對濕度場均更加均勻。

關(guān)鍵詞：流場；冷藏庫；均勻送風(fēng)管道；氣流組織；不均勻系數(shù)

劉曉菲，南曉紅.裝設(shè)均勻送風(fēng)管道對冷藏庫氣流流場特性的改善[J].農(nóng)業(yè)工程學(xué)報(bào)，2016，32（01）：91-96.doi：10.11975/j.issn.1002-6819.2016.01.012 http://www.tcsae.org

Liu Xiaofei, Nan Xiaohong.Improvement on characteristics of air flow field in cold storage with uniform air supply duct [J].Transactions of the Chinese Society of Agricultural Engineering（Transactions of the CSAE）, 2016, 32（01）: 91－96.（in Chinese with English abstract）doi：10.11975/j.issn.1002-6819.2016.01.012 http://www.tcsae.org

0　引言

冷藏庫作為果蔬貯藏的重要設(shè)施，庫內(nèi)溫濕度場的均勻性對果蔬的貯藏品質(zhì)有較重要的影響，庫內(nèi)氣流組織的不均勻會導(dǎo)致果蔬貯藏品質(zhì)的下降并降低冷藏庫企業(yè)經(jīng)濟(jì)效益。所以，如何使得庫內(nèi)溫、濕度均勻合理地保持在最為適宜的貯藏條件是冷藏庫設(shè)計(jì)和運(yùn)行時必須要考慮的問題。在庫內(nèi)設(shè)置配風(fēng)系統(tǒng)能夠獲得更為理想的均勻流場，在冷風(fēng)機(jī)送風(fēng)口裝設(shè)均勻送風(fēng)管道是一種主要方式[1]。

冷藏庫內(nèi)冷卻、加濕設(shè)備的設(shè)計(jì)參數(shù)、運(yùn)行參數(shù)及貨物質(zhì)量、品種以及擺放形式等均會對庫內(nèi)溫濕度場產(chǎn)生影響，庫內(nèi)溫濕度場的分布較為復(fù)雜，國內(nèi)外學(xué)者對庫內(nèi)溫度場及相對濕度場的影響因素做了大量的研究工作。文獻(xiàn)[2-4]建立綜合傳熱傳質(zhì)的三維CFD（computational fluid dynamics）計(jì)算模型，模擬研究了果蔬呼吸熱、孔隙率及堆貨尺寸對庫內(nèi)溫濕度場的影響；劉永娟通過建立三維CFD穩(wěn)態(tài)模型研究了冷風(fēng)機(jī)送風(fēng)速度及貨物擺放形式對庫內(nèi)速度場及溫度場的影響[5]；文獻(xiàn)[6-7]分析了3種不同風(fēng)機(jī)擺設(shè)形式對庫內(nèi)氣流場的影響，得出風(fēng)機(jī)背對吹風(fēng)式為最佳擺設(shè)形式；文獻(xiàn)[8-9]通過建立包含冷庫內(nèi)部對流、冷庫門空氣幕射流及室外環(huán)境風(fēng)場流動在內(nèi)的三維耦合數(shù)值模型研究了空氣幕性能的影響因素及空氣幕運(yùn)行時庫內(nèi)流場的分布情況；文獻(xiàn)[10-11]采用理論計(jì)算和CFD模擬相結(jié)合的方法分析了熱壓作用下的冷庫門的冷風(fēng)滲透率；謝晶等人用計(jì)算流體力學(xué)技術(shù)對冷庫開關(guān)門過程進(jìn)行了非穩(wěn)態(tài)數(shù)值模擬，分析了冷庫開關(guān)門過程中庫內(nèi)外溫度場的變化，研究了室外空氣對庫內(nèi)溫度場的影響情況[12]。

公開發(fā)表的文章鮮有涉及到冷藏庫內(nèi)裝設(shè)與不裝設(shè)均勻送風(fēng)管道時庫內(nèi)速度場、溫度場及相對濕度場分布情況的研究。因此，為了確定在庫內(nèi)裝設(shè)均勻送風(fēng)管道后溫度場和相對濕度場的均勻性及合理性的改善程度，本文以陜西省延安市某200 t的蘋果冷藏庫為研究對象，采用紊流模型和Simple算法對比研究了庫內(nèi)有無均勻送風(fēng)裝置時速度場、溫度場及相對濕度場的分布情況。研究結(jié)果為冷庫的建設(shè)和提高果蔬貯藏品質(zhì)提供理論指導(dǎo)。

1　三維耦合數(shù)學(xué)求解模型

1.1物理模型

本文研究的冷藏庫容量為200 t，冷卻設(shè)備采用落地式冷風(fēng)機(jī)，送風(fēng)裝置為均勻送風(fēng)管道，均勻送風(fēng)管道的總長度為16.1 m，首段風(fēng)管斷面尺寸為2 m×0.6 m，末端風(fēng)管斷面尺寸為0.4 m×0.6 m，均勻送風(fēng)管道的兩側(cè)面上以0.617 m的間隔均勻設(shè)置50個圓錐形送風(fēng)口，風(fēng)口的直徑為0.18 m。貨物采用兩堆放置，貨物堆距兩側(cè)墻面0.6 m，貨物堆間距2.2 m，距冷庫前后側(cè)墻面各2 m，距地面0.2 m。計(jì)算模型見圖1。

圖1　設(shè)有均勻送風(fēng)裝置的冷藏庫模型Fig.1 Physical model of cold storage that equipped with uniform air supply duct

作為對比，現(xiàn)將該冷藏庫設(shè)計(jì)另外設(shè)置一套冷卻設(shè)備，依據(jù)冷卻面積為主、風(fēng)量作為參考的原則[13]，為該冷藏庫選擇一臺吊頂式冷風(fēng)機(jī)作為冷卻設(shè)備，該冷風(fēng)機(jī)的參數(shù)如表1所示。冷風(fēng)機(jī)距兩側(cè)墻為6.1 m，距背面墻為0.5 m。庫內(nèi)貨物堆放形式不變。計(jì)算模型見圖2。

圖2　冷卻設(shè)備為吊頂式冷風(fēng)機(jī)的冷藏庫模型Fig.2 Physical model of cold storage without uniform air supply duct

表1　冷卻設(shè)備參數(shù)Table 1 Air-cooler parameters

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

1.2.1控制方程

計(jì)算域內(nèi)流體流動與傳熱的控制方程包括連續(xù)性方程、動量方程、能量方程和k-ε兩方程模型方程，聯(lián)立這些方程可以得到方程組的通用形式如下：

式中φ、Γφ、Sφ分別是通用變量、廣義擴(kuò)散系數(shù)和廣義源項(xiàng)。對應(yīng)于不同的通用變量，式中各項(xiàng)參數(shù)的具體表達(dá)式可參見文獻(xiàn)[14]。

1.2.2貨物區(qū)的處理

針對蘋果冷藏庫，將庫內(nèi)貨物看作多孔介質(zhì)，水果間隙的空氣對應(yīng)于多孔介質(zhì)中的流體。

1）貨物區(qū)多孔介質(zhì)阻力

當(dāng)庫內(nèi)空氣受冷風(fēng)機(jī)送風(fēng)作用向果蔬內(nèi)滲流時要受到水果對其的阻力，這部分阻力在計(jì)算時作為流體動量方程的源項(xiàng)被考慮，其主要表現(xiàn)為黏性阻力和慣性阻力：

式中μ代表流體的動力黏度，Pa·s；α為多孔介質(zhì)的滲透率，m2；C2是慣性阻力系數(shù)，1/m；ρ為流體的密度，kg/m3；vi代表速度在X、Y和Z方向的分量，m/s；vmag是流體的速度，m/s。

依據(jù)Ergun方程[15]可以得到黏性阻力系數(shù)和慣性阻力系數(shù)：

式中Dp為多孔介質(zhì)顆粒的平均粒徑，m；ε為多孔介質(zhì)的孔隙率。

2）果蔬呼吸熱

果蔬在庫內(nèi)貯藏期間一直發(fā)生呼吸作用，經(jīng)過呼吸作用產(chǎn)生的能量有45%轉(zhuǎn)化為生物能維持果蔬的生命活動，另外55%的能量則以呼吸熱的方式釋放[16]。果蔬呼吸熱在計(jì)算時作為能量方程的源項(xiàng)被考慮。針對所研究的蘋果冷藏庫，果蔬的呼吸產(chǎn)熱速率可由文獻(xiàn)[17]查得。

1.2.3邊界條件

合理地設(shè)置邊界條件對問題的研究起著非常重要的作用，針對本文的研究對象，邊界條件進(jìn)行如下的設(shè)置：

1）冷風(fēng)機(jī)送風(fēng)口：速度入口，設(shè)有均勻送風(fēng)裝置的冷藏庫送風(fēng)速度為5.24 m/s，無均勻送風(fēng)裝置的冷藏庫送風(fēng)速度為8.48 m/s，送風(fēng)溫度為271.15 K，送風(fēng)空氣中水蒸氣的質(zhì)量分?jǐn)?shù)為0.002 9（相對濕度為90%）。

2）冷風(fēng)機(jī)回風(fēng)口：自由出流，適用于在計(jì)算完成前無法確定壓強(qiáng)和速度的情況，在計(jì)算時fluent軟件默認(rèn)將除壓強(qiáng)以外的所有流動參數(shù)的法向梯度都設(shè)為0。

3）冷庫壁面熱邊界條件

外墻，第三類邊界條件，外表面與庫外空氣對流換熱系數(shù)為23 W/m2·K，庫外空氣溫度299.25 K，墻體厚度668 mm，墻體材料導(dǎo)熱系數(shù)0.26 W/m·K；

屋頂，第三類邊界條件，屋頂與庫外空氣對流傳熱系數(shù)為23 W/m2·K，庫外空氣溫度299.25 K，屋面厚度314mm，屋頂導(dǎo)熱系數(shù)為0.2 W/m·K；

地面，第二類邊界條件，熱流密度為3 W/m2。

2　計(jì)算模型的試驗(yàn)驗(yàn)證

為驗(yàn)證本文建立的數(shù)值模型的可靠性，對所研究的冷藏庫（冷卻設(shè)備為落地式冷風(fēng)機(jī)加設(shè)均勻送風(fēng)管道）的庫內(nèi)空氣速度和溫度場進(jìn)行了測試，試驗(yàn)測試的條件為：均勻送風(fēng)管道送風(fēng)口的送風(fēng)速度為5.24 m/s，送風(fēng)溫度為-2℃。試驗(yàn)測試用的主要儀器有熱電風(fēng)速儀ZRQFF30J（北京協(xié)亞電子有限公司）、紅外線測溫儀EC-8339（上海舒佳電氣有限公司）。熱電風(fēng)速儀的測速精度為（±0.01）m/s，測速范圍為0.05～30 m/s，測溫精度為（±0.1）℃，測溫范圍為-10～50℃。紅外線測溫儀的測溫精度為（±1）℃，測溫范圍為-32～600℃。試驗(yàn)測試共布置16個測點(diǎn)，測點(diǎn)位置如圖3所示。

圖3　試驗(yàn)測點(diǎn)位置示意圖（距地面1 m的水平面）Fig.3 Positions of measuring points （Distance between measuring points and floor is 1 m）

圖4顯示了測試結(jié)果與模擬計(jì)算結(jié)果的對比情況，從圖中可以看到，距冷風(fēng)機(jī)回風(fēng)口較遠(yuǎn)的測點(diǎn)的空氣的溫度較低；而由于回風(fēng)不斷卷吸周圍果蔬呼吸熱量，距離回風(fēng)口越近的測點(diǎn)溫度越高。測試值與計(jì)算值偏差在15%以內(nèi)的測點(diǎn)有13個，占總測點(diǎn)數(shù)的81.25%，測試值與計(jì)算值偏差在15%以上的測點(diǎn)有3個，占總測點(diǎn)數(shù)的18.75%?？梢?，模擬結(jié)果和現(xiàn)場測試值差別較小，表明本文所建立的數(shù)值模型是可靠的。

圖4　庫內(nèi)測點(diǎn)溫度測試值與模擬值的對比Fig.4 Comparison of predicted temperatures with experimental temperatures in cold storage

3　氣流組織評價指標(biāo)

目前，常用的氣流組織評價指標(biāo)主要包括3方面：送風(fēng)有效性的描述參數(shù)、污染物排除有效性的描述參數(shù)及能量有效利用和熱舒適的描述參數(shù)[18]。對于冷藏庫這類特殊的建筑，人們關(guān)心的主要是對果蔬貯藏品質(zhì)影響較大的因素，如風(fēng)速、溫度及相對濕度等。為了能夠更加合理直觀地評價庫內(nèi)的氣流組織效果，本文采用不均勻系數(shù)及空氣分布特性指標(biāo)對果蔬冷藏庫內(nèi)氣流組織進(jìn)行評價。

3.1溫度不均勻系數(shù)

果蔬貯藏時，人們?yōu)樘岣呃鋷斓慕?jīng)濟(jì)效益，庫內(nèi)一般擺放較多的貨物，容積利用系數(shù)較大，而果蔬又在不斷地進(jìn)行呼吸作用及蒸騰作用，因而庫內(nèi)不同區(qū)域溫度、相對濕度可能存在較大差異，而果蔬貯藏對風(fēng)速、溫度及相對濕度等參數(shù)的均勻性要求較高，所以對庫內(nèi)空氣溫度、相對濕度的均勻性評價顯得尤為重要。鑒于此，本文引入“不均勻系數(shù)”來評價庫內(nèi)氣流組織的均勻性，包括溫度不均勻系數(shù)和相對濕度不均勻系數(shù)，溫度不均勻系數(shù)的概念及計(jì)算公式如下。

溫度不均勻系數(shù)kt、相對濕度不均勻系數(shù)kR（概念和計(jì)算公式均類似于溫度不均勻系數(shù)）都是無量綱數(shù)。kt、kR的值越小，表示溫度場及相對濕度場的均勻性越好。

3.2空氣分布特性指標(biāo)

為整體評價冷庫內(nèi)部環(huán)境滿足果蔬貯藏條件的程度，引入“空氣分布特性指標(biāo)（air diffusion performance index, ADPI）”。通過在冷庫內(nèi)部建立測點(diǎn)，然后測試各個測點(diǎn)的風(fēng)速、溫度及相對濕度的值，則同時滿足果蔬貯藏溫度及相對濕度要求的測點(diǎn)數(shù)與總測點(diǎn)數(shù)的比值即為空氣分布特性指標(biāo)的值[19]。ADPI的值越大，說明滿足果蔬貯藏要求的測點(diǎn)比例越大，則冷庫內(nèi)部的氣流組織效果就越好。其計(jì)算公式如下：

為了使氣流組織評價指標(biāo)的計(jì)算結(jié)果更為客觀準(zhǔn)確，本文在庫內(nèi)貨物區(qū)依長寬高3個方向以1 m為步長均勻取點(diǎn)，共建立780個測點(diǎn)。

4　計(jì)算結(jié)果分析及評價

4.1計(jì)算結(jié)果分析

4.1.1溫度場對比與分析

從圖5a和圖5b可以看出，有均勻送風(fēng)裝置時庫內(nèi)貨物區(qū)的溫度范圍在272.5～273.3 K，貨物堆內(nèi)最大溫差為0.8 K，貨物區(qū)的溫度基本滿足蘋果的貯藏條件（273.15 K，并允許上下波動0.5 K）[1]；當(dāng)庫內(nèi)無均勻送風(fēng)裝置時，貨物區(qū)的溫度范圍為272～273.2 K，貨物區(qū)內(nèi)最大溫差為1.2 K，無均勻送風(fēng)裝置的冷藏庫內(nèi)貨物區(qū)內(nèi)最大溫差較有均勻送風(fēng)裝置時高33%。且無均勻送風(fēng)管道時，庫內(nèi)貨物區(qū)溫度偏低。經(jīng)過以上分析可以發(fā)現(xiàn)，有均勻送風(fēng)裝置時庫內(nèi)貨物區(qū)溫度場更均勻且更接近貯藏條件。

圖5　不同冷卻設(shè)備形式下x=7.67 m截面溫度分布圖Fig.5 Temperature distribution at plane of x=7.67 m under two different conditions

4.1.2相對濕度場對比與分析

從圖6a和6b可以看出，當(dāng)庫內(nèi)有均勻送風(fēng)裝置時貨物區(qū)相對濕度變化范圍為87%～92%，最大相對濕度差值為5%；當(dāng)庫內(nèi)無均勻送風(fēng)裝置時，貨物區(qū)相對濕度變化范圍為87%～94%，最大差值為7%，無均勻送風(fēng)裝置的冷藏庫內(nèi)貨物區(qū)相對濕度最大差值比有均勻送風(fēng)裝置高28.57%。因此，當(dāng)庫內(nèi)設(shè)有均勻送風(fēng)裝置時貨物區(qū)相對濕度更加均勻。

圖6　不同冷卻設(shè)備形式下x=7.67 m截面相對濕度分布圖Fig.6 Relative humidity distribution at plane of x=7.67 m under two different conditions

4.2氣流組織評價指標(biāo)的計(jì)算及分析

在冷藏庫內(nèi)貨物區(qū)以1 m為步長沿3個坐標(biāo)軸方向均勻建立780個測點(diǎn)，調(diào)取每個測點(diǎn)的溫度、風(fēng)速及相對濕度的值，然后按照公式5、公式6及公式7進(jìn)行氣流組織評價指標(biāo)的計(jì)算，計(jì)算結(jié)果如表2所示。

表2　氣流組織評價指標(biāo)計(jì)算結(jié)果Table 2 Value of air distribution evaluation indexes

從表格2可以看出，裝設(shè)均勻送風(fēng)裝置的冷庫和不設(shè)均勻送風(fēng)裝置的冷庫庫內(nèi)溫度和相對濕度均能滿足蘋果貯藏的條件，然而庫內(nèi)貨物區(qū)的溫度和相對濕度的不均勻度有明顯差別：裝設(shè)均勻送風(fēng)裝置的冷藏庫內(nèi)貨物區(qū)溫度不均勻系數(shù)比不設(shè)均勻送風(fēng)裝置的冷藏庫內(nèi)貨物區(qū)溫度不均勻系數(shù)小20%，裝設(shè)均勻送風(fēng)裝置時貨物區(qū)相對濕度不均勻系數(shù)比不設(shè)均勻送風(fēng)裝置時小22.73%，裝設(shè)均勻送風(fēng)裝置時貨物區(qū)滿足規(guī)定風(fēng)速、溫度和相對濕度的測點(diǎn)數(shù)與總測點(diǎn)數(shù)的比值（空氣分布特性指標(biāo)，ADPI）比不設(shè)均勻送風(fēng)裝置時大11.13%。因此，裝設(shè)均勻送風(fēng)裝置時貨物區(qū)溫度場和相對濕度場會更加均勻。

5　結(jié)論

本文以延安市某蘋果冷藏庫為研究對象，建立了庫內(nèi)氣體流動與傳熱的三維數(shù)學(xué)求解模型，依據(jù)冷卻面積、送風(fēng)量及送風(fēng)溫度相同的原則，為該冷庫重新設(shè)置一套冷卻設(shè)備后，對比研究了2套不同冷卻設(shè)備的庫內(nèi)溫濕度場的分布情況，主要得出如下結(jié)論：

1）在對冷藏庫內(nèi)氣流組織進(jìn)行評價時，通過引入溫濕度的平均值、不均勻系數(shù)及ADPI指標(biāo)能夠更加合理直觀地評價冷藏庫內(nèi)溫濕度場的分布情況。

2）采用均勻送風(fēng)裝置時，庫內(nèi)貨物區(qū)的溫度不均勻系數(shù)比無均勻送風(fēng)裝置時小20%，貨物區(qū)的相對濕度不均勻系數(shù)比無均勻送風(fēng)裝置時小22.73%，ADPI指標(biāo)比無均勻送風(fēng)裝置時大11.13%。設(shè)置均勻送風(fēng)裝置時庫內(nèi)溫度場與濕度場的均勻性有較大改善。

[參考文獻(xiàn)]

[1]張祉祐.冷藏與空氣調(diào)節(jié)[M].北京：機(jī)械工業(yè)出版社，1995.

[2] Chourasia M K, Goswami T K.Steady state CFD modeling of airflow, heat transfer and moisture loss in a commercial potato cold store[J].International Journal of Refrigeration, 2007, 30（4）: 672-689.

[3] Chourasia M K, Goswami T K.Simulation of effect of stack dimensions and stacking arrangement on cool-down characteristics of potato in a cold store by computational fluid dynamics [J].Biosystems Engineering, 2007, 96（4）: 503-515.

[4] Ho H Son, Rosario L, Rahman M M.Numerical simulation of temperature and velocity in a refrigerated warehouse [J].International Journal of Refrigeration, 2010, 33（5）: 1015-1025.

[5]劉永娟.貨物不同擺放形式下冷庫內(nèi)氣流組織的模擬研究[D].西安：西安建筑科技大學(xué), 2011.Liu Yongjuan.Simulation of Cold Store′s Air Flow Under Different Layout Forms of Goods[D].Xi′an: Xi′an University of Architecture and Technology, 2011.（in Chinese with English abstract）

[6]謝晶，瞿曉華，徐世瓊.冷藏庫內(nèi)氣體流場數(shù)值模擬與驗(yàn)證[J].農(nóng)業(yè)工程學(xué)報(bào)，2005，21（2）：11-16.Xie Jing, Qu Xiaohua, Xu Shiqiong.Numerical simulation and verification of airflow in cold-store[J].Transactions of the Chinese Society of Agricultural Engineering（Transaction of the CSAE）, 2005, 21（2）:11-16.（in Chinese with English abstract）

[7]湯毅，謝晶，王金鋒，等.CFD預(yù)測風(fēng)機(jī)擺設(shè)形式對冷庫的影響[J].食品與機(jī)械，2012，28（1）：124-128.Tang Yi, Xie Jing, Wang Jinfeng, et al.CFD simulation of air flow distribution with different arrangements of fan though a cold store[J].Food & Machinery, 2012, 28（1）: 124-128.（in Chinese with English abstract）

[8]南曉紅，何媛，劉立軍.冷庫門空氣幕性能的影響因素[J].農(nóng)業(yè)工程學(xué)報(bào)，2011，27（10）：334-338.Nan Xiaohong, He Yuan, Liu Lijun.Impact factors of air curtain performance in cold store[J].Transactions of the Chinese Society of Agricultural Engineering（Transaction of the CSAE）, 2011, 27 （10）: 334－338.（in Chinese with English abstract）

[9]繆晨，謝晶.冷庫空氣幕流場的非穩(wěn)態(tài)數(shù)值模擬及驗(yàn)證[J].農(nóng)業(yè)工程學(xué)報(bào), 2013, 29（7）: 246-253.Miao Chen, Xie Jing.Unsteady state numerical simulation and verification of flow field of air curtain in cold stores [J].Transactions of the Chinese Society of Agricultural Engineering （Transaction of the CSAE）, 2013, 29（7）: 246－253.（in Chinese with English abstract）

[10]何媛，南曉紅.三維CFD模型預(yù)測熱壓作用下冷庫門的冷風(fēng)滲透率[J].農(nóng)業(yè)工程學(xué)報(bào)，2008，24（6）：26-30.He Yuan, Nan Xiaohong.Theoretical study and CFD prediction of infiltration rate through a cold store entrance under thermal pressure[J].Transactions of the Chinese Society of Agricultural Engineering（Transaction of the CSAE）, 2008, 24（6）: 26-30.（in Chinese with English abstract）

[11]翟玉玲，鐘曉暉，勾昱君，等.小型冷庫開關(guān)門過程的數(shù)值模擬[J].河北理工大學(xué)學(xué)報(bào)（自然科學(xué)版），2011，33（2）：144-149.Zhai Yuling, Zhong Xiaohui, Gou Yujun, et al.Numerical simulation of switch door in small cold storage[J].Journal of Hebei Polytechnic University（Natural Science Edition）, 2011, 33 （2）: 144-149.（in Chinese with English abstract）

[12]謝晶，吳天.小型冷庫開關(guān)門過程溫度場的數(shù)值模擬[J].上海水產(chǎn)大學(xué)學(xué)報(bào)，2006，15（3）：333-339.Xie Jing, Wu Tian.Numerical simulation on temperature field in the doorway of a mini type cold store[J].Journal of Shanghai Fisheries University, 2006,15（3）: 333-339.（in Chinese with English abstract）

[13] GB 50072-2010，冷庫設(shè)計(jì)規(guī)范[S].

[14]陶文銓.數(shù)值傳熱學(xué)[M].西安：西安交通大學(xué)出版社，2005.

[15]毛迪凡.孔隙介質(zhì)滲流基本方程的改進(jìn)[D].武漢：中國地質(zhì)大學(xué)，2012.Mao Difan.The Improvement of the Porous Media Seepage Equation[D].Wuhan: China University of Geosciences, 2012.（in Chinese with English abstract）

[16]趙家祿，黃清華.小型果蔬氣調(diào)庫[M].北京：科學(xué)出版社，2000.

[17]張祉祐.氣調(diào)貯藏和氣調(diào)庫[M].北京：機(jī)械工業(yè)出版社，1994.07.

[18]孟彬.體育館氣流組織評價方法[D].北京：清華大學(xué)，2004.Meng Bin.Evaluation Method of Air Distribution in Gymnasium [D].Beijing: Tsinghua University, 2004.（in Chinese with English abstract）

[19]李先庭，趙彬.室內(nèi)空氣流動數(shù)值模擬[M].北京：機(jī)械工業(yè)出版社，2009.

Improvement on characteristics of air flow field in cold storage with uniform air supply duct

Liu Xiaofei, Nan Xiaohong※
（School of Environmental and Municipal Engineering, Xi′an University of Architecture and Technology, Xi′an 710055, China）

Abstract:The cold storage of fruit and vegetable has higher request for the uniformity of temperature, relative humidity and air velocity.The air-cooler with uniform air supply duct can meet the condition of fruit and vegetable storage better.The authors chose a cold storage of apples as a research object and explored the velocity, temperature and relative humidity fields under two different conditions: one used the distribution system and the other did not.Then a combination of CFD （computational fluid dynamics）calculation and experimental study was carried out to further understand the complicated phenomenon of air flow and heat transfer in cold storage.A three dimension physical and numerical model for flow, heat and mass transfer inside the air-cooler（with uniform air supply duct）-cargo-air in the cold storage was developed.The k-ε turbulent model was considered as CFD model.The apple zone was regarded as porous medium zone.When the inlet temperature was 271.15K, inlet velocity was 8.48 m/s and inlet relative humidity was 90%, the velocity, temperature and relative humidity distributions were simulated by CFD software.The simulation results were validated by experimental results.The main instruments for experiment included thermal anemometer ZRQF-F30J and infrared thermometer EC-8399.The total amount of 16 measurement points were tested at following air-cooler inlet parameters: 271.15 K, 8.48m/s and relative humidity 90% , and the air supply system was a uniform air supply duct.The relative difference of the temperature between simulation value and experimental data at 13 measurement points was less than 15%, occupying 81.25% of total measurement points.The numerical results were in good agreement with the measured results.Simulation results showed that the maximum difference in temperature was 0.8 K for the cold storage equipped with a uniform air supply duct and 1.2 K the cold storage without uniform air supply duct.The maximum difference in temperature for the cold storage with uniform air supply duct was 33.33% lower than that without the uniform air supply duct.The maximum difference in relative humidity was 5% for the cold storage equipped with a uniform air supply duct, and 7% when the cold storage without uniform air supply duct.The maximum difference in relative humidity for the cold storage with uniform air supply duct was 28.75% lower than that without the uniform air supply duct.Further, the non-uniformity coefficient and air diffusion performance index（ADPI）were introduced to assess the fields.The total amount of 780 measurement points were selected, which were uniformly distributed and the step length of three coordinate directions was 1 m in the studied cold storage.Both of the non-uniformity coefficient and air diffusion performance index（ADPI）was computational by gathered velocity, temperature and relative humidity of every measurement points.When the uniform air supply duct was installed in the cold storage, the non-uniformity coefficient of air velocity was 22.65% lower than the cold storage without uniform air supply duct.The non-uniformity coefficient of temperature was 20% lower than the cold storage without uniform air supply duct.The non-uniformity coefficient of relative humidity was 22.73% lower than the cold storage without uniform air supply duct.The value of ADPI was 11.13% higher than the cold storage without uniform air supply duct.Thus, the air velocity, temperature and relative humidity fields were more uniform in the cold storage which had the uniform air supply duct.

Keywords:flow fields; cold storage; uniform air supply duct; air distribution; non-uniformity coefficient

通信作者：※南曉紅，女，陜西西安人，教授，主要從事制冷技術(shù)領(lǐng)域的科研與教學(xué)工作研究。西安西安建筑科技大學(xué)環(huán)境與市政工程學(xué)院，710055。Email: nanxh@xauat.edu.cn

作者簡介：劉曉菲，男，河南平頂山人，研究方向?yàn)槔鋷靸?nèi)溫濕度場的數(shù)值模擬。西安西安建筑科技大學(xué)環(huán)境與市政工程學(xué)院，710055。Email:1126558853@qq.com

基金項(xiàng)目：“十一五”國家科技支撐計(jì)劃（2008BAJ08B07-1）

收稿日期：2015-08-09

修訂日期：2015-11-12

中圖分類號：TB61+1

文獻(xiàn)標(biāo)志碼：A

文章編號：1002-6819（2016）-01-0091-06

doi：10.11975/j.issn.1002-6819.2016.01.012