空调换热器结霜管理

空调换热器结霜管理摘要采用空气源热泵冷热水机组的动态数学模型对空气侧换热器的结霜工况进行了模拟。模拟中同时考虑了结霜的密度和厚度随时间的变化，首次提出了结霜密度随时间的变化关系式。计算了不同工况下的结霜速度、霜的密度、霜的厚度随时间的变化。将模拟结果与实验数据进行了比较，进一步验证了所建模型的正确性。关键词空气源热泵冷热水机组结霜动态模拟1前言空气源热泵冷热水机组作业中央空调的冷热源有很多优势，如冬夏共用，设备利用率高；省去了锅炉房和一套冷却水系统；机组可安装在室外，节省了机房的建筑面积；不污染环境等。因此该机组在气候适宜地区的中小型建筑中得到了广泛地应用。但机组在冬季运行时，当空气侧换热器表面温度低于周围空气的露点温度且低于0时，换热器表面就会结霜。结霜后换热器的传热效果急剧恶化，严重时机组会停止运行。因此换热器结霜是影响机组应用和发展的主要问题，研究机组在结霜工况下的工作性能具有十分重要的意义。2结霜模型的建立霜的积累速率是由进出室外换热器空气湿度的变化决定的：式中：-空气的质量流量，kg/s；di，d0-分别为空气进、出换热器的含湿量，kg/kg。由于霜的多孔性和分子扩散作用，在表面温度低于0的换热器上沉降为霜的水分一部分用以提高霜层的厚度，一部分用以增加霜的密度【1】，即式中用于霜密度变化的结霜量变化率由下式确定【2】：式中：-换热器的全热交换量，W；iSV-水蒸气的升华潜热，J/kg；fr-霜的导热系数，W/对于每一个时间步长t，霜密度的变化和厚度的变化为：式中：At-换热器的总换热面积，m2；t-霜层的厚度，m。3模型的求解我们对空气侧换热器后个换热单元在不同工况下的结霜情况进行了模拟计算，该单元的结霜情况可以反映出整个换热器的结霜情况。空气侧换热器由160个这种换热单元组成。计算的换热器单元结构参数见表1，计算工况见表2。换热器单元的结构参数表1管材铜管径0.15风向管排数4迎风管排数20管间距S125.4管间距S222mm翅片材料铝片型波纹片片厚0.2mm片间距2.0mm翅化系数17.8单根管长16m分液路数10在求解结霜的动态模型时，必须考虑结霜的密度和厚度随时间的变化，但在以往的结霜量计算中，均未同时考虑结霜的密度和厚度随时间的变化。如.-【3】建立了空冷器上结霜的数学模型，并将模拟结果与实验数据进行了比较。由于假设霜层均匀分布，且霜的厚度随时间线性增加，而霜的密度不随时间变化，使得模拟霜的厚度比实验值大20%30%。计算工况表2工况编号空气温度1A0651062-1355050.0096B0751062-1355050.0096C0851062-1355050.00962D-4651062-1355050.00816E-4751062-1355050.00816F-4851062-1355050.00816本文根据一些实验数据和结霜密度的变化规律，首次提出了结霜密度随时间的变化关系式，并认为在刚开始结霜时，结霜量度要是增加霜的厚度，而密度变化很小。随着时间的推移，霜厚度的增加变缓，而密度变化增加，而且霜的密度随时间呈抛物线规律变化。由稳态模型和公式式中为结霜的时间，min。为验证所建的换热器结霜模型正确性，将模拟结果与实验数据进行了比较，我们采用文献【4】中的实验数据。实验是日本工业标准实验工况1.585-7.53.3模拟工况085-132.5由于实验工况与模拟工况换热器的换热面积不同，因此单纯地比较结霜量的变化是没有实际意义的。为此提出了单位换热面积结霜量的概念，即结霜量与总换热面积之比。实验工况与模拟工况的单位换热面积结霜量变化见图1。由图可见，模拟值与实验工况的条件略有差异造成的，因为模拟工况的蒸发温度比实验工况低，且迎面风速小，而蒸发温度越低，结霜量越多；迎面风速越低，结霜量也越多。这两方面的因素造成了模拟值略大于实验值。通过比较进一步验证了所建模型的正确性。图2为空气温度一定3T.Senshu,H.Yasuda.HeatPumpPerformanceunderFrostingConditions:partI:-HeatandMassTransferonCross-finnedTubeHeatExchangersunderFrostingConditions.ASHRAETrans,1990,96(1):324-3294T.S.Gatchilov.翅片式空气冷却器表面结霜的特性，第十五届国际制冷大会论文译丛（中册），制冷学报编辑部，1981：356362。