辐射供冷空调技术于20世纪80年代兴起于欧洲,在中国起步较晚[1]。许多研究表明,辐射空调系统比变风量空调系统更为节能[2-4]。在目前可持续发展的战略背景下,辐射供冷系统成为一个绿色低能耗建筑首先考虑的空调系统,具有良好的发展前景[5]。主动式冷梁技术作为辐射供冷的空调末端技术,在中国发展迅速。由一些案例中得知,主动式冷梁系统仍有节能潜力尚未挖掘[6-8]。近年来, 人们对冷梁的诱导器结构和运行系统做了大量的研究和改进[9-12],但目前在对冷梁的研究中尚无针对主动式冷梁建立仿真模型。本文在简化假设的前提下,对主动式冷梁诱导送风规律和传热规律进行了分析,为主动式冷梁建立了仿真工作模型。对主动式冷梁样品风侧诱导性能和换热器换热性能进行了实验研究,将实验结果与理论计算进行对比,并对二者的差异性和一致性进行了分析。
主动式冷梁工作原理如图 1所示。以供冷工况为例,经热湿处理后的一次风送入诱导器,在诱导作用下产生二次风(室内热空气),二次风在换热器处换热后温度降低,和一次风混合后送入室内,达到降低室内温度的目的[13-14]。基于此,将主动式冷梁分为一次风诱导过程和二次风的换热过程两部分分别研究。
一次风送入主动式冷梁静压箱后,经喷嘴高速喷出,通过诱导作用卷吸室内二次风通过换热器发生换热作用,该诱导作用发生在冷梁的诱导器处,静压箱则起到均压的作用。如图 2所示,冷梁喷嘴直径为d,喷嘴中心距侧板的距离为l1,喷嘴中心距挡板的距离为l2。
以单个喷嘴诱导送风为研究对象,过喷嘴中心轴线且垂直于图面的平面将射流流场分为两部分,分别记为部分Ⅰ和部分Ⅱ,射流流场与侧板和挡板搭接时的射程分别为s1和s2,喷嘴出口处的风量为Q1,部分Ⅰ和部分Ⅱ的风量分别记为QⅠ和QⅡ。由无限空间圆断面淹没射流理论,分别计算部分Ⅰ和部分Ⅱ的流场风量。
考虑到主动式冷梁诱导腔的诱导作用是非常复杂的,在对此处诱导送风理论进行研究过程中进行适当的简化,同时, 保证可靠性和准确度。假设条件如下:
1) 忽略侧板及挡板对射流流场的影响,将受限空间射流理想化为自由射流,射流流场断面为圆形。
2) 当射流与挡板和侧板刚刚搭接后,或相邻喷嘴射流流场边界层刚好搭接后,就不再产生诱导效果。
3) 忽略相邻射流对流场的影响,即相邻射流互不混合,搭接后紧密流动。
4) 忽略表冷器侧的阻力对二次风量的影响。
5) 假设射流不射入换热器上方区域,即射流完全射入冷梁挡板和侧板所围城的夹缝区域。
诱导送风过程中诱导能力的大小关系到二次风量Q2的大小,诱导能力越强,卷吸的二次风量越多,相同条件下产生的换热量越大。本文引入诱导比的概念来衡量诱导能力的大小,定义主动式冷梁的诱导比为二次风量与一次风量的比值
单排多喷嘴诱导射流,如图 3(a)所示,设喷嘴等距离分布,间距为l3,相邻喷嘴搭接时射流的射程为s3。图 3(b)为射流流场沿射流方向的正视图,流场截面被分为如图所示4部分,分别命名为Ⅰ′、Ⅱ′、Ⅲ′和Ⅳ′。Ⅰ′部分对应的射程为s1,Ⅱ′部分对应的射程为s2,Ⅲ′和Ⅳ′两部分为对称结构,可视为一体,对应射程s3,由此将射流诱导过程分为三部分。分别计算以上四部分射流的诱导比,然后取均值,即认为是受限空间内多喷嘴射流诱导比。
sn为无限空间淹没紊流射流条件下起始段核心长度[15]。对于部分(i=1, 2, 3, 4),当si≤sn时,部分的理论诱导比可按式(2)[15]计算
当si > sn时,部分的理论诱导比可按式(3)[15]计算:
式中:a为紊流系数,由实验决定,是表示射流流动结构的特征系数。
综上,主动式冷梁理论诱导比可表示为
主动式冷梁表冷器供冷工况的供水温度应高于室内回风温度对应的露点温度,以防止结露的发生,即表冷器处于干工况运行模式,整个供冷过程不产生冷凝水。主动式冷梁目前多采用管翅式表冷器作为换热元件,铜管铝片结构能够使换热器达到较高的换热效率,室内二次风由于一次风的诱导作用,经过管翅式换热器后温度逐渐降低,与一次风混合后送入室内。
二次风与表冷器换热过程可看作单相流体外掠翅片管束强制对流换热模型。翅片管束可以视作圆管外加等截面环肋模型处理。对于实际的管翅式换热器,管间距为s1,排间距为s2,换热器铜管铝片通过铜管胀管方式紧密连接。设tf1 > tf2,则在稳态传热工况时,通过表冷器的传热量可表示为如下关系式:
通过管内壁换热量
通过铜管和管外包裹铝片层壁的换热量
通过肋片换热量(肋与肋间换热之和)
结合肋片效率η[16],式(7)可改写为
式(8)中,翅片间铝片与空气换热的表面传热系数h′2可通过单向流体外掠管束传热实验关联式计算得到,翅片与空气换热的表面传热系数h″2可通过单相流体外掠平板传热计算式计算得到,对于主动式冷梁运行工况而言,h′2计算值小于h″2,且翅片面积远大于翅间铝片的面积,因此,可忽略h′2A′2,将式(8)进行简化处理
联立式(5)、式(6)和式(8′),消去tw1和tw2得到
上式以沿管长方向单位长度的换热器为基准。
将诱导送风模型与表冷器干工况换热模型相结合,即得到主动式冷梁换热模型。连接诱导送风模型与表冷器干工况换热模型的纽带为二次风。由对诱导送风模型和表冷器干工况换热模型的分析可知,当主动式冷梁结构参数一定时、主动式冷梁输入的一次风量一定时,可通过式(1)和式(4)计算得到主动式冷梁诱导的二次风量;通过二次风量可以计算得到主动式冷梁表冷器翅片表面对流换热系数,已知供水参数和表冷器详细结构参数的情况下,可以通过式(9)计算得到主动式冷梁水侧理论换热量。由此可得到主动式冷梁水侧换热量与一次风量的关系式。
主动式冷梁工作时,表冷器迎面风速不大,一般不超过1 m/s,翅片表面流场处于层流。由层流外掠平板传热计算式可计算得到翅片与空气换热的表面传热系数h″2,带入式(9)得到主动式冷梁一次风量与水侧换热量关系式:
在实际计算中tf1取供回水温度的算术平均值,tf2取二次风进出换热器的平均温度。
对某主动式冷梁样品进行了热工性能测试。测试内容包括静压箱静压、一次风量、二次风量、水侧换热量等。实验过程中,设定5种冷梁静压箱静压值分别为20、40、60、80、100 Pa,在以上几种静压条件下,测定冷梁其他热工参数值。实验测试的主动式冷梁结构参数详见表 1。
根据本文建立的诱导送风模型,计算得到的主动式冷梁诱导比理论值与实测值如表 2所示。
由表 2得知,测试冷梁的诱导比实测值与理论值相差不多,模型计算值与实验测试值符合性良好,则可将其作为主动式冷梁工作模型的验证,对结果进行相应的分析。由于实验过程中未对经过表冷器后的空气温度进行测试,根据热力学第一定律,本文通过水侧换热量计算得到经过表冷器后的空气温度,进而求得二次风进出换热器平均温度tf2。根据供水温度及流量可求得冷水与铜管换热的对流表面传热系数,根据换热器结构参数可求得铜管及铝片的热阻。在计算肋片侧热阻的过程中,将表冷器等效为环肋模型,即将翅片等面积分到每根铜管上,建立圆管环肋换热模型进行计算。将求得的参数和一次风量及理论诱导比等参数带入式(10),即可求得主动式冷梁工作时的理论换热量。
实验分别测试了冷梁在4种喷嘴间距时的换热量,在不同静压箱静压和换热温差实验条件下,共进行了300组实验,选取其中静压箱静压值为40 Pa与80 Pa的数据进行分析。图 4为测试冷梁喷嘴间距不同时,在不同换热温差下的理论制冷量与实际制冷量对比图。
由图 4可以看出实际与理论制冷量曲线在不同喷口间距情况下走势出入不大,证明了建立的仿真模型在实际应用中具有一定的准确性。但理论值与实际值之间的误差则需要做进一步的误差分析。
使用完全的实验数据与模型计算结果进行比对,模型计算值与实验值的相对误差保持在15%以内,各实验工况下计算值与实验值误差分布规律如图 5所示。由图 5可知,模型计算值与实验值的最大相对误差在-10%~10%范围内的点占总数据点的比例为77%,说明模型计算值与实验测试值差异性不大。
采用Bland-Altman方法对各工况实验值和模型计算值进行一致性分析,如图 6所示。Bland-Altman方法是一种比较常用的用于分析模拟数据和实验数据一致性的方法,它可以把两种不同方式得到的数据的定性和定量分析结合到一起,用图表的形式形象的反映出来。Bland-Altman方法不仅包含了数据得到的过程的随机误差,同时也将系统误差纳入到了考虑的范围内,能够比较全面准确地反映不同方法得到的数据的一致性。限于篇幅,文中未列出各工况具体的实验条件及实验结果。
由图 6可以明显地看出两种数据的一致性较好,大部分数值点落在在95%置信区间范围内。经计算,超出95%的置信区间的数值点占总数值点的比例为5.67%。对于这样的模型计算结果,是可以认为模型计算值和实验测试值具有较强的一致性的,说明主动式冷梁工作模型在文中所列约束条件下具有较强的适用性。
在简化条件下建立了诱导送风模型和表冷器换热模型,在两个模型的基础上最终建立起了主动式冷梁工作模型。实验测试了某主动式冷梁样品的热工特性,验证了所建模型的准确性和适用性。