管式蒸發冷却器冷却性能的数值求解.docVIP

管式蒸发冷却器冷却性能的数值求解 山东建筑大学 刘乃玲 济南同圆建筑设计研究院 邵东岳 摘要：根据热湿交换理论建立了管式蒸发冷却器的数学模型，采用四阶Runge-Kutta方法对该模型进行了数值求解，并将计算结果与实验数据进行了对比，结果表明计算值与实验值吻合较好，证明用数值分析的方法研究类似的工程问题是可行的。分析了运行参数对蒸发冷却器冷却性能的影响，并得到了一些有价值的结论。 关键词：管式蒸发冷却器；数学模型；数值求解 引言 管式蒸发冷却器是一种常用的冷却设备，与板翅式间接蒸发冷却器相比，具有布水均匀，流道较宽，不会产生堵塞现象，流动阻力小等优点，因而管式蒸发冷却器逐渐被人们重视起来。闭式冷却塔就是管式蒸发冷却器的一种应用形式。它的主要原理是喷淋水在盘管外壁上蒸发以冷却管内流体，并利用风机及时把产生的水蒸汽带走，以提高冷却效果。闭式冷却塔与开式冷却塔相比由于被冷却工质在盘管中流动，避免了冷却工质与空气的直接接触，管内工质不易被污染。常规机械制冷空调系统中若采用闭式冷却塔，在过渡季节可以将冷却水与冷冻水系统进行切换，从而实现免费供冷，可大大节省制冷机的运行费用[1,2]。因此闭式冷却塔在空调节能方面有着广阔的应用前景可分为顺流逆流两种形式[3]，计算模型微元控制体如图1所示： 管内流体失去的热量为： (1) 即 (2) 式中：mc为管内流体流量，kg/s；cpc为管内流体的比热，J/(kg·K)，tc为管内流体的温度,℃；K为从传热管内流体到管外喷淋水的总传热系数，W/（m2·℃）；tw为管外流体的温度，℃；A为蒸发冷却器的断面面积，m2； z为距蒸发冷却器底部的距离，m；为蒸发冷却器单位容积的传热管的传热面积，m2/m3。 空气得到的热量为： (3) 即 (4) 式中： ma为管外空气流量，kg/s；为管外喷淋水温度下饱和湿空气的焓，J/kg；为管外湿空气的焓，J/(kg·K)；hm为从管外喷淋水向空气的总传质系数，kg/(m3·s)；为蒸发冷却器单位容积的传热管的有效气液界面面积，m2/m3。 管外冷却水失去的热量为： (5) 取图1中虚框内的部分为计算控制体，所得能量方程为： (6) 略去二阶无穷小，整理得 (7) 式中： mw为管外喷淋水流量，kg/s；为管外喷淋水的焓，J/kg；为管外喷淋水温度，℃。 由于蒸发冷却器喷淋水的蒸发率很小[4,5]，因此可近似取，即。将式（1）、(2)、(4)代入式(7)整理得 (8) 联立式（2）、（4）、（8）可得到方程组 (9) 为使方程有唯一解，需确定方程组的边界条件，由于喷淋水是循环的，因此微分方程组满足，方程组(9)的边界条件为： (10) 令，，，，则微分方程组(9)变为： (11) 数学模型的求解 方程组(11)为线性齐次微分方程组，本文采用四阶Runge-Kutta方法求解[6]。方程组(11)可改写为如下形式： (12) 对于顺流式，沿冷却盘管高度方向逐层计算，每层间隔取管间距。分层示意如图2所示。 根据计算原理则有如下的方程组 (13) 其中 (14) (15) (16) (17) 从上面的方程组可知只有当前一层的所有变量已知时，才能求得下一层的所有变量，对于顺流式，计算第一层时未知，因此需假定该变量的初值，然后采用迭代法依次计算各层变量的值，当时，计算所得的各层变量值即为最终解。 逆流式冷却器数学模型的建立及计算方法与顺流式类似，计算结果可参阅文献[7]。 数学模型的实验验证 本文以Hasan[8]所研究的冷却器模型为对象进行了理论计算，并将计算结果与Hasan的实验结果进行了对比。Hasan研究模型的运行参数及结构参数如表1所示。 表1 蒸发冷却器参数表 运行参数 结构参数 /℃ /(kg/s) /℃ /% /(kg/s) /(kg/s) /m /m /m 23.96 0.40 26.19 47 1.27