串列双圆柱流致振动数值模拟

基于开源程序OpenFOAM和动网格技术,利用切应力平衡法建立水流作用下的海上风电基础局部冲刷数学模型。通过模拟结果与实验数据的对比发现,所建立的冲刷数学模型能够合理反映圆柱型单桩基础周围的水流结构,冲刷深度与实验结果吻合较好。     

基于高雷诺数的并联双圆柱绕流研究

文章采用基于边界瞬时涡量守候的拉格朗日方法(IVCBC),结合并列双圆柱的特点,建立了双圆柱绕流数值计算模型。对高雷诺数下Re=6×10~4,间隙率为T/D=1.1~7的并联圆柱双圆柱二维绕流特性进行了研究。提出了双圆柱间隙中点的速度区别宽窄尾流的新方法。分别讨论了流体力系数随间距增加的特点,脉动流体力的特点,尾流特征以及斯托哈尔数特征。研究发现:区别宽窄尾流的新方法是可靠的;在双圆柱尾流附近有五种尾流模型存在;同时发现了在宽窄尾流的频率,存在一个中间频率。     

Numerical simulation of piggyback pipelaying under current loadings

For piggyback pipelaying operations, current-induced force and its effect on the piggyback pipe have not been thoroughly studied. In the present study, an improved method in hydrodynamic load calculation and structural modelling is proposed to simulate the pipelaying of a piggyback pipeline. In order to obtain the mean drag and lift force coefficients for the piggyback pipeline subjected to different inflow angles, two-dimensional Computational Fluid Dynamics (CFD) simulations are performed by modelling the piggyback pipeline as two cylinders attached to each other without gap. Then, the acquired force coefficients are used to calculate the hydrodynamic loads through a user-defined function in OrcaFlex based on a cross-flow principle approach. The interaction between the pipeline and the piggyback cable is modelled using two types contact elements which are ring penetrator and non-penetrating contact. The present proposed method is compared with other two widely used engineering methods based on (1) the equivalent diameter and (2) two separate cylinders without accounting for hydrodynamic interaction, in terms of the top tension, and the bending moments at Hang-off Clamps (HOC) and sagbend of the pipeline. The comparison shows that the two widely used engineering methods are not always conservative in force and response predictions. Hence, it is important to consider the hydrodynamic and structural interactions between the piggyback cable and the pipeline. With different current directions, the bending moments at the HOC predicted by the present method vary from 40% lower to 100% higher than those predicted by the two widely used engineering methods.     

Empirical model of the wake-induced lift force on a cylinder with low mass ratio

The vortex-induced vibrations (VIV) of two flexible circular cylinders in a tandem configuration were studied numerically for spacing ratios ranging from 6 to 18 and the reduced velocities ranging from 2.35 to 12.59. The VIV response amplitude, response frequency, fluid force, pressure distribution and vortex structure of the tandem cylinders with different spacing ratios under different reduced velocities were compared. The results indicate that there is a great difference between the lift forces on the downstream and upstream cylinders. The lift coefficient of the downstream cylinder undergoing the wake-induced vibrations (WIV) is larger than that of the upstream cylinder, and the dominant frequency curves of the lift coefficients of the upstream and downstream cylinders separate. It can be found that the length and intensity of the wake are quite different under different reduced velocities and spacing ratios, and the reattachment positions between the wake and the downstream cylinder are different, which leads to a great change in the flow around the downstream cylinder and have a great effect on the wake-induced lift force on the downstream cylinder. Considering these factors, an empirical model for the wake-induced lift force on a cylinder with low mass ratio was proposed and verified.     

带有变化率的同轴圆柱体受到的波浪力研究(英文)

Wave diffraction of two concentric porous cylinders with varying porosity was studied by using an analytical method based on eigenfunction matching.The fluid domain around the cylinders is divided into three sub-domains and in each sub-domain an eigenfunction expansion of the velocity potential is obtained by satisfying the Laplace equation,the boundary conditions on the free surface and on the sea bed.The unknown coefficients of eigenfunction expansions are determined by boundary conditions on the porous hulls.In the paper,the boundary conditions are based upon the assumption that the flow in the porous medium is governed by Darcy’s law.Two porous-effect parameters applied on two porous cylinders are functions of the vertical coordinate instead of the constant.Wave loading on the outer and inner cylinder is presented in the numerical results.     

四柱结构近场干涉的粘流分析（英文）

The hydrodynamic analysis of multi-floating bodies is important and widely used in marine engineering. In this study, we systematically simulated the wave diffraction problem of a fixed vertical four-cylinder structure in regular waves in the time domain in a viscous numerical wave tank. The hydrodynamic interaction of waves with a bottom-mounted structure consisting of four vertical cylinders arranged at the corners of a square shows a complicated interference phenomenon. In this paper, we illustrate and analyze the run-up around the structure and the corresponding wave forces. To investigate the viscous effect on the near-trapping phenomenon, we pay particular attention to investigating the waves near-trapped inside the four-cylinder structure,and make a comparative study of the viscous-and inviscid-flow solutions with the experimental measurements. The results show that the maximum wave elevation occurs on the inner side of the leeside cylinder, and that the wave elevations on the outer side of the cylinders are lower than those on the inner side. We can conclude that viscosity has an obvious damping effect on wave elevations inside the structure. The cylinders show a tendency to drift apart from each other when the near-trapping phenomenon occurs.     

Numerical prediction of the modal response of flexible cylinders in cross-flow with a current dependent form of damping

A quadratic eigenvalue problem (QEP) was posed in order to study the dynamics of flexible cylinders in cross-flow, simulating slender offshore structures such as risers, catenaries or tendons. The Euler–Bernoulli equation was used to model the structure assuming a fluid loading model, and yielding a quadratic eigenvalue problem that included a form of damping dependent not only on the structural damping itself, but also on the free stream velocity and the fluid force coefficients. We solved the QEP using the finite element method. We also derived a simplified analytical solution in this work for comparison with the QEP, however this solution does not consider changes in tension along the length of the cylinder as the QEP does. In our study, the QEP solutions were first validated against the simplified analytical solution, and also against a well-known experimental dataset obtained in 2003, in which a flexible circular cylinder model was used to model the dynamics of a riser undergoing multi-mode vortex-induced vibrations.     

简单Green函数法求解三维水动力系数

本文提出用简单Green函数1/r求解无航速的三维物体在水面、水下或在无界域中作各向振荡时水动力系数的数值方法。对于无穷远处扩散条件转换成辐射面条件的问题作了专门的考虑,即用圆柱将物体包围在其内部,该圆柱将流场分成内、外两部分,外部域中速度势的特征函数展开式与内部域中用Green公式求得的解相匹配。 根据本计算方法编制的程序对一些典型情况作了试算,其中包括在水面振荡的半球、方盒和圆盘的水动力系数曲线;水面及水底对圆球水动力系数的影响等。为考核本方法,与其它计算方法得到的结果作了比较。对方盒的三维结果与切片法作了比较。     

两圆柱体在波浪载荷作用下的运动响应分析

基于势流理论和波浪辐射/绕射理论,利用水动力分析软件(AQWA)研究两圆柱体在波浪作用下的运动响应特性,给出不同入射角与不同频率规则波作用下,小圆柱体单独运动时的六自由度响应幅值算子(RAOs)及运动响应,并与两个并排圆柱体六自由度响应幅值算子(RAOs)及运动响应作比较。结果表明:大圆柱体的干扰将使小圆柱体的纵摇、横摇、横荡、垂荡运动加剧,对小圆柱体漂移力也有较大影响,但对小圆柱体的纵荡与首摇运动影响不大。计算结果与传统方法计算值相近,证明此软件合理,可为两船在波浪中的干扰模拟计算提供参考。     

Exciting forces for a wave energy device consisting of a pair of coaxial cylinders in water of finite depth

Two coaxial vertical cylinders-one is a riding hollow cylinder and the other a solid cylinder of greater radius at some distance above an impermeable horizontal bottom,were considered.This problem of diffraction by these two cylinders,which were considered as idealization of a buoy and a circular plate,can be considered as a wave energy device.The wave energy that is created and transferred by this device can be appropriately used in many applications in lieu of conventional energy.Method of separation of variables was used to obtain the analytical expressions for the diffracted potentials in four clearly identified regions.By applying the appropriate matching conditions along the three virtual boundaries between the regions,a system of linear equations was obtained,which was solved for the unknown coefficients.The potentials allowed us to obtain the exciting forces acting on both cylinders.Sets of exciting forces were obtained for different radii of the cylinders and for different gaps between the cylinders.It was observed that changes in radius and the gap had significant effect on the forces.It was found that mostly the exciting forces were significant only at lower frequencies.The exciting forces almost vanished at higher frequencies.The problem was also investigated for the base case of no plate arrangement,i.e.,the case having only the floating cylinder tethered to the sea-bed.Comparison of forces for both arrangements was carried out.In order to take care of the radiation of the cylinders due to surge motion,the corresponding added mass and the damping coefficients for both cylinders were also computed.All the results were depicted graphically and compared with available results.     

Vortex-induced vibrations of piggyback pipelines near the horizontal plane wall in the upper transition regime

Slender subsea structures like pipelines, jumpers and umbilicals when exposed to currents may experience vortex-induced vibrations (VIV), which can shorten their fatigue life and increase the risk of structural failure. In the present study, flow around different configurations of a piggyback pipeline close to a flat seabed has been investigated using the two-dimensional (2D) Unsteady Reynolds-Averaged Navier Stokes (URANS) equations with the k − ω Shear Stress Transport (SST) turbulence model. The Reynolds number (based on the free stream velocity and large cylinder diameter) is equal 3.6 × 10⁶ corresponding to upper-transition regime. The drag forces acting on the cylinders and base pressure coefficient value are well predicted by the present simulations, while the other hydrodynamic quantities (root-mean-square lift coefficient, Strouhal number) are predicted reasonably well as compared to published experimental data. The piggyback pipeline in the present study is modeled as two circular cylinders with a diameter ratio d/D = 0.2 (D denoting diameter of the large cylinder, d is diameter of the small cylinder). These two cylinders are clamped together at a distance G/D = 0.2. The two rigidly coupled cylinders are elastically supported and free to vibrate in two degrees of freedom. The effects on the vibration amplitudes and hydrodynamic forces are analyzed. The flow structures around the cylinders are investigated to explain the variations in observed structural responses. Depending on the angular position (α) of the small cylinder, the lock-in regime is narrower (α = 0°) or significantly wider (α = 180°) when compared to that of a single cylinder.     

块体形状对水流拖曳力的影响

在计算块体受到的水流拖曳力时,常借鉴针对球体所提出的公式,然而它并不能客观地反映出块体形状对拖曳力产生的影响。文中将块体形状概化为长方体,通过水槽试验测量了水流作用下不同尺度的大颗粒块体受到的水平拖曳力,研究了块体形状对水流拖曳力的影响,并提出了考虑块体形状影响的拖曳力计算公式。     

圆柱绕流研究进展及展望

圆柱绕流问题广泛存在于水利工程、建筑工程和环境工程之中。通过对前人圆柱绕流研究的成果整理和分析,将其分为解析研究、物理模型研究和数值模拟研究3个阶段,总结了单圆柱、双圆柱以及圆柱群绕流在流动形态和圆柱受力特性方面的研究进展。分析发现,现有研究主要集中于简单条件下单柱和双柱的受力特性和尾流流场结构,今后需进一步开展柱群在海洋立管、高桩码头和环境工程等复杂条件下的研究,开发尾流形态的主动控制技术。     

Effect of mass ratio on hydrodynamic response of a flexible cylinder

The effect of the mass ratio on the flow-induced vibration (FIV) of a flexible circular cylinder is experimentally investigated in a towing tank. A Tygon tube with outer and inner diameters of 7.9 mm and 4.8 mm, respectively, was employed for the study. The tube was connected to a carriage and towed from rest to a steady speed up to 1.6 m/s before slowing down to rest again over a distance of 1.6 m in still water. Reynolds number based on the cylinder’s outer diameter was 800–13,000, and the reduced velocity (velocity normalized by the cylinder’s natural frequency and outer diameter) spanned from 2 to 25. When connected, the cylinder was elongated from 420 mm to 460 mm under an axial pre-tension of 11 N. Based on the cylinder’s elongated length, the aspect ratio (ratio of the cylinder’s length to outer diameter) was calculated as 58. Three mass ratios (ratio of the cylinder’s structural mass to displaced fluid mass, m*) of 0.7, 1.0, and 3.4 were determined by filling the cylinder’s interior with air, water, and alloy powder (nickel-chromium-boron matrix alloy), respectively. An optical method was adopted for response measurements. Multi-frequency vibrations were observed in both in-line (IL) and cross-flow (CF) responses; at high Reynolds number, vibration modes up to the 3^rd one were identified in the CF response. The mode transition was found to occur at a lower reduced velocity for the highest tested mass ratio. The vibration amplitude and frequency were quantified and expressed with respect to the reduced velocity. A significant reduced vibration amplitude was found in the IL response with increasing mass ratios, and only initial and upper branches existed in the IL and CF response amplitudes. The normalized response frequencies were revealed to linearly increase with respect to the reduced velocity, and slopes for linear relations were found to be identical for the three cases tested.     

全非线性方法对于沉浸双圆柱定常前行运动引起自由表面流的研究（英文）

The free surface flow generated by twin-cylinders in forced motion submerged beneath the free surface is studied based on the boundary element method. Two relative locations, namely, horizontal and vertical, are examined for the twin cylinders. In both cases, the twin cylinders are starting from rest and ultimately moving with the same constant speed through an accelerating process. Assuming that the fluid is inviscid and incompressible and the flow to be irrotational, the integral Laplace equation can be discretized based on the boundary element method. Fully-nonlinear boundary conditions are satisfied on the unknown free surface and the moving body surface. The free surface is traced by a Lagrangian technique. Regriding and remeshing are applied, which is crucial to quality of the numerical results. Single circular cylinder and elliptical cylinder are calculated by linear method and fully nonlinear method for accuracy checking and then fully nonlinear method is conducted on the twin cylinder cases, respectively. The generated wave elevation and the resultant force are analysed to discuss the influence of the gap between the two cylinders as well as the water depth. It is found that no matter the kind of distribution, when the moving cylinders are close to each other, they suffer hydrodynamic force with large absolute value in the direction of motion. The trend of force varying with the increase of gap can be clearly seen from numerical analysis. The vertically distributed twin cylinders seem to attract with each other while the horizontally distributed twin cylinders are opposite when they are close to each other.     

Experimental investigation on hydrodynamics of floating cylinder in oscillatory and steady flows by forced oscillation test

In this paper the hydrodynamic characteristics of a floating cylinder are investigated via forced oscillation experiments in towing tank. The effects of Keulegan–Carpenter number, Reynolds number, reduced velocity and overtopping on hydrodynamics of the floating cylinder in oscillatory and steady flow are studied. The results show a considerably difference of the hydrodynamic characters between the floating and the fully immerged cylinders due to the influences of free surface. The growth of the reduced velocity, a proven notable effect on hydrodynamics, will lead to the increase of added mass coefficient and the decrease of drag coefficient. Meanwhile the overtopping, a particular phenomenon for the floating cylinder, render the added mass coefficients reach up to 3.6 while for the drag coefficient small influences were made.     

Vortex-induced vibration (VIV) dynamics of a tensioned flexible cylinder subjected to uniform cross-flow

This paper presents the experimental results of a study on the effects of axial applied tension on the vibration amplitude, the suppression of vibration, hydrodynamic force coefficients and in-line (IL) and cross-flow (CF) frequency responses during vortex-induced vibration of a horizontally mounted flexible cylinder with a low mass ratio (cylinder’s mass/mass of displaced water), low bending-stiffness, and high aspect ratio (length/diameter 200) in the subcritical Reynolds number regime (Re = 1000–16000). The effect of tension is studied by applying four different tensions. It was revealed that higher applied tensions, which reduce the vibration amplitude, could significantly raise the hydrodynamic lift force coefficient. In addition, higher applied tensions generate narrower lock-in bandwidths. After the highest vibration amplitude and during the region of lower vibration amplitudes, within the first lock-in region (in the first mode of vibration), power spectral densities show broad bandwidth, while within other regions and higher modes they appear narrow-banded. The ratio of the dominant IL to CF frequency is approximately equal to 2.0, except for the lower reduced velocities, where the ratio values reach 3.83 for the highest tension accompanied by widening of the region in which this ratio is over 2.0. This ratio is 2.76 for the lowest applied tension with a narrower region.     

3D numerical modeling of wave forces on tandem fixed cylinders using the BEM

In this paper a 3D numerical model was developed to study the complicated interaction between waves and a set of tandem fixed cylinders.The fluid was considered to be inviscid and irrotational.Therefore,the Helmholtz equation was used as a governing equation.The boundary element method（BEM） was adopted to discretize the relevant equations.Open boundaries were used in far fields of the study domain.Linear waves were generated and propagated towards tandem fixed cylinders to estimate the forces applied on them.Special attention was paid to consideration of the effect on varying non-dimensional cylinder radius and distance between cylinders,ka and kd on forces and trapped modes.The middle cylinder wave forces and trapped modes in a set of nine tandem cylinders were validated utilizing analytical data.The comparisons confirm the accuracy of the model.The results of the inline wave force estimation on n tandem cylinders show that the critical cylinder in the row is the middle one for odd numbers of cylinders.Furthermore the results show that the critical trapped mode effect occurs for normalized cylinder radiuses close to 0.5 and 1.0.Finally the force estimation for n tandem cylinders confirms that force amplitude of the middle cylinder versus normalized separation distance fluctuates about that of a single cylinder.     

基于多目标遗传算法的复式转叶舵机结构优化

针对舰船转叶舵机运动时舵叶与水动力之间存在强力位耦合的问题,设计一种新型复式液压摆动缸转舵机构。为减小结构整体占地面积,基于多目标遗传算法(Multi-Objective Genetic Algorithm, MOGA)建立以复式摆动缸壳体内径、力矩解耦缸转子内径、驱动缸转子轮毂直径、动静叶片高度、动静叶片宽度和壳体厚度几何尺寸为设计变量的多目标优化模型。在使用遗传算法(Genetic Algorithm, GA)求解非线性约束条件下的多变量优化问题时,通过设计罚函数,将适应度函数与约束条件分离以处理约束条件,根据设定的准则择优筛选最优尺寸参数。通过有限元软件ANSYS对摆动缸壳体、力矩解耦缸转子和驱动缸转子进行应力应变分析。结果表明,在满足结构强度、刚度和安全因数的前提下,复式摆动缸结构尺寸最优,验证优化后结构的可靠性。     

Numerical simulation with a macroscopic CFD method and experimental analysis of wave interaction with fixed porous cylinder structures

Porous structures have been widely applied in the coastal and ocean engineering due to their wave energy dissipation mechanism. The macroscopic computational fluid dynamics (CFD) approach where the quadratic pressure drop condition of porous surface is introduced to model the wave interaction with porous cylinders. A series of CFD simulations of waves interacting with a single porous cylinder and the combined structure of a porous cylinder with a concentric inner solid column are performed, with corresponding tank tests conducted. The CFD method is compared with experiments, linear potential model, and the quadratic BEM (boundary element method) model. The effects of porosity and porous cylinder radius on wave force and wave heights inside porous cylinder are analyzed to evaluate the performance of porous shell reducing wave loads and wave surface elevation, and the wave force variation with incident wave amplitudes are also investigated. The results demonstrate that the established CFD model is reliable for engineering analysis and thereby being of great significance for reference purpose in the CFD simulations of waves interacting with porous structures.