带三根绊线圆柱气动力特性研究.docVIP

带三根绊线圆柱气动力特性研究 　　摘要： 在亚临界雷诺数下（Re=4×104），采用大涡模拟方法对带三根绊线圆柱绕流场进行了数值模拟，研究了圆柱表面风压分布、平均气动力和Strouhal数（St）随绊线位置角（β）的变化规律，分析了带绊线圆柱的流场流态与气动性能之间的相关性，探讨了分离泡对圆柱气动性能的影响机制。研究表明：三根绊线会对圆柱的气动性能造成很大影响，其最大平均阻力系数较带单根绊线圆柱与光圆柱分别提升了32%及50%左右，最大平均升力系数则比单绊线圆柱高37%；带三绊线圆柱的St数低于光圆柱，并会随着β的增加单调递减；受绊线的影响，气流往往会在绊线位置提前分离，从而导致带三根绊线圆柱的尾流宽度较光圆柱更宽，并导致圆柱的阻力系数增大；绊线在β=50°附近时，在绊线分离的剪切层会再附到圆柱表面，形成单侧分离泡，这会导致圆柱表面出现局部强负压，并受到平均升力的作用。 　　关?I词： 涡激振动； 气动控制； 流场特性； 圆柱； 大涡模拟 　　中图分类号： V211； TU352.1文献标志码： A文章编号： 1004-4523（2018）02-0291-09 　　DOI：10.16385/j.cnki.issn.1004-4523.2018.02.012 　　1概述 　　在圆柱表面缠绕螺旋线是一种常用的气动控制措施。螺旋线可破坏圆柱旋涡脱落的规律性，因而可抑制或减小圆柱的涡激振动[1]。缠绕螺旋线的气动措施也被用来减少斜拉桥拉索的风雨激振的发生[2-3]。螺旋线的减振效果取决于螺旋线的数量、直径和间距等参数[4]，缠绕螺旋线圆柱也有存在阻力增大的问题[5]。 　　尽管在圆柱表面螺旋线的气动措施已得到广泛应用，然而螺旋线确切的减振机理尚未被澄清[5]。这主要是因为螺旋状表面会造成圆柱绕流呈现复杂的三维特性。为了降低问题的复杂性，研究者常将缠绕螺旋线圆柱的三维绕流问题简化为带绊线圆柱的二维问题进行研究[6-14]。图1为带三根螺旋线圆柱和带三根直绊线圆柱的示意图。研究表明，带单根绊线圆柱的气动力系数、表面风压系数和Strouhal数（St）等对绊线的位置、直径和雷诺数非常敏感[6-10，13]。基于绊线位置，Nebres和Batill[10]，Igarashi[11]等将带单根绊线圆柱的流场结构分为4种不同的流场流态。Fage和Warsap[15]，James和Truong[16]则对带两根对称绊线圆柱的气动特性进行了试验研究。Alam[17]研究了对带两根对称绊线双圆柱的尾流干扰问题。 　　图1缠绕三根螺旋线及带三根直绊线圆柱示意图 　　Fig.1Schematic of the cylinder with three helical wires and three tripping wires 　　此外，Weave[18]，Hara和Shimada[19]等的研究表明：表面缠绕3或4根螺旋线或螺旋状箍条对圆柱涡激振动的抑制效果最好。然而以往文献主要集中在对带单根绊线或两根对称绊线圆柱的研究，尚未见到对带三根绊线圆柱的研究。研究带三根绊线圆柱的气动性能和流场机理有助于理解多螺旋线对圆柱旋涡脱落和涡激振动的控制机理。 　　本文通过大涡模拟（LES）的方法，在亚临界雷诺数下（Re=4×104），对带单根和三根绊线圆柱绕流场进行了数值模拟，研究了圆柱表面风压分布、平均气动力和St数等气动性能随绊线位置的变化规律，从流场角度分析了绊线对圆柱绕流场流态结构和气动性能的作用机理，并探讨了螺旋线抑制圆柱卡门涡脱的机理。 　　2数值方法和计算模型〖*2〗2.1控制方程和亚格子模型在大涡模拟方法中，大尺度涡通过滤波后的Navier-Stokes方程直接求解，而小尺度的涡则采用亚格子尺度模型（SGS）模拟。与雷诺平均法（RANS）相比，大涡模拟方法可更好地模拟流场中的湍流旋涡，经过滤波函数的滤波，可得到大尺度涡的不可压缩Navier-Stokes方程：ixi=0（1） 　　it+ijxj=-1ρxi+μ2ixjxj-τijxj（2）式中i为滤波后速度，ρ为流体密度，为滤波后压力，μ为流体动力黏度，τij为亚格子应力张量。 　　亚格子应力τij采用Smagorinsky-Lilly的亚格子尺度模型，具有以下形式：τij-13τkkδij=-2μtij（3） 　　ij=12ixj+jxi（4）式中μt为亚格子尺度的湍动黏度，且μt=（CSΔ）2||（5） 　　||=2ijij（6） 　　Δ=（ΔxΔyΔz）13（7）式中Δi代表沿坐标轴i方向的网格尺寸；CS为Smagorinsky常数，取0.1。本文数值模拟采用Fluent软件实现。 　　2.2计算模型和计算参数 　　为了验证本文采用的计算方法和计算参数的正确