不同CO2濃度下雞蛋呼吸仿真分析及驗(yàn)證

王嬌嬌，王巧華，曹 芮，謝娟娟

不同CO2濃度下雞蛋呼吸仿真分析及驗(yàn)證

王嬌嬌1，王巧華2，曹 芮1，謝娟娟1

（1. 周口師范學(xué)院機(jī)械與電氣工程學(xué)院，周口 466001；2. 華中農(nóng)業(yè)大學(xué)工學(xué)院，國家蛋品加工技術(shù)研發(fā)分中心，武漢 430070）

貯存環(huán)境中CO2濃度對雞蛋呼吸作用及新鮮度變化具有重要影響，為了進(jìn)一步考察在不同CO2濃度下雞蛋呼吸釋放CO2的擴(kuò)散及新鮮度變化情況，該研究測定貯存在溫度為25 ℃、相對濕度為65%、CO2體積分?jǐn)?shù)分別為1.5%、3.0%、4.5%及空氣（對照組）環(huán)境下雞蛋呼吸和新鮮度，利用FLUENT軟件完成不同CO2濃度下雞蛋第1天呼吸釋放CO2的擴(kuò)散過程仿真。結(jié)果表明，利用Fluent軟件計算所得氣體速度值和試驗(yàn)值基本一致，模擬值與實(shí)測值相對誤差在4%～9%。不同CO2濃度下雞蛋釋放CO2的擴(kuò)散過程符合重氣擴(kuò)散的特點(diǎn)，其擴(kuò)散方向受到CO2濃度的影響，且隨著貯存環(huán)境中CO2濃度的增加，雞蛋呼吸釋放CO2的擴(kuò)散量和擴(kuò)散速度逐漸變小，呼吸得到抑制。再通過分析貯存過程中不同CO2體積分?jǐn)?shù)下雞蛋呼吸強(qiáng)度和新鮮度變化得出，對照組1.5%CO2與3.0%CO2下的雞蛋呼吸強(qiáng)度和新鮮度具有顯著差異（<0.05），當(dāng)貯存環(huán)境中CO2體積分?jǐn)?shù)達(dá)到3.0%時，對雞蛋呼吸強(qiáng)度的抑制效果不再明顯，且新鮮度變化保持不變。3.0%CO2和4.5%CO2體積分?jǐn)?shù)下貯存雞蛋20 d，其新鮮度等級仍在AA級以上，綜合經(jīng)濟(jì)因素，3.0%CO2保鮮雞蛋效果較優(yōu)。該研究可為雞蛋呼吸和氣調(diào)保鮮技術(shù)提供理論及數(shù)據(jù)基礎(chǔ)。

貯藏；呼吸；雞蛋；CO2濃度；新鮮度；仿真

0 引 言

在貯存過程中，雞蛋通過蛋殼表層的氣孔與外界進(jìn)行CO2、水分等物質(zhì)交換，其中，氣體交換是由氣體濃度在蛋殼內(nèi)外分布的不均勻所產(chǎn)生的，其驅(qū)動力是濃度驅(qū)動，即氣體由高濃度向低濃度方向擴(kuò)散[1]。這種氣體交換指的就是雞蛋的呼吸作用，它是雞蛋的一種生理代謝活動，也是生命存在的重要標(biāo)志[2]。因此，貯存環(huán)境中的CO2濃度對雞蛋與外界CO2氣體交換具有重要影響。

CO2濃度對雞蛋的呼吸和新鮮度均有重要的影響。高濃度CO2貯存雞蛋，抑制雞蛋內(nèi)部CO2的排出，從而延緩貯存期間雞蛋哈夫單位的降低，該方法可以有效抑制雞蛋品質(zhì)劣變，延長雞蛋的保鮮期[3-5]。Rocculi等[6]通過使用空氣、100% N2、以及100% CO2包裝雞蛋進(jìn)行試驗(yàn)對比，結(jié)果得出用100% CO2包裝有助于提高雞蛋的蛋白品質(zhì)。袁曉龍等[7]通過采用不同濃度的CO2、N2、O2以及空氣的配比包裝雞蛋，結(jié)果發(fā)現(xiàn)配比為60%～100%的CO2包裝保鮮效果較好。劉美玉等[8]通過試驗(yàn)比較分析了在25 ℃和4 ℃貯存條件下不同氣體配比對雞蛋保鮮效果的影響，結(jié)果得出在25℃下，采用50%CO2、7%～11%O2以及39%～43%N2的氣體配比包裝雞蛋30 d，雞蛋的新鮮度等級仍保持AA級。黃群等[9]通過試驗(yàn)得出當(dāng)貯存環(huán)境中CO2體積分?jǐn)?shù)大于5%時，能顯著抑制S-卵白蛋白生成。馬逸霄等[10]對0.03%（空氣）、5%CO2、10%CO2及20%CO2貯存下雞蛋呼吸和品質(zhì)的變化情況，結(jié)果得出5%CO2保存雞蛋時，其呼吸被抑制，10%CO2時抑菌效果較優(yōu)。綜上可知貯存環(huán)境中CO2濃度是研究雞蛋呼吸強(qiáng)度和新鮮度變化的關(guān)鍵因素，且保鮮雞蛋所需的CO2濃度在逐漸降低。鑒于此，在實(shí)際生產(chǎn)過程中，若采用CO2氣體包裝，其體積分?jǐn)?shù)仍需進(jìn)一步精確調(diào)配。

雞蛋在貯藏過程中內(nèi)部的氣體會通過蛋殼擴(kuò)散到外界，其內(nèi)部發(fā)生的生理變化存在復(fù)雜性、多樣性，不能直觀觀察，因此雞蛋呼吸強(qiáng)度是其內(nèi)部變化的最直觀體現(xiàn)，也是影響貯存保鮮的關(guān)鍵因素[11]，同時有研究表明雞蛋呼吸強(qiáng)度與品質(zhì)之間具有密切聯(lián)系[12]，鑒于此，本文旨在研究雞蛋呼吸釋放CO2的量，探尋能夠抑制雞蛋呼吸實(shí)現(xiàn)品質(zhì)保鮮較佳的CO2氣體濃度。

FLUENT是目前流行的CFD仿真軟件，通常用于流體分析[13-15]。Perianu等[16]利用有限元方法模擬了蛋殼的動態(tài)特性和蛋殼內(nèi)流體的作用力。Sellés等[17]利用有限元法描述了幾種雞蛋類型的不同應(yīng)力模式，并分析了自由落體沖擊時雞蛋殼上的應(yīng)力。Fabbri等[18]利用有限元模型反演法確定CO2在雞蛋各組分中的擴(kuò)散量。因此，本文首先利用FLUENT軟件對貯存第1天的雞蛋呼吸釋放CO2的擴(kuò)散情況進(jìn)行了仿真模擬，重點(diǎn)分析了不同時間CO2的擴(kuò)散情況、不同CO2濃度下雞蛋釋放CO2擴(kuò)散情況。然后，根據(jù)不同CO2濃度下雞蛋呼吸強(qiáng)度與新鮮度的變化，綜合得出抑制雞蛋呼吸強(qiáng)度且實(shí)現(xiàn)保鮮的CO2濃度，以期為雞蛋貯存保鮮提供參考。

1 材料與方法

1.1 材料

試驗(yàn)樣本由武漢九峰雞場提供，選取新鮮海蘭粉殼雞蛋280枚，將其隨機(jī)分為4組，在溫度為25 ℃、相對濕度為65%、CO2體積分?jǐn)?shù)分別為1.5%、3.0%、4.5%及對照組（空氣）條件下80 L的培養(yǎng)箱中，貯存20d，每個貯存條件下取10枚雞蛋進(jìn)行編號，每2d測量其呼吸強(qiáng)度，并且每2d隨機(jī)取樣5枚測定雞蛋的哈夫單位。

1.2 儀器與設(shè)備

HWS-250型恒溫恒濕培養(yǎng)箱，武漢格雷思科技開發(fā)有限公司；CB210型二氧化碳培養(yǎng)箱，上海丙林電子科技有限公司；呼吸測定儀SY-1022，武漢鑫星星科學(xué)儀器有限公司；電子游標(biāo)卡尺，上海美耐特實(shí)業(yè)有限公司；SPX智能型生化培養(yǎng)箱，寧波江南儀器廠；JY3002稱重儀，上海精密科學(xué)儀器有限公司。

1.3 方法

1.3.1 雞蛋呼吸釋放CO2速度測定

雞蛋呼吸測定裝置如圖1所示。試驗(yàn)前將橡膠管連接呼吸室的一端，置于室外空曠處進(jìn)行校準(zhǔn)21 min。校準(zhǔn)完成后，將橡膠管重新連接到呼吸室，待數(shù)據(jù)穩(wěn)定后，將已編號的雞蛋稱質(zhì)量，放入1 L的呼吸室中密封，記錄此時的CO2初始濃度0，溫度0，運(yùn)行一段時間后，記錄CO2最終濃度，溫度T。

雞蛋的呼吸強(qiáng)度根據(jù)公式（1）計算得出[21]。

本試驗(yàn)采用的呼吸室大小1 L，雞蛋呼吸釋放CO2摩爾濃度根據(jù)公式（2）計算。

質(zhì)量流量根據(jù)式（3）計算。

式中為質(zhì)量流量，g/s；t為擴(kuò)散時間，s；為呼吸室體積，L。

速度根據(jù)式（4）計算。

1.3.2 雞蛋新鮮度測定

測量完雞蛋呼吸強(qiáng)度，用電子天平測雞蛋質(zhì)量，將雞蛋打破，用數(shù)顯游標(biāo)卡尺測雞蛋蛋白高度，測定3次取平均值，利用公式（5）計算樣本的哈夫單位（Hough Unit, HU），以其表征雞蛋新鮮度。

式中H為蛋白高度，mm。

1.3.3 雞蛋呼吸數(shù)值模擬

1）幾何模型構(gòu)建

雞蛋長軸尺寸范圍為53.29～60.42 mm，雞蛋短軸尺寸范圍為41.27～44.6 mm，蛋形指數(shù)在1.25～1.36之間。根據(jù)公式（6）繪制雞蛋模型二維圖，如圖2所示。

注：，分別為雞蛋的短軸和長軸；，分別為雞蛋邊緣曲線的橫坐標(biāo)軸和縱坐標(biāo)軸；10o為蛋形角。

Note:andare the short axis and long axis of the egg respectively;andare the abscissa axis and ordinate axis of the egg edge curve respectively; 10°is the egg shape angle.

圖2 雞蛋模型二維示意圖

Fig.2 Two dimensional diagram of egg model

2）幾何模型網(wǎng)格劃分

為了方便觀察雞蛋釋放CO2的過程，在雞蛋外部建立流體域，利用ICEM CFD軟件對模型的流體域進(jìn)行網(wǎng)格劃分，考慮到直觀觀察及便于觀察，在雞蛋外部設(shè)置為受限空間，其為126.43 mm×108.66 mm×108.66 mm的長方體，可以全方位觀察雞蛋釋放CO2的濃度分布情況。對模型流體區(qū)域進(jìn)行六面體結(jié)構(gòu)網(wǎng)格劃分，劃分后的網(wǎng)格內(nèi)沒有重疊，網(wǎng)格質(zhì)量為0.7，劃分后的網(wǎng)格數(shù)為431 884，符合計算要求，如圖3所示。

將劃分好的模型導(dǎo)入Fluent 15.0中，根據(jù)CFD理論方法設(shè)置相應(yīng)的參數(shù)以及邊界條件，對雞蛋呼吸釋放CO2的擴(kuò)散過程進(jìn)行數(shù)值模擬。

3）擴(kuò)散過程中模擬模型的基本假設(shè)條件

雞蛋內(nèi)部發(fā)生著復(fù)雜的化學(xué)變化，但是具體是何種化學(xué)變化是未知的，所以對其進(jìn)行完全的數(shù)學(xué)模擬是非常困難的，通過試驗(yàn)反復(fù)驗(yàn)證，CO2從雞蛋內(nèi)部擴(kuò)散到外部的速率在30 min內(nèi)是基本不變的，為了便于對CO2的擴(kuò)散過程進(jìn)行仿真，做如下基本假設(shè)[19]：①假設(shè)擴(kuò)散速度和面積始終保持某一值，不隨時間變化；②假設(shè)在CO2擴(kuò)散過程中，不發(fā)生化學(xué)、物理變化以及相變反應(yīng)，只發(fā)生擴(kuò)散；③假設(shè)CO2氣體與空氣完全作為理想氣體；④假設(shè)CO2的擴(kuò)散方向不隨時間、地點(diǎn)和高度變化；⑤ 假設(shè)在擴(kuò)散過程中，環(huán)境溫度在整個擴(kuò)散過程中不發(fā)生變化；⑥假設(shè)CO2在擴(kuò)散的過程中，外界環(huán)境參數(shù)改變時，不影響重力加速度，重力加速度始終為某一值。

4）擴(kuò)散過程中的基本方程

雞蛋內(nèi)部CO2向外界擴(kuò)散的基本方程除了連續(xù)方程、動量方程以及能量守恒方程外還有組分方程，其具體擴(kuò)散基本方程可表述如下。

連續(xù)方程

動量守恒方程

能量守恒方程

組分輸送方程

5）湍流模型的選擇

目前，在自然環(huán)境、工程及試驗(yàn)裝置中的流動一般為湍流流動，其數(shù)值模擬主要有直接數(shù)值模擬DNS（Direct Numerical Simulation）、雷諾平均方法RANS（Reynolds Average Navier-Stokes）以及大渦模擬LES（Large Eddy Simulation）三類[20]。其中，DNS在從理論分析上較另外兩種模擬更為準(zhǔn)確，但是對計算機(jī)的硬件要求較高，所以適用于模擬低雷諾數(shù)流動的情況。LES是利用濾波函數(shù)將大小尺度渦分離，對大尺度渦和小尺度渦進(jìn)行不同處理。大渦模擬從N-S（Navier-Stokes）方程出發(fā)，能夠模擬湍流過程中的一些細(xì)節(jié)，然而計算工作量仍然很大，難以適用于情況和條件發(fā)生復(fù)雜變化的流動情況。RANS是流體計算較為常用的方法，其理論依據(jù)為氣體流動理論，結(jié)合試驗(yàn)數(shù)據(jù)進(jìn)行數(shù)值模擬進(jìn)而得到模擬結(jié)果。

Fluent 軟件中提供的RANS模型包括-模型、-模型和雷諾應(yīng)力模型。選取湍流模型時一般考慮計算機(jī)能力、精度要求、流體的可壓性、流體域的運(yùn)動狀態(tài)以及能否滿足特定的求解問題等，本文選用標(biāo)準(zhǔn)-模型。

6）邊界條件的設(shè)定

本試驗(yàn)?zāi)M的是三維受限空間內(nèi)雞蛋釋放CO2的擴(kuò)散過程，所以選擇空間為三維空間模型。本試驗(yàn)采用單精度求解器、非耦合隱式求解法以及非穩(wěn)態(tài)計算模式，其中非穩(wěn)態(tài)計算方法選用一階隱式計算法。采用單相流模型，允許有少量的熱交換過程，所以設(shè)置開啟能量方程選項(xiàng)。在整個擴(kuò)散過程中，由于CO2密度比空氣大，所以選擇重力參數(shù)，考慮到雞蛋在貯存過程中，均為尖端向下，所以根據(jù)三維模型結(jié)構(gòu)，在方向上設(shè)置一個9.8 m/s2的重力加速度。另外考慮到氣體浮力影響，設(shè)置全浮力參數(shù)打開。

本試驗(yàn)采用初始化流場的邊界條件，從雞蛋表面進(jìn)行計算，設(shè)置擴(kuò)散開始時的速度為0.08 m/s，溫度為298.15 K，每10 s自動保存一次數(shù)據(jù)文件，迭代時間步長設(shè)為1，迭代總步數(shù)為1 800次，實(shí)現(xiàn)連續(xù)計算30 min內(nèi)的CO2的瞬態(tài)擴(kuò)散過程。

1.3.4 數(shù)據(jù)統(tǒng)計分析方法

應(yīng)用Excel 2007和SPSS Statistics 19.0軟件對試驗(yàn)結(jié)果進(jìn)行統(tǒng)計分析。通過Duncan多重比較法進(jìn)行差異顯著性檢驗(yàn)；＜0.05表示差異性顯著。

2 結(jié)果與分析

2.1 不同CO2環(huán)境下雞蛋呼吸釋放CO2擴(kuò)散分析

在貯存過程中，雞蛋新鮮度逐漸下降，內(nèi)部生理活動逐漸減弱，釋放CO2的量逐漸降低，由此可知貯存第1天的雞蛋最為新鮮，呼吸釋放CO2的量最大，后面逐漸降低[21]，因此本文模擬貯存在CO2濃度分別為空氣狀態(tài)、1.5%、3.0%、4.5%下雞蛋第1天呼吸釋放CO2的擴(kuò)散情況，其擴(kuò)散過程仿真如圖4所示。

由于CO2的密度比空氣大，在擴(kuò)散過程中，應(yīng)考慮到重力的影響，由圖4可知，雞蛋釋放CO2的擴(kuò)散過程符合重氣擴(kuò)散的特點(diǎn)，根據(jù)重氣擴(kuò)散的特征，可以把雞蛋呼吸釋放CO2的擴(kuò)散過程大致分為4個階段：初始階段，下降階段，混合階段以及湍流階段[22]。

2.1.1 初始階段

CO2濃度為空氣狀態(tài)時，雞蛋30 min內(nèi)釋放CO2的量為0.03 cm3/m3；CO2體積分?jǐn)?shù)為1.5%時，雞蛋30 min內(nèi)釋放CO2的量為0.01 cm3/m3；CO2體積分?jǐn)?shù)為3.0%時，雞蛋30 min內(nèi)釋放CO2的量為0.002 cm3/m3；CO2體積分?jǐn)?shù)為4.5%時，雞蛋30 min內(nèi)釋放CO2的量為0.001 9 cm3/m3。當(dāng)擴(kuò)散時間為60 s時，CO2擴(kuò)散屬于初始階段，附著在雞蛋表層，4個CO2濃度場的CO2分布沒有明顯的區(qū)別；當(dāng)擴(kuò)散時間為300 s時，貯存在CO2濃度為空氣狀態(tài)、1.5%CO2及3.0%CO2雞蛋釋放的CO2層明顯增厚，由此表明雞蛋上方氣云的厚度隨擴(kuò)散速度的增大而增大，這一現(xiàn)象與胡世明[23]所分析的釋放速度對重氣擴(kuò)散影響的相關(guān)研究一致。

2.1.2 下降階段

當(dāng)擴(kuò)散時間為600 s時，由于CO2重力的影響，4個CO2濃度場下雞蛋釋放的CO2分布均具有向下方沉積的趨勢，準(zhǔn)備進(jìn)入下降階段。當(dāng)擴(kuò)散時間為1 200 s時，4個CO2濃度場下雞蛋釋放的CO2分布均已進(jìn)入下降階段。CO2體積分?jǐn)?shù)為3.0%和4.5%中，CO2層分布均完全向下沉積，然而CO2濃度為空氣狀態(tài)和1.5%中，CO2層沉積并不完全，仍然具有上下浮動的趨勢，且浮動較為明顯。由此表明擴(kuò)散量越大，擴(kuò)散源周圍的氣體濃度越大，這一現(xiàn)象與司凡[24]所研究的天然氣氣體量大，泄露擴(kuò)散濃度大一致。

2.1.3 混合階段

當(dāng)擴(kuò)散時間為1 400 s時，CO2濃度為空氣狀態(tài)、3.0%CO2和4.5%CO2中，雞蛋釋放CO2的擴(kuò)散為逐漸向下沉積過程，為混合階段，然而CO2體積分?jǐn)?shù)為1.5%中，雞蛋釋放CO2的擴(kuò)散分布仍為下降階段，正在逐漸向下沉積，逐漸進(jìn)入混合階段。

2.1.4 湍流階段

當(dāng)擴(kuò)散時間為1 800 s時，CO2體積分?jǐn)?shù)為3.0%和4.5%中，雞蛋呼吸釋放CO2的擴(kuò)散分布已經(jīng)進(jìn)入湍流階段，CO2云團(tuán)隨著空氣流動的方向進(jìn)行擴(kuò)散，進(jìn)而混合達(dá)到完全擴(kuò)散，無明顯區(qū)別。CO2體積分?jǐn)?shù)為1.5%中，雞蛋釋放CO2的擴(kuò)散分布已進(jìn)入湍流階段，但擴(kuò)散明顯慢于另外2個CO2濃度?？諝鉅顟B(tài)下的雞蛋呼吸釋放CO2的擴(kuò)散分布已經(jīng)進(jìn)入湍流階段，附著在蛋殼表層的云團(tuán)明顯厚于其它3個濃度場的。結(jié)果表明擴(kuò)散量和擴(kuò)散速度對擴(kuò)散方向具有一定的影響。

2.2 模擬結(jié)果驗(yàn)證

上述對數(shù)值模擬所獲得的結(jié)果進(jìn)行了定性的分析與比較，本節(jié)通過試驗(yàn)驗(yàn)證模型的準(zhǔn)確性。試驗(yàn)測定貯存在CO2濃度分別為空氣狀態(tài)、1.5%、3.0%、4.5%下雞蛋第1天呼吸釋放CO2的速度，與利用仿真計算所得速度值進(jìn)行比較分析，如表1所示。

表1 氣體擴(kuò)散試驗(yàn)數(shù)據(jù)與仿真對照

利用Fluent軟件計算所得氣體速度值和試驗(yàn)值基本一致，模擬值與實(shí)測值相對誤差在4%～9%之間，說明-模型較適合雞蛋呼吸釋放CO2擴(kuò)散過程分析。

2.3 不同CO2濃度下雞蛋呼吸強(qiáng)度與新鮮度的變化

不同CO2濃度下雞蛋呼吸強(qiáng)度變化如表2所示，在25 ℃下，CO2體積分?jǐn)?shù)為1.5%、3.0%及4.5%時，雞蛋的呼吸強(qiáng)度均在1 mg/(kg·h)以下，當(dāng)CO2體積分?jǐn)?shù)為1.5%時，對雞蛋呼吸強(qiáng)度有一定的抑制作用，而CO2濃度為3.0%時抑制效果更為顯著（<0.05），此時繼續(xù)增加CO2體積分?jǐn)?shù)為4.0%時，抑制效果不再明顯改善（>0.05）。

根據(jù)表2分析可知，對照組（空氣）下貯存的雞蛋呼吸強(qiáng)度和新鮮度隨著貯存時間的延長而下降，且下降速度較快。隨著貯存環(huán)境中CO2體積分?jǐn)?shù)的增加，雞蛋呼吸強(qiáng)度逐漸減弱，新鮮度下降逐漸變緩，當(dāng)貯存環(huán)境中CO2體積分?jǐn)?shù)增至3.0%再增加CO2濃度，雞蛋新鮮度保持不變（>0.05）。綜合結(jié)果表明在CO2體積分?jǐn)?shù)為3.0%和4.5%的條件下，雞蛋呼吸強(qiáng)度均較小且變化情況較一致，雞蛋的新鮮度變化情況趨于一致，貯存20 d后，雞蛋的新鮮度均在新鮮等級AA以上。因此，綜合經(jīng)濟(jì)因素，3.0%CO2保鮮雞蛋效果較優(yōu)。

呼吸作用是雞蛋采集后內(nèi)部生理活動的外在表征，也是影響其貯運(yùn)效果的重要因素[25]。貯存環(huán)境中CO2的濃度是雞蛋呼吸及保鮮研究的熱點(diǎn)，那么雞蛋自身呼吸釋放CO2的量是相關(guān)研究的重要前提。根據(jù)檢測雞蛋呼吸釋放CO2量進(jìn)行貯存環(huán)境中CO2濃度的調(diào)控，有利于精確調(diào)配抑制雞蛋呼吸實(shí)現(xiàn)保鮮的CO2濃度。以上研究結(jié)果豐富了雞蛋呼吸相關(guān)理論，有助于禽蛋銷售前貯運(yùn)過程的保鮮研究。

表2 不同CO2濃度下雞蛋的呼吸強(qiáng)度與新鮮度

注：同行相同指標(biāo)不同小寫字母表示差異顯著（<0.05）。

Note: There is a significant differences in different lowercase letters for the same index (< 0.05).

3 結(jié) 論

1）一般情況下，雞蛋在受限空間內(nèi)釋放CO2的量逐漸增加，開始時，氣體在雞蛋表層進(jìn)行擴(kuò)散，然后由于重力作用，向下方沉積，最后再隨空氣流動逐漸上升，直至與空氣完全混合；不同CO2體積分?jǐn)?shù)下，雞蛋呼吸釋放CO2的擴(kuò)散過程基本一致，但雞蛋釋放CO2的量和速度有所不同，所以在相同時間點(diǎn)CO2分布具有一定差異，然而，當(dāng)CO2體積分?jǐn)?shù)增至3.0%后，雞蛋呼吸釋放的CO2分布無明顯區(qū)別。利用Fluent軟件計算所得氣體速度值和試驗(yàn)值基本一致，模擬值與實(shí)測值相對誤差在4%～9%之間，說明-模型較適合雞蛋呼吸釋放CO2擴(kuò)散過程分析。

2）根據(jù)不同CO2體積分?jǐn)?shù)下雞蛋新鮮度變化分析得知，隨著CO2體積分?jǐn)?shù)增加，貯藏20 d內(nèi)雞蛋哈夫單位下降速度逐漸變緩慢，當(dāng)CO2體積分?jǐn)?shù)達(dá)到3.0%時，再增加CO2體積分?jǐn)?shù)至4.5%，哈夫單位下降速度無明顯改善（>0.05），雞蛋新鮮度變化一致。因此，綜合經(jīng)濟(jì)因素，3.0%CO2保鮮雞蛋效果較優(yōu)。

[1]Leleu S, Herman L, Heyndrickx M, et al. Effects on Salmonella shell contamination and trans-shell penetration of coating hens' eggs with chitosan[J]. International Journal of Food Microbiology, 2011, 145(1): 43-48.

[2]王嬌嬌，王巧華，馬美湖，等. 雞蛋外部品質(zhì)與其呼吸強(qiáng)度的相關(guān)關(guān)系[J]. 食品科學(xué)，2018，39(1)：42-46. Wang Jiaojiao, Wang Qiaohua, Ma Meihu, et al. Correlation between egg external quality and respiration intensity[J]. Food Science, 2018, 19(1): 42-46. (in Chinese with English abstract)

[3]Pasquali F, Manfreda G, Olivi P, et al. Modified-atmosphere packaging?of?hen?table?eggs:?Effects?on?pathogen?and?spoilage?bacteria[J]. Poultry Science, 2012, 91(12): 3253–3259.

[4]Rocculi P, Cocci E, Sirrli F, et al. Modified atmosphere packaging of hen table eggs: Effects on functional properties of albumen[J]. Poultry Science, 2011, 90(8): 1791-1798.

[5]Jia F, Yan W J, Yuan X L, et al. Modified atmosphere packaging of eggs: Effects on the functional properties of albumen[J]. Food Packaging and Shelf Life, 2019, 22:100377-100377.

[6]Rocculi P, Tylewicz U, P?kos?awska A, et al. MAP storage of shell hen eggs, Part 1: Effect on physico-chemical characteristics of the fresh product[J]. Lwt Food Science & Technology, 2009, 42(3): 758–762.

[7]袁曉龍，高婧嫻，杜穎，等. 氣調(diào)儲藏中不同平衡氣體對雞蛋保鮮品質(zhì)的影響[J]. 食品工業(yè)科技，2014，35(2)： 300-303，307. Yuan Xiaolong, Gao Qianxian, Du Ying, et al. Effect of the balance gas of modified atmosphere packaging on the shelf-life of eggs[J]. Science and Technology of Food Industry, 2014, 35(2): 300-303, 307. (in Chinese with English abstract)

[8]劉美玉，連海平，任發(fā)政. 貯藏溫度和氣調(diào)包裝對雞蛋保鮮效果的影響[J]. 農(nóng)業(yè)工程學(xué)報，2011，27(12)：378-382. Liu Meiyu, Lian Haiping, Ren Fazheng. Effects of storage temperature and modified atmosphere packaging on egg preservation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2011, 27(12): 378-382. (in Chinese with English abstract)

[9]黃群，馬美湖，金永國，等. 貯存條件對蛋清S-卵白蛋白形成的影響[J].農(nóng)業(yè)工程學(xué)報，2012，28(05)：288-292. Huang Qun, Ma Meihu, Jin Yongguo, et al.Effect of storage condition on S-ovalbumin formation in albumen[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(05): 288-292. (in Chinese with English abstract)

[10]馬逸霄，王巧華，王賢波，等. 二氧化碳濃度對雞蛋呼吸強(qiáng)度及品質(zhì)的影響[J]. 分析化學(xué)，2020，48(12)：1728-1736. Ma Yixiao, Wang Qiaohua, Wang Xianbo, et al. Effect of concentration of carbon dioxide on respiration intensity and quality of eggs[J]. Chinese Journal of Analytical Chemistry, 2020, 48(12): 1728-1736. (in Chinese with English abstract)

[11]王巧華，張濤，馬美湖. 基于非損傷微測技術(shù)監(jiān)測貯期雞蛋氧呼吸規(guī)律[J]. 農(nóng)業(yè)工程學(xué)報，2014，30(5)：255-261. Wang Qiaohua, Zhang Tao, Ma Meihu. Oxygen breathing regularity of eggs in storage period monitoring by non-invasive micro-test technique[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(5): 255-261. (in Chinese with English abstract)

[12]王嬌嬌，王巧華，祝志慧，等. 雞蛋內(nèi)部品質(zhì)與其呼吸強(qiáng)度的關(guān)系[J]. 食品科學(xué)，2018，39(23)：36-40. Wang Jiaojiao, Wang Qiaohua, Zhu Zhihui, et al. Correlation between egg internal quality and respiration intensity[J]. Food Science, 2018, 39(23): 36-40. (in Chinese with English abstract)

[13]Guo Y, Wang G, Qian D, et al. Accident safety analysis of flow blockage in an assembly in the JRR-3Mresearch reactor using system code RELAP5 and CFD code FLUENT[J]. Annals of Nuclear Energy, 2018, 122: 125–136.

[14]Shen R, Jiao Z, Parker T, et al. Recent application of Computational Fluid Dynamics (CFD) in process safety and loss prevention: A review[J]. Journal of Loss Prevention in the Process Industries, 2020, 67(104252):1-80.

[15]Conradi M, Sánchez-Moyano J E, Galotti A, et al. CO2leakage simulation: Effects of the decreasing pH to the survival and reproduction of two crustacean species[J]. Marine Pollution Bulletin, 2019, 143:33-41.

[16]Perianu C, Ketelaere B D, Pluymers B, et al. Finite element approach for simulating the dynamic mechanical behaviour of a chicken egg[J]. Biosystems Engineering, 2010, 106(1): 79-85.

[17]Sellés A G, Marcè-Nogué J, Vila B, et al. Computational approach to evaluating the strength of eggs: Implications for laying in organic egg production[J]. Biosystems Engineering, 2019, 186: 146-155.

[18]Fabbri A, Cevoli C, Cocci E, et al. Determination of the CO2mass diffusivity of egg components by ?nite element model inversion[J]. Food Research International, 2011, 44(1): 204-208.

[19]岳佳淦. 氯氣泄漏擴(kuò)散數(shù)值模擬的研究[D]. 鄭州：鄭州大學(xué)，2016. Yue Jinfu. Studies on the Numerical Simulation of Chlorine Leakage and Diffusion[D]. Zhengzhou: Zhengzhou University, 2016. (in Chinese with English abstract)

[20]郭嘉. 小流量工況混流式水輪機(jī)壓力脈動研究[D]. 杭州：浙江大學(xué)，2013. Guo Jia. Simulation of Pressure Fluctuations under Part-load Operating Condition in A Francis Turbine[D]. Hangzhou: Zhejiang University, 2013. (in Chinese with English abstract)

[21]王嬌嬌，王巧華，曹芮. 不同貯存條件對雞蛋呼吸強(qiáng)度的影響[J]. 食品科技，2020，45(4)：26-30. Wang Jiaojiao, Wang Qiaohua, Cao Rui. Effects of different storage conditions on the respiration of eggs[J]. Food Science and Technology, 2020, 45(4): 26-30. (in Chinese with English abstract)

[22]崔啟筆. 基于FLUENT的氯氣泄漏擴(kuò)散的數(shù)值模擬與應(yīng)急預(yù)防的研究[D]. 焦作：河南理工大學(xué)，2009. Cui Qibi. Research on Numerical Simulation and Emergency Prevention of Chlorine Leakage and Diffusion Based on FLUENT[D]. Jiaozuo: Henan Polytechnic University, 2009. (in Chinese with English abstract)

[23]胡世明. 氣體釋放源的三維瞬態(tài)重氣擴(kuò)散模型及數(shù)值研究[D]. 北京：北京化工大學(xué)，2000. Hu Shiming. Three-Dimensional Transient “HeavyGas” Dispersion Model and Numerical Study[D]. Beijing: Beijing University of Chemical Technology, 2000. (in Chinese with English abstract)

[24]司凡. LNG氣化站氣相管道天然氣泄漏擴(kuò)散的數(shù)值模擬與分析[D]. 武漢：華中科技大學(xué)，2013. Si Fan. Numerical Simulation and Analysis of Natural Gas Leakage and Diffusion from LNG Vaporizing Station Gas Pipeline[D]. Wuhan: Huazhong University of Science and Technology, 2013. (in Chinese with English abstract)

[25]劉美玉，連海平，任發(fā)政. 不同貯藏溫度對雞蛋呼吸強(qiáng)度及品質(zhì)的影響[J]. 食品科學(xué)，2011，32(06)：270-274.

Liu Meiyu, Lian Haiping, Ren Fazheng. Effect of Storage Temperature on Respiratory Intensity and Quality of Eggs[J]. Food Science, 2011, 32(06): 270-274. (in Chinese with English abstract)

Simulation analysis and verification of egg respiration under different CO2concentrations

Wang Jiaojiao1, Wang Qiaohua2, Cao Rui1, Xie Juanjuan1

(1.,,466001,; 2.,,430070,)

Egg quality depends strictly on the storage environment after acquisition. The concentration of CO2has an important influence on the respiration and freshness of eggs during storage. This study aims to further observe the diffusion and freshness changes of CO2released from egg respiration under different CO2concentrations. The respiration intensity and freshness of stored eggs were measured at the temperature of 25°C, the relative humidity of 65%, the CO2volume fraction of 1.5%, 3.0%, 4.5%, and air (control group). The egg model was built by UG software. An ICEM CFD software was utilized to establish the fluid domain of the egg released carbon dioxide diffusion, and then divide it into a structural grid. A FLUENT software was selected to complete the diffusion simulation of the first day breathing release CO2from eggs with different CO2concentrations. The results show that the law of heavy gas diffusion was followed by the diffusion process of CO2released from eggs with different CO2volume fractions. Four stages were divided into: the initial, descending, mixing, and turbulent flow. The diffusion time in each stage was affected by the concentration of external CO2. Specifically, the diffusion amount and speed of CO2released by egg respiration gradually decreased, with the increase of CO2concentration in the storage environment. Egg respiration was inhibited, as the outside CO2concentration increased. There was a significant difference in the respiration intensity and freshness of eggs stored in the external environment with the CO2volume fraction of 1.5%, 3.0%, and 4.5% (<0.05) during the entire storage. The inhibitory effect on the respiration intensity of the egg was no longer obvious, and the freshness remained unchanged (>0.05) when the CO2volume fraction reached 3.0%. The freshness level of eggs was still above the AA level, when stored under the environment of CO2volume fraction of 3.0%, and 4.5% for 20 d. Taking economic factors into account, the CO2volume fraction of 3.0% was better to achieve the optimal freshness of eggs during storage. This study provides new ideas and data support for modified atmosphere preservation of eggs.

storage; respiration; eggs; CO2concentration; freshness; simulation

王嬌嬌，王巧華，曹芮，等. 不同CO2濃度下雞蛋呼吸仿真分析及驗(yàn)證[J]. 農(nóng)業(yè)工程學(xué)報，2021，37(6)：302-308.doi：10.11975/j.issn.1002-6819.2021.06.037 http://www.tcsae.org

Wang Jiaojiao, Wang Qiaohua, Cao Rui, et al. Simulation analysis and verification of egg respiration under different CO2concentrations[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(6): 302-308. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2021.06.037 http://www.tcsae.org

2020-09-19

2021-02-26

國家自然科學(xué)基金面上項(xiàng)目（31871863）；河南省科技廳科技發(fā)展計劃項(xiàng)目（212102110002）。

王嬌嬌，博士，講師，研究方向?yàn)檗r(nóng)產(chǎn)品加工與智能檢測。Email：1058864657@qq.com。

10.11975/j.issn.1002-6819.2021.06.037

TS253.2

A

1002-6819(2021)-06-0302-07