Finite element analysis of second order wave radiation by a group of cylinders in the time domain

A finite element based numerical method is employed to analyze the wave radiation by multiple or a group of cylinders in the time domain. The nonlinear free surface and body surface boundary conditions are satisfied based on the perturbation method up to the second order. The first- and second-order velocity potential problems at each time step are solved through a Finite Element Method （FEM）. The matrix equation of the FEM is solved through iteration and the initial solution is obtained from the result at the previous time step. The three-dimensional （3-D） mesh required is generated based on a two-dimensional （2-D） hybrid mesh on a horizontal plane and its extension in the vertical direction. The hybrid mesh is generated by combining an unstructured grid away from cylinders and two structured grids near the cylinder and the artificial boundary. The fluid velocity on the free surface and the cylinder surface are calculated by using a differential method. Results for various configurations including the cases of two cylinders and four cylinders and a group of eighteen cylinders are obtained to show the joint influences of cylinders on the first- and second- order waves and forces, including the effects of spacing ratios and wave frequency on the second order waves and the mean force, in particular.

A SECOND ORDER STATIONARY SOLUTION OF THE THREE DIMENSIONAL WAVE DIFFRACTION POTENTIAL

A theory on the second order wave diffraction by a three dimensional body fixed in a regular sea has been developed in the present paper. By regarding the sinusoidal disturb potential as a stationary solu- tion of an initial value problem, and using Laplace transformation method and Tauberian theorem, the boundary value problems of stationary solution of the first and second order diffraction potential have been de- rived in this paper. Furthermore, the explicit solution of the second order stationary diffraction potential has been obtained with the method of the double Fourier transformation. It is found that the asymptotic behaviour of the second order stationary solution at far field is dependent on two wave systems, the first is 'free wave', travelling independently of the first order wave system, the other is 'phase locked waves', which accompany the first order waves. At the same time, the radiation conditions of the second order diffraction problems are derived. We also find that one can still pursue a steady state formulation with the inclusion of Rayleigh damping. Finally, as an example, the second order wave forces upon a fixed vertical cir- cular cylinder have been calculated, and the numerical results agree well with the experimental data.

A second order random wave model for predicting the power performances of a wave energy converter

The power performances of a point absorber wave energy converter(WEC)operating in a nonlinear multidirectional random sea are rigorously investigated.The absorbed power of the WEC Power-Take-Off system has been predicted by incorporating a second order random wave model into a nonlinear dynamic filter.This is a new approach,and,as the second order random wave model can be utilized to accurately simulate the nonlinear waves in an irregular sea,avoids the inaccuracies resulting from using a first order linear wave model in the simulation process.The predicted results have been systematically analyzed and compared,and the advantages of using this new approach have been convincingly substantiated.

Effects of Second-Order Sum-and Difference-Frequency Wave Forces on the Motion Response of a Tension-Leg Platform Considering the Set-down Motion

This paper presents a study on the motion response of a tension-leg platform(TLP) under first-and second-order wave forces, including the mean-drift force, difference and sum-frequency forces. The second-order wave force is calculated using the full-field quadratic transfer function(QTF). The coupled effect of the horizontal motions, such as surge, sway and yaw motions, and the set-down motion are taken into consideration by the nonlinear restoring matrix. The time-domain analysis with 50-yr random sea state is performed. A comparison of the results of different case studies is made to assess the influence of second-order wave force on the motions of the platform. The analysis shows that the second-order wave force has a major impact on motions of the TLP. The second-order difference-frequency wave force has an obvious influence on the low-frequency motions of surge and sway, and also will induce a large set-down motion which is an important part of heave motion. Besides, the second-order sum-frequency force will induce a set of high-frequency motions of roll and pitch. However, little influence of second-order wave force is found on the yaw motion.

SECOND-ORDER INTERACTION OF IRREGULAR WAVES WITH A TRUNCATED COLUMN

A complete semi-analytical solution is obtained for second-order diffraction of plane bichromatic waves by a fixed truncated circular column.The fluid domain is divided into interior and exterior regions.In the exterior region,the second-order velocity potential is expressed in terms of‘locked-wave’and‘free-wave’ components,both are solved using Fourier and eigenfunction expansions.The re- sulting‘locked wave’potential is expressed by one-dimensional Green's integrals with oscillating integrands.In order to increase computational efficiency,the far-field part of the integrals are carried out analytically.Solutions in both regions are matched on the interface by the potential and its normal derivative continuity conditions.Based on the present approach,the sum-and difference-frequency potentials are efficiently evaluated and are used to generate the quadratic transfer functions which correlates the incident wave spectrum with second-order forcing spectrum on the column.The sum-frequency QTFs for a TLP column are present,which are compared for some frequency pairs with those from a fully numerical procedure.Satisfactory agreement has been obtained.QTF spectra for a case study TLP column,generated using the semi-analytical solution are presented.Also given are the results for nonlinear wave field around the column.

Second-order random wave solutions for interfacial internal waves in N-layer density-stratified fluid

This paper studies the random internal wave equations describing the density interface displacements and the velocity potentials of N-layer stratified fluid contained between two rigid walls at the top and bottom. The density interface displacements and the velocity potentials were solved to the second-order by an expansion approach used by Longuet-Higgins （1963） and Dean （1979） in the study of random surface waves and by Song （2004） in the study of second- order random wave solutions for internal waves in a two-layer fluid. The obtained results indicate that the first-order solutions are a linear superposition of many wave components with different amplitudes, wave numbers and frequencies, and that the amplitudes of first-order wave components with the same wave numbers and frequencies between the adjacent density interfaces are modulated by each other. They also show that the second-order solutions consist of two parts： the first one is the first-order solutions, and the second one is the solutions of the second-order asymptotic equations, which describe the second-order nonlinear modification and the second-order wave-wave interactions not only among the wave components on same density interfaces but also among the wave components between the adjacent density interfaces. Both the first-order and second-order solutions depend on the density and depth of each layer. It is also deduced that the results of the present work include those derived by Song （2004） for second-order random wave solutions for internal waves in a two-layer fluid as a particular case.

Effects of Second Order Scattered Wave in Stratified Ocean

By using the theories on Stokes multicolored water waves and taking the two- layer ocean as a basic model of stratified ocean, the paper analyzes the problems related to the effects of the nonlinear water wave on offshore structures. A mathematical expression is presented to describe second order wave radiation conditions. Using integral principle, the analytical integral solutions are given to evaluate second order scattered wave loads on general vertical circular cylinders in the two-layer ocean, and the special recurrence formulas for infinite integrals over free and stratified surfaces are derived.

Effects of Second-Order Difference-Frequency Wave Forces on Floating Wind Turbine Under Survival Condition

In this paper, the effects of second-order difference-frequency wave forces on the global motion of an offshore wind turbine system with a large displacement under the survival condition are studied. In this case, the hydrodynamic force is the main force because the blades are feathered to reduce the lifting force. The first-order hydrodynamic forces are calculated by WADAM, while the second-order wave forces are calculated by a customized MATLAB module. Then the hydrodynamic coefficients are transferred to the wind turbine analytical code FAST. Through the comparisons of dynamic responses between the first- and second-order numerical models, it is found that the second-order wave forces significantly influence the motion of floating wind turbine under the survival condition. Moreover, neglecting the second-order force significantly underestimates the tension forces in the mooring lines. © 2017, Tianjin University and Springer-Verlag Berlin Heidelberg.

Second-Order Solutions for Random Interfacial Waves Under Steady Uniform Currents

In the present research, the study of Song (2004) for random interfacial waves in two-layer fluid is extended to the case of fluids moving at different steady uniform speeds. The equations describing the random displacements of the density interface and the associated velocity potentials in two-layer fluid are solved to the second order, and the wave-wave interactions of the wave components and the interactions between the waves and currents are described. As expected, the extended solutions include those obtained by Song (2004) as one special case where the steady uniform currents of the two fluids are taken as zero, and the solutions reduce to those derived by Sharma and Dean (1979) for random surface waves if the density of the upper fluid and the current of the lower fluid are both taken as zero.

Second-order random interfacial wave solutions for two-layer fluid with a free surface

A previous study （Song. 2004. Geophys Res Lett, 31 （15）：L15302） of the second-order solutions for random interracial waves is extended in a constant depth, two-layer fluid system with a rigid lid is extended into a more general case of two-layer fluid with a top free surface. The rigid boundary condition on the upper surface is replaced by the kinematical and dynamical boundary conditions of a free surface, and the equations describing the random displacements of free surface, density-interface and the associated velocity potentials in the two-layer fluid are solved to the second order using the same expansion technology as that of Song （2004. Geophys Res Lett, 31 （15）：L15302）. The results show that the interface and the surface will oscillate synchronously, and the wave fields to the first-order both at the free surface and at the density-interface are made up of a linear superposition of many waves with different amplitudes, wave numbers and frequencies. The second-order solutions describe the second-order wave-wave interactions of the surface wave components, the interface wave components and among the surface and the interface wave components. The extended solutions also include special cases obtained by Thorpe for progressive interracial waves （Thorpe. 1968a.Trans R Soc London, 263A：563～614） and standing interracial waves （Thorpe. 1968b. J Fluid Mech, 32：489-528） for the two-layer fluid with a top free surface. Moreover, the solutions reduce to those derived for random surface waves by Sharma and Dean （1979.Ocean Engineering Rep 20） if the density of the upper layer is much smaller than that of the lower layer.

Establishment of Numerical Wave Flume Based on the Second-Order Wave-Maker Theory

With growing computational power, the first-order wave-maker theory has become well established and is widely used for numerical wave flumes. However, existing numerical models based on the first-order wave-maker theory lose accuracy as nonlinear effects become prominent. Because spurious harmonic waves and primary waves have different propagation velocities, waves simulated by using the first-order wave-maker theory have an unstable wave profile. In this paper, a numerical wave flume with a piston-type wave-maker based on the second-order wave-maker theory has been established. Dynamic mesh technique was developed. The boundary treatment for irregular wave simulation was specially dealt with. Comparisons of the free-surface elevations using the first-order and second-order wave-maker theory prove that second-order wave-maker theory can generate stable wave profiles in both the spatial and time domains. Harmonic analysis and spectral analysis were used to prove the superiority of the second-order wave-maker theory from other two aspects. To simulate irregular waves, the numerical flume was improved to solve the problem of the water depth variation due to low-frequency motion of the wave board. In summary, the new numerical flume using the second-order wave-maker theory can guarantee the accuracy of waves by adding an extra motion of the wave board. The boundary treatment method can provide a reference for the improvement of nonlinear numerical flume.

P-and S-wavefield simulations using both the firstand second-order separated wave equations through a high-order staggered grid finite-difference method

In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this, especially when the velocity field is complex. A useful approach in multi-component analysis and modeling is to directly solve the elastic wave equations for the pure P- or S-wavefields, referred as the separate elastic wave equa- tions. In this study, we compare two kinds of such wave equations： the first-order （velocity-stress） and the second- order （displacement-stress） separate elastic wave equa- tions, with the first-order （velocity-stress） and the second- order （displacement-stress） full （or mixed） elastic wave equations using a high-order staggered grid finite-differ- ence method. Comparisons are given of wavefield snap- shots, common-source gather seismic sections, and individual synthetic seismogram. The simulation tests show that equivalent results can be obtained, regardless of whether the first-order or second-order separate elastic wave equations are used for obtaining the pure P- or S-wavefield. The stacked pure P- and S-wavefields are equal to the mixed wave fields calculated using the corre- sponding first-order or second-order full elastic wave equations. These mixed equations are computationallyslightly less expensive than solving the separate equations. The attraction of the separate equations is that they achieve separated P- and S-wavefields which can be used to test the efficacy of wave decomposition procedures in multi-com- ponent processing. The second-order separate elastic wave equations are a good choice because they offer information on the pure P-wave or S-wave displacements.

The second order spectrum of two-dimensional random waves

On the basis of the second order solution of two-dimensional random gravity waves in finite uniform depth,which is rederived by a perturbation expansion method,the analytical expression of the second order spectrum is strictly deduced, and for infinite depth, the correct form of the kernel function is given for each octant do main.In fact,the present study improves and corrects the generally accepted results obtained by Tick

Statistical Distribution of Depth-Integrated Local Horizontal Momentum for Second-Order Random Ocean Waves in Finite Water Depth

Based on the second order random wave solutions of water wave equations in finite water depth, statistical distributions of the depth integrated local horizontal momentum components are derived by use of the characteristic function expansion method. The parameters involved in the distributions can be all determined by the water depth and the wave number spectrum of ocean waves. As an illustrative example, a fully developed wind generated sea is considered and the parameters are calculated for typical wind speeds and water depths by means of the Donelan and Pierson spectrum. The effects of nonlinearity and water depth on the distributions are also investigated.

AN ANALYTICAL SOLUTION OF SECOND－ORDER WAVE FORCE ON A VERTICAL CIRCULAR CYLINDER

So far many investigations have been made on nonlinear wave diffraction problemfor a large-diameter vertical circular cylinder. However,there are still some problemsworthy to be further discussed. It includes that the second order radiation condition isnot very clear and the inhomogeneous term of second order free surface boundarycondition makes the calculation of second order wave force either not easier to beperformed accurately due to its slowly decaying with radial distance or toocomplicated for practical application. In this paper, the second order radiationcondition is posed of the circumferential Fourier components of second orderpotential, instead of the second order potential. It is found that the circumferenatialFourier cormponents of second order potential have to satisfy Sommerfeld radiationcondition. By means of the mathematical formulae derived in this paper, theinhomogeneous term of second order free surface boundary condition were simplifiedand then an exact expression of second order wave force was obtained, which issimpler in form and easier to be used in practical calculation. The calculation resultsagree well with some experimental data.

Second-Order Potentials and Forces for Two-Dimensional Diffraction Problem of Finite Water Depth

-In this paper, an analytical solution in the outer region of finite water depth is derived for the second-order diffraction potential, which gives a clear physical meaning of the wave transmission and reflection characteristics in the far field. A numerical method-simple Green's function technique-for calculating the second-order diffraction potential in the inner region is also described. Numerical results are provided for the second-order wave forces on a semi-submerged cylinder. It is found that the contribution of second-order diffraction potential to second-order wave forces is important. The effect of water depth and submerged depth on the wave force is also discussed.

浮式风力机若干特征动力学问题综述

总结浮式风力机类型及其对应的特征动力学问题,针对浮式风力机气动荷载、水动荷载的计算方法以及结构动力学、控制动力学典型问题进行论述。讨论了气动—水动—结构—伺服耦合分析的难点,重点分析了二阶波浪力、畸形波等非线性波浪荷载、流荷载及涡激运动对浮式风力机特征动力响应的影响。阐述了浮式风力机动力学研究的试验方法、数值仿真方法、样机测试方法,并对模型试验技术的相似理论、气动模型的实现和难点以及数值仿真的频域方法、时域方法和分析工具进行了归纳对比。研究表明:浮式风力机多场、多体耦合动力分析机理及相关技术仍不成熟,气动荷载、高阶非线性波浪荷载耦合模型的建立是动力学问题研究的重点,数值仿真及模型试验是浮式风力机动力响应研究的主要方法,样机测试技术的积累将促进设计标准的完善及浮式风电的产业化发展。

二阶Stokes波作用下多孔弹性海床动力响应解析解

基于二阶Stokes波理论,研究了非线性海洋波作用下多孔弹性海床的动力响应问题。利用Biot固结理论以及Verruijt储存方程,在笛卡尔坐标下建立了海床的控制方程。通过将海洋波函数拓展到复数域,采用严格的数学推导求得了海床动力响应的解析解,并将该解与既有解进行对比验证。最后分析了波浪和海床特征参数对海床有效应力、剪切应力和孔压分布的影响,结果表明波浪和海床的参数对海床动力响应具有显著的影响。其中波浪周期和水深对二阶Stokes波的二阶项有较大影响,渗透系数和剪切模量则分别影响海床动力响应的变化速率和幅值。

二阶波浪力对浮式风机系统动力响应的影响研究

文中以OC4-DeepCwind海上浮式风机系统为例,考虑浮式基础运动与气动力之间的耦合,研究了二阶波浪力对浮式风机系统动力响应的影响.结果表明:二阶波浪力增加了浮式基础纵荡和纵摇运动的标准差,而浮式基础纵荡和纵摇对风机气动力有较大的影响,从而使得风轮的推力及功率波动都增大.因此,在风力机系统控制设计及风力机叶片结构设计时,需要充分考虑二阶波浪力的影响.此外,二阶和频力和塔柱一阶固有频率接近,激发了塔柱自振,在设计塔柱时应避开二阶和频力.

有限水深波浪作用下弱惯性塑料块体漂移规律

塑料漂浮垃圾对海洋环境影响巨大,在近岸的传播过程主要受波浪作用。针对塑料漂浮垃圾在近岸的运动特征,以往研究尚不深入。本文采用物理模型实验,对塑料漂浮物在有限水深波浪作用下的漂移规律开展研究,探讨了弱惯性塑料块体水平漂移速度与漂浮物性状以及波陡之间的关系。实验结果表明,塑料块体的漂移量受斯托克斯漂移与欧拉回流的共同影响,与二阶拉格朗日漂移理论吻合良好。在漂浮物尺寸远小于波长的情况下,漂浮物尺寸以及密度的改变对漂移没有显著的影响;漂浮物的漂移量与波陡的平方成正比,根据本文的实验及以往公开发表的实验数据,修正了经验公式,为塑料漂浮物在近岸的传播与预测提供了有益借鉴。