WEC-透空式防波堤集成装置能量转换效率研究

基于成本共享理念的波浪能发电装置(WEC)的开发与设计为降低建设成本提供了新的研究思路,应用计算流体力学方法对在透空箱式防波堤前附加垂荡浮子的集成装置进行数值模拟研究,主要研究了流体黏性和非线性PTO系统以及浮子形状对此类集成装置能量转换效率的影响。结果表明,此类集成装置可以获得较高的波能捕获宽度比(CWR),最高可达0.7,可以较好地利用反射波波能。共振区间上受黏性影响较大,相比于线性无黏理论解,CWR下降明显;实现了2种非线性PTO系统的模拟,有待于进一步的优化未获得更高的CWR值;浮子形状优化效果明显,采用圆底形浮子受黏性影响更小,可以获得更大的CWR值。此种集成模型的模拟和研究可以为新型防波堤设计和现有防波堤的改造提供思路和参考。     

Simulation of Random Waves and Associated Laminar Bottom Shear Stresses

This work presents a new approach for simulating the random waves in viscous fluids and the associated bottom shear stresses. By generating the incident random waves in a numerical wave flume and solving the unsteady two-dimensional Navier-Stokes equations and the fully nonlinear free surface boundaiy conditions for the fluid flows in the flume, the viscous flows and laminar bottom shear stresses induced by random waves axe determined. The deterministic spectral amplitude method implemented by use of the fast Fourier transform algorithm was adopted to generate the incident random waves. The accuracy of the numerical scheme is confirmed by comparing the predicted wave spectrum with the target spectrum and by comparing the nanlerical transfer function between the shear stress and the surface elevation with the theoretical transfer function. The maximum bottom shear stress caused by random waves, computed by this wave model, is compared with that obtained by Myrhaug＇ s model （1995）. The transfer function method is also employed to determine the maximum shear stress, and is proved accurate.     

Viscous free surface effect on coastal embankment hydrodynamics

A finite-differnece method was used to calculate the nonlinear hydrodynamic pressures acting on the coastal embankment faces by seismic-wave actions. The nonlinearity of free surface flow, convective acceleration, viscosity and surface tension of fluid are included in the analysis. The kinematic and dynamic free surface boundary conditions are employed for calculating the horizontal fluid velocity, pressure at the free surface and the surface profile of the fluid. The time-dependent water surface is transformed to the horizontal plane, and the flow field is mapped onto a rectangular, making it convenient to model the complex sea bottom geometry and the wavy water surface by the finite-difference method. Fully nonlinear and weakly nonlinear dynamic free surface conditions are used and compared. The effects of surface tension of fluid are also discussed. The nonslip boundary condition is applied on the most part of the interface between fluid and solid face, except the region near the intersection between free surface and wall face. The numerical results are presented for various water depths and ground motion intensities, and their associate viscous effects on coastal embankment hydrodynamics are discussed.     

A Numerical Study on Wave-Mud Interaction

1 .Introduction The dynamics of soft mud under surface water waves is of great importance to the sedimentationprocesses in approach channels and harbors ,and has long been drawing attention. Advancements innumerous engineering applications inthe shallowco…     

Dynamics of arrays of floating point-absorber wave energy converters with inter-body and bottom slack-mooring connections

If point absorbers are employed in the extensive exploitation of the offshore wave energy resource, they should be deployed in arrays, the distance between the elements in the array being possibly tens of meters. In such cases, it may be more convenient that the array is spread moored to the sea bottom through only some of its elements, located in the periphery, while the other array elements are prevented from drifting and colliding with each other by connections to adjacent elements. An array of identical floating point absorbers located at the grid points of an equilateral triangular grid is considered in the paper. A spread set of slack-mooring lines connect the peripheric floaters to the bottom. A weight is located at the centre of each triangle whose function is to pull the three floaters towards each other and keep the inter-body mooring lines under tension. The power take-off system (PTO) is a linear damper activated by the buoy heaving motion. The whole system–buoys, moorings and power take-off systems–is assumed linear, so that a frequency domain analysis may be employed. Hydrodynamic interference between the oscillating buoys in array is accounted for. Equations are presented for a set of three identical point absorbers. This is then extended to more complex equilateral triangular grid arrays. Results from numerical simulations, with regular and irregular waves, are presented for the motions and power absorption of hemispherical converters in arrays of three and seven elements and for different mooring parameters and wave incidence angles. Comparisons are given with the unmoored and independently-moored buoy situations.     

南黄海和东海“人工水母”投放试验

1984—1986年,作者在南黄海和东海投放大量“人工水母”测量底层流。本文介绍了“人工水母”的投放、回收和漂流概况,根据这一实验结果,参考有关黄、东海海流研究的成果,绘出了调查试验海区的底层流模式。     

Effect of wave nonlinearity on fatigue damage and extreme responses of a semi-submersible floating wind turbine

Floating wind turbine has been the highlight in offshore wind industry lately. There has been great effort on developing highly sophisticated numerical model to better understand its hydrodynamic behaviour. A engineering-practical method to study the nonlinear wave effects on floating wind turbine has been recently developed. Based on the method established, the focus of this paper is to quantify the wave nonlinearity effect due to nonlinear wave kinematics by comparing the structural responses of floating wind turbine when exposed to irregular linear Airy wave and fully nonlinear wave. Critical responses and fatigue damage are studied in operational conditions and short-term extreme values are predicted in extreme conditions respectively. In the operational condition, wind effects are dominating the mean value and standard deviation of most responses except floater heave motion. The fatigue damage at the tower base is dominated by wind effects. The fatigue damage for the mooring line is more influenced by wind effects for conditions with small wave and wave effects for conditions with large wave. The wave nonlinearity effect becomes significant for surge and mooring line tension for large waves while floater heave, pitch motion, tower base bending moment and pontoon axial force are less sensitive to the nonlinear wave effect. In the extreme condition, linear wave theory underestimates wave elevation, floater surge motion and mooring line tension compared with fully nonlinear wave theory while quite close results are predicted for other responses.     

Viscosity and nonlinearity effects on the forces and waves generated by a floating twin hull under heave oscillation

Nonlinear hydrodynamics of a twin rectangular hull under heave oscillation is analyzed using numerical methods. Two-dimensional nonlinear time-domain solutions to both inviscid and viscous problems are obtained and the results are compared with linear, inviscid frequency-domain results obtained in [26] to quantify nonlinear and viscous effects. Finite-difference methods based on boundary-fitted coordinates are used for solving the governing equations in the time domain [2]. A primitive-variables based projection method [6] is used for the viscous analysis and a mixed Eulerian–Lagrangian formulation [11] for inviscid analysis. The algorithms are validated and the order of accuracy determined by comparing the results obtained from the present algorithm with the experimental results of Vugt [22] for a heaving rectangle in the free surface. The present study on the twin-hull hydrodynamics shows that at large and non-resonant regular frequencies, and small amplitude of body oscillation, the fluid viscosity does not significantly affect the wave motion and the radiation forces. At low frequencies however the viscosity effect is found to be significant even for small amplitude of body oscillation. In particular, the hydrodynamic force obtained from the nonlinear viscous analysis is found to be closer to the linear inviscid force than the nonlinear inviscid force to the linear inviscid force, the reason for which is attributed to the wave dampening effect of viscosity. Since the wave lengths generated at smaller frequencies of oscillation are longer and therefore the waves could have a more significant effect on the dynamic pressure on the bottom of the hulls which contribute to the heave force, the correlation between the heave force and the wave elevation is found to be larger at smaller frequencies. Because of nonlinearity, the wave radiation and wave damping force remained nonzero even at and around the resonant frequencies – with the resonant frequencies as determined in [26] using linear potential flow theory. As to be expected, the nonlinear effect on the wave force is found to be significant at all frequencies for large amplitude of oscillation compared to the hull draft. The effect of viscosity on the force, by flow separation, is also found to be significant for large amplitude of body oscillation.     

不同地形下单浮子集成装置水动力性能研究

为探讨不同地形特征对单浮子集成装置水动力特性的影响,本文基于黏性计算流体动力学理论,应用Star-CCM+软件建立二维数值波浪水槽,研究了规则波作用下受不同地形影响的单浮子式防波堤与波能转换装置集成系统的水动力性能.结果表明:在低频区,不同地形下单浮子集成系统的波能转换效率均大于无地形影响下单浮子集成系统的波能转换效率...     

均匀流中直立圆柱体绕流三维数值模拟

研究直立贯底圆柱体的三维粘性绕流问题。以不可压缩Navier-Stokes方程为控制方程,采用有限体积法和SIMPLE算法,建立了数值模拟方法。考察在不同水平和垂直断面上,圆柱体绕流产生的尾涡和流动速度场的分布特性,成功地数值模拟了直立贯底圆柱体绕流场的三维特性。结果表明,在考虑重力影响的情况下,直立圆柱体周围的流动具有明显的三维特性,而且沿圆柱体轴向不同断面上的尾涡分布是不相同的。     

Sound-pressure level calculations using the RRA algorithm fordepth-dependent speeds of sound valid at turning points and focal points

Sound-pressure level (SPL) calculations are made along individual ray paths for arbitrary, one-dimensional, depth-dependent speeds of sound using an enhanced version of the RRA (recursive ray acoustics) algorithm. The SPL calculations are valid (i.e., finite) at turning points and focal points and do not require the use of Airy functions. The SPL calculations include the effects of frequency-dependent volume attenuation and frequency-dependent attenuation due to surface and bottom reflections. The ocean surface and bottom are treated as boundaries between viscous fluid media. Although the ocean surface is modeled as a planar boundary, the bathymetry is an arbitrary function of horizontal range. Sound speed versus depth and bathymetric data are represented by orthogonal function expansions. Computer simulation results from preliminary test cases are presented     

两层粘性流体中圆柱体受迫振荡数值模拟

研究两层粘性流体中无限长水平圆柱体的受迫振荡问题。在湍流模式下,采用VOF方法追踪两层流体的内界面,基于动网格技术模拟圆柱体的运动边界,对均匀流中横向振荡圆柱体的绕流场进行了数值模拟。计算受迫振荡圆柱体的升力系数、阻力系数随时间的演化曲线和圆柱体的尾涡分布,以及圆柱体的受迫振荡激发两层流体内界面的扰动,并与均匀流体的情况进行了比较分析。研究表明,流体的两层分层效应对受迫振荡圆柱体的升阻力系数和尾涡分布特性都有显著影响,在水下输油气管道涡激振动特性的工程评估中,应考虑流体的密度分层效应。     

Quantitative evaluation of a passively leveled ocean bottom seismometer

A problem in the use of ocean bottom seismometers is the difficulty in leveling the sensors while ensuring good coupling to the seafloor. We have investigated the coupling characteristics of the seismic sensors in the new ONR ocean bottom seismometer. In the deployable sensor package for that instrument, a three-component seismometer set is suspended on a 2-axis passive leveling gimbal and is immersed in a viscous fluid. We report tests, conducted in a seismic vault, comparing the output of a gimbaled seismometer set to that of a set rigidly coupled to the ground. Our results show that the degree to which the gimbaled set is coupled to ground motion is a function of the viscosity of the coupling fluid. The coherence between the two sensor sets is poor (<0.4) at some frequencies within the band of interest (0.15 to 20 Hz) and on some components when the viscosity of the coupling fluid is comparatively low (14 Pa-s or 0.16 kSt kinematic viscosity). In addition, the outputs of some components over portions of this frequency band are attenuated and are phase-shifted relative to the outputs of the set rigidly coupled to the ground. Coherence and phase response similarity improve as the viscosity of the coupling fluid is increased. With a coupling fluid viscosity of 980 Pa-s (10 kSt), coherence and phase agreement between the two sensor sets is good (>0.9) across nearly the entire band of interest on all three components. A simple analytical model of the gimbaled seismometer set as a damped, driven, compound-pendulum provides a basis for understanding the test results.     

Viscous Effects on Wave Forces on A Submerged Horizontal Circular Cylinder

Numerical simulations are carried out for wave action on a submerged horizontal circular cylinder by means of a viscous fluid model, and it is focused on the examination of the discrepancies between the viscous fluid results and the potential flow solutions. It is found that the lift force resulted from rotational flow on the circular cylinder is always in anti-phase with the inertia force and induces the discrepancies between the results. The influence factors on the magnitude of the lift force, especially the correlation between the stagnation-point position and the wave amplitude, and the effect of the vortex shedding are investigated by further examination on the flow fields around the cylinder. The viscous numerical calculations at different wave frequencies showed that the wave frequency has also significant influence on the wave forces. Under higher frequency and larger amplitude wave action, vortex shedding from the circular cylinder will appear and influence the wave forces on the cylinder substantially.     

Quadratic Dissipation Effect on the Moonpool Resonance

This paper adopted a semi-analytical method based on eigenfunction matching to solve the problem of sharp resonance of cylindrical structures with a moonpool that has a restricted entrance. To eliminate the sharp resonance and to measure the viscous effect, a quadratic dissipation is introduced by assuming an additional dissipative disk at the moonpool entrance. The fluid domain is divided into five cylindrical subdomains, and the velocity potential in each subdomain is obtained by meeting the Laplace equation as well as the boundary conditions. The free-surface elevation at the center of the moonpool, along with the pressure and velocity at the restricted entrance for first-order wave are evaluated. By choosing appropriate dissipation coefficients, the free-surface elevation calculated at the center of the moonpool is in coincidence with the measurements in model tests both at the peak period and amplitude at resonance. It is shown that the sharp resonance in the potential flow theory can be eliminated and the viscous effect can be estimated with a simple method in some provided hydrodynamic models.     

Dynamics of gravity-capillary waves on the surface of a nonuniformly heated fluid

The method of slowly varying amplitudes is used to determine the reflection of surface gravity-capillary waves from regions of regular surface convection in the approximation of an ideal and homogeneous viscous fluid. A model of a quasi-steady distribution of temperature is employed. The conditions of both weak (nonresonant) interaction, when the wave passes almost completely through the convective zone, and strong resonant interaction are revealed. When interaction is strong, standing waves form in front of the convective region and the spectral composition of the wave packet changes. Viscosity and the accompanying boundary layers lead to the broadening of forbidden zones.     

Numerical modelling of a point-absorbing wave energy converter in irregular and extreme waves

Based on the Navier-Stokes (RANS) equations, a three-dimensional (3-D) mathematical model for the hydrodynamics and structural dynamics of a floating point-absorbing wave energy converter (WEC) with a stroke control system in irregular and extreme waves is presented. The model is validated by a comparison of the numerical results with the wave tank experiment results of other researchers. The validated model is then utilized to examine the effect of wave height on structure displacements and connection rope tension. In the examined cases, the differences in WEC’s performance exhibited by an inviscid fluid and a viscous fluid can be neglected. Our results also reveal that the differences in behavior predicted by boundary element method (BEM) and the RANS-based method can be significant and vary considerably, depending on wave height.     

基于SPH方法的开孔沉箱比尺效应研究

由于在前壁上设置了尺寸较小的孔,开孔沉箱受流体黏性力作用显著,依照弗劳德数相似准则设计模型存在比尺效应。为揭示比尺效应,建立了模拟波浪与开孔沉箱相互作用的光滑粒子流体动力学(SPH)模型。其中流体运动由连续性方程和Navier-Stokes方程控制,固壁边界由改进的动力边界粒子施加。模型收敛性通过分析不同粒子分辨率下的波浪反射系数得到,模型精度通过比较计算与理论波浪反射系数证明。使用经过验证的SPH模型,计算并比较了不同几何比尺和开孔率下开孔沉箱附近的涡量场、箱体外侧的波面时程曲线和波浪反射系数。结果表明,随着模型几何比尺的减小,开孔沉箱受到偏大的流体黏性力,致使更多波能在湍流运动中耗散,进而减小了波浪反射系数并降低了箱体外侧的波面高度。     

浮式两级缓冲防撞系统及浮体结构构型研究

以浮式两级缓冲防撞系统为研究对象,介绍其工作原理,利用动力学分析软件LS-DYNA构建了碰撞模型并进行了仿真计算,依据仿真结果完成了以阻尼索组与浮体结构构型为主要内容的防撞系统设计。研究确定了阻尼索组的布置位置、材料及数量,为浮体内部结构构型设计了7种比选方案,分别对5 000 t船舶撞击浮体前梁和侧梁的动态过程进行了仿真分析,对比了不同方案及不同板梁厚度下的接触力。研究结果表明:浮体内部构造若选用板梁厚度12 mm的方案6设计,则碰撞发生时浮体前梁和侧梁分别可以稳定地提供10.4 MN和7.9 MN的均值力,以及11.5 MN和8.9 MN的峰值力,设计浮体能够满足整体刚度要求,防撞系统符合双保护原则。防撞系统及浮体结构构型设计为下一步全面开展船舶撞击防撞系统后的动力学响应分析奠定了基础。     

Vertical structure of the quasi-two-dimensional velocity field of a viscous incompressible flow and the problem of nonlinear friction

An approximate theory is constructed to describe quasi-two-dimensional viscous incompressible flows. This theory takes into account a weak circulation in the vertical plane and the related divergence of the two-dimensional velocity field. The role of the nonlinear terms that are due to the interaction between the vortex and potential components of velocity and the possibility of taking into account the corresponding effects in the context of the concept of bottom friction are analyzed. It is shown that the nonlinear character of friction is a consequence of the three-dimensional character of flow, which results in the effective interaction of vortices with vertical and horizontal axes. An approximation of the effect of this interaction in quasi-two-dimensional equations is obtained with the use of the coefficient of nonlinear friction. The results based on this approximation are compared to the data of laboratory experiments on the excitation of a spatially periodic fluid flow.