冲击式速冻机中送风孔板对鱼体冻结效率的影响.pdf
《冲击式速冻机中送风孔板对鱼体冻结效率的影响.pdf》由会员分享,可在线阅读,更多相关《冲击式速冻机中送风孔板对鱼体冻结效率的影响.pdf(9页珍藏版)》请在咨信网上搜索。
1、第 50 卷第 4 期渔 业 现 代 化Vol.50 No.42023 年 8 月FISHERY MODERNIZATIONAug.2023DOI:10.3969/j.issn.1007-9580.2023.04.003收稿日期:2022-12-21基金项目:国家重点研发计划项目(2019YFD0901802)作者简介:欧阳杰(1983),男,副研究员,研究方向:水产品加工技术与装备。E-mail:ouyangjie 通信作者:赵日晶(1989),男,副教授,研究方向:制冷与低温工程。E-mail:zhaorijing.4 冲击式速冻机中送风孔板对鱼体冻结效率的影响欧阳杰1,3,沈 建1,3,
2、郭文华2,孙鑫宇2,赵日晶2,黄 东2(1 中国水产科学研究院渔业机械仪器研究所,上海 200092;2 西安交通大学能源与动力工程学院,陕西 西安 710049;3 海洋食品精深加工关键技术省部共建协同创新中心,大连工业大学,辽宁 大连 116034)摘要:为探明冲击式速冻机中送风孔板对鱼体冻结效率的影响,以小黄鱼为研究对象,通过数值模拟和试验验证,研究了送风孔板的排布方式、送风孔板孔径和鱼体中心距上下送风孔板的距离对鱼体冻结效率的影响,分析了不同送风孔板排布方式、送风孔板孔径及送风距离下鱼体的温度分布及中心温度变化情况。结果显示:上下送风孔板交错排布避免了两侧风的对冲损耗,使得冻结速率最高
3、、冻结时间最短;送风孔板孔径的增加有效减小了第二阶段的鱼体冻结和第三阶段的鱼体降温时间;上下孔板距鱼体中心的距离越短,冻结时间越短。研究表明:筛选出适宜的送风孔板结构和参数,可为水产品冷冻提供新的工艺指导,为冲击式速冻机的优化设计提供技术参考。关键词:送风孔板;冻结特性;冲击式速冻;数值模拟;小黄鱼;水产品加工 中图分类号:S969;TK263.4 文献标志码:A 文章编号:1007-9580(2023)04-0021-009 相比于传统的速冻机,冲击式速冻机有着冻结时间短、能耗低的优点,逐渐成为速冻行业的主流设备1-2。气流喷射是冲击式速冻的关键技术,主要是因为气流通过不同类型的喷口结构后速
4、度增加,可有效缩短冻品的冻结时间3-6。喷口的结构会影响气流组织的均匀性和换热效率,实际生产中,比较常见的喷口结构主要有送风孔板、条缝喷口、圆漏斗喷口等7。国外关于不同类型的喷口结构对传热特性的影响研究较多。Bijarchi 等8研究发现喷嘴的喷射速度和对流换热速率随雷诺数的增加而增加,喷射距离的增加虽然导致了对流换热速率的减小,但可使换热较为均匀;Xu 等9 研究发现螺旋喷射有助于提升目标面的对流换热系数,且换热效率随螺旋角的增加而增加;Tepe 等10-11研究发现延伸喷射孔有助于增加带有肋片的粗糙表面和水平表面的努赛尔数(Nu)。Yang 等12研究了圆形、喷嘴和 2D 轮廓喷嘴下的喷射
5、对流换热特性,发现在喷嘴出口和自由表面之间存在一个 Nu 峰值的滞止区,圆形喷嘴和 2D 轮廓喷嘴在此达到最大湍流强度;Fenot 等13研究发现对流换热受阻塞比(Dint/D)和喷射高度(H/D)的影响。Chang 等14研究发现带有凹槽的表面有助于提升传热性能;Yadav等15研究了带有 3 种不同孔口直径下的送风孔板的传热特性,发现 Nu 随喷嘴直径的减小而增大;Ingole 等16研究了圆形倾斜射流换热特性,发现平均 Nu 取决于射流的倾角,且平均 Nu 随射流倾角的减小而减小;Ahmed 等17研究了一种新型的反向射流冲击传热,发现将喷嘴延伸至横流通道能显著降低横流的影响、提高平均
6、Nu;Youn等18研究发现优化送风孔板的排列可以提升换热效率,尤其拓展喷口可以减弱横流来提升传热效率;Makatar 等19研究了送风孔板阵列下冲击射流的传热特性,得到了关于喷射距离、喷射间距和喷射雷诺数(Re)的关系式。国内也有较多关于喷口形状、结构对冻机内部换热特性影响的研究,优化形成了相关的工艺和结构参数20-26。上述研究多集中在喷口结构对于速冻机内部流场和温度场的影响,而关于喷口结构对于冻品(尤其是鱼类)冻结效率的研究较少。渔 业 现 代 化2023 年本研究以小黄鱼为研究对象,建立了冻品内部的冻结模型,通过数值模拟和试验研究送风孔板的排布方式、送风孔板孔径和鱼体中心距上下送风孔板
7、的距离对鱼体冻结效率的影响,筛选出适宜的送风孔板结构和参数,为水产品冷冻提供新的工艺指导,为冲击式速冻机的优化设计提供技术支撑。1 材料与方法1.1 材料与仪器冰鲜小黄鱼,由浙江兴业集团有限公司提供,每条质量 80100 g,泡沫箱加冰运输至实验室。冲击式速冻机(MJN700-0818,冰轮环境技术股份有限公司);热线风速仪(TESTO-425,德图仪器);热电偶(T 型,精度0.1,开普森电子有限公司)。1.2 方法1.2.1 数值模拟1.2.1.1 物理模型本研究主要集中于送风孔板对于鱼体冻结过程的影响,不考虑速冻机的运转情况,因此只对单个鱼体进行建模,按照实际的冲击式速冻设备对模型尺寸进
8、行设计。将上下表面定义为上下两侧的送风口,四周定义为出风口。送风孔板分别置于上下送风口与鱼体之间的位置,其中上送风孔板距离鱼体中心为 h1(h1可取 80 mm、100 mm),下送风孔板距离鱼体中心为 h2(h2可取 35 mm、55 mm),如图 1a 所示。送风孔板的孔径为 D(D可取 6 mm、8 mm 和 10 mm),孔间距 26 mm。鱼体的尺寸细节如图 1b 所示。图 1 鱼体及周围流体域模型和鱼体尺寸细节图Fig.1 Fish and surrounding fluid domain model and detailed size of fish1.2.1.2 数学模型鱼体表
9、面受到高速冷空气喷射对流换热,鱼体内部伴随着凝结相变过程。假设空气为不可压缩牛顿型流体、空气的物性为常数。鱼体内部凝结过程中体积膨胀和黏性耗散忽略不计。则描述数学模型的连续性、动量和能量方程如下:连续性方程:ux+vx+wx=0(1)x 方向的动量方程:ut+(uu)x+(uv)y+(wu)z=-px+x(+t)2ux()+y(+t)uy+vx()+z(+t)uz+wx+Sx(2)y 方向的动量方程:vt+(uv)x+(vv)y+(wv)z=-py+x(+t)vx+uy()+y(+t)2vy()(+22第 4 期欧阳杰等:冲击式速冻机中送风孔板对鱼体冻结效率的影响 z(+t)vz+vx()+S
10、y(3)z 方向的动量方程:wt+(uw)x+(vw)y+(ww)z=-py+x(+t)wx+uy()+y(+t)wy+vz()+z(+t)2wz()+Sz(4)能量方程:Ht+(uH)x+(vH)y+(wH)z=xP+tT()Tx+yP+tT()Ty+zP+tT()Tz(5)k-方程:(uk)x+(vk)y+(wk)z=x+tk()kx+y+tk()ky+z+tk()kz+Gk-(6)(u)x+(v)y+(w)z=x+tk()x+y+tk()y+z+tk()z+k(c1Gk-c2)(7)式(2)(4)中,S 表示源项,定义如下27:Sx=(1-)23+AmuSy=(1-)23+AmvSz=(
11、1-)23+Amw(8)式(8)中,为液相分数,其与温度的关系如下所示,其中 Ts和 TL分别为固相和液相温度28。=0 T Ts(T-Ts)/(Tl-Ts)Ts T Tl1 T Tl(9)总焓 H 等于显热 h 和潜热 H 之和,即H=h+Hh=href+TTrefcPdTH=L(10)1.2.1.3 边界条件及模型参数上下采用速度入口条件,其中上入口的速度为 2 m/s,下入口速度为 1 m/s,送风温度为236.15 K。其中上下送风风速是根据风机风量和风口面积确定的,此冲击式速冻机中每台风机的风量为 12 000 m3/h。四周采用自由出流条件。鱼体的初始温度为 292.15 K。模拟
12、在瞬态条件下选用有限体积法求解,时间步长为 0.5 s。采用焓-多孔度方法29模拟鱼体冻结过程和追踪固液相面,固液糊状区常数 Am取值为 105 kg/(m3s)。利用 PISO 算法耦合速度和压力,利用 PRESTO 算法处理压力修正方程,动量和能量方程通过 QUICK 差分法处理。动量、能量和压力关联方程的松弛因子分别为 0.7、0.5、0.3。鱼体的物性参数为:=998 kg/m、=8.910-4 Pa、Cps|Cpl=2 100 4 179 J/(kgK)、ks|kl=2.2|0.6 W/(mK)、Ts|TL=271.65 K、L=335 000 J/kg、=2.110-4 K-1。为
13、研究网格数和时间步长对计算精度的影响,选取 4 种不同网格数(包括 100 000、200 000、400 000 和 800 000)的结构化六面体网格和 3 种不同的时间步长(0.1、0.5 和 1.0 s)进行分析。表 1 给出了上下送风孔板对称排布,孔径为 8 mm时,不同网格数和不同时间步长情况下计算的冻结时间。可以看出,网格数由 100 000 增加到400 000,冻结时间偏差从 7.6%降低到 1.5%;时间步长由 1 s 减小到 0.1 s,冻结时间偏差由32渔 业 现 代 化2023 年3.8%降低到 0.8%;故同时考虑计算精度和速度,网格数采用 400 000、时间步长
14、采用 0.5 s。表 1 网格数和时间步长对计算精度的影响Tab.1 Effects of grid number and time step on accuracy ofnumerical simulation网格数冻结时间/s偏差/%时间步长/s冻结时间/s偏差/%100 0006107.61.06303.8200 0006353.80.56500.8400 0006501.50.1655基准800 000660基准1.2.2 试验验证设计上下冲击式试验台,测试了小黄鱼中心温度随时间的变化,如图 2 所示。该试验台主要由离心风机、静压箱、上下送风孔板和网带组成。冷空气被离心风机吸入静压箱内
15、,通过上下孔板后喷射到小黄鱼的表面。上送风孔板距离小黄鱼50 mm,下送风孔板距离小黄鱼 20 mm,且上下送风孔板对称。利用热线风速仪测量距离上孔板120 mm 和距离下孔板 60 mm 处的风速。数值模拟中上下风速设置与其保持一致。利用热电偶测量小黄鱼的中心温度,利用温度采集仪检测鱼体中心温度随着时间的变化。待鱼体中心降到-18时,视为冻结完成。图 2 上下冲击式试验台及小黄鱼冻品Fig.2 Impingement quick-freezing test and yellow fish 小黄鱼中心温度随时间变化的试验和模拟值的对比如图 3 所示。KTQ#,图 3 小黄鱼中心温度随时间变化的
16、试验和模拟值的对比 Fig 3 Comparison of fish center temperature with time between numerical simulation and experimental results 图 3 对比了试验测试和数值模拟的小黄鱼中心温度随时间的变化,可以看出试验测试小黄鱼的冻结时间为 649 s,数值模拟得到小黄鱼的冻结时间为 600 s,模拟相对偏差为 7.55%,验证了数值模拟方法的准确性。2 结果与讨论2.1 不同送风孔板的排布方式对冻结特性的影响送风孔板的排布方式对鱼体的冻结特性有较大的影响,研究了只有上送风孔板、上下送风孔板对称和上下送
17、风孔板交错 3 种不同的排布方式下的鱼体冻结特性。图 4 表示 3 种不同排布方式下鱼体中心温度随时间的变化。从图 4 中可以看出,鱼体的速冻分为 3 个阶段,第 1 阶段是初温到冰点的显热释放阶段(初温271.65 K)。第 2 阶段是鱼体从冰点降到中心温度为 268.15 K,此过42第 4 期欧阳杰等:冲击式速冻机中送风孔板对鱼体冻结效率的影响程80%水分被冻结成冰晶。第3 阶段是鱼体中心温度由 268.15 K 降到 255.15 K 30。只有上送风孔板、上下孔板对称和上下孔板交错布置时鱼体冻结(从初始温度 278.15 K 冷冻到 255.15 K)所需的时间分别为 800 s、6
18、50 s 和 540 s。上下孔板交错布置时的速冻时间最短,比只有上送风孔板和上 下 孔 板 对 称 布 置 时 减 少 了 32.50%和16.92%。上下孔板交错排布时冻结时间的减小主要体现在第 2 阶段,鱼体内部水分冻结为冰晶的时间缩短,这一现象可通过鱼体内部的液相分数随时间的变化来表征。KTQ#,图 4 3 种不同排布方式下鱼体中心温度随时间的变化Fig.4 Evolution of fish center temperature with time for three different orifice plate arrangement图 5 表示 3 种不同排布方式下鱼体内部液相
19、分数随时间的变化。液相分数为 1.0 时表示鱼体内部的水分没有冻结生成冰晶,为 0.0 时表示水分全部冻结成为冰晶。鱼体中心温度由冰点降到268.15 K,只有上送风孔板、上下孔板对称和上下孔板交错排布所需的时间分别为 643 s、509 s 和402 s。结果表明上下孔板交错布置有助于鱼体快速通过最大冰晶生成阶段。1.00.80.60.40.20.0KT,图 5 3 种不同排布方式下鱼体内部液相分数随时间的变化Fig.5 Evolution of fish liquid fraction with time for three different orifice plate arrangem
20、ent图 6 为 3 种不同排布方式下鱼体温度和液相分数云图,选取的时刻分别为 90 s、240 s、390 s 和540 s。从图中可以看出 3 种排布方式下鱼体的冻结都是先由左右两端到中心,然后转变为从四周到中心。相比于只有上送风和上下对称排布时,上下交错排布的转变较快。由于上下交错送风避免了风的对冲损耗,从而冻结速率加快。/K图 6 3 种不同排布方式下液相分数和鱼体温度分布云图Fig.6 Liquid fraction and fish temperature colour for three different orifice plate arrangement52渔 业 现 代 化
21、2023 年 图 7 为 3 种不同排布方式下的速度矢量图和迹线图。从速度矢量图可以看出,冷风经过孔板加速后直接喷射到鱼体表面。只有上送风孔板时,鱼体的上表面受到了冷空气的喷射,冷量只从上表面传递到下表面,导致了冻结时间较长。在送风孔板上下对称布置时,鱼体的上下表面都受到了冷空气的喷射,冷量是从鱼体上下表面到中心传递,冻结速率加快。但上下孔板的对称布置导致了冷风的对冲,鱼体的上表面出现了涡旋,导致了冷量的损耗31-32。上下交错布置避免了冷风的对冲,充分利用了上下孔板的喷射作用,有效地提升冻结速率、减小冻结时间。EME-FE3/J7DMPDJUZ$POUPVS图 7 3 种不同排布下的速度矢量
22、图和迹线图Fig.7 Vector and streamline colour for three different orifice plate arrangement2.2 不同送风孔板孔径对冻结效率的影响从上面的分析可得,上下孔板交错排布时的冻结效率最高、冻结时间最短。因此选择此种排布形式下,研究不同送风孔板孔径(6 mm、8 mm和 10 mm)对冻结特性的影响。图 8 表示了不同送风孔板孔径下鱼体中心温度和液相分数随时间的变化。从图 9 a 可以看出,鱼体冻结时间(初始温度 278.15 K 冷冻到 255.15 K)随孔板孔径的增加而增加。孔板孔径为 6 mm、8 mm 和 10
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 冲击 速冻 送风 冻结 效率 影响
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。