第32卷 第23期 2016钲 农业工程学报 Transactions of the Chinese Society of Agricultural Engineering Vb1.32 No.23 Dec.2016 213 12月 蒸发式过冷水制冰液滴蒸发结晶的模拟 马善军 ,李鹏辉2,孔令健2,李少华 ,韩吉田2※ (1.华北电力大学能源动力与机械工程学院,北京102206;2.山东大学能源与动力工程学院,济南250061) 摘要:蒸发式过冷水制冰系统中液滴传热传质特性研究具有重要意义。为了分析蒸发式过冷水制冰系统中液滴的热质 传输特性,该文将单个液滴冷却结晶过程分为全液相区、固.液两相区和全固相区3个不同区域,建立液滴在大空间内蒸 发降温结晶过程的数学模型,通过数值求解方法研究其在大空间蒸发室内降温结晶过程的温度变化特性,进一步分析液 滴和空气的入口参数与液滴整个蒸发结晶特性的关系,为蒸发式过冷水制冰系统的结构优化设计和运行管理提供参考依据。 关键词:蒸发;冷冻;冰;蒸发过冷水;冰蓄冷;冰浆;液滴;降温结晶 doi.10.119758.issn.1002-6819.2016.23.029 中图分类号:TK02 文献标志码:A 文章编号:1002—6819(2016)一23—0213—05 马善军,李鹏辉,孑L令健,李少华,韩吉田.蒸发式过冷水制冰液滴蒸发结晶的模拟[J].农业工程学报,2016,32(23): 213--217.doi:10.11975 ̄.issn.1002-6819.2016.23.029 http://www.tcsae.org Ma Shanjun,Li Penghui,Kong Lingjian,Li Shaohua,Han Jitian.Modeling of water droplet in super-cooling water evaporative system for ice slurry producfion[ ̄.Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSA ̄), 2016,32(23):213—217.(in Chinese Wim English abstract)doi:10.11975 ̄.issn.1002—6819.2016.23.029 http://www.tcsae.org 0引 言 冰蓄冷技术是实现电力“移峰填谷”和缓解尖峰电 力供需矛盾的有效途径之一,对于解决由于制冷空调导 致的用电矛盾具有重要意义。由于冰水混合物(冰浆) 具有流动性好和相变潜热大等优点,特别适合于大规模 和高能量密度的冷量储存和输送过程,基于冰浆技术的动 态冰蓄能技术己成为目前蓄冷领域的重要发展方向之--[ 1。 蒸发式过冷水制冰是一种冰浆制备新方法。其工作 原理为:水通过压力喷嘴后被雾化成小水滴,低温的湿 空气中的水滴由于蒸发吸热而被过冷到0℃以下,过冷 水滴在解冷装置的解冷作用下即可变成冰浆。该方法克 服了传统过冷水法制冰易出现冰堵等问题,且具有较高 的制冰效率[24】。因此,蒸发式过冷水冰浆制备技术被认 为是一种具有良好工程应用前景的制冰技术。 为了发展蒸发式过冷水冰浆制备技术,国内外已对 低温湿空气中液滴冷却结晶过程进行了研究。Kim等【5】 研究了蒸发式过冷水制冰系统中液滴全液相区的降温特 性,得到了其降温过程的模拟结果。Shin 6J对水雾蒸发 法制取球形冰颗粒过程进行了理论与试验研究,其发展 论分析和试验研究,并建立了其数学模型,所得模拟结 果与试验数据吻合较好。Satoh等[8]对真空室中水滴蒸发 过程的冷却/冻结现象进行了试验研究,发展了预测液滴 结晶过程的数学模型。B6d6carrats 9J研究了利用过冷水 流制取冰浆的方法,得到了过冷度、流量和制冷温度等 参数对结晶过程的影响。Heneghan等[1 0]对气泡加速纯过 冷水核化的液体.冰晶非均匀核化现象进行了研究。章学 来等[11]研究了不同种类成核添加剂对降低蓄冰成核过冷 度的影响,王世清等[1 】对基于热管的自然冷源制冰技术 方案进行了初步的研究。李福乐[1 3]研究了基于空气喷射 解除过冷的冰浆连续制备方法,其研究结果表明空气喷 射法可以有效解除水结冰过程中的过冷度。马军等[1 】和 李鹏辉【I5】对蒸发式过冷水法制冰系统中液滴的热质传输 特性进行了研究,通过将液滴的蒸发结晶过程分为全液相 区、固一液两相区和全固相3个不同区域,发展了描述液滴 冷却结晶过程的数学模型,对其进行了数值模拟研究。 本文在已有研究工作的基础上,以蒸发式过冷水冰 浆制备系统中垂直下降的液滴为研究对象,分别按照全 液相、固一液两相和全固相3个不同区域,建立大空间蒸 发室内单个液滴降温结晶过程的数学模型,研究不同液 滴和空气初始参数等对于整个液滴结晶过程的影响,得 到蒸发结晶过程中液滴参数的变化特性,为蒸发式过冷 水制冰系统的有效设计和可靠运行提供参考。 的扩散控制蒸发模型的预测结果与试验数据吻合良好。 Strub等[J7】对低温湿空气中水滴的冷却结晶特性进行了理 收稿日期:2016.03.31 修订日期:2016.08.27 芋金项目:国事自然科学基金 (2013013110006) ’ ¨。 ;教育部博士点基金 1 数学模型 。 作者简介:马善军,博士生,高级工程师,主要从事换热器设计、材料腐蚀 防护和传热传质方面研宄。北京华北电力大学能源动力与机械工程学院, 102206。Email:mashanjun@cgdc.COrn.cn 为了建立蒸发式过冷水冰浆制备过程中液滴蒸发结 ※通信作者:韩吉田,教授,博士,博士生导师,主要从事燃料电池、总能 系统和多相流方面研究。济南山东大学能源与动力工程学院,250061。 Email:jthan@sdu.edu.cn 晶的数学模型,假定液滴在低温低湿大空间中经历由降 温、结晶而达到结冰的过程,不考虑其解除过冷环节的 影响,并做如下简化假设:1)空气和水蒸气均是理想气 体;2)液滴形状为球形;3)液滴与空气为逆向流动,且 ● 214 农业工程学报(http://www.tcsae.org) 2016年 不计液滴辐射换热影响。过冷水滴蒸发结晶示意图见图1。 数;Pwa、Pwd、p 和pwd分别是空气和液滴表面的水蒸气 密度和压力, 和乃分别为空气和液滴的温度, 为水 ........ ▲ , ......●● 空气环境 Air environment 注: 和 分别为空气、液滴速度;Dd为液滴直径。 Note: and ua indicate speed of air and liquid droplet,respectively;Dd indicates diameter of liquid droplet. 图1过冷水滴蒸发结晶示意图 Fig.1 Schematic of evaporative crystallization of super-cooling water droplet 1.1 全液相区 全液相区内液滴与空气间的对流传热传质而导致液 滴的温度下降和直径减小。假设该区域液滴内部温度为 均匀分布,则其运动方程可表示为 警一( +G (1) FD: 1 ・ .cD.兀. = 1 .g-尢. (2) ’ b Ur=Ud—Ua,dH=Ud・dt (3) md ÷丌・ ・Pa (4) P: (5) 式中 、 和G分别为粘滞阻力、浮力和重力, 、 和 分别为空气、液滴速度和液滴与空气的相对速度,m 和Dd分别为液滴的质量和直径, 和 分别为液滴和空 气的密度, 和 分别为阻力系数和雷诺数,dH为计 算时间步长dt内液滴的下降高度,Va为空气的运动粘度。 当Re<1000,阻力系数 可由下式计算[8] cD= (1+0.125Re0.72) 传质方程可表示为 粤 .hm.(‰一Pwd) (6) Sh:2+0.6RP .& (7) So= , = ㈦ 根据理想气体状态方程PV=m・RM・T(P、 m、 I 和 分别是理想气体的压力、体积、质量、通用气体 常数和绝对温度),可得 P wa一 fl 一一 1 j (9) 式中h 为传质系数, 和&分别是舍伍德数和施密特 的二元扩散系数。D 可由下式确定 D’,:2-2×1 f,量 ×273) \f 1J (10) 式中 为空气压力。 传热方程可表示为 ( + = 1兀・ -Pw・ …) =2+0.6Re2.pr3 (12) h:蔓Nu (13) Dd 式中h、Nu和 分别为液滴表面的对流换热系数、努射 尔特数和普朗特数, 为水的密度, 为空气导热系数, c 和htg分别为水的定压比热和汽化潜热。 1.2固一液两相区 在固.液两相区,液滴的过冷状态由于过冷解除装置 的作用而被破坏。因此,液滴的温度可在瞬间升高到0℃ 而在液滴表面形成薄冰层[ 。 假设该区固相中具有温度梯度,其外表面温度为 , 液相水的温度为冰熔点温度 固相的初始质量为 :三 . 。.( 一 ) : 生 (14) 式中 为冰的密度,D 为水的直径,h 为水的固化潜热。 液相质量为 m =÷兀・Pw・ (15) 液滴总质量为 md= +mw (16) 液滴体积为 一 -1 一.6 Dj: Ps + Pw (1L/J 7) 固.液两相区内液滴的运动方程为 d :一( + )+G (18) 由于固.液两相区发生的是固相冰的升华过程,其传 质方程可表示为 .,.h msub.( 一Pwd) (19) md= + = 7c・[pw・D 3+ ・( 一 )] (20) 式中 为固相冰升华过程的传质系数,可由式(8)确定。 根据能量平衡原理,假定液滴结晶过程中的放热量 等于由固相导出的热量,且等于固相外表面的对流放热 量,则可得到下述的传热方程 警h/s--- ̄.札,㈡一 嚣≥ 第23期 马善军等:蒸发式过冷水制冰液滴蒸发结晶的模拟 215 勘札丝D ̄-Dw一 j(22 为温度对半径的微分, 5 4 式中, 为水的固化潜热, 旨0 1 ,a暑景鲁8 醚赠 和 分别为冰的升华潜热和导热系数。 1.3全固相区 当液相直径减小到0时固.液相变过程结束而进入全 固相区。假定在全固相区内发生冰的升华过程,且可用 固相区的平均温度表示其温度,可得到下述的运动和传 蓍 4 5 0 l 2 3 4 -Dropping height/m 下降高度 注:D 为液滴入口直径。 热传质方程 :一( + )+G (23) :兀. .h2 .( 一 ) (24) 一。 = c 曙 詈 式(25)中, 是冰的定压比热。 上面得到的式(1)到式(25)即为蒸发式过冷水制 冰系统中液滴降温结晶过程的完整数学模型。 2模拟结果与分析 髌艇 利用上述建立的数学模型,在时间步长df内通过计 算其参数值可得到液滴降温结晶过程中的参数变化关系。 2.1与试验结果比较 为了检验上述所建数学模型的可靠性,在液滴入口 直径为1.8 min,入口温度为6℃,空气温度为一12℃, 空气相对湿度为50%,空气流速为1.5 m/s的条件下,利 用该模型模拟了蒸发式过冷水制冰过程中液滴温度随时 间的变化关系,并与相同条件下文献[7]的研究结果进行 了比较(图2)。可以看出,该模型的模拟结果与文献[7] 的研究结果吻合较好,证明了该模型的可靠性。 O 5O 100 150 时间Time/s 图2本文模拟结果与文献【7]结果对比 Fig.2 Comparison of obtained simulation results by the developed model、Ⅳith data in reference[7] 2.2液滴结晶模拟结果与分析 蒸发式过冷水制冰系统中水经压力喷嘴雾化成为小 直径的液滴。下面通过液滴参数随下降高度变化的模拟 结果分析液滴结晶过程的参数变化关系。 图3是液滴入口温度和速度分别为2℃和10 m/s,空 气的温度、相对湿度和流速分别为0℃、50%和0的条件 下,不同入口直径液滴结晶过程中液滴温度随着下降高 度变化的模拟结果。由图3可以看出,入口液滴直径越 小,其温度下降速度越快,结晶过程需要的时间和下降 高度越短,降温结晶过程的速度也越快。因此,为了提 高液滴的结晶速度,在工程应用中应通过采用合适的喷 嘴参数来减小液滴的初始直径。 Note:Da indicates inlet diameter ofliquid droplet. 图3不同入口直径液滴温度随下降高度的变化 Fig.3 Temperature variation of liquid droplets with different inlet diameter along with dropping height 图4是液滴入口直径为200/xm,速度为10 m/s,空 气温度为0℃,相对湿度为50%,空气流速为0的条件 下,不同入口温度液滴结晶过程中液滴温度随着下降高 度变化的模拟结果。由图4可知,随着液滴入口温度的 升高其温度下降速度也提高,因为液滴入口温度的升高 导致其与空气间的传热增大。由于降低液滴的入口温度 需要消耗冷量,因此,液滴入口温度的选择需要综合考 虑液滴达到结晶所需过冷度的时间和系统总能耗等因素 的影响。 毯赠壤避 0 l Z 3 4 Droppi下降高度 ngheight/m 注: 为液滴入口温度。 Note:乃indicates inlet temperature ofliquid rdoplet. 图4不同入口温度液滴温度随下降高度的变化 Fig.4 Temperature variation of liquid rdoplets with different inlet temperature along with dropping height 图5是液滴入口直径为200pm,入口温度为2℃,速 度为10 m/s,空气相对湿度为50%,空气流速为0的条 件下,不同空气温度下液滴温度随着下降高度变化的模 拟结果。由图5可知,空气温度越低,沿着下降高度液 滴降温越快,这是由于空气温度越低,空气与液滴的温 差则越大,从而导致液滴温度降低得越快。 下降高度 Droppingheight/m 注: 为空气温度。 Note: indicates air temperature. 图5 不同空气温度下液滴温度随下降高度的变化 Fig.5 Temperature variation of liquid rdoplets wiht different air temperature along with rdopping height 216 农业工程学报(http://www.tcsae.org) 2016矩 0 1 2 3 4 下降高度 Dropping h ̄ight/m 注: 为空气速度。 Note: indicates air speed. 图7不同空气速度下液滴温度随下降高度的变化 Fig.7 Temperature variation of liquid droplets with diferent air speed along with rdopping height 3结论 本文在分析蒸发式过冷水制冰系统单个液滴蒸发结 晶过程中传热传质的基础上,建立了该过程的数学模型, 对液滴温度随系统参数的变化特性进行了模拟研究,得 到了以下结论: 1)减小喷入蒸发室的液滴直径可提高其蒸发结晶速 度,实际应用中应尽可能减少喷入蒸发室的液滴直径。 2)空气湿球温度决定了液滴的最小过冷温度。实际 应用中可通过选择较高的空气温度达到降低能耗的目的。 3)增大空气流速有利于系统制取冰浆产量的提高, 但空气速度过高会携带走较小粒径的液滴。因此,需要 根据液滴直径的大小选择合适的空气速度。 【参考文献】 [1]张小松,陈瑶,殷勇高,等.流态冰制取技术研究进展及 实验初探[J].东南大学学报(自然科学版),2013,43(6): 1343——1352. Zhang Xiaosong,Chen Yao,Yin Yonggao,et a1.Research progress and tentative experimental study on ice slurry producing technology[J].Journal of Southeast University (Natural Science Edition),20 1 3,43(6):1 343—1 352.(in Chinese with English abstract) [2】李秀伟,张小松.蒸发式过冷水制流态冰方法[J】.东南大 学学报(自然科学版),2009,39(2):269--275. Li Xiuwei,Zhang Xiaosong.Evaporative super-cooled water method for ice-slurry producing[J].Journal of Southeast Universiyt(Natural Science Edition),2009,39(2):269--275. (in Chinese wiht English abstract) [3]何国庚,吴锐,柳飞.冰浆生成技术研究进展及实验初探 [J].建筑热能通风空调,2006,25(4):22—27. He Guogeng,Wu Rui,Liu Fei.The review and experiment on dynamic ice slurry generated technology[J].Building Energy&Environment,2006,25(4):22—27.(in Chinese with English absrtact) [4】曲凯阳,江亿.高性能过冷水连续制冰系统[J].太阳能学 报,2002,23(3):317—321. Qu Kaiyang,Jiang Yi.The dynamic ice—making with higher performances[J].Acta Energiae Solaris Sinica,2002,23(3): 317—321.(in Chinese wiht English abstract) [5]Kim B S,Shin H T,Lee Y P,et a1.Study on ice slurry production by water spray[J].International Journal of Refrigeration,2001,24(2):176—184. 【6]Shin H T,Lee Y P,Jurng J.Spherical-shaped ice particle production by spraying water in a vacuum chamber[J]. Applied Thermal Engineering,2000,20(5):439--454. [7】Strub M,Jabbour O,Stnrb F,et a1.Expeirmental study and modeling of the crystallization of a water droplet[J]. International Journal ofRefirgeration,2003,26(1):59—68. [8]Satoh I,Fushinobu K,Hasihmoto Y.Freezing of a water rdoplet due to evaporation—heat trnasfer domin ̄ing the evaporation—freezing phenomena and hte effect of boiling on freezing characteristics[J]. International Journal of Refrigeration,2002,25(2):226--234. [9】B6d6carrats J,David T,Castaing-Lavignottes J.Ice slurry production using supercooling phenomenon[J].International Journal ofRefirgeration,2010,33(1):196--204. A F,Haymet A D J.Liquid-to-crystal heterogeneous nucleation:bubble accelerated nucleation ofpure supercooled water[J].Chemical Physics Letters,2003,368(s 1-2):177—182. 1]章学来,孟祥来,赵群志,等.冰蓄冷空调蓄冰成核[J].化 工学报,2014,65(1 1):4321--4324. Zhang Xuelai,Meng Xianglai,Zhao Qunzhi,et a1.Ice nucleation of ice thermal storage air—condiiton[J].CIESC Journal,2014,65(1 1):4321--4324.(in Chinese wiht English absrtact) 【12】王世清,宋庆武,张岩,等.基于热管的自然冷源制冰技 术方案的初步试验【J].农业工程学报,2008,24(1 1):23O一 232. Wang Shiqing,Song Qingwu,Zhang Yan,et a1.Preliminary study on technology for ice—making using natural cold resource based on heat pipe[J].Transactions of hte Society of Agricultural Engineering(Transactions of hte CSAE),2008, [10]Heneghan [1第23期 马善军等:蒸发式过冷水制冰液滴蒸发结晶的模拟 217 24(1 11:230--232.(in Chinese with English abstract) [13] 李福乐.基于空气喷射解除过冷的冰浆连续制备系统的研 ofawaterdropletinside avacuum chamber[J].Refrigerationand Air Conditioning,2012,26(3):213--218.(in Chinese with English absrtact) [15]李鹏辉.蒸发式过冷水法制冰浆的模拟研究[D】.济南:山 东大学,20t3: Li Penghui.Simulation study on Water Evaporative 究[D】.济南:山东大学,2013. Li Fule.Study on Ice Slurry Production System Based on Air Jet Lifting Supercooled State[D].Jinan:Shandong University, 2013.(in Chinese with English absrtact) 【14] 马军,赵红霞,韩吉田.水滴在真空室内结晶过程的模拟 『J].制冷与空调,2012,26(3):213--218. Ma Jun,Zhao Hongxia,Han Jitian.Modeling Crystallization Supercooling Water Method for Ice Slurry Production[D]. Jinan:Shandong University,2013.(in Chinese wih tEnglish absrtact) Modeling of water droplet in super・-cooling water evaporative system for ice slurry production Ma Shanjun ,Li Penghui ,Kong Lingjian ,Li Shaohua ,Han 3itian ※ f1.School ofEnergy,Power and Mechanical Engineering,Ⅳ0砌China Electric Power University,Beijing 102206,China; 2.SchoolofEnergyandPowerEngineering,ShandongUniversity,Jinan 250061,China) Abstraet:Ice storage technology is considered as one of the most promising options to achieve the so.called‘'peak load leveling of power system”and to relieve the contradiction between the supply and demand of peak power currently encountered in China.Among the ice.producing methods being developed around the world.the evaporative supercooling water ice.making one is a new and efficient way for ice slurry production of large scale without ice blockage.Therefore.it is of signiicafnt importance to investigate the heat and mass仃ansfer characteristies during the cooling and crystallization process ofwater droplets in the evaporative supercooling water ice—making chamber to the development ofpracticalice.making system. In order to analyze the heat and mass仃ansfer characteristies of water ropldets in an evaporative super-cooling water system for ice slurry production.we proposed a mathematical model for the cooling and crystallization process of a single water droplet falling in the evaporation chamber with 1arge space.which comprehensively took into account the three different zones of entire liquid phase.solid-liquid interphase and entire solid phase during the cooling and crystallization process of a water roplet.The developed matdhematical model was then validated by use of he ttheoretical and experimental results presented in he availtable literature and satisfactory agreement was achieved in between the mode1 simulation results and the research ones reported in the literature.indicating the correc仃1ess of the mathematical mode1.The parameter variations of he water dtroplet with changes in some of the key system operating variables,such as the inlet size and temperature of water droplet,the flow rate and relative humidity of the cold air,were numerically determined by solving the developed mathematical mode1.Effects of he itnlet size and temperature of water droplet,inlet temperature,flow rate and relative humidiy of tthe cold air in the evaporation chamber on the entire cooling and crystallization process were analyrzed in detail.The results indicated mat the cold air temperature,relative humidiy tnd alfow rate had signiifcant effect on he coolting and crystallization process of water roplets.And dchange in the air temperature had the most obvious impact on the cooling rate of water droplets.Both of he aitr temperature and relative humidiy dettermined he mitnimum temperature fi.e.。the wet.bulb temperature)that he dtroplets can reach in the evaporative supercooling water chamber for ice slurry production.Air velociy wast one of the key factors influencing the ice crystal production.It can be GOnGluded from the obtained simulation results that decrease in the inlet diameter of water droplets entering the evaporative supercooling chamber was he1D 1 to enhance the speed of evaporative and crystallization process of water droplets.As a result.it was recommended to incorporate water ropldets with diameters as small as possible in practical applications.In order to reduce the energy consumption for ice—making,it was important to employ a higher air temperature in meeting the wet—bulb temperature requirement.Whi1e increasing the cold air speed can improve the ice slurry output of ice—making system,it may result in increase in the carrying volme of uwater droplets.Therefore。the cold air speed should be determined according to the corresponding diameters of water droplets entering the supercooling evaporative chamber.The conelusions obtained n ithis PaDer were of significance to the e伍cient structural design and practical operation of the super.cooling water evaporative system for ice slurry production of large scale. Keywords:evaporation;freezing;ice;evaporative super-cooling water;ice cold storage;ice slurry;water droplet;cooling and crystallization