# A Kinetic Model for Predicting the Relative Humidity in Modified Atmosphere Packaging and Its Application inLentinula edodes Packages

Authors: Li-xin Lu; Ya-li Tang; Su-yue Lu
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

Adjusting and controlling the relative humidity (RH) inside package is crucial for ensuring the quality of modified atmosphere packaging (MAP) of fresh produce. In this paper, an improved kinetic model for predicting the RH in MAP was developed. The model was based on heat exchange and gases mass transport phenomena across the package, gases heat convection inside the package, and mass and heat balances accounting for the respiration and transpiration behavior of fresh produce. Then the model was applied to predict the RH in MAP of freshLentinula edodes (one kind of Chinese mushroom). The model equations were solved numerically using Adams-Moulton method to predict the RH in model packages. In general, the model predictions agreed well with the experimental data, except that the model predictions were slightly high in the initial period. The effect of the initial gas composition on the RH in packages was notable. In MAP of lower oxygen and higher carbon dioxide concentrations, the ascending rate of the RH was reduced, and the RH inside packages was saturated slowly during storage. The influence of the initial gas composition on the temperature inside package was not much notable.


---

## Body

## 1. Introduction

Modified atmosphere packaging (MAP) is an important technology which can efficiently maintain quality and extend shelf life of produce. However, most polymeric films used in MAP have lower water vapor permeation rates relative to larger respiration rates of fresh produce; excessively high RH inside package may occur, and it may lead to undesirable results such as moisture condensation, microbial growth, and shortened shelf life [1–4]. Some methods were adopted to reduce the RH inside the package for extending shelf life of produce [5–8].Research on the mechanism of the RH change in MAP is very important for proper control of the RH and the quality of fresh produce. The RH in MAP is influenced by many factors such as respiration and transpiration of produce, the permeability of the packaging films, storage conditions, and packaging technology. Because of the RH inside the package with dynamic changes during storage, and the difficulty to accurately determine the respiration and transpiration process, little work has been done on predicting the RH in MAP.It is critical to establish the kinetic model of produce respiration and transpiration in order to predict and control the change of the RH in MAP. Transpiration rate model in MAP system has not been well developed mostly because transpiration phenomena on product surface due to respiratory heat has not been well acknowledged at changing O2 and CO2 environment [9]. A model was developed to predict the change of temperature and RH in MAP containing blueberry and moisture absorbent [9, 10]. A simple model to predict the quality changes of apples inside the packages was proposed [11]. Some effect factors and kinetic processes were neglected in the above models, such as the temperature difference of the fresh produce, the package headspace and the surrounding, gases heat convection inside the package, and heat change caused by gas transmission across the package.In MAP of nonrespiring foods, several mathematical models for gas diffusion and sorption processes and predicting the temperature change inside package were developed. Simpson proposed a model for gas diffusion and gas sorption processes [12] and then developed a model to predict consequences of temperature abuse on shelf life of pacific hake packaged in a MAP system [13]. The model included simultaneous heat and mass transfer phenomena coupled with a predictive shelf-life model and was solved using an explicit finite difference scheme. Simulated results were in good agreement with experimental results. Recently, these predictive models were successfully applied to study and determine the optimal conditions for packaging scallops in a modified atmosphere system [14].In this study, an improved kinetic model for predicting the RH in MAP of produce was developed. The model was based on heat exchange and gases mass transport phenomena across the package, gases heat convection inside the package, and mass and heat balances accounting for the respiration and transpiration behavior of fresh produce. To verify the model, a set of testing systems was established to measure the gas concentrations, RH, and temperature in MAP. The RH and temperature were valuated with experimental data of freshLentinula edodes(one kind of Chinese mushroom) packaged in modified atmospheres.


## 2. The Improved Kinetic Model for Predicting the RH in Map of Produce

Produce released respiratory energy during respiration in MAP. When released, respiratory heat partly is absorbed by produce leading to the rising of self temperature, and part of the heat is absorbed by free water of produce to vaporize water vapor. Meanwhile, the rest of respiratory heat is released to package.According to law of heat conservation, a heat balance relation inside package can be represented as(1)Qr=Qp+Qt+Qg+Qf+Qe,
where Qr is the respiratory heat of produce (J·h−1), Qp is the heat absorbed by produce (J·h−1), Qt is the heat absorbed by water vapor during evaporation (J·h−1), Qg is the heat absorbed by gas inside package (J·h−1), Qf is the heat absorbed by package (J·h−1), and Qe is the heat change caused by gas transmission across package (J·h−1).


### 2.1. Respiratory Heat of Produce

Considering the respiration heat of produce as the only internal heat, it could be described using the following respiratory equation [9]
(2)C6H12O6+6O2⟶6CO2+H2O+2816kJ
Respiration rate is generally described as O2 consumption rate RO2 or CO2 evolution rate RCO2. Respiration rate of produce could be estimated using Michaelis-Menten type of equations [15]: (3a)RO2=Vm1[O2]Km1+(1+[CO2]/Ku1)[O2],(3b)RCO2=Vm2[O2]Km2+(1+[CO2]/Ku2)[O2], where RO2, RCO2 are O2 consumption rate, CO2 evolution rate, respectively (mL·kg·h−1), [O2], [CO2] are O2, CO2 concentration inside package, respectively (%), Vm1, Vm2 are maximum O2, CO2 consumption rate, respectively (mL·kg·h−1), Km1, Km2 are Michaelis constants of O2, CO2 consumption, respectively (%), and ku1, ku2 are inhibition constants in O2, CO2 consumption, respectively (%).The quantities of O2 and CO2 inside package are  dynamically changeable in storage. On the basis of gas mass balance across the package, the change of O2, CO2 concentration inside package can be obtained [16] by
(4)d[O2]dt=(100)[PO2·Ap·patm([O2]out/100-[O2]in/100)z-RO2·Wp(([O2]in/100))z]×V-1,d[CO2]dt=(100)×[Pco2·Ap·patm([CO2]out/100-[CO2]in/100)z+RCO2·Wp(([CO2]in/100))z]×V-1,
where PO2, PCO2 are film permeabilities to O2, CO2, respectively (mL·m·m−2·h−1·atm−1), Ap is surface area of package (m2), z is the thickness of package film (m), Wp is the mass of produce (kg), patm is the pressure of 1 atm, and V is the free volume of package (mL).Many researches indicate that respiratory quotientRQ (RCO2/RO2) changes between 0.7 and 1.3. Therefore, respiration rates R of produce can be described by the average of O2 consumption and CO2 evolution rate [9, 11]:
(5)R=RCO2+RO22.Based on (2)–(5), the respiratory heat of produce in unit time can be represented as [9, 11]
(6)Qr=(28166)RWp.


### 2.2. Heat Absorbed by Produce

The heat absorbed by produce results in the change of temperature of produce. The mass of produce may decrease because of evaporation. Considering it as insignificant to the initial mass of produce, the loss of mass is negligible in this study. So the heat absorbed by produce in unit time can be represented as [17]
(7)Qp=CpWpdTpdt,
where Cp is the specific heat of surface tissue of produce (J·kg−1·k−1), Tp is the temperature on produce surface (k).


### 2.3. Heat Absorbed by Water Vapor during Evaporation

It is known that the heat absorbed by free water of unit mass to vaporize is definite in a certain condition. In unit time, the heat absorbed by water vapor during evaporation is [9, 11]
(8)Qt=Lmλ,
where Lm is the rate of water evaporation from produce to package (kg·h−1), λ is coefficient of evaporation heat (J·kg−1).


### 2.4. The Heat Absorbed by Gas inside Package

The heat absorbed by gas inside package consists of the following processes.(1)
Gas around produce inside package absorbs heat and results in rise in temperature.(2)
The difference in temperature of gas between near and far away from the produce makes convection of gas inside package and results in exchange of heat.


#### 2.4.1. Heat Absorbed by Gas around the Produce

The heat absorbed by gas around the produce in unit time is expressed by [17]
(9)Qg1=WsCsdTsdt,
where Ws is the mass of gas around the produce (kg), Cs is the specific heat of gas  (J·kg−1·k−1), and Ts is the temperature of gas around produce (k).


#### 2.4.2. Heat Exchanged by Convection of Gas inside Package

In the closed container, convection is one natural heat-releasing convection caused by the difference in temperature between gas near produce surface and free space inside package. Hence, Nusselt and Grashof rules could be used to describe the process.It is well known that Nusselt numberNu is the function of Grashof number Gr and Prandtl number Pr on the condition of natural heat-releasing convection. To simpler geometric shapes, the following equations could be applied to deal with these problems in engineering [17]:
(10)Nu=HLk=c(GrI·Pr)m,
where H is the coefficient of heat-exchange convection (J·h−1·m·k−2), L is the qualitative length of Grashof number (m), k is the coefficient of heat conduction of gas (J·h−1·m·k−1), c is the proportion constant, GrL is the Grashof number when the qualitative length equals L, and the value of exponent m, according to experience, is 1/4 for laminar flow and 1/2 for overfall flow [17].To simplify analysis, a series of experimental formulas have been proposed. For example, on the condition of the natural convection inside the vertical or horizontal closed cylinder,Nu can be expressed as [17]
(11)Nu≡HDk=0.55·(GrL·Pr)1/4,
where D is the diameter of cylinder (m).The coefficients of heat-exchanging convectionH can be obtained from (10)
(12)H=ck(GrL·Pr)mL.Therefore, the exchanged convective heat in unit time can be expressed as(13)Qg2=HA(Ts-Th),
where Th is the temperature far away from the heat-releasing surface (k), namely, the temperature of headspace gas inside package.As a result, in unit time, the heat absorbed by gas inside package is written as(14)Qg=Qg1+Qg2.


### 2.5. Heat Absorbed by Package

The heat absorbed by package is caused by the heat transfer produced by the difference in temperatures between inside and outside of the package. Rigid package consists of the packaging bracket and sealed film, while soft package is made up of films. For universality, rigid package is discussed in this study.The heat absorbed by a rigid package in unit time can be represented as(15)Qf=Qft+Qfm=WtCtdTtdt+WmCmdTmdt,
where Qft, Qfm are the heat absorbed by packaging bracket, sealed film, respectively (J·h−1), Wt, Wm is the mass of bracket, package film, respectively (kg), Ct, Cm is the specific heat of bracket material, sealed film, respectively (J·kg−1·k−1), and Tt, Tm is the temperature of bracket, sealed film, respectively (k).


### 2.6. Heat Change Caused by Gas Transmission across the Package

Gas exchanges exist across the packaging film because of the difference of gas pressures inside and outside the package during storage. It results in the heat change inside package.The gas energy is the sum of the internal energy and pressure static energy [18]. The internal energy of perfect gas, determined only by temperature, is produced by heat movement of gas molecules. Average movement energy Ek of a single gas molecule has a direct ratio relationship with absolute temperature [18]:
(16)Ek=i2kBT=i2n·Rc·T,
where i is the free degree of gas molecule, kB is Boltzmann constant (1.38×10-23 J·k−1), Rc is gas constant  (8.314 J·mol−1·k−1), T is the absolute temperature of gas (k), and n is the mol number of gas (mol).The pressure static energyEs of gas can be described as [18]
(17)Es=nRcT.The free degrees of CO2, H2O, and O2 molecules are 6, 12, and 5, respectively [18]. Generally, the gas environment in MAP contains lower O2 concentration, higher CO2 and water vapor concentrations relative to those of outside package. the energy which CO2 and water vapor permeate out of package in unit time can be described as
(18)Qg-out=(dnCO2dt×62+dnH2Odt×122)RcTh+(dnCO2dt+dnH2Odt)RcTh=(4dnCO2dt+7dnH2Odt)RcTh,
where nCO2, nH2O are the mol number of CO2, water vapor permeating out of package, respectively.Meanwhile, the energy which O2 permeates into package in unit time can be denoted as
(19)Qg-in=72dnO2dtRcTout,
where nO2 is the mol number of O2 permeating into package, Tout is storage temperature of package (k).The heat change caused by gas transmission across the package is(20)Qe=Qg-in-Qg-out.In terms of the results above, the heat balance relation in MAP can be represented as(21)RWp(Q6)=WpCpdTpdt+Lmλ+HA(Ts-Th)+WsCsdTsdt+WtCtdTtdt+WmCmdTmdt+72dnO2dtRcTout-(4dnCO2dt+7dnH2Odt)RcTh.


### 2.7. Predicting the RH in MAP

There are two sources that cause the change of quantity of water vapor inside package:(1)
the mass of water vapor produced by transpiration of produce,(2)
the mass of water vapor penetrated out or into packages across the package film.


#### 2.7.1. The Mass Change of Water Vapor Caused by Transpiration of Produce

The mass of water vapor produced by transpiration of produce in unit time, namely, transpiration rate of water vapor, can be obtained according to (21) as follows:
(22)Lm=[Q6RWp-WpCpdTpdt-HA(Ts-Th)-WsCsdTsdt-WtCtdTtdt-WmCmdTmdt-72dnO2dtRcTout+(4dnCO2dt+7dnH2Odt)RcTh]×λ-1.


#### 2.7.2. The Mass Change of Water Vapor Penetrated out or into Packages across the Package Film

The mass of water vapor permeating across the package in unit time is represented as [9]
(23)dmH2Odt=[PH2O·Ap·(pH2Oout-pH2Oin)z][0.018patmRcTh],
where mH2O is the mass of water vapor (kg), PH2O is the package film permeability to water vapor (mL·m·m−2·h−1·atm−1), and pH2Oout,pH2Oin are the water vapor pressures outside and inside package, respectively (atm).


#### 2.7.3. The Change of the RH inside Package

The changing rate of mass of water vapor inside the package can be described as(24)dHindt=Lm-dmH2O/dtWa,
where Hin is the absolute humidity inside package ((kg·H2O)·(kg−1·dry·H2O)), Wa is the mass of dry air inside package (kg).The RH inside package can be denoted as(25)RHin=HinHsp,
where Hsp is the absolute humidity of saturated water vapor at corresponding temperature ((kg·H2O)·(kg−1·dry·H2O)).


## 2.1. Respiratory Heat of Produce

Considering the respiration heat of produce as the only internal heat, it could be described using the following respiratory equation [9]
(2)C6H12O6+6O2⟶6CO2+H2O+2816kJ
Respiration rate is generally described as O2 consumption rate RO2 or CO2 evolution rate RCO2. Respiration rate of produce could be estimated using Michaelis-Menten type of equations [15]: (3a)RO2=Vm1[O2]Km1+(1+[CO2]/Ku1)[O2],(3b)RCO2=Vm2[O2]Km2+(1+[CO2]/Ku2)[O2], where RO2, RCO2 are O2 consumption rate, CO2 evolution rate, respectively (mL·kg·h−1), [O2], [CO2] are O2, CO2 concentration inside package, respectively (%), Vm1, Vm2 are maximum O2, CO2 consumption rate, respectively (mL·kg·h−1), Km1, Km2 are Michaelis constants of O2, CO2 consumption, respectively (%), and ku1, ku2 are inhibition constants in O2, CO2 consumption, respectively (%).The quantities of O2 and CO2 inside package are  dynamically changeable in storage. On the basis of gas mass balance across the package, the change of O2, CO2 concentration inside package can be obtained [16] by
(4)d[O2]dt=(100)[PO2·Ap·patm([O2]out/100-[O2]in/100)z-RO2·Wp(([O2]in/100))z]×V-1,d[CO2]dt=(100)×[Pco2·Ap·patm([CO2]out/100-[CO2]in/100)z+RCO2·Wp(([CO2]in/100))z]×V-1,
where PO2, PCO2 are film permeabilities to O2, CO2, respectively (mL·m·m−2·h−1·atm−1), Ap is surface area of package (m2), z is the thickness of package film (m), Wp is the mass of produce (kg), patm is the pressure of 1 atm, and V is the free volume of package (mL).Many researches indicate that respiratory quotientRQ (RCO2/RO2) changes between 0.7 and 1.3. Therefore, respiration rates R of produce can be described by the average of O2 consumption and CO2 evolution rate [9, 11]:
(5)R=RCO2+RO22.Based on (2)–(5), the respiratory heat of produce in unit time can be represented as [9, 11]
(6)Qr=(28166)RWp.


## 2.2. Heat Absorbed by Produce

The heat absorbed by produce results in the change of temperature of produce. The mass of produce may decrease because of evaporation. Considering it as insignificant to the initial mass of produce, the loss of mass is negligible in this study. So the heat absorbed by produce in unit time can be represented as [17]
(7)Qp=CpWpdTpdt,
where Cp is the specific heat of surface tissue of produce (J·kg−1·k−1), Tp is the temperature on produce surface (k).


## 2.3. Heat Absorbed by Water Vapor during Evaporation

It is known that the heat absorbed by free water of unit mass to vaporize is definite in a certain condition. In unit time, the heat absorbed by water vapor during evaporation is [9, 11]
(8)Qt=Lmλ,
where Lm is the rate of water evaporation from produce to package (kg·h−1), λ is coefficient of evaporation heat (J·kg−1).


## 2.4. The Heat Absorbed by Gas inside Package

The heat absorbed by gas inside package consists of the following processes.(1)
Gas around produce inside package absorbs heat and results in rise in temperature.(2)
The difference in temperature of gas between near and far away from the produce makes convection of gas inside package and results in exchange of heat.


### 2.4.1. Heat Absorbed by Gas around the Produce

The heat absorbed by gas around the produce in unit time is expressed by [17]
(9)Qg1=WsCsdTsdt,
where Ws is the mass of gas around the produce (kg), Cs is the specific heat of gas  (J·kg−1·k−1), and Ts is the temperature of gas around produce (k).


### 2.4.2. Heat Exchanged by Convection of Gas inside Package

In the closed container, convection is one natural heat-releasing convection caused by the difference in temperature between gas near produce surface and free space inside package. Hence, Nusselt and Grashof rules could be used to describe the process.It is well known that Nusselt numberNu is the function of Grashof number Gr and Prandtl number Pr on the condition of natural heat-releasing convection. To simpler geometric shapes, the following equations could be applied to deal with these problems in engineering [17]:
(10)Nu=HLk=c(GrI·Pr)m,
where H is the coefficient of heat-exchange convection (J·h−1·m·k−2), L is the qualitative length of Grashof number (m), k is the coefficient of heat conduction of gas (J·h−1·m·k−1), c is the proportion constant, GrL is the Grashof number when the qualitative length equals L, and the value of exponent m, according to experience, is 1/4 for laminar flow and 1/2 for overfall flow [17].To simplify analysis, a series of experimental formulas have been proposed. For example, on the condition of the natural convection inside the vertical or horizontal closed cylinder,Nu can be expressed as [17]
(11)Nu≡HDk=0.55·(GrL·Pr)1/4,
where D is the diameter of cylinder (m).The coefficients of heat-exchanging convectionH can be obtained from (10)
(12)H=ck(GrL·Pr)mL.Therefore, the exchanged convective heat in unit time can be expressed as(13)Qg2=HA(Ts-Th),
where Th is the temperature far away from the heat-releasing surface (k), namely, the temperature of headspace gas inside package.As a result, in unit time, the heat absorbed by gas inside package is written as(14)Qg=Qg1+Qg2.


## 2.4.1. Heat Absorbed by Gas around the Produce

The heat absorbed by gas around the produce in unit time is expressed by [17]
(9)Qg1=WsCsdTsdt,
where Ws is the mass of gas around the produce (kg), Cs is the specific heat of gas  (J·kg−1·k−1), and Ts is the temperature of gas around produce (k).


## 2.4.2. Heat Exchanged by Convection of Gas inside Package

In the closed container, convection is one natural heat-releasing convection caused by the difference in temperature between gas near produce surface and free space inside package. Hence, Nusselt and Grashof rules could be used to describe the process.It is well known that Nusselt numberNu is the function of Grashof number Gr and Prandtl number Pr on the condition of natural heat-releasing convection. To simpler geometric shapes, the following equations could be applied to deal with these problems in engineering [17]:
(10)Nu=HLk=c(GrI·Pr)m,
where H is the coefficient of heat-exchange convection (J·h−1·m·k−2), L is the qualitative length of Grashof number (m), k is the coefficient of heat conduction of gas (J·h−1·m·k−1), c is the proportion constant, GrL is the Grashof number when the qualitative length equals L, and the value of exponent m, according to experience, is 1/4 for laminar flow and 1/2 for overfall flow [17].To simplify analysis, a series of experimental formulas have been proposed. For example, on the condition of the natural convection inside the vertical or horizontal closed cylinder,Nu can be expressed as [17]
(11)Nu≡HDk=0.55·(GrL·Pr)1/4,
where D is the diameter of cylinder (m).The coefficients of heat-exchanging convectionH can be obtained from (10)
(12)H=ck(GrL·Pr)mL.Therefore, the exchanged convective heat in unit time can be expressed as(13)Qg2=HA(Ts-Th),
where Th is the temperature far away from the heat-releasing surface (k), namely, the temperature of headspace gas inside package.As a result, in unit time, the heat absorbed by gas inside package is written as(14)Qg=Qg1+Qg2.


## 2.5. Heat Absorbed by Package

The heat absorbed by package is caused by the heat transfer produced by the difference in temperatures between inside and outside of the package. Rigid package consists of the packaging bracket and sealed film, while soft package is made up of films. For universality, rigid package is discussed in this study.The heat absorbed by a rigid package in unit time can be represented as(15)Qf=Qft+Qfm=WtCtdTtdt+WmCmdTmdt,
where Qft, Qfm are the heat absorbed by packaging bracket, sealed film, respectively (J·h−1), Wt, Wm is the mass of bracket, package film, respectively (kg), Ct, Cm is the specific heat of bracket material, sealed film, respectively (J·kg−1·k−1), and Tt, Tm is the temperature of bracket, sealed film, respectively (k).


## 2.6. Heat Change Caused by Gas Transmission across the Package

Gas exchanges exist across the packaging film because of the difference of gas pressures inside and outside the package during storage. It results in the heat change inside package.The gas energy is the sum of the internal energy and pressure static energy [18]. The internal energy of perfect gas, determined only by temperature, is produced by heat movement of gas molecules. Average movement energy Ek of a single gas molecule has a direct ratio relationship with absolute temperature [18]:
(16)Ek=i2kBT=i2n·Rc·T,
where i is the free degree of gas molecule, kB is Boltzmann constant (1.38×10-23 J·k−1), Rc is gas constant  (8.314 J·mol−1·k−1), T is the absolute temperature of gas (k), and n is the mol number of gas (mol).The pressure static energyEs of gas can be described as [18]
(17)Es=nRcT.The free degrees of CO2, H2O, and O2 molecules are 6, 12, and 5, respectively [18]. Generally, the gas environment in MAP contains lower O2 concentration, higher CO2 and water vapor concentrations relative to those of outside package. the energy which CO2 and water vapor permeate out of package in unit time can be described as
(18)Qg-out=(dnCO2dt×62+dnH2Odt×122)RcTh+(dnCO2dt+dnH2Odt)RcTh=(4dnCO2dt+7dnH2Odt)RcTh,
where nCO2, nH2O are the mol number of CO2, water vapor permeating out of package, respectively.Meanwhile, the energy which O2 permeates into package in unit time can be denoted as
(19)Qg-in=72dnO2dtRcTout,
where nO2 is the mol number of O2 permeating into package, Tout is storage temperature of package (k).The heat change caused by gas transmission across the package is(20)Qe=Qg-in-Qg-out.In terms of the results above, the heat balance relation in MAP can be represented as(21)RWp(Q6)=WpCpdTpdt+Lmλ+HA(Ts-Th)+WsCsdTsdt+WtCtdTtdt+WmCmdTmdt+72dnO2dtRcTout-(4dnCO2dt+7dnH2Odt)RcTh.


## 2.7. Predicting the RH in MAP

There are two sources that cause the change of quantity of water vapor inside package:(1)
the mass of water vapor produced by transpiration of produce,(2)
the mass of water vapor penetrated out or into packages across the package film.


### 2.7.1. The Mass Change of Water Vapor Caused by Transpiration of Produce

The mass of water vapor produced by transpiration of produce in unit time, namely, transpiration rate of water vapor, can be obtained according to (21) as follows:
(22)Lm=[Q6RWp-WpCpdTpdt-HA(Ts-Th)-WsCsdTsdt-WtCtdTtdt-WmCmdTmdt-72dnO2dtRcTout+(4dnCO2dt+7dnH2Odt)RcTh]×λ-1.


### 2.7.2. The Mass Change of Water Vapor Penetrated out or into Packages across the Package Film

The mass of water vapor permeating across the package in unit time is represented as [9]
(23)dmH2Odt=[PH2O·Ap·(pH2Oout-pH2Oin)z][0.018patmRcTh],
where mH2O is the mass of water vapor (kg), PH2O is the package film permeability to water vapor (mL·m·m−2·h−1·atm−1), and pH2Oout,pH2Oin are the water vapor pressures outside and inside package, respectively (atm).


### 2.7.3. The Change of the RH inside Package

The changing rate of mass of water vapor inside the package can be described as(24)dHindt=Lm-dmH2O/dtWa,
where Hin is the absolute humidity inside package ((kg·H2O)·(kg−1·dry·H2O)), Wa is the mass of dry air inside package (kg).The RH inside package can be denoted as(25)RHin=HinHsp,
where Hsp is the absolute humidity of saturated water vapor at corresponding temperature ((kg·H2O)·(kg−1·dry·H2O)).


## 2.7.1. The Mass Change of Water Vapor Caused by Transpiration of Produce

The mass of water vapor produced by transpiration of produce in unit time, namely, transpiration rate of water vapor, can be obtained according to (21) as follows:
(22)Lm=[Q6RWp-WpCpdTpdt-HA(Ts-Th)-WsCsdTsdt-WtCtdTtdt-WmCmdTmdt-72dnO2dtRcTout+(4dnCO2dt+7dnH2Odt)RcTh]×λ-1.


## 2.7.2. The Mass Change of Water Vapor Penetrated out or into Packages across the Package Film

The mass of water vapor permeating across the package in unit time is represented as [9]
(23)dmH2Odt=[PH2O·Ap·(pH2Oout-pH2Oin)z][0.018patmRcTh],
where mH2O is the mass of water vapor (kg), PH2O is the package film permeability to water vapor (mL·m·m−2·h−1·atm−1), and pH2Oout,pH2Oin are the water vapor pressures outside and inside package, respectively (atm).


## 2.7.3. The Change of the RH inside Package

The changing rate of mass of water vapor inside the package can be described as(24)dHindt=Lm-dmH2O/dtWa,
where Hin is the absolute humidity inside package ((kg·H2O)·(kg−1·dry·H2O)), Wa is the mass of dry air inside package (kg).The RH inside package can be denoted as(25)RHin=HinHsp,
where Hsp is the absolute humidity of saturated water vapor at corresponding temperature ((kg·H2O)·(kg−1·dry·H2O)).


## 3. Predicting the RH in Map ofLentinula Edodes

### 3.1. Materials and Method

#### 3.1.1. FreshLentinula edodes and Packaging Film

Field-grownLentinula edodes used in this study were grown in Wuxi, China, and harvested at the stage of maturity in March, 2008. Lentinula edodes were then transported to the laboratory within one hour. They were selected for being similar in size and regular in shape.The low-density polyethylene, 32 µm thick (LDPE) (Shenda Flexible Packaging Material Co., Ltd., China), was used for packaging film. The film permeabilities to O2, CO2, and H2O were measured by gas permeability testers (Languang Instrument Ltd., China, Model BTY-B1,TSY-T, test precision: 5%) under 8°C. The film permeabilities to O2, CO2, and H2O are 5.66×10-13, 16.12×10-13, and 4.32×10-14 mL·m·m−2·h−1·atm−1, respectively.


#### 3.1.2. Testing System

As shown in Figure1, the testing system consisted of four parts that were testing jar, intelligent measurement system of RH-temperature, adjustment system of RH-temperature, and gas concentration analysis system. In order to measure the temperature, RH, and gas concentrations inside the metal jar, one RH-temperature sensor (RH: 10–100±2.0%, temperature: 0–60±0.5°C) and one gas pump were fixed on the upper portion of the jar and connected to intelligent measurement instrument of RH-temperature and gas concentration analysis instrument (Illinois Ltd., USA, model 6600, O2 range: 0.001–100±0.2%, CO2 range: 0.1–100±2.0%), respectively.Figure 1
Schematic of testing system (A, intelligent measurement system of RH-temperature; B, adjust system of RH-temperature; C, testing jar; D, gas concentration analysis instrument; 1, samples; 2, sensor of RH-temperature; 3, sealed film/metal cap; 4, pump of gases; 5, jar body).The testing jar can be closed with metal cap or packaging film in order to conduct the airtight or permeability packaging testing. After packaging, the testing jars were placed in adjuster of relative humidity-temperature (Qingshen Ltd., Taiwan, model THS-D7C-100AS, temperature:0–100±0.5°C, RH: 20–100±2.0%) to be stored in a testing condition.


#### 3.1.3. Estimation of Respiration Model Parameter Values

The samples and the jars were pretreated at 8°C, 80% RH for 4 h. The normal air packages were applied to pack570±10 g of samples using testing jars, which were closed with metal caps. The packages were stored at 8°C, 80% RH for 48 h. The concentrations of O2 and CO2 were measured every 2 h. The parameter values of respiration rates were calculated according to closed system method. The experiment was simultaneously repeated three times.Respiration rate data from above experiments were used for estimation of respiration model parameters values. Equations (3a) and (3b) were linearized as follows for the purpose of fitting the data:(26a)1RO2=1Vm1+Km1Vm11[O2]+1Ku1Vm1[CO2],(26b)1RCO2=1Vm2+Km2Vm21[O2]+1Ku2Vm2[CO2].The multiple linear regression analysis for estimating the parameter values of respiration model was performed using the statistical software SAS.


#### 3.1.4. Packaging  and Storage of the Samples

First, normal air packages were applied to pack480±10, 700±10 g of samples using test jars, which were closed with LDPE films of constant area 0.177 m2.In the second test, modified initial gas composition (3.5% O2/11.0% CO2/85.5% N2) packages were applied to pack 700±10 g of samples using testing jars. Gas mixtures were conducted using a gas mix machine (Desen Packaging Machine Ltd., China, model MA-35, CO2, O2, and N2 range: 0–90±1%). The machine used cylinders of compressed oxygen, carbon dioxide, and nitrogen to prepare the gas mixtures of any required proportions of oxygen, carbon dioxide, and nitrogen. Before packaging, the testing jars were flushed with the required gas composition (initial RH: 31%, temperature: 8°C) for five minutes through the lower valve of the jars, and then the samples were placed in jars. In the meantime, the flush was kept; finally the jars were closed with LDPE films.The samples, testing jars were conditioned at 8°C, 80% RH for 4 h before packaging. The packaging operation was conducted in a cold room at 8°C. All packages were stored at 8°C, 80% RH for 40 h, with periodic sampling to measure the RH and temperature. Three replicates were used for each condition.


### 3.2. Model Simplification and Numerical Analysis

#### 3.2.1. Simplification of the Kinetic Model

In terms of the above experiment conditions, the kinetic model can be simplified as follows.(1) Since the headspace of package is small, thermal equilibrium between the produce and the headspace is assumed to be reached within a short time. Consequently, the heat caused by gas convection inside the package could be neglected, and the heat absorbed by gas inside the package in unit time can be simplified as [18]
(27)Qg=WaCadTindt,
where Ca is the specific heat of dry air inside the package (J·kg−1·k−1), Tin is the temperature of gas inside the package (k).(2) Usually the temperature change inside package of produce is not notable [9]; its effect on the permeability of packaging film and the respiration model parameters can be ignored.


#### 3.2.2. Numerical Analysis

Based on the experimental conditions and corresponding parameters, the numerical analysis of the simplified kinetic model and corresponding equations was conducted using the Adams-Moulton method to obtain the change of RH inside package.


## 3.1. Materials and Method

### 3.1.1. FreshLentinula edodes and Packaging Film

Field-grownLentinula edodes used in this study were grown in Wuxi, China, and harvested at the stage of maturity in March, 2008. Lentinula edodes were then transported to the laboratory within one hour. They were selected for being similar in size and regular in shape.The low-density polyethylene, 32 µm thick (LDPE) (Shenda Flexible Packaging Material Co., Ltd., China), was used for packaging film. The film permeabilities to O2, CO2, and H2O were measured by gas permeability testers (Languang Instrument Ltd., China, Model BTY-B1,TSY-T, test precision: 5%) under 8°C. The film permeabilities to O2, CO2, and H2O are 5.66×10-13, 16.12×10-13, and 4.32×10-14 mL·m·m−2·h−1·atm−1, respectively.


### 3.1.2. Testing System

As shown in Figure1, the testing system consisted of four parts that were testing jar, intelligent measurement system of RH-temperature, adjustment system of RH-temperature, and gas concentration analysis system. In order to measure the temperature, RH, and gas concentrations inside the metal jar, one RH-temperature sensor (RH: 10–100±2.0%, temperature: 0–60±0.5°C) and one gas pump were fixed on the upper portion of the jar and connected to intelligent measurement instrument of RH-temperature and gas concentration analysis instrument (Illinois Ltd., USA, model 6600, O2 range: 0.001–100±0.2%, CO2 range: 0.1–100±2.0%), respectively.Figure 1
Schematic of testing system (A, intelligent measurement system of RH-temperature; B, adjust system of RH-temperature; C, testing jar; D, gas concentration analysis instrument; 1, samples; 2, sensor of RH-temperature; 3, sealed film/metal cap; 4, pump of gases; 5, jar body).The testing jar can be closed with metal cap or packaging film in order to conduct the airtight or permeability packaging testing. After packaging, the testing jars were placed in adjuster of relative humidity-temperature (Qingshen Ltd., Taiwan, model THS-D7C-100AS, temperature:0–100±0.5°C, RH: 20–100±2.0%) to be stored in a testing condition.


### 3.1.3. Estimation of Respiration Model Parameter Values

The samples and the jars were pretreated at 8°C, 80% RH for 4 h. The normal air packages were applied to pack570±10 g of samples using testing jars, which were closed with metal caps. The packages were stored at 8°C, 80% RH for 48 h. The concentrations of O2 and CO2 were measured every 2 h. The parameter values of respiration rates were calculated according to closed system method. The experiment was simultaneously repeated three times.Respiration rate data from above experiments were used for estimation of respiration model parameters values. Equations (3a) and (3b) were linearized as follows for the purpose of fitting the data:(26a)1RO2=1Vm1+Km1Vm11[O2]+1Ku1Vm1[CO2],(26b)1RCO2=1Vm2+Km2Vm21[O2]+1Ku2Vm2[CO2].The multiple linear regression analysis for estimating the parameter values of respiration model was performed using the statistical software SAS.


### 3.1.4. Packaging  and Storage of the Samples

First, normal air packages were applied to pack480±10, 700±10 g of samples using test jars, which were closed with LDPE films of constant area 0.177 m2.In the second test, modified initial gas composition (3.5% O2/11.0% CO2/85.5% N2) packages were applied to pack 700±10 g of samples using testing jars. Gas mixtures were conducted using a gas mix machine (Desen Packaging Machine Ltd., China, model MA-35, CO2, O2, and N2 range: 0–90±1%). The machine used cylinders of compressed oxygen, carbon dioxide, and nitrogen to prepare the gas mixtures of any required proportions of oxygen, carbon dioxide, and nitrogen. Before packaging, the testing jars were flushed with the required gas composition (initial RH: 31%, temperature: 8°C) for five minutes through the lower valve of the jars, and then the samples were placed in jars. In the meantime, the flush was kept; finally the jars were closed with LDPE films.The samples, testing jars were conditioned at 8°C, 80% RH for 4 h before packaging. The packaging operation was conducted in a cold room at 8°C. All packages were stored at 8°C, 80% RH for 40 h, with periodic sampling to measure the RH and temperature. Three replicates were used for each condition.


## 3.1.1. FreshLentinula edodes and Packaging Film

Field-grownLentinula edodes used in this study were grown in Wuxi, China, and harvested at the stage of maturity in March, 2008. Lentinula edodes were then transported to the laboratory within one hour. They were selected for being similar in size and regular in shape.The low-density polyethylene, 32 µm thick (LDPE) (Shenda Flexible Packaging Material Co., Ltd., China), was used for packaging film. The film permeabilities to O2, CO2, and H2O were measured by gas permeability testers (Languang Instrument Ltd., China, Model BTY-B1,TSY-T, test precision: 5%) under 8°C. The film permeabilities to O2, CO2, and H2O are 5.66×10-13, 16.12×10-13, and 4.32×10-14 mL·m·m−2·h−1·atm−1, respectively.


## 3.1.2. Testing System

As shown in Figure1, the testing system consisted of four parts that were testing jar, intelligent measurement system of RH-temperature, adjustment system of RH-temperature, and gas concentration analysis system. In order to measure the temperature, RH, and gas concentrations inside the metal jar, one RH-temperature sensor (RH: 10–100±2.0%, temperature: 0–60±0.5°C) and one gas pump were fixed on the upper portion of the jar and connected to intelligent measurement instrument of RH-temperature and gas concentration analysis instrument (Illinois Ltd., USA, model 6600, O2 range: 0.001–100±0.2%, CO2 range: 0.1–100±2.0%), respectively.Figure 1
Schematic of testing system (A, intelligent measurement system of RH-temperature; B, adjust system of RH-temperature; C, testing jar; D, gas concentration analysis instrument; 1, samples; 2, sensor of RH-temperature; 3, sealed film/metal cap; 4, pump of gases; 5, jar body).The testing jar can be closed with metal cap or packaging film in order to conduct the airtight or permeability packaging testing. After packaging, the testing jars were placed in adjuster of relative humidity-temperature (Qingshen Ltd., Taiwan, model THS-D7C-100AS, temperature:0–100±0.5°C, RH: 20–100±2.0%) to be stored in a testing condition.


## 3.1.3. Estimation of Respiration Model Parameter Values

The samples and the jars were pretreated at 8°C, 80% RH for 4 h. The normal air packages were applied to pack570±10 g of samples using testing jars, which were closed with metal caps. The packages were stored at 8°C, 80% RH for 48 h. The concentrations of O2 and CO2 were measured every 2 h. The parameter values of respiration rates were calculated according to closed system method. The experiment was simultaneously repeated three times.Respiration rate data from above experiments were used for estimation of respiration model parameters values. Equations (3a) and (3b) were linearized as follows for the purpose of fitting the data:(26a)1RO2=1Vm1+Km1Vm11[O2]+1Ku1Vm1[CO2],(26b)1RCO2=1Vm2+Km2Vm21[O2]+1Ku2Vm2[CO2].The multiple linear regression analysis for estimating the parameter values of respiration model was performed using the statistical software SAS.


## 3.1.4. Packaging  and Storage of the Samples

First, normal air packages were applied to pack480±10, 700±10 g of samples using test jars, which were closed with LDPE films of constant area 0.177 m2.In the second test, modified initial gas composition (3.5% O2/11.0% CO2/85.5% N2) packages were applied to pack 700±10 g of samples using testing jars. Gas mixtures were conducted using a gas mix machine (Desen Packaging Machine Ltd., China, model MA-35, CO2, O2, and N2 range: 0–90±1%). The machine used cylinders of compressed oxygen, carbon dioxide, and nitrogen to prepare the gas mixtures of any required proportions of oxygen, carbon dioxide, and nitrogen. Before packaging, the testing jars were flushed with the required gas composition (initial RH: 31%, temperature: 8°C) for five minutes through the lower valve of the jars, and then the samples were placed in jars. In the meantime, the flush was kept; finally the jars were closed with LDPE films.The samples, testing jars were conditioned at 8°C, 80% RH for 4 h before packaging. The packaging operation was conducted in a cold room at 8°C. All packages were stored at 8°C, 80% RH for 40 h, with periodic sampling to measure the RH and temperature. Three replicates were used for each condition.


## 3.2. Model Simplification and Numerical Analysis

### 3.2.1. Simplification of the Kinetic Model

In terms of the above experiment conditions, the kinetic model can be simplified as follows.(1) Since the headspace of package is small, thermal equilibrium between the produce and the headspace is assumed to be reached within a short time. Consequently, the heat caused by gas convection inside the package could be neglected, and the heat absorbed by gas inside the package in unit time can be simplified as [18]
(27)Qg=WaCadTindt,
where Ca is the specific heat of dry air inside the package (J·kg−1·k−1), Tin is the temperature of gas inside the package (k).(2) Usually the temperature change inside package of produce is not notable [9]; its effect on the permeability of packaging film and the respiration model parameters can be ignored.


### 3.2.2. Numerical Analysis

Based on the experimental conditions and corresponding parameters, the numerical analysis of the simplified kinetic model and corresponding equations was conducted using the Adams-Moulton method to obtain the change of RH inside package.


## 3.2.1. Simplification of the Kinetic Model

In terms of the above experiment conditions, the kinetic model can be simplified as follows.(1) Since the headspace of package is small, thermal equilibrium between the produce and the headspace is assumed to be reached within a short time. Consequently, the heat caused by gas convection inside the package could be neglected, and the heat absorbed by gas inside the package in unit time can be simplified as [18]
(27)Qg=WaCadTindt,
where Ca is the specific heat of dry air inside the package (J·kg−1·k−1), Tin is the temperature of gas inside the package (k).(2) Usually the temperature change inside package of produce is not notable [9]; its effect on the permeability of packaging film and the respiration model parameters can be ignored.


## 3.2.2. Numerical Analysis

Based on the experimental conditions and corresponding parameters, the numerical analysis of the simplified kinetic model and corresponding equations was conducted using the Adams-Moulton method to obtain the change of RH inside package.


## 4. Results and Discussion

### 4.1. Respiration Model Parameter Values ofLentinula edodes

The respiration rate data obtained from the closed system experiments were applied to fit (26a) and (26b). The estimated parameter values were presented in Table 1. The regression equation fitted the data well with the coefficient of determination R2>0.95.Table 1
Parameter estimates of respiration model forLentinula edodes at 8°C.
Gas
V
m (mL·kg−1·h−1)
K
m (%)
K
u (%)
R
2
O2
122.52
26.73
−88.30
0.9824
CO2
110.81
30.11
−83.23
0.9575


### 4.2. Verification of the Kinetic Model

The model predictions were compared to the results of three sets of experiments. The experiments varied in product weight and initial gas composition.Figures2 and 3 compared the model predictions and experiment data for 480,700 g, Lentinula edodes and normal air packages, respectively. The results showed that the RH inside packages rose rapidly. In the packages with the masses of 480, 700 g, the RH inside packages reached saturation (100% RH) 12 h, 10 h later, respectively. The main reason was that the gas exchange inside and outside the package was slow because of the initial gas of air and the respiration and transpiration rates were a bit higher. As the mass of produce increased, the quality of vaporized moisture was enhanced, the headspace decreased, and the rising rate of the RH quickened. Meanwhile, the results showed good agreement between theoretical prediction and experimental results with a little higher theoretical result at the beginning of package and the time for RH to reach saturated state shortened about 1 h mainly because of the simplification in theoretical analysis and neglects of influences of some factors on energy changes.Figure 2
Relative humidity in packages oflentinula edodesduring storage (mass: 480 g; initial gas composition : normal air).Figure 3
Relative humidity in packages ofLentinula edodes during storage (mass: 700 g; initial gas composition: normal air).Figure4 compared the model prediction and experiment data for 700 g Lentinula edodes and modified initial gas composition packages. Comparing with the package of same mass with initial gas of air, the increasing rate of the RH inside package was slowed down greatly. About 22 h later, the RH reached saturation. The main reason is that modified atmosphere of low CO2 and high O2 decreased the respiration rates of produce and furthermore the transpiration rates were slowed down. This shows that initial gas composition could influence the RH inside package evidently. Similarly, the theoretical prediction agreed well with the experimental results.Figure 4
Relative humidity in packages ofLentinula edodes during storage (mass: 700 g; initial gas composition: 3.5% O2/11.0% CO2/85.5% N2).The temperature change in three experimental conditions was shown in Figure5. As a whole, the change of temperature inside packages was not much notable. During storage of 40 h, corresponding to the experimental condition, the highest temperatures inside packages rose approximately 3.3°C, 4.3°C, and 2.0°C, respectively. Meanwhile, temperature inside packages decreased a little along with the prolonging of storage. This situation showed that the respiratory heat of sample was mainly absorbed to vaporize water vapor and result in the rapid rise of RH inside packages. The slightly increased temperatures justified our simplification of the kinetic model.Figure 5
Temperature in packages ofLentinula edodes during storage.


## 4.1. Respiration Model Parameter Values ofLentinula edodes

The respiration rate data obtained from the closed system experiments were applied to fit (26a) and (26b). The estimated parameter values were presented in Table 1. The regression equation fitted the data well with the coefficient of determination R2>0.95.Table 1
Parameter estimates of respiration model forLentinula edodes at 8°C.
Gas
V
m (mL·kg−1·h−1)
K
m (%)
K
u (%)
R
2
O2
122.52
26.73
−88.30
0.9824
CO2
110.81
30.11
−83.23
0.9575


## 4.2. Verification of the Kinetic Model

The model predictions were compared to the results of three sets of experiments. The experiments varied in product weight and initial gas composition.Figures2 and 3 compared the model predictions and experiment data for 480,700 g, Lentinula edodes and normal air packages, respectively. The results showed that the RH inside packages rose rapidly. In the packages with the masses of 480, 700 g, the RH inside packages reached saturation (100% RH) 12 h, 10 h later, respectively. The main reason was that the gas exchange inside and outside the package was slow because of the initial gas of air and the respiration and transpiration rates were a bit higher. As the mass of produce increased, the quality of vaporized moisture was enhanced, the headspace decreased, and the rising rate of the RH quickened. Meanwhile, the results showed good agreement between theoretical prediction and experimental results with a little higher theoretical result at the beginning of package and the time for RH to reach saturated state shortened about 1 h mainly because of the simplification in theoretical analysis and neglects of influences of some factors on energy changes.Figure 2
Relative humidity in packages oflentinula edodesduring storage (mass: 480 g; initial gas composition : normal air).Figure 3
Relative humidity in packages ofLentinula edodes during storage (mass: 700 g; initial gas composition: normal air).Figure4 compared the model prediction and experiment data for 700 g Lentinula edodes and modified initial gas composition packages. Comparing with the package of same mass with initial gas of air, the increasing rate of the RH inside package was slowed down greatly. About 22 h later, the RH reached saturation. The main reason is that modified atmosphere of low CO2 and high O2 decreased the respiration rates of produce and furthermore the transpiration rates were slowed down. This shows that initial gas composition could influence the RH inside package evidently. Similarly, the theoretical prediction agreed well with the experimental results.Figure 4
Relative humidity in packages ofLentinula edodes during storage (mass: 700 g; initial gas composition: 3.5% O2/11.0% CO2/85.5% N2).The temperature change in three experimental conditions was shown in Figure5. As a whole, the change of temperature inside packages was not much notable. During storage of 40 h, corresponding to the experimental condition, the highest temperatures inside packages rose approximately 3.3°C, 4.3°C, and 2.0°C, respectively. Meanwhile, temperature inside packages decreased a little along with the prolonging of storage. This situation showed that the respiratory heat of sample was mainly absorbed to vaporize water vapor and result in the rapid rise of RH inside packages. The slightly increased temperatures justified our simplification of the kinetic model.Figure 5
Temperature in packages ofLentinula edodes during storage.


## 5. Conclusions

An improved kinetic model was developed to predict the RH in MAP of produce based on heat exchange and gases mass transport phenomena across the package, gases heat convection inside the package, and mass and heat balances accounting for the respiration and transpiration behavior of fresh produce. The experiments of MAP containing freshLentinula edodes were conducted to obtain the RH and temperature inside packages. The results showed that(1)
as a whole, the model predictions agreed well with experimental results,(2)
initial gas composition had notable influence on the RH inside packages,(3)
the change of temperature inside package and effect of initial gas composition on temperature were not much notable.



---

*Source: 304016-2013-09-19.xml*