Vortex-induced vibration of two parallel risers: Experimental test and numerical simulation

The vortex-induced vibration of two identical rigidly mounted risers in a parallel arrangement was studied using Ansys-CFX and model tests.The vortex shedding and force were recorded to determine the effect of spacing on the two-degree-of-freedom oscillation of the risers.CFX was used to study the single riser and two parallel risers in 2–8D spacing considering the coupling effect.Because of the limited width of water channel,only three different riser spacings,2D,3D,and 4D,were tested to validate the characteristics of the two parallel risers by comparing to the numerical simulation.The results indicate that the lift force changes significantly with the increase in spacing,and in the case of 3D spacing,the lift force of the two parallel risers reaches the maximum.The vortex shedding of the risers in 3D spacing shows that a variable velocity field with the same frequency as the vortex shedding is generated in the overlapped area,thus equalizing the period of drag force to that of lift force.It can be concluded that the interaction between the two parallel risers is significant when the risers are brought to a small distance between them because the trajectory of riser changes from oval to curve 8 as the spacing is increased.The phase difference of lift force between the two risers is also different as the spacing changes.     

Combined action of uniform flow and oscillating flow around marine riser at low Keulegan-Carpenter number

With the increase of petroleum and gas production in deep ocean, marine risers of circular cylinder shape are widely used in the offshore oil and gas platform. In order to research the hydrodynamic performance of marine risers, the dynamic mesh technique and User-Defined Function （UDF） are used to simulate the circular cylinder motion. The motion of a transversely oscillat-ing circular cylinder in combination of uniform flow and oscillating flow is simulated. The uniform flow and oscillating flow both are in x direction. SIMPLE algorithm is used to solve the Navier-Stokes equations. The User-Defined Function is used to control the cylinder transverse vibration and the inlet flow. The lift and drag coefficient changing with time and the map of vorticity isolines at different phase angle are obtained. Force time histories are shown for uniform flow at Reynolds number （Re） of 200 and for the com-bination of uniform and oscillating flows. With the increase of amplitude of oscillating flow in combined flow, the change of lift am-plitude is not sensitive to the the change of cylinder oscillating frequency. Lift amplitude increases with the increase of oscillating flow amplitude in the combined flow, but there is no definite periodicity of the lift coefficient. The drag and inertia force coefficients change when the maximum velocity of the oscillating flow increases in the combined flow. The vortex shedding near the circular cylinder shows different characteristics.     

Investigation of a new vortex-induced vibration suppression device in laboratory experiments

In order to mitigate vortex-induced vibration (VIV) of marine risers, especially to eliminate the phenomenon of frequency ’lock-in’, a new suppression device of crescent-shaped flow spoiler was designed with seven different layout schemes. VIV model tests with six flow levels were conducted in a large wind-wave-current flume. In all cases, vibration responses in both in-line and cross-flow cases were measured. With the installation of suppression devices vibration frequency evolution of a riser was analyzed by Morlet wavelet transform. The principle of VIV suppression was interpreted through vibration characteristics. Fatigue life of the riser was calculated by the Palmgren-Miner rule. Compared with a bare riser, vibration of an outfitted riser with suppression devices disturbed the steady flow, the vibration amplitudes in the two flow directions were reduced, and the riser fatigue life was improved.     

Fatigue life assessment of top tensioned risers under vortex-induced vibrations

The fatigue life of top tensioned risers under vortex-induced vibrations (VIVs) with consideration of the effect of internal flowing fluid on the riser is analyzed in the time domain. The long-term stress histories of the riser under VIVs are calculated and the mean stresses, the number of stress cycles and amplitudes are determined by the rainflow counting method. The Palmgren-Miner rule for cumulative damage theory with a specified S-N curve is used to estimate the fatigue life of the riser. The corresponding numerical programs numerical simulation of vortex-induced vibrations (NSVIV) which can be used to calculate the VIV response and fatigue life of the riser are compiled. Finally the influences of the riser’s parameters such as flexural rigidity, top tension and internal flow velocity on the fatigue life of the riser are analyzed in detail and some conclusions are drawn.     

Synergetic Analysis and Possible Control of Vortex-Induced Vibrations in a Fluid-Conveying Steel Catenary Riser

This work aimed to demonstrate possibilities for both active and passive control of the vortex-induced vibration and fatigue life of steel catenary risers via an analysis of the self-organization and evolution of the structural vibration based on synergetic theory. An analysis of the complex interrelated and synergistic relationship between the order parameter and the fast variable was performed, and the master equation of the nodal displacements was established as the order parameter for the evolution of the riser’s structural vibration. Passive control methods include modifying the structure’s elastic modulus, the internal fluid velocity, the top tension and the structural damping ratio, while an active control involves adjusting the external flow rate. Optimized parameters were obtained by analyzing the non-steady state solution of the master equation. The results show that the fatigue life greatly increases as the riser’s elastic modulus decreases. In contrast, the fatigue life decreases with an increase of the internal fluid velocity. With an increase of the top tension, the vibration amplitudes and the number of modes may decrease, resulting in fewer bending stress cycles and a longer fatigue life. Furthermore, the structural damping ratio should be as large as possible. Finally, an active and passive control of the riser structure’s response to vortex-induced vibration and its fatigue life can be achieved by carefully modifying the parameters mentioned above. The results may provide a theoretical framework for engineering practice concerning the design and control of steel catenary riser structures which are affected by vortex-induced vibration.     

Parametric Stability Analysis of Marine Risers with Multiphase Internal Flows Considering Hydrate Phase Transitions

This study investigates the effects of multiphase internal flows that consider hydrate phase transitions on the parametric stability of marine risers.A numerical model of the multiphase internal flow that considers a hydrate phase transition is established.The model first solves the flow parameters and subsequently obtains the natural frequencies of risers with different gas intake ratios.The stability charts of marine risers with different gas intake ratios are plotted by applying Floquet theory,and the effects of the gas intake ratio on the instability and vibration response of the risers are identified.The natural frequency increases with an increase in the gas intake ratio;thus,instability zones move to higher frequency ranges in the stability charts.As the increasing gas intake ratio reduces the damping effect of the Coriolis force,the critical amplitude of the heave in the unstable region decreases,especially when hydrodynamic damping is not considered.As a result,higher-order unstable regions are excited.When in an unstable region,the vibration response curve of a riser with a high gas intake ratio excited by parametric resonance diverges quickly due to parametric resonance.     

Numerical investigation of vortex-induced vibrations for flow past a circular cylinder

The objective of the present investigation is to study the vortex-induced vibrations （VIV） for flow past a circular cylinder. The turbulent flow is simulated by using a 2-D standard k-ε model incorporating the finite volume method （FVM） and the Semi-Implicit Method for the Pressure Linked Equations （SIMPLE） algorithm on non-orthogonal boundary-fitted collocated grids. The wall boundaries are approximated with wall functions. In the numerical cases, the turbulent wake patterns are studied by plotting the streamlines and the turbulent kinetic energy contours. The pressure distributions are investigated. Analyses of the vortex-induced force coefficients and the structural vibrations are carried out. The variations of the Strouhal number with the Reynolds number and of the vortex-induced force coefficients with the reduced velocity are obtained. The results show that this numerical approach is feasible and efficient in investigating the VIV problem for a circular cylinder.     

3种附属管排布方式对主管路水动力特性的影响

【目的】研究均匀流多管束干涉流动下圆柱受迫振动的水动力特性。【方法】基于SSTκ-ω模型,在亚临界雷诺数下(Re=1×10⁵)对多管束共振强迫的涡激振动问题进行二维数值模拟,比较与分析三种典型附属管排布方式对主管路流体动力学特征的影响。【结果与结论】采用模型3下的附属管排布方式可在较大范围的振幅比下(A_y/D=0.1~0.8)有效改善主管路水动力特性:1)有效降低主管路上平均升力系数的幅值;2)抑制在单管路系统中出现的脉动升力系数突变衰减。同时,由于多管束对流动产生干涉效应,主管路上表现的尾迹涡度随着振动幅度的增大而呈现出不同的模式。此外,功率谱密度分析发现,多管束系统相较于单圆柱系统,在频率比为1时,模型2与模型3的共振“锁定”状态得到改善。     

Characteristics from a hydrodynamic model of a trapezoidal artificial reef

Flume experiments and numerical simulation were conducted to characterize the hydrodynamics of a trapezoid artificial reef. Measurements in particle image velocimetry were conducted to observe the formation of upwelling and vortices; and forces for the reef model were measured by load cell. The results of flume experiments agree well with the numerical data. In addition, the flow structure around a reef combining trapezoidal and cubic blocks was simulated numerically under two deployment schemes, showing a more complicated flow structure than that of a stand-alone reef. Relationship between drag coefficient and Reynolds number suggest that the degree of turbulence can be assessed from the value of drag coefficient downstream from the reef. The role of the reef in water flow is to reduce flow velocity and generate turbulence.     

Deep-water riser fatigue monitoring systems based on acoustic telemetry

Marine risers play a key role in the deep and ultra-deep water oil and gas production. The vortex-induced vibration (VIV) of marine risers constitutes an important problem in deep water oil exploration and production. VIV will result in high rates of structural failure of marine riser due to fatigue damage accumulation and diminishes the riser fatigue life. In-service monitoring or full scale testing is essential to improve our understanding of VIV response and enhance our ability to predict fatigue damage. One marine riser fatigue acoustic telemetry scheme is proposed and an engineering prototype machine has been developed to monitor deep and ultra-deep water risers’ fatigue and failure that can diminish the riser fatigue life and lead to economic losses and eco-catastrophe. Many breakthroughs and innovation have been achieved in the process of developing an engineering prototype machine. Sea trials were done on the 6th generation deep-water drilling platform HYSY-981 in the South China Sea. The inclination monitoring results show that the marine riser fatigue acoustic telemetry scheme is feasible and reliable and the engineering prototype machine meets the design criterion and can match the requirements of deep and ultra-deep water riser fatigue monitoring. The rich experience and field data gained in the sea trial which provide much technical support for optimization in the engineering prototype machine in the future.     

Preliminary Study on the Water Drag of Schlegel＇s Black Rockfish Sebastes schlegeli

Drag forces acting on Schlegel‘s black rockfish Sebastes schlegeli were studied. A new drag force transducer was designed and used to measure the water drag on Schlegel‘s black rockfish in a vertical recirculating flume tank. Fourteen individuals were investigated, yielding two mean drag coefficients referred to the cross-sectional area and volume^2/3 respectively at water velocities ranging from 0.3 to 1.0 ms^-1 The drag coefficients can be used for estimating the drag forces acting on Sebastes schlegeli in water.     

Damage localization of marine risers using time series of vibration signals

Based on dynamic response signals a damage detection algorithm is developed for marine risers. Damage detection methods based on numerous modal properties have encountered issues in the researches in offshore oil community. For example, significant increase in structure mass due to marine plant/animal growth and changes in modal properties by equipment noise are not the result of damage for riser structures. In an attempt to eliminate the need to determine modal parameters, a data-based method is developed. The implementation of the method requires that vibration data are first standardized to remove the influence of different loading conditions and the autoregressive moving average (ARMA) model is used to fit vibration response signals. In addition, a damage feature factor is introduced based on the autoregressive (AR) parameters. After that, the Euclidean distance between ARMA models is subtracted as a damage indicator for damage detection and localization and a top tensioned riser simulation model with different damage scenarios is analyzed using the proposed method with dynamic acceleration responses of a marine riser as sensor data. Finally, the influence of measured noise is analyzed. According to the damage localization results, the proposed method provides accurate damage locations of risers and is robust to overcome noise effect.     

CFD-based numerical analysis of a variable cross-section cylinder

Using ANSYS-CFX, a general purpose fluid dynamics program, the vortex-induced vibration (VIV) of a variable cross-section cylinder is simulated under uniform current with high Reynolds numbers. Large eddy simulation (LES) is conducted for studying the fluid-structure interaction. The vortex shedding in the wake, the motion trajectories of a cylinder, the variation of drag and lift forces on the cylinder are analyzed. The results show that the vortices of variable cross-section cylinder are chaotic and are varying along the cylinder. In places where cross-sections are changing significantly, the vortices are more irregular. The motion trail of the cylinder is almost the same but irregular. The drag and lift coefficients of the cylinder are varying with the changes of diameters.     

THE DYNAMIC EFFECTS OF A FINITE DEPTH AMBIENT FLUID ON BOTTOM EDDY

Scaling analysis shows that if o(ε~2, β_1ε, γε) ~o (δ), frontal geostrophic dynamics governs the behav-ior of an isolated bottom eddy in a finite depth ambient fluid; and that the ambient flow induced bybottom eddy migration satisfies quasi-geostrophic dynamics. This two layer model includes the impor-tant processes of advection of bottom eddy due to ambient flow, baroclinic instability, and form dragintroduced by Rossby waves. The numerical results show that the three processes enhance theinstability and alter the migration speed of the bottom eddy, and that the form drag induces asignificant meridional drift of the eddy.     

Numerical simulation of the transport and diffusion of dissolved pollutants in the changjiang (Yangtze) river estuary

Based on a coupled hydrodynamic-ecological model for regional and shelf seas (COHERENS), a three-dimensional baroclinic model for the Changjiang (Yangtze) River estuary and the adjacent sea area was established using the sigma-coordinate in the vertical direction and spherical coordinate in the horizontal direction. In the study, changing-grid technology and the “dry-wet” method were designed to deal with the moving boundary. The minimum water depth limit condition was introduced for numerical simulation stability and to avoid producing negative depths in the shallow water areas. Using the Eulerian transport approaches included in COHERENS for the advection and dispersion of dissolved pollutants, numerical simulation of dissolved pollutant transport and diffusion in the Changjiang River estuary were carried out. The mass centre track of dissolved pollutants released from outlets in the south branch of the Changjiang River estuary water course has the characteristic of reverse current motion in the inner water course and clockwise motion offshore. In the transition area, water transport is a combination of the two types of motion. In a sewage-discharge numerical experiment, it is found that there are mainly two kinds of pollution distribution forms: one is a single nuclear structure and the other is a double nuclear (dinuclear) structure in the turbid zone of the Changjiang River estuary. The rate of expansion of the dissolved pollutant distribution decreased gradually. The results of the numerical experiment indicate that the maximum turbid zone of the Changjiang River estuary is also the zone enriched with pollutants. Backward pollutant flow occurs in the north branch of the estuary, which is similar to the backward salt water flow, and the backward flow of pollutants released upstream is more obvious.     

A nonlinear vortex induced vibration model of marine risers

With the exploitation of oil and gas in deep water, the traditional vortex induced vibration (VIV) theory is challenged by the unprecedented flexibility of risers. A nonlinear time-dependent VIV model is developed in this paper based on a VIV lift force model and the Morison equation. Both the inline vibration induced by the flow due to vortex shedding and the fluid-structure interaction in the transverse direction are included in the model. One of the characteristics of the model is the response-dependent lift force with nonlinear damping, which is different from other VIV models. The calculations show that the model can well describe the VIV of deepwater risers with the results agreeing with those calculated by other models.     

Particle image velocimetry and numerical simulations of the hydrodynamic characteristics of an artificial reef

This article reports a particle image velocimetry study and the comparative results of a numerical simulation into the hydrodynamic characteristics around an artificial reef.We reveal the process of flow separation and vortex evolution,and compare the force terms generated by our artificial reef model.The numerical simulation agrees well with experimental results,showing the applicability of computational fluid dynamics to the hydrodynamics of an artificial reef.Furthermore,we numerically simulate the hydrodynamics of the reef model for seven velocities.The results show that the drag coefficient is approximately 1.21 in a self-modeling region for Reynolds numbers between 2.123×104and 9×104.Therefore,the upwelling height and current width of the flow field do not change significantly when the inflow velocity increases.Our study indicates that computational fluid dynamics can be applied to study the hydrodynamics of an artificial reef and offer clues to its construction.     

Two-dimensional numerical simulation of flow around three-stranded rope

Three-stranded rope is widely used in fishing gear and mooring system. Results of numerical simulation are presented for flow around a three-stranded rope in uniform flow. The simulation was carried out to study the hydrodynamic characteristics of pressure and velocity fields of steady incompressible laminar and turbulent wakes behind a three-stranded rope. A three-cylinder configuration and single circular cylinder configuration are used to model the three-stranded rope in the two-dimensional simulation. The governing equations, Navier-Stokes equations, are solved by using two-dimensional finite volume method. The turbulence flow is simulated using Standard κ-ε model and Shear-Stress Transport κ-ω(SST) model. The drag of the three-cylinder model and single cylinder model is calculated for different Reynolds numbers by using control volume analysis method. The pressure coefficient is also calculated for the turbulent model and laminar model based on the control surface method. From the comparison of the drag coefficient and the pressure of the single cylinder and three-cylinder models, it is found that the drag coefficients of the three-cylinder model are generally 1.3–1.5 times those of the single circular cylinder for different Reynolds numbers. Comparing the numerical results with water tank test data, the results of the three-cylinder model are closer to the experiment results than the single cylinder model results.     

A STUDY ON THE RELATIONSHIPS OF THE NUTRIENTS NEAR THE CHANGJIANG RIVER ESTUARY WITH THE FLOW OF THE CHANGJIANG RIVER WATER

Investigations from August, 1985 to July , 1986 showed that the high concentration area of PO4-P , SiO3-Si and NO3-N gradually reduced with the reduction of the area of the Changjiang River diluted water from summer, autumn to winter , and that the seasonal distributions and variations of the nutrients concentrations were mainly controlled by the river flow and were also related to the growth and decline of phytoplankton . The conservation of SiO3-Si and NO3-N in the estuary in the flood season was poorer than that in the dry season .. The behaviour of PO4-P in the estuary shows that aside from -biological removal, buffering of PCU-P is possible in the estuary . The highest monthly average concentrations and annual average concentrations in the river mouth were respectively 0.88 and 0.57 umol/L for PO4-P,191.5 and 96.2 umol/L for SiO3-Si, and 81.6 and 58.6 umol/L for NOs-N . The Changjiang's annual transports of PO4-P , SiO3-Si and NO3-N to the sea were about 1.4×104tons , 204.4×104 tons and 63.6×104     

Distribution of borehole temperature at four high-altitude alpine glaciers in Central Asia

The distribution of borehole temperature at four high-altitude alpine glaciers was investigated. The result shows that the temperature ranges from -13.4℃ to -1.84℃, indicating the glaciers are cold throughout the boreholes. The negative gradient （i.e., the temperature decreasing with the increasing of depth） due to the advection of ice and climate warming, and the negative gradient moving downwards relates to climate warming, are probably responsible for the observed minimum temperature moving to lower depth in boreholes of the Gyabrag glacier and Miaoergou glacier compared to the previously investigated continental ice core borehole temperature in West China. The borehole temperature at 10m depth ranges from -8.0℃ in the Gyabrag glacier in the central Himalayas to -12.9℃ in the Tsabagarav glacier in the Altai range. The borehole temperature at 10 m depth is 3-4 degrees higher than the calculated mean annual air temperature on the surface of the glaciers and the higher 10 m depth temperature is mainly caused by the production of latent heat due to melt-water percolation and refreezing. The basal temperature is far below the melting point, indicating that the glaciers are frozen to bedrock. The very low temperature gradients near the bedrock suggest that the influence of geothermal flux and ice flow on basal temperature is very weak. The low temperature and small velocity of ice flow of glaciers are beneficial for preservation of the chemical and isotopic information in ice cores.