空气源热泵冷热水机组空气侧换热器结霜工..
详细内容
摘要 采用空气源热泵冷热水机组的动态数学模型对空气侧换热器的结霜工况进行了模拟。模拟中同时考虑了结霜的密度和厚度随时间的变化，首次提出了结霜密度随时间的变化关系式。计算了不同工况下的结霜速度、霜的密度、霜的厚度随时间的变化。将模拟结果与实验数据进行了比较，进一步验证了所建模型的正确性。
关键词 空气源热泵冷热水机组 结霜 动态模拟
1 前言
空气源热泵冷热水机组作业中央空调的冷热源有很多优势，如冬夏共用，设备利用率高；省去了锅炉房和一套冷却水系统；机组可安装在室外，节省了机房的建筑面积；不污染环境等。因此该机组在气候适宜地区的中小型建筑中得到了广泛地应用。但机组在冬季运行时，当空气侧换热器表面温度低于周围空气的露点温度且低于0℃时，换热器表面就会结霜。结霜后换热器的传热效果急剧恶化，严重时机组会停止运行。因此换热器结霜是影响机组应用和发展的主要问题，研究机组在结霜工况下的工作性能具有十分重要的意义。
2．结霜模型的建立
霜的积累速率 是由进出室外换热器空气湿度的变化决定的：
（1）
式中： ----空气的质量流量，kg/s；
di，d0----分别为空气进、出换热器的含湿量，kg/kg。
由于霜的多孔性和分子扩散作用，在表面温度低于0℃的换热器上沉降为霜的水分一部分用以提高霜层的厚度，一部分用以增加霜的密度[1]，即
（2）
式中用于霜密度变化的结霜量变化率 由下式确定[2]：
（3）
式中： ----换热器的全热交换量，W；
iSV----水蒸气的升华潜热，J/kg；
λfr----霜的导热系数，W/（m・K）；
R----水蒸气的气体常数，461.9/（kg・K）；
TS----霜表面的温度，K；
pV----水蒸气的分压力，Pa；
vV,vi----分别为水蒸气、冰的比容，kg/m3。
ρfr，ρi----分别为霜、冰的密度，m3/kg；
DS----霜表面水蒸气的扩散系数，m2/s。
而霜的密度ρfr与换热器表面的温度、空气的温度、相对湿度、流速和结霜的时间等有关，结霜时间越长，霜的密度越大。计算时，先假设一个初始密度，由下式计算霜的导热系数，再计算霜密度和厚度的变化。
（4）
对于每一个时间步长Δt，霜密度的变化和厚度的变化为：
（5）
（6）
式中：At----换热器的总换热面积，m2；
δt----霜层的厚度，m。
3 模型的求解
我们对空气侧换热器后个换热单元在不同工况下的结霜情况进行了模拟计算，该单元的结霜情况可以反映出整个换热器的结霜情况。空气侧换热器由160个这种换热单元组成。计算的换热器单元结构参数见表1，计算工况见表2。
换热器单元的结构参数 表1
空气源热泵冷热水机组空气侧换热器结霜工况的动态模拟 : 图5和图6为动态工况下霜厚度随时间的变化。图5为空气温度一定（0℃）时，不同相对湿度（65%、75%、85%）下霜厚度的变化。由图可见，随着时间的增加，霜的厚度迅速增加，而且相对湿度越大，霜厚度增加越快。在该计算工况下，霜厚度在到0.5mm左右时，应开始融霜。
图5 不同相对湿度下霜厚度的变化
图6 不同温度下直厚度的变化
图6为相对湿度一定（75%）时，不同空气温度（0℃、-4℃）下霜厚度的变化。由图可见，0℃，75%工况（工况B）下，运行60分钟左右就需要融霜，而-4℃、75%工况（工况E）下，则运行115分钟时才需融霜。
显然，空气源热泵冷热水机组除霜控制方法常用的时间控制法和时间-温度控制法是不符合霜厚度随时间的变化规律的。如当机组设定的固定除霜时间按工况C确定时，那么工况B和工况A将会出现不必要的除霜，从而影响了机组的效率。同样，许多生产厂家虽采用时间-温度控制法，但还是采用统一固定的除霜启动值和除霜时间值，因此由于空气温度、相对湿度的不同，结霜的厚度不同，除霜效果也就不一样。结霜规律的正确预测，才是保证除霜效果良好的前提。
4 结论
空气侧面换热器结霜过程中，不仅霜的厚度发生变化，霜的密度也在发生变化，本文同时考虑了结霜的密度和厚度随时间的变化，并根据一些实验数据和结霜密度的变化规律，首次提出了用于霜密度变化的结霜量变化率随时间的变化关系式，并认为在刚开始结霜时，结霜量主要是增加霜的厚度，而密度变化很小。随着时间的推移，霜的厚度增加减缓，而密度变化增加，而且霜的密度随时间呈抛物线规律变化。
分别计算了不同工况下的结霜速率、霜的密度、霜的厚度随时间的变化。计算结果表明，在不同的工况下，空气侧换热器的结霜情况是不同的。在空气温度一定时，相对温度越大，结霜越严重，融霜的时间间隔越短；在空气相对湿度一定时，0℃工况的结霜比-4℃工况的结霜严重。而且计算出了不同工况下融霜的时间间隔，为采取有效的除霜控制方法提供了依据。
将模拟结果与实验数据进行了比较，两者吻合很好，进一步验证了所建模型的可靠性。
参考文献
1 S.N. Kondepudi, D.L.O Neal. Performance of Finned-Tube Heat Exchangers under Frosting Conditions: I. Simulation Model. Int. J. of Refrig. 1993, 16 (3) :175～180
2 姚杨，姜益强，马最良，空气源热泵冷热水机组空气侧换热器结霜规律的研究，哈尔滨工业大学学报，2002（6）
3 T. Senshu, H. Yasuda. Heat Pump Performance under Frosting Conditions: part I: -Heat and Mass Transfer on Cross-finned Tube Heat Exchangers under Frosting Conditions. ASHRAE Trans, 1990, 96 (1): 324-329
4 T.S.Gatchilov. 翅片式空气冷却器表面结霜的特性，第十五届国际制冷大会论文译丛（中册），制冷学报编辑部，1981：356～362。
空气源热泵冷热水机组空气侧换热器结霜工况的动态模拟 :