基于CFD技術(shù)的預(yù)冷參數(shù)對(duì)盒裝荔枝微環(huán)境溫濕度的影響

郭嘉明，蔣易宏，林濟(jì)誠(chéng)，張霄丹，蔡威，林國(guó)鵬，劉東峰，曾志雄，呂恩利

基于CFD技術(shù)的預(yù)冷參數(shù)對(duì)盒裝荔枝微環(huán)境溫濕度的影響

郭嘉明，蔣易宏，林濟(jì)誠(chéng)，張霄丹，蔡威，林國(guó)鵬，劉東峰，曾志雄，呂恩利*

（華南農(nóng)業(yè)大學(xué) 工程學(xué)院，廣州 510642）

深入研究外環(huán)境參數(shù)對(duì)荔枝包裝盒內(nèi)環(huán)境的影響，獲得荔枝包裝盒內(nèi)外環(huán)境溫度濕度規(guī)律。建立盒裝荔枝果實(shí)和外流場(chǎng)的三維模型。結(jié)合荔枝果實(shí)和包裝物理特性，對(duì)預(yù)冷過(guò)程中包裝荔枝果實(shí)內(nèi)外環(huán)境溫濕度變化進(jìn)行數(shù)值分析，獲得包裝內(nèi)流場(chǎng)的分布情況。在達(dá)到二分之一冷卻時(shí)間之前，提高風(fēng)速可以縮短預(yù)冷時(shí)間，并且能夠加快荔枝果實(shí)的傳熱；在達(dá)到二分之一冷卻時(shí)間之后，風(fēng)速超過(guò)6 m/s后，荔枝果實(shí)溫度變化不顯著，相鄰風(fēng)速間果實(shí)溫度標(biāo)準(zhǔn)差的差值較小，標(biāo)準(zhǔn)差趨于穩(wěn)定；包裝內(nèi)相對(duì)濕度主要受荔枝果實(shí)溫度的影響；預(yù)冷初始階段選擇較低環(huán)境溫度，有利于加快荔枝果實(shí)與空氣的傳熱，促進(jìn)果實(shí)溫度的降低，但是荔枝果實(shí)間的溫度差異隨著環(huán)境溫度降低而增大。在預(yù)冷期間，外環(huán)境參數(shù)溫度和風(fēng)速對(duì)包裝內(nèi)環(huán)境溫濕度有較大的影響。

荔枝；預(yù)冷；數(shù)值模擬

荔枝生產(chǎn)于高溫高濕的夏季，采后易腐爛褐變[1]，不耐貯運(yùn)，通過(guò)預(yù)冷可以有效延長(zhǎng)荔枝保鮮周期、減小荔枝質(zhì)量損失[2-3]。冷庫(kù)預(yù)冷是常見(jiàn)預(yù)冷方式之一，其原理是冷風(fēng)經(jīng)過(guò)果蔬時(shí)，與果蔬進(jìn)行對(duì)流換熱，從而降低果蔬的溫度[4]。冷庫(kù)預(yù)冷通過(guò)調(diào)節(jié)風(fēng)速、溫度、相對(duì)濕度等多目標(biāo)參數(shù)，可延長(zhǎng)荔枝保鮮周期。然而，環(huán)境參數(shù)波動(dòng)或制冷風(fēng)機(jī)高風(fēng)速易導(dǎo)致荔枝失水褐變加速[5]。包裝可以減小環(huán)境參數(shù)波動(dòng)對(duì)荔枝品質(zhì)的影響，但會(huì)阻礙內(nèi)外環(huán)境的熱質(zhì)傳遞，導(dǎo)致內(nèi)外參數(shù)存在差異。因此，對(duì)帶包裝荔枝果實(shí)內(nèi)外環(huán)境參數(shù)變化進(jìn)行研究，有利于降低預(yù)冷能耗、提高預(yù)冷效果。通過(guò)試驗(yàn)可以研究預(yù)冷過(guò)程中包裝荔枝果實(shí)內(nèi)外環(huán)境參數(shù)變化，但傳感器很難準(zhǔn)確地顯示出包裝內(nèi)各參數(shù)的平均值，包裝內(nèi)流場(chǎng)變化復(fù)雜，并且多目標(biāo)參數(shù)的調(diào)控需要花費(fèi)大量的成本和人力。

計(jì)算流體力學(xué)（Computational Fluid Dynamics, CFD）是通過(guò)計(jì)算機(jī)模擬流體流動(dòng)、傳熱、化學(xué)反應(yīng)等物理現(xiàn)象的技術(shù)。近年來(lái)，CFD已在預(yù)冷研究領(lǐng)域得到廣泛運(yùn)用[6-8]，其卓越的可視化能力和數(shù)值模擬的高精度性質(zhì)備受認(rèn)可[9]，能夠有效地深入分析氣流特性和溫度分布，為研究對(duì)象提供更為全面的解析。郭嘉明等[10]建立了荔枝果實(shí)預(yù)冷二維數(shù)值模型，對(duì)不同風(fēng)速和不同溫度的預(yù)冷過(guò)程進(jìn)行了數(shù)值模擬，研究了荔枝果實(shí)的降溫特性，但沒(méi)對(duì)包裝荔枝內(nèi)外環(huán)境參數(shù)對(duì)荔枝降溫效果的影響進(jìn)行研究，并且三維數(shù)值模型具有與實(shí)際更接近、邊界條件簡(jiǎn)便等優(yōu)點(diǎn)，有利于開(kāi)展不同參數(shù)對(duì)比研究[11]。本文以PET盒包裝荔枝果實(shí)預(yù)冷過(guò)程為研究對(duì)象，建立PET盒裝荔枝果實(shí)和外流場(chǎng)的三維模型，運(yùn)用ICEM進(jìn)行非結(jié)構(gòu)化網(wǎng)格劃分，并采用FLUENT對(duì)模型進(jìn)行數(shù)值模擬，獲得PET盒內(nèi)外環(huán)境溫濕度規(guī)律，為荔枝果實(shí)預(yù)冷過(guò)程優(yōu)化提供一定的參考。

1 材料與方法

本文采用已有的保鮮試驗(yàn)平臺(tái)，廂體基于壓差原理進(jìn)行運(yùn)作。該廂體的外尺寸：長(zhǎng)為1.9 m、寬為1.1 m、高為1.5 m，如圖1所示，整體被隔熱材料聚氨酯板所覆蓋。該廂體被劃分為保鮮區(qū)和制冷區(qū)，位于風(fēng)機(jī)后方的回風(fēng)道產(chǎn)生負(fù)壓，使保鮮區(qū)內(nèi)的氣體被推入制冷區(qū)，在經(jīng)過(guò)蒸發(fā)器的冷卻后，冷卻的氣體通過(guò)開(kāi)孔隔板重新流回至保鮮區(qū)，從而形成一個(gè)封閉式循環(huán)。通風(fēng)風(fēng)速可以通過(guò)調(diào)整回風(fēng)道的風(fēng)速進(jìn)行改變，從而實(shí)現(xiàn)通風(fēng)效果的調(diào)控。

1.風(fēng)機(jī)；2.蒸發(fā)器；3.壓力區(qū)；4.開(kāi)孔隔板； 5.回風(fēng)道；6.架子；7.盒裝荔枝；8.保鮮區(qū)。

本文主要研究塑料盒包裝的內(nèi)外環(huán)境溫差特性，通過(guò)CATIA建立了實(shí)驗(yàn)配置的計(jì)算模型，精確模擬整個(gè)試驗(yàn)平臺(tái)需要大量的計(jì)算和時(shí)間成本，因此對(duì)模型進(jìn)行簡(jiǎn)化。將荔枝堆放在塑料盒內(nèi)，單個(gè)產(chǎn)品被離散地建模。塑料盒具體結(jié)構(gòu)及外觀如圖2所示，材質(zhì)為PET塑料，開(kāi)孔率為7.5%。其中，換氣口A、B尺寸（長(zhǎng)×寬）為26 mm×7 mm，換氣口C、D尺寸（長(zhǎng)×寬）為48 mm×7 mm，換氣孔直徑為8 mm。荔枝果實(shí)被建模為直徑300 mm的球體和各向同性的連續(xù)介質(zhì)。共23顆荔枝，并且荔枝果實(shí)被擺放為雙層，默認(rèn)它們大小均勻一致，理想地布置在塑料盒內(nèi)。本文研究了通風(fēng)風(fēng)速、預(yù)冷溫度對(duì)盒裝荔枝內(nèi)外環(huán)境溫差的影響。

圖2 三維計(jì)算模型

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

塑料箱內(nèi)荔枝果實(shí)的不穩(wěn)定傳熱和傳質(zhì)在強(qiáng)制空氣冷卻過(guò)程中非常復(fù)雜，不僅包括以對(duì)流為主導(dǎo)的顯熱交換，還包括通過(guò)熱傳導(dǎo)、輻射、蒸騰和呼吸等進(jìn)行的潛熱傳遞。因此，為了簡(jiǎn)化計(jì)算且仍然能正確描述實(shí)驗(yàn)系統(tǒng)，做出了以下假設(shè)[12-13]：

1）忽略荔枝果皮與果肉、果肉與果核之間傳熱差異的影響，即整個(gè)果實(shí)的熱物理參數(shù)是統(tǒng)一的，并且認(rèn)為成熟或衰老不會(huì)影響荔枝果實(shí)的呼吸和蒸騰作用。

2）忽略輻射傳熱。

3）箱體內(nèi)空氣為牛頓不可壓縮流體，符合Boussinesq假設(shè)。

4）各物質(zhì)的熱物理性質(zhì)參數(shù)在模擬過(guò)程中是恒定的。

在預(yù)冷過(guò)程中，每個(gè)域的流動(dòng)和傳熱特性由它們的控制方程描述，如下所示[14]。

質(zhì)量守恒方程：

式中：f為流體密度，g/m3；為速度矢量，m/s；m為質(zhì)量源項(xiàng)，kg/(m3·s)。

動(dòng)量守恒方程：

能量方程[10]：

荔枝呼吸熱模型[10]：

式中：s為呼吸熱，W/m3；p為荔枝果實(shí)平均溫度，K；為荔枝密度，kg/m3。

3 CFD模型及參數(shù)

3.1 網(wǎng)格獨(dú)立性研究

文中運(yùn)用ICEM對(duì)模型結(jié)構(gòu)進(jìn)行網(wǎng)格劃分，將計(jì)算域劃分為四面體非結(jié)構(gòu)化網(wǎng)格，如圖3所示。采取4種不同數(shù)量的網(wǎng)格，用于本節(jié)的網(wǎng)格獨(dú)立性的研究[15]，在保證數(shù)值計(jì)算結(jié)果滿足精度的前提下，減少網(wǎng)格數(shù)量能有效縮短計(jì)算時(shí)間。表1顯示了4種網(wǎng)格數(shù)量，網(wǎng)格尺寸由粗糙到精細(xì)，相鄰的2套網(wǎng)格數(shù)量大概相差1倍。在相同速度（=4 m/s）的情況下，采用瞬態(tài)求解器對(duì)模型進(jìn)行求解，到1 800 s時(shí)停止迭代，得到荔枝果實(shí)的平均溫度。表1列出了4套網(wǎng)格的果實(shí)平均溫度以及相對(duì)誤差。

圖3 網(wǎng)格模型結(jié)構(gòu)

表1 不同網(wǎng)格數(shù)量的對(duì)比

Tab.1 Comparison of the number of different meshes

從表1中可以明顯看出Ⅲ號(hào)網(wǎng)格與Ⅳ號(hào)網(wǎng)格之間的相對(duì)誤差明顯小于Ⅰ號(hào)網(wǎng)格、Ⅱ號(hào)網(wǎng)格與Ⅳ號(hào)網(wǎng)格之間的相對(duì)誤差?？紤]到計(jì)算時(shí)間問(wèn)題以及保證計(jì)算結(jié)果精度，本文數(shù)值模擬選擇Ⅲ號(hào)網(wǎng)格，即網(wǎng)絡(luò)數(shù)量為286萬(wàn)，其相對(duì)誤差在可接受范圍內(nèi)。

3.2 邊界條件及初始條件設(shè)置

雷諾數(shù)作為仿真的重要邊界參數(shù)，并通過(guò)雷諾數(shù)判斷流體類型。對(duì)模型的雷諾數(shù)進(jìn)行求解。

式中：為雷諾數(shù)；為流體的密度（空氣的密度為1.225 kg/m3）；為特征長(zhǎng)度，m；為入口速度，m/s；為黏度系數(shù)，一般為1.789 4×10?5。

3.2.1 入口邊界條件

速度入口邊界條件用于定義入口速度。在恒定入口溫度（溫度為3 ℃）的模擬中，入口速度分別設(shè)置為2、4、6、8、10 m/s。在恒定入口速度（速度為4 m/s）的模擬中，入口溫度分別設(shè)置為?1、1、3、5 ℃。入口相對(duì)濕度恒為90%。水力直徑計(jì)算入口為長(zhǎng)方形，湍流強(qiáng)度計(jì)算見(jiàn)式（6）[16]。

3.2.2 出口邊界條件

模型出口為壓力出口，其中湍流強(qiáng)度和水力直徑設(shè)置與入口一致。

3.2.3 壁面條件

縱觀中國(guó)文學(xué)發(fā)展史可以得知：“五四”以前的文學(xué)作品中，農(nóng)民是沒(méi)有地位的，即便是出現(xiàn)過(guò)描寫農(nóng)民起義的《水滸傳》，但由于歷史和階級(jí)的局限性，作者沒(méi)有塑造出能完全反映農(nóng)民社會(huì)地位和精神面貌的人物形象，作品中雖然刻畫了農(nóng)民出身的宋江等系列人物，但他們并沒(méi)有脫離原始的宗法意識(shí)形態(tài)的藩籬。至于其他作品中歪曲，甚至丑化農(nóng)民形象的比比皆是。

將外流場(chǎng)壁面設(shè)為熱流量邊界條件，即將外界對(duì)保鮮試驗(yàn)平臺(tái)的影響通過(guò)設(shè)置壁面熱流條件添加進(jìn)模型。

將網(wǎng)格文件導(dǎo)入FLUENT軟件后，進(jìn)行穩(wěn)態(tài)求解器求解，待流場(chǎng)穩(wěn)定，即穩(wěn)態(tài)收斂后（殘差收斂標(biāo)準(zhǔn)為10?3），以穩(wěn)態(tài)收斂值作為瞬態(tài)求解的初始值。初始化環(huán)境、荔枝果實(shí)溫度為23 ℃。本研究中，采用SIMPLE算法來(lái)處理壓力-速度耦合問(wèn)題，對(duì)梯度進(jìn)行計(jì)算時(shí)運(yùn)用了最小二乘法。在計(jì)算壓力、動(dòng)量和能量時(shí)，采用了二階離散方案，而湍動(dòng)能和湍能耗散率的求解則采用了一階離散方案。在文中，將半冷卻時(shí)間和八分之七冷卻時(shí)間作為預(yù)冷結(jié)束的標(biāo)準(zhǔn)，并采用以下計(jì)算方法[17]：

模型中具體的物性參數(shù)如表2所示

表2 物性參數(shù)設(shè)置

Tab.2 Setting of physical property parameters

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

4.1 模型驗(yàn)證

為了驗(yàn)證簡(jiǎn)化后的模型，將模擬結(jié)果與試驗(yàn)測(cè)量值進(jìn)行對(duì)比。試驗(yàn)在保鮮試驗(yàn)臺(tái)上進(jìn)行，調(diào)節(jié)隔板孔出口處的空氣溫度為–1～2 ℃，調(diào)節(jié)回風(fēng)道的風(fēng)速為5 m/s。放入2個(gè)包裝盒進(jìn)行試驗(yàn)，并且每盒包裝布置2路PT100傳感器（德國(guó)賀利氏鉑電阻，量程為–50～200 ℃，精度為±0.1 ℃）測(cè)量包裝盒內(nèi)環(huán)境溫度。使用葉輪風(fēng)速儀（Testo410i，量程為0.4～30 m/s，精度為（±0.2+0.02）m/s，Testo SE & Co. KGaA）測(cè)出開(kāi)孔隔板后面的風(fēng)速，以模擬入口的風(fēng)速。制冷前，維持廂體環(huán)境溫度和荔枝果實(shí)溫度為（20±0.5）℃；關(guān)閉廂門后，開(kāi)啟制冷系統(tǒng)2 h后停止預(yù)冷。溫度模擬值和試驗(yàn)值如圖4所示，總體而言，試驗(yàn)值與模擬值差異較小。通過(guò)計(jì)算，模擬值和試驗(yàn)值的平均誤差為2.80%、均方根誤差為3.20%，證明所建模型擬合程度較優(yōu)。

圖4 模擬溫度和試驗(yàn)溫度對(duì)比

4.2 風(fēng)速對(duì)包裝內(nèi)環(huán)境的影響

圖5a給出了塑料板條箱內(nèi)的荔枝果實(shí)的平均中心溫度，顯示了不同風(fēng)速的相應(yīng)半冷卻時(shí)間和八分之七冷卻時(shí)間。荔枝果實(shí)的預(yù)冷時(shí)間隨風(fēng)速增大不斷縮短，但風(fēng)速超過(guò)6 m/s后只有相對(duì)較小的差異，說(shuō)明在風(fēng)速達(dá)到臨界點(diǎn)時(shí)，風(fēng)速對(duì)預(yù)冷效果和預(yù)冷時(shí)間的影響減小。半冷卻時(shí)間之前，除了2 m/s風(fēng)速外，其他風(fēng)速下的荔枝果實(shí)的溫度在此階段快速下降。這說(shuō)明超過(guò)一定風(fēng)速后對(duì)半冷卻時(shí)間之前的冷卻速率沒(méi)有顯著影響，這可能是因?yàn)樵诖穗A段，荔枝果實(shí)和冷空氣之間溫差較大（從圖5a可以看出），在較低風(fēng)速下也會(huì)進(jìn)行大量的對(duì)流換熱。綜上，在半冷卻時(shí)間之前，采用較低風(fēng)速也可以滿足現(xiàn)階段產(chǎn)品冷卻的需求，并且還可以節(jié)省能源。

圖5b給出了不同風(fēng)速下PET包裝盒內(nèi)相對(duì)濕度的變化。從圖5b中可以看出，風(fēng)速越大，相對(duì)濕度上升得越快，并且到八分之七預(yù)冷時(shí)間的相對(duì)濕度越高。在八分之七預(yù)冷時(shí)間時(shí)，最大風(fēng)速和最小風(fēng)速的平均相對(duì)濕度相差2%。對(duì)比圖5a和圖5b可以看出，在預(yù)冷過(guò)程中進(jìn)行加濕處理，包裝內(nèi)的平均相對(duì)濕度上升，與溫度下降趨勢(shì)相反。這說(shuō)明風(fēng)速對(duì)包裝內(nèi)相對(duì)濕度有一定的影響，但包裝內(nèi)荔枝果實(shí)溫度對(duì)相對(duì)濕度的影響較大。

圖6顯示了預(yù)冷30 min時(shí)PET塑料盒和荔枝果實(shí)的溫度分布。從圖6中可以看出，隨著風(fēng)速的增加，荔枝果實(shí)的冷卻效果越好（高溫區(qū)域減小），尤其是當(dāng)風(fēng)速增大到4 m/s時(shí)，荔枝果實(shí)高溫區(qū)域驟降。說(shuō)明當(dāng)風(fēng)速在2～4 m/s時(shí)，增加風(fēng)速能夠加快荔枝果實(shí)的冷卻，然而風(fēng)速進(jìn)一步超過(guò)6 m/s后，荔枝果實(shí)溫度變化不顯著。

從表3可以看出不同風(fēng)速下八分之七冷卻時(shí)間的溫度標(biāo)準(zhǔn)差，其中樣本為每個(gè)荔枝果實(shí)的中心溫度。標(biāo)準(zhǔn)差隨著風(fēng)速的增大在減少，說(shuō)明風(fēng)速高，有利于熱量交換，減小盒內(nèi)荔枝果實(shí)溫度差異性，提高冷負(fù)荷利用效率，每個(gè)荔枝都能有較好的降溫效果。然而相鄰風(fēng)速間標(biāo)準(zhǔn)差的差值逐漸減小，6 m/s和8 m/s風(fēng)速間標(biāo)準(zhǔn)差的差值為0.099，小于0.1。說(shuō)明過(guò)度增加風(fēng)速不能有效地提高預(yù)冷效果，反而會(huì)增加能耗。綜上可得，6 m/s的風(fēng)速能夠有助于荔枝果實(shí)的傳熱與降溫。

4.3 溫度對(duì)包裝內(nèi)溫度的影響

圖7顯示了不同環(huán)境溫度下PET盒內(nèi)外環(huán)境溫差。從圖7中可以看出，4種環(huán)境溫度的溫差下降速度趨于一致，并且最終都趨向于0，說(shuō)明溫差下降速度不受環(huán)境溫度的影響。然而，隨著環(huán)境溫度的升高，開(kāi)始階段的溫差有明顯的區(qū)別，環(huán)境溫度越高，開(kāi)始階段溫差越小。這有可能是因?yàn)槔笾Τ跏紲囟认嗤?，且荔枝不斷產(chǎn)生呼吸熱，造成盒內(nèi)溫度較高，而環(huán)境溫度越高就更加接近包裝內(nèi)溫度。因此，在預(yù)冷開(kāi)始階段選擇較低環(huán)境溫度有利于提高PET盒內(nèi)外環(huán)境溫差，加快荔枝果實(shí)與空氣的傳熱，促進(jìn)果實(shí)溫度的降低。

環(huán)境溫度對(duì)果實(shí)中心溫度標(biāo)準(zhǔn)差的影響如表4所示。從表4中看出，隨著環(huán)境溫度的升高，荔枝果實(shí)中心溫度標(biāo)準(zhǔn)差在減小。說(shuō)明整體環(huán)境溫度越高，果實(shí)中心溫度的差異性就越小，并且環(huán)境溫度不是越低越好。

圖5 風(fēng)速對(duì)荔枝果實(shí)溫度、相對(duì)濕度、內(nèi)外環(huán)境溫差的影響

圖6 不同風(fēng)速下包裝內(nèi)的溫度分布

表3 風(fēng)速對(duì)荔枝中心溫度標(biāo)準(zhǔn)差的影響

Tab.3 Effect of wind speed on standard deviation of litchi center temperature

圖7 溫度對(duì)包裝內(nèi)外環(huán)境溫差的影響

表4 溫度對(duì)荔枝中心溫度標(biāo)準(zhǔn)差的影響

Tab.4 Effect of temperature on standard deviation of litchi center temperature

5 結(jié)語(yǔ)

以PET盒包裝荔枝果實(shí)為研究對(duì)象，建立三維數(shù)值模型進(jìn)行數(shù)值模擬，通過(guò)改變通風(fēng)風(fēng)速和通風(fēng)溫度，研究包裝盒內(nèi)外環(huán)境溫差對(duì)荔枝果實(shí)降溫的影響。經(jīng)研究，得出以下結(jié)論：

1）在半冷卻時(shí)間之前，風(fēng)速對(duì)荔枝果實(shí)冷卻效果不顯著，可以使用較低風(fēng)速進(jìn)行預(yù)冷，以減少能耗。

2）在半冷卻時(shí)間之后，使用6 m/s的風(fēng)速進(jìn)行預(yù)冷，能夠快速對(duì)荔枝果實(shí)進(jìn)行降溫，減少預(yù)冷時(shí)間。

3）風(fēng)速增加，有助于熱量釋放并提高冷負(fù)荷的利用效率，并且有利于加濕。然而過(guò)度增加風(fēng)速不能有效地提高預(yù)冷效果，會(huì)增加能耗。

4）預(yù)冷開(kāi)始階段，選擇較低環(huán)境溫度能夠加快荔枝果實(shí)的降溫，但是會(huì)增大果實(shí)間的溫度差異。

[1] 魏華, 李小東. 可釋放SO2保鮮包裝袋在荔枝保鮮中的應(yīng)用[J]. 包裝工程, 2020, 41(3): 57-61.

WEI Hua, LI Xiao-dong. Application of Fresh-Keeping Packaging Bags with Releasable SO2in Litchi Preservation[J]. Packaging Engineering, 2020, 41(3): 57-61.

[2] HAN Qiang, GAO Hai-yan, CHEN Hang-jun, et al. Precooling and Ozone Treatments Affects Postharvest Quality of Black Mulberry (Morus Nigra) Fruits[J]. Food Chemistry, 2017, 221: 1947-1953.

[3] 甘玲, 陳明林, 呂恩利. 不同預(yù)冷方式對(duì)荔枝預(yù)冷效果的影響[J]. 現(xiàn)代農(nóng)業(yè)裝備, 2020, 41(5): 8-13.

GAN Ling, CHEN Ming-lin, LYU En-li. Effects of Different Precooling Methods on Precooling Effect of Litchi[J]. Modern Agricultural Equipment, 2020, 41(5): 8-13.

[4] 呂盛坪, 呂恩利, 陸華忠, 等. 不同預(yù)冷方式對(duì)荔枝貯藏品質(zhì)的影響[J]. 現(xiàn)代食品科技, 2014, 30(3): 157-162.

LYU Sheng-ping, LYU En-li, LU Hua-zhong, et al. Effect of Different Precooling Methods on the Storage Quality of Litchi[J]. Modern Food Science and Technology, 2014, 30(3): 157-162.

[5] 郭嘉明. 荔枝氣調(diào)保鮮熱質(zhì)耦合機(jī)理與優(yōu)化研究[D]. 廣州: 華南農(nóng)業(yè)大學(xué), 2016: 23-54.

GUO Jia-ming. Study on Coupling Mechanism and Optimization of Heat and Mass for Modified Atmosphere Preservation of Litchi[D]. Guangzhou: South China Agricultural University, 2016: 23-54.

[6] GUO Jia-ming, WEI Xin-yu, LI Bin, et al. Characteristic Analysis of Humidity Control in a Fresh-Keeping Container Using CFD Model[J]. Computers and Electronics in Agriculture, 2020, 179(4): 105816.

[7] WANG Da, LAI Yan-hua, JIA Bin-guang, et al. The Optimal Design and Energy Consumption Analysis of Forced Air Pre-Cooling Packaging System[J]. Applied Thermal Engineering, 2019, 165(2): 114592.

[8] GAO Yu-ping, SHAO Shuang-quan, TIAN Shen, et al. Energy Consumption Analysis of the Forced-Air Cooling Process with Alternating Ventilation Mode for Fresh Horticultural Produce[J]. Energy Procedia, 2017, 142: 2642-2647.

[9] AMBAW A, MUKAMA M, OPARA U L. Analysis of the Effects of Package Design on the Rate and Uniformity of Cooling of Stacked Pomegranates: Numerical and Experimental Studies[J]. Computers and Electronics in Agriculture, 2017, 136: 13-24.

[10] 郭嘉明, 呂恩利, 陸華忠, 等. 盒裝荔枝果實(shí)降溫特性數(shù)值分析與驗(yàn)證[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2016, 47(5): 218-224.

GUO Jia-ming, LYU En-li, LU Hua-zhong, et al. Numerical Analysis and Verification on Characteristics of Temperature Decreasing of Litchi Fruits with Packages[J]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(5): 218-224.

[11] HO S H, 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.

[12] 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.

[13] WANG X, FAN Z, LI B, et al. Variable air Supply Velocity of Forced-Air Precooling of Iceberg Lettuces: Optimal Cooling Strategies[J]. Applied Thermal Engineering, 2021, 187: 116484.

[14] GONG Y, CAO Y, ZHANG X. Forced-Air Precooling of Apples: Airflow Distribution and Precooling Effectiveness in Relation to the Gap Width between Tray Edge and Box Wall[J]. Postharvest Biology and Technology, 2021, 177: 111523.

[15] Guo JM, LV EL, LU HZ, et al., Numerical Simulation of Gas Exchange in Fresh-Keeping Transportation Containers with a Controlled Atmosphere[J]. Food Science and Technology Research, 2016. 22(4): 429-441.

[16] 朱紅鈞, 林元華, 謝龍漢. FLUENT流體分析及仿真實(shí)用教程[M]. 北京: 人民郵電出版社, 2010.

ZHU Hong-jun, LIN Yuan-hua, XIE Long-han. Practical Course of FLUENT Fluid Analysis and Simulation[M]. Beijing: Posts & Telecom Press, 2010.

[17] WANG Xi-fang, LI Bao-guo, FAN Zhong-yang. Packaging Optimization and Integral Evaluation of Cooling Efficiency of Apples[J]. Food Science and Technology Research, 2019, 25(5): 677-686.

Effect of CFD Technology Based Pre-cooling Parameters on Temperature and Humidity of Micro-environment for Litchi with Packages

GUO Jia-ming, JIANG Yi-hong, LIN Ji-cheng, ZHANG Xiao-dan, CAI Wei, LIN Guo-peng, LIU Dong-feng, ZENG Zhi-xiong, LYU En-li*

(College of Engineering, South China Agricultural University, Guangzhou 510642, China)

The work aims to study the effect of environmental parameters on the environment inside the litchi packaging box in depth to obtain the temperature and humidity pattern of the environment inside and outside the packages. A three-dimensional model of the litchi fruit in the box and the external flow field was established. Combining the physical characteristics of the litchi fruit and the package, the numerical analysis of the temperature changes inside and outside the packaged litchi fruit during the pre-cooling process was carried out, and the distribution of the flow field in the package was obtained. Before one-half the cooling time, increasing the wind speed shortened the pre-cooling time and accelerated the heat transfer to the litchi fruit. After reaching one-half of the cooling time and wind speeds above 6 m/s, the variation of litchi fruit temperature was not significant. The difference in standard deviation of fruit temperature between adjacent wind speeds was small and the standard deviation tended to stabilize. The relative humidity inside the package was mainly affected by the temperature of the litchi fruit. Selection of low ambient temperature in the initial stage of pre-cooling was conducive to accelerating the heat transfer between litchi fruit and air, and promoting the reduction of fruit temperature. But the temperature difference between litchi fruit increased as the ambient temperature decreased. During pre-cooling, the external ambient parameters temperature and wind speed have a large effect on the ambient temperature and humidity inside the package.

litchi; pre-cooling; numerical simulation

S667.1

A

1001-3563(2023)17-0059-07

10.19554/j.cnki.1001-3563.2023.17.008

2022-12-21

廣東省2019年省級(jí)農(nóng)業(yè)科技創(chuàng)新及推廣項(xiàng)目（2023KJ101）；農(nóng)產(chǎn)品保鮮物流共性關(guān)鍵技術(shù)研發(fā)創(chuàng)新團(tuán)隊(duì)項(xiàng)目（2023KJ145）；國(guó)家自然科學(xué)基金項(xiàng)目（31901736，31971806）；國(guó)家荔枝龍眼產(chǎn)業(yè)體系項(xiàng)目（CARS-32-11）

責(zé)任編輯：曾鈺嬋