流速对盘管内冰浆流动及换热特性的影响研究.docx

1、doi：10.3969/j.issn.2095-4468.2019.0i.207流速对盘管内冰浆流动及换热特性的影响研究姬长发二王展荣，姬晨阳，李美晨，陈柳，鱼欣媛(西安科技大学能源学院，陕西西安710054)摘要冰浆作为一种性能优越的潜热蓄能和输送介质，在蓄冷空调中具有极大的应用潜力。但由于冰浆输运、储存过程中受到诸多因素的影响呈现复杂的流动特征，如分层、聚集等，容易堵塞管道，目前仍未大规模应用。文章模拟研究了一定冰浆含冰率以及管外径条件下入口速度对盘管内冰浆流动及换热特性的影响。结果表明：冰浆在盘管的流动过程中，弯管部位都会出现与单相液体流动样的二次流，速度越大、二次流现象越明显。入口速度

2、对冰浆在盘管内流动与换热特性影响较大，流速过小，管道顶部冰粒子聚集易形成滞留层，不利于冰浆流动，但有利于冰粒子的相变释放潜热；流速过大时，管道内冰粒子均匀悬浮，流动状态好，但冰浆在管道内持续的时间短，不利于相变释放潜热，针对本文条件，冰浆在末端设备流动入口速度不应超过0.5m/s,比较理想的入口速度为0.25m/s左右。关键词蓄冷空调：冰浆；盘管；流动；换热StudyonInfluenceofFlowVelocityonFlowandHeatTransferCharacteristicsofIceSlurryinCoilTubesJIChangfaWANGZhanrong,JIChenyang

3、,LIMeichcn,CHENLiu,YUXinyuan(SchoolofEnergyResources,XianUniversityofScienceandTechnology,Xian,Shaanxi710054,China)AbstractAsakindoflatentheatstorageandtransportmediumwithsuperiorperformance,iceslurryhasgreatapplicationpotentialinstorageairconditioning.However,duetotheinfluenceofmanyfactorsintheproc

4、essoficeslurrytransportationandstorage,i(presentscomplexflowcharacteristics,suchasstratificationandaggregation,whicharceasytoblockthepipeline.Sofar,ithasnotbeenappliedinalargescale.Inthispaper,theeffectsoftheinletvelocityontheflowandheattransferoficeslurryinthetubeundertheconditionoficecontentanddia

5、meterof(heiceslurryarestudied.Theresultsshowthatintheflowprocessoftheiceslurryinthecoil,thesecondaryflowofthepipeappearsassameasthesinglephaseliquidflow;thegreaterthevelocityis,themoreobviousthesecondaryflowphenomenonis.Theinletvelocityhasgreatinfluenceontheflowandheattransfercharacteristicsoficeslu

6、rryinthecoil.Whentheflowvelocityistoosmall,theaccumulationoficeparticlesatthetopofthepipeiseasytoformastagnantlayer,whichisnotconducivetotheflowoficepulp,butitisbeneficialtorelease(helatentheatoftheiceparticles.Whentheflowvelocityistoolarge,theiceparticlesinthepipearesuspendedevenlyandtheflowstateis

7、good,butthedurationoftheiceslurryinthepipelineisshort,whichisnotconducivetoreleasethelatentheatofthephasechange.Accordingtotheconditionsofthispaper,theinletspeedoficeslurryinterminalequipmentshouldnotexceed().5m/s,andtheidealinletspeedisabout0.25m/s.KeywordsColdstorageairconditioning;Iceslurry;Coil;

8、Flow;Heattransfero弓I言网昼夜负荷差值巨大。为了缓解白天高峰时期的电网压力，蓄冷空调技术对电力负荷的“移峰填谷集中空调应用过程中的负荷不平衡性使得电作用显现叫冰浆具有较高的热输送能力，在区域本姬长发(1963-),男，教授，博士。研究方向：制冷空调节能技术。联系地址：西安市雁塔路58号，邮编：710054c联系电话E-mail：jicf。基金项目：国家自然青年科学基金(No.51404191)。木论文选自2018年第十届全国制冷空调新技术研讨会。的冷却系统中，如果使用含冰率为5%的冰浆，在相同流速下，其冷却能力超过传统系统的两倍以上，且运输、储存无

9、需改变泵送设备、配电网络或蓄积罐12气冰浆制冷系统在一些发达国家，如美国、口本和丹麦等国家有较广泛应用1501。由于冰浆流动不同于一般的两相流动，它在输送的长管道内可能会出现分层，通过聚集而凝结成冰块，甚至会发生冰堵现象。所以冰浆的输运过程是一个急需解决的问题，国内外围绕冰浆的制取、输送及换热做了很多方面的研究。KITANOVSKI等采用了数值模拟结合实验的方法研究了刮削式制冰浆法。ROMAN等圆用过冷法进行了冰浆制取实验。华中科技大学的肖凯网对冰浆的制取方式进行了实验研究。清华大学的曲凯阳3）成功实现过冷水连续制取冰浆。浙江大学王炜I研究了一种双直接接触式冰浆制取系统。河南科技大学梁坤峰等借

10、鉴了流化床顺利进行了流体冰的制取实验。天津商业大学刘圣春等3】设计出一台小型的风冷式冰浆制取装置。华中科技大学的杨丽媛等电对冰浆在三通中的流动进行了数值模拟。中南大学徐爱祥等I以对间接式冰浆制取过程中冰层生长速率与冰层厚度变化及其影响因素进行数值模拟研究。大连理工大学王继红等成分别对竖直管道和水平90。弯管内冰浆流体流动特性进行数值模拟研究。河南科技大学的杜军恒廿引建立了混合物模型来模拟研究冰浆在水平直管，90。弯管及离心泵内的流动状态，杨帆等Ml对动态冰浆在水平直管内的流动压降进行了研究。上述研究中，针对盘管内冰浆流动特性的研究较少，本文利用数值模拟的方法研究冰浆在风机盘管中的流动特性，分析

11、流速对其流动特性的影响。1物理模型的建立、网格划分及边界条件设定1.1几何模型盘管的几何模型如图1所示，设计时完全参考常用空调末端设备，以家用分体式空调柜机蒸发器盘管为依据。设计盘管为3排铜管串联，管间距为20mm,铜管的管壁厚为0.35mm,直管部分管长为30nim。肋片为0.15mm的厚铜片增加散热效率，肋片间距设为1cm。含冰率为10%、盘管管径都为10min,其他条件不变，分别对入口速度为0.1m/s、0.25m/s0.5m/s、0.75m/s和1m/s的冰浆在管内流动进行了数值模拟。图1管径为10mm盘管模型图1.2网格划分网格划分在Fluent模拟中起着至关重要的影响，网格划分的效

12、果将直接决定模拟计算的速度快慢及效果。在本模拟中，由于冰浆流体具有粘性，本文首先对边界层进行了划分，然后以四边形非结构性网格形式对面网格进行了划分，最后整体网格划分为以六面体网格为主，在适当位置上包括楔形网格形式的非结构性网格。文献11对具有类似结构的弯管模型进行了网格独立性验证，网格数最多为34万时己经独立，本文在计算时网格数量约为40万，盘管划分的网格如图2所示。魂枷臾QWE顺顾箫加血油皿血f疽.训。购血祯阿皿琪（口贰明心血现血H棚M顷血此0.9n?/h、1.1m3/ho结果如图3、图4所示。换热量的计算公式如式(1)所示。Qw=WCpw(tw2-twl)+W(A/PF)L(1)式中：Qw

13、冰浆换热量，kw；W冰浆流量，kg/s；Cpw水的定压比热，KJ/(KgLC)；ZwH，w2分别为冰浆进出口温度，C；1PF浆融化的含冰率，%；L冰的融化潜热，335KJ/Kgo图3盘管流体进出口温差模拟值与实验值对比R制屈口王1.1流量/(mJ/h)图4换热量模拟值与实验值对比由图3和图4可以看出，在同等情况下，冰浆在盘管中的换热情况的模拟结果与实验结果有相同的规律性及数值吻合性。模型的建立比较符合实际情况、较准确，说明使用此种方法数值模拟冰浆在盘管内流动及换热结果接近实际并具有可参考性。4模拟结果及分析4.1入口流速对管道内冰粒分布的影响图5为不同流速条

14、件卜.管内冰粒子分布图。(a)速度为0.1m/s管道剖面冰粒子分布(b)速度为0.25m/s管道剖面冰粒子分布(c)速度为0.5m/s管道剖面冰粒子分布(d)速度为0.75m/s管道剖面冰粒子分布(e)速度为1m/s管道剖面冰粒子分布图5不同流速条件下管内冰粒子分布从图5(a)可以看出，当流速为0.1m/s时，除刚进入管道内部分，其他部分未融化的冰粒子儿乎全部聚集在管子顶部，特别是刚进入管道部分聚集程度较高。由于冰品粒了的密度低于水溶液的密度，出现冰粒了悬浮在管道顶部而滞留的现象，这与华中科技大学何国臾教授等mi总结的冰浆流动的3层模型完全吻合，冰浆流体在管道内可能呈现出3种分层：1)部分冰晶

15、粒子在顶部滞留，形成滞留层：2)往下形成冰粒子的整体移动层；3)最下层为粒于分布不均的混合溶液层。滞留粒子层会导致冰浆流动阻力增大，不利于冰浆在管道内的输送。从图5(b)可以看出，当流速为0.25m/s时，未融化的冰粒子并没有完全悬浮在管道顶部而滞留，只有部分冰粒子聚集靠近管道顶部，形成冰粒子维体移动层。从图5(c)和图5(d)可以看出，当流速大于0.5m/s时，未融化的冰粒子几乎没有出现聚集在管道顶部而悬浮的现象，而以小结块儿的形式分布在管道中间。从图5(e)可以看出，冰浆流速达到1m/s时，冰浆流动时冰晶颗粒与水相容得比较好，没有明显结块或者冰晶粒子聚集现象，冰浆流动属于混合溶液层。这是因

16、为，当冰浆流体流速比较大时，所形成的紊流强度也会比较大，这种情况下颗粒受到的紊流切应力远大悬浮力，因此颗粒在管道内基本对比图5(a)至图5(e)可以看出，当平均流速减小时，降低了紊流的影响，在管道顶部出现粒子高浓度区，若平均流速进一步减小，才可能在顶部出现滞留粒子层。只有冰浆在流速为0.1m/s的时候会出现在弯道后滞留层小部分聚集。4.2入口流速对换热的影响表1为不同流速的冰浆进出含冰率及融化了的含冰率对比。从表中可以看出流速为0.1m/s的冰浆融化了的含冰率最高，说明换热效果最好，随着流速的增加，出口含冰率增加，融化了的含冰率明显减少，即冰浆在盘管中流动无法很好地释放潜热，虽然冰粒子分布均匀

17、，流动阻力小，却失去了在末端设备流动换热的意义，所以高流速流动比较适合冰浆输送运输。成均匀悬浮状态。表I不同流速冰浆融化含冰率根据模拟得到的数据，结合冰浆换热量的计算公式得出不同流速条件下盘管的进出口温差、换热量及压降变化分别如图6、图7和图8所示。冰浆流速/(m/s)0.751入口含冰率APR10%10%10%10%10%出口含冰率/IPF26.02%7.37%8.36%8.87%9.23%融化含冰/AIPF3.98%2.63%1.64%1.13%0.77%p.、代噩口壬归-冬p.、代噩口壬归-冬流速/(m/s)图6不同速度流体进出口温差1(H)9(1SCI7(16(15

18、(10.20.4000.8I冰浆流速/(m/s)图7不同速度冰浆总换热量对比图8不同速度冰浆流动压降对比从图6中可以看出，在管道半径、冰浆流体含冰率不变的条件下，冰浆的进出口温度差随着流动速度的增加而减小，但减少程度小流速增加程度。流速0.1m/s的冰粒融化量最大；当流速为0.75m/s以上时温差接近2C,冰粒融化量较小。从图7可以看出，冰浆的总换热量随着流速的增加而增加，虽然含冰率随着流速的增大而减小，但其减小程度远小于流速增大的程度。从图8可以看出，在管道半径、冰浆流体含冰率不变的条件下，冰浆的压降随着流动速度的增加而增大，速度越高，增大幅度越大；这是由于，随着冰浆流体流速的增大，形成的紊

19、流强度也随之增大，由此产生的能量损耗也增加，从而压力损失也增大。5结论1) 冰浆在盘管流动过程中，弯管部位都会出现与单相液体流动一样的二次流，速度越大，二次流现象越明显。2) 入口速度对冰浆在盘管中的流动换热有很大影响，流速过小，冰浆中冰粒子聚集在管道顶部形成滞留层现象比较严重，不利于冰浆流动及换热，但冰浆流经管道持续的时间较长，换热时间长，有利于冰粒子的相变释放潜热；流速过大时，冰粒子在管道内基本成均匀悬浮状态，是冰浆最好的流动状态；但冰浆流经管道持续的时间太短，很多冰粒子来不及相变释放潜热就流出了盘管，换热效果不好;由于流体流速大，冰浆的压降也较大，所以冰浆在末端设备流动需要控制入口速度不

20、应超过0.5m/s,本文得出的比较理想的入口速度为0.25m/s左右。参考文献：1 梁坤峰，任幌乐，贾雪迎，等.冰蓄冷辐射空调系统运行优化与能耗分析J.制冷技术,2014,34(6):68-74.2 李叶，刘圣春,饶志明，等.冰浆储存与融化及流动换热特性研究现状及展望J.制冷技术，2012,40(1):55-60.|3|李新，侯予，张兴群.冰浆流动及传热特性研究进展J.制冷与空调,2007.8(7):15-18.|4)简夕忠.二元冰浆流动与传热特性研究D.南京：南京航空航天大学,2008.5 GRAYBF,JOHNSONCA,SCHOENAUGJ,ctal.Energyconsumptiona

21、ndeconomicevaluationofthermalstorageandrecoverysystemsforalargecommercialbuildingfJ.ASHRAETransactions,1988,94(1):412-424.6 WANGMJ,KUSUMOTON.Iceslurrybasedthermalstorageinmullifunctionalbuildings(J).HeatandMassTransfer.2001,37:597-604.71KITANOVSKIA,POREDOSA,Concentrationdistributionandviscosityofice-slurryinheterogeneousflowJ.InternationalJournalofRefrigeration,2002,25:827-835.81ROMAND.G