牛肉胴体预冷库温度时空分布模拟研究

　　摘要：该文应用计算流体动力学（computational fluid dynamics，CFD）技术，建立数学模型求解牛胴体预冷库中的温度场，模拟研究牛胴体预冷过程中温度的时空分布。结果表明，牛胴体上部分比下部分降温快，正对冷风机的牛胴体降温速率明显高于非正对冷风机的牛胴体；每个胴体平均温度都低于7℃，时间为44h，胴体中心最高温度低于7℃，需预冷72h左右；牛胴体温度离散系数呈上凸曲线变化，峰值为0.5，且前两排牛胴体温度离散系数比后两排先18h达到峰值。模拟结果可为合理掌控预冷时间减少能源消耗、提升预冷效果改善牛肉品质提供参考，同时也为不同预冷需求的预冷库结构优化、预冷设备合理分布提供理论参考与借鉴。
 

　　关键词：计算流体动力学；预冷；牛胴体；温度场；牛肉品质
 

　　我国是牛肉生产和消费大国，2019年我国牛肉产量仅次于美国、巴西和欧盟；牛肉消费量位于世界第二位，仅次于美国。随着人们生活水平的提高，牛肉需求量越来越大的同时，大家对牛肉品质要求也随之提高。牛胴体预冷排酸对牛肉品质有重要意义。牛被宰杀后体温升高，胴体中心温度可达40℃，成为滋生细菌的温床，低温预冷能够很好地抑制细菌生长繁殖，保持牛肉品质。另外，牛在宰杀的时候，因为恐惧等因素，导致肉质坚硬，干躁，缺少弹性，且肌肉中的糖元无氧酵解产生乳酸，不仅影响牛肉品质，还不利于人体健康。低温预冷排酸，可将乳酸分解成二氧化碳、水和酒精然后挥发。由于低温预冷排酸经历了较为充分的解僵、成熟过程，不仅肉质柔软有弹性、味道鲜美，而且安全营养。因此，牛宰后胴体预冷排酸是生产高品质牛肉的第一个环节也是必要环节，通过预冷改善牛肉品质成为研究热点。
 

　　然而，国内预冷库大多是根据经验来设计，导致预冷阶段牛胴体温度分布不均匀，浪费能源的同时还影响牛肉品质。利用数值模拟技术研究牛肉胴体预冷库温度时空分布，使温度场在预冷库中分布更加均匀，对改善牛肉品质具有重要意义。计算流体动力学（computational fluid dynamics，CFD）是基于计算机技术的一种数值计算工具，用于求解流体的流动和换热问题。CFD技术被国内外研究者广泛应用于研究低温冷却过程中的温度场分布。
 

　　文献分别模拟了苹果、番茄、蓝莓、茎状蔬菜等冷却时的温度变化情况，为改善果蔬冷却效果提供了参考。刘泽勤等以压差原理的果蔬仓库为研究对象，对稳态下的库内温度分布情况进行模拟仿真，讨论了果蔬间距的改变对库内果蔬温湿度的影响。赵春江等建立了求解短距离运输的冷藏车车厢内温度场分布的计算模型，模拟分析了不同边界条件和货物不同堆栈方式对车厢内温度场分布的影响。赵时等建立了冷藏车厢的仿真模型，利用CFD技术研究梯级送风对空仓时冷藏车厢内温度场的影响。国内虽然针对果蔬包装箱、冷藏车、冷库中温度场的研究比较多，但目前利用CFD技术对牛胴体预冷库研究的很少。而国外对于牛胴体预冷库研究大多是一次预冷一个或几个牛胴体，并没有研究实际上的一次预冷几十上百个胴体，无法体现预冷时温度场的均匀性，在提高牛胴体预冷效果方面缺少实用价值。为此，本文以吊挂有76个牛胴体二分体的实际预冷库为研究对象，利用CFD技术进行模拟计算，分析不同预冷时间不同位置的温度分布，对比不同时间不同位置的降温快慢、冷却均匀性，客观评价预冷效果。
 

　　1 材料与方法

　　1.1 物理模型

　　预冷库物理模型根据实地测量北京某畜牧有限公司牛肉胴体预冷库得出。牛胴体二分体模型参考实际尺寸并做适当简化，应用3dmax进行三维建模绘制。预冷库空间尺寸（长×宽×高）为1.6m×5.25m×3.8m，4个冷风机空间尺寸（长×宽×高）都为1.6m×0.6m×0.65m，风机距天花板0.35m，距墙壁0.6m。预冷库进风口为圆形，位于冷风机前部，直径0.5m，风速6m/s，温度－5℃。预冷库出风口为圆形，位于冷风机背面，直径0.5m。本次模拟实验一次性预冷76个二分体。4根吊挂轨道，每根吊轨吊挂19个二分体。吊轨间距1m，相邻吊钩间距0.8m。牛胴体二分体长2.3m、宽0.75m、厚0.35m。预冷库三维图（省略了吊挂牛胴体的吊钩）及牛胴体二分体模型如图1所示。

图1 预冷库结构示意图

　　1.2 网格划分

　　本文采用软件ICEM CFD 15.0进行网格划分，整个预冷库及内部结构采用六面体非结构化网格进行划分。全局网格最大尺寸不超过0.2m，壁面及牛胴体最大网格尺寸不超过0.1m。网格数量为1 621 719个单元，划分结果如图2所示。非结构网格如果质量太差，在Fluent软件计算过程中容易导致计算发散或者收敛困难。网格质量正交性指标分布从0到1，越接近0网格质量越差。利用mesh quality对网格的正交性进行质量检测，结果表明网格质量均大于0.35，平均网格质量为0.77，表示网格质量良好。

a-预冷库整体网络；b-牛胴体网络

图2 网格划分

　　1.3 数学模型

　　建立数学模型求解牛胴体二分体预冷过程的温度场分布情况时，为了达到预期效果的同时尽量减少计算时间，研究对实际模型进行简化，做出如下假设：空气为不可压缩气体且符合Boussinesq假设；预冷库内管道、铁架等对流场的影响忽略不计；冷风机内部温度场变化对整个预冷库内温度场无影响；忽略墙壁与外界的热交换；全部牛二分体胴体是同时进入预冷库的；牛胴体的物理参数不随温度的变化而改变，忽略牛胴体预冷时的重量损失。

　　预冷库中空气、牛胴体的具体热物理性能参数如表1所示。

 

　　图5显示的是不同时间下胴体表面的温度分布云图，可以看出，正对风机的牛胴体，降温速率明显高于非正对风机的牛胴体。在条件允许范围内，建议定时改变风机吹风方向，或者定时调换正对风机与非正对风机牛胴体的位置，使牛胴体冷却更均匀。

　　2.2 预冷时间

　　依据《牛羊屠宰与分割车间设计规范》（GB 51225-2017），牛胴体冷却后中心温度不应高于7℃。所以当胴体中心平均温度达到7℃时，胴体基本达到预期预冷温度，所用时间为预冷时间。图6是所有牛胴体表面和胴体内部平均温度随预冷时间变化的曲线。由图6可知，所有胴体表面平均温度达到7℃需要18h。所有胴体内部平均温度达到7℃，需要40h。通过2.1节图3的温度云图可知，要使胴体中心最高温度低于7℃，则需要72h左右。

 

　　图7是预冷库中平均温度低于7℃的牛胴体个数随预冷时间增长而增加的曲线。从图7可看出，预冷28h后开始有1个胴体平均温度低于7℃，从30h到44h，平均温度低于7℃胴体个数呈线性增长趋势。38h后超过一半胴体平均温度低于7℃，这是通常所说的二分之一预冷时间。42h后有超过八分之七的胴体平均温度低于7℃，这是通常所说的八分之七预冷时间。44h后每个胴体平均温度都低于7℃。因为28h后陆续有胴体平均温度低于7℃，建议在28h后，将先预冷好的牛胴体陆续移出预冷库进行下一步加工处理。
 

　　2.3 冷却均匀性

　　温度离散系数用来评估冷却均匀性，温度离散系数越大，冷库内牛胴体温度数据离散程度越大，温度分布越不均匀；反之，温度离散系数越小，冷库内牛胴体温度数据离散程度越小，温度分布越均匀。采用热力学温度计算温度离散系数，计算如公式（1）所示：

　　式中：Cv为温度离散系数，%；n为牛胴体数量；Tave为n个牛胴体的平均温度；Ti为牛胴体i的平均温度。
 

　　以离冷风机由近到远分成第1排、第2排、第3排、第4排牛胴体，图8是4排牛胴体和所有牛胴体的温度离散系数曲线图。从图8可以看出，各排和整体温度离散系数都是先升高再降低。当温度离散系数达到峰值后，随着冷却时间的增加，预冷库内胴体间温度差异性越来越小，温度离散系数越来越小。整体温度离散系数峰值为0.5，各排牛胴体温度离散系数变化趋势相近，但后两排比前两排晚18h达到峰值。前两排在12h达到峰值，后两排在30h分别达到峰值。这可能是因为前两排牛胴体离冷风机近，牛胴体降温相对快于后两排，随着温度的降低，前两排胴体间温度差异先减小，温度离散系数先降低。

　　3 结论

　　风机朝一个方向吹风，导致牛胴体上部分比下部分降温快，非正对冷风机的牛胴体冷却效果不如正对冷风机的牛胴体。在条件允许范围内，建议改进风机，让风机朝不同方向发散吹风或定时改变风机吹风方向，或者定时调换正对风机与非正对风机牛胴体的位置，减少牛胴体温度离散系数，使牛胴体冷却更均匀。
 

　　冷却28h后，开始有牛胴体平均温度低于7℃；44h后，每个胴体平均温度都低于7℃；要使胴体中心最高温度低于7℃，则需要72h左右。建议在28h后，将先预冷好的牛胴体陆续移出预冷库进行下一步加工处理。
 

　　对于不同冷库大小、不同初始温度、不同预冷胴体大小及个数，可参考本研究，适当选择风机（包括风机个数、风速、送风温度等），适当调整预冷时间。