鼓风式机械通风冷却塔空气动力特性数值模拟研究.docVIP

鼓风式机械通风冷却塔空气动力特性数值模拟研究 李红莉、赵顺安 （中国水利水电科学研究院，北京 100038） 摘要：冷却塔的空气动力特性是确保冷却塔高效的设计关键。本文采用Fluent软件，对内陆核电重要厂用水系统鼓风式机械通风冷却塔在不同塔型条件下的空气动力特性进行了数值模拟计算，获得了冷却塔塔型与冷却塔阻力及填料断面风速分布均匀性的关系。计算结果表明，填料安装高度对鼓风式机械通风冷却塔整塔阻力影响不大，但会影响填料断面上风速分布均匀性，填料安装高度越低，风速分布越均匀，填料阻力系数越小影响越明显；出口收缩段的高度越高，整塔阻力越小，风速分布越均匀。 关键词：鼓风式冷却塔；塔型；阻力系数；风速均匀性 1研究背景 核电厂的重要厂用水系统关乎着核电厂的安全，核电厂特别是内陆核电厂的重要厂用水系统一般采用鼓风式机械通风冷却塔，以适应核电对冷却塔的安全性和抗震性能的要求。这种设计在法国核电较为常见，美国大多数内陆核电厂采用每台机组配2台机械通风冷却塔的标准设计[1]。我国目前还没有内陆核电厂，但是随着我国内陆核电建设的快速发展，鼓风式机械通风冷却塔将会越来越多的被用于内陆核电。 内陆核电重要厂用水系统鼓风式机械通风冷却塔不仅在通风方式上有别于常规的抽风式机械通风冷却塔，在塔型结构布置上也有明显差异，目前我国的相关设计规范对此并没有明确的设计原则与方法[2～4]。为了解塔内气流特性并对塔型进行优化，需要通过相关的研究来确定其空气动力特性。实践证明，通过物理模型试验来研究冷却塔空气动力特性是一个十分有效的手段，但是由于鼓风式机械通风冷却塔模型本身的复杂性及塔型的不固定性，使模型试验研究受到局限，而利用流体力学软件Fluent则十分方便有效。 本文利用Fluent软件建立了内陆核电鼓风式机力塔空气动力计算的三维数学模型，研究了不同塔型条件下塔内气流分布及阻力特性，并对塔型进行了优化计算，提高了鼓风式机械通风冷却塔的配风均匀性，增强了冷却塔的整体性能，为确保内陆核电厂用水系统冷却塔安全稳定运行提供了十分有利的条件。 2数学模型及计算方法 2.1 空气流场控制方程 塔内外流场为等温、不可压流动，其控制方程包括连续方程、动量方程，并选用双方程湍流模式对方程进行封闭，各方程等可写为统一形式［4］： （1） 式中：为空气密度，kg/m3；为空气流速，m/s。各控制方程的变量、扩散系数项与源项如下表1。 表1 控制方程中各变量代表参数 控制方程 连续方程 1 0 0 动量方程（流速） ，， 湍能方程 耗散方程 其中生成项，为空气分子粘性系数；为压力；为紊流粘性系数，由动能和紊动耗散率求出：，为经验常数；和分别为和的紊流普朗特数。 2.2 边界条件 空气动力计算整个计算域分为塔内与塔外两个区域。对于塔外区域的边界：底部为固壁无滑移边界条件；四周及顶部采用压力出口边界条件；塔壳采用固壁边界条件。塔内区域：进风口及塔的出口都设置成内部边界；填料区域设置成多孔介质边界条件，并根据实测填料阻力系数设置各向阻力系数；风机采用Fluent风扇边界条件，根据风机性能参数，拟合出风压p与风速v之间的关系曲线。 （2） 2.3 冷却塔阻力系数及风速分布均匀性计算 鼓风式机械通风冷却塔，气流经由风机鼓入塔内，依次经过进风口，穿过雨区、填料等，并经由出口排入到大气中，气流经过各部分的阻力为该区域前后断面的全压差，可表示为阻力系数与该区域气流动压之积： （3） 则气流经过各区域的阻力系数可表示为： （4） 式中：为气流经过某区域前后断面的全压差，Pa；为空气密度，kg/m3；为填料断面平均风速，m/s。 对于鼓风式机械通风冷却塔来说，良好的进风条件和塔内气流流态，特别是填料断面处风速分布状况是保证其充分进行热质交换的重要条件，均匀顺畅的气流能够最大效率的发挥热水与空气的换热潜能，提高冷却塔换热效率。为此定义一个变量来表征风速分布的均匀性，可表示为： （5） 式中：为填料断面风速分布均布系数； 为填料断面各点风速，m/s；n为点的个数。 3计算结果及分析 图1为本文所研究的鼓风式机械通风冷却塔立面布置图，单元塔格平面尺寸为9