High aspect ratio (L/D) riser VIV prediction using wake oscillator model

A two-dimensional (2-D) vortex-induced vibration (VIV) prediction model for high aspect ratio (L/D) riser subjected to uniform and sheared flow is studied in this paper. The nonlinear structure equations are considered. The near wake dynamics describing the fluctuating nature of vortex shedding is modeled using classical van der Pol equation. A new approach was applied to calibrate the empirical parameters in the wake oscillator model. Compared the predicted results with the experimental data and computational fluid dynamic (CFD) results. Good agreements are observed. It can be concluded that the present model can be used as simple computational tool in predicting some aspects of VIV of long flexible structures.     

VIV response of a long flexible riser fitted with strakes in uniform and linearly sheared currents

An experimental investigation was conducted on a flexible riser with and without various strake arrangements. The aim of the present work was to further improve the understanding of the response performance of the vortex-induced vibration (VIV) of a riser with helical strakes. Two current profiles, including uniform and linearly sheared flows, were simulated. The uniform current was simulated by towing the riser model in one direction using the towing carriage, and the linearly sheared current was simulated by fixing one end of the riser and using a driven cantilever to traverse a circular arc. Based on the modal superposition method, the displacement responses were obtained from the measured strain. Strakes with different heights and pitches were analysed, and response parameters such as the displacement response and fatigue damage were studied. The results of the bare model test show that the lock-in phenomenon displays multi-order characteristics, and the VIV displacement decreases with an increased order of the lock-in regime. The results of the straked model test indicate that the response characteristics of a bare riser can be quite distinct from those of a riser with helical strakes, and the response performance depends closely on the geometry of the strake configuration.     

抑制涡激振动的螺旋列板设计参数研究

基于水池模型实验结果和工程设计经验,结合国内外试验数据,着重分析用于抑制海洋立管涡激振动的螺旋列板几何参数(鳍高和螺距)及覆盖率对立管涡激振动的影响;并对水动力直径和水动力系数的选取对预报涡激振动的影响进行了分析,进而提出了适合于海洋立管工程应用的螺旋列板几何和设计参数选取的建议,为螺旋列板工程应用、海洋立管强度和疲劳设计提供参考。     

A prediction technique for the transverse vortex-induced oscillations of tensioned risers

In this paper, a model for the prediction of the vortex-induced transverse oscillations of tensioned risers is proposed.The model, based on the peak amplitudes of the oscillations recorded for cylindrical bodies in a uniform current, has been extended to multimodal oscillations in a sheared flow condition through a modal priority technique. The difference between this model and other similar ones lies in the following features:

the hydrodynamic damping, which has been introduced in a more general form without ‘a priori’ approximations;

the modal priority, which is based on the local amplitudes of the active vibration modes, rather than on their frequencies.

Some comparisons with experimental and numerical results available in the literature are presented and discussed.     

The effect of yaw angle on VIV suppression for an inclined flexible cylinder fitted with helical strakes

Experimental studies were carried out to investigate the response features of an inclined flexible bare cylinder as well as a straked cylinder in a towing tank, with the main purpose of further improving the understanding of the effect of yaw angle on vortex-induced vibration (VIV) suppression. Four yaw angles (a = 0°, 15°, 30°, 45°), which is defined as the angle between the cylinder axis and the plane orthogonal to the oncoming fluid flow, were tested. The cylinder model was towed along the tank to generate a uniform fluid flow. The towing velocity was in the range of 0.05–1.0 m/s with an interval of 0.05 m/s. The corresponding Reynolds number ranged from 800 to 16000. The strakes selected for the experiments had a pitch of 17.5D and a height of 0.25D, which is generally considered as the most effective configuration for VIV suppression of a flexible cylinder in water. The experimental results indicate that VIV suppression effectiveness of the inclined flexible straked cylinder is closely related to the yaw angle. The displacement amplitudes are significantly suppressed in both cross-flow (CF) and in-line (IL) directions at a = 0°. However, with increasing yaw angle, the suppression efficiencies of the CF and IL displacement amplitudes gradually decrease. In addition, the CF dominant frequencies of the straked cylinder obviously deviate from those of the bare cylinder at a = 0° and 15°. This deviation is substantially alleviated with increasing yaw angle. The IL dominant frequencies show less dependency on the yaw angle. Similar trends are also observed on the dominant modes of vibration and the mean drag coefficients.     

Analysis and prediction of vortex-induced vibrations of variable-tension vertical risers in linearly sheared currents

Many studies have tackled the problem of vortex-induced vibrations (VIV) of a vertical riser with a constant tension and placed in uniform currents. In this study, attention is focused on the cross-flow VIV modelling, time-domain analysis and prediction of variable-tension vertical risers in linearly sheared currents. The partial-differential equation governing the riser transverse motion is based on a flexural tensioned-beam model with typical pinned-pinned supports. The hydrodynamic excitation model describing the modulation of lift force is based on a distributed van der Pol wake oscillator whose nonlinear equation is also partial-differential due to the implementation of a diffusion term. The variation of empirical wake coefficients with system parameters and the water depth-dependent Reynolds number is introduced. Based on the assumed Fourier mode shape functions obtained by accounting for the effect of non-uniform tension, the Galerkin technique is utilized to construct a low-dimensional multi-mode model governing the coupled fluid-riser interaction system due to VIV. Numerical simulations in the case of varying sheared flow profiles are carried out to systematically evaluate riser nonlinear dynamics and highlight the influence of fluid-structure parameters along with associated VIV aspects. In particular, the effects of shear and tensioned-beam (tension versus bending) parameters are underlined. Some comparisons with published experimental results and observations are qualitatively and quantitatively discussed. Overall parametric analysis and prediction results may be worthwhile for being a new benchmark against future experimental testing and/or numerical results predicted by an alternative model and methodology.     

深水顶张式立管全耦合动力响应及参数敏感性分析

通过建立双套管与油管的全耦合运动控制方程,使用Newmark-β法,运用Matlab软件进行编程,对模型进行非线性时域动力耦合响应分析,分析在一定波浪环境,不同扶正器布置间距下,立管各层相同节点处耦合动力响应,研究各层管之间发生碰撞和摩擦的最大扶正器布置间距;分析相同扶正器布置间距下,不同顶张力系数对耦合动力响应的影响。结果表明:在波浪作用下,全耦合动力响应分析比等效管更能展现各层管的实际动力响应;扶正器布置间距与顶张力系数的取值不同,会影响各层立管的耦合响应与相对位移,取值不当时,立管会发生碰撞。     

深水钻井隔水管避台撤离动力与长度优化

钻井平台作业过程中遭遇强台风时,可能由于种种原因导致隔水管未被完全回收,平台只能下挂隔水管实施撤离。撤离过程中隔水管受到强烈的海流载荷作用,使平台的撤离航向与航速受到严重限制。提出隔水管避台撤离分析方法,建立软硬两种悬挂模式下撤离隔水管有限元模型,进行不同航向与航速下的隔水管悬挂撤离分析,研究隔水管悬挂撤离作业窗口,并对隔水管悬挂长度进行优化。分析表明,隔水管软悬挂撤离作业窗口较大,推荐平台采用软悬挂模式进行撤离,如不具备软悬挂实施条件,提前回收部分隔水管并将其余部分硬悬挂在平台上进行撤离也是一种可行的方案。     

A new approach for predicting the extreme response of structures

In this article a method to estimate the mean upcrossing intensity, μ(u), of a stochastic process is proposed. It is assumed that the stochastic process is a sum of a Gaussian process and a second-order correction term. The method is based on the two-dimensional Saddlepoint approximation. The accuracy of the method is tested on processes having analytical solutions for μ(u). Numerical examples are given where the stochastic process represents (i) the horizontal response of a floating offshore structure in a Gaussian sea, and (ii) the response of a structure subjected to a Gaussian wind velocity process. In addition, the estimates are compared to empirical upcrossing intensities of simulated responses. For case (ii), the obtained μ(u) estimates are compared to estimates obtained by numerical integration.     

错流超滤技术在海水胶体态铀、钍、镭同位素和有机碳研究中的应用

建立了由预过滤装置、蠕动泵、中空纤维超滤膜(AmiconH10P10-20,标称截留分子量10KDa)和连接管组成的错流超滤系统,利用荧光标记的40KDa葡聚糖和已知放射性活度的234Th示踪剂评估了超滤膜的截留和吸附性质,探讨了234Th在超滤过程中的渗透行为,考查了该系统用于实际海水样品时铀、钍、镭同位素和有机碳的质量平衡状况.结果表明,10Kda中空纤维超滤膜对40Kda葡聚糖具有良好的截留效率(85%),而吸附损失率为18%.铀、钍、镭同位素和有机碳在超滤过程中均达到极佳的质量平衡,回收率R=95%～98%,优于大多数文献报道的值.234Th在超滤过程中的渗透行为可以很好地用渗透模型加以描述.研究组分胶体态含量占“溶解”态含量的份额大小顺序如下:钍同位素、有机碳、镭同位素约等于铀同位素,这与钍为强颗粒活性元素、铀和镭为水溶性元素的地球化学性质相吻合.     

Experiments with flexible shrouds to reduce the vortex-induced vibration of a cylinder with low mass and damping

Experiments employing a low-mass-damping cylinder have been conducted to determine the vortex-induced vibration (VIV) response of four suppressors of the flexible-shroud family. The VIV suppressors were inspired in the concept of the Ventilated Trousers (VT), a flexible shroud composed of a flexible net fitted with three-dimensional bobbins. Reynolds number varied between 5 × 10³ and 25 × 10³, while reduced velocity varied from 2 to 26. The VIV dynamic response showed that the VT suppressed the peak amplitude of vibration down to 40% of that of a bare cylinder. Other flexible shrouds also achieved suppression, but not as efficiently. Drag was reduced during the VIV synchronization range, but remained above the value for a bare static cylinder thereafter. Spectral analysis of displacement and lift revealed that, depending on the geometry and distribution of the bobbins, the flexible shroud can develop an unstable behavior, capturing energy from the wake and sustaining vibrations for higher reduced velocities. PIV measurements of the wake revealed that the entrainment flow through the mesh is necessary to extend the vortex-formation length of the wake; this mechanism only occurs for the VT mesh.     

Physical Modeling and Parametric Study on Two-Degree-of-Freedom VIV of A Cylinder near Rigid Wall

Unlike most previous studies on the transverse vortex-induced vibration(VIV) of a cylinder mainly under the wallfree condition (Williamson & Govardhan,2004),this paper experimentally investigates the vortex-induced vibration of a cylinder with two degrees of freedom near a rigid wall exposed to steady flow.The amplitude and frequency responses of the cylinder are discussed.The lee wake flow patterns of the cylinder undergoing VIV were visualized by employing the hydrogen bubble technique.The effects of the gap-to-diameter ratio (e0/D) and the mass ratio on the vibration amplitude and frequency are analyzed.Comparisons of VIV response of the cylinder are made between one degree (only transverse) and two degrees of freedom (streamwise and transverse) and those between the present study and previous ones.The experimental observation indicates that there are two types of streamwise vibration,i.e.the first streamwise vibration (FSV) with small amplitude and the second streamwise vibration (SSV) which coexists with transverse vibration.The vortex shedding pattem for the FSV is approximately symmetric and that for the SSV is alternate.The first streamwise vibration tends to disappear with the decrease of e0/D.For the case of large gap-to-diameter ratios (e.g.e0/D = 0.54～1.58),the maximum amplitudes of the second streamwise vibration and transverse one increase with the increasing gapto-diameter ratio.But for the case of small gap-to-diameter ratios (e.g.e0/D = 0.16,0.23),the vibration amplitude of the cylinder increases slowly at the initial stage (i.e.at small reduced velocity V,),and across the maximum amplitude it decreases quickly at the last stage (i.e.at large Vr).Within the range ofthe examined small mass ratio (m＜4),both streamwise and transverse vibration amplitude of the cylinder decrease with the increase of mass ratio for the fixed value of V,.The vibration range (in terms of Vr ) tends to widen with the decrease of the mass ratio.In the second streamwise vibration region,the vibration frequency of the cylinder with a small mass ratio (e.g.mx = 1.44) undergoes a jump at a certain Vr,.The maximum amplitudes of the transverse vibration for two-degree-of-freedom case is larger than that for one-degree-of-freedom case,but the transverse vibration frequency of the cylinder with two degrees of freedom is lower than that with one degree of freedom (transverse).     

A methodology for structural analysis and optimization of riser connection joints

In riser structural systems, the region where the riser is connected to the platform is critical due to the development of high stress levels. To reduce the stress concentration in this region, bend stiffeners and stress joints are used in order to provide a gradual stiffness transition between the riser and the platform. The present paper presents an optimization design approach for bend stiffeners and stress joints, using a slender beam procedure for the connection joint analysis and an evolution strategies multi-objective optimization algorithm. This approach produces similar results when compared to a finite element analysis of the complete riser, but with significant reduction of computational costs. Numerical examples are presented, showing the efficiency and robustness of the suggested methodology.     

A wavelet for predicting the time-domain response of vertically tethered systems

The response wavelets for tether tension and the displacement of a package on the end of a tether are derived from a linear, frequency-domain model of vertically tethered systems using inverse fourier transform techniques. These wavelets are convolved with a record of ship displacement to predict the instantaneous tension and displacement to an accuracy of better than 20% for a typical system. Although snap loads are non-linear and the response cannot be calculated while the tether is slack, the onset of zero tension — a precursor to snap loads — is predicted by the wavelet. Thus, extreme and potentially hazardous conditions are foreseeable with the wavelet and this information can be used in real-time to guide operations at sea.     

Bounded and unbounded coupled transverse response of parametrically excited vertical marine risers and tensioned cable legs for marine applications

The paper deals with the non-linear dynamic response in the transverse direction of vertical marine risers or a tensioned cable legs subjected to parametric excitation at the top of the structure. The dynamic model contains both elastic and bending effects. The analytical approach reveals that the dynamic lateral response is governed by effects originated from the coupling of modes in transverse direction. The mathematical model is being treated numerically by retaining a sufficient number of transverse modes. Numerical results are given for specific case studies and refer both to the time histories of the lateral response for all modes of motion, and to the corresponding power spectral densities obtained through FFT. The numerical predictions are suitably plotted and discussed. The calculations concern both the undamped and the damped dynamic system. The damping in the system is a non-linear Morison type term, which describes the effect of the hydrodynamic drag. Both coupled and uncoupled equations are treated and points as well as regions of coupled and uncoupled stability and instability are defined. It is shown that the impacts originated from the coupling, evaluate new instabilities for the respective undamped system. The numerical results obtained through FFT of the time histories, provide qualitative conclusions for the features of the dynamic response for the modes of motions considered. Special attention has been paid to the effect of the hydrodynamic drag for the parametric excitation frequencies that guide the dynamic system to lie within a region of coupled instability.     

海水中胶体有机碳的测定──高温燃烧法和紫外/过硫酸钾法的比较

利用切向超滤(cross-flow ultrafiltration)技术对海水中胶体有机碳(COC)和真溶解有机碳(UOC)进行了分离,并分别用高温燃烧法和紫外/过硫酸钾法对胶体有机碳和真溶解有机碳进行了测定。测定结果表明,切向超滤能定量分离海水中的胶体有机碳,高温燃烧法和紫外/过硫酸钾法对胶体有机碳和真溶解有机碳的测定结果没有明显的系统误差,表明这两种方法的氧化效率基本相同,胶体粒子的存在,对紫外/过硫酸钾法的氧化效率没有显着影响。     

A finite differences formulation for the linear and nonlinear dynamics of 2D catenary risers

A finite differences (FD) solution method is proposed for the numerical treatment of the dynamic equilibrium problem of 2D catenary risers. The method is based on the so-called Box approximation, which in the scope of the present contribution is applied to the complete nonlinear model as well as to the reduced linearized formulation. The application of the Box method transforms the original governing systems into convenient sets of algebraic equations, which in turn are solved efficiently by the relaxation method. Extensive numerical calculations are presented that describe the dynamic behaviour of the structure and evaluate the amplification in loading due to the dynamic components. The effect of the geometric nonlinearities is assessed through comparative calculations that concern both mathematical formulations examined in the present, i.e. the complete nonlinear, and the reduced linearized model. Special attention is paid to the heave excitations as they amplify significantly the magnitudes of the loading components.     

A simple method to determine breaker height and depth for different deepwater wave height/length ratios and sea floor slopes

Existing, easily applicable methods to calculate the depth and height of breaking waves are hampered by two obstacles. First, the breaker depth is usually required to compute its height, and vice versa. Second, the equations take into account either the deepwater height to wavelength ratio or the sea floor slope, but not both. A simple iterative procedure is therefore proposed which incorporates both elements. For fully developed waves breaking over a nearly horizontal bottom, the breaker height and depth are also direct functions of the deepwater wavelength.     

A spectral response approach for detecting dominant phytoplankton size class from satellite remote sensing

An important goal in ocean colour remote sensing is to accurately detect different phytoplank- ton groups with the potential uses including the validation of multi-phytoplankton carbon cycle models; synoptically monitoring the health of our oceans, and improving our understanding of the bio-geochemical interactions between phytoplankton and their environment. In this paper a new algorithm is developed for detecting three dominant phytoplankton size classes based on distinct differences in their optical signatures. The technique is validated against an independent cou- pled satellite reflectance and in situ pigment dataset and run on the 10-year NASA Sea viewing Wide Field of view Sensor (SeaWiFS) data series. Results indicate that on average 3.6% of the global oceanic surface layer is dominated by microplankton, 18.0% by nanoplankton and 78.4% by picoplankton. Results, however, are seen to vary depending on season and ocean basin.     

A finite element approach for predicting the full resistance profile of a spudcan deeply penetrating in dense sand overlying clay

This paper presents a finite element approach to calculate the full resistance profile of a spudcan deeply penetrating in dense sand overlying clay, in which a potential for an installing spudcan to experience a sudden uncontrolled punch-through failure exists. A modified Mohr-Coulomb model characterized by incorporating a four-phase variation of the mobilized strength and dilation parameters with an equivalent accumulated plastic strain is developed and tested for the overlying dense silica sand. An extended Tresca model is used for the strain softening of the underlying clay. A series of large deformation finite-element (LDFE) analyses are carried out, varying the strength and dilation parameters as well as the spudcan geometries. A fairly good performance of the present approach is verified by validating against groups of centrifuge tests data, allowing the numerical study to be extended parametrically. The four-phase variation of the mobilized strength and dilation parameters involved in the progressive failure of the upper dense sand is parametrically studied and extended to cover the range of sand relative densities that are of practical interest. Additionally, comparisons with the typical existing LDFE analyses using both simple and sophisticated constitutive models are carried out. It shows that the present approach performs fairly well to calculate the full resistance profile of a spudcan deeply penetration in both thin and thick dense sand overlying clay, especially the peak and post-peak resistance, within around 5% of the corresponding centrifuge tests results.