入射涡与圆柱相互作用的数值模拟

本文用快速涡方法对入射涡与圆柱的相互作用进行了数值模拟,观察到了入射涡在圆柱表面上诱导的二次分离和三次分离现象。二次涡的产生,与入射涡配对,改变了它们的运动轨迹。二次涡是入射涡“回跳”现象的主要原因。本文还对不同入射涡强度及相互位置作了计算,并分析了不同参数对涡运动轨迹的影响。这些现象与涡的无粘圆柱绕流有着本质的差异。     

涡激诱导并列双圆柱碰撞数值模拟研究

圆柱类结构物的涡激振动是工程中较为常见的一种现象,如果圆柱结构物之间的距离较小, 就会产生涡激诱导碰撞现象,而涡激碰撞会比涡激振动对结构物疲劳破坏产生更严重的威胁.采用浸入边界法模拟流体中的动边界问题,避免了传统贴体网格方法在求解流体中存在固体间碰撞问题时出现数值求解不稳定问题,采用有限元方法对圆柱的运动和碰撞进行求解,通过数据回归方法建立了流体流动条件下的润滑模型,对不同间隙比下涡激诱导并列双圆柱振动及碰撞过程进行了数值模拟, 数值结果表明,如果两圆柱产生了碰撞将会有连续的碰撞发生, 碰撞时出现了多阶频率,振动主频率要比无碰撞时大, 两圆柱碰撞时的相对速度比自由来流速度小;当两圆柱相互接近时, 随着涡环分离角度的逐渐倾斜, 横向流体力先逐渐减小,当两圆柱间涡环开始相互影响发生挤压时, 横向流体力开始逐渐增大;当两圆柱开始反弹时, 两圆柱间形成了低压区, 改变了横向流体阻力的方向,使两圆柱又产生了接近运动,如此反复从而产生了碰撞后横向流体力和圆柱速度的振荡现象.      

横向振荡圆柱绕流的格子Boltzmann方法模拟

基于格子Boltzmann方法(LBM)对不可压横向振荡圆柱绕流问题进行了数值研究. 与传统的求解宏观的N-S方程的数值方法不同, LBM求解此类问题不需要采用动网格, 而且不需要对网格进行特殊处理, 从而节约了计算成本. 结果显示, 当振荡频率增加到相应的静止圆柱绕流的自然涡脱落频率附近时, 圆柱后最新形成的集中涡距离柱体越来越近, 直到达到一个极限位置. 随后, 集中涡突然转向圆柱体另一侧脱落. 当振荡频率接近于静止圆柱的自然涡脱落频率时, 发生频率同步的现象. 随着振荡频率远离自然涡脱落频率, 同步现象消失. 在几种次谐振荡和超谐振荡下, 尾流区的涡脱落频率仍为相应的静止圆柱绕流的自然涡脱落频率.      

涡激诱导并列双圆柱碰撞数值模拟研究

圆柱类结构物的涡激振动是工程中较为常见的一种现象,如果圆柱结构物之间的距离较小,就会产生涡激诱导碰撞现象,而涡激碰撞会比涡激振动对结构物疲劳破坏产生更严重的威胁.采用浸入边界法模拟流体中的动边界问题,避免了传统贴体网格方法在求解流体中存在固体间碰撞问题时出现数值求解不稳定问题,采用有限元方法对圆柱的运动和碰撞进行求解,通过数据回归方法建立了流体流动条件下的润滑模型,对不同间隙比下涡激诱导并列双圆柱振动及碰撞过程进行了数值模拟,数值结果表明,如果两圆柱产生了碰撞将会有连续的碰撞发生,碰撞时出现了多阶频率,振动主频率要比无碰撞时大,两圆柱碰撞时的相对速度比自由来流速度小;当两圆柱相互接近时,随着涡环分离角度的逐渐倾斜,横向流体力先逐渐减小,当两圆柱间涡环开始相互影响发生挤压时,横向流体力开始逐渐增大;当两圆柱开始反弹时,两圆柱间形成了低压区,改变了横向流体阻力的方向,使两圆柱又产生了接近运动,如此反复从而产生了碰撞后横向流体力和圆柱速度的振荡现象.     

二次涡与卡门涡形成过程

本文用离散涡位流理论与边界层理论相结合的方法,同时考虑到边界层及后剪切层分离的离散涡模型,研究并计算了圆柱非定常运动过程中二次涡的发生、发展与影响。研究了对称主涡发生运动不稳定性的直接原因和条件,分析卡门涡的形成过程。本文提出,在二次涡的影响下对称主涡从有涡量补充的非自由涡变为没有涡量补充的自由涡,这是对称主涡发生运动不稳定性的条件,在小扰动影响下可诱发成卡门涡。 本文考虑不可压、大雷诺数、层流情况。所得结果与流场显示结果十分相似。     

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

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

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

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

圆柱固结刚性分离盘结构振动风洞实验研究

分离盘作为常用的流动控制装置之一,对于固定圆柱旋涡脱落具有很好的抑制作用。但是,对于结构振动,固结刚性分离盘是否有类似的抑制作用及控制机制尚待深入研究。本文通过风洞实验,利用非接触式激光位移传感器,研究圆柱固结短分离盘后在"从静止开始"和"连续增加风速"这两种不同的初始条件下的结构响应。结果表明,附加分离盘后,与典型的涡激振动(Vortex-induced vibration, VIV)类似,结构的响应虽然仍然限定在一段风速范围内,但是其幅值比单圆柱VIV幅值大得多。当分离盘长度为0.25D(D为圆柱直径)时,结构响应与圆柱VIV类似,分为两个分支。当分离盘长度为0.4D、0.5D时,结构响应存在明显的迟滞现象,同时,响应会在达到最大振动幅值后突然消失。连续增加风速时,振动随约化速度增加而增强。与典型的涡激振动类似,可以发现明显的频率锁定现象,与典型驰振特性不同。通过FFT分析尾流速度频谱,可以发现,附加分离盘后,频谱中出现大量的高阶谐频成分,它们的出现可能与旋涡脱落模式有关。     

纳米尺度圆柱绕流尾迹区流动形式模拟研究

采用非平衡分子动力学模拟方法,对微尺度低{Re}数下的圆柱绕流问题进行了研究,模拟结果表明:当{Re}<12时,圆柱下游形成对称、无分离的定常流;当{Re}>20时,圆柱下游形成周期性交替出现的对称涡;当12相似文献   

空气中平面激波同圆柱的非定常干涉

平面激波在圆柱上的反射与绕射是目前国内外研究激波非线性复杂现象的典型课題。近年来,Takayama(1987)在入射激波马赫数M_i=1.3及2.6的条件下,得到了在圆柱周围的激波绕射与反射的光测实验结果。众所周知,激波在圆柱上的反射与绕射是一非定常现象,为了了解和分析整个非定常过程的发展和变化,应得到圆柱后较远区域的结果。本文利用两次曝光全息激光干涉仪不仅得到了圆柱周围的波系图案,而且也得到了激波脱离圆柱后驻点后约10倍圆柱直径的范围内的流动形态。同时也给出了不同M_i条件下的三波点轨迹和入射激波与反射激波之间夹角ω随距离变化的实验结果。     

VORTICES INCIDENT UPON AN OSCILLATING CYLINDER: FLOW STRUCTURE AND LOADING

Interaction of an incident vortex street with an oscillating cylinder is addressed using high-image-density particle image velocimetry and simultaneous force measurements. This approach reveals that the timing of the incident vortices relative to the cylinder motion controls the large-scale vortex formation in the near-wake, and thereby the phase shift between the loading on the cylinder and its motion. As a consequence, it is possible to change the sign of the fluid-dynamic work done by the fluid on the cylinder. The incident vortices dramatically shorten the formation length of vortices in the near-wake and yield values of lift coefficient up to a factor of five larger than that for an isolated cylinder subjected to controlled oscillations in the absence of incident vortices. These alterations of the wake structure and loading occur in conjunction with globally locked-on patterns of incident and shed vortices with respect to the cylinder oscillation. Different states of global lock-on are attainable for different values of timing of the incident vortices.     

Study of sub-critical Taylor vortex flow between eccentric rotating cylinders by torque measurements and visual observations

This study of the flow between a rotating inner and stationary outer cylinder revealed secondary vortex flows at Taylor numbers well below the value at which Taylor vortices are formed. A distinct sequence of transitional stages has been observed. Torque measurements, visual observations, and varying the eccentricity have been used to study sub-critical vortices with zero axial flow     

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.     

圆柱尾迹涡的三维演化及结构特征

应用无粘涡丝运动学理论和局部诱导近似(LIA)方法,以Lagrange观点跟踪涡丝在背景流场中运动,用数值方法研究了中等Re数(≈10~3)下圆柱分离尾迹中Kármán涡和涡辫区涡丝的三维演化的机制和动力学过程,及其结构特征。背景流场考虑为尾迹时间平均速度流场和Kármán涡街流场。初始展向小扰动为指数形式和谐波形式。结果指出:Kármán涡和涡辫区中的涡丝具有展向不稳定性,形成流向涡量。在尾迹的初期输运过程中,表现出有序的大、小尺度涡结构。并进一步分析了其产生的机理。     

Large eddy simulation of a vortex ring impinging on a three-dimensional circular cylinder

A vortex ring impacting a three-dimensional circular cylinder is studied using large eddy simulation (LES) for a Reynolds number Re = 4 × 10⁴ based on the initial translation speed and diameter of the vortex ring. We have investigated the evolution of vortical structures and identified three typical evolution phases. When the primary vortex closely approaches to the cylinder, a secondary vortex is generated and its segment parts move inward to the primary vortex ring. then two large-scale loop-like vortices are formed to evolve in opposite directions. Thirdly, the two loop-like vortices collide with each other to form complicated small-scale vortical structures. Moreover, a series of hair-pin vortices are generated due to the stretching and deformation of the tertiary vortex. The trajectories of vortical structures and the relevant evolution speeds are analyzed. The total kinetic energy and enstrophy are investigated to reveal their properties relevant to the three evolution phases.     

Unsteady flow around impacting bluff bodies

The flow resulting from the collision without rebound of generic bluff bodies with a wall in a still viscous fluid is investigated both computationally and experimentally. Emphasis is on the case of a circular cylinder impact (two-dimensional geometry), but comparisons with the flow generated by the impact of a sphere (axisymmetric geometry) are included. For normal cylinder impacts, the two counter-rotating vortices forming behind the body during its motion continue their trajectory towards the wall after the collision, leading to the generation of opposite-signed secondary vorticity at the cylinder and wall surfaces. Secondary vortices forming from this vorticity at higher Reynolds numbers exhibit a short-wavelength three-dimensional instability. Comparison with the sphere impact reveals significant differences in the scales of the vortices after the collision, due to the additional vortex stretching acting in the axisymmetric geometry. This leads to a delay in the onset of three-dimensionality and to a different instability mechanism. Oblique cylinder impacts are also considered. For increasing impact angles, the wall effect is gradually reduced on one side of the cylinder, which favours the roll-up of the secondary vorticity and increases the rebound height of the vortex system.     

Numerical analysis of a flowfield produced by a pair of rectilinear vortices approaching a circular cylinder

Two kinds of discrete vortex methods and a direct method to solve the Navier-Stokes equations were carried out for the flow-field induced by a pair of rectilinear vortices approaching a circular cylinder. The major features of the flowfield observed by smoke visualization, i.e., the large scale secondary vortices and the subsequent rebounding trajectories of the primary vortices were well simulated by both numerical methods.     

Resonance in flow past oscillating cylinder under subcritical conditions

Flow around an oscillating cylinder in a subcritical region are numerically studied with a lattice Boltzmann method(LBM). The effects of the Reynolds number,oscillation amplitude and frequency on the vortex wake modes and hydrodynamics forces on the cylinder surface are systematically investigated. Special attention is paid to the phenomenon of resonance induced by the cylinder oscillation. The results demonstrate that vortex shedding can be excited extensively under subcritical conditions, and the response region of vibration frequency broadens with increasing Reynolds number and oscillation amplitude. Two distinct types of vortex shedding regimes are observed. The first type of vortex shedding regime(VSR I) is excited at low frequencies close to the intrinsic frequency of flow, and the second type of vortex shedding regime(VSR II)occurs at high frequencies with the Reynolds number close to the critical value. In the VSR I, a pair of alternately rotating vortices are shed in the wake per oscillation cycle,and lock-in/synchronization occurs, while in the VSR II, two alternately rotating vortices are shed for several oscillation cycles, and the vortex shedding frequency is close to that of a stationary cylinder under the critical condition. The excitation mechanisms of the two types of vortex shedding modes are analyzed separately.