1、间冷式冰箱空气循环系统能量流动特性研究王坦1， 彭玲2， 向东1， H.F.TENG1（1.清华大学 机械工程学院精密超精密制造装备及控制北京市重点实验室， 北京 100084；2.四川长虹电器股份有限公司， 四川 绵阳621000 ）来稿日期： 2018-03-17基金项目： 国家高技术研究发展计划 （863 计划 ） （2013AA041305 ）作者简介： 王坦，（1992- ） ， 男， 山东人， 硕士研究生， 主要研究方向： 机电产品绿色设计理论及应用、 模块化设计；彭玲，（1971- ） ， 女， 四川人， 高级工程师， 主要研究方向： 家电产品绿色制造、 新材料及工业工程领域1
2、引言空气循环系统的能量消耗在冰箱总能耗中占有较大比例，对冰箱优化设计， 需要对空气循环系统进行研究分析。目前对于冰箱空气循环系统的研究主要集中在对冰箱腔室及风道的仿真分析方面。文献1用 FLUENT 仿真了网状和平板状搁物架对冰箱冷藏室内的温度场和流场的影响。文献2运用 CFD分析手段分析优化了某款风冷冰箱风道系统的风扇及风道结构，提升了风道的性能。文献3通过对冰箱冷藏室进行 CFD 建模仿真， 提高了冰箱冷藏室内部温度分布的均匀性。 文献4对冰箱风道的流场进行仿真， 改善了风道系统流场分布， 提高了冰箱出风口量， 有效降低了噪音。文献5对风道进行了数值模拟及相关实验，优化风道结构参数。文献6
3、通过二维数值仿真的方法对冰箱自然对流换热进行了研究。文献7通过数值仿真的方法得出了冰箱箱内空气传热的规律。 文献8-9采用了有限元仿真与实验验证相结合的方法， 对冰箱内部的空气流场和温度场进行了研究。文献10通过计算仿真的方法， 对冰箱保温层及内部空气的温度分布和导热问题进行了深入研究。 以上建模方法对于冰箱腔室及风道的仿真分析相对独立， 且关注点主要集中在流场和温度场的分析， 并不直接反映空气循环系统的能量流动和消耗。为此， 首先将间冷式冰箱空气循环系统划分为多个模块， 利用 CFD 仿真方法求出每个模块的接口参数，并通过计算确定空气循环系统能量流动的分布规律， 为系统节能和性能优化提供理论
4、依据。2 空气循环系统对于间冷式冰箱而言，空气在循环系统中按照蒸发器腔送风风道冷藏室回风风道蒸发器腔的方向循环流动。为了直观描述间冷式冰箱空气循环系统的能量流动特性，提出了如图 1 所示的空气循环系统建模流程， 该流程首先对空气摘要： 能量在间冷式空气循环系统中的流动规律对产品的能耗与性能具有重要的影响。为此， 首先针对空气循环系统流固耦合的特征， 提出了体积元法和局部划分法相结合的模块划分方法； 在此基础上， 利用 CFD 方法仿真计算得到相应的接口参数； 然后， 运用数值计算得出系统中各部分间的能量流动特性。最后以某型号间冷式冰箱空气循环系统为例， 得出了其能量流动特性并验证了这里方法的有
5、效性。关键词： 间冷式冰箱； 空气循环系统； CFD 方法； 能量流动特性中图分类号： TH16； TH122文献标识码： A文章编号： 1001-3997 （2018 ） 08-0001-04Research on Energy Flow Characteristics of Air CirculationSystem of Indirect-Cool RefrigeratorWANG Tan1， PENG Ling2， XIANG Dong1， H.F.TENG1（1.Beijing Key Lab of Precision/Ultra-Precision Manufacturing Eq
6、uipments and Control， Tsinghua University， Beijing 100084，China； 2.Sichuan Changhong Electroic Co. Ltd， Sichuan Mianyang 621000， China ）Abstract： The energy flow routine in each part of air circulation system plays a key role in affecting energy consumption andthe performance of indirect-cool refrig
7、erator. Therefore，the modeling approach，both the volume element method and localdivision method，is proposed. On this basis，CFD method was employed to calculate the interface parameters of the system.Then，the numerical calculation method was applied to exhibit energy flow characteristics among each p
8、art of the system.Finally，an instance of a certain type of indirect-cool refrigerator air circulation system is given to draw a description to theenergy flow characteristics， which is synchronously a solid proof to the validation of method in this research.Key Words： Indirect-Cool Refrigerator； Air
9、Circulation System； CFD Method； Energy Flow CharacteristicsMachineryDesignManufacture机械设计与制造第 8 期2018 年 8 月1ChaoXing循环系统进行模块划分；然后运用 CFD 仿真方法计算出各个模块间的接口参数； 最后， 求解系统各模块的能量流动特性。 该流程主要环节的详细介绍， 如图 1 所示。模块划分CFD 方法求解接口参数计算能量流动体积元划分法局部划分法图 1 冰箱空气循环系统建模流程Fig.1 The Modeling Processes of Air Circulation System
10、3 系统建模根据间冷式冰箱空气循环系统流固耦合的特性，需要对其进行模块划分。3.1 基于体积元法的流体模块划分在间冷式冰箱的空气循环系统中， 空气在蒸发器腔和冷藏室循环流动， 在蒸发器处放出热量， 从保存的物品和隔热层中吸收热量。流体在循环流动时， 温度、 压力、 流速等物理特性会随着流动发生较大范围的改变。为解决整体描述中误差较大的情况，提出体积元划分法。冷藏室内由于搁架的存在， 上下两个部分相互之间特性参数差异较大， 所以以搁架及其延伸面为分界线， 将冷藏室内的空气分为上下两个模块。而送风风道和回风风道内的流体区域， 根据其与周围零部件的接口特征， 分为位于冷藏室后部隔热层内的上送风风道内
11、空气模块和上回风风道内空气模块， 以及位于蒸发器腔后部隔热层内的下送风风道内空气模块和下回风风道内空气模块。具体划分结果， 如图 2、 表 1 所示。送风风道上冷藏室下冷藏室上风道下风道蒸发器腔搁架冷藏室图 2 冰箱内空气流通区域的模块划分Fig.2 Module Division of Air Circulation Area表 1 体积元划分法划分结果Tab.1 The Results of Volume Element Method编号名称编号名称E1蒸发器腔空气E5冷藏室下层空气E2下送风风道内空气E6上回风风道内空气E3上送风风道内空气E7下回风风道内空气E4冷藏室上层空气3.2 基
12、于局部划分法的固体模块划分空气循环系统的风道和箱体是由固态的隔热层组成的， 同样存在体积大、 特性参数差异大的现象。 对于这种情况， 根据隔热层的位置、 厚度、 换热面积、 导热系数， 将整个隔热层划分为局部特性一致的多个模块。具体划分结果， 如表 2 所示。表 2 局部划分法划分结果Tab.2 The Result of Local Classification Method编号名称编号名称E8蒸发器腔后壁面E14冷藏室后壁面E9蒸发器腔前壁面E15冷藏室前壁面E10蒸发器腔左侧壁面E16冷藏室左侧壁面E11蒸发器腔右侧壁面E17冷藏室右侧壁面E12蒸发器腔顶部壁面E18冷藏室顶部壁面E13
13、蒸发器腔底部壁面E19冷藏室底部壁面采用体积元划分法和局部划分法， 间冷式冰箱空气循环系统被划分为表 1 和表 2 所示的 19 个模块， 各模块功能明确， 特性参数和接口参数形式简单、 易获取， 而且模块之间的能量流动易于表达和计算， 从而有利于用能量流动的方式更加准确地表达出系统消耗能量实现功能/性能的过程。3.3 接口参数计算为了计算空气循环系统能量流动大小， 需要首先获取以上划分模块的特性参数和接口参数。 其中一部分特性参数如结构参数是已知的， 还有一部分， 如空气循环系统中流场的流量、 温度、压强等， 需要通过 CFD 仿真方法计算获得。3.4 计算模块间的能量流动根据各模块间的接口
14、特性，可以按照下面几类能量流动形式计算能量流动大小。3.4.1 空气模块之间的能量流动在间冷式冰箱工作时， 空气在模块之间循环流动， 当空气从一个模块流出进入另一个模块时， 空气所携带的机械能和内能从前一个模块输出， 输入了后一个模块。其形式可由式 （1 ） 表示。EM=f （m， T， p， u ） =mp+mu22+U（1 ）式中： m流量； T温度； p压强； u流速； U内能； 空气的压力能； mu22空气的运动能。3.4.2 空气模块- 隔热层模块- 环境空气间的能量流动空气循环系统内的空气、 环境空气与壁面的换热过程属于受迫对流换热， 如图 3 所示。考虑如下条件：（1 ） 只研究
15、稳态状态下的温度场， 忽略时间项的影响；（2 ） 环境空气温度恒定， 边界条件为第一类边界条件；（3 ） 空气为牛顿流体， 湿度为常量；（4 ） 箱体内空气流动形式为稳定层流和非边界层型流动；（5 ） 空气在壁面上满足无滑移边界条件；（6 ） 忽略空气的粘性耗散。空气模块-隔热层模块-环境空气间热量交换的能量流动可以用式 （2 ） 计算：内壁面与箱内空气的对流传热系数 1Tr环境温度 Te外壁面与箱外空气的对流传热系数 2壁厚 h， 壁面材料导热系数 图 3 壁面传热模型Fig.3 The Wall Heat Transfer Model第 8 期王坦等： 间冷式冰箱空气循环系统能量流动特性研
16、究2ChaoXingEc=Q=F k T （W ）（2 ）其中， k111+12+h（W/m2 K ）式中： F换热面积；T箱内空气温度 Tr和环境温度 Te之差；1、 2内外壁面与空气的对流换热系数；h壁面厚度；壁内热传导系数。3.4.3 冷藏室空气模块与搁架之间的能量流动在冷藏室内， 搁架将腔内空间分为不同的空气模块。在稳态下， 空气温度与搁架温度均保持恒定， 故搁架输入输出的能量流动在数值上一致。 空气模块与搁架之间的能量流动的计算方式与空气模块与壁面的能量流动计算方式类似。3.4.4 蒸发器腔空气模块与蒸发器之间的能量流动在蒸发器腔内，能量以强制对流换热的方式从蒸发器腔内空气流入蒸发器
17、。 冰箱内部气体模块与蒸发器的换热过程属于受迫对流换热。基于 3.4.2 中的假设， 可利用公式 （1 ） 进行能量流动量的计算。除去以上几种能量流动形式外， 一些能量流动如辐射换热等由于数值较小不予考虑。运用以上方法， 可以计算空气循环系统各模块间能量流动的方向及大小。下面以具体某型号冰箱为例。3.5 各部分间的能量流动分布根据以上模块划分方法，间冷式冰箱空气循环系统能量流动方式， 如图 4 所示。 其中， a 箭头表示空气间能量流动， b 箭头表示环境与壁面间的能量流动， c 箭头表示固体模块与空气模块间的能量流动。Ec， in， 18EM3Ec， in， 14EM5Ec， 14， 5Ec
18、， 14， 4Ec， 20， 4Ec， 18， 4EM4Ec， 15， 4Ec， in， 15Ec， 5， 20Ec， 19， 4Ec， in， 19EM1Ec， in， 8Ec， 8， 1Ec， I， evaEc， in， 9Ec， 9， IEM7Ec， in， 17Ec， 17， 4EM2Ec， in， 12EM6Ec， 12， 1Ec， in， 10Ec， 10， 1Ec， 13， 1Ec， in， 13Ec， 11， 1Ec， in， 11Ec， 16， 4Ec， in， 16bbbabaaabbabaabbabbb图 4 间冷式冰箱空气循环系统能量流动Fig.4 The Energy
19、 Flow of Air Circulation System4 实验验证4.1 参数设置使用某型号的间冷式冰箱作为研究对象。其部分结构参数，如表 3 所示。4.2 接口参数计算与实验测试为了计算空气循环系统各模块的接口参数， 设定：（1 ） 以空气循环系统的整个工作周期为研究对象， 可以认为系统处于稳态， 各腔室、 蒸发器以及内部的温度场不随时间变化。（2 ） 环境温度为 25， 蒸发器温度为-40， 空气质量流量为0.03kg/s。基于所设定的边界条件， 利用 Fluent 软件进行仿真计算， 得出空气流场仿真以及温度场分布情况， 其中， 温度场分布情况， 如图 5 所示。表 3 冰箱结构
20、参数Tab.3 The Structural Parameters of Refrigerator序号模块名称长宽（mm ）换热面积（mm2）厚度（mm ）密度（kg/m3）导热系数（W/m ）1蒸发器腔后壁面40739015873084300.022蒸发器腔前壁面40739015873060300.023蒸发器腔左侧壁面3907027300116300.024蒸发器腔右侧壁面3907027300116300.025蒸发器腔顶部壁面407702849080300.026蒸发器腔底部壁面407702849036.5300.027冷藏室后壁面540270.5146070167.5300.028冷藏
21、室前壁面540270.5146070100300.029冷藏室左侧壁面445270.5 120372.553300.0210冷藏室右侧壁面445270.5 120372.553300.0211冷藏室顶部壁面54044524030036.5300.0212冷藏室底部壁面54044524030060300.02图 5 空气温度场的仿真结果Fig.5 Simulation Result of the Ai Temperature Field为判断仿真得出的接口参数的准确性， 在蒸发器附近等便于测量的位置分别设置温度传感器， 利用冰箱测试系统及软件采集稳态下的温度数据， 如图 6 所示。并与仿真结果中
22、对应处的温度进行对比， 如表 4 所示。图 6 冰箱样品测试平台接口及软件Fig.6 The Interface of Refrigerator Test Software表 4 温度仿真结果与测试结果对比Tab.4 The Comparison of TemperatureSimulation and Test Result位置仿真结果 （K ）测试结果 （K ）误差%蒸发器腔出口242.5241.80.29蒸发器腔入口244.7245.7-0.41冷藏室入口243.1242.40.29冷藏室出口244.4245.2-0.33冷藏室搁架处244.1244.00.04通过对比仿真结果与测试结果
23、，各点温度的仿真结果与测试结果的误差均在 0.41%以下。说明 CFD 仿真方法能够对冰箱内温度场进行较为准确的预测， 得出的间冷式冰箱空气循环系统的接口参数是可信的。机械设计与制造No.8Aug.20183ChaoXing4.3 模块间能量流动计算与分析根据系统中各部分的特性参数和上节仿真得出接口参数， 可对空气循环系统中各模块间的能量流动进行计算。例如， 由空气流场仿真结果得到各模块的接口参数， 如表5所示。可计算蒸发器腔空气模块到下送风风道空气模块的能量流动EM1。已知空气流量为0.03kg/s， 由表 5 中数据可知， 蒸发器腔出口处温度 T=242.5K、 压力p=101325Pa、
24、 空气密度为 =1.4572kg/m2、 流速为 u=7.1608445m/s，则由式 （1 ） 可得， 该处的能量流动大小为EM1=11053.92W。表 5 仿真得到的接口参数Tab.5 The Interface Parameters from Simulation位置温度 （T ）压力 （p ）密度 （ ）流速 （u ）蒸发器腔入口244.71015751.457.21蒸发器腔出口242.51013251.467.16冷藏室入口243.21025251.477.10冷藏室出口244.41018251.457.18冷藏室搁架处244.11020251.460.40若计算蒸发器腔前壁面模块
25、处的能量流动 Ec， 9， 1， 已知换热面积 F=0.15873m2， 环境温度 T0=298.15K， 蒸发器腔温度 T1=243.6K，壁面厚度 h=0.06m， 导热系数 =0.05， 则由式 （2 ） 可计算出该处的能量流动大小为 Ec， 9， 1=7.22W。经过计算， 可获得的各模块间的能量流动值， 如表 6 所示。结合图 5 所示的能量流动示意图， 可以直观地描述空气循环系统在完成功能时的能量流动方向与大小。表 6 能量流动计算结果Tab.6 The Calculation Results of Energy Flow能量流动值 （W ）能量流动值 （W ）能量流动值 （W ）
26、EM111053.92Ec， 12， 12.13Ec， 16， 42.74EM211062.68Ec， 13， 10.97Ec， 17， 42.74EM311073.63Ec， 8， 10.52Ec， 18， 415.93EM411102.26Ec， 8， 20.64Ec， 19， 415.75EM511111.37Ec， 8， 72.16Ec， 20， 40.51EM611115.34Ec， in， 83.32Ec， 5， 200.51EM711121.12Ec， in， 144.99Ec， 14， 51.17Ec， 1， eva78.14Ec， 14， 62.16Ec， 15， 51.71
27、Ec， 9， 17.22Ec， 14， 32.63Ec， 16， 52.71Ec， 10， 10.64Ec， 14， 41.19Ec， 17， 52.71Ec， 11， 10.64Ec， 15， 41.73表 6 显示蒸发器腔空气与外界环境的能量交换 Ec， 1， eva、 冷藏室顶部壁面与冷藏室空气的能量交换 Ec， 18， 4以及冷藏室底部壁面与冷藏室空气的能量交换 Ec， 19， 4相对较多，可加强相应模块的保温设计。5 结论对冰箱空气循环系统进行能量建模是实现冰箱性能设计优化的基础。为了得出冰箱空气循环系统中各部分间的能量流动，运用模块划分、 CFD 仿真及数值计算等方法提出了一种有效
28、的建模方法。利用所提出的能量模型， 对某款冰箱进行建模分析， 得出其蒸发器外壁、 冷藏室上下避免处能量交换为 78.14W、 15.93W 及15.75W， 能量交换相对较大， 需要对其进行进一步的保温设计。参考文献1吴小华， 张璟， 宋春节冰箱冷藏室温度场和流场的仿真与优化 J 北京石油化工学院学报， 2006， 14 （3 ） ： 8-11（Wu Xiao-hua， Zhang Jing， Song Chun-jie.Simulation of TemperatureField and Flow Field Inside Refrigerator and Its Structure Opt
29、imizationJ .Journal of Beijing Institute of Petro-Chemical Technology， 2006， 14（3 ） ： 8-11 ）2钟明， 王瑶， 王冬祥基于 CFD 的冰箱风道优化研究 J 轻工标准与质量， 2015 （4 ） ： 53-57（Zhong Ming， Wang Yao， Wang Dong-xiang.Study on Optimization ofRefrigerator Wind Tunnel Based on CFD J .Standard & Quality ofLight Industry， 2015 （4 ）
30、： 53-57 ）3白连社， 王瑶， 关鹏基于 CFD 的风冷冰箱风道系统研究 C .中国家用电器技术大会论文集， 2012 （S1 ） ： 122-126（Bai Lian-she， Wang Yao， Guan Peng.Researches on the Air Duct Systemof the Air-cooled Refrigerator Based on CFD C .Electrical Equipment，2012 （S1 ） ： 122-126 ）4刘永辉， 王晶， 李标基于 CFD 技术的大容量风冷冰箱流场分析及改进设计 J 流体机械， 2014， 42 （7 ） ： 6
31、9-73（Liu Yong -hui， Wang Jing， Li Biao.Numerical Analysis of Air FlowField in Air-cooled Refrigerator with Large Capacity and its DesignImprovement J .Fluid Machinery， 2014， 42 （7 ） ： 69-73 ）5陈成， 刘志峰， 罗建基于 Fluent 的风冷冰箱的风道结构优化设计 J 日用电器， 2013 （5 ） ： 30-35（Chem Cheng， Liu Zhi-feng， Luo Jian.The Optimiz
32、ation Design of Air-cooled Refrigerator Duct Based on Fluent J .Electrical Appliances， 2013（5 ） ： 30-35 ）6凌长明， 陶文铨.冰箱内非稳态自然对流的二维数值模拟 J .西安交通大学学报， 1995 （10 ） ： 35-41.（Ling Chang-ming， Tao Wen-quan.Two-Dimensional Numerical Simul-ation of the Transient Natural Convection in Refrigerator J .Journal ofX
33、i an Jiaotong University， 1995 （10 ） ： 35-41. ）7周湘江， 凌长明.家用冰箱内传热的非稳态数值研究 J .中南工学院科技通讯， 1999， 15 （1 ） ： 21-24.（Zhou Xiang-jiang， Ling Chang-mingJ .Numerical Study of UnsteadyHeat Transfer in the Fridgerator， 1999， 15 （1 ） ： 21-24. ）8丁国良， Oellrich L. R. 冰箱箱内空气温度场与流场的优化研究 J .制冷学报， 1998 （1 ） ： 24-29.（Di
34、ng Guo-liang， Oellrich L.R.Researches on the Optimization of Air Flowand Temperature Field inside a Refrigerated Chamber J .RefrigerationJournal， 1998 （1 ） ： 24-29. ）9丁国良， Oellrich L. R. 冰箱中空气流场和温度场的有限元模拟 J .上海交通大学学报， 1998 （07 ） ： 20-24.（Ding Guo-liang， Oellrich L.R.Simulation of Air Flow and Temper
35、atureField inside a Refrigerated Chamber with Finite Element MethodJ .Journal of Shanghai Jiaotong University， 1998 （7 ） ： 20-24. ）10俞炳丰， 葛军， 王志刚.冷冻箱内温度场、 速度场的计算模拟 J .制冷学报， 1997 （1 ） ： 53-55.（Yu Bing-feng， Ge Jun， Wang Zhi-gang.Computerized Simulation ofTemperature and Velocity Field in Freezer J .Refrigeration Journal，1997 （1 ） ： 53-55. ）机械设计与制造No.8Aug.20184ChaoXing
