立交盘静态混合器流动状态的计算流体力学研究.docxVIP

立交盘静态混合器流动状态的计算流体力学研究 静态混合器可以通过不同的流管结构在非常宽的雷诺数范围内混合水流，而不需要机械可动件的流管电路结构。两个或两个以上的混合机制是将两个或两个以上的流通过这些因素反复提取，形成完整的混合。自静态混合器的引入以来，国内外已经应用了几种静态混合器，其中最广泛的是美国的凯斯-楠木、瑞士的苏尔士smx和smv，以及日本的东利高效风格。 随着计算机技术的发展和计算流体力学方法 (computational fluid dynamics, CFD) 的日趋成熟, 越来越多的学者采用CFD方法对静态混合器的流体力学行为进行研究.无论是Kenics型还是SMX和SMV型静态混合器, 国内外学者都有相关报道.大多数报道的结果显示静态混合器的流场由雷诺数、混合元件的结构、摆放形式及个数等条件决定.Zalc 等在中低雷诺数下, 计算了SMX型静态混合器中的流动特性.认为在Re小于1时, 流动情况不随雷诺数变化;而由于惯性力的存在, 在较高的Re下, 情况完全不同.Hobbs等对Kenics静态混合器进行模拟, 将研究重点放在对一个混合元件内部流场的变化上.Song等主要考察了Kenics静态混合器的压降情况, 并得出摩擦因子与雷诺数和混合元件长径比的关联式.还有一些学者的数值模拟针对其他形式的静态混合器, 例如Mickaily-Huber等、Wadley等、Heniche等分别对SMRX、SMV、KMV型静态混合器进行了研究.研究普遍将对静态混合器的流体力学性能的充分认识作为传热、混合等其他方面研究的起点.立交盘静态混合器是一种新型的静态混合器, 对这种设备的研究还仅局限在工程应用所关心的一些实验研究上, 缺乏细致的理论研究, 尤其是基础的流场分布的研究.目前, 只有少量报道对立交盘静态混合器的压降和混合效果进行实验测定, 但是无法得知该结构对流体流动的影响以及内部流动的具体状态.笔者采用CFD方法对牛顿流体在立交盘中的流动状况进行模拟, 并着重混合元件内部流场的细节变化, 发现结构变化对流场的影响;并进一步计算了压强变化情况, 对结构改进和优化提供了指导依据. 1 态混合器的结构 立交盘静态混合器的结构及其流道分布见图1, 盘内设有从中央区到周边区和从周边区到中央区的通道, 并相互交叉, 结构上与立交桥相似, 因此而得名.这种静态混合器由3种不同形式的叶片组成:贴着管路壁面的折弯平面组成的外叶片, 位于中央的扇形内叶片以及起连接作用的四棱形叶片.立交盘外轮廓与圆管内壁吻合, 内轮廓线将流经的流体分为两大部分:中央区流体和周边区流体.中央区流体进入到立交盘后受到其中两片内叶片阻隔而分成两股, 这两股流体在前行中会遇到另外两个内叶片, 再次受阻后从周边区流出;周边区流体流经立交盘的过程中, 受外叶片剥离和引导逐渐进入到中央区流出.这样的两个通道就实现了对流体的分割以及交叉换位. 2 圆管及压降计算 立交盘静态混合器的三维稳态流场和压强由CFD软件FLUENT模拟, 采用非结构化网格结构划分装有静态混合器的圆管计算域, 总共包括544 463个四面体和103 013个节点.对速度场的研究集中在两个立交盘的管路上 (图2) .计算中所用到的流动介质和物理模型的参数见表1.而对于压降的计算, 为了减小计算误差, 选取了尽量长的计算区域, 总共计算了10个立交盘的压降情况, 并且盘与盘之间紧密排列.所有的湍流计算均选用标准κ-ε模型. 3 结果与讨论 3.1 轴向速度对比 为详细研究立交盘结构所引起的流场变化, 在两个盘上选取了6个测试截面, 见图2:截面1为未进入立交盘, 距入口20 mm处, 以立交盘元件进口为起始位置, 截面1、2、3、4、5和6在轴向距离上分别为-2L/5、L/5、L/2、4L/5、7L/5和3L.得到了6个截面在Re为87时的三维速度场. 对于轴向速度 (图3) , 截面1为典型的层流速度分布.对于截面2, 由于立交盘混合器的叶片对流体的拉伸作用, 可明显看出圆心周围的高速区被拉长变形;边界层内的低速流体受混合器叶片剥离并逐渐进入中央区, 其流速大小在这一过程中得到提升.随着轴向速度的增加, 混合器流场变化更加复杂, 截面3处, 叶片的十字交叉结构将整个圆形流道分成大小相等的4个部分, 各部分内流速分布类似.在截面4的分布图中, 截面3图中的上下两个区逐渐靠拢并联成一体, 左右两个区也被拉长了;同样, 接触到混合器叶片的边界层被剥离, 到此, 进入混合器的全部边界层更新完毕, 完成一次周边区和中央区流体的换位.截面5的流速分布显示, 高速流体由中心向外分散.第二个元件的出口截面6和第一个元件的出口截面5的速度分布很相似, 体现了大多数静态混合器所具有的“周期相似”特性.从轴向速度的大小看, 不