Further study on vortex shedding flow by a circular cylinder

In the present paper the mechanism involved in vortex shedding flows is investigated in detail. In the early stage of the unsteady separated flow the interaction between secondary vortices and primary vortices is quite strong. In the later stage of the flow, corresponding to the vortex shedding the recirculating flow region on each side of the aft body goes through such a cycle: growth-contraction-growth, the secondary separation occurs periodically rather than continuously. The reduction of circulation is taken into account in three cases with different decay factors to study its influence on the prediction of main flow characteristics. Results show that to simulate vortex shedding flow it is necessary to include the reduction of circulation to bring the calculated results into good agreement with experiments. An improved discrete vortex model is suggested in which both the secondary separation and the reduction are incorporated. The processes of vortex shedding, the forces prediction and other flow characteristics are given and some discussions are made. Porject is supported by National Natural Science Foundation of China.     

单窄条控制件对横向振荡柱体尾流 2P模式 旋涡脱落的改变

横向强迫振荡柱体尾流控制是柱体涡激振动控制的基础,在海洋、土木等工程中具有重要意义. 横向强迫振荡柱体尾流中存在一种锁频旋涡脱落模式,即在一个振荡周期内柱体上、下侧各脱落旋转方向相反的一对涡,称为2P模式. 本文将相对宽度b/D=0.32的窄条控制件置于横向强迫振荡柱体下游,对振幅比A/D=1.25, 无量纲振频f_e D/V_∞=0.22,雷诺数Re=1 200的2P模式旋涡脱落进行干扰,并通过改变控制件位置,研究旋涡的变化规律. 采用二维大涡模拟和实验验证方法进行研究,在控制件位置范围0.8≤X/D≤3.2, 0.4≤Y/D≤3.2内,得到了2P, 2S, P+S和另外6种新发现的旋涡脱落模式,并对各模式旋涡的形成过程作了详细描述. 在控制件位置平面上给出了各旋涡模式的存在区域,画出了旋涡脱落强度的等值线图,并发现在一个相当大的区域内,旋涡脱落强 度可减小一半以上,尾流变窄. 发现柱体大幅振荡引起的横向剪切流在旋涡生成中起关键作用. 探讨了控制件对横向剪切流的影响,分析了控制件在每种旋涡模式形成中的作用机制.      

方柱的涡脱落频率及端部影响的风洞实验

本文报告了用压力传感器测量方形柱体涡脱落的实验方法和结果。实验在大气边界层风洞进行。方柱的截面为正方形，高宽比为５。压力传感器在测试前经过严格的标定。测量结果表明，随着来源速度沿高度的变化，方柱的涡脱落频率也随高度而变，而且，在柱体端部附近，涡脱频率和Ｓｔ数发生明显的抑制。和二维柱体一样，在一定的Ｒｅ数范围内，三维柱体的Ｒｏ数与Ｒｅ数也存在一种线性关系。     

湍流涡致振动频率特性

用数值模拟方法对固定圆柱湍流涡脱落频率与弹性圆柱湍流涡致振动频率特性进行了研究,湍流计算模型采用标准κ-ε模型,压力泊松方程提法基于非交错网格系统.研究结果表明:固定圆柱湍流绕流涡脱落频率基本不随雷诺数而变,对于同一固有频率弹性圆柱,涡振频率基本不随雷诺数而变;对于某一固定雷诺数流动涡振频率在一定范围内与系统固有频率有关.     

Control of vortex shedding behind circular cylinder for flows at low Reynolds numbers

It has been observed by researchers in the past that vortex shedding behind circular cylinders can be altered, and in some cases suppressed, over a limited range of Reynolds numbers by proper placement of a second, much smaller, ‘control’ cylinder in the near wake of the main cylinder. Results are presented for numerical computations of some such situations. A stabilized finite element method is employed to solve the incompressible Navier–Stokes equations in the primitive variables formulation. At low Reynolds numbers, for certain relative positions of the main and control cylinder, the vortex shedding from the main cylinder is completely suppressed. Excellent agreement is observed between the present computations and experimental findings of other researchers. In an effort to explain the mechanism of control of vortex shedding, the streamwise variation of the pressure coefficient close to the shear layer of the main cylinder is compared for various cases, with and without the control cylinder. In the cases where the vortex shedding is suppressed, it is observed that the control cylinder provides a local favorable pressure gradient in the wake region, thereby stabilizing the shear layer locally. Copyright © 2001 John Wiley & Sons, Ltd.     

Vortex shedding from a circular cylinder with a slit and concave rear surface

Vortex shedding from short circular cylinders with a slit was studied using a flow visualization and amplitude spectrum analysis of a thermoanemometry probe signal. It was found that a circular cylinder with a slit and concave rear surface produces stronger vortices than other bluff cylinders but that these vortices are very vulnerable to the end wall conditions. It was established that two small splitter plates (tails) fixed directly behind the cylinder at the end walls effectively isolate the vortices shed from the cylinder from the end wall boundary layer effects. For this arrangement a perfect regularity of vortex shedding and almost constant Strouhal number were achieved in the Reynolds number test range of about 250 to 43,000.On a leave from Technical University, 60965 Poznan, Piotrowo 3, Poland.     

钝体尾流控制机理及方法研究进展

首先从涡脱落生成理论出发对钝体尾流控制方法进行了分类,并简单介绍了国内尾流控制研究情况. 之后介绍了我们用窄条或小方柱取代小圆柱后,对Strykowsky和Sreenivasan 控制方法的改进及其在高雷诺数下对圆柱和方柱尾流涡脱落的有效抑制情况, 并探讨了控制件钝度对抑制效果的影响.第3部分用实验数据对各个涡脱落生成模型做了分析与检验, 指出控制件方法的机理与改变钝体分离位置、减小钝体背压吸力、改变流动的展向相关性、 防止钝体两侧剪切层相互作用等无关,而与钝体近尾流速度剖面的局部修正及其稳定性的改变有关. 最后简单介绍了控制件方法今后研究工作展望及其工程应用前景.     

流向振荡柱体尾流控制研究进展

详细介绍了近几年采用尾部喷射、隔离板和小窄条控制件等3 种方法对流向振荡柱体绕流旋涡脱落的抑制情况. 在研究范围内存在非锁频和3种锁频旋涡脱落模式. 风洞实验表明, 尾部喷射对这4 种模式都有抑制效果,窄条控制件对非锁频和2种锁频模式具有抑制效果, 而隔离板仅对非锁频和1 种锁频模式有效. 在不同流动和振荡条件下找出了每种方法的有效控制区, 研究了减阻和减少脉动力的效果, 并探讨了控制机理.      

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

圆柱绕流涡脱落诱发较大的振动和声,如何有效地抑制值得关注.利用大涡模拟技术求解了Navier-Stokes方程,得到了涡脱落频率,升力脉动幅值及平均阻力系数.计算表明二维模拟不能体现流动基本特征,三维计算与实验吻合较好.为了抑制涡脱落,在直径为D的圆柱表面装入间距为1D,直径为0.0167D的O型环.通过升力、速度谱分析以及柱向横截面流场分析可知,在光滑圆柱外表面加入O型环能诱发流体边界层分离,有效地抑制涡脱落现象,升力脉动和观测点速度脉动幅值几乎完全消失,阻力系数也略微降低,适合在实际工程中采用.     

隔离板对流向振荡圆柱尾流旋涡脱落的抑制

用隔离板对直径为D, 沿流向振荡的圆柱后涡脱落进行抑制. 隔离板放于圆柱尾流中心线上,控制参数包括隔离板长度L/D以及隔离板前缘到柱体振荡中心的距离G/D. 实验的雷诺数范围Re=V_∞D/v=1.01×10⁴～1.69×10⁴,柱体折减振频范围f_eD/V_∞=0～0.03, 柱体振幅固定为A/D=0.2. 风洞烟线显示和热线测量结果表明:当 G/D位于一个有效区域内时,可有效抑制振荡柱体尾流的旋涡脱落. 该有效区的大小随着隔离板板长的增大而增大, 随着Re数和圆柱振荡频率的增大而减小.     

外部激励Stuart-Landau模型的数值研究

将非线性常微分方程组周期解的求解看作一个边值问题 ,运用Newton迭代构造求解这组方程的数值方法。利用上述方法求得了激励Stuart Landau方程的周期解 ,研究了圆柱振动对圆柱后Karman涡街的抑制现象 ,和振动的频率锁定现象 ,证明了激励Stuart Landau方程描写钝体尾迹动力系统的有效性     

基于浸入边界-格子Boltzmann通量求解法的椭圆柱流动特性分析

基于浸入边界-格子Boltzmann通量求解法,开展了雷诺数Re=100不同几何参数下单椭圆柱及串列双椭圆柱绕流流场与受力特性对比研究。结果表明,随长短轴比值的增加,单椭圆柱绕流阻力系数先减小后缓慢上升,最大升力系数则随长短轴比值的增大而减小;尾迹流动状态从周期性脱落涡到稳定对称涡。间距是影响串列圆柱及椭圆柱流场流动状态的主要因素,间距较小时,串列圆柱绕流呈周期性脱落涡状态,而椭圆柱则为稳定流动;随着间距增加,上下游圆柱及椭圆柱尾迹均出现卡门涡街现象,且串列椭圆柱临界间距大于串列圆柱。串列椭圆柱阻力的变化规律与圆柱的基本相同,上游平均阻力大于下游阻力;上游椭圆柱阻力随着间距的变大先减小,下游随间距的变大而增加,当间距达到临界间距时上下游阻力跃升,随后出现小幅度波动再逐渐增加,并趋近于相同长短轴比值下单柱体绕流的阻力。     

Flow past a sphere in density-stratified fluid

Stratified flow past a three-dimensional obstacle such as a sphere has been a long-lasting subject of geophysical, environmental and engineering fluid dynamics. In order to investigate the effect of the stratification on the near wake, in particular, the unsteady vortex formation behind a sphere, numerical simulations of stratified flows past a sphere are conducted. The time-dependent Navier–Stokes equations are solved using a three-dimensional finite element method and a modified explicit time integration scheme. Laminar flow regime is considered, and linear stratification of density is assumed under Boussinesq approximation. The effects of stratification is implemented by density transport without diffusion. The computed results include the characteristics of the near wake as well as the effects of stratification on the separation angle. Under increased stratification, the separation on the sphere is suppressed and the wake structure behind the sphere becomes planar, resembling that behind a vertical cylinder. With further increase in stratification, the wake becomes unsteady, and consists of planar vortex shedding similar to von Karman vortex streets.     

Experimental study on the wake behind tapered circular cylinders

The flow around tapered cylinders can act as basic models for numerous bluff body flows with a spanwise variation of either the body shape or the inflow conditions. The well-known vortex street is influenced by strong three-dimensional effects from the spanwise variation of the shedding frequency, namely oblique vortex shedding and vortex dislocations. Stereo-PIV was chosen to study these phenomena, since it allows analyzing planes with the full three-component, instantaneous velocity fields and local, time-dependent variations in the same setting. Hence, detailed aspects of the vortex dislocation phenomenon are presented. Single vortex dislocation events are presented through the local variation of the three measured velocity components u, v and w. Longer time-series reveal both period and location of these dislocation events, as well as quantity and sizes of the cells of constant shedding velocity in between them. The influence of the Reynolds number and the cylinder aspect ratio on the vortex cells could be shown. The analysis of the vortex shedding behavior shows good agreement with previously published results. At the same time, the applied PIV technique provides more spatial information than point-based measurements and offers insight into a Reynolds number range that is currently out of reach of Direct Numerical Simulations.     

流向强迫振荡圆柱绕流的涡脱落模态分析

通过求解采用ALE方法描述的运动坐标系Navier-Stokes方程组,分析均匀来流下雷诺 数为150的静止和流向振荡的圆柱绕流. 主要研究了强迫振荡频率和较大振幅比 (A/D=0.3-1.2)对圆柱升力、阻力变化特性以及涡脱落模态的影响. 研究表 明,流向振荡圆柱绕流存在多种涡脱落模态,如对称S以及反对称A-I, A-III, A-IV等多种形式;比较研究结果,拓展了各模态下对应的锁定区域,并将其分为5个 子区;A-I模态中圆柱受力较以前所知更复杂;通过分析计算结果,发现最大加速度 比Af_{c}^{2}/Df_{s0}^{2}可能是涡脱落模态(尤其是对称S模态)最有效的控制参数.     

串列双圆柱绕流问题的数值模拟

本文运用有限体积方法,对绕串列放置的双圆柱的二维不可压缩流动进行了数值计算。为研究两圆柱不同间距对圆柱相互作用和尾流特征的影响,选取间距比Ｌ／Ｄ（Ｌ为两圆柱中心间的距离,Ｄ为圆柱直径）在１．５～５．０之间每隔０．５共八个有代表性的间距进行了计算模拟。计算均在Ｒｅ＝２００条件下进行。计算结果表明：对该绕流问题,流动特征在很大程度上取决于间距的大小。且间距存在一临界值,间距比从小于临界值变化到大于临界     

均匀来流中横向振动圆柱近迹涡结构的数值模拟

从涡量流函数形式的N-S方程出发,在不同的振动频率、振幅及Re数下数值研究了均匀来流中横向振动圆柱粘性统流的涡脱泻现象。着重探讨了近迹复杂的涡结构及其非定常演化过程,以及它们对物体受力特性的影响,并首次成功地模拟了近年来实验研究中所发现的一些重要的流动现象,如相位“开关”现象(phase switch phenomena)及复杂涡结构(complex vortex structure),等等。通过数值模拟,不仅能够再现实验研究中所发现的一些重要的流动现象,还可进一步预示某些新的流动现象,使数值计算起到与实验研究相辅相成的作用。     

Induced parallel vortex shedding from a circular cylinder at Re of O(10) by using the cylinder end suction technique

In the present work, the objective is to attempt to induce parallel vortex shedding at a moderately high Reynolds number (=1.578 × 10⁴) by using the cylinder end suction method, and measure the associated aerodynamic parameters.We first measured the aerodynamic parameters of a single circular cylinder without end suction, and showed that the quantities measured are in good agreement with equivalent data in the published literature. Next, by using different amount of end suction which resulted in increasing the cylinder end velocity by 1%, 2% and 2.5%, we were able to show that the above corresponded to the situation of under suction, optimal suction and over suction, respectively. With optimal suction, we demonstrated that the end suction method works at Re = 1.578 × 10⁴. The shape of the primary vortex shed became straighter than when there is no end suction, and parameters like cylinder surface pressure distribution, drag force per unit span, as well as vortex shedding frequency all showed negligible spanwise variation. Further careful analyses showed that when compared to the naturally existing curved vortex shedding, with parallel vortex shedding the mid-span drag per unit span became slightly smaller, but the drag averaged over the cylinder span became slightly larger. For cylinder surface pressure, it was found that cylinder end effects mainly influenced the surface pressure in the angular ranges −180°  β < −60° and 60° < β  180°. Without end suction, the cylinder surface pressure in the above ranges was found to increase (become less negative) slightly with |z/d|, but such increase disappeared when optimal end suction was applied, and the cylinder surface pressure distribution became spanwise location independent. As for the vortex shedding frequency (Strouhal number), although the Strouhal number showed spanwise variation when there is no end suction and negligible spanwise variation when optimal suction was applied, the difference between the spanwise averaged Strouhal number was quite negligible. With under suction, the spanwise dependence of various aerodynamic parameters existed, but was found to be not as significant as when no end suction was applied at all. With over suction, the flow situation was found to be practically no change from the optimal suction situation.     

A numerical exploration of secondary separation in starting flow around a flat plate by the discrete vortex model

In the present work, a further numerical simulation of the starting flow around a flat plate normal to the direction of motion in a uniform fluid has been made by means of the discrete vortex method. The secondary separation occurring at rear surface of the plate is explored, and predicted approximately using Thwait's method. The calculated results show that in the early stages of the flow secondary separation does occur. The evolution of flow field, the vortex growing process and the characteristics of secondary vortices have been described. The time dependent drag coefficients, the vorticity shed from the edges and rear surface, and the separation positions are calculated as well as the distributions of velocity and pressure on the plate. In the case of flow normal to the plate, the calculated secondary vortices are weak. Their existence will change the local velocity distributions and affect pressure distributions. However, the effect on drag coefficient is negligible.