冰浆流动过程中冰晶粒径演化规律的实验研究.pdf

第 44 卷 第 5 期2023 年 10 月冰浆流动过程中冰晶粒径演化规律的实验研究Vol.44,No.5October,2023文章编号:0253-4339(2023)05-0131-08doi:10.3969/j.issn.0253-4339.2023.05.131冰浆流动过程中冰晶粒径演化规律的实验研究蔡玲玲1 米 沙1 刘志强2(1 中国石油大学(北京)克拉玛依校区 重质油国家重点实验室 克拉玛依 834000;2 中南大学能源科学与工程学院 长沙 410083)摘 要 本文实验研究了流速为 1.02.5 m/s、质量含冰率为 10%30%、初始粒径为 350550 m 条件下的冰晶粒径演化规律,分析了流速、冰晶颗粒浓度及初始粒径对冰晶粒径演化规律的影响。结果表明:过冷法制备的冰晶粒径分布符合对数正态分布;在湍流条件下,出口平均粒径明显小于进口,表明破碎比团聚对冰晶粒径演化的影响显著,使大粒径冰晶颗粒减少,小粒径冰晶颗粒增多;大粒径的冰浆颗粒更容易破碎成小粒径颗粒,并且随着流速的增大冰晶粒径破碎作用更显著。当质量含冰率为10%,速度由 1.0 m/s 增至 2.5 m/s 时,进出口平均粒径差(d32)由 16.1 m 增至 34.6 m,这对冰浆的安全输运具有重要的指导价值。关键词 冰浆;粒径演化;团聚;破碎中图分类号:TB61+1;TB64文献标识码:AExperimental Study of Evolution of Ice Crystal Size in Turbulent Flow of Ice SlurryCai Lingling1 Mi Sha1 Liu Zhiqiang2(1.State Key Laboratory of Heavy Oil Processing,China University of Petroleum-Beijing at Karamay,Karamay,834000,China;2.School of Energy Science and Engineering,Central South University,Changsha,410083,China)Abstract Experiments were conducted with flow rates between 1.0 m/s and 2.5 m/s,ice mass fractions between 10%and 30%,and initial particle sizes between 350 m and 550 m to study the effects of parameters on the evolution of ice crystal size in ice slurry isothermal flow.The results indicated that the ice particle-size distribution of the ice slurry prepared using a supercooling method presented a lognormal distribution.Under turbulent conditions,the Sauter mean diameter at the outlet was significantly smaller than that at the inlet.This indicates that the effect of breakage on the ice particle-size distribution was more significant than that of aggregation,resulting in an increase in small ice particles and a decrease in large ice particles.Large ice particles are more likely to break into smaller particles;the effect of breakage becomes more significant with an increase in the flow rate.With an ice mass fraction of 10%,the velocity increased from 1.0 m/s to 2.5 m/s,and the Sauter mean diameter difference between the inlet and outlet(d32)increased from 16.1 m to 34.6 m.These findings have significance in the safe transportation of ice slurries.Keywords ice slurry;particle-size evolution;aggregation;breakage 收稿日期:2022-07-05;修回日期:2022-10-21 冰浆是一种由载流溶液和冰晶颗粒组成的固-液两相流体1,具有良好的流动性、蓄冷密度高、无污染等优点,因此被广泛应用于建筑和矿井的蓄冷空调系统2、运动后冷却3、食品预冷4-5等领域。冰浆在实际应用过程中,涉及冰浆的管内输送环节,其呈现的流体动力学特征是确保冰浆安全输运的重要理论基础。由于冰晶颗粒与载流体间存在密度差,不同的流动工况冰浆呈现出不同的流动型态,主要包括均质流、非均质流、移动床流和固定床流6。A.Bordet等7采用高速摄像机拍摄并分析了冰晶粒径在 150350 m 范围的冰浆在水平圆管内的流型特征,研究发现改变冰浆流速、含冰率等条件,冰浆会呈现均匀流、分层流和环状流动。D.Vuarnoz 等8实验发现冰浆在均匀流时,速度呈对称分布,而在非均质性流动时,速度呈非对称分布,因而可通过速度分布判断冰浆流型特点。P.Doron 等9指出对于冰浆浓度分布特性的研究应考虑冰晶粒径分布(particle size dis-tribution,PSD)及冰浆流动型态的影响。同时,在冰浆流动过程中,冰晶颗粒将发生团聚、破碎、磨损等动力学行为,使冰晶粒径分布发生变化,从而对冰浆流动力学特性产生影响,甚至导致“冰堵”。131第 44 卷 第 5 期2023 年 10 月制 冷 学 报Journal of RefrigerationVol.44,No.5October,2023目前人们对冰晶颗粒间的相互作用规律和冰晶粒径演化机理的研究多集中于冰晶的形成10-12及冰浆的存储阶段13-15,对冰浆流动与换热过程冰晶颗粒动力学原理和冰晶粒径演化及分布规律的研究较少。Xu Dan 等16-17以数群平衡(population balance model,PBM)为基础,基于经典颗粒团聚、破碎动力学理论,开展了关于管内冰浆流动特性以及冰晶粒径演化特性的数值模拟研究。在该模型的基础上,邹良旭等18对冰浆在水平管内的流动换热特性及冰晶粒径演化规律进行了数值模拟研究。虽然他们采用该模型成功实现了冰晶生长(融化)、团聚、破碎等动力学过程的数值仿真,但由于目前还未有关于冰浆管内流动过程冰晶粒径演化规律的实验研究,该动力学模型中的关键参数采用了化工领域颗粒动力学特性进行计算,导致该模型在计算冰晶颗粒动力学行为时难以用实验数据验证其准确性。本文主要对水平管内冰浆流动过程冰晶粒径演化特性展开实验研究,研究不同流速、质量含冰率和初始粒径条件下冰晶颗粒的团聚、破碎等动力学行为对冰晶粒径分布特性的影响规律。研究结果为冰浆流动特性研究提供了理论基础,可为冰浆安全输运提供建议和指导。1 实验1.1 实验系统 实验系统包括过冷法冰浆制备及存储系统、流动测试管路、冰晶图像采集系统以及数据采集系统等,如图 1 所示。图 1 实验系统Fig.1 Experimental system 过冷法冰浆制备及存储系统包括制冷循环、载冷剂循环、过冷水循环以及控制系统,用于冰浆的制备及存储。制冷循环主要由压缩机、蒸发器、冷凝器和膨胀阀等设备组成,以 R22 为制冷剂进行循环制冷,为系统提供冷量。载冷剂循环主要包括两个板式换热器、乙二醇泵、缓冲罐以及膨胀罐等,以乙二醇为载冷剂将制冷循环产生的冷量持续、稳定地传递给过冷循环。过冷水循环主要包括过冷换热器、超声波促晶器、储冰槽、过滤器、循环泵等设备。储冰槽内的水溶液由循环泵输送至过冷换热器内,与载冷剂循环中的乙二醇进行热量交换后处于过冷状态,过冷水溶液在超声波促晶器内受到扰动后解除过冷状态发生相变形成晶核,生成的冰浆经管道输送回储冰槽内。为了防止储冰槽内的冰晶颗粒随水溶液进入过冷换热器内诱发结冰,在储冰槽内设置过滤器实现冰水分离。此外,为了防止储冰槽内的冰浆分层,在储冰槽内设置直流电机驱动搅拌器实现冰浆的均匀存储。测试管段为长度为 1.2 m,内径为 14 mm,壁厚为 1 mm(测量精度0.1 mm)的铜管。在测试段的前、后端均设置了 0.6 m 的入口段和出口段,以减小测试段前后两端的弯头对测试段流体流动的影响,确保冰浆在进入测试段时达到充分发展。测试管段和循环管段外均包裹了多层绝热材料,以减小周围环境与管内冰浆的换热。测试管段管材及绝热材料的物231第 44 卷 第 5 期2023 年 10 月冰浆流动过程中冰晶粒径演化规律的实验研究Vol.44,No.5October,2023性参数如表 1 所示。实验过程中冰浆流过测试管后的最大温升为 0.2 K,此时冰浆所处的流动工况可近似视为等温流动。表 1 管材及绝热材料物性参数Tab.1 Physical parameters of pipe and insulation materials材料密度/(kg/m3)导热系数/W/(m K)比热容/J/(kg K)铜8 900401385绝热材料950.0341.2 实验参数及测量方法 1)压降压降的测量采用麦克传感器的 MDM484A 型差压变送控制器。在测试管段的进、出口两端分别设置了取压口,通过螺纹将导压管与取压口进行连接。为避免气体及冰晶颗粒进入导压管从而影响压差计的测量结果,取压口设置在测试管段的侧面。2)温度温度测量采用上海自动化仪表三厂定制的插入式 Pt100 温度传感器。在测试管段的进、出口位置各布置一个插入式温度传感器,安装时在测温点预留三通,通过螺纹将插入式温度传感器固定在三通上,并确保热电偶的感温点在管道中心。3)流速采用麦克传感器的 MFE600 型电磁流量计测量流体的体积流量,通过体积流量、流速及管道截面积的关系得到流体的截面平均流速。为了避免流体在管道中冲刷产生的电势对测量结果产生误差,需对电磁流量计外壳单独设置接地。4)质量含冰率采用称重法测量冰浆的质量含冰率。在测试管段的进、出口各设置一个取样口,测量时从取样口取样并称重,随后采用不锈钢细筛将样品中的冰颗粒进行分离并称重。被测冰浆的质量含冰率可以通过测得的样品分离前后的质量计算得到。冰晶分离时采用 150 目(筛孔尺寸约为 100 m)不锈钢细筛网。为避免取样及分离过程冰晶融化引起较大的测量误差,实施多次少量的取样方法,并且取样时用到的烧杯需要提前预冷至 0 以下。5)粒径冰晶粒径的测量采用显微拍摄的方法。通过测试管段进、出口处的取样口进行取样,迅速转移至显微拍摄仪进行图像采集。冰晶图像采集前调整好显微镜的放大倍数并进行标定。用显微拍摄方法获得冰晶粒径的显微图像,再利用 ImageJ 软件对显微图像进行处理获得冰晶粒径。图 2 所示为冰晶图像采集系统的实物图。图像采集系统主要由显微拍摄仪、计算机、图像处理及分析软件等组成。图 2 图像采集系统Fig.2 Image acquisition system为了避免在图像采集时冰晶发生融化对粒径测量产生影响,取样前对载玻片、镊子等器具进行预冷处理,同时将显微拍摄仪置于恒温箱内使其温度低于-1。显微拍摄得到的冰晶图像需通过 ImageJ 软件进行处理以获取冰晶颗粒粒径,由冰晶的二维投影面积 Ai计算当量圆直径作为冰晶的特征尺寸 di。其中,下标 i 表示第 i 个冰晶颗粒。di=4Ai/(1)冰晶的平均粒径以索特平均直径 d32进行描述,可通过每个冰晶的粒径及冰晶颗粒的总数 N 计算得到。d32=Ni=1d3i/Ni=1d2i(2)为减小抽样误差,在实验过程中每组工况的样本取样均保证样本容量300 个冰晶,进行有意义的统计分析。实验所需的主要测量仪器如表 2 所示。表 2 主要测量仪器Tab.2 Main measuring instruments主要仪器型号测量范围精度电磁流量计MFE60008 m3/h0.5%FS压差计MDM484A060 kPa0.5%FS温度传感器Pt100-50200 0.2%FS电子秤05 kg0.01 g采集卡DAM-30551.3 实验过程及工况 过冷法制备的冰浆含冰率较低,质量含冰率 约为 2.5%19。为了获得实验所需的,冰浆被循环连续制取,本实验台的过冷法冰浆制备机组制备出 331第 44 卷 第 5 期2023 年 10 月制 冷 学 报Journal of RefrigerationVol.44,No.5October,202330%的冰浆至少需要2 h。因此,实验中对制冰时间分别为 2、2.5、3、3.5 h 储冰槽内的冰晶颗粒粒径进行了测量。由实验数据可知,冰晶颗粒粒径分布约在101 000 m 范围内,采用离散概率分布表示冰晶颗粒粒径分布规律,对测得的粒径进行均匀分区统计,每个区间的宽度为 20 m。实验测得不同制冰时刻无添加剂冰浆冰晶颗粒粒径分布如图 3 所示,并根据实验数据拟合得到了正态分布、对数正态分布、威布尔分布以及伽马分布曲线。由图 3 可知,对数正态分布概率密度函数能够较好地预测冰晶粒径分布的峰值以及最小和最大粒径概率密度,故本文采用对数正态分布描述冰晶粒径演化规律。随着制冰时间的增加,冰浆冰晶粒径分布曲线明显变矮变胖,峰值逐渐减小,即小粒径冰晶颗粒数量不断减少,大粒径冰晶颗粒数量不断增加。一方面,因为制备的冰浆溶液含有一定的过冷度还未完全释放,随着制冰时间的增大,冰晶颗粒在该过冷度的驱动下不断长大。另一方面,由于存储过程中冰晶颗粒的团聚作用明显比破碎作用强烈15,随着时间推移冰晶颗粒尺度不断增大。实验时为了防止储冰槽内的冰浆分层,在储冰槽内设置直流电机驱动搅拌器实现冰浆的均匀存储。测量时,调好冰浆的流量和质量含冰率,待系统稳定后在测试管进、出口处取样进行粒径测量,即得到流动前后颗粒的粒径分布。实验的有效变量分别为流图 3 不同时刻无添加剂冰浆冰晶颗粒粒径分布Fig.3 Ice crystal size distribution of pure water ice slurry at different time速 u、质量含冰率 和初始冰晶粒径 d32,i。通过改变这 3 个变量,对不同操作工况下冰晶粒径演化规律进行测量,实验工况如表 3 所示。表 3 实验工况Tab.3 Working conditions of the experiments工况u/(m/s)/%d32,i/mI1.510350II1.510450III1.510550IV1.520450V1.530450VI2.010450VII2.510450 由于实验过程中较难精准控制进口冰浆质量含冰率和粒径分布,因此在实验数据分析时,将 和d32,i分为 8%12%、18%22%、28%32%和 335 365 m、435465 m、535565 m 三组,并分别记为 10%、20%、30%和 350、450、550 m。431第 44 卷 第 5 期2023 年 10 月冰浆流动过程中冰晶粒径演化规律的实验研究Vol.44,No.5October,20232 实验结果与讨论2.1 流速影响 不同流速下进、出口冰晶颗粒粒径分布如图 4 所示。由图 4(a)可知,与进口冰晶颗粒粒径分布相比,出口冰晶颗粒粒径分布曲线向左移动、峰值向左上方移动。即大粒径冰晶颗粒的数量降低,小粒径冰晶颗粒的数量增多。随 u 的增大,出口冰晶颗粒粒径分布曲线左移幅度及峰值均略有增大。由图 4(b)可知,进口处 200 m 以下颗粒的数量密度约为 25%;u 为1.5 m/s 时,200 m 以下颗粒的数量密度升至 38%以上;u 升至 2.5 m/s 时,出口处 200 m 以下颗粒的数量密度升至约 45%。图 4 不同流速下进、出口冰晶颗粒粒径分布Fig.4 Ice crystal size distribution at inlet and outlet underdifferent flow rates冰浆在管内流动时,一方面,由于湍流作用冰晶颗粒相互碰撞并融合在一起形成较大粒径的冰晶颗粒,导致冰晶颗粒数量减少,大粒径颗粒增多;另一方面,较大粒径的冰晶颗粒受流体剪切、压力梯度等作用发生破碎,导致冰晶颗粒的数量增多,且主要集中于小粒径范围。由此可知,为 10%、d32,i为 450 m的冰浆在管内流动时破碎作用大于团聚作用。在相同 和 d32,i下,u 的增大使管内冰晶颗粒的破碎作用增强,最小破碎粒径逐渐减小。这是因为 u 增大增强了流体剪切作用,同时湍流强度增强,颗粒碰撞概率增大,导致冰晶颗粒破碎作用增强。根据实验结果,计算获得不同流速下出口截面冰晶的平均粒径 d32,o及进、出口平均粒径的变化量d32,如图 5 所示。随着 u 的增大,d32,o由 431 m 降至 413 m,d32由16.1 m 增至34.6 m,d32增大了两倍多。由此可知,为 10%、d32,i为 450 m 的冰浆在管内流动过程中破碎作用大于团聚作用,并且在相同的 和 d32,i下,随着 u 的增大管内冰晶颗粒的破碎作用增强。图 5 流速对出口截面冰晶平均粒径 d32,o的影响Fig.5 Effect of flow rates on d32,o2.2 含冰率影响 不同质量含冰率下进、出口冰晶颗粒粒径分布如图 6 所示。由图 6(a)可知,出口冰晶颗粒粒径分布曲线与进口冰晶颗粒粒径分布相比向左移动、峰值增大。随着 的增大,冰晶颗粒粒径分布曲线左移幅度增大,峰值也增大。但当 由 10%增至 20%时变化微小,而当 增至 30%时,出口冰晶颗粒粒径分布曲线明显左移,且峰值也有较为显著的增大。进一步由图 6(b)中累计数量密度变化可知,由 10%升至20%时,出口处 200 m 以下颗粒的数量密度由 38%升至 41%;增至 30%时,200 m 以下颗粒的数量密度升至 51%,小粒径颗粒数量增多,大粒径颗粒数量减小。质量含冰率对出口截面冰晶平均粒径的影响如图 7 所示。由进、出口冰晶平均粒径分析可知,u 为1.5 m/s、d32,i为 450 m 的冰浆在管内流动过程中,随着 的增大,d32,o由 426 m 降至 418 m,d32由21 m 增至 30 m,d32增大了近 1.5 倍。实验结果说明,对于 u 为 1.5 m/s、d32,i为 450 m 的冰浆其在管内流动时破碎作用依然大于团聚作用,且随着 的增大,颗粒之间的碰撞概率增大,破碎作用增大效果更显著。531第 44 卷 第 5 期2023 年 10 月制 冷 学 报Journal of RefrigerationVol.44,No.5October,2023图 6 不同质量含冰率下进、出口冰晶颗粒粒径分布Fig.6 Ice crystal size distribution at inlet and outlet under different ice mass fractions图 7 质量含冰率对出口截面冰晶平均粒径的影响Fig.7 Effect of ice mass fractions on d32,o2.3 初始冰晶粒径影响 不同初始粒径进出口冰晶颗粒粒径分布如图 8所示。与进口冰晶颗粒粒径分布相比,所有初始粒径条件下的出口冰晶颗粒粒径分布曲线向左移动、峰值向左上方移动。即大粒径冰晶颗粒数量减少,小粒径冰晶颗粒数量增多。随着 d32,i的增大,d32,o分布曲线左移幅度增大。进一步由图 8(b)可知,d32,i为 350 m 时,进口处 200 m 以下颗粒的数量密度约为45%,出口处 200 m 以下颗粒的数量密度约升至53%,小粒径颗粒占比提高了 17.8%;d32,i为 450 m时,进口处 200 m 以下颗粒的数量密度约为 25%,出口处 200 m 以下颗粒的数量密度约升至 38%,小粒径颗粒占比提高了 52%;d32,i为 550 m 时,进口处 200 m 以下颗粒的数量密度约为 12%,出口处200 m 以下颗粒的数量密度约升至 24%,小粒径颗粒占比提高了 100%。即随着 d32,i的增大,冰浆在管内流过后出口处小粒径颗粒数量占比的增量越多,d32,i较大的冰晶颗粒破碎作用越大。图 8 不同初始粒径进、出口冰晶颗粒粒径分布Fig.8 Ice crystal size distribution at inlet and outlet under different initial particle size由进、出口冰晶平均粒径可知,为 10%、u 为1.5 m/s 的冰浆在管内流动过程中,随着 d32,i的增大,d32,o由 346 m 增至 507 m,d32由 10 m 增至41 m,d32增大了 4 倍多。结果说明,对于 为10%、u 为 1.5 m/s 的冰浆,其在管内流动时破碎作用大于团聚作用,并且随着 d32,i的增大,破碎作用显著增强,即大粒径的冰晶颗粒更容易发生破碎形成小粒径的冰晶颗粒。可以发现,湍流条件下冰晶颗粒的破碎作用占主导,但不同变量对于冰晶粒径分布的影响程度有所差异。改变流速对冰晶粒径分布影响最大,流速越大,破碎作用越显著;冰浆的质量含冰率越高,颗粒之间的碰撞概率增大,破碎作用增大效果更显著;初始冰晶粒径越大,越容易发生破碎。631第 44 卷 第 5 期2023 年 10 月冰浆流动过程中冰晶粒径演化规律的实验研究Vol.44,No.5October,2023图 9 初始粒径对出口冰晶平均粒径的影响Fig.9 Effect of initial particle size on d32,o3 结论 本文以无添加剂冰浆为实验介质,搭建了冰晶粒径测试以及冰浆流动实验平台,研究了初始冰晶粒径、流速以及质量含冰率对冰晶粒径演化规律的影响,得到如下结论:1)实验制备获得的无添加剂冰浆的冰晶颗粒形状多为椭圆形或近似圆形,对数正态分布可较好地预测无添加剂冰浆冰晶颗粒粒径分布。2)湍流运动条件下冰晶颗粒粒径分布曲线向小粒径方向移动,表明颗粒破碎作用大于其团聚作用,即湍流时破碎作用占主导。3)大粒径的冰晶颗粒更容易发生破碎形成小粒径的冰晶颗粒,且流速的增大可显著减小冰晶颗粒的平均粒径。当质量含冰率为 10%,速度由 1 m/s 增至2.5 m/s 时,d32由 16.1 m 增至 34.6 m,d32增大了两倍多。本文受中国石油大学(北京)克拉玛依校区科研启动基金(XQZX20220013),新 疆 维 吾 尔 自 治 区 自 然 科 学 基 金(2022D01F39)和湖南省教育厅科学基金(20C0651)项目资助。(The project was supported by Research Foundation of China University of Petroleum-Beijing at Karamay(No.XQZX20220013),Natural Science Foundation of Xinjiang Uygur Autonomous Region(No.2022D01F39)and Scientific Research Foundation of Hunan Provincial Education Department(No.20C0651).)符号说明A 投影面积,m2d 当量直径,md32 索特平均直径,md32,i 进口截面(初始)平均直径,md32,o 出口截面平均直径,md32 进、出口平均直径变化量,mN 颗粒数量,个u 流速,m/s 质量含冰率,%参考文献1 EGOLF P W,KAUFFELD M.From physical properties of ice slurries to industrial ice slurry applications J.International Journal of Refrigeration,2005,28(1):4-12.2 KAUFFELD M,GUND S.Ice slurryhistory,current technologies and future developments J.International Journal of Refrigeration,2019,99:264-271.3 NAKAMURA M,NAKAMURA D,YASUMATSU M,et al.Effect of ice slurry ingestion on core temperature and blood pressure response after exercise in a hot environmentJ.Journal of Thermal Biology,2021,98:102922.4 LI Xiang,WU Weidong,LI Kun,et al.Experimental study on a wet precooling system for fruit and vegetables with ice slurryJ.International Journal of Refrigeration,2022,133:9-18.5 ZHANG Roujia,CHENG Zhiming,DING Fuyuan,et al.Improvements in chitosan-based slurry ice production and its application in precooling and storage of Pampus argenteusJ.Food Chemistry,2022,393:133266.6 WANG J,BATTAGLIA F,WANG S,et al.Flow and heat transfer characteristics of ice slurry in typical components of cooling systems:a review J.International Journal of Heat and Mass Transfer,2019,141:922-939.7 BORDET A,PONCET S,POIRIER M,et al.Flow visualizations and pressure drop measurements of isothermal ice slurry pipe flowsJ.Experimental Thermal and Fluid Science,2018,99:595-604.8 VUARNOZ D,SARI O,EGOLF P W,et al.Ultrasonic velocity profiler UVP-XW for ice-slurry flow characterisa-tionC/Proceedings of the third international symposium on ultrasonic Doppler method for fluid mechanics and fluid engineering.Lausanne,Switzerland,2002.9 DORON P,BARNEA D.Flow pattern maps for solid-liquid flow in pipesJ.International Journal of Multiphase Flow,1996,22(2):273-283.10 SHIMADA W,FURUKAWA S.Rapid growth of ice crystal dendrite tips in dilute solution of trehaloseJ.Journal of Crystal Growth,2018,493:25-29.11 COLUCCI D,FISSORE D,BARRESI A A,et al.A new mathematical model for monitoring the temporal evolution of the ice crystal size distribution during freezing in pharmaceutical solutions J.European Journal of Pharmaceutics and Biopharmaceutics,2020,148:148-731第 44 卷 第 5 期2023 年 10 月制 冷 学 报Journal of RefrigerationVol.44,No.5October,2023159.12 余云霞,刘志强,唐艺芳,等.水平直管内冰浆流动过程冰晶粒径演化研究J.制冷学报,2020,41(2):41-47.(YU Yunxia,LIU Zhiqiang,TANG Yifang,et al.Investigation on the evolution of ice crystal size during ice slurry flowing in horizontal straight pipesJ.Journal of Refrigeration,2020,41(2):41-47.)13 刘曦,林淑娴,李岁,等.二氧化硅纳米颗粒对冰浆中冰晶粒径分布及存储演化特性的影响J.化工学报,2017,68(3):870-878.(LIU Xi,LIN Shuxian,LI Sui,et al.Effect of nanosilica on size distribution and evolution of ice crystal particles during storage of ice slurry J.CIESC Journal,2017,68(3):870-878.)14 赵腾磊,刘志强,徐爱祥,等.冰浆存储过程冰晶粒径演化数值模拟J.中南大学学报(自然科学版),2014,45(10):3651-3656.(ZHAO Tenglei,LIU Zhiqiang,XU Aixiang,et al.Numerical simulation of evolution of ice crystal size distribution during storageJ.Journal of Central South University(Science and Technology),2014,45(10):3651-3656.)15 徐爱祥,刘志强,赵腾磊,等.冰浆存储过程中冰晶粒径动力学演化影响因素J.中南大学学报(自然科学版),2015,46(8):3138-3144.(XU Aixiang,LIU Zhiqiang,ZHAO Tenglei,et al.Factors influencing dynamics evolution of ice crystals during ice slurry storageJ.Journal of Central South University(Science and Technology),2015,46(8):3138-3144.)16 XU Dan,LIU Zhiqiang,CAI Lingling,et al.A CFD-PBM approach for modeling ice slurry flow in horizontal pipesJ.Chemical Engineering Science,2018,176:546-559.17 CAI Lingling,LIU Zhiqiang,MI Sha,et al.Investigation on flow characteristics of ice slurry in horizontal 90 elbow pipe by a CFD-PBM coupled modelJ.Advanced Powder Technology,2019,30(10):2299-2310.18 邹良旭,马非,孟昭男,等.基于群体平衡模型的冰浆流动与传热特性数值研究J.上海交通大学学报,2019,53(12):1459-1465.(ZOU Liangxu,MA Fei,MENG Zhaonan,et al.Numerical investigation of flow and heat transfer of ice slurry based on population balance model J.Journal of Shanghai Jiao Tong University,2019,53(12):1459-1465.)19 KAUFFELD M,KAWAII M,EGOLF P W.Handbook on ice slurries-fundamental and engineering M.Paris:International Institute of Refrigeration,2005.通信作者简介米沙,男,讲师,中国石油大学(北京)克拉玛依校区重质油国家重点实验室,15111154203,E-mail:misha 。研究方向:浆体储能关键技术研究与应用。About the corresponding authorMi Sha,male,lecturer,State Key Laboratory of Heavy Oil Processing,China University of Petroleum-Beijing at Karamay,86-15111154203,E-mail:misha .Research fields:research and application of the key technology of ice slurry energy storage.831