等价湿球温度法分析传热传质过程.doc

板式蒸发空气冷却器传热传质模型分析

摘 要： 根据热力学和传热学理论，建立了板式间接蒸发式空冷器传递过程的微分方程组。引入空气湿球比热和对流换热系数，推导得到了基于空气湿球温度的数学传递模型。应用四阶-龙格库塔法求解了高炉循环水用板式蒸发冷却器基于湿球温度的传递模型，得到了热流体、喷淋水、空气干湿球温度的变化规律和各流体的换热强度情况。该分析模型为进一步分析蒸发式空冷器的热工性能和优化设计提供了理论依据。
关键词：传热传质；板式蒸发空冷器；传递模型；湿球温度

蒸发式空冷器利用自然环境中空气的干湿球温差取得冷量来冷却热流体，由于它比空气冷却换热效果好，比水冷节水等显著优点，因此在石油、化工、冶金、电力等领域中有着广泛的应用。蒸发式空冷器冷却侧空气与喷淋水之间传热和传质过程同时进行，相互耦合，传递机理相当复杂。目前针对管式和翅片管式蒸发空冷器的研究成果较多[1-3]，然而这些成果大多是基于Merkel理论，空气和喷淋水的热湿交换以焓差为推动力建立数学模型。本文建立了基于湿球温度的板式蒸发式空冷器的数学模型，并用四阶-龙格库塔法对一个实例进行数值求解，分析了各流体的温度变化和换热规律。
1 传递过程的基本方程组
图1 逆流板式蒸发空冷器示意图
空气
喷淋水shui
热流体
本文以板式蒸发空冷器为研究对象，热流体与空气逆流，喷淋水与空气并流，流体组合如图1所示。对蒸发式空冷器的传热传质和流动过程做如下假设：
1. 空冷器内传热传质和流动过程处于稳态，忽略外壳的散热损失。
2. 各流体热物性参数为常数。流体的状态参数仅沿流动方向变化。
3. 喷淋水膜均匀分布，忽略水膜波动和水膜厚度对传递和流动的影响，忽略水膜的蒸发损失；水膜在传热壁面上完全润湿。
4. 忽略空气中离散水珠对传热传质的影响。
喷淋水
热流体

空气
水膜
壁面
dz
图2 微元面积传递过程
5. 湿空气的传热传质过程符合刘易斯关系式，即，为喷淋水与空气的传质系数，kg·m-2·s-1；为喷淋水与空气的传热系数，W·m-2·K-1；为空气比热容J·kg-1·℃-1。
取微元面积Bdz进行传热传质分析。热流体侧的能量守恒方程为：
（1）
式中：为热流体的比热容，J·kg-1·℃-1；为热流体的质量流量，kg·s-1；为空冷器板片数；为空冷器宽度，m；为热流体温度，℃；为喷淋水温度，℃；为热流体到水膜的总传热系数，W·m-2·℃-1。
因为水膜很薄，可认为气液界面的湿空气的饱和温度等于水膜的温度，则空气侧水蒸气的质量守恒方程为：
（2）
式中：为空气（干空气）的质量流量，kg·s-1；为湿空气的含湿量，kg·kg(干空气)-1；为水膜温度对应的饱和含湿量，kg·kg(干空气)-1。
湿空气的能量守恒方程为：
（3）
式中：为湿空气的比焓，J·kg-1；为0℃时水的汽化潜热，J·kg-1；为湿空气的温度，℃。
把式和式（2）代入式（3）并化简得到：
（4）
式中：为干空气的比热容，J·kg-1·℃-1；为水蒸汽的比热容，J·kg-1·℃-1。
因为，因此式（4）可化简为：
（5）
式（5）表明，空气的干球温度变化主要取决于空气和水膜之间的显热交换，而潜热交换对空气干球温度的变化影响很小。
由于忽略水膜的蒸发损失，所以水膜的能量守恒方程为：
（6）
式中：为水的比热容，J·kg-1·℃-1；为喷淋水的质量流量，kg·s-1。
综上，板式蒸发空冷器传递过程的基本微分方程组由式（1）、（2）、（5）、（6）组成。
2 基于湿球温度的传递模型
根据ASHRAE手册[4]，空气湿球温度在0~200℃范围内，有：

其中：为与对应的饱和蒸汽压，C1=-5.8002206×103，C2=1.3914993，C3=-4.8640239×10-2，C4=4.1764768×10-5，C5=-1.4452093×10-8，C6=6.5459673。

其中，为与对应的饱和含湿量，为大气压力。则湿空气饱和线的斜率b为:
= （7）
引入湿空气的湿球比热容，和湿球对流换热系数[5]，则：

（8）
（9）
其中，,为水的比热容。
在数值计算的步长内假设湿空气饱和含湿量与湿球温度成线性关系，即：，a为待定系数，则有：。
由于，所以可得：
（10）
式（10）重组并微分，把式（2）和（5）代入式（10），整理化简得到：
（11）
把式（11）减去式（5），并把式（8）和（9）代入得到：

（12）
同理把式（8）和式（9）代入式（6）并进一步化简可得到：
（13）
综上，基于空气湿球温度的板式蒸发空冷器传递过程的控制微分方程组由式（1）、（11）、（12）和（13）组成。根据已知条件、边界条件和相应的传热传质准则关联式，联立方程组（1）、（7）、（8）、（9）、（11）、（12）、（13），可求解热流体、喷淋水和空气干湿球温度沿流动方向的变化规律以及换热强度情况。对于逆流的板式蒸发空冷器，喷淋水循环利用，则边界条件为：
，，；，；。
3 传热传质系数的确定
热流体到水膜之间的总传热系数由式(14)确定：
（14）

式中：为热流体侧的对流换热系数，W·m-2·K-1；、和分别为热流体侧的污垢热阻、壁面热阻和水膜侧的污垢热阻，m2·℃·W-1；为水膜的对流换热系数，W·m-2·K-1。
热流体侧的对流换热系数由式(15)确定[6]：
（15）
式中：，为热流体通道当量直径，，为通道间距，m。
喷淋水膜的对流换热由式（16）确定[7]：
（16）
式中： ，为水膜厚度，m；为水的动力粘度，m2·s-1；为
水的密度，kg·m-3；为重力加速度，m2·s1；为水的导热系数，W·m-1·K-1； ，为单位长度喷淋水质量流量，kg·m-1·s-1。
冷却侧空气的对流换热系数计算式由式（17）确定[2]：
（17）
式中：为空气密度，kg·m-3；为通道内空气流速，m·s-1；为空气通道的当量直径，，为空气通道间距，m。根据假设，则空气的传质系数为。
4 实例计算和分析
已知某高炉工艺循环冷却水采用板式蒸发空冷器，其结构参数：长×宽×高为1.2 m×1 m×0.8 m，工艺冷却水通道宽度为2.5 mm, 质量流量为10 kg·s-1，进口温度45℃；空气通道宽度为5 mm进口空气干球温度为35 ℃，湿球温度为28℃，质量流量为4.5 kg,干空气/s；喷淋水采用循环水，质量流量2 kg/s。
图4 各流体沿流程的热流密度分布
图3 各流体温度分布
用四阶-龙格库塔法求解方程组（1）、（7）、（8）、（9）、（11）、（12）、（13），解得各流体沿流动方向的温度变化如图3所示，各流体沿流动方向的换热情况如图4所示。图3与图4结果表明：工艺冷却水温度沿流动方向降低，工艺冷却水放热强度取决于其与喷淋水之间的温差。喷淋水温度是冷却器两侧换热平衡的结果，喷淋水温度呈现先减小后增大的变化趋势，由于其循环利用，所以进出口水温相等；喷淋水温度先降低的原因是由于空气在入口处相对湿度较小，喷淋水传递给空气的显热和潜热总合大于热流体侧传递给喷淋水的热量，喷淋水放热导致喷淋水温降低；随着空气相对湿度增加，喷淋水传递给空气的潜热也在减少，当热流体侧传递给喷淋水的热量大于喷淋水传递给空气的热量时，喷淋水吸热，温度升高。空气湿球温度沿流动方向一致升高，湿球温度上升是由于空气与喷淋水之间的显热和潜热交换导致了空气焓增加；而空气干球温度变化较为复杂，空气干球温度的升降取决于干球温度与喷淋水温度的比较，当喷淋水温高于空气干球温度时，空气吸收显热，干球温度升高，反则降低。
4 结语
(1) 引入空气湿球比热和对流换热系数，建立了基于空气湿球温度的板式蒸发空气冷却器传热传质过程的数学模型微分方程组。
(2）该数学模型的微分方程组形式与空气-空气换热器的基本微分方程组一致，所以该分析模型为进一步分析蒸发式空冷器的热工性能模拟和优化设计提供了理论依据。

参考文献：
[1] R. L. Webb．A unified theoretical treatment for thermal analysis of cooling towers, evaporative condensers and fluid coolers．ASHRAE trans., 1984, vol. 90, Part 2B: 398-415
[2] N. J. Stoitchkov, G I. Dimitrov. Effectiveness of cross-flow plate heat exchanger for indirect evaporative cooling. Int. J. Refrig., 1998, 21(6): 463–471
[3] Guo, X. C.，Zhao, T. S. A parametric study of an indirect evaporative air cooler. International Communications in Heat and Mass Transfer，1998，Vol.25（2）：217-226

[4] ASHRAE Handbook 1989-Fundamental. American Society of Heating, Refrigerating, and Air-Conditioning Engineers, Inc., Atlanta, 1989
[5] J. E. Braun，S. A. Klein，J. W. Mitchell. Effectiveness models for cooling towers and cooling coils. ASHRAE Transactions, 1989, 95(2): 164-174
[6] 杨世铭编. 传热学（第二版）. 北京：高等教育出版社，1992.
[7] Wilke, W. Wa¨rmeu¨bergang am Rieselfilm. VDI-Forschungsheft-490, Du¨sseldorf, 1962, p. 36