一种V型百叶窗内部流场的数值模拟.pdf

一种V型百叶窗内部流场的数值模拟 湖南大学 马 扬 龚光彩 广东省建筑科学研究院 杨仕超 摘 要 采用k-e紊流模型对一种V型百叶窗过滤器内部流场进行了数值模 拟，采用控制容积积分法离散控制方程，对流项为二阶迎风格式，应用SIMPLE 算法进行求解。得出了几种不同过滤风量下的流场矢量图，了解了该种过滤器内部 流场的分布规律。其中采用的一种处理格栅障碍物的新方法具有一定的通用性，适 合用来处理空调风口等。获得的模拟程序具有实用价值，只需稍加改动部分子程序 即可用来模拟空调器风口，空调系统中散流器的出风等情况，还可用于该过滤器内 部的两相流的数值模拟。 关键词 百叶窗过滤器数值模拟k-e紊流模型 计算流体力学CFD是指运用先进的软件和计算技术模拟与分析流场结构、传质与传热 及相关的现象。CFD现已成为工程设计与环境分析工程师的一个有力的工具。它能在工程 设计实施前或产品生产前预测方案的实际效果与产品的性能【1]。百叶窗式过滤器广泛应用于 众多领域，可是却没有对它的内部流场进行过数值模拟。本文研究的过滤器是一种由12组 V型挡板相互耦合而成的，在它内部有很多弯道，当含尘气流通过时，在通道内方向多次变 化，使尘粒产生离心力脱离气流流线而撞到板壁上，并在挡板上由上而下运动至贮灰槽而分 离。它的局部具体结构见下图1。本文对该过滤器内部流场进行了数值模拟。并形成了一套 具有一定实用价值的数值模拟程序。 1物理模型和数值计算方法 1．1计算区域的简化、网格生成及障碍物的处理 实施过滤器内部流动数值模拟的困难主要有 两方面：一是百叶窗的挡板排列紧密，必须生成十 分细密的网格才能实现有效的空间离散，对计算 八八 机的内存容量要求高；二是挡板障碍物的处理。 由于本文研究的过滤器是由12组相同的挡板耦合 ／V＼ 的，通过每个弯道的气流具有相同的规律，所以 本文只模拟两组挡板，即图1中左半部，这样可 图1过滤器截面结构图 以减小运算量。 障碍物的处理和网格划分是联系在一起的。由于矩形障碍物的处理方法已经很成熟，就 是要使障碍物处的内部各个内点的速度以及其他各量保持一定值Ⅲ，本文就通过微观改变挡 板障碍物的形状成为矩形，以达到可利用原来的处理方法。具体方法如图2所示：挡板与z 489 =0．00143 m，这样可以用一个个小的矩形体代替原本的挡板。而且，由于实际的挡板表 面并不是完全的光滑的，用这样的小矩形体代替更符合实际情况，并不会影响计算结果 的真实可靠性。由于速度沿z轴变化不大，在Z轴的网格可以适当稀疏一些，取 △z=O．005ITI。为了更好的展现气流经过过滤器时的流场变化，计算区域在进、出风口要 适当的扩大一些，本文在进、出风口的计算区域各延长了0．02In。 图2删格划分 1．2流场模拟和边界条件 在本文研究中，气相流动模拟采用肛e紊流模型‘3|。方程组如下： 流体连续方程：鑫(融)=o (1) (2) 流体i方向动量方程：矗(鹏¨一差+矗[p。(罄+罄)]+Apgi 流体湍能方程：毒(腭走)=毒(戋差)+G—IDe (3) (4) 流体湍能耗散率方程：蠹(融e)=矗(尝蠢)+詈(c-G—czPe) 流体能量方程：矗(丹^)=鑫[等差]一qr (5)