基于COMSOL Multiphysics的甜樱桃气调包装过程的呼吸数值模拟研究
|
摘要
以甜樱桃呼吸及其自发气调过程为研究对象,基于COMSOL Multiphysics软件建立其三维反应工程和传质扩散模型模拟,并优化其薄膜透气率参数。使用动态果蔬呼吸速率测定仪测定其呼吸速率,并与已有实验结果进行对比。基于酶动力学原理CO_2作为O_2的非竞争性抑制剂的米氏(Michaelis-Menten)呼吸的偏微分方程用来计算甜樱桃氧气的消耗和二氧化碳的产生,并用CO_2作为O_2的非竞争性抑制剂起到减缓呼吸的作用;Arrhenius表达式用于呼吸作用反应速率和反应级数的计算;可逆的反应工程表达式用于可逆反应的逆反应速率的计算;稀物质传递方程用O_2、CO_2的在罐体内、PE/LDPE薄膜内和空气中的气体扩散的计算;通过参数化扫描对薄膜气体透过率进行系列赋值,反向求解罐内的气体浓度以获得适合甜樱桃气体氛围的薄膜气体透过率。结果证明,动态果蔬呼吸速率测定仪的测定结果、已有实验结果和模型仿真的气体氛围变化结果较为一致。仿真构建的0～36 h的气体浓度分布和总通量分布,O_2在36 h的浓度为0.39%～0.44%;CO_2在36 h的浓度为19.06%～19.1%;O_2和CO_2浓度差在10 h最大,分别为888 mL/(m~3·h)和959 mL/(m~3·h)。查看参数化扫描结果,氧气扩散系数超过8.7145×10~(-11) m~2/s,二氧化碳扩散系数超过1.4815×10~(-11) m~2/s时,气调保鲜效果较好。下一步,将继续探索其他适宜使用气调的水果蔬菜的相关呼吸参数,可以依靠优化此模型较快的得出其最适透过率参数或扩散系数。
Based on COMSOL Multiphysics software,a three-dimensional reaction engineering and mass transfer model of sweet cherry were established to simulate its breathing and MAP (modified atmosphere preservation) process,and its membrane permeability parameters were optimized. Respiratory rate of sweet cherry was measured by dynamic respirator and compared with previous experimental results. Based on the principle of enzymatic kinetics,the partial differential equation of Michaelis-Menten respiration was used to calculate oxygen consumption and carbon dioxide production in sweet cherries,with CO_2 be used as a non-competitive inhibitor of oxygen to slow down respiration. Arrhenius expression was used to calculate the reaction rate and reaction order of respiration. Reversible reaction engineering expression was used to calculate the reversible reaction rate and order of respiration. The dilute matter transfer equation was used to calculate the gas diffusion of O_2 and CO_2 in the tank,PE/LDPE film and air. The gas permeability of the film was evaluated by a series of parametric scanning,and the gas concentration in the postharvest physiological tank was inversely calculated to obtain the film gas permeability which was suitable for sweet cherry. The results showed that the results of measurement instrument,previous experimental results and the simulation results were in good agreement. The simulated gas concentration distribution and total flux distribution for 0～36 h were 0.39%～0.44% for O_2 at 36 h,19.06%～19.1% for CO_2 at 36 h,and the maximum difference between CO_2 and CO_2 concentration at 10 h was 888 mL/(m~3·h)and 959 mL/(m~3·h),respectively. When the diffusion coefficient of O_2 and CO_2 exceeded 8.7145×10~(-11) m~2/s and 1.4815×10~(-11) m~2/s respectively,the effect of MAP was relatively good,based on parametric scanning result. The relevant breathing parameters of fruits and vegetables,which were suitable for MAP,would be explored,and the optimal transmittance parameters or diffusion coefficients would be attained by this model.
Based on COMSOL Multiphysics software,a three-dimensional reaction engineering and mass transfer model of sweet cherry were established to simulate its breathing and MAP (modified atmosphere preservation) process,and its membrane permeability parameters were optimized. Respiratory rate of sweet cherry was measured by dynamic respirator and compared with previous experimental results. Based on the principle of enzymatic kinetics,the partial differential equation of Michaelis-Menten respiration was used to calculate oxygen consumption and carbon dioxide production in sweet cherries,with CO_2 be used as a non-competitive inhibitor of oxygen to slow down respiration. Arrhenius expression was used to calculate the reaction rate and reaction order of respiration. Reversible reaction engineering expression was used to calculate the reversible reaction rate and order of respiration. The dilute matter transfer equation was used to calculate the gas diffusion of O_2 and CO_2 in the tank,PE/LDPE film and air. The gas permeability of the film was evaluated by a series of parametric scanning,and the gas concentration in the postharvest physiological tank was inversely calculated to obtain the film gas permeability which was suitable for sweet cherry. The results showed that the results of measurement instrument,previous experimental results and the simulation results were in good agreement. The simulated gas concentration distribution and total flux distribution for 0～36 h were 0.39%～0.44% for O_2 at 36 h,19.06%～19.1% for CO_2 at 36 h,and the maximum difference between CO_2 and CO_2 concentration at 10 h was 888 mL/(m~3·h)and 959 mL/(m~3·h),respectively. When the diffusion coefficient of O_2 and CO_2 exceeded 8.7145×10~(-11) m~2/s and 1.4815×10~(-11) m~2/s respectively,the effect of MAP was relatively good,based on parametric scanning result. The relevant breathing parameters of fruits and vegetables,which were suitable for MAP,would be explored,and the optimal transmittance parameters or diffusion coefficients would be attained by this model.
引文
[1]施俊凤,薛梦琳,王春生,等.甜樱桃采后生理特性与保鲜技术的研究现状与进展[J].保鲜与加工,2009(6):7-10.
[2]Pasquariello M S,Patre D D,Mastrobuoni F,et al.Influence of postharvest chitosan treatment on enzymatic browning and antioxidant enzyme activity in sweet cherry fruit[M].Postharvest Biology and Technology,2015,109:45-56.
[3]姜齐永,尤艳莉,周志才,等.自发气调保鲜“红灯”樱桃及数学模型的建立[J].食品研究与开发,2014,35(17):119-121.
[4]王宝刚,李文生,侯玉茹,等.甜樱桃物流及气调箱贮藏期间的品质变化[J].食品工业科技,2017,38(4)323-326.
[5]谢晶,刘晓丹.利用酶动力学拟合在密闭条件下香菇呼吸速率方程及米氏方程[J].中国食品学报,2006,6(1):96-100.
[6]万森,李喜宏,薛婷,等.温度对菠菜呼吸特性的影响[J].农业机械,2012(12):133-136.
[7]王亮,周建伟,徐剑,等.基于COMSOL Multiphysics模型的液态罐头食品热杀菌过程模拟研究[J].中国食品学报,2014(11):76-82.
[8]徐轶,徐青.基于COMSOL Multiphysics的渗流有限元分析[J].武汉大学学报:工学版,2014(2):165-170.
[9]周建伟,袁剑,王亮,等.基于COMSOL Multiphysics模型的瓶装液态食品巴氏灭菌过程的数值模拟研究[J].中国食品学报,2015(9):158-164.
[10]王刚,安琳.COMSOL Multiphysics工程实践与理论仿真-多物理场与数值仿真技术[M].北京:电子工业出版社,2012:76-80.
[11]Padmavati R,Anandharamakrishnan C.Computational fluid dynamics modeling of the thermal processing of cannedpine-apple slices and titbits[J].Food and Bioprocess Technology,2013,6(4):882-895.
[12]Pasquariello M S,Patre D D,Mastrobuoni F,et al.Influence of postharvest chitosan treatment on enzymatic browning and antioxidant enzyme activity in sweet cherry fruit[M].PostharvestBiology and Technology,2015,109:45-56.
[13]Jiang F,Weng J,Jia M,et al.Microstructural model in COMSOL packages with simulation to aging behavior of paper materials[J].Cellulose,2018,11:35-41.
[2]Pasquariello M S,Patre D D,Mastrobuoni F,et al.Influence of postharvest chitosan treatment on enzymatic browning and antioxidant enzyme activity in sweet cherry fruit[M].Postharvest Biology and Technology,2015,109:45-56.
[3]姜齐永,尤艳莉,周志才,等.自发气调保鲜“红灯”樱桃及数学模型的建立[J].食品研究与开发,2014,35(17):119-121.
[4]王宝刚,李文生,侯玉茹,等.甜樱桃物流及气调箱贮藏期间的品质变化[J].食品工业科技,2017,38(4)323-326.
[5]谢晶,刘晓丹.利用酶动力学拟合在密闭条件下香菇呼吸速率方程及米氏方程[J].中国食品学报,2006,6(1):96-100.
[6]万森,李喜宏,薛婷,等.温度对菠菜呼吸特性的影响[J].农业机械,2012(12):133-136.
[7]王亮,周建伟,徐剑,等.基于COMSOL Multiphysics模型的液态罐头食品热杀菌过程模拟研究[J].中国食品学报,2014(11):76-82.
[8]徐轶,徐青.基于COMSOL Multiphysics的渗流有限元分析[J].武汉大学学报:工学版,2014(2):165-170.
[9]周建伟,袁剑,王亮,等.基于COMSOL Multiphysics模型的瓶装液态食品巴氏灭菌过程的数值模拟研究[J].中国食品学报,2015(9):158-164.
[10]王刚,安琳.COMSOL Multiphysics工程实践与理论仿真-多物理场与数值仿真技术[M].北京:电子工业出版社,2012:76-80.
[11]Padmavati R,Anandharamakrishnan C.Computational fluid dynamics modeling of the thermal processing of cannedpine-apple slices and titbits[J].Food and Bioprocess Technology,2013,6(4):882-895.
[12]Pasquariello M S,Patre D D,Mastrobuoni F,et al.Influence of postharvest chitosan treatment on enzymatic browning and antioxidant enzyme activity in sweet cherry fruit[M].PostharvestBiology and Technology,2015,109:45-56.
[13]Jiang F,Weng J,Jia M,et al.Microstructural model in COMSOL packages with simulation to aging behavior of paper materials[J].Cellulose,2018,11:35-41.