利用O型环抑制圆柱绕法流涡脱落的数值研究.pdfVIP

第Z5卷第4期 计算力学学报 V01．25。No．4 Chinese of Mechanics 2008 2008年8月 Journal August Computational 利用O型环抑制圆柱绕流涡脱落的数值研究 黄翌宇’， 蒋伟康 (上海交通大学振动、冲击、噪声研究所机械系统与振动国家重点实验室，上海200240) 擒要：圆柱绕流涡脱落诱发较大的振动和声，如何有效地抑制值得关注。利用大涡模拟技术求解了Navier- Stokes方程，得到了涡脱落频率，升力脉动幅值及平均阻力系数。计算表明二维模拟不能体现流动基本特征。三 维计算与实验吻合较好。为了抑制涡脱落．在直径为D的圆柱表面装入间距为1D，直径为0．0167D的0型环。 通过升力、速度谱分析以及柱向横截面流场分析可知，在光滑圆柱外表面加入。型环能诱发流体边界层分离， 有效地抑制涡脱落现象，升力脉动和观测点速度脉动幅值几乎完全消失，阻力系数也略微降低，适合在实际工程 中采用。 关键词：圆柱绕流，大涡模拟；O型环；涡脱落抑制 中图分类号：0351．2 文献标识码：A 1 引 言 2数值模拟方法 圆柱绕流涡脱落是工程中广泛存在的一种现 2．1涡脱落频率表达 象，如空调通风系统的管道、换热器，因此对它的研 在圆柱绕流问题中，涡脱落会引起升力和阻力 究引起了广大学者的关注[1刁]。但就目前而言，圆 脉动从而诱发圆柱振动或诱发气动声。旋涡脱落 频率可以表达为 柱绕流的涡脱落抑制机理尚未完全掌握。数值计 f—s，U／D (1) 算方面，雷诺平均模型通常不能对复杂、高雷诺数 式中S，是相应的Strouhal数，【厂为自由来流速度， 流动进行有效地求解，而直接湍流模拟又太耗时， D为圆柱直径。 因此大涡模拟的应用势在必行[4]。通常涡脱落都 升力和阻力系数的计算表达式为 采用被动控制，给圆柱表面加突起或凹槽，圆柱后 ————’————～‘，O， 加分离平板，使流动呈流线型等都是较常用的方 法[5’6]。而对于圆柱表面加O型环控制涡脱落的式中FL和FD分别为作用于圆柱上总的升力和阻 力，H为圆柱体的高度。 方法非常简单适用，但仅在Lim[73的试验中提及 2．2模型的建立 到，并且只分析了阻力的减少效应。由于0型环 带有O型环的圆柱体如图1所示，O型环间距 与圆柱直径相比非常小，在周围需要细致的网格才 X10DX 为1D，直径为0．0167D。计算区域为18D 能得到较为满意的结果，这就大大增加了数值计算 2D，圆柱附近网格俯视图如图2所示。光滑圆柱共 的难度。本文通过大涡模拟方法求解了三维 计网格85440