层流条件下并列三圆柱涡激振动响应与尾流形态

应用基于嵌入式压强-力迭代的高精度浸入边界法研究等间距并列三圆柱涡激振动。其中,雷诺数Re=100,间距比T/D=2.0～5.0,圆柱质量比m*=2.0,折合流速Ur=2.0～10.0,忽略振动系统的阻尼且三圆柱仅横向振动。研究发现,圆柱的振动响应随折合流速的增加呈现初始响应分支和下端响应分支两种模式;振幅响应出现不连续现象,且随着间距比的增加,该不连续现象对应的折合流速增加;尾流模式与间距比和折合流速密切相关。共发现六种尾流形态,分别为窄宽窄尾流、不规律尾流、反相同步尾流、调制尾流、同相同步尾流和偏斜尾流。总结并绘制了尾流形态在参数空间[Ur,T/D]内的分区图。     

正弦振荡来流下柔性立管涡激振动发展过程

在风浪流的作用下,海洋浮式结构物将带动悬链线立管在水中作周期性往复运动,从而在立管运动方向上产生相对振荡来流,这种振荡来流将激励立管悬垂段发生“间歇性” 的涡激振动. 在海洋工程水池中对一个4m 长的立管微段进行模型试验研究,以探索相对振荡来流作用下立管涡激振动产生的机理及其发展的物理过程. 试验通过振荡装置带动模型作正弦运动来模拟不同最大约化速度U_Rmax、不同KC（Keulegan-Carpenternumber）的相对振荡来流,利用光纤应变片测量立管涡激振动响应. 结合模态分析方法处理试验数据得到位移响应时历,继而提出相对振荡来流下柔性立管涡激振动发展的3 个阶段:建立阶段、锁定阶段以及衰减阶段. 并进一步总结了最大约化速度U_Rmax,KC 对涡激振动发展过程的影响规律. 最终获得不同最大约化速度U_Rmax下,涡激振动各发展阶段随KC 所占时间分布比例图.      

串列双圆柱流致振动的数值模拟及其耦合机制

对雷诺数Re= 100 条件下串列双圆柱的流致振动进行了数值模拟. 圆柱的质量比m*均为2.0,间距比L/D 为2.0 5.0. 考虑两种工况:(a) 上游圆柱固定,下游圆柱可沿横流向自由振动;(b) 上、下游圆柱均可沿横流向自由振动. 结果表明:无论上游圆柱静止或者振动,下游圆柱横向振幅明显大于单圆柱的. 工况(b) 的下游圆柱最大振幅要大于工况(a) 的,这是由于两圆柱均振动时,两圆柱之间耦合作用增强,上游圆柱的尾流和下游圆柱的振动之间“相互调节” 作用显著. 对工况(b) 的下游圆柱振动和间隙流之间的作用机制进行了详细的研究,发现当上游圆柱脱落的自由剪切层重新附着于下游圆柱上并且完全从间隙之间通过时,下游圆柱的振幅最大.      

低雷诺数下串列三圆柱涡激振动中的弛振现象及其影响因素

对间距比为1.2和雷诺数为100的串列三圆柱涡激振动进行数值模拟, 发现在某个折合流速之后, 三圆柱的响应均呈现为随着折合流速增大而增大的弛振现象, 平衡位置偏移、低频振动以及旋涡脱落与圆柱运动之间的时机三个因素共同决定了弛振现象的出现. 进一步的研究发现, 串列三圆柱的弛振现象仅出现在质量比不大于2.0和雷诺数不大于100的工况下. 当质量比较大时, 串列三圆柱的平衡位置固定不变, 且圆柱的振动不规律, 使得旋涡脱落与圆柱运动的时机处于变化之中. 当雷诺数较高时, 最上游圆柱的平衡位置在折合流速较大时回到初始位置, 不再参与对圆柱振动的调节, 使得圆柱的振动响应不再规律, 旋涡脱落与圆柱运动的时机也一直处于变化之中.      

低雷诺数下串列三圆柱涡激振动中的弛振现象及其影响因素

对间距比为1.2和雷诺数为100的串列三圆柱涡激振动进行数值模拟,发现在某个折合流速之后,三圆柱的响应均呈现为随着折合流速增大而增大的弛振现象,平衡位置偏移、低频振动以及旋涡脱落与圆柱运动之间的时机三个因素共同决定了弛振现象的出现.进一步的研究发现,串列三圆柱的弛振现象仅出现在质量比不大于2.0和雷诺数不大于100的工况下.当质量比较大时,串列三圆柱的平衡位置固定不变,且圆柱的振动不规律,使得旋涡脱落与圆柱运动的时机处于变化之中.当雷诺数较高时,最上游圆柱的平衡位置在折合流速较大时回到初始位置,不再参与对圆柱振动的调节,使得圆柱的振动响应不再规律,旋涡脱落与圆柱运动的时机也一直处于变化之中.     

倾斜圆柱结构涡激振动研究进展

圆柱结构涡激振动现象在生活中十分常见, 如海洋工程中的管道、土木工程中的高耸建筑、桥梁斜拉索, 核工程中的热交换器等频繁受到涡激振动影响, 诱发结构的疲劳损伤, 甚至破坏失效. 现阶段, 人们对垂直来流作用下圆柱结构涡激振动机理已有较为全面的认识. 然而, 当圆柱倾斜置于流场中, 结构后缘的尾流形态与垂直放置差异显著, 结构与流体的耦合作用机理更为复杂. 为简化倾斜圆柱涡激振动问题, 提出了不相关原则, 来流速度被分解为垂直圆柱结构轴向和平行圆柱结构轴向的两个速度分量, 仅考虑垂直结构轴向速度分量的影响, 忽略平行结构轴向速度分量的影响. 近年来, 针对倾斜圆柱涡激振动及不相关原则的适用性, 出现了大量实验和数值模拟研究成果. 为了深化对倾斜圆柱结构涡激振动相关机理的认知, 本文全面阐述了倾斜圆柱结构涡激振动响应规律、尾迹流场模式和流体力特性等方面的研究进展, 分析了不相关原则的适用范围, 探讨了倾斜圆柱结构涡激振动抑制措施, 并对今后该领域的研究进行了力所能及的展望.      

虚拟边界法研究正交双圆柱及串列双圆球绕流

把Goldstein等人提出的虚拟边界法推广到三维情况，研究了 Re=150时不同间距下正交双圆柱绕流，和Re=250时不同间距下串列双 圆球绕流流场. 对于正交双圆柱绕流，当间距比大于3，下游圆柱对上游圆柱尾流的影响只 限定在下游圆柱的尾流所扫过的范围之内；当间距比小于等于3，下游圆柱对上游圆柱尾流 的影响扩大，下游圆柱尾流扫过区上下出现两排三维流向二次涡结构. 对于串列圆球绕流， 研究发现，在小间距比(L/D≈ 1.5)的情况下，由于上下游圆球尾流区的相互抑 制消除了压力不稳定性，整个流场呈现稳 态轴对称特征；间距比为2.0时，周向压力梯度诱发出流体的周向输运，流场呈现稳态非对 称性，但流场中存在特定的对称面；间距比增大到2.5后，绕流场开始周期振荡，原有的对 称面依旧存在；在间距比3.5时下游圆球下表面的涡结构强度有所减弱，导致占优频率发生 交替；间距比增至7.0时，整个流场恢复稳态特征，两圆球尾部同时出现双线涡，这时流场 对称面的位置发生了变动.     

不同剪切率来流作用下柔性圆柱涡激振动数值模拟

采用浸入边界法对细长柔性圆柱在线性剪切流条件下的涡激振动进行三维数值模拟。细长柔性圆柱振动采用三维索模型模拟,其两端铰接,质量比为6,长细比为50,无量纲顶张力为496。来流为线性剪切流,剪切率从0到0.024变化,最大雷诺数为250。研究发现:剪切流作用下柔性立管横流向振动表现为驻波模式,而顺流向振动表现为行波-驻波混合模式。随着剪切率增大,振动频谱呈现多频响应,振动能量逐渐向低频转移。阻力系数平均值随着展向变化,脉动阻力系数和升力系数的均方根值均表现为“双峰”模式。流固能量传递系数沿立管轴向的分布表明,振动激励区集中于高流速区,而振动阻尼区多位于低流速区。剪切率较小时,圆柱的泻涡为平行交叉模式;剪切率较大时,圆柱的泻涡为倾斜泻涡模式,且由于泻涡频率沿立管轴向变化,尾流发生涡裂现象,形成泻涡频率不同的胞格结构。      

串列双圆柱绕流下游圆柱两自由度涡致振动研究

数值研究了串列双圆柱绕流下游圆柱两自由度涡致振动问题，研究发现：(1) 双自由度的圆柱振幅峰值及出现振峰的频率比都比单自由度的大；(2) 尾流圆柱中的升力远大于均匀来流的，而阻力却相反；(3) 下游圆柱的位移响应对于频率比的变化没有均匀来流中的"敏感"；(4) 尾流中，在频率比1.16和0.87之间，出现了明显的"拍"现象，即圆柱的振幅响应包含不同的频率，而在均匀来流中，并无明显的"拍"现象. 采用ALE方法，计算网格采用H-O非交错网格系统，结合分块耦合方法. N-S方程的对流项和扩散项分别采用三阶迎风紧致格式和四阶中心紧致格式离散. 圆柱振动采用弹簧柱体阻尼器模型，柱体的振动方程采用龙格-库塔法求解. 通过模拟柱体和流体之间的非线性耦合作用，成功地捕捉到了"拍"和"相位开关"等现象.     

流致振动能量收集的钝头体几何设计研究

流致振动蕴含着可观的能量, 通过能量收集技术可将其转化为电能. 为提高低速流场中能量转化效率, 本文实验研究了不同截面下钝头体以及它们的宽厚比(W/T)对流致振动能量收集特性的影响, 并通过计算流体动力学(computational fluid dynamic, CFD)仿真分析了尾流特性. 流致振动能量收集装置由压电悬臂梁和不同截面的钝头体构成. 首先搭建了流致振动能量收集风洞实验平台, 钝头体的截面分别设置为矩形、三角形和D形, 宽厚比分别设定为1, 1.3, 1.8和2.5. 然后利用实验方法分析不同形状钝头体的宽厚比(W/T)对位移响应和电压响应的影响规律. 最后通过计算流体动力学模拟揭示实验结果的内在力学机理. 实验结果表明, 当钝头体截面为矩形时, 增大宽厚比可以显著提高电压输出峰值; 当钝头体为三角形和D形时, 增加宽厚比将使系统呈现“驰振”→“驰振 + 涡激振动”→“涡激振动”响应特性变化趋势, 提高了低风速时的能量收集效果. CFD结果解释了实验现象, 即随着宽厚比增加, 钝头体尾流会产生更加强劲的涡街, 显著提高流致振动能量收集效果. 相关结果可优化流致振动能量收集装置结构, 为提高低速流场的能量收集效果提供理论和实验依据.      

NUMERICAL INVESTIGATION ON FLOW-INDUCED VIBRATION OF TWO CYLINDERS IN TANDEM ARRANGEMENTS AND ITS COUPLING MECHANISMS

对雷诺数Re= 100 条件下串列双圆柱的流致振动进行了数值模拟. 圆柱的质量比m^*均为2.0,间距比L/D 为2.0 5.0. 考虑两种工况:(a) 上游圆柱固定,下游圆柱可沿横流向自由振动;(b) 上、下游圆柱均可沿横流向自由振动. 结果表明:无论上游圆柱静止或者振动,下游圆柱横向振幅明显大于单圆柱的. 工况(b) 的下游圆柱最大振幅要大于工况(a) 的,这是由于两圆柱均振动时,两圆柱之间耦合作用增强,上游圆柱的尾流和下游圆柱的振动之间“相互调节” 作用显著. 对工况(b) 的下游圆柱振动和间隙流之间的作用机制进行了详细的研究,发现当上游圆柱脱落的自由剪切层重新附着于下游圆柱上并且完全从间隙之间通过时,下游圆柱的振幅最大.     

Wake flow pattern modified by small control cylinders at low Reynolds number

Passive wake control behind a circular cylinder in uniform flow is studied by numerical simulation for Re_D ranging from 80 to 300. Two small control cylinders, with diameter d/D=1/8, are placed at x/D=0.5 and y/D=±0.6. Unlike the 1990 results of Strykowski and Sreenivasan, in the present study, the vortex street behind the main cylinder still exists but the fluctuating lift and the form drag on the main cylinder reduces significantly and monotonously as the Reynolds number increases from 80 to 300. Obstruction of the control cylinders to the incoming flow deflects part of the fluid to pass through the gap between the main and control cylinders, forming two symmetric streams. These streams not only eliminate the flow separation along the rear surface of the main cylinder, they also merge toward the wake centerline to create an advancing momentum in the immediate near-wake region. These two effects significantly reduce the wake width behind the main cylinder and lead to monotonous decrease of the form drag as the Reynolds number increases. As the Reynolds number gets higher, a large amount of the downstream advancing momentum significantly delays the vortex formation farther downstream, leading to a more symmetric flow structure in the near-wake region of the main cylinder. As the Reynolds number increases from 80 to 300, both increasing symmetry of the flow structure in the near-wake and significant delay of the vortex formation are the main reasons for the fluctuating lift to decrease monotonously.     

Vortex induced vibration and vortex shedding characteristics of two side-by-side circular cylinders of different diameters in close proximity in steady flow

This paper presents results of a numerical study of vortex-induced vibrations of two side-by-side circular cylinders of different diameters in steady incompressible flow. The two-dimensional Reynolds-averaged Navier–Stokes equations with a SST k–ω turbulence model are solved using the Petrov–Galerkin finite element method and the Arbitrary-Lagrangian–Eulerian scheme. The diameter ratio of the two cylinders is fixed at 0.1 and the mass ratio of both cylinders is 5.0. Both cylinders are constrained to oscillate in the transverse direction only. The Reynolds number based on the large cylinder diameter and free stream velocity is fixed at 5000. The effects of the reduced velocities of the cylinders on the vibration amplitude and vortex shedding regimes are investigated. It is found that for the range of parameters considered, collision between the two cylinders is dependent on the difference of the reduced velocities of the cylinders. Presence of the small cylinder in the proximity of the large one appears to have significant effects on the vortex shedding regime and vibration amplitude of the large cylinder.     

柔性圆柱涡激振动流体力系数识别及其特性

涡激振动是诱发海洋立管、浮式平台系泊缆和海底悬跨管道等柔性圆柱结构疲劳损伤的重要因素.目前,海洋工程中用于柔性圆柱涡激振动预报的流体力系数主要来源刚性圆柱横流向受迫振动的实验数据,存在一定缺陷和误差.本文综合考虑横流向与顺流向振动耦合作用,建立了柔性圆柱涡激振动流体力模型,运用有限元法和最小二乘法确定升力系数、脉动阻力系数和附加质量系数.为了准确识别柔性圆柱涡激振动流体力系数,设计并开展了拖曳水池模型实验,实验用柔性圆柱模型的质量比为1.82,长径比为195.5.通过与刚性圆柱流体力系数对比,深入分析了柔性圆柱流体力系数的特性.结果表明:柔性圆柱在一阶模态控制区,流体力系数随约化速度变化趋势与刚性圆柱大致相似;二阶模态控制区,升力系数和脉动阻力系数显著增大;附加质量系数在响应频率较低时与振动位移的相关性增强;当响应频率较低时,振动位移较大区域为能量耗散区,当响应频率较高时,振动位移较大区域为能量输入区.     

串列布置三圆柱涡激振动频谱特性研究

对串列三圆柱体双自由度涡激振动问题进行了数值计算, 并分析了雷诺数、固有频率比和约化速度对串列三圆柱体结构动力响应及频谱特性的影响. 研究发现: 雷诺数、频率比对上游圆柱的振幅和流体力系数的影响较小. 中游圆柱频率锁定区域随着雷诺数的增大而增大, 其动力响应受上游圆柱尾流的影响较大, 但频率比的影响较小. 同时, 流体力系数在约化速度较小时受雷诺数和频率比的影响较大. 另外, 下游圆柱的振幅和流体力系数受雷诺数及频率比的影响较大. 雷诺数、频率比和约化速度对圆柱流体力系数能量谱密度(PSD)曲线中主峰幅值、频谱成分及波动性的影响较大. 流体力系数PSD曲线波动性的增强, 导致圆柱运动轨迹会从"8"字形转变成不规则形状. 当频率比为2.0时, 上游圆柱尾流出现P$+$S模式, 导致其发生非对称运动, 且升、阻力系数PSD曲线主峰重合. 最后, 激励荷载平均功率值随约化速度的变化趋势与对应的结构动力响应的变化类似. 在同一约化速度区间内, 结构振动响应的强弱与位移的平均功率值成正比. 对不同约化速度区间内的升力系数功率谱密度分析时, 振动频率比($f_{s}/f_{n, y})$对结构振动响应的影响更大.      

带有新型涡激振动抑制罩的圆柱体的水动力特性

通过模型实验和数值模拟计算,研究了带有涡激振动抑制罩的圆截面柱体的水动力特性.模型实验主要测试了柱体上附加谐波型和类圆锥型涡激振动抑制罩的单摆结构在不同流速下发生涡激振动的性质;数值模拟则针对谐波型和圆锥型扰动,在雷诺数Re为102到105范围内,研究其水动力参数,如阻力、升力和涡脱落频率等,随扰动波长和波动强度的变化.模型实验结果表明,在直圆柱开始发生共振的流速下,带抑制罩的柱体的振幅显著降低,而在更高流速下则显著增大.数值模拟结果表明,谐波型和圆锥型扰动具有相似的水动力特性;且在不同Re时,阻力、升力和涡脱落频率具有相似的变化规律;随波动强度的增大,阻力一般逐渐增大,升力则在多数情况下先减小而后增大,而涡脱落频率一般逐渐减小.      

On characteristics of two-degree-of-freedom vortex induced vibration of two low-mass circular cylinders in proximity at low Reynolds number

Two-degree-of-freedom vortex induced vibrations (VIVs) of two identical spring-supported circular cylinders in proximity with the mass ratio of 2 and zero damping at Re of 100 are numerically studied. Totally 20 arrangements of cylinders are investigated combining four stagger angles and five normalized center-to-center spacings. Results show that the in-line vibration amplitude is comparable to the transverse one for most arrangements and usually accompanies irregular cylinder trajectories. Extremely slender figure-8 cylinder trajectories usually seen in single-cylinder VIVs exist only for the tandem arrangements. Arranging the trailing cylinder to vibrate near the wake boundary of the leading cylinder enhances the possibility of irregular trajectories and impacts of both cylinders. Impact between cylinders must occur in cases with irregular cylinder trajectories; however, irregular cylinder trajectories could be found in impact-free cases. The stagger angle significantly changes the attribute of the transverse vibration frequency, toward either the single-cylinder VIV frequency or natural structure frequency in still fluid. The major transverse vibration frequency and the natural structure frequency in still fluid are decoupled for all the side-by-side arrangements and some far spaced tandem arrangements and highly correlated for non-tandem and non-side-by-side arrangements. The time-averaged impact frequency increases with decreasing normalized center-to-center spacing for most combinations of stagger angle and reduced velocity. Apart from the side-by-side arrangements, high-frequency impacts occur when the trailing cylinder is initially located in or near the wake zone of the leading cylinder. The mechanism of trailing cylinder chopping the gap-flow vortices plays an important role in determining the near-wake vortex structures for all non-side-by-side arrangements.     

Characteristics and suppression of flow-induced vibrations of two side-by-side circular cylinders

A free-vibration experiment was conducted to examine flow-induced vibration (FIV) characteristics of two identical circular cylinders in side-by-side arrangements at spacing ratio T^⁎ (=T/D)=0.1–3.2, covering all possible flow regimes, where T is the gap spacing between the cylinders and D is the cylinder diameter. Each of the cylinders was two-dimensional, spring mounted, and allowed to vibrate independently in the cross-flow direction. Furthermore, an attempt to suppress flow-induced vibrations was undertaken by attaching flexible sheets at the rear stagnation lines of the cylinders. Based on the vibration responses of the two cylinders, four vibration patterns I, II, III and IV are identified at 0.1≤T^⁎<0.2, 0.2≤T^⁎≤0.9, 0.9<T^⁎<2.1 and 2.1≤T^⁎≤3.2, respectively. Pattern I is characterized by the two cylinders vibrating inphase^, with the maximum amplitudes occurring at the same reduced velocity Ur=10.47 almost two times that (Ur=5.25) for an isolated cylinder. Pattern II features no vibration generated for either cylinder. Pattern III exemplifies the occurrence of the maximum vibration amplitude of a cylinder at a smaller Ur and that of the other cylinder at a higher Ur compared to its counterpart in an isolated cylinder. Pattern IV represents each cylinder response resembling an isolated cylinder response; the vibrations of the two cylinders are, however, coupled inphase or antiphase. Linking maximum vibration amplitudes to fluctuating lift forces acting on fixed cylinders reveals that fluid–structure interactions between two fixed cylinders and between two elastic cylinders are not the same, even though vibration is not significant. As such, while two fixed cylinders generate narrow and wide wakes at 0.2≤T^⁎<1.7, two elastic cylinders do the same for a longer range of T^⁎ (0.2≤T^⁎<2.1). The flexible sheets effectively suppress FIV of the two cylinders in patterns III and IV, and reduce the vibration amplitude in pattern I. For the effectively controlled cases (patterns III and IV), the flexible sheet of each cylinder folds into a semicircle at the base, forming two free edges.     

TWO-DIMENSIONAL VORTEX-INDUCED VIBRATION OF CABLE SUSPENSIONS

Resonant responses of suspended elastic cables driven by a steady current are investigated. Phenomenological fluid force models for alternate vortex-shedding are coupled with the nonlinear partial differential equations of cable motion. Decoupled cross-flow and in-line vortex-induced vibrations (VIV) are examined first using linearized and nonlinear cable models. The linearized cable model predicts well the basic characteristics of VIV and the nonlinear cable model captures the hysteresis often observed in experiments. Next, coupled cross-flow and in-line vibrations are evaluated by considering two principal coupling mechanisms: (i) cable structural nonlinearities, and (ii) coupled fluid lift and drag. Attention is focused on a “worst-case” resonant response where the natural frequencies for cable modes in the cross-flow and in-line directions are in the same 1:2 ratio as the excitation frequencies associated with lift and drag. The inclusion of cable structural nonlinearities alone leads to coupled responses that differ qualitatively (i.e., in number and stability of periodic motions) when compared to those of the decoupled model. The inclusion of coupled fluid lift and drag produces non-planar “figure eight” motions of the cable cross-section that exhibit similar characteristics to those previously measured on spring supported cylinders.