柱状茭白真空预冷热质传递分析

奥鹏网院作业

柱状茭白真空预冷热质传递分析柱状茭白真空预冷热质传递分析

王 宁，阚安康*，黄孜沛，陆家希，潘雨芹

（上海海事大学商船学院，上海 201306）

摘 要：茭白味道鲜美、营养价值高，但其含水率高达90%以上，易腐变质，货架期短；为延长其货架期，本实验对柱状茭白的真空预冷进行研究，综合考虑了非稳态传热和质量守恒，建立了柱状茭白真空预冷的热质传递数学模型，并对茭白内部热质传递进行计算流体动力学模拟预测；模拟了真空室内压力变化和茭白质量变化及温度分布情况，并将模拟结果与实验数据进行了对比分析。结果表明，模拟结果与实验数据相吻合：模拟与实验压力误差为7.4%，质量损失误差在8%以内，内外最大平均温度误差仅3.7%。该模型可以预测柱状茭白真空预冷过程中热质传递的变化，对柱状蔬菜真空环境下的水分迁移特性研究和食品工业真空预冷工艺的改进具有借鉴意义。

关键词：真空预冷；柱状茭白；热质传递；计算流体动力学模拟

真空预冷是食品保鲜的重要技术之一，也是冷链物流过程中重要的支撑技术，与其他冷却方法相比其冷却时间十分短暂，有着明显的优势[1-6]。其原理是真空条件可使材料中水分快速蒸发，蒸发的水蒸气带走材料潜热并向外界释放，从而达到使材料迅速降温的目的[7-9]。在真空预冷过程中，食品温度分布和失水率一直是业界关注的问题。温度分布不均可能导致低温损伤，对食品品质产生一定负面影响，进而导致经济损失[10]。数学模型是有效的模拟工具，可以用来从理论上描述真空预冷过程，有利于改进设备、优化工艺和控制产品品质[11-13]；针对肉类、蔬菜和木材等不同材料真空预冷过程的模拟，国内外学者已建立了众多数学模型，用于预测系统压力、温度和产品失水率等各种参数。Torres等[14]提出了一种理想固体多孔介质真空干燥模型的数值求解方法，研究介质内部水蒸气和空气的动力变化情况。Salehi等[15]建立真空干燥数学模型，研究了红外-真空联合干燥对柠檬片的干燥动力学、水分扩散性、表面收缩和颜色变化的影响。结果表明，数学模型和实验结果较为吻合。金听祥等[16]为研究肉类在真空冷却中的温度分布和压力的变化情况，建立圆柱形肉块水分迁移的真空冷却数值模型，结果表明模拟结果和实验数据基本一致，温度最大误差在5%以内。Bahrani等[17]建立了淀粉材料真空降压过程中热质传递的物理模型，解释了淀粉加工过程中热量和质量传递的变化机理。Sun Dawen等[18]使用计算流体软件建立了多孔材料真空冷却的非稳态数值模型，并通过多孔食品数值模拟和实验结果的对比，成功验证了真空冷却过程中多孔食品的热质传递。Drummond等[19]对熟牛肉浸泡式真空冷却过程建立了数学模型，研究牛肉尺寸、孔隙率大小与工艺参数对熟牛肉样品在真空冷却期间的冷却时间、总质量损失、压力和温度分布的影响，并成功进行了模型验证。贺素艳等[20]对甘蓝真空预冷过程进行了理论模拟，考虑装置泄漏面积，通过实验计算出了甘蓝实际传质系数，这为其他真空预冷系统的泄漏面积计算提供了新方法。王雪芹等[21]以卷心菜为研究对象，建立了真空冷却过程中球形果蔬的非稳态传热传质模型，并考虑了蒸发、对流和辐射等因素对温降的影响，同时通过实验成功验证了模型的准确性。韩志等[22]提出了真空预冷过程中球形食品的传热模型，并利用UDF自定义编程，引入食品材料的物性参数，使用FLUENT软件对食品真空预冷过程进行了数值模拟。阚安康[23]、陈雪梅[24]等建立了多孔蔬菜真空预冷传热传质数值模型，通过模拟数据和实验结果的对比，验证了柱状蔬菜降温速率和质量损失变化情况。He Zhenbin等[25]建立了木材真空干燥过程中各时间点的热质耦合传递模型，并对模型进行了验证，结果表明，水分含量高于纤维饱和点木材的平均干燥速率约为水分含量低于纤维饱和点木材的2.91 倍。Demarchi等[26]建立了玫瑰果皮真空干燥过程中水分、温度和抗坏血酸含量的变化模型，对热量和质量损失进行了较好的预测。Song Xiaoyan等[27]利用数学模型模拟真空预冷中切叶蔬菜的传热过程，并对鲜切青菜真空预冷时的表面温度分布及其叶部和叶柄的平均表面温度进行了预测，在红外热成像技术的帮助下，通过理论和实验数据的对比，成功地验证了数学模型。

在当前的教学模式下，教师在讲授相关理论知识时，普遍采用的是视频文件、幻灯片、文档等形式，而学校往往因为各类因素的限制，导游课程教学与实训难以与实际结合起来，学生与社会需求的人才普遍脱节，导致学生社会适应能力下降。

茭白味道鲜美、营养价值高，但其含水率高达90%以上，易腐烂变质，货架期短[28]。为解决茭白保鲜问题，在流入消费市场前，可以对茭白进行真空预冷处理。然而真空条件下蔬菜失水的现象不可避免，真空预冷会对茭白品质产生影响。因此，本实验从理论上分析真空下柱状蔬菜热质传递的过程，并以柱状茭白为研究对象，建立柱状茭白热质传递数值模型，通过实验数据与模拟结果对比来检验模型的准确性，这对茭白保鲜和延长茭白市场货架期具有一定的理论和指导意义，同时也为同类柱状蔬菜在真空预冷保鲜过程中提供相应借鉴。

1 材料与方法1.1 材料

新鲜茭白为市售。

1.2 仪器与设备

ZLG0.1A型真空预冷试验机 上海善如水保鲜科技有限公司；HZY-B1000电子天平 福州华志科学仪器有限公司。

真空预冷装置系统如图1所示，主要由真空泵和压缩机组成，通过配备T型热电偶和压力传感器的数据采集系统，对蔬菜内部温湿度和真空室内压力进行连续测量和采集。真空室内尺寸500 mm×410 mm×490 mm；最低真空度为400 Pa，显示精度为1 Pa；预冷最低温度-5 ℃，在-5～40 ℃之间可任意设定，显示精度0.1 ℃，相对湿度传感器精度0.5%。

     

图1 真空预冷系统图
Fig. 1 Schematic diagram of the vacuum cooling device

1.3 方法

1.3.1 热质传递模型的建立

真空预冷过程主要分为3 个阶段。第一阶段是恒温和降压阶段，在此阶段样品温度基本不变，真空室内压力逐渐下降。当压力下降到材料初始温度对应的饱和压力（闪点）时，材料中的自由水开始吸热蒸发，材料即迅速降温，降温效果明显，此为真空预冷的第二阶段。随着水分的不断蒸发和向外排放，材料温度将迅速下降到设定的温度工况。第三阶段是压力恢复阶段，此阶段预冷装置进气阀打开，外界空气流入真空室内直至内外压力平衡[29]。综上，主要的热质传递发生在第二阶段，故本研究的热质传递模型就是基于这一阶段所建立。

多孔材料真空预冷中的各部位压强条件和汽化强度是相同的，其温降速率亦相同；在真空预冷过程中，热量传递的方式有相变热（即材料水分蒸发吸收的热量）、热传导（即材料内部的导热）、对流传热（即材料外表面与真空室内气体、腔体内壁的热交换）、辐射热（即真空室内材料与气体、腔体内壁之间的换热量）。辐射换热量与材料自身温度、黑度和辐射系数有关，在本实验模拟中热量小，可被忽略。

为了简化计算，在建立柱状茭白数学模型时作如下假设：1）真空预冷装置密闭性好，不考虑外界空气的渗透；2）辐射换热量少，忽略不计；3）模型初始温度和水分含量分布均匀，模拟过程中材料的热物性参数恒定不变且各向同性；4）材料吸热产生的水蒸气完全被捕水器捕获。

1.3.1.1 传热模型

根据以上假设条件，茭白真空预冷过程的非稳态传热控制方程可用公式（1）表示。考虑到新鲜的茭白是一个有机生命体，进行的呼吸作用会产生呼吸热，故传热过程涉及的内热源可由公式（2）表示。

   

式中：ρf为茭白密度/（kg/m3）；cp为茭白比热容/（kJ/（kg·K））；λ为茭白导热系数/（W/（m·K））；T为温度/K；t为时间/s；x、r、θ分别为柱坐标系；φ为内热源/（kJ/kg），由蒸发潜热和呼吸热两部分组成；hfg为水的汽化潜热/（kJ/kg）；fv为茭白中水蒸气的产生速率/（kg/（m·s））；fbh为茭白自身的呼吸热/（kJ/（m·s））。

以上事例均为同一航向才能使用TCPA进行对遇、追越判断，不同航向则需心算大致的TCPA，并留有3～5分钟的安全余量。所有的追越应在被追越船同意并配合下小心追越。无论如何，都必须特别谨慎地驾驶，直到追越完成为止。

根据公式（3）～（5）确定内热源。

   

式中：ε为茭白的孔隙率；M为相对分子质量；R0为通用气体常数，取为8.314 J/（mol·K）；da为多孔介质的茭白孔径/m；μ为水蒸气的动力黏度/（kPa·s）；P为真空室内环境压力/Pa；ξx、ξr和ξθ分别为x、r、θ所对应的水蒸气迁移阻力因子。

1.3.1.2 传质模型

真空预冷过程中水分达到闪点后，将从茭白内部吸热蒸发，茭白单位体积的蒸发速率可用公式（6）[16]表示。

2）培训机构层面：为学习者提供功能完善的网络研修平台，并对学习者定期做学习分析，及时对学习者进行反馈。

   

式中：D为茭白的直径/m；hm为茭白的沸腾系数，取为8.4×10-7 kg/（Pa·m2·s）；Ps为茭白在温度T下对应的饱和压力/Pa；Pv为真空室内的压力/Pa。

Ps与茭白温度T的关系可用公式（7）表示。

   

真空室内压力Pv随时间的变化关系用公式（8）表示。

基于此，提出基于词向量与CNN(Convolutional Neural Network)结合的中文文本语义相似度计算方法WV-CNN，并应用到中文答案选择系统中，选取2种典型的传统语义相似度计算模型与其进行对比，实验结果表明：WV-CNN模型效果良好，对当前的中文文本语义相似度计算有很大的指导意义。

国外对修建过鱼设施的共识度高、起步早，技术较为成熟。如美国在哥伦比亚河干支流修建了16座鱼道，成功地解决了鲑科鱼类的洄游问题。国内关于修建过鱼设施仍存争议，有学者认为，过鱼对象主要是为了解决上溯繁殖的鱼类洄游问题，葛洲坝工程修建后，已没有必须上溯繁殖的洄游性鱼类。但从遗传多样性保护层面来看，阻隔会影响坝上、坝下种群的基因交流，造成种群的遗传分化，加剧物种的濒危，建造过鱼设施还是有一定的必要性。

   

式中：S表示真空泵的抽气速率/（m3/s）；Vf表示真空室的容积/m3。

1.3.1.3 初始条件和边界条件

李站长说，我这里好说，顶多搞点罚款过去，要是被市林业公安的发现了，可能我也脱不了干系的。何泽说，李站长尽管放心，市里的工作，由买主去做，负责没有障碍。

大中型轴流泵站由于叶轮直径大，一般要求运行转速较低。而低转速交流异步电动机的功率因数较低，因此,目前多采用同步电动机作为主电机。 同步电动机的特点是：其运行转速恒定不变，除非改变其电源频率f或电机的磁极对数p。因此，轴流泵站反向发电运行时，其运行方式一般有：同转速运行、变极降速运行、变频运行3种方式。

初始条件如式（9）所示。

   

式中：T0为茭白初始温度/K；P0为茭白所处的真空室内初始压力/Pa。

边界条件如式（10）所示。

   

式中：ω为茭白真空预冷产生的水蒸气质量分数/%。

茭白表面边界条件如式（11）所示。

   

式中：n、qr和qc分别表示茭白的表面法向量、辐射热和对流换热/（kJ/kg），其中数值求解时辐射热和对流换热可忽略，即qr＝0，qc＝0。

1.3.1.4 茭白的热物性参数

茭白的导热系数按照公式（12）[30]求取。

   

式中：ωw、ωp、ωc、ωf和ωa分别表示茭白中水分、蛋白质、碳水化合物、脂肪和灰分的质量分数/%。

在进行屋顶施工作业过程中，必须合理的应用具有保温性能的新型材料。想要阻隔顶层以及外界环境；就要凸显顶面种植工艺的特征优势，在保障防水性能的基础上，着重分析建筑的承载范围，合理进行顶层种植，不仅会提升资源利用效率，还让屋顶具备了自动调节温度的功能，使得夏季降温冬季保暖的功效。同时，在建筑顶层还可以合理应用太阳能技术，在供暖以及制冷系统中应用太阳能技术，进而避免冬季过于寒冷、夏季过于炎热的问题。

老师的即时评价不仅仅只是评价功能，还应有启发功能。在教学过程中老师要用启发性的评价语言引导学生动脑思考，让学生在探寻学习方法过程中发现问题，及时纠正自己的行为偏差。通过渗透教育，培养学生良好的学习习惯。

茭白的比热容按照公式（13）[30]求取。

   

实验中孔隙率取0.63，茭白多孔介质的平均孔径取0.002 5 m[18]，其他指标参数值如表1所示，其中，茭白质量由电子天平测定；各成分质量分数委托上海海洋大学食品学院实验测定；直径由直尺测量；密度由质量与体积比求平均值；初始温度、初始压力由真空预冷设备自带传感器测定。将各数据代入FLUENT软件求解。

表1 真空预冷过程中的水分迁移模型参数
Table 1 Model parameters for moisture migration during vacuum cooling

     

指标 数值初始质量（m）/kg 0.139 32灰分质量分数（ωa）/% 0.52碳水化合物质量分数（ωc）/% 5.9脂肪质量分数（ωf）/% 0.2蛋白质量分数（ωp）/% 1.2水分质量分数（ωw）/% 91.1直径（D）/m 0.04密度（ρf）/（kg/m3） 628导热系数（λ）/（W/（m·K）） 0.55比热容（cp）/（kJ/（kg·K）） 4.05初始温度（T0）/K 297初始压力（P0）/Pa 101 325

1.3.2 验证实验

1）将购买的茭白进行筛选，选择半径2 cm、长18 cm的茭白适当去皮并称质量，确保其几何大小与模型尺寸参数一致。2）将温度传感器的两个探头分别插进茭白中心、表面位置（图2A）；茭白样品初始温度为24 ℃，最终温度设置为2 ℃。3）启动设备，制冷机和真空泵同时运行，直至茭白中心处最终温度达到设定值，系统每隔10 s自动采集压力、温度及相对湿度等数据并记录。4）重复1～3过程。5）实验测试工作结束，关闭实验系统。

由于结核病属于人畜共患的疾病，需要加强对卫生方面的监督，严格按照卫生检验制度执行，尤其是畜产品，如牛乳产品等，目的是减少消化道感染[4]。每年定期对从事牛养殖的工作人员进行健康检查，尤其是人畜共患的传染性疾病，若结核病检测结果为阳性，需要将其立即调离岗位，并给予有效的治疗，促进牛养殖业的健康发展，同时保障牛乳产品的质量。人员在对病牛饲养环境进行消毒灭菌时必须佩戴一定的防护措施，防止细菌感染。

     

图2 茭白的测试点位置及装置
Fig. 2 Test locations and equipment for Zizania latifolia

2 结果与分析2.1 模型网格划分

真空预冷中的材料温度和压力都是随时间而不断变化的，因此这是一个非稳态过程。本实验中柱状模型使用隐式差分法来求解茭白真空预冷中的热质传递微分方程。首先划分网格，建立x、r、θ柱坐标系，再将计算区域和时间坐标分别均匀划分为（a-1）、（b-1）份，因此空间节点为a 个（1≤a≤A），而时间节点b（0≤b≤n）可由n×Δi表示，Δi为时间间隔。模型中a＝0位置为茭白的中心，令其为第一节点，茭白表面为最大节点A。基于能量守恒定律、质量守恒定律以及提出的柱状数学模型，针对茭白内某一微元体Δr/2（△r为空间两相邻节点的距离）处的热质平衡，可得茭白内部节点的有限差分方程，如式（14）所示。

   

式中：Tab＋1和Tab分别表示同一个空间节点在两相邻时间下的温度/K；Tab和Tba-1分别表示在相同时间下的相邻空间节点的温度/K；Δr为空间两相邻的节点距离。

在r＝0，即a＝0时，茭白中心节点的差分方程如式（15）所示。

   

式中：T0b＋1和T0b分别表示中心节点在两相邻时间下的温度/K；T1b和T0b分别表示在相同时间下的中心及其相邻空间节点的温度/K。

在r＝R，即a＝A时，茭白表面节点的差分方程如式（16）所示。

   

式中：Rw为水的气体常数/（J/kg）；D0为在温度273 K下一个标准大气压下水蒸气在空气中的扩散系数（0.22×104 m2/s）；l为传质阻力层厚度/m；Paw为真空室内环境水蒸气分压/Pa；aw为水分活度。

吡嗪类具有坚果和焙烤等香气，是酱油酱香味的来源。本试验所检出吡嗪类化合物的相对含量较高，以样品7最为显著，以该种方式添加产香对酱油的酱香味贡献最大。

基于上述数学模型分析，建立茭白柱状数值模型，ICEM-CFD（the integrated computer engineering and manufacturing code for computational fluid dynamics）网格划分如图3所示。通过FLUENT软件，根据上述各参数的设置，对柱状蔬菜的真空预冷过程进行模拟求解。

咳嗽跟呼吸的空气有很大关系。如果空气太脏，或者是太干燥，宝宝就咳得厉害。所以，在家里护理宝宝时，若天气好，应多开门窗通气，使屋里的空气保持清新；若天气不好，可以使用空气净化器来改善室内的空气；若空气干燥，要用加湿器增加室内湿度，使用加湿器使室内湿度保持在40%～50%，湿度太大也不行，容易使房间滋生霉菌。也可以选择让宝宝吸入水蒸汽的方式缓解咳嗽。睡前在浴室内放会儿热水，待蒸汽充满浴室，把宝宝抱进去尽可能多待一些时间，让呼吸道通过多吸入一些水蒸汽获得充分的滋润，这个方法也有助于缓解鼻塞和咳嗽。还可以用妈妈们蒸脸的蒸汽机让宝宝的呼吸道滋润，但蒸汽机里不能使用自来水或矿泉水，要使用蒸馏水。

     

图3 柱状模型网格划分
Fig. 3 Cylinder model meshing

2.2 模拟茭白温度变化     

图4 茭白真空预冷后的内部温度云图
Fig. 4 Temperature profiles inside Zizania latifolia after vacuum cooling

真空预冷过程中茭白横向和纵向切面的温度变化云图如图4所示。在模拟初始时刻（0 s），茭白的起始温度为297 K，此时预冷设备未启动，茭白温度与真空室内环境之间没有温度差，故茭白温度保持不变。当进行300 s后，在真空泵的抽吸作用下，真空室内气体被连续排除，茭白外表面与环境之间发生强制对流，室内压力和流体温度也因此降低；当茭白表面温度达到其压力对应饱和水的蒸发温度时（闪点），茭白外层水分开始吸热蒸发，其温度由内到外依次迅速降低，从图4b能够看出茭白的温度分布是以芯部为中心的圆弧圈。随着模拟的进行，茭白外表面温度在735 s时最先到达设定值（275 K），此时茭白芯部温度在内外温差驱动下主要依赖导热缓慢下降，这是因为茭白外层水分的闪发，使得整个茭白周围充满水蒸气，即增加了茭白内层水分的闪发阻力。在真空泵的持续抽吸作用下，真空室内的水蒸气不断被排出，茭白内层水分在压力差的驱动下向外层迁移，同时伴随着吸热蒸发，导致茭白芯部的温度也迅速降低。800 s时，柱状茭白真空预冷热质传递过程模拟结束，茭白停止吸热蒸发，其外表面和内部温度基本相同。

2.3 真空室内压力和温度变化验证实验结果     

图5 茭白压力和温度模拟值与实验值的对比
Fig. 5 Comparison of simulated and experimental values for pressure in the vacuum chamber and internal and surface temperature of Zizania latifolia

真空预冷过程中实验数据和模拟数据的压力变化曲线如图5a所示。在真空预冷初始阶段，真空室内的模拟压力与实验压力均以较高的压降速率下降，但前者数值稍低；在预冷中期，茭白内部水分到达闪点后吸热蒸发，真空室内压降放缓，到预冷后期两者曲线基本吻合并保持水平不变。这是因为模拟真空预冷过程的真空泵功率是维持不变的平均功率，因此，在降压初始阶段，真空泵工作效率高，抽气速率快，降压速率高，且不考虑外界空气的渗入；然而在验证实验过程中，真空泵的初始功率低，真空室内压力下降速率较小，而且外界的空气会向真空室内渗透，抽气速率也相对较慢。随着真空泵的连续工作，验证实验中的真空泵功率趋于恒定，200 s后模拟值和实验值的压力误差逐渐缩小，最大偏差不超过100 Pa。真空预冷过程结束时，室内压力均保持恒定不变，模拟压力为635 Pa，实验压力为686 Pa，模拟与实验压力误差仅为7.4%，两者压力变化曲线基本吻合。

茭白表面和中心位置的温度模拟数据与实验数据的曲线变化如图5b所示。在预冷初始阶段，真空室内空气被真空泵连续抽吸并向外排出，室内压力持续下降，直至达到茭白内部水分闪点，这一阶段主要是茭白与真空室内的空气发生对流换热，换热量微弱，茭白温度几乎没有下降。随后达到闪点的茭白外层水分在压力差作用下最先吸热蒸发，温度急剧下降，而这时茭白内层主要依赖导热降温。随着预冷过程的进行，茭白温度不断下降，当达到设定工况时，其内外层温度基本相同。真空预冷过程结束时，茭白终点处中心温度实验值比模拟值高，表面温度实验值比模拟值低，但茭白内外平均温度的模拟值要比实验值偏高，两者在整个真空预冷过程的最大误差仅为3.7%，这是因为验证实验所用的茭白都是经过机械鲜切处理的，这种机械损伤既破坏了茭白的组织结构，又改变了茭白内部水分迁移的通透性；因此，在真空预冷时，实验中的茭白水分更易流失，温度下降也更快。

2.4 真空室内质量损失和相对湿度变化验证实验结果     

图6 茭白质量（a）和相对湿度（b）模拟值与实验值的对比
Fig. 6 Comparison of simulated and experimental values for mass (a)and relative humidity (b) inside Zizania latifolia

构成蔬菜组织结构的水分包括结合水和自由水，结合水是指通过化学键与溶质分子、非水组织牢固结合的一类水，因此，在真空预冷过程出现损耗的水都是蔬菜组织结构和细胞中可以自由移动的自由水。在预冷过程中，茭白内部水分（自由水）达到闪点时会大量吸热蒸发并向外界排出，因此茭白的质量也会降低。如图6a所示，在真空预冷前期，验证实验与模拟茭白质量损失误差在1%以内；但在真空预冷后期，两者质量误差逐渐增大。真空预冷过程结束时验证实验茭白质量损失为6.3%，而模拟数据显示的质量损失为5.8%，两者间的误差仅为7.9%。一方面，是由于茭白在预冷实验中期达到闪点后内部水分大量吸热蒸发，蒸发的水蒸气包围整个茭白，阻碍了水分进一步的蒸发损失；另一方面，到实验后期，真空泵的连续抽吸使得真空室内气体被迅速排至外界，在压力差驱动下，外层水分的蒸发增加了茭白中的孔隙数量，结果内层水分向外转移能力得到提升的同时也能顺利吸热蒸发。

目前，医院统计工作在各医院中还没有引起足够的重视，医院的统计数据也没有在工作中发挥应有作用，对此，医院要充分认识到统计数据的重要性，改变现有的统计工作态度和方法，并积极创新工作方法，提升数据统计质量。

如图6b所示，在预冷起步阶段，茭白所处环境的相对湿度为66%，在预冷230 s时，真空室内相对湿度第一次达到最低值（24%）；而在220～400 s，真空室内相对湿度出现上升的态势并达到41%；400 s以后相对湿度又开始缓慢降低，第800秒时，真空室内相对湿度第二次达到最低值（19%），随后相对湿度维持不变。这是因为预冷初期真空室内气体被真空泵迅速抽吸排至外界，此时茭白的内部水分尚未达到闪点，真空室内相对湿度也因此第一次迅速到达最低点；随着预冷过程的进行，茭白内部水分在230 s后到达闪点，水分吸热快速蒸发，大量水蒸气包围茭白附近区域，使得真空室内水蒸气分压力升高，相对湿度也随之上升；随着水蒸气的生成与排出，茭白内外层的自由水含量也不断减少，所以400 s后真空室内相对湿度又开始缓慢下降；800 s后茭白内部水分基本完全蒸发，茭白温度达到设定值，室内相对湿度也第二次达到最低点并维持不变。

2.5 误差分析

该数学模型能有效模拟柱状茭白热质传递过程，但此模型存在误差，主要原因有：1）实验所用的茭白都是经过机械鲜切处理的，而且实验中热电偶需要插孔，这些无疑破坏了茭白的组织结构，进而对实验数据产生影响。2）由于茭白的上下端、内外部的密度和形态各有差异，因此对于建立的柱状数值模型，茭白难以真正实现物质分布均匀，整体一致。

3 结 论

本实验以茭白为柱状蔬菜的典型代表，建立了真空预冷热质传递模型，并分析了柱状茭白真空预冷实验的真空室内压力、相对湿度和茭白温度分布以及质量损失的情况，通过模拟和验证实验数据的对比，成功验证了数学模型。模拟与验证实验结果表明：1）在真空预冷初始阶段，真空室内降压速率快，模拟压力比实验压力稍低，茭白与真空室内空气温度基本相同，对流换热微弱，茭白温度变化不明显；茭白水分到达闪点后吸热蒸发，真空室内压降放缓，而降温速率增加；真空预冷过程中，茭白外表面部分降温速率高于中心处部分，模拟与验证实验结果中压力和温度变化曲线基本吻合。2）在真空预冷过程中，受真空室内水蒸气分压影响，相对湿度出现两次最低值，第一次在闪点之前，为24%，第二次在闪点之后，为19%。3）在真空预冷初期，模拟质量损失数据和验证实验数据相差不明显，误差在1%以内；真空预冷结束时，真空室内压力模拟值为635 Pa，实验值为686 Pa，模拟数值与实验数值误差仅为7.4%；温度模拟值与实验值的最大平均温度误差仅为3.7%。实验结果的质量损失为6.3%，模拟结果的质量损失为5.8%，两者误差在8%以内。

参考文献：

[1] 付艳武, 高丽朴, 王清, 等. 蔬菜预冷技术的研究现状[J]. 保鲜与加工, 2015, 15(1): 58-63. DOI:10.3969/j.issn.1009-6221.2015.01.013.

[2] 李靖, 谢如鹤, 刘广海, 等. 冷藏运输用新型低温相变材料及装备的研制[J]. 制冷学报, 2018, 39(4): 35-40. DOI:10.3969/j.issn.0253-4339.2018.04.032.

[3] 严灿, 刘升, 贾丽娥, 等. 蔬菜冷链物流技术研究进展[J]. 食品与机械,2015, 31(4): 260-265. DOI:10.13652/j.issn.1003-5788.2015.04.064.

[4] 吴欣蔚, 朱志伟, 孙大文. 三种形态叶菜真空预冷过程的比较研究[J]. 现代食品科技, 2017, 33(9): 134-139. DOI:10.13982/j.mfst.1673-9078.2017.9.019.

[5] ZHU Zhiwei, LI Ying, SUN Dawen, et al. Developments of mathematical models for simulating vacuum cooling processes for food products: a review[J]. Critical Reviews in Food Science and Nutrition,2019, 59(5): 715-727. DOI:10.1080/10408398.2018.1490696.

[6] 郑先章. 关于减压贮藏技术及理论主流观点的商榷[J]. 农业工程学报, 2017, 33(14): 9-18. DOI:10.11975/j.issn.1002-6819.2017.14.001.

[7] 陈颖, 刘宝林, 宋晓燕. 荷兰豆真空预冷及其对贮藏品质的影响[J].食品科学, 2013, 34(6): 276-279. DOI:10.7506/spkx1002-6630-201306061.

[8] ZHANG Zhihua, DRUMMOND L, SUN Dawen. Vacuum cooling in bulk of beef pieces of different sizes and shape: evaluation and comparison to conventional cooling methods[J]. Journal of Food Engineering, 2013,116(2): 581-587. DOI:10.1016/j.jfoodeng.2012.12.036.

[9] 宋晓燕, 刘宝林. 食品真空冷却技术研究进展[J]. 食品科学, 2014,35(11): 319-324. DOI:10.7506/spkx1002-6630-201411062.

[10] SUN Dawen, ZHENG Liyun. Vacuum cooling technology for the agrifood industry: past, present and future[J]. Journal of Food Engineering,2006, 77(2): 203-214. DOI:10.1016/j.jfoodeng.2005.06.023.

[11] CAMPA Ñ ONE L A, ZARITZKY N E. Mathematical modeling and simulation of microwave thawing of large solid foods under different operating conditions[J]. Food and Bioprocess Technology, 2010, 3(6):813-825. DOI:10.1007/s11947-009-0249-0.

[12] AZARPAZHOOH E, RAMASWAMY H S. Modeling and optimization of microwave osmotic dehydration of apple cylinders under continuous-flow spray mode processing conditions[J]. Food and Bioprocess Technology, 2012, 5(5): 1486-1501. DOI:10.1007/s11947-010-0471-9.

[13] CAMPA Ñ ONE L A, PAOLA C A, MASCHERONI R H. Modeling and simulation of microwave heating of foods under different process schedules[J]. Food and Bioprocess Technology, 2012, 5(2): 738-749.DOI:10.1007/s11947-010-0378-5.

[14] TORRES S S, JOMAA W, PUIGGALI J R, et al. Multiphysics modeling of vacuum drying of wood[J]. Applied Mathematical Modelling, 2011, 35(10): 5006-5016. DOI:10.1016/j.apm.2011.04.011.

[15] SALEHI F, KASHANINEJAD M. Modeling of moisture loss kinetics and color changes in the surface of lemon slice during the combined infrared-vacuum drying[J]. Information Processing in Agriculture,2018, 5(4): 516-523. DOI:10.1016/j.inpa.2018.05.006.

[16] 金听祥, 张海川, 李改莲, 等. 熟肉真空冷却过程中水分迁移理论分析和实验[J]. 农业工程学报, 2008, 24(8): 309-312. DOI:10.3321/j.issn:1002-6819.2008.08.069.

[17] BAHRANI S A, MONTEAU J Y, REZZOUG S A, et al. Physicsbased modeling of simultaneous heat and mass transfer intensification during vacuum steaming processes of starchy material[J]. Chemical Engineering and Processing: Process Intensification, 2014, 85: 216-226. DOI:10.1016/j.cep.2014.09.005.

[18] SUN Dawen, HU Zehua. CFD simulation of coupled heat and mass transfer through porous foods during vacuum cooling process[J].International Journal of Refrigeration, 2003, 26(1): 19-27.DOI:10.1016/s0140-7007(02)00038-5.

[19] DRUMMOND L, SUN Dawen. Evaluation of the immersion vacuum cooling of cooked beef joints: mathematical simulation of variations in beef size and porosity and pressure reduction rates[J]. Innovative Food Science & Emerging Technologies, 2012, 16: 205-210. DOI:10.1016/j.ifset.2012.06.001.

[20] 贺素艳, 张光存, 王翠苹, 等. 果蔬真空预冷过程中传质系数的实验研究[J]. 青岛大学学报(工程技术版), 2011, 26(3): 86-90.DOI:10.3969/j.issn.1006-9798.2011.03.018.

[21] 王雪芹, 刘宝林, 闫静文. 卷心菜真空冷却过程的CFD传热传质模拟[J]. 制冷学报, 2012, 33(1): 65-68. DOI:10.3969/j.issn.0253-4339.2012.01.065.

[22] 韩志, 谢晶, 潘迎捷. 球性食品真空冷却非稳态过程的数值分析[J].真空科学与技术学报, 2010, 30(6): 657-661. DOI:10.3969/j.issn.1672-7126.2010.06.16.

[23] 阚安康, 韩厚德. 柱状蔬菜冷藏运输前真空预冷热质迁移数学模型及实验[J]. 上海海事大学学报, 2010, 31(3): 58-62. DOI:10.3969/j.issn.1672-9498.2010.03.011.

[24] 陈雪梅, 韩厚德, 阚安康. 多孔蔬菜在真空冷却过程中热质流失的研究[J]. 安徽农学通报, 2011, 17(5): 143-145. DOI:10.3969/j.issn.1007-7731.2011.05.069.

[25] HE Zhenbin, QIAN Jing, QU Lijie, et al. Simulation of moisture transfer during wood vacuum drying[J]. Results in Physics, 2019, 12:1299-1303. DOI:10.1016/j.rinp.2019.01.017.

[26] DEMARCHI S M, IRIGOYEN R M T, GINER S A. Vacuum drying of rosehip leathers: modelling of coupled moisture content and temperature curves as a function of time with simultaneous timevarying ascorbic acid retention[J]. Journal of Food Engineering, 2018,233: 9-16. DOI:10.1016/j.jfoodeng.2018.03.027.

[27] SONG Xiaoyan, LIU Baolin, JAGANNATHAN G K. Mathematical simulation on the surface temperature variation of fresh-cut leafy vegetable during vacuum cooling[J]. International Journal of Refrigeration, 2015, 65: 228-237. DOI:10.1016/j.ijrefrig.2015.12.009.

[28] 颜建春, 胡志超, 吴朋来, 等. 热板-微波联合真空冷冻干燥茭白工艺优化[J]. 农业工程学报, 2017, 33(1): 262-270. DOI:10.11975/j.issn.1002-6819.2017.01.036.

[29] AUGUSTO C M, RIBEIRO J B, GASPAR A R, et al. A mathematical model describing the two stages of low-pressure-vaporization of free water[J]. Journal of Food Engineering, 2012, 112(4): 274-281.DOI:10.1016/j.jfoodeng.2012.05.013.

[30] SWEAT V E. Thermal properties of foods[M]// RAO M A, RIZIL S S H.Engineering properties of foods. New York: Marcel Dekker, Inc.,1986: 49-87.

Analysis of Heat and Mass Transfer of Cylindrical Pieces of Zizania latifolia during Vacuum Precooling

WANG Ning, KAN Ankang*, HUANG Zipei, LU Jiaxi, PAN Yuqin
(Merchant Marine College, Shanghai Maritime University, Shanghai 201306, China)

Abstract: Zizania latifolia has delicious taste and high nutrition value, but is perishable with a short shelf life due to the water content of over 90%. In order to extend the shelf life, a mathematical model to describe heat and mass transfer in cylindrical pieces of Zizania latifolia during vacuum cooling was developed taking into consideration both transient heat transfer and mass conservation, and using computational fluid dynamics, simulation and prediction of heat and mass transfer in Zizania latifolia were conducted. Pressure variations inside the vacuum chamber, material mass variations,and temperature distribution were also simulated and compared with the experimental data. The results indicated that the simulation results agreed with the experimental data. The deviation between simulated and experimental values for pressure,mass loss and average temperature difference between the inside and the outside was 7.4%, less than 8% and only 3.7%,respectively. The model proved to be useful to predict the change of heat and mass transfer in cylindrical pieces of vegetables during vacuum cooling, which will provide useful references to study the water migration characteristics of cylindrical pieces of vegetables under vacuum environment and improve the vacuum cooling process in the food industry.

Keywords: vacuum precooling; cylindrical pieces of Zizania latifolia; heat and mass transfer; computational fluid dynamics simulation

收稿日期：2019-05-07

基金项目：上海市自然科学基金面上项目（15ZR1419900）

第一作者简介：王宁（1990—）（ORCID: 0000-0002-6505-3862），男，硕士研究生，研究方向为多孔介质传热传质。E-mail: ning12578@126.com

*通信作者简介：阚安康（1981—）（ORCID: 0000-0002-5333-7702），男，高级工程师，博士，研究方向为多孔介质传热传质。E-mail: ankang0537@126.com

DOI:10.7506/spkx1002-6630-20190507-054

中图分类号：TS205.9；TB79

文献标志码：A

文章编号：1002-6630（2020）11-0050-07

引文格式：王宁, 阚安康, 黄孜沛, 等. 柱状茭白真空预冷热质传递分析[J]. 食品科学, 2020, 41(11): 50-56. DOI:10.7506/spkx1002-6630-20190507-054. http://www.spkx.net.cn

WANG Ning, KAN Ankang, HUANG Zipei, et al. Analysis of heat and mass transfer of cylindrical pieces of Zizania latifolia during vacuum precooling[J]. Food Science, 2020, 41(11): 50-56. (in Chinese with English abstract) DOI:10.7506/spkx1002-6630-20190507-054. http://www.spkx.net.cn