二维圆柱涡激振动的数值模拟

采用动网格的动态层铺模型和滑移网格模型以及用户自定义接口编程,将计算结构响应的Newmark-β代码嵌入FLU-ENT软件,建立了2D圆柱横流向涡激振动(VIV)的计算模型。在雷诺数6×103~2×104或对应的约化速度Vr=3~10范围内,结合k-ω湍流模型,并且以静止圆柱为初始状态,首先计算了质量比m*=12.8、255.5和结构阻尼比ζ=0.0%、0.5%几种情况下圆柱的VIV。通过与经典的Feng-type图示比较发现,在小"组合参数"(m*ζ)情况下,计算结果与实验的"初始"支和最大响应吻合较好,同时观察到"初始"支与"下"支之间的"跳跃"现象。通过另外两组工况的计算发现,m*与ζ对圆柱振幅的影响不同,前者是完全非线性的,而后者几乎是线性的。     

单自由度圆柱涡激振动特性及其影响因素分析

对二维刚性圆柱涡激振动进行了系统理论和数值模拟研究,对比了不同因素下圆柱振幅、锁定区间、升阻力系数等关键参数和特征,探究了影响圆柱涡激振动特性的因素。通过线性化分析推导出了锁定区间柱体响应的频率比、振幅与流体荷载和系统物性参数之间的定性关系。为检验理论分析结果,对相关经典实验的结果进行了整理,并用RANS模型开展了相应的数值模拟。理论分析、实验和数值模拟三者的结论相符,可以确认以下结论:当圆柱处在锁定状态时,其振幅大小由组合因子决定。组合因子越小,圆柱的振幅越大,反之则越小,但是组合因子对振幅的影响幅度有限;柱体锁定区间范围的大小由质量比决定,质量比越小,锁定范围越大,反之则越小,质量比对锁定区间范围的大小有显著影响,影响的主要范围是在约化速度较大一端。     

并列双圆柱涡激振动的经验性模型研究

建立一种新的预报并列双圆柱涡激振动响应的经验性模型,根据特定间距比条件下旋涡脱落频率出现分支的现象,提出以两个具有不同固有频率的尾流振子来共同描述结构的近壁尾涡动力特性,同时两个振子均满足van der Pol方程,进而得到结构振子和流体振子的耦合方程组。使用该模型分别对中高质量比和低质量比的并列圆柱涡激振动问题进行数值计算,结果表明,结构的位移响应和最大振动幅值等变化规律与实验结果趋势一致,数值基本吻合。     

结构与尾流非线性耦合涡激振动预测模型

预测圆柱涡激振动的尾流振子模型中,通常采用线性的耦合模型,例如位移或者速度、加速度耦合来表征结构对尾流的作用。三种线性模型在预测圆柱锁频阶段的动力特性时存在差异,而且适用范围也受质量比的限制。提出了考虑结构与尾流动力非线性耦合的模型,该模型基于加速度耦合并结合速度耦合进行修正,适用范围不受质量比的影响;与实验结果的对比表明该模型可以更合理地给出锁频区域以及结构位移响应和尾流升力。最后,利用新模型讨论了质量比对锁频阶段结构振动幅值、尾流升力及频率比的影响;结果表明,随着质量比的增大,结构锁频区域变窄,结构振幅和尾流升力幅值减小。     

分离盘控制圆柱涡激振动的数值模拟研究

圆柱涡激振动问题一直以来备受关注,分离盘作为涡激振动抑制装置得到广泛研究。分离盘长度L与圆柱直径D之比L/D是影响抑制效果的主要因素。运用有限体积法结合RANS方程与一定的湍流模式离散和求解流场,通过编写自定义程序,使用动网格模拟结构物的运动带来的流域边界的变化,针对弹性支撑的圆柱及附加长度为0.5 D的分离盘模型,在约化速度Ur为2.5~13的情况下,对涡激振动及其抑制进行研究。结果表明:分离盘可以抑制甚至消除圆柱涡激振动,99%以上的振幅被抑制;锁定区始点被推后,锁定区变窄;附加分离盘的圆柱阻力和升力被抑制;其斯特鲁哈数(St)稍高于单圆柱St但差别不大。     

低雷诺数下圆柱涡激振动的二维有限元数值模拟

采用有限元方法求解原始变量的二维不可压粘性流体的N-S方程,计算了雷诺数从90到150范围内圆柱绕流引起的涡激振动,完整地再现了流固耦合系统从不共振到频率锁定,再到脱离锁定的过程,成功地预测到了涡激振动的“拍”和“锁定”现象,并与A nagnostopou los和B earam an的试验结果进行了比较。计算涡激振动时用ALE方法分析圆柱和流体的耦合作用,圆柱振动被简化为质量-弹簧-阻尼系统。     

不同运动自由度组合串列双圆柱涡激振动数值模拟

基于自主研发的紧致插值曲线CIP (Constrained Interpolation Profile)方法数学模型,对均匀来流条件时不同运动自由度组合下的串列双圆柱涡激振动问题开展二维数值模拟。模型针对雷诺数Re=100,质量比m*=2的串列双圆柱涡激振动问题,选取上、下游圆柱不同运动自由度组合工况进行模拟。重点分析圆柱的升阻力系数、运动位移随折合速度Ur变化的响应。研究表明:当上游圆柱双自由度运动时,随着下游圆柱运动自由度的增加,下游圆柱对上游圆柱涡激振动响应的影响减弱;当下游圆柱双自由度运动时,随着上游圆柱运动自由度的增加,上游圆柱对下游圆柱涡激振动响应的影响变强。研究结果表明圆柱运动自由度组合形式对串列双圆柱涡激振动的影响不可忽略。     

质量比对柔性立管涡激振动影响实验研究

质量比是影响海洋立管涡激振动的一个重要因素.通过在室内物理实验中使立管模型内部分别充填空气、水和沙来改变立管的质量比,从而研究质量比对柔性细长立管涡激振动的影响.实验结果表明:在相同流速下,质量比大的立管模型所激起的模态更高.在低约化速度区域,空管和水管的涡激振动响应频率与涡脱落频率相同,沙管的响应频率则与自振频率更接近,三种质量比立管的响应位移较接近;在高约化速度区域,三种质量比的立管模型的响应频率处于自振频率和涡脱落频率之间,但空管的响应频率随约化速度的增大而不断增大,同一流速下,质量比大的立管模型响应位移小,其中空管的涡激振动响应一直处于大振幅的锁定状态下.共振区域对应约化速度的范围随着质量比增大而减小.     

不同固有频率比条件下圆柱双自由度涡激振动特性的数值研究

采用RANS方法,结合SST k-ω湍流模型,对不同顺流向与横顺流向固有频率的比值(即固有频率比,fnx/fny)条件下低质量比圆柱体的双自由度涡激振动进行了二维数值模拟。圆柱体的质量比为2.6,雷诺数范围为2 500~18 750,相应的约化速度范围为2~15,包括了经典试验中出现的整个锁定范围。通过研究发现,固有频率比是影响振动特性的重要参数,随着固有频率比的增加,响应幅值逐渐降低且向更高的约化速度偏移;在低约化速度范围内,固有频率比对顺流向和横流向振动之间的相位差以及升力频率有较大影响,从而得到各种不同偏向的8字形轨迹;最后对不同固有频率比条件下的尾涡模式进行了讨论,给出了对应不同约化速度时的尾涡模式。     

大长细比柔性杆件涡激振动实验

涡激振动(vortex-induced vibration,VIV)是导致深海细长柔性立管发生疲劳破坏的重要因素。采用实验观测手段研究了长细比为1 750的柔性立管多模态涡激振动特性。实验中,通过采用拖车拖拉立管模型在水池中匀速行进来模拟均匀流作用下的涡激振动响应。利用光纤光栅传感器测量立管模型在横流向(cross-flow,CF)和顺流向(in-line,IL)的应变,进而通过模态分解的方法,获得立管模型涡激振动的位移。在此基础上,研究了CF以及IL方向的响应频率、位移标准差的平均值和最大值等随流速的变化规律,并分析了立管模型上测点的运动轨迹及其影响因素。     

Prediction Model for Vortex-Induced Vibration of Circular Cylinder with Data of Forced Vibration

A model based on the data from forced vibration experiments is developed for predicting the vortex-induced vibrations (VIV) of elastically mounted circular cylinders in flow. The assumptions for free and forced vibration tests are explored briefly. Energy equilibrium is taken into account to set up the relationship between the dynamic response of self-excited oscillations and the force coefficients from forced vibration experiments. The gap between these two cases is bridged straightforwardly with careful treatment of key parameters. Given reduced mass m and material damping ratio ζ of an elastically mounted circular cylinder in flow, the response characteristics such as amplitude, frequency, lock-in range, added mass coefficient, cross-flow fluid force and the corresponding phase angle can be predicted all at once. Instances with different combination of reduced mass and material damping ratio are compared to investigate their effects on VIV. The hysteresis phenomenon can be interpreted reasonably. The predictions and the results from recent experiments carried out by Williamson's group are in rather good agreement.     

A Novel Wake Oscillator Model for Vortex-Induced Vibrations Prediction of A Cylinder Considering the Influence of Reynolds Number

It is well known that the Reynolds number has a significant effect on the vortex-induced vibrations (VIV) of cylinders. In this paper, a novel in-line (IL) and cross-flow (CF) coupling VIV prediction model for circular cylinders has been proposed, in which the influence of the Reynolds number was comprehensively considered. The Strouhal number linked with the vortex shedding frequency was calculated through a function of the Reynolds number. The coefficient of the mean drag force was fitted as a new piecewise function of the Reynolds number, and its amplification resulted from the CF VIV was also taken into account. The oscillating drag and lift forces were modelled with classical van der Pol wake oscillators and their empirical parameters were determined based on the lock-in boundaries and the peak-amplitude formulas. A new peak-amplitude formula for the IL VIV was developed under the resonance condition with respect to the mass-damping ratio and the Reynolds number. When compared with the results from the experiments and some other prediction models, the present model could give good estimations on the vibration amplitudes and frequencies of the VIV both for elastically-mounted rigid and long flexible cylinders. The present model considering the influence of the Reynolds number could generally provide better results than that neglecting the effect of the Reynolds number.     

Numerical Simulation on Flow Past Two Side-by-Side Inclined Circular Cylinders at Low Reynolds Number

A series of three-dimensional numerical simulations is carried out to investigate the effect of inclined angle on flow behavior behind two side-by-side inclined cylinders at low Reynolds number Re=100 and small spacing ratio T/D=1.5 (T is the center-to-center distance between two side-by-side cylinders, D is the diameter of cylinder). The instantaneous and time-averaged flow fields, force coefficients and Strouhal numbers are analyzed. Special attention is focused on the axial flow characteristics with variation of the inclined angle. The results show that the inclined angle has a significant effect on the gap flow behaviors behind two inclined cylinders. The vortex shedding behind two cylinders is suppressed with the increase of the inclined angle as well as the flip-flop gap flow. Moreover, the mean drag coefficient, root-mean-square lift coefficient and Strouhal numbers decrease monotonously with the increase of the inclined angle, which follows the independent principle at small inclined angles.     

Experimental Study on Coupled Cross-Flow and in-Line Vortex-Induced Vibration of Flexible Risers

In this work, we study the coupled cross-flow and in-line vortex-induced vibration （VIV） of a fixedly mounted flexible pipe, which is free to move in cross-flow （ Y- ） and in-line （ X- ） direction in a fluid flow where the mass and natural frequencies are precisely the same in both X- and Y-direction. The fluid speed varies from low to high with the corresponding vortex shedding frequency varying from below the first natural frequency to above the second natural frequency of the flexible pipe. Particular emphasis was placed on the investigation of the relationship between in-line and cross-flow vibration. The experimental results analyzed by using these measurements exhibits several valuable features.     

Numerical Simulation of Self-Excited and Forced Vibration of Circular Cylinders in Current

Numerical simulations of a low-mass-damping circular cylinder which can oscillate freely at transverse and streamwise directions are presented in this work.The Navier-Stokes equations are solved with finite volume method,and large eddy simulation of vortex is also performed in the calculation.In order to implement dynamic mesh,overlapping grids are generated to lessen the computation for mesh field itself.Self-excited vibrations are firstly calculated to obtain the average amplitudes and frequencies of the target circular cylinder in the current flow situation,and then forced oscillations are implemented with parameters obtained in vortex-induced vibrations previously.With slight amplitude modulation,time series of displacements in vortex-induced vibrations are essentially harmonic.Regarding the fluid force,which are larger in forced oscillations than those in corresponding self-excited cases because the fluid subtracts energy from the forced cylinders.The phase angles between forces and displacements are 0° and 180° for self-excited case and forced case respectively.In vortex-induced vibrations,the interactions between fluid and structure produce some weakly energetic vortices which induce the modulations of amplitude and frequency.     

The Effect of Internal Fluid on the Response of Vortex-Induced Vibration of Marine Risers

Based on Iwan‘s wake oseillator model developed with the classical van der Pol equalion, the differential equation for the response of the vortex-induced vihration (VIV) of the riser considering the effect of the internal flowing fluid and the external marine environmental condition is derived. The effect of the intermal flowing fluid on the response of VIV of the riser is studied by means of the Finite Element Method. The results show that the effect of the internal fluid velocity on the VIV of the riser is strong when the natural frequency of the riser is close to the vortex shedding frequency. In addition, the increase of the top tension can decrease the sensitivity of the riser to the internal fluid velocity.     

Stress Wave Propagation Analysis on Vortex-Induced Vibration of Marine Risers

To analyze the stress wave propagation associated with the vortex-induced vibration (VIV) of a marine riser, this paper employed a multi-signal complex exponential method. This method is an extension of the classical Prony''s method which decomposes a complicated signal into a number of complex exponential components. Because the proposed method processes multiple signals simultaneously, it can estimate the "global" dominating frequencies (poles) shared by those signals. The complex amplitude (residues) corresponding to the estimated frequencies for those signals is also obtained in the process. As the signals were collected at different locations along the axial direction of a marine riser, the phenomena of the stress wave propagation could be analyzed through the obtained residues of those signals. The Norwegian Deepwater Program (NDP) high mode test data were utilized in the numerical studies, including data sets in both the in-line (IL) and cross-flow (CF) directions. It was found that the most dominant component in the IL direction has its stress wave propagation along the riser being dominated by a standing wave, while that in the CF direction dominated by a traveling wave.     

陡波形立管涡激振动疲劳损伤参数敏感性分析

基于柔性杆理论和尾流振子模型计算陡波形立管的涡激振动响应;综合使用S-N曲线法、雨流计数法、Palmgren-Miner线性累积疲劳理论对立管涡激振动导致的疲劳损伤进行计算分析.并以MATLAB为平台编写相应计算程序,将本文计算得到的静力分析结果、固有频率和疲劳损伤分别与专业海工计算软件O rcaFlex和已发表文献进...     

Tension and Drag Forces of Flexible Risers Undergoing Vortex-Induced Vibration

This paper presents the results of an experimental investigation on the variation in the tension and the distribution of drag force coefficients along flexible risers under vortex-induced vibration (VIV) in a uniform flow for Reynolds numbers (Re) up to 2.2×10⁵. The results show that the mean tension is proportional to the square of the incoming current speed, and the tension coefficient of a flexible riser undergoing VIV can be up to 12. The mean drag force is uniformly and symmetrically distributed along the axes of the risers undergoing VIV. The corresponding drag coefficient can vary between 1.6 and 2.4 but is not a constant value of 1.2, as it is for a fixed cylinder in the absence of VIV. These experimental results are used to develop a new empirical prediction model to estimate the drag force coefficient for flexible risers undergoing VIV for Reynolds number on the order of 10⁵, which accounts for the effects of the incoming current speed, the VIV dominant modal number and the frequency.