Experimental study and numerical simulation on vortex-induced vibration of flexible riser

In order to study the effect of internal flow on vortex-induced vibration of flexible riser, the experiment on the vortex-induced vibration of flexible riser transporting fluid in the current was conducted in the physical oceanography laboratory of Ocean University of China. Considering the internal flowing fluid and external marine environment, the dynamic response of the flexible riser was measured. The corresponding numerical simulation was performed using the wake oscillatory model considering the extensibility of the riser system. Both the experiment and the numerical simulation indicated that with the increase of internal flow speed, the response amplitude increases, while the response frequency decreases.     

Numerical and Physical Investigation on Vortex-Induced Vibrations of Marine Risers

1 .Introduction Risers ,in deep water and strong current environments ,are prone to vortex-induced vibrations(VIV) .Vortex-inducedforces may excite the risersintheir normal modes of transverse direction.Un-der the“lock-in”condition,large resonant oscill…     

Experiments with a steel catenary riser model in a towing tank

A model test with a steel catenary riser (SCR) was conducted in a towing tank. The main purpose of the experiment was to gain further understanding of the global SCR dynamic behavior with vortex-induced vibration (VIV). To this end, a large model scale factor (250) was considered. The mass and stiffness of the riser model are very low, and the Reynolds number in the experiment ranged from 400 to 600. The experimental results demonstrated the influence of traveling waves on the cross-flow response of the riser model, which was verified by several analytical methods. This study provided important indications that the sources of the traveling waves (power-in regions) seem to change position along the riser length over time. Finally, the experimental design, data processing procedure, experimental set-up, and results are described in detail.     

波浪作用下三角型人工鱼礁水动力特性数值模拟与实验验证

人工鱼礁在波浪作用下水动力特性的研究对于人工鱼礁的工程安全与设计具有非常重要的意义。基于有限体积法,采用推板造波形式,通过利用VOF方法求解波面的方式建立了模拟人工鱼礁与波浪相互作用的三维数值波浪水槽。基于该数值模型研究了波浪作用下三角型镂空人工鱼礁的受力情况,并与物理模型实验结果进行比较,结果显示模拟得到的人工鱼礁受力和波浪形态均与实验结果吻合良好。根据数值模拟结果,拟合出波浪作用下人工鱼礁的水动力系数。单体三角型镂空人工鱼礁的速度力系数Cd随着Kc数、Re数的增加呈现减小的趋势;惯性力系数Cm则随着Kc数、Re数的增加呈现波动趋势。数值模拟结果显示在一个波浪周期内,人工鱼礁周围产生了较强的上升流和回流,其内部产生了明显的涡旋结构。研究结果为人工鱼礁的设计优化提供了理论依据。     

多频参激—涡激联合作用下顶张力立管分岔分析

复杂海洋环境中,顶张力立管顶部平台受到不同分量叠加的波浪作用,导致其轴向张力与多个频率相关,系统产生多频参激与涡激的联合振动,其动力学行为变得更为复杂,尤其是分岔和混沌等非线性振动特性有待进一步分析。首先,基于欧拉—伯努利梁理论,引入范德波尔尾流振子,建立和推导受多频参激—涡激联合作用下立管的动力学模型及其运动微分方程;接着,利用伽辽金法得到离散后的常微分方程组,采用多尺度法得到系统共振响应的调谐方程;最后,通过数值算例探究系统在多源和多频载荷联合作用下的共振响应。结果表明：对于多频参数激励,其中一个参数激励幅值在立管振动中起主导作用,而另一个会引起分岔的位置、数量和性质发生改变;随质量—阻尼参数增大,系统共振响应幅值整体先增大后减小,各分岔点出现的位置整体先向泄涡频率较大处移动,随后向泄涡频率较小处移动;质量—阻尼参数的变化还可能诱使系统发生倍周期分岔和混沌现象,危害结构安全。此外,采用直接数值积分求解微分方程组,验证近似分析结果,两者吻合较好。     

张力腿平台涡激运动模型试验与数值分析

为揭示张力腿平台涡激运动响应规律,采用数值模拟与模型试验相结合的方法研究了全水深系泊张力腿平台涡激运动响应。根据张力腿平台主尺度参数,按照几何相似制作了水池试验模型,在满足运动相似和动力相似的条件下开展了均匀流、剖面流模型试验,并将模型试验结果与数值模拟结果进行了对比,验证了数值模拟与模型试验结果的一致性。分析结果表明在均匀流作用下张力腿平台涡激运动的锁定区间在5.5<U_r<8.5,来流角对涡激运动影响较大。剖面流作用下平台涡激运动规律与均匀流作用基本一致,但剖面流造流引起的能量分散,使平台在XY平面的运动轨迹变得不规律,系泊系统提供的回复力对涡激运动幅值有抑制作用,来流角和流速对张力腿或立管模态影响明显。论文得到结论对于张力腿平台的工程设计有一定借鉴意义。     

网衣对波浪传播影响的三维数值模拟

为研究网衣对波浪传播的影响,采用多孔介质模型模拟网衣,建立了模拟网衣与波浪相互作用的三维数值波浪水槽模型。基于该数值模型研究了波浪经过网衣作用后波高衰减的变化规律,并与物理模型试验结果进行比较。通过对比数值模拟和物理模型试验的结果,证明了该三维数值模型模拟网衣对波浪传播影响的可行性。     

深海立管相互碰撞的响应面法研究

针对深海立管在理想环境和外部水环境中发生碰撞的情况,从立管碰撞最大径向位移和应力峰值两方面展开研究,充分利用ABAQUS仿真得到的数据,基于数值模拟和回归分析,提出四因素三水平的立管碰撞响应面分析方法。采用Design-Expert软件的BOX设计,对立管的非线性碰撞过程中的碰撞应力峰值和最大径向位移等两个试验指标进行回归分析并分别建立二次响应面回归模型。通过将响应面分析所得最优化解与有限元计算结果的对比,最终验证了响应面分析方法及所得模型的准确性。     

四叉树网格下的椭圆型缓坡方程数值模型研究

波浪是近岸海域关键的水动力因素之一。考虑到近岸地形复杂、波浪演化显著的特点,建立了四叉树网格体系下的椭圆型缓坡方程数值模型,采用有限体积法对模型进行数值离散,应用GPBiCG(m, n)算法求解离散后的控制方程。模型中根据波浪波长布局计算网格,生成多层次四叉树网格,对复杂计算域有较好的适应性,并且在离散和方程求解中无需引入形函数、不产生复杂的交叉项,节约了存储空间和计算时间。将模型成功应用于物理模型实验及Acapulco海湾的波浪场数值模拟,结果表明该模型能够准确、高效地模拟近岸波浪场,可为近岸波浪场的模拟提供一定的理论和技术支持。     

Simulating offshore sand waves

Sand waves form a prominent regular pattern in the offshore seabeds of sandy shallow seas and pose a threat to a range of offshore activities. A two-dimensional vertical (2DV) flow and morphological simulation model describing the behaviour of these sand waves has been developed. The simulation model contains the 2DV shallow water equations, with a free water surface and a general bed load formula. The water movement is coupled to the sediment transport equation with a seabed evolution equation. The domain is non-periodic in both directions. The spatial discretisation is performed by a spectral method based on Chebyshev polynomials. A fully implicit method is chosen for the discretisation in time. Firstly, we validate the simulation model mathematically by reproducing the results obtained using a linear stability analysis for infinitely small sand waves. Hereby, we investigate a steady current situation induced by a wind stress applied at the sea surface. The bed forms we find have wavelengths in the order of hundreds of metres when the resistance at the seabed is relatively large. The results show that it is possible to model the initial evolution of sand waves with a numerical simulation model. This paper forms the necessary first step to investigate the intermediate term behaviour of sand waves.     

Numerical simulations using momentum source wave-maker applied to RANS equation model

For Navier-Stokes equation model using the VOF scheme, Lin and Liu (Lin, P. and Liu, P.L.-F. (1999). Internal wave-maker for Navier-Stokes equations models. J. Waterw. Port Coast. Ocean Eng., 125 (4), 207–215.) developed an internal wave-maker method for which a mass source function of the continuity equation was used to generate target wave trains. Using this internal wave-maker method, various numerical experiments have been conducted without any problems due to waves reflected by a wave-maker. In this study, an internal wave-maker method using a momentum source function was proposed. Various numerical simulations in two and three dimensions were performed using the momentum source wave-maker applied to the RANS equation model in a CFD code, FLUENT. To verify their applicability in 2 dimensions, the computational results obtained using the momentum source wave-maker in a channel of constant depth were compared with the results obtained by using the mass source wave-maker and with the analytical solutions. And the results of the present numerical simulations of hydraulic experiments, which represent nonlinear waves on a submerged shoal and breaking waves on a plane beach, were compared with measurements. The comparisons showed good agreements between them. To see their applicability into 3-dimensional cases, the present results in a basin of constant depth were compared with the analytical solutions, and they agreed well with each other. In addition, vertical variation of longshore current was presented by using the 3-dimensional simulation results.     

圆形数值波浪水池斜向和多向入射波浪模拟研究

为了研究斜向入射波浪,基于三维不可压缩两相流模型,开发了一套圆形数值波浪水池数值模型。在圆形波浪水池中,通过源项造波法成功生成了任意入射方向的波浪,并且利用人工摩擦项模拟阻尼区以数值耗散反射波浪。模型基于嵌入式多块网格体系,采用FVM法(finite volume method)离散Navier-Stokes方程,VOF法(volume of fluid)追踪自由水面。试验结果表明,斜向入射波浪的模拟结果与理论值基本一致,圆形波浪水池在模拟斜向入射波浪时,有效区域的面积较传统波浪水池显著增大,而且有效区域受波浪入射角度的影响也较小。同时,通过叠加多列斜向入射波浪,模拟出了多向交叉波列,并通过与理论结果对比,发现其具有较高的精度。     

Numerical analysis of the forces exerted on offshore structures by ship waves

Ship-generated waves can contribute to the fatigue of offshore structures. This paper presents a numerical model for evaluating the forces exerted on a nearby fixed structure by ship-generated waves. The ship waves were modeled using Michell's thin-ship theory (Wigley waves), and the forces were computed using a boundary element method in the time domain. The simulation was validated by comparing its results with those of frequency-domain methods reported in the literature. It was then applied to calculate the forces exerted on a hemisphere by ship waves varying with the ship's speed, dimensions and distance from the hemisphere to the ship's path. Our results indicate that the ship waves have enormous effects on offshore structures and are not neglectable.     

基于WRF-SWAN耦合模式的台风“威马逊”波浪场数值模拟

获取高分辨率的风场数据和气压场数据是精确模拟台风浪的基础,采用经验公式构建台风风场和气压场对海浪模式进行驱动,无法反映台风影响下海气动力过程,难以提供高精度的风场、气压场数据。本文基于中尺度大气模式WRF(Weather Research and Forecasting model)和第三代海浪模式SWAN(Simulating WAves Nearshore model),构建了南中国海地区大气—海浪实时双向耦合模式,针对超强台风"威马逊"进行数值模拟。将数值模拟结果与现场观测结果及卫星高度计观测结果进行对比验证,验证结果表明,本文建立的WRF-SWAN耦合模式在对台风"威马逊"影响下的南中国海台风浪的模拟中展现出较高的模拟精度,揭示了台风风场分布和台风浪分布在空间上的"右偏性"不对称分布特征及其形成机制。基于WRF和SWAN建立的大气-海浪实时双向耦合模式能够准确模拟台风动力过程以及台风浪的时空分布特征,可以推广用于南中国海地区台风浪的模拟分析。     

基于格子Boltzmann方法的斜坡堤越浪数值模拟研究

准确确定越浪量对于斜坡堤设计有重要意义。利用格子Boltzmann方法(LBM),并采用主动吸收式速度入口造波、出流边界消波、VOF方法追踪自由表面以及静态Smagorinsky模型模拟紊流运动,建立二维数值波浪水槽,对光滑斜坡堤上规则波与不规则波越浪进行数值模拟。模拟结果与试验值及其他数值模型结果比较表明,二维LBM数值波浪水槽具有模拟斜坡堤越浪的能力,但对于破碎较为剧烈的越浪过程模拟,该模型还存在一定的不足,未来可从提高自由表面模型精度等方面进一步改善其性能。     

钢悬链线立管与海床相互作用小尺寸三维试验研究

围绕钢悬链线立管(SCR)与海床的相互作用,在水箱内开展三维试验研究,研究在不同模拟运动激励下SCR触地点的应力状态。针对当前模拟试验中,全尺寸试验耗资巨大,且试验环境难以控制,缩尺试验大多模拟立管二维运动等现状,提出一套简单易行的三维管土作用试验装置,通过横向、纵向、垂向三个方向轨道位置的合理布置,使得立管可在单向、二维耦合和三个方向同时运动,对模型立管的顶端、底端的边界条件进行处理,通过驱动器在顶端施加位移,模拟在周期运动作用下,立管触地区与土的相互作用,在三维空间内研究立管的力学特性。由此指导立管的整体设计与分析,对保证SCR在深水油气开采中的安全可靠性,具有非常重要的意义。     

A nonlinear viscoelastic bend stiffener steady-state formulation

Bend stiffeners are essential components of a flexible riser system, employed to ensure a smooth transition at the upper connection and to protect the riser against over bending and from accumulation of fatigue damage. The highly nonlinear rate dependent behavior of these structures directly affects the integrity assessment of the riser in one of its most critical regions, the top connection. A steady-state formulation (disregarding inertial forces) and numerical solution procedure is developed in this work employing the perturbation method for a nonlinear viscoelastic bend stiffener large deflection beam model subjected to harmonic loading conditions. For stochastic loading conditions, the response is calculated employing the superposition principle by summing up the steady-state result of a number of individual frequency components. A time domain formulation is also derived employing the state-variable approach for the numerical solution of the resulting hereditary integral in the governing equations. A case study is presented for the top connection system of a 4″ ID flexible riser using relaxation and tensile experimental data obtained from a typical class of bend stiffener polyurethane. Harmonic and stochastic input loading conditions are employed for time and frequency domain model comparison/validation and to assess loading history and frequency influence in the curvature response.     

Transparent boundary condition for simulating nonlinear water waves by a particle method

Open boundaries are important when simulating water waves. In this study, a transparent boundary condition at an open boundary was developed for simulating nonlinear water waves propagating to a distant area using the Moving Particle Semi-implicit method. The novelty of this study is that the technique of wave analysis used in the experiment was introduced into the particle simulation to absorb incident waves; the simulation cost was reduced by employing inflow and outflow regions instead of a long dissipation region. Incident waves in front of the boundary were evaluated using Fourier analysis, and the particles on the transparent boundary were forced to move at the velocity of the analytical solution for Stokes waves in order to absorb the incident waves. The analysis was restricted to periodic waves. Wave propagation was simulated for two wave periods using the developed transparent boundary condition. The results showed that this transparent boundary transmitted the incident waves with small reflection and the simulation cost was lower than that for wave damping by a conventional highly viscous region.     

Numerical Investigation on Severe Slugging in A Catenary Riser

A mathematical model is presented to study the main characteristics of severe slugging in a downward inclined pipeline followed by a catenary riser. In this model, both simplified transient model and phase distribution model are included so that the flow characteristics of each stage for severe slugging can be accurately reproduced, especially for blowout stage. The results show that the flow features of severe slugging can be simulated by the proposed mathematical model. A good agreement between the experimental data and the numerical results is observed. The model can predict the transient fluctuation of many respects, such as the superficial gas velocity at the bottom of the riser and the average velocity at the outlet of the riser.