部分负荷运行时风冷热泵机组特性模拟分析.docVIP

部分负荷运行时风冷热泵机组特性模拟分析 　　【摘要】建立了热泵空调器制冷系统风冷热泵在结霜工况下稳态分布参数模型和结霜过程分布参数模型，编制了适于计算机上运行的matlab模拟程序；模拟了室外换热器结霜引起风机流量下降及结霜过程中与制冷循环参数的相互影响；并以压缩机的台数为调节参数，改变机组输出，根据模拟结果分析了热泵空调器部分负荷时结霜工况下机组的性能，并对比了一种优化方案与常规方案结霜工况的性能，通过对结霜工况下机组满负荷和部分负荷时性能的数值模拟,结果表明该优化方案提高了热泵机组部分负荷运行时的COP值，有很好的节能效果。 　　【关键词】热泵;结霜;部分负荷;性能 　　 　　0.前言 　　大型集中空调的热泵机组在部分负荷运行时，一般采用多台分级调节的控制方案（例如美国某公司及其上海分公司完成的上千个大型集中空调项目中，有86%的项目采用了两台或两台以上的热泵机组组合运行方案），受传统调节控制模式的制约，制冷系统配件(如热力膨胀阀)存在缺陷，机组运行的安全要求(回油)不到满足等等。因此为了使机组具备调节能力，常常由多台压缩机组成一台机组，每个压缩机构成一个相对独立的制冷回路，运行时保持机组的冷凝压力恒定(或在一个较小的范围内变化,根据负荷情况启停压缩机，启动后的压缩机都满负荷运行。这种多台分级调节方法虽然能满足部分负荷的使用要求，但冷量的调节不是连续的，而且还造成了部分机组闲置，设备的潜能不能得到充分发挥。众所周知，机组实际运行性能系数COP是由机组内部(制冷循环、压缩机效率、传热效果等)因素和外部条件共同影响的，而通常冬季部分负荷运行时，外部条件较满负荷时优越,环境温度相对较高。在不考虑传热性能变化影响的情况下，可充分利用环境温度的升高来增加蒸发温度温度、缩小压缩机的高低压力差以改善系统性能系数。本文提出一种新的优化运行方案，该方案从制冷系统入手，利用已有设备,通过增大实际运行机组的传热面积来提高热泵机组部分负荷运行时的COP值以实现节能，且可以实现能量的连续调节。 　　1.空气源热泵空调器系统结霜数学模型 　　1.1压缩机模型、冷凝器模型、电子膨胀阀模型 　　压缩机模型、冷凝器模型、电子膨胀阀模型详见[1]. 　　1.2 蒸发器结霜模型 　　1.2.1 蒸发器模型空气侧控制方程 　　没有结霜现象时，空气侧的控制方程为： 　　(C.△t+r.△d)dM=a(t-t).F.dx（2.1） 　　式中：F＝A+A 　　 C――空气定压比热,kJ/(kg#8226;K) 　　 r――水蒸气汽化潜热,kJ/kg 　　 M――空气质量流量,kg/s 　　 a――空气侧换热系数,W/(k#8226;m) 　　 ――析湿系数 　　 F――每米管长翅片表面积,m 　　 t――空气温度， ℃ 　　 t――翅片表面温度， ℃ 　　 d――空气湿度,kg/kg 　　 n――翅片效率 　　结霜后，霜层在翅片表面沉积,增加了换热热阻,空气侧的换热方程变为： 　　(C.△t+r.△d)dM=K(t-t).F.dx （2.2） 　　式中：K――霜层表面的换热系数，W/(k#8226;m) 　　结霜后空气与管壁之间的总换热系数 　　K= （2.3） 　　式中： 　　 a――空气侧换热系数,W/(k#8226;m) 　　 F――每米管长结霜的翅片的外表面积,m/m 　　 ――霜的厚度， m 　　 ――霜层的导热系数，W/(k#8226;m) 　　1.2.2 霜层参数的计算 　　(1)霜的密度与导热系数 　　霜的密度与导热系数采用Hosuda和Uznhashi[2]的经验公式，公式中将霜的特性与外界空气的参数连系起来： 　　 =340t+85v（2.4） 　　式中： 　　 ――霜层密度, kg/m 　　 t――肋片表面温度,℃ 　　 v――空气比容,m/kg 　　 = （2.5） 　　（2）霜的厚度 　　在时间 内增加的霜层的厚度： 　　△= （2.6） 　　（3）翅片效率 　　当翅片表面结霜时，由于霜的绝热效应，使其表面温度趋于一致，因此，翅片效率有所提高。本文采用Sanders[3](1974)的结霜工况下的翅片效率模型： 　　 =（2.7） 　　式中： 　　 m：结霜量； 　　 m= 　　 h：当量肋高 　　 b：常数,是拟合直线关系式h=a+bT的斜率，由湿空气物性插值求出。 　　（4）结霜后蒸发器进口风速的变化 　　换热器翅片表面结霜后其空气阻力增加，进而影风机的工作特性。 　　a. 空气流动阻力 　　本文采用Hosada及Uzahashi给出的结霜工况下阻力经验公式[4]。 　　P=5.88×10N 　　 ×.v（2.8）