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

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

深海细长柔性立管涡激振动数值分析方法研究进展

深海细长柔性立管的涡激振动与刚性圆柱的涡激振动问题相比,既有联系又有区别.当立管的长细比(L/D)较大时,会展现出一些刚性圆柱涡激振动没有的现象.本文介绍了近几十年来深海开发的细长柔性立管涡激振动的研究进展,主要侧重于数值计算分析研究领域,对目前常用的数值分析方法进行了归纳和总结.在此基础之上,详细介绍了目前理论研究和实际工程中的热点问题,诸如流场剖面的影响、立管顶端浮体运动的耦合、涡激振动抑制措施、振动多模态特性等.     

深海采矿系统中悬臂式立管涡激振动分析

不同于传统的海洋立管, 深海采矿系统中的垂直提升管道可以被视为一个底部无约束的柔性悬臂式立管, 工作过程中同样面临涡激振动和柔性变形问题. 本文采用一种无网格离散涡方法和有限元耦合的准三维时域求解数值模型, 系统性地研究了不同流速下悬臂式立管的涡激振动问题. 结果表明: 悬臂式立管的横向振动模态阶数随折合速度增加而增大, 在一定折合速度范围内主导振动模态保持不变; 当主导模态转变时, 对应的横向振幅会发生突降, 但是当新的高阶模态被激发后, 立管振幅随来流速度增加而再次逐渐增大; 在相同的振动模态下, 立管底部位移均方根值随折合速度线性增加, 主导振动频率在模态转变时会出现跳跃现象; 特别地, 本文讨论了三阶主导模态下悬臂式立管的振动响应, 无约束的立管底部呈现出较大的振动能量, 且振幅的驻波特征随折合速度增加而逐渐增强; 本文比较了两端铰支立管与悬臂式立管的涡激振动响应特征, 两者在振幅和主导振动频率两方面均表现出了相同的变化趋势.      

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

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

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

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

海洋柔性结构涡激振动的流固耦合机理和响应

对近几十年来国内外在涡激振动的基础研究包括机理认识和动响应分析等方面的进展进行了论述,尤其针对海洋油气平台中的立管、隔水管等细长柔性结构的涡激振动.描述了涡激振动这种典型的非线性流固耦合现象所具有的特征,包括自激、自限制、展向相关、尾迹水动力与结构动力的流固耦合等及其主要影响参数.介绍了目前常用的结构响应预测方法和相关实验.通过讨论当前理论研究和实际工程中的热点问题,诸如多模态宽带振动、浮体运动与水下立管的耦合、响应抑制措施、双向振动、高雷诺数下的大尺度物理实验等,对今后该领域的研究方向进行了力所能及的展望.     

剪切流场中含内流立管横向涡激振动特性

立管是海洋工程中输送油气或其他矿产资源的必备结构, 外部洋流引起的立管涡激振动影响着立管的疲劳寿命, 危害深海资源开发. 本文基于欧拉?伯努利梁方程, 结合半经验时域水动力模型, 建立剪切流与内流耦合作用下海洋立管涡激振动预报模型, 运用有限元方法和Newmark-β逐步积分法求解方程, 首先将数值模拟结果与实验数据进行对比, 验证模型正确性. 然后, 运用此模型, 对剪切流作用下含内流的顶张立管在不同内流速度和密度下的横向涡激振动响应特性进行研究, 主要分析了立管的横向振动模态、振动频率以及均方根位移等涡激振动参数随内流速度和密度等参数的变化规律. 结果表明, 在剪切流场中, 含内流海洋立管在横向上表现出多模态多频率的涡激振动;立管横向振动的最大均方根位移随内流速度和密度的增大而增大, 特别是当内流速度较大时, 横向最大均方根位移增大明显;立管横向振动的主导频率随内流速度和密度的增大而减小, 并且内流密度的增大同样会引起模态转换和频率转换.      

海洋柔性结构涡激振动的流固耦合机理和响应

对近几十年来国内外在涡激振动的基础研究包括机理认识和动响应分析等方面的进展进行了论述,尤其针对海洋油气平台中的立管、隔水管等细长柔性结构的涡激振动.描述了涡激振动这种典型的非线性流固耦合现象所具有的特征,包括自激、自限制、展向相关、尾迹水动力与结构动力的流固耦合等及其主要影响参数.介绍了目前常用的结构响应预测方法和相关实验.通过讨论当前理论研究和实际工程中的热点问题,诸如多模态宽带振动、浮体运动与水下立管的耦合、响应抑制措施、双向振动、高雷诺数下的大尺度物理实验等,对今后该领域的研究方向进行了力所能及的展望.     

Fluid-solid coupling and dynamic response of vortex-induced vibration of slender ocean cylinders

对近几十年来国内外在涡激振动的基础研究包括机理认识和动响应分析等方面的进展进行了论述,尤其针对海洋油气平台中的立管、隔水管等细长柔性结构的涡激振动.描述了涡激振动这种典型的非线性流固耦合现象所具有的特征,包括自激、自限制、展向相关、尾迹水动力与结构动力的流固耦合等及其主要影响参数.介绍了目前常用的结构响应预测方法和相关实验.通过讨论当前理论研究和实际工程中的热点问题,诸如多模态宽带振动、浮体运动与水下立管的耦合、响应抑制措施、双向振动、高雷诺数下的大尺度物理实验等,对今后该领域的研究方向进行了力所能及的展望.     

两端铰接的细长柔性圆柱体涡激振动响应特性数值研究

目前细长柔性圆柱体涡激振动响应的研究方法主要包括实验方法、计算流体动力学方法以及半经验模型方法. 鉴于实验方法研究成本较高、计算流体动力学方法计算时间较长,本文基于尾流振子模型对线性剪切来流下两端铰接的细长柔性圆柱体涡激振动响应特性进行了半经验模型方法研究. 先建立了柔性圆柱体结构振子以及尾流振子之间的耦合模型,紧接着基于二阶精度中心差分格式对耦合模型先离散后迭代进行求解. 对不同剪切参数下柔性圆柱体涡激振动响应的振动波长、振动频率、振动位移以及响应频率随时间的变化特性等参数进行了分析. 分析结果表明:圆柱体的涡激振动响应由驻波和行波混合组成. 当无量纲弯曲刚度较小时,在圆柱体两端附近,驻波占主导;而在圆柱体中间段附近,行波占主导. 当无量纲弯曲刚度较大时,在圆柱体整个长度区间上均为驻波占主导. 随着剪切参数的增大,振动位移以及振动波长均逐渐减小,而振动频率和频率带宽均逐渐增大.      

Drag amplification of long flexible riser models undergoing multi-mode VIV in uniform currents

Measurement data of long flexible riser models are used to study the drag amplification due to multi-mode vortex induced vibration (VIV). The riser model was towed in a towing tank, and vibration along its length was measured indirectly at a set of discrete points using accelerometers, along with the total drag and lift forces on the model. These data were used to fit expressions for the drag amplification taking into account the spatial variation of the vibration amplitude along the riser length. The results were compared with the existing empirical expression used in the VIV analysis tool SHEAR7 as well as other test data. It was found that the expression in SHEAR7 agrees well with the test data. However, for larger vibration amplitudes, it appears to under-predict the drag amplification.     

DNS-DERIVED FORCE DISTRIBUTION ON FLEXIBLE CYLINDERS SUBJECT TO VORTEX-INDUCED VIBRATION

We use direct numerical simulation (DNS) based on spectral methods to simulate turbulent flow past rigid and flexible cylinders subject to vortex-induced vibrations (VIV). We present comparisons of amplitude, and lift and drag forces, at Reynolds number 1000 for a short and a long cylinder, and we examine differences between a traveling wave response and a standing wave response. The DNS data suggest that the often-used empirical formula proposed by Skop, Griffin & Ramberg in 1977 overpredicts the drag coefficient. We propose an appropriate modification and present preliminary results that indicate that low-dimensional modeling may be an accurate and efficient approach in predicting forces in VIV. Given the lack of any benchmark experiments in VIV currently, the DNS results presented here, both distributions as well as span- and time-averaged quantities, should be helpful to experimentalists and modelers.     

Optimal lift force coefficient databases from riser experiments

Significant past effort has gone into understanding the complicated flow–structure interaction problem of vortex-induced vibration (VIV) of long flexible cylindrical structures (e.g., risers, mooring lines, tendons, conductors) in the ocean environment. However, major challenges persist with regard to riser VIV modeling and response prediction. The existing prediction schemes are based on a number of hypotheses, experimental facts and data like strip theory, energy balance, correlation length and, most importantly, the use of lift force coefficient databases. Recent advances in observing the VIV motions on experimental risers with high confidence shows that some of these assumptions may not be valid. One important source of the discrepancies between theoretical estimates and experimental observations arise from the use of experimentally obtained lift coefficient databases. These databases were obtained under the laboratory conditions of limited Reynolds number, and under the assumption that the cross-flow motions are not influenced by restraining the in-line motions. In this paper we develop a method to improve the modeling capability of riser VIV by extracting empirical lift coefficient databases from field riser VIV measurements. The existing laboratory-based lift coefficient databases are represented in a flexible parameterized form using a set of carefully chosen parameters. Extraction of the lift coefficient parameters is posed as an optimization problem, where the objective is to minimize the error between the prediction using a theoretical model and the experimental data. Application of the method to data from the Norwegian Deepwater Programme experiments shows that the new optimal databases significantly reduce the error in estimating the riser VIV cross-flow response.     

带整流罩隔水管流场特性的大涡模拟研究

为了深入揭示整流罩抑制涡激振动的机理，采用CFD软件FLUENT，结合大涡模拟方法，分别对附加整流罩前后隔水管模型的绕流流场进行模拟，将附加整流罩前后隔水管周围的流场特性、升力系数、阻力系数及其频谱特性等进行对比分析．结果表明，整流罩能有效抑制隔水管尾流区的漩涡脱落，隔水管背部的压力升高，升力系数和阻力系数显著降低，从而抑制了隔水管的涡激振动．     

细长锥边界层绊线转捩风洞自由飞试验

通过在半锥角θ_c=10°细长锥面上布置一定数量的人工绊线,促使细长锥表面边界层在相应轴向位置上发生层流向湍流转变的固定转捩,采用运动自由度不受约束的风洞自由飞试验技术研究边界层转捩对高超声速细长锥再入体无控自由飞行下的运动特性和气动特性影响规律,并与以往无人工绊线的细长锥风洞自由飞试验结果作对比.试验马赫数Ma=5.0,通过改变风洞前室总压P_0实现两个雷诺数的模拟,以模型长为特征尺寸自由流雷诺数分别为0.84×10~6和1.68×10~6.结果表明:当自由流雷诺数Re=0.84×10~6时,人工绊线尚不足以促使边界层发生转捩,有绊线的细长锥气动特性与无绊线基本一致,动稳定导数大于零;当自由流雷诺数Re=1.68×10~6时,人工绊线促使边界层发生固定转捩,细长锥的动稳定导数小于零,细长锥自由飞行动稳定.     

亚临界雷诺数下圆柱受迫振动的数值研究

在Navier-Stokes方程和k-ω湍流模型的基础上,利用流线迎风有限元方法结合ALE动网格技术对亚临界雷诺数下的圆柱受迫振动问题开展了数值模拟研究。本文的数值模型成功模拟了Re=5000条件下,圆柱发生受迫振动时尾迹区内的2S,2P和P+S尾流模式;对Re=10000情况下,无量纲振幅分别为0.3,0.4,0.5的圆柱受迫振动问题开展了数值模拟,分析了给定振幅条件下圆柱受力随振动频率的变化关系以及受迫振动的锁定区间。以上数值计算结果与Gopalkrishnan (1993)的实验结果基本符合。研究结果表明,二维数值模型能够基本正确地反映出圆柱发生受迫振动时的涡激振动特性以及有关的受力变化趋势,为今后进一步开展三维数值分析工作奠定了基础。     

A harmonic balance technique for the reduced-order computation of vortex-induced vibration

We present a harmonic balance (HB) method to model frequency lock-in effect during vortex-induced vibration (VIV) of elastically mounted circular cylinder and a flexible riser section in a freestream uniform flow. The fluid flow and structure are coupled by a fixed-point iteration process through a frequency updating algorithm. By minimizing the structural residual in the standard least-square norm, the convergence of HB-based fixed-point algorithm is achieved for a range of reduced velocity. To begin with, the HB solver is first assessed for a periodic unsteady flow around a stationary circular cylinder. A freely vibrating circular cylinder is then adopted for the reduced-order computation of VIV at low Reynolds numbers of Re=100 and 180 with one- and two-degrees-of-freedom. The coupled VIV dynamics and the frequency lock-in phenomenon are accurately captured. The results show that the HB solver is able to predict the amplitude of vibration, frequency and forces comparable to its time domain counterpart, while providing a significant reduction with regard to overall computational cost. The proposed new scheme is then demonstrated for a fully-coupled three dimensional (3D) analysis of a linear-elastic riser section undergoing vortex-induced vibration in the lock-in range. The results reveal the 3D effects through isosurfaces of streamwise vorticity blobs distributed over the span of flexible riser section. In comparison to time domain results, the 3D flow-structure interactions are accurately predicted while providing a similar speed up rate that of 2D simulations. This further corroborates that the HB solver can be extended to 3D flow-structure dynamics without compromising efficiency and accuracy.     

Numerical study on the suppression of the vortex-induced vibration of an elastically mounted cylinder by a traveling wave wall

In the present paper, the commercial CFD code “Fluent” was employed to perform 2-D simulations of an entire process that included the flow around a fixed circular cylinder, the oscillating cylinder (vortex-induced vibration, VIV) and the oscillating cylinder subjected to shape control by a traveling wave wall (TWW) method. The study mainly focused on using the TWW control method to suppress the VIV of an elastically supported circular cylinder with two degrees of freedom at a low Reynolds number of 200. The cross flow (CF) and the inline flow (IL) displacements, the centroid motion trajectories and the lift and drag forces of the cylinder that changed with the frequency ratios were analyzed in detail. The results indicate that a series of small-scale vortices will be formed in the troughs of the traveling wave located on the rear part of the circular cylinder; these vortices can effectively control the flow separation from the cylinder surface, eliminate the oscillating wake and suppress the VIV of the cylinder. A TWW starting at the initial time or at some time halfway through the time interval can significantly suppress the CF and IL vibrations of the cylinder and can remarkably decrease the fluctuations of the lift coefficients and the average values of the drag coefficients; however, it will simultaneously dramatically increase the fluctuations of the drag coefficients.