基于CFD轴流通风机三维建模及其内部流场的探讨.docVIP

— —2009年第2期 摘要：运用FLUENT 软件对轴流通风机进行了三维实体建模，并对其内部三维粘性流场进行了数值模拟，得到了内部流动的一些规律，并进行了性能分析和预测，为风机的内部结构优化提供了依据和方向。 关键词：轴流式通风机；数值模拟；性能中图分类号：TH432.1 文献标识码：B 文章编号：1006-8155（2009）02-0016-04Discussion of 3-D Modeling for Axial -flow Fan and Its Internal Flow Field Based on CFD Abstract ：The 3-D solid modeling for axial-flow fan and the numerical simulation for its internal 3-D viscous flow field are carried out by FLUENT software. Some regulars of internal flow are obtained and performance analysis and prediction are carried out. The results have offered important foundation and direction for optimizing internal structure of fan. Key words ：axial-flow fan ；numerical simulation ；performance 0引言 风机内部流场的流动状况直接决定风机的 性能[1] ， 因此，弄清风机内部流场对改进和开发高性能的风机产品有着非常重要的意义。传统的风机设计方法已不能满足高效、节能和降低噪声的要求，本文根据实际轴流通风机的结构参数，利用CFD 进行计算模拟和优化轴流通风机的内部流场，从而为轴流通风机的内部结构优化设计提供了依据。 1 数学模型的建立 1.1 实体建模和网格划分 对流体机械内部流场进行数值模拟的第一 步是在FLUENT 前处理软件GAM BIT 中完成几何建模和生成计算网格。 作为三维实体模型，风机的叶型是通过多个截面的参数来定义，在建模过程中，采用由点到线再到面再到体的建模方法。将每个截面的几何尺寸转换成三维坐标后，通过GAM BIT 中的Turbo 模块来生成风机实体。首先要编写一个后缀为tur 的文件，定义叶片的一些关键参数：如旋转轴、叶片顶端和外壳的间隙，叶片数目，截面数目及各个截面的详细坐标等。编写好tur 文件后，将其导入到GAM BIT 中，就可以得到叶片各截面的轮廓图，见图1。用GAM BIT 工具中的Turbo 命令组定义轮毂和机壳的进口以及各个截面的进口，就可得到图2所示的各个截面同心圆示意图。适当调整进口和出口的中心线位置，定义叶片的数目，就可以生成一个透平体单位，见图3。生成透平体图后，定义轮毂、机壳、进口、出口、压力面、吸力面等关键表面；定义叶片根部截面和顶端截面的网格间距；然后就可以生成透平体的网格，见图4。1.2 计算模型的选取 本风机模型包含旋转的动边界和静止不动的静边界，因此将计算区域划分成定子和转子两个子区域。文中旋转叶轮和静止机壳之间的耦合采用了M RF 模型。生成约70万四面体非 刘 松 吴树福 张立祥/安徽理工大学 基于CFD 轴流通风机 三维建模及其内部流场的探讨 收稿日期：2008－11－26 淮南市232001 16 — —结构网格，经检查，网格倾斜度均在0.85以下，符合计算要求。 1.3边界条件和相关设置 操作条件：设工作压力为一个大气压，即所 求的压力值均为相对于大气压的表压力。 气流进口：采用质量入流边界条件（M ass-Flow Inlet ），作为整个流体区域气流的进 口，入流参数随不同的工况变化。 气流出口：取压力出口边界条件（Pressure Outlet ），作为整个流体区域气流的出口。定义出口压力相对大气压力为0， 即没有附加的压力作用。固壁边界：风道内壁、龙骨、支撑及叶片均取壁面边界条件（wall ）。在具体的设置中选择无滑移固壁边界条件。 将叶片的吸力面（Suction ）和压力面（Pressure ）以轮毂定义为旋转壁面条件，以便在后面的计算中得到旋转的力矩，以计算风机的全效率。 在旋转壁面条件的定义中，按照M RF 的要求，将旋转壁面的旋转速度定义为相对速度（Relative Velocity ），并且相对周围流体速度为0。1.4 数值计算方法 采用SIM PLE 算法和Realizable k -ε湍流模型，非耦合隐式求解器，求解三维时均雷诺N-S 方程[2]，具体的理论方程不在赘述。 2 计算结果