果蔬氣調(diào)庫(kù)快速降氧時(shí)間的影響因素

周 博， 南曉紅

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

氣調(diào)貯藏是目前最先進(jìn)的果蔬保鮮貯藏方式，它在冷藏的基礎(chǔ)上增加氣體成分調(diào)節(jié)，采用專(zhuān)門(mén)的降氧設(shè)備，使庫(kù)內(nèi)的氧氣濃度迅速降低至規(guī)定的指標(biāo)，有效抑制果蔬的呼吸作用，達(dá)到盡量延長(zhǎng)果蔬貯藏時(shí)間的目的。為了使果蔬的貯藏環(huán)境盡快達(dá)到氣調(diào)指標(biāo)，作者對(duì)氣調(diào)庫(kù)快速降氧時(shí)間的影響因素進(jìn)行了研究。

影響氣調(diào)庫(kù)降氧時(shí)間的因素有：制氮機(jī)制氮能力和工作參數(shù)、貯藏果蔬的多少、果蔬的品種以及氣調(diào)庫(kù)的氣密性等。由于影響降氧時(shí)間的因素較多，而且降氧過(guò)程持續(xù)時(shí)間長(zhǎng)，如果采用實(shí)驗(yàn)研究將會(huì)消耗大量的財(cái)力、物力以及人力。計(jì)算流體力學(xué)（CFD）是基于計(jì)算機(jī)技術(shù)的一種數(shù)值計(jì)算工具，克服了傳統(tǒng)理論分析法在對(duì)象簡(jiǎn)化和計(jì)算求解方面的不足，突破了試驗(yàn)過(guò)程人力物力消耗以及試驗(yàn)周期長(zhǎng)等諸多限制[1]。因此，關(guān)于冷庫(kù)內(nèi)氣流組織的研究多采用CFD的方法。目前，國(guó)內(nèi)外利用CFD軟件關(guān)于冷庫(kù)內(nèi)氣流組織的研究多集中于普通冷庫(kù)內(nèi)氣流速度場(chǎng)和溫度場(chǎng)[2-11]，關(guān)于氣調(diào)庫(kù)的研究則較少。文獻(xiàn)[12-14]運(yùn)用CFD軟件對(duì)影響氣調(diào)庫(kù)內(nèi)氣流組織的多個(gè)設(shè)計(jì)參數(shù)進(jìn)行數(shù)值分析，為合理設(shè)計(jì)氣調(diào)庫(kù)內(nèi)的氣流組織和優(yōu)化氣調(diào)庫(kù)設(shè)計(jì)提供了參考依據(jù)。但是對(duì)氣調(diào)庫(kù)降氧過(guò)程的研究不夠深入，而且研究成果難以推廣到其它氣調(diào)庫(kù)。

作者利用CFD軟件建立了通用的氣調(diào)庫(kù)內(nèi)氣體組分濃度數(shù)值計(jì)算模型，該計(jì)算模型適用于求解氣調(diào)庫(kù)不同階段的氣體組分濃度分布和變化情況。為了使模型盡量與實(shí)際情況相一致，作者通過(guò)編譯User-Defined Function（UDF）程序來(lái)解決冷風(fēng)機(jī)送風(fēng)口邊界的各氣體組分質(zhì)量分?jǐn)?shù)和蘋(píng)果的呼吸強(qiáng)度隨庫(kù)內(nèi)氣體組分濃度的變化而改變這兩個(gè)問(wèn)題，并將計(jì)算結(jié)果與文獻(xiàn) [15]（其實(shí)驗(yàn)氣調(diào)庫(kù)與本研究對(duì)象為同一氣調(diào)庫(kù)）的實(shí)驗(yàn)結(jié)果進(jìn)行對(duì)比。

1 物理模型

主要研究氣調(diào)庫(kù)內(nèi)快速降氧過(guò)程，這一過(guò)程是通過(guò)充入高純度的氮?dú)庵脫Q庫(kù)內(nèi)氣體使氧氣濃度降至所需要的值來(lái)實(shí)現(xiàn)的。作者以西安某蘋(píng)果氣調(diào)庫(kù)為研究對(duì)象，氣調(diào)庫(kù)模型見(jiàn)圖1。庫(kù)體尺寸（長(zhǎng)寬高）為8.0 m×4.6 m×6.5 m，貨物堆放方式為一堆，尺寸（長(zhǎng)寬高）為6.0 m×3.6 m×5.5 m。該氣調(diào)庫(kù)氮?dú)獬渥⒐艿篮蛶?kù)內(nèi)氣體出口管道直徑均為0.15 m。

圖1 氣調(diào)庫(kù)物理模型Fig.1 Physicalmodel of CA room

2數(shù)學(xué)模型的建立

2.1 紊流模型的選擇

采用k-ε湍流模型和組分輸運(yùn)模型，該模型的控制方程包括連續(xù)性方程、動(dòng)量方程、能量方程、組分濃度方程、k方程及ε方程。模型控制方程的通用形式如下：

式中，Φ為通用變量，Γ為廣義擴(kuò)散系數(shù)，S為廣義源項(xiàng)。

2.2 氣體區(qū)數(shù)學(xué)模型的建立

進(jìn)行計(jì)算時(shí)，對(duì)氣調(diào)庫(kù)內(nèi)氣體區(qū)作出如下假設(shè)：氣調(diào)庫(kù)內(nèi)氣體為牛頓流體；氣體在流動(dòng)過(guò)程中是不可壓縮的；庫(kù)內(nèi)氣體不參與輻射；氣體物性參數(shù)為常數(shù)。

k-ε湍流模型求解氣體區(qū)控制方程用到的參數(shù)見(jiàn)表1。的建立

表1 k－ε模型中的參數(shù)Table 1 Coefficients of the k－ε model

2.3 貨物區(qū)數(shù)學(xué)模型的建立

氣調(diào)庫(kù)內(nèi)貨物堆放方式與普通冷藏庫(kù)不同，通常采用高堆滿(mǎn)放方式，這樣使庫(kù)內(nèi)氣體所占容積相對(duì)較小，有助于縮短降氧時(shí)間。針對(duì)貨物的特點(diǎn)，將蘋(píng)果視為球形固體顆粒，蘋(píng)果間存在間隙，因此將貨物區(qū)視為多孔介質(zhì)處理。

2.3.1 動(dòng)量方程 由于貨物箱壁的阻礙以及貨物區(qū)多孔結(jié)構(gòu)造成極大的流動(dòng)阻力，貨物區(qū)中的氣體流動(dòng)極其微弱，對(duì)氣調(diào)庫(kù)內(nèi)主流區(qū)的影響可以忽略，因此對(duì)貨物區(qū)中氣體流速進(jìn)行修正為u=v=w=0。

2.3.2 能量方程 本文研究的蘋(píng)果堆滿(mǎn)足文獻(xiàn)[17]中多孔介質(zhì)的假設(shè)，因此貨物區(qū)的能量傳遞控制方程為：

式中，Γ為熱擴(kuò)散系數(shù)，S為果蔬呼吸引起的源項(xiàng)，計(jì)算公式如下：

式（3）中：

式（3）～（5）中，下標(biāo) g 和 s分別代表氣體、固體；Qr為呼吸熱，W/kg，ξ為多孔介質(zhì)的孔隙率，可以取為 ξ=0.35[16]，λe為貨物區(qū)的當(dāng)量導(dǎo)熱系數(shù)。Qr數(shù)值是由呼吸強(qiáng)度的大小確定。呼吸強(qiáng)度是指單位質(zhì)量的果蔬在單位時(shí)間內(nèi)釋放出的二氧化碳摩爾數(shù)或者吸收的氧氣摩爾數(shù)，用WCO2或WO2表示，單位為mol/(kg·h)。在貯藏過(guò)程中，果蔬的呼吸強(qiáng)度是貯藏溫度T、貯藏環(huán)境中氧氣和二氧化碳組分濃度的函數(shù)。在本研究中，蘋(píng)果進(jìn)入氣調(diào)庫(kù)時(shí)已經(jīng)進(jìn)行預(yù)冷處理，溫度為蘋(píng)果適宜的貯藏溫度0℃。根據(jù)文獻(xiàn)[18]，經(jīng)過(guò)實(shí)驗(yàn)數(shù)據(jù)回歸后得到的蘋(píng)果呼吸強(qiáng)度計(jì)算公式（當(dāng)時(shí)t=0℃）如下：

式中η表示呼吸能量轉(zhuǎn)換成熱量的轉(zhuǎn)換系數(shù)，一般認(rèn)為呼吸能量轉(zhuǎn)換為熱量的部分約為55%[19]。

2.3.3 氣體組分濃度方程 貨物區(qū)的氣體組分濃度方程可以用以下控制方程描述：

根據(jù)果蔬有氧呼吸原理，則呼吸作用產(chǎn)生的呼吸熱為：

式中Γ為多孔介質(zhì)中的質(zhì)量擴(kuò)散系數(shù)，S為果蔬的呼吸作用所引起的源項(xiàng)。

果蔬呼吸作用消耗氧氣產(chǎn)生二氧化碳，這樣就導(dǎo)致了氧氣組分濃度方程和二氧化碳組分濃度方程源項(xiàng)的產(chǎn)生。根據(jù)文獻(xiàn)[18]，蘋(píng)果在健康呼吸條件下，吸收的氧氣摩爾數(shù)和釋放的二氧化碳摩爾數(shù)比值是十分接近于1的。因此，根據(jù)式子（6）確定氧氣組分濃度方程的源項(xiàng)SO2：

式子（9）中，ρS為蘋(píng)果密度，取 790 kg/m3；WO2為氧氣的摩爾質(zhì)量，為0.032 kg/mol，MO2為二氧化碳的摩爾質(zhì)量，為 0.044 kg/mol。

2.4 邊界條件的設(shè)定

將冷風(fēng)機(jī)送風(fēng)口和氮?dú)膺M(jìn)口均設(shè)置為速度入口，參數(shù)設(shè)置根據(jù)冷風(fēng)機(jī)型號(hào)和制氮機(jī)型號(hào)以及氮?dú)夤軓酱_定，出口邊界均設(shè)置為 outflow。壁面采用無(wú)滑移邊界，溫度按第三類(lèi)邊界條件處理，各氣體組分在壁面無(wú)通量。

2.5 UDF程序功能

編譯的UDF程序主要包含兩個(gè)方面：一是冷風(fēng)機(jī)送風(fēng)口邊界的各氣體組分的質(zhì)量分?jǐn)?shù)賦值，二是貨物區(qū)的氧氣、二氧化碳組分方程的源項(xiàng)以及能量方程的源項(xiàng)。

冷風(fēng)機(jī)送風(fēng)口邊界賦值UDF程序的功能為：每次計(jì)算時(shí)將冷風(fēng)機(jī)回風(fēng)口邊界各組分的平均質(zhì)量分?jǐn)?shù)計(jì)算后賦值給冷風(fēng)機(jī)送風(fēng)口邊界。這種處理方法是因?yàn)椋撼涞笛蹼A段氣調(diào)庫(kù)內(nèi)的氧氣、二氧化碳和氮?dú)饨M分的體積分?jǐn)?shù)是時(shí)刻變化的，而庫(kù)內(nèi)氣體通過(guò)冷風(fēng)機(jī)進(jìn)行強(qiáng)制對(duì)流時(shí)，經(jīng)過(guò)冷風(fēng)機(jī)的各組分質(zhì)量分?jǐn)?shù)并不會(huì)改變，因此該UDF程序的處理方法符合冷風(fēng)機(jī)處真實(shí)的傳質(zhì)情況。該UDF程序?yàn)檠h(huán)流動(dòng)問(wèn)題中入口邊界條件的確定提供了一種解決方法。

貨物區(qū)的源項(xiàng)UDF程序主要是描述蘋(píng)果呼吸強(qiáng)度的變化引起的貨物區(qū)氣體組分濃度方程和能量方程源項(xiàng)的變化。充氮降氧階段庫(kù)內(nèi)的O2和CO2體積分?jǐn)?shù)是時(shí)刻變化的，而蘋(píng)果的呼吸強(qiáng)度隨著周?chē)h(huán)境中O2和CO2體積分?jǐn)?shù)變化而改變，這樣就導(dǎo)致了蘋(píng)果吸收的O2量、釋放的CO2量以及產(chǎn)生的呼吸熱都在時(shí)刻變化。根據(jù)文獻(xiàn)[18]的蘋(píng)果呼吸強(qiáng)度擬合公式編寫(xiě)的貨物區(qū)源項(xiàng)UDF程序可以真實(shí)的描述蘋(píng)果在氣體調(diào)節(jié)過(guò)程中的生理活動(dòng)對(duì)周?chē)h(huán)境的影響。

3 結(jié)果與分析

3.1 充氮降氧過(guò)程氧氣體積分?jǐn)?shù)計(jì)算值與實(shí)驗(yàn)值對(duì)比

蘋(píng)果的最佳氣調(diào)貯藏參數(shù)為O2體積分?jǐn)?shù)3%，而文獻(xiàn)[20]指出庫(kù)內(nèi)O2體積分?jǐn)?shù)降至比規(guī)定值高出2～3個(gè)百分點(diǎn)時(shí)，即可停止降氧。作者對(duì)一實(shí)際氣調(diào)庫(kù)的降氧過(guò)程進(jìn)行了數(shù)值計(jì)算，其制氮機(jī)工作參數(shù)是制氮量為30 m3/h，氮?dú)饧兌葹?5%。庫(kù)內(nèi)氣體組分測(cè)點(diǎn)坐標(biāo)為X=3.2 m，Y=0.4 m，Z=5.0 m，為了進(jìn)行對(duì)比，在計(jì)算過(guò)程中在貨物區(qū)中心位置設(shè)置測(cè)點(diǎn)坐標(biāo)為X=2.3 m，Y=4.5 m，Z=2.85 m。為了驗(yàn)證氣調(diào)庫(kù)內(nèi)氣體流動(dòng)與傳熱傳質(zhì)數(shù)學(xué)模型及數(shù)值計(jì)算結(jié)果的正確性和可靠性，作者以相關(guān)文獻(xiàn)[15]的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行驗(yàn)證，得到了充氮降氧階段庫(kù)內(nèi)氣體組分濃度測(cè)點(diǎn)處的氧組分濃度計(jì)算值與實(shí)驗(yàn)值的對(duì)比，結(jié)果見(jiàn)圖2。計(jì)算值與實(shí)驗(yàn)值的相對(duì)偏差最大值為0.068，最小值為0.005 6，平均值為0.027，數(shù)值計(jì)算結(jié)果與實(shí)驗(yàn)值有較好的一致性。

3.2 充氮降氧過(guò)程結(jié)束時(shí)O2體積分?jǐn)?shù)分布規(guī)律

圖3為充氮降氧階段結(jié)束后，x=2.3 m截面O2體積分?jǐn)?shù)場(chǎng)分布。圖中可以看出氣體區(qū)的O2的組分由于風(fēng)機(jī)的強(qiáng)制對(duì)流作用使其分布比較均勻，貨物區(qū)O2體積分?jǐn)?shù)分布則呈現(xiàn)出明顯的梯度，這是由于貨物區(qū)的O2體積分?jǐn)?shù)變化主要是由于擴(kuò)散作用，根據(jù)文獻(xiàn)[21]高體積分?jǐn)?shù)側(cè)與低體積分?jǐn)?shù)側(cè)的組分在空間上的分布是均勻遞減的，從而導(dǎo)致了這種體積分?jǐn)?shù)分層現(xiàn)象。

圖2 氣調(diào)庫(kù)內(nèi)測(cè)點(diǎn)處氧氣體積分?jǐn)?shù)變化計(jì)算值與實(shí)驗(yàn)值的對(duì)比Fig.2 Comparison of simulated values and experimental values of oxygen specie concentration at measuring point

圖3 x=2.3 m截面O2組分濃度分布Fig.3 Distribution of O2specie concentration at the plane of x=2.3 m

通過(guò)上述分析可知：果蔬氣調(diào)庫(kù)內(nèi)氧氣體積分?jǐn)?shù)受制氮機(jī)工作參數(shù)和果蔬呼吸作用的綜合影響，并且?guī)靸?nèi)氧氣分布并不均勻，因此難以通過(guò)簡(jiǎn)單的數(shù)學(xué)計(jì)算獲得準(zhǔn)確的降氧時(shí)間，運(yùn)用CFD軟件建立三維數(shù)值計(jì)算模型進(jìn)行求解是必要且可行的。

3.3 氮?dú)饧兌葘?duì)降氧時(shí)間的影響

由于降低庫(kù)內(nèi)氧氣濃度是通過(guò)向庫(kù)內(nèi)充入高純度氮?dú)鈱?shí)現(xiàn)的，那么氮?dú)饧兌染蜁?huì)影響降氧速度。

對(duì)于任何已經(jīng)選定型號(hào)的制氮機(jī)或制氮系統(tǒng)，其流量和純度是呈反比的，即提高氮?dú)饧兌葎t氮?dú)饬髁繒?huì)下降。本研究中制氮機(jī)工作參數(shù)見(jiàn)表2。

表2 制氮機(jī)工作參數(shù)Table 2 Nitrogen machine operating parameters

圖4為不同氮?dú)饧兌认碌玫降臍怏w區(qū)測(cè)點(diǎn)位置的降氧曲線(xiàn)?？梢钥闯?，不同的氮?dú)饧兌葘?duì)應(yīng)著最快的降氧區(qū)間，當(dāng)氧氣體積分?jǐn)?shù)從21%降至8.8%之間任意一點(diǎn)時(shí)，95%純度的氮?dú)馑钑r(shí)間最短；而降至8.8%～6%之間任意一點(diǎn)時(shí)，96%純度的氮?dú)馑钑r(shí)間最短；而降至6%～4.4%之間任意一點(diǎn)時(shí)，97%純度的氮?dú)馑钑r(shí)間最短；而降至4.4%～2.9%之間任意一點(diǎn)時(shí)，98%純度的氮?dú)馑钑r(shí)間最短；而降至2.9%～1%之間任意一點(diǎn)時(shí)，99%純度的氮?dú)馑钑r(shí)間最短。

圖4 不同氮?dú)饧兌瘸渥⑶闆r的O2體積分?jǐn)?shù)變化Fig.4 Changing value of O2concentration of different purity nitrogen filled cases

3.4 不同種類(lèi)果蔬對(duì)降氧時(shí)間的影響

不同種類(lèi)或品種的果蔬呼吸強(qiáng)度相差很大，呼吸強(qiáng)度影響著果蔬的生理變化，同樣也影響著貯藏環(huán)境的變化。作者在已有的果蔬呼吸強(qiáng)度相關(guān)研究基礎(chǔ)上，選擇了與蘋(píng)果貯藏指標(biāo)接近但呼吸強(qiáng)度卻相差較大的黃桃作為對(duì)比，研究了蘋(píng)果和黃桃的呼吸作用對(duì)充氮降氧過(guò)程中降氧時(shí)間的影響，其中黃桃呼吸速率擬合公式由文獻(xiàn)[16]得到。

圖5為貯藏貨物分別為黃桃和蘋(píng)果時(shí)，庫(kù)內(nèi)氣體區(qū)測(cè)點(diǎn)位置的氧氣濃度變化對(duì)比。貨物堆體積相同，采用純度為96%、流量為26 m3/h的氮?dú)獬淙搿?/p>

根據(jù)相關(guān)研究的成果分別計(jì)算蘋(píng)果和黃桃的呼吸速率，得到溫度為0℃時(shí)，處于空氣中黃桃的呼吸速率為 11.8 mL/（kg·h）， 蘋(píng)果為 1.77 mL/（kg·h），3%氧氣、3%二氧化碳的氣體條件下黃桃的呼吸速率 6.4 mL/（kg·h），蘋(píng)果為 0.65 mL/（kg·h）。 根據(jù)計(jì)算結(jié)果可以看出，黃桃在同樣貯藏條件下呼吸速率是遠(yuǎn)大于蘋(píng)果的。從圖4可以看出，在同樣的氮?dú)饧兌群土髁肯拢捎邳S桃的呼吸強(qiáng)度比蘋(píng)果大，降氧時(shí)間有明顯縮短，蘋(píng)果需要15 h，黃桃只需要11.5 h。由于呼吸強(qiáng)度大的果蔬氧氣消耗速度和二氧化碳生成速度都較快，因此制氮設(shè)備和二氧化碳洗滌裝置進(jìn)行選型時(shí)，應(yīng)該根據(jù)果蔬呼吸強(qiáng)度的大小匹配合適的型號(hào)。

圖5 同樣條件下黃桃和蘋(píng)果的降氧時(shí)間對(duì)比Fig.5 Comparison of the reducing oxygen time of peaches and apples under the same conditions

4 結(jié)語(yǔ)

作者以西安某蘋(píng)果氣調(diào)庫(kù)為研究對(duì)象，建立了氣調(diào)庫(kù)內(nèi)氣體流動(dòng)、傳熱與傳質(zhì)的三維耦合數(shù)學(xué)求解模型，對(duì)氣調(diào)庫(kù)的快速降氧過(guò)程進(jìn)行了數(shù)值計(jì)算，并對(duì)數(shù)學(xué)模型及數(shù)值計(jì)算結(jié)果的正確性、可靠性進(jìn)行了驗(yàn)證。經(jīng)過(guò)驗(yàn)證，本研究得到的數(shù)值計(jì)算結(jié)果與實(shí)驗(yàn)值有較好的一致性，因此建立的數(shù)學(xué)模型可以較真實(shí)地反映氣調(diào)庫(kù)內(nèi)氧氣組分濃度的變化情況。

不同的氮?dú)饧兌葘?duì)應(yīng)著不同的最快降氧區(qū)間，對(duì)于蘋(píng)果氣調(diào)而言，制氮機(jī)的制氮純度為96%時(shí)可以在最短時(shí)間內(nèi)把氧氣體積分?jǐn)?shù)降至6%，不僅可以使果蔬快速進(jìn)入氣調(diào)效果，還可以節(jié)約降氧動(dòng)力的消耗。呼吸強(qiáng)度差別較大的果蔬降氧時(shí)間有明顯差異，根據(jù)果蔬呼吸強(qiáng)度的大小匹配合適型號(hào)的制氮設(shè)備和二氧化碳洗滌裝置不僅可以使果蔬達(dá)到更好的貯藏效果，還可以使設(shè)備投資更加合理。本研究對(duì)于獲取實(shí)際氣調(diào)庫(kù)降氧時(shí)間及選擇合理的制氮機(jī)、改進(jìn)氣調(diào)工藝具有重要參考價(jià)值。

[1]ZHAO Changqing，F(xiàn)U Zetian，LIU Xue，et al.Temperature monitoring and warming system in cool food chain transportation[J].Microcomputer Information，2010，26（6）：27-28.（in Chinese）

[2]DELELE M A，SCHENK A，RAMON H，et al.Evaluation of a chicory root cold store humidification system using computational fluid dynamics[J].Journal of Food Engineering，2009，94（1）：110-121.

[3]PAULL R.Effect of temperature and relative humidity on fresh commodity quality[J].Postharvest Biology and Technology，1999，15（3）：263-277.

[4]LOPEZ A，PIQUE M，ROMERO A，et al.Influence of cold-storage conditions on the quality of unshelled walnuts[J].International Journal of Refrigeration，1995，18（8）：544-549.

[5]TASHTOUSH B.Natural losses from vegetable and fruit products in cold storage[J].Food Control，2000，11（6）：465-470.

[6]DEVRES Y，BISHOP C.Computer model for weight loss and energy conservation in a fresh-produce refrigerated store[J].Applied Energy，1995，50（2）：97-117.

[7]CHOURASIA M，GOSWAMI T.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.

[8]CHOURASIA M，GOSWAMI T.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.

[9]JI Amin.Numerical simulation and experimental Verification of temperature Variation tluring controlled atomospherecold storage offruitsandVegetables[J].Transactions of the Chinese Society of Agricultural Engineering，2006，22（5）：24-27.（inChinese）

[10]WANG Yizhong，HU Chunyuan，CHEN Shaohui，et al.Study on numerical simulation of temperature and humidity fields in a refrigerated vehicle[J].Storage and Process，2010，10（3）：26-29.（in Chinese）

[11]LIU Bin，YANG Zhao，LI Xihong，et al.Experiment on air disrtibution and storage effects in mini-cold storage[J].Journal of Tianjin University，2005，38（10）：897-900.（in Chinese）

[12]李敏華.氣調(diào)庫(kù)制氮系統(tǒng)的改進(jìn)設(shè)計(jì)及其氣流組織的數(shù)值模擬[D].廣州：華南理工大學(xué)，2004.

[13]瞿曉華.果蔬氣調(diào)庫(kù)的氣體流場(chǎng)分析及優(yōu)化設(shè)計(jì)[D].上海：上海水產(chǎn)大學(xué)，2004.

[14]李勇祁.裝配式氣調(diào)庫(kù)濃度分布及CAD研究[D].上海：上海水產(chǎn)大學(xué)，2001.

[15]HU Hao，YUAN Xiuling，SUN Dawei，et al.Three dimensional dynamic simulation and experiment investigation about temperature and oxygen density variationprocess of the fruit CA storage[J].Journal of Refrigeration，2009，94：110-122.（in Chinese）

[16]劉穎.果蔬快速降氧氣調(diào)貯藏環(huán)境的熱動(dòng)力學(xué)研究[D].上海：上海交通大學(xué)，2004.

[17]林瑞泰.多孔介質(zhì)傳熱傳質(zhì)引論[M].北京：科學(xué)出版社，1995.

[18]LIANG Dawei，ZENG Danling，AO Yue，et al.The mathematical model of the storage of fruits and vegetables by adjusting the atmosphere inside the system[J].Journal of Engineering Thermophysics，1990，11（1）：5-8.（in Chinese）

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

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

[21]連之偉.熱質(zhì)交換原理與設(shè)備[M].北京：中國(guó)建筑工業(yè)出版社，2006.