并列双圆柱流致振动的不对称振动和对称性迟滞研究

对雷诺数Re = 100 间距比s/D = 2.5 和5.0 的并列双圆柱流致振动进行了数值模拟研究, 其中圆柱质量比m = 2.0, 折合流速U_r 在2.0~10.0 之间, 两圆柱仅能做横流向振动. 研究发现, 当间距比s/D = 2.5 时, 在折合流速4.4 < U_r< 4.8区间内, 两圆柱流致振动响应出现不对称振动现象, 在折合流速4.4 < U_r< 4.8 区间内, 两圆柱流致振动响应出现对称性迟滞现象; 而当间距比s/D = 2.5时, 圆柱流致振动响应与单圆柱涡激振动响应相似, 没有出现不对称振动和对称性迟滞现象. 在不对称振动区间内, 两圆柱的升、阻力参数也出现了不相等的情况. 此外, 当两圆柱不对称振动时, 圆柱间隙流稳定地偏斜向其中的一个圆柱; 相应地, 尾涡也出现了宽窄不等的模式. 窄尾流圆柱的振幅和升、阻力均较宽尾流圆柱的大. 通过对比不对称振动现象发生前后的尾涡模式, 对新现象的产生机制进行了阐述.      

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

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

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

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

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

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

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

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

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

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

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) 的下游圆柱振动和间隙流之间的作用机制进行了详细的研究,发现当上游圆柱脱落的自由剪切层重新附着于下游圆柱上并且完全从间隙之间通过时,下游圆柱的振幅最大.     

两交错排列弹性管的流致振动特性研究

为了深入研究反应堆结构中诸如燃料棒、蒸汽发生器传热管束等的流致振动问题,利用有限体积法离散大涡模拟的流体控制方程,利用有限元方法离散结构动力学方程,并结合动网格技术建立了模拟双向流固耦合作用的三维数值模型,实现了弹性管与弹性管之间以及弹性管与湍流流场之间的交互作用。利用该模型分别研究了节径比P/D为1.2、1.6、2.0、3.0、4.0的两交错排列弹性管在横向湍流作用下的振动响应及流场特性。通过分析发现:两交错排列弹性管的运动轨迹和尾涡结构与其间距、流速密切相关;两交错管的临界节径比为2.0;P/D=3.0时双管的临界流速Upr=2.2,P/D=1.6时双管的临界流速Upr=4.0。     

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

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

不同控制角下附加圆柱对圆柱涡激振动影响

在弹性支撑的圆柱周围布置直径更小圆柱会影响剪切层发展以及旋涡脱落,进而改变其涡激振动状态.通过不同的布置形式和附加小圆柱个数可以实现对圆柱涡激振动的促进或抑制.激励更大幅值的振动可以更好地将水流动能转化为可利用的机械能或电能,抑制其振动则可以实现对海洋平台等结构物的保护.采用基于迭代的嵌入式浸入边界法对前侧对称布置两个小圆柱的圆柱涡激振动进行数值模拟研究,系统仅做横向振动,其中基于主圆柱直径的雷诺数为100,质量比为2.0,折合流速为3~11.小圆柱与主圆柱的直径比为0.125,间隙比为0.125.结果表明,在研究的控制角范围内(30°~90°),附加小圆柱可以很大程度上改变圆柱涡激振动的状态.当控制角较小(30°)时,附加小圆柱对主圆柱的振动起抑制作用;当控制角为45°~60°时,圆柱的振动分为涡振和弛振两个阶段,在弛振阶段,圆柱振幅随折合流速增加而持续增加;当控制角较大(75°~90°)时,附加小圆柱的促进作用随着控制角增加而减小.进一步地,结合一个周期内不同时刻旋涡脱落以及圆周压强分布,解释了附加小圆柱对主圆柱涡激振动的作用机制.应用能量系数对圆柱系统的进一步分析发现,弛振阶段由流体传递到主圆柱的能量系数随折合流速的增加逐渐下降,旋涡结构的改变是产生这种变化的直接原因.      

Flow-induced vibrations of two circular cylinders in tandem with shear flow at low Reynolds number

The flow-induced vibrations of two elastically mounted circular cylinders subjected to the planar shear flow in tandem arrangement are studied numerically at Re=160. A four-step semi-implicit Characteristic-based split (4-SICBS) finite element method is developed under the framework of the fractional step method to cope with the vortex-induced vibration (VIV) problem. For the computational code verification, two benchmark problems are examined in the laminar region: flow-induced vibration of an elastically mounted cylinder having two degrees of freedom and past two stationary ones in tandem arrangement. Regarding the two-cylinder VIVs in shear flow, the computation is conducted with the cylinder reduced mass M_r=2.5π and the structural damping ratio ξ=0.0. The effects of some key parameters, such as shear rate (k=0.0, 0.05, 0.1), reduced velocity (U_r=3.0–18.0) and spacing ratio (L_x/D=2.5, 3.5, 4.5, 8.0), are demonstrated. It is observed that the shear rate and reduced velocity play an important role in the VIVs of both cylinders at various center-to-center distances. Additionally, in comparison with the single cylinder case, a further study indicated that the gap flow has a significant impact on such a dynamic system, leading it to be more complex. The results show that, the performances of ‘dual-resonant’ are discovered in the shear flow. A valley is formed in transverse oscillation amplitude of DC for each spacing ratio when U_r is about 6.0. For the X–Y trajectories of the circular cylinders, figure-eight, figure-O and oval shape are obtained. Finally, the interactions between cylinders are revealed, together with the wake-induced vibration (WIV) mechanism underlying the oscillation characteristics of both cylinders exposed to shear flow. Besides, the “T+P” wake pattern is discovered herein.     

Response and wake patterns of two side-by-side elastically supported circular cylinders in uniform laminar cross-flow

Vortex-induced vibrations (VIV) of two side-by-side elastically supported circular cylinders in a uniform flow with the Reynolds number of 100 are numerically investigated by using the immersed boundary method. The cylinders are constrained to oscillate in the cross-flow direction with a center-to-center spacing ratio $T/D$ ranging from 2 to 5. The structural damping is set to zero to enable large vibration amplitudes in the range of reduced velocity $U_r=3-10$. It is found that the proximity of the cylinders does not have a significant impact to the lock-in region and cylinder responses, except at a small spacing ratio of $T/D=2$. The critical spacing ratio is determined as $T/D=4$ and beyond that the interaction between the cylinders is negligible. The following six near-wake patterns are observed; the irregular pattern, in-phase flip-flopping pattern, out-of-phase flip-flopping pattern, in-phase-synchronized pattern, antiphase-synchronized pattern and the biased antiphase-synchronized pattern. These patterns are plotted in a plane of $U_r$ and $T/D$, together with approximate borderlines to distinguish one region from the others. The time histories, spectral features and wavelet transform contours of drag and lift forces are presented to elucidate the mechanisms of the in-phase and out-of-phase flip-flopping phenomena. It is established that the in-phase flip-flopping stems from the long-short near-wake pattern and its low-frequency flip-over, whereas the out-of-phase pattern originates from the large vortex shedding from the fictitious bluff-body with an augmented characteristic length.     

Free-stream turbulence effects on vortex-induced vibration of two side-by-side elastic cylinders

The effect of free-stream turbulence on vortex-induced vibration of two side-by-side elastic cylinders in a cross-flow was investigated experimentally. A turbulence generation grid was used to generate turbulent incoming flow with turbulence intensity around 10%. Cylinder displacements in the transverse direction at cylinder mid-span were measured in the reduced velocity range 1.45<U_r0<12.08, corresponding to a range of Reynolds number (Re), based on the mean free-stream velocity and the diameter of the cylinder, between Re=5000–41 000. The focus of the study is on the regime of biased gap flow, where two cylinders with pitch ratio (T/D) varying from 1.17 to 1.90 are considered. Results show that the free-stream turbulence effect is to enhance the vortex-induced force, thus to restore the large-amplitude vibration associated with the lock-in resonance. However, the enhancement is significant at a different Strouhal number (St) for different pitch ratios. When the spacing between two cylinders is relatively small (1.17<T/D<1.50), the enhancement is significant at St≈0.1. When the spacing is increased, the Strouhal number at which the enhancement is significant shifts to St≈0.16. This enlarges the range of reduced velocity to be concerned. An energy analysis showed that free-stream turbulence feeds energy to the cylinder at multiple frequencies of vortex shedding. Therefore, the lock-in region is still of main concern when the approach flow is turbulent.     

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.     

Flow interference between a stationary cylinder and an elastically mounted cylinder arranged in proximity

The flow interference between two circular cylinders, one stationary and the other free to oscillate in the transverse direction, are studied numerically at Re=150. The incompressible Navier–Stokes equation in two space dimensions, an assumption that is expected to be valid at the considered Re, is solved by the characteristic-based-split (CBS) finite element method using the T4/C3 MINI triangular element. The center-to-center spacing between the two cylinders is fixed at 4D, where D is the cylinder diameter. The angle of incident flow, α, with respect to the line through the two cylinder centers, varies within the range from α=0° to 90°. For the elastically mounted cylinder, the reduced mass considered is M_r=2.0; the structural damping coefficient is assigned to be zero, which encourages high amplitude oscillations. For each α, the computations are conducted for a wide range of reduced velocities, U_r. The flow interference is examined by scrutinizing (i) the frequency characteristics of the vortex shedding and oscillation; (ii) the dynamic response of the oscillating cylinder, including the amplitude of displacement, the drag and lift force characteristics and the phase relationship between the lift and the displacement series; and (iii) the flow response in terms of the instantaneous vorticity field. It was found that the flow interference type is significantly affected by the angle of the incident flow. As the cylinder is oscillated outside of the region of the wake behind the stationary cylinder (α≥30°), it behaves similarly to its isolated counterpart. In contrast, if the cylinder is partially or entirely submerged within the upstream wake (α<30°), then both the flow and body responses are substantially modified due to the vigorous interaction between the upstream wake and the oscillating cylinder; the response therefore belongs to the wake-induced regime. The U_r range associated with the higher amplitude response is significantly shifted toward a higher U_r. The maximum vibration amplitude builds up to a significantly higher level, even increasing the U_r far beyond the resonance regime. In general, the wake flow associated with the wake-induced vibration (WIV) regime appears to be more unperiodic than does that corresponding to the vortex-induced oscillation regime. It is also revealed that both the vortex-cylinder and the shear layer-cylinder interaction mechanisms are responsible for the characteristics of the responses of oscillating cylinder. The larger momentum required for the higher oscillation amplitude is obtained from the duration of the energy transfer from the fluid to the cylinder, which is ascribed to the phase lag between the lift force and the cylinder displacement.     

Wake-induced vibrations of an elastically mounted cylinder located downstream of a stationary larger cylinder at low Reynolds numbers

Two-dimensional numerical simulations of flow past two unequal-sized circular cylinders in tandem arrangement are performed at low Reynolds numbers (Re). The upstream larger cylinder is stationary, while the downstream cylinder has both one (transverse-only) and two (transverse and in-line) degrees of freedom (1-dof and 2-dof, respectively). The Re, based on the free stream velocity U_∞ and the downstream cylinder diameter d, varies between 50 and 200 with a wide range of reduced velocities U_r. The diameter of the upstream cylinder is twice that of the downstream cylinder, and the center-to-center spacing is 5.5d. In general, for the 1-dof case, the calculations show that the wake-induced vibrations (WIV) of the downstream cylinder are greatly amplified when compared to the case of a single cylinder or two equal-sized cylinders. The transverse amplitudes build up to a significantly higher level within and beyond the lock-in region, and the U_r associated with the peak amplitude shifts toward a higher value. The dominant wake pattern is 2S mode for Re=50 and 100, while with the increase of Re to 150 and 200, the P+S mode can be clearly observed at some lower U_r. For the 2-dof vibrations, the transverse response characteristics are similar to those presented in the corresponding 1-dof case. The in-line responses are generally much smaller, except for several significant vibrations resulting from in-line resonance. The obvious in-line vibration may induce a C (chaotic) vortex shedding mode for higher Re (Re=200). With regard to the 2-dof motion trajectories, besides the typical figure-eight pattern, several odd patterns such as figure-double eight and single-looped trajectories are also obtained due to the wake interference effect.     

Three-dimensional numerical simulations of cross-flow around four cylinders in an in-line square configuration

This paper presents a numerical study of three-dimensional (3-D) laminar flow around four circular cylinders in an in-line square configuration. The investigation focuses on effects of spacing ratio (L/D) and aspect ratio (H/D) on 3-D flow characteristics, and the force and pressure coefficients of the cylinders. Extensive 3-D numerical simulations were performed at Reynolds number of 200 for L/D from 1.6 to 5.0 at H/D=16 and H/D from 6 to 20 at L/D=3.5. The results show that the 3-D numerical simulations have remedied the inadequacy of 2-D simulations and the results are in excellent agreement with the experimental results. The relation between 3-D flow patterns and pressure characteristics around the four cylinders is examined and discussed. The critical spacing ratio for flow pattern transformation was found to be L/D=3.5 for H/D=16, while a bistable wake pattern was observed at L/D=1.6 for the same aspect ratio. Moreover, a transformation of flow pattern from a stable shielding flow pattern to a vortex shedding flow pattern near the middle spanwise positions of the cylinders was observed and was found to be dependent on the aspect ratio, spacing ratio, and end wall conditions. Due to the highly 3-D nature of the flows, different flow patterns coexist over different spanwise positions of the cylinders even for the same aspect ratio. It is concluded that spacing ratio, aspect ratio, and the no-slip end wall condition have important combined effects on free shear layer development of the cylinders and hence have significant effects on the pressure field and force characteristics of the four cylinders with different spacing ratios and aspect ratios.