Hydrodynamic forces on partly buried tandem, twin pipelines in current

This study extends the investigations of the forces on a cylinder, laid on, or partly buried in the bed with a parallel twin dummy cylinder nearby and without it, and were determined by measuring the pressure distribution on the cylinder in the case of a steady current. The pressure distribution around the cylinder was measured by using pressure transducers. The forces on the cylinder were calculated by the integration of the measured pressures on the surface of the cylinder. Force coefficients were obtained for the ranges of Re=0.8×10⁴–1.5×10⁴ for the burial depth to diameter RATIO=0:0.7. The distance between axis of the measurement and dummy cylinders to diameter ratio (x/D) was 2, 1.5 and 1. The dummy cylinder was replaced downstream and upstream of the measurement cylinder.     

Hydrodynamic forces on a partly buried cylinder exposed to combined waves and current

The present study extends the investigations of the hydrodynamic forces on a cylinder, laid on, or partly buried in the bed. They were determined by measuring the pressure distribution on the cylinder surface in the case of steady current, waves and coexisting flow. The pressure distribution around the cylinder was measured by using pressure transducers, which were replaced in the cylinder. Force coefficients were obtained for the ranges of Re=0.8×10⁴–1.5×10⁴, for steady current, low KC numbers (KC<5) for wave alone case and, for current-to-wave velocity RATIO=0, 3, 6 and infinity (current) for coexisting flow. The forces were also determined for the various burial-depth-to-diameter ratios between 0 and 0.7 values of the cylinder.     

不同倒角半径四柱体绕流数值模拟及水动力特性分析

为研究四柱体布置情况下倒角半径变化对柱体绕流水动力特性的影响,使用Fluent软件,采用大涡模拟方法研究了在雷诺数Re=3 900下6种不同倒角半径的柱体在方形四柱体布置时的三维流场。在模型分析验证有效后,分析了柱体后方瞬时流场、水动力参数、时均流场的变化情况。分析结果表明:随着倒角半径的增大,上游柱体的平均阻力系数逐渐减小,下游柱体的平均阻力系数除了在R~+=0.1处增幅很大以外,其余均随倒角半径变大而平稳变大;各柱体的升力系数均方根变化趋势基本相同;R~+=0.1、0.5时,上下游两柱体的升力系数曲线相位相反,而在R~+=0.2、0.3和0.4时,上下游两柱体的升力系数曲线相位相同。     

Flow induced vibrations of a sharp edge square cylinder in the wake of a circular cylinder

The response of an oscillating circular cylinder at the wake of an upstream fixed circular cylinder was classified by different researchers as galloping, wake induced galloping or wake induced vibration. Furthermore it is already known that a sharp edge square cylinder would undergo galloping if it is subjected to uniform flow. In this study the influence of the wake of a fixed circular cylinder on the response of a downstream square cylinder at different spacing ratios (S/D = 4, 8, 11) is experimentally investigated. The subject appears not to have received previous attention. The lateral displacements, lift forces and the pressure data from gauges mounted in the wake of the oscillating cylinder are recorded and analyzed. The single degree of freedom vibrating system has a low mass-damping parameter and the Reynolds number ranges from 7.7 × 10² to 3.7 × 10⁴.In contrast to that for two circular cylinders in tandem arrangement, the freely mounted downstream square cylinder displays a VIV type of response at all spacing ratios tested. There is no sign of galloping or wake induced galloping with the square cylinder. With increase at the spacing ratio the cross-flow oscillations decrease. It is shown that the vortices arriving from the upstream fixed circular cylinder play a major role on the shedding mechanism behind the downstream square cylinder and cause the square cylinder to shed vortices with frequencies above Strouhal frequency of the fixed square cylinder (St = 0.13). The VIV type of oscillations in the downstream square cylinder is most probably caused by the vortices newly generated behind the square cylinder.     

Three-dimensional flow around a square cylinder near a wall

In this study, two- and three-dimensional numerical simulations were performed to investigate the effect of the flow structure in the wake of a square cylinder placed near a plane wall by applying a fully implicit finite-difference method to the Navier-Stokes equations. The gap ratio between the cylinder and the wall, G/D, was varied from 0.2 to 4 for the Reynolds numbers of 175, 185 and 250. The role of the 3D structure on the lift and drag coefficients and Strouhal number was investigated. The results were compared with those of the 2D numerical simulations. The deviations of the 3D flow structure of the cylinder-wall pair from that of a single cylinder were also reported. At Re=250, B type secondary vortices were determined in the wake region. At Re=175 and 185, transition from A type vortex to fully periodic B type vortices was observed when the cylinder was brought closer to the wall.     

Numerical investigation on the suppression of VIV for a circular cylinder by three small control rods

Laminar flow past a circular cylinder with 3 small control rods is investigated by numerical simulation. This study is concerned with the suppression efficacy of vortex induced vibration by small control rods located around a main cylinder. The effects of the attack angle and rod-to-cylinder gap ratio on the hydrodynamics and vibration responses of the main cylinder are investigated. The attack angle of α = 45° is performed as the critical angle for VIV suppression of 3 control rods. The 3 control rods have no effect on VIV suppression when the attack angle is less than the critical angle. The 3 control rods have an excellent VIV suppression efficacy when the attack angle is larger than the critical angle. The transverse vibration frequency of the cylinder with 3 control rods is less than that for an isolated cylinder for all the configurations. The numerical results for the configurations of α = 45° & 60°, G/D = 0.6–1.2 show excellent suppression efficient among the cases investigated in this study. The best suppression efficient is found at α = 45°, G/D = 0.9 for 3 control rods. 2 rods in behind of the main cylinder perform more efficient than that of 1 rod in front for VIV suppression as the gap ratio of G/D less than 1.0.     

Wave forces on, and water-surface fluctuations around a vertical cylinder encircled by a perforated square caisson

The wave forces and moments on and the water surface fluctuations around a vertical circular cylinder encircled by a perforated square caisson were experimentally investigated. The porosity of the outer square caisson was varied from 4.24 to 14.58%. The in-line wave forces on the inner vertical cylinder are influenced by changing the porosity of the outer caisson, whereas the variations in the water surface fluctuations are less influenced in this porosity range. The in-line moment on the vertical cylinder is relatively less sensitive when the porosity is increased from 4.24 to 8.75%, but varies substantially when it is increased from 8.75 to 14.58%. The force and moment ratio (i.e. the ratio of the force or moment on the vertical cylinder, when it is encircled by the perforated caisson to the force or moment on the cylinder without any protection around it) reduces with increased wave height, H, and wave length, L, whereas the wave height ratio (ratio of the wave height at a point in the vicinity of the structure to the incident wave height) is less sensitive for the varying H and L. A new non-dimensional parameter, p^1.5 (D/L)/(H/d), is introduced to predict the in-line force and moment on the inner vertical cylinder, where d is local water depth, D is the diameter of the inner cylinder and p is the porosity of the outer caisson in percentage. Simple predictive equations for forces, moments and water surface fluctuations are provided.     

Comparison of flow patterns in the near wake of a circular cylinder and a square cylinder placed near a plane wall

The flow patterns in the near wake of a cylinder (either circular or square in shape, D=25 mm) placed in the proximity of a fully developed turbulent boundary layer (thickness δ=0.4D) are investigated experimentally using particle image velocimetry (PIV). The effects of changing the gap height (S) between the cylinder bottom and the wall surface, over the gap ratio range S/D=0.1–1.0, have been investigated. The results show that both the ensemble-averaged and instantaneous flow fields are strongly dependent on S/D. The flow patterns for the two types of cylinders share many similarities with respect to the change in S/D, such as the reduced recirculation length and increased velocity fluctuation in the near wake with increasing S/D, as well as the trend of suppression of vortex shedding at small S/D and onset of vortex shedding at large S/D. However, developments of the shear layers, in terms of wake width, flow curvature, etc., are considerably different for these two types of cylinders. In general, the wake development and momentum exchange for the square cylinder are slower those for the circular cylinder at the same gap ratio. Correspondingly, it is shown that the periodic vortex shedding is delayed and weakened in the case of square cylinder, as compared to that of the circular cylinder at the same S/D.     

张力腿平台绕流数值模拟分析

在一定来流条件下,张力腿平台(tension leg platform,简称TLP)的立柱后缘出现周期性的交替旋涡脱落,致使立柱受到垂直于来流方向的升力和平行于来流方向的阻力作用,导致TLP产生大幅度往复运动,显著增加平台结构和系泊系统的负载。目前,关于单柱、多柱结构绕流问题的研究较多,但对于TLP绕流特性的研究较少,机理尚存不明确的地方。为研究TLP的绕流力变化情况和流场特征,开展了数值模拟分析。利用计算流体动力学数值模拟软件,基于雷诺平均(RANS)法和分离涡模拟(DES)法对TLP绕流场进行仿真分析,揭示了TLP的绕流特性。结果表明,在3种来流角度和多个折合速度V_r下,TLP绕流的流体力系数存在明显差异,升力系数时域曲线呈现脉动性。TLP的上、下游立柱间存在明显的相互作用,影响了旋涡的形成与发展。TLP的旋涡脱落大多集中在平台固有频率附近,且在V_r=7,来流角度为0°时,升力系数频谱峰值最大,旋涡脱落集中。     

Experimental investigation of the forces and motion responses of two interfering VIV circular cylinders at various tandem and staggered positions

Experimental results are presented in the paper of two elastically supported rigid circular cylinders subjected to steady flows in a flume. The two cylinders were initially placed at various tandem and staggered positions with one in the wake of the other when subject to the steady flows. The in-line centre-to-centre distance varied from 2 to 5 diameters whilst the cross-flow distance from 0 to 2 diameters. The nominal Reynolds numbers were in the sub-critical regime and ranged from 1.12 × 10⁴ to 5.52 × 10⁴, and the nominal reduced velocities from 1.78 to 8.77. The damping ratio of the test set-up is low at 0.003 which gives a combined mass-damping parameter of 0.0046. Both the cylinders were free to respond in both the in-line and the cross-flow directions. The cylinder motion was measured simultaneously with the hydrodynamic loading in the two directions. It was found that the motion trajectories of the downstream cylinder show qualitative difference depending upon whether it is in tandem with the upstream cylinder or in the wake with a transverse offset. The VIV response of the downstream cylinder is dependent upon the reduced velocity of the upstream cylinder and its own reduced velocity based upon the actual mean wake velocity. The drag amplification of the downstream cylinder in the wake appears to be fundamentally different from that of a single VIV cylinder in isolation. Furthermore, unlike the two fixed cylinders in cross flow, the downstream cylinder undergoing VIV no longer experiences a marked non-zero mean lift. The upstream cylinder is largely unaffected by the downstream cylinder when the initial spacing is greater than 3 diameters. On the other hand, the motion response of and the fluid loading on the downstream cylinder are strongly influenced by the upstream cylinder in the spacing range tested.     

Assessment of Risk Due to Debris Flow and Its Application to a Marine Environment

The purpose of this study is to evaluate the behavior and mechanism of a debris flow on various slopes through numerical simulation. The numerical simulation consisted of using equations related to mass conservation and momentum conservation in order to consider erosion and deposition, and the Finite Difference Method was applied. As the inflow water discharge in the upstream of the channel increases, the curve of the water discharge exhibits instability and, as time passes, the fluctuation of the high water discharge continues. In regions where the mountain areas and the ocean are connected, it is deduced that the high level of sediment concentration can greatly affect the environment surrounding the ocean. The numerical model of this study was applied in Kangwon Province of South Korea. The results show that when the debris flow reaches downstream, the flow discharge and water flow depth increase. Erosion occurs more than deposition and much of the sediment runs off downstream. The result of the simulation performed at point of sediment discharge runoff is 114,216 m³. This study will provide useful information in predicting disasters caused by debris flow and in planning for various countermeasures to prevent debris-flow-related disasters.     

Numerical simulation of local scour around two pipelines in tandem using CFD–DEM method

In the field of offshore oil and gas engineering, the arrangement of multiple pipelines are becoming more common, the spacing between the pipelines and the incoming stream velocity will significantly affect the scouring process around the pipelines. In this study, the effect of space ratio (G/D) and the stream velocity on the scouring process around two pipelines in tandem are investigated using the coupled approach of computational fluid dynamics (CFD) and discrete element method (DEM). Here G is the spacing between the pipelines and D is the diameter of the pipeline. Specifically, the effect of space ratio and the stream velocity are discussed by simulating the gap ratio (G/D) between two pipelines ranging from 1 to 3 with an interval of 1, under the stream velocity U = 0.5,1 and 2 m/s, The results indicate that when G/D ≤ 2, the equilibrium scour depth below the upstream pipeline (S1) is slightly larger than that under the downstream pipeline (S2), S1 and S2 slightly increase as the gap ratio increases. Whereas for G/D > 2, the equilibrium scour depth beneath the upstream pipeline is slightly smaller than that under the downstream pipeline, S1 and S2 slightly decrease as the gap ratio increases. Furthermore, the scour depths are highly dependent on and positively related to the incoming stream velocity, the equilibrium bed profiles are similar under the same incident stream velocity with different gap ratios.     

Dynamic interaction between two parallel submarine pipelines considered vortex-induced vibration and local scour

Abstract

The mechanism of local scour under two vibrating pipelines is investigated numerically in this research. A sediment scour model is adopted to estimate the motion of sediment. The general moving objects model, which is dynamically coupled with fluid flow, is set up to simulate the vortex-induced vibration (VIV) of the pipeline. The sediment scour model and pipeline vibration model are verified with the previous experimental results and show good agreement. Then, the coupling effects between the pipeline vibration and the local scour are investigated numerically. The effects of G/D (the ratio of the distance between the two pipelines to the diameter of the pipelines) on the local scour and the VIV of the pipeline are examined. The results indicate that the maximum scour depth under the vibrating pipelines is much larger than the scour depth under the fixed pipelines. Due to the shadowing effect of the upstream pipeline, the maximum scour depth under the upstream pipeline is deeper than that under the downstream pipeline. The pipeline vibration magnitude is closely related to the strength of the vortex that sheds behind the pipeline. The effect of G/D on the shape and strength of the vortices that shed behind the pipelines is significant.     

Patterns and rates of fluid flow during burial diagenesis of the Middle Jurassic Lincolnshire Limestone, eastern England

A phase of ferroan burial calcite from the Middle Jurassic Lincolnshire Limestone exhibits a systematic spatial arrangement of oxygen isotopic characteristics. Mean δ¹⁸O values of the ferroan calcites from each of 15 core and outcrop localities over a study area 25 × 25 km were obtained. These values show a marked depletion from west to east across the study area of approximately 3‰, such that the oxygen isotopic composition of the ferroan calcites can be contoured. The systematic change in oxygen isotopic composition across the study area is believed to have recorded the thermal gradient in the limestone during ferroan calcite precipitation. This thermal gradient can be partially attributed to approximately 200 m of differential burial of the Lincolnshire Limestone across the study area during the Chalk deposition, with a maximum burial of 550 m to the east of the area at this time. A component of up-dip fluid flow (from east to west) through the formation is required to generate the temperature enhancements above those predicted for conduction alone by simple differential burial. Using a finite-difference step computer program, rates of fluid flow during ferroan calcite precipitation are calculated to be approximately 25 m/year. This rate of fluid flow is considerably greater than rates usually predicted for buried sedimentary basins. The causes of such rapid, probably relatively short-lived flow-rates may be the sudden dewatering of adjacent shales, the release of overpressure within the formation of interest, seismic pumping, or fluid circulation round a supracrustal convective loop.     

Numerical simulation of flow around a circular cylinder close to a flat seabed at high Reynolds numbers using a k–ε model

High Reynolds number flows around a circular cylinder close to a flat seabed have been computed using a two-dimensional standard high Reynolds number k–ε turbulence model. The effects of gap to diameter ratio, Reynolds number and flat seabed roughness for a given boundary layer thickness of the inlet flow upstream of the cylinder have been investigated. Hydrodynamic quantities and the resulting bedload transport have been predicted, and the vortex shedding mechanisms have been investigated. Predictions of hydrodynamic quantities around a cylinder located far away from the bed (so that the effect of the bed is negligible) are in satisfactory agreement with published experimental data and numerical results obtained for the flow around an isolated cylinder. Results for lower Reynolds number flows have also been computed for comparison with the high Reynolds number flow results. Overall it appears that the present approach is suitable for design purposes at high Reynolds numbers which are present near the seabed in the real ocean.     

Heave motion suppression of a Spar with a heave plate

Vortex shedding flow of an oscillating vertical cylinder with a disk attached at its keel is considered. This configuration is of interest for the offshore oil and gas industry. A finite difference method is employed to solve the incompressible Navier–Stokes equations in the primitive-variables formulation. Test cases were used to guide selection of the size of flow domain, numerical parameters, and to verify that the resultant method was both convergent and accurate. Numerical simulations have shown that the geometry configurations of the cylinder and disk, such as aspect ratio of the disk t_d/D_d and diameter ratio, D_d/D_c have significant influence on the vortex shedding modes and associated hydrodynamic properties, e.g. hydrodynamic damping and added mass coefficients. These in turn affect the performance in heave motion control of the structures.     

Geometry of scour hole around, and the influence of the angle of attack on the burial of finite cylinders under combined flows

Experiments were conducted to investigate the geometry of the scour hole and flow structure around short cylinders under the action of waves alone (WA) and combined flows (CF). The study is aimed at better understanding the dynamics of isolated objects on a sandy floor under oscillatory flows as occurs in shallow water regions in coastal areas. Flow velocities within the fluid core were recorded and 3D mapping of the bottom was performed with sub-aquatic acoustic sensors. Experiments were conducted for cylinder Reynolds wave number and Keulegan-Carpenter number within the ranges 10⁴R_e1.7×10⁵ and 2KC71, respectively. The present experimental evidence shows that the geometric characteristics of the scour hole (length and width) depend primarily on the Keulegan-Carpenter number (KC) and the cylinder aspect ratio (a_r=L_c/D). The effect of variation in the angle of attack of the flow with respect to the cylinder main axis was also investigated. Initial orientations of zero and ninety degrees were found to be stable while cylinders with intermediate initial orientations tended to orientate their main axes perpendicular to the flow direction. The final angle of orientation was found to be primarily a function of the Shields parameter, θ, and the initial angle of attack, α_i.     

Wave interaction with a submarine pipeline in clayey soil due to random waves

The wave pressure and uplift force due to random waves on a submarine pipeline (resting on bed, partially buried and fully buried) in clayey soil are measured. The influence of various parameters viz., wave period, wave height, water depth, burial depth and consistency index of the soil on wave pressures around and uplift force on the submarine pipeline was investigated. The wave pressures were measured at three locations around the submarine pipeline (each at 120° to the adjacent one). It is found that the wave pressure and uplift force spectrum at high consistency index of the soil is smaller compared to that of low consistency index. Just burying the pipeline (e/D=1.0) in clayey soil reduces the uplift force to less than 60% of the force experienced by a pipeline resting on the seabed (e/D=0.0) for I_c=0.33.     

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.     

集成于海上风电钢管桩基础上的升降式网箱结构波浪场特性研究

对安装在海上风电钢管桩基础上的升降式网箱结构的波浪场特性进行研究,掌握升降式网箱结构内部及结构后方水体的运动特征及速度场变化情况。基于OpenFOAM软件包开发了波、流与孔隙介质结构相互作用的数值计算模型,将网箱结构按等效阻力简化为多孔介质结构,开展升降式网箱结构的波浪场特性研究。研究结果表明：网箱结构对流体具有一定的阻流作用,网箱结构内部速度场得到一定程度的减小,网箱结构背浪侧也有一定的速度衰减区域;对比网箱结构顶部不同潜深条件下的网箱结构内部流场特征,网箱结构顶部潜深在1/4D~1/2D（D为水深）范围内网箱结构内部流场速度最小、流场最为稳定、速度分布均匀,网箱结构向浪侧前方和背浪侧后方流场波动较小。所得结论表明在钢管桩基础上安装升降式网箱结构时需要关注的网箱结构对流场特性的影响,充分考虑网箱结构阻力对流速的作用,掌握极端工况下升降式网箱结构保持优良养殖环境需要下潜的高度范围,以保障升降式网箱的安全。