换热器核心算法的设计与实现.docVIP

（3）换热器模型 换热器模型是整个系统仿真软件的核心模块，从热力学的角度看，换热器直接决定了制冷系统的能量平衡，对系统性能有着至关重要的影响。由于换热器流路设计、风道结构、翅片形式、传热管内为两相流动换热等各种复杂因素，要纯粹在理论上建立一个精度较高的换热器模型是异常困难的，因此我们通常都是依赖于半经验的方法来建立制冷系统中的换热器模型。 换热器模型目前主要有两个典型代表，分别是以美国Purdue大学为代表的ACMODEL，以及以美国国家橡树岭实验室为代表的MARK模型。 在ACMODEL仿真模型中，针对换热器每个流路的铜管进行分段计算，即每个流路的铜管分为多个管段，针对每个管段计算传热量、流动阻力以及每个管段内的制冷剂质量。上一个管段的制冷剂出口状态条件作为下一个管段的制冷剂入口条件，可以根据制冷剂的参数（比如流量、干度）等确定管内的流态，从而选择相应的换热关联式，如图2-1所示。上游两个管段的空气出口状态条件的平均值作为下游对应管段的空气入口状态条件，如图2-2所示。 图1对每个流程的铜管进行分段计算 图2 空气侧分段计算的进口状态参数确定 对于ACMODEL仿真模型，首先需要利用输入文件清晰的定义好分配头、流路的逻辑关系，它可以指定每个流路的制冷剂流量和风量的分布，因此可以非常灵活的依据试验数据来建立各种修正系数。同时，由于每个管段的制冷剂状态都可以计算确定，比如干度和密度，所以非常容易就可以算出每个管段的制冷剂质量，从而能够模拟制冷剂充灌量对制冷系统性能的影响。理论上该方案无疑是最合理和准确的，但由于换热和流动的计算关联式本身就不精确，特别是两相流动状态下的关联式，通常文献报道的误差都达到了20%，在具体的应用场合下误差就更大，而且空气在换热器内部的流场分布也存在一定的不可知性，同时该方案的程序开发的工作量巨大。 在MARK换热器仿真模型中，他们把整个换热器作为一个整体来进行计算，只是采用了滑动界面模型来区分单相区（汽相或液相）和两相区，对于不同的区域选用不同的换热关联式，通过模型的不断迭代计算，确定每个区域的换热面积，从而计算出相应的换热量和出口状态。该算法只能对风侧和制冷剂侧的整个传热系数增加修正因子，无法模拟复杂的换热器流路设计。 Mark模型尽管理论上存在一定的缺陷，但根据BO SHEN的研究，经过修正后其仿真精度并不比ACMODEL模型差，而且可以大大减少整个软件的开发工作量，因此本项目的工作主要是在Mark模型的基础上进行。为了在Mark模型的算法基础上能够模拟充灌量对系统性能的影响，我们的思路是，以某台机型的实验数据作为基准点，制冷剂充灌量对性能的影响可以折算为内部有效气体容积变化对压力的影响以及换热器中液体段长度变化对换热的影响。 4 冷凝器模型与算法 冷凝器模型的计算流程如图所示，其中涉及的重要计算包括：等。下面将计算详细介绍。 图 （1）冷凝器进、出口参数 冷凝器进风温度等于室外干球温度，即已知。空气出口温 （11） 式中Va，c为冷凝器的风量，ρa和Cpa则分别为冷凝器冷却空气的密度和热比容。 对于制冷剂侧的状态参数计算，可参照第二章的屋顶机热力学分析。 （2）冷凝器结构计算 图6 套片管结构示意图 近年来，为了强化空气侧的换热，出现了各种形状的强化传热翅片，例如波纹片、百叶窗片、开缝片等等，约可提高空气侧换热系数20～60%。对于该种套片型的换热器，其基本几何参数计算可按下述方法进行。 每米管长翅片的外表面积 ： （12） 式中的S1和S2分别为换热器的管间距和排间距，db为管基外径，Sf为翅片间距。 每米管长翅片间表面积： （13） 式中δf为翅片厚度。 每米管长的总外表面积： （14） 每米管长平均直径处的表面积： （15） 式中，其中di为换热管的内径。 每米管长内表面积 ： （16） 肋化系数： （17） 冷凝器最窄流通面与迎风面之比： （18） 换热器的迎面风速wf： （19） 式中的qair为冷凝器风量，L1和L2则分别为整张翅片的宽和高。 在换热器最小截面处的最大速： （20） （3）空气侧换热系数和压降计算 对于平直翅片，可以采用下式计算空气侧的换热系数ao： 、、、为系数及指数，可分别参见表3和4。为空气的导热系数， L1为翅片宽度，和Rea则为当量直径和最大风速雷诺数，可分别用下式计算：