轴对称湿式冷却塔水汽相变换热机理数值模拟研究.pdf

轴对称湿式冷却塔水汽相变换热机理的数值 模拟研究 詹杰民 ，余凌晖，徐士倩,岑显荣1 1 2 1 1 （中山大学应用力学与工程系，510275） 2 ( 山东电力工程咨询院 济南 250013 ） 摘要：本文通过FLUENT软件，采用二维轴对称模型和多相流的欧拉模型对冷却塔水汽相变换热的散热机理进行 研究。针对被广泛接受的设计规范中的案例，本文采用相对密的网格进行计算和可靠性的验证，把填料层的流体区 域设为多孔介质条件，以模拟空气通过填料层时的压力损失。根据Merkel方程，水和空气的热量交换分接触传热和 相变传热两部分，接触传热部分使用FLUENT里现有的模型计算，相变传热是分别对空气加上质量源项（不加能量源 项是因为FLUENT能量方程里空气的焓不包括水蒸汽的信息，水蒸汽的部分另外用自定义标量（UDS）表示），而对水 则加上质量源项和能量源项实现的。湿空气的密度根据理想气体定律采用用户自定义函数（UDF）改写，并考虑重力 的影响。通过添加求解自定义标量（UDS）方程，把环境湿度的影响加进计算模型里。根据不同空气干球温度下的饱 和空气含湿量的不同，对不同条件下的环境相对湿度进行模拟。本文系统地研究了不同环境湿度在冷却塔降温效果 的方面影响。 关键词：轴对称模型，欧拉模型，Merkel焓差法,UDF,UDS 引言 随着人们对冷却塔研究的不断深入，关于冷却塔相变换热机理的研究正逐渐受到学者们的关注， 其研究手段也从传统的模型试验发展到利用先进的 CFD 技术对冷却塔的空气动力场和温度场进行数 值模拟。 在国外，Gan等人[1-3]基于Merkel方程建立了一系列的冷却塔热工数值模型，这些模型都是在一 维稳态的假设条件下对雨区水滴及填料层液膜两种特殊的流动传热过程进行了简化，但它们在研究 [4]等人使用PHOENICS软件建立了三维稳态的 自然风对冷却塔的影响方面却显得无能为力。Spalding 冷却塔CFD热工模型，发现即使在均匀的来流条件下，空气流动在塔内的分布也是不均匀的。Rafat Al-Waked和Masud Behnia[5]用数值方法研究了在不同环境和侧风条件下自然通风蒸汽冷却塔(NDWCT) 内传热和传质过程。他们采用了标准的三维k–ε紊流模型。在进行数值模拟时，对空气采用了欧拉 法，而对于液体的水则采用了拉格朗日法。通过计算塔内各部分的换热量发现，风速增大时雨区的 传热和传质会增强，因此，雨区的存在能够有效地降低横向风的不利影响。在国内，周兰欣[6]等利 用FLUENT软件和Poppe模型，研究了淋水密度、进塔水温和环境条件等因素对冷却塔热力特性的影响。 王凯、孙奉仲[7]等通过实验与数值计算相结合的方法，研究了自然风对冷却塔的影响规律。本文通 过FLUENT软件，采用二维轴对称模型和多相流的欧拉模型对冷却塔水汽相变换热的散热机理进行研 究。 1．模型与参数 下面是冷却塔运行系统示意图（图 1），热水从淋水系统洒下，与冷空气接触进行质量和热量的 交换，经过填料层和雨区，冷却后的水最后从集水池中流出。 2．计算原理 2.1 基本控制方程： 连续性方程 在FLUENT里我们采用Eulerian多相流模型，空气 作为连续相（主相），而水则作为次相。在Eulerian 模型中，q相的连续性方程为[8] ∂ v n (1) (α ρ ) +∇⋅(α ρ v ) ∑(m =−m ) +S q q q q q pq qp q ∂t p 1 其中， p 相为除 q 相以外的其他相，αq 是q 相的体 v q m 是p 相到q 相的质量 积分数；