Difference-frequency wave loads on a large body in multi-directional waves

The second-order difference-frequency wave forces on a large three-dimensional body in multi-directional waves are computed by the boundary integral equation method and the so-called FML formulation (assisting radiation potential method). Semi-analytic solutions for a bottom-mounted vertical circular cylinder are also developed to validate the numerical method. Difference-frequency wave loads on a bottom-mounted vertical cylinder and stationary four legs of the ISSC tension-leg platform (TLP) are presented for various combinations of incident wave frequencies and headings. These force quadratic transfer functions (QTF) can directly be used in studying slowly varying wave loads in irregular short-crested seas described by a particular directional spectrum. From our numerical results, it is seen that the slowly varying wave loads are in general very sensitive to the directional spreading function of the sea, and therefore wave directionality needs to be taken into account in relevant ocean engineering applications. It is also pointed out that the uni-directionality of the sea is not necessarily a conservative assumption when the second-order effects are concerned.     

波浪方向对两矩形驳船间窄缝内流体共振影响的数值研究

应用二阶非线性三维势流时域模型模拟两个并靠矩形驳船间的窄缝内流体共振现象。模型中,速度势分解为入射势和未知的散射势,其中散射势通过求解边界积分方程计算得到。采用四阶预报校正的方法对自由水面边界条件在时间上进行数值积分。本文研究了波浪方向对窄缝内一阶、二阶波高的影响。计算结果表明,入射波浪的角度不会影响窄缝内流体的共振频率,同时窄缝内流体发生共振时,最大的共振波高总是发生在窄缝中点处。然而,共振时对应的最大波高随着入射角度的变化而变化。同时,入射波浪的频率与窄缝内流体共振频率的相对大小导致窄缝内最大波高位置向上游或者下游移动。     

Numerical analysis on ship maneuvering coupled with ship motion in waves

This paper considers a numerical analysis of ship maneuvering performance in the presence of incident waves and resultant ship motion responses. To this end, a time-domain ship motion program is developed to solve the wave–body interaction problem with the ship slip speed and rotation, and it is coupled with a modular-type 4-DOF maneuvering problem. In this coupled problem, the second-order mean drift force, which can play an important role in the ship maneuvering trajectory, is estimated by using a direct pressure integration method. The developed method is validated by observing the second-order mean drift force, and planar trajectories in maneuvering tests with and without the presence of incident waves. The comparisons are made for two ship models, Series 60 with block coefficient 0.7 and the S-175 containership, with existing experimental data. The maneuvering tests observed in this study include a zig-zag test in calm water, and turning tests in calm water and in regular waves. The present results show a fair agreement of overall tendency in maneuvering trajectories.     

A novel decomposition of the quadratic transfer function (QTF) for the time-domain simulation of non-linear wave forces on floating bodies

In the present study, a novel method is proposed for the separation of the second-order sum- and difference-frequency wave forces—that is, quadratic transfer functions (QTFs)—on a floating body into three components due to wave–wave, wave–motion, and motion–motion action. By applying the new QTF components, the second-order wave forces on a floating body can be strictly computed in the time domain. In this work, the boundary value problems (BVPs) corresponding to the three kinds of QTF components were derived, and non-homogeneous boundary conditions on the free surface and the body surface were obtained. The second-order diffraction potentials were determined using the boundary integral equation method. In the solution procedure, the highly oscillatory and slowly converging integral on the free surface was evaluated in an accurate and effective manner. Furthermore, the application of the QTF components in the time domain was demonstrated. The second-order exciting forces in the time domain were divided into three parts. Each part of these forces was computed via a two-term Volterra series model based on the incident waves, the first-order motion response, and the QTF components. This method was applied to several numerical examples. The results demonstrated that this decomposition yields satisfactory results.     

Nonlinear wave loads on a vertical cylinder

Nonlinear contributions due to elevation of the free surface, the dynamic head, and the second-order velocity potential on the wave loads are presented in closed-form expressions. Such nonlinearities resulting from large-amplitude ocean waves are associated with the irrotational flow interacting with a fixed bottom-mounted vertical cylinder piercing the surface. These are expressed in the form of dynamic, waterline and quadratic forces all of which depend on the square of the wave amplitude. The appropriate modifications are made to both the classical Morison equation and the well-known linear diffraction theory of MacCamy and Fuchs for accounting the second-order effects.A limited comparative study is performed to verify the present theoretical derivations. In general, satisfactory agreements have been obtained with the test results from various laboratory studies by different researchers. However, under certain environmental conditions, some discrepancies still exist with the measured results.     

A numerical study of second-order springing of an elastic body using higher-order boundary element method (HOBEM)

Second-order springing on an elastic body with forward speed is analyzed by numerical simulations. The boundary-value problem for the velocity potential is solved by means of the direct time-domain higher-order boundary element method (HOBEM). The free-surface boundary condition in the boundary-value problem is approximated on the mean surface up to second order by use of perturbation and Taylor-series expansion methods. The body boundary condition for an elastic body is derived with various quantities which are redefined in the generalized mode. These variables such as mode shape, normal vector, etc. are obtained by using directional derivative and continuum mechanics, and the same mathematical expressions are used to obtain several second-order generalized forces. To validate the numerical results, the second-order hydrodynamic force on the bottom-mounted rigid/elastic cylinders without forward speed is compared with other semi-analytic results. The property of second-order forces on an elastic ship is studied by changing the flexural rigidity and forward speed with elastic response. It is confirmed that the second-order velocity potential is important for a body with forward speed and investigation should be made more on numerical methods for accurate computation of the second-order velocity-potential force with forward speed.     

层化海洋中二阶多色波对可渗透结构的绕射问题

可渗透结构具有使波浪作用减弱的效应,而海水的层化及水波的非线性使结构的波绕射产生多层复杂机制。将可渗透结构应用于复杂海况条件中,海水的层化性、波浪的非线性及结构的透空性构成了波绕射的一个十分复杂的数学问题。该问题存在理论研究的必要性,而文章则着重探讨其数学分析的可能性。通过引入二层海的层化海模式及Stokes二阶波的非线性波模式,给出了二阶多色波对透空结构的波绕射的定解问题提法,提出了复合形式的二阶多色波辐射条件式及可渗透结构的二阶物面条件式,应用特征函数解法与积分法推导了多色波对结构绕射的一阶势解与二阶作用的耦合积分解式,并讨论了解式所涉及无穷积分的算法。     

A method for second-order diffraction potential from an axisymmetric body

The paper provides a detailed analysis for the second-order diffraction of monochromatic waves. For the second-order potential on the free surface, the paper proposed a forward prediction method for computing the integration on the free surface. By this method we only need to run the infinity integration on the free surface directly for a few points; a one-step quadrature can then be applied successively outward from the body for potentials at other points. For wave diffraction from a body of revolution with a vertical axis, the paper derives a new integral equation, which can cancel the leading singularity in the derivative of ring Green's functions automatically. To obtain accurate results, different approaches are also used to deal with singularities in the ring Green's functions in the integration on both the body surface and free surface. The method has been implemented for bodies of revolution with vertical axes, but the theory is also available for arbitrary bodies.A numerical examination is made to validate the numerical code by comparing the second-order forces and moments on uniform and truncated cylinders and second-order diffraction potentials on the free surface with some published results. The comparison shows that the present results are in good agreement with those published. The method is also used to compute the second-order wave elevation around uniform and truncated cylinders.     

Second-order vertical and horizontal wave forces on a circular dock

A second-order potential solution is presented for the diffraction of a nonlinear progressive wave in finite-depth water, incident on a fixed circular dock. The usual perturbation analysis is used to produce first- and second-order subproblems. The mathematical method is based on the assumption that inner and outer solutions exist and these are matched by the requirements of continuity for mass flux and pressure between adjacent regions. It is shown that the solutions for the second-order problem can be derived in the same manner as in the first-order theory.     

Analytical method of radiation by a wave energy device with dual rectangular floating bodies

The system with one floating rectangular body on the free surface and one submerged rectangular body has been applied to a wave energy conversion device in water of finite depth. The radiation problem by this device on a plane incident wave is solved by the use of an eigenfunction expansion method, and a new analytical expression for the radiation velocity potential is obtained. The wave excitation force is calculated via the known incident wave potential and the radiation potential with a theorem of Haskind employed. To verify the correctness of this method, an example is computed respectively through the bound element method and analytical method. Results show that two numerical methods. are in good agreement, which shows that the present method is applicable. In addition, the trends of hydrodynamic coefficients and wave force are analyzed under different conditions by use of the present analytical method.     

Nonlinear hydroelastic analysis of a moored floating body

By integration of the second-order fluid pressure over the instantaneous wetted surface, the generalized first- and second-order fluid forces used in nonlinear hydroelastic analysis are obtained. The expressions for coefficients of the generalized first- and second-order hydrodynamic forces in irregular waves are also given. The coefficients of the restoring forces of a mooring system acting on a flexible floating body are presented. The linear and nonlinear three-dimensional hydroelastic equations of motion of a moored floating body in frequency domain are established. These equations include the second-order forces, induced by the rigid body rotations of large amplitudes in high waves, the variation of the instantaneous wetted surface and the coupling of the first order wave potentials. The first-order and second-order principal coordinates of the hydrelastic vibration of a moored floating body are calculated. The frequency characteristics of the principal coordinates are discussed. The numerical results indicate that the rigid resonance and the coupling resonance of a moored floating body can occur in low frequency domain while the flexible resonance can occur in high frequency domain. The hydroelastic responses of a moored box-type barge are also given in this paper. The effects of the second-order forces on the modes are investigated in detail.     

非线性载荷对船中弯矩响应影响

给出舰船结构水弹性分析的通用公式,分析典型高速船的水弹性响应,将总响应分为线性和非线性,前者用线性势流理论计算,后者用船体—流体系统脉冲响应函数的卷积积分和水动力非线性修正得到。计算表明线性解答过高地估计了共振点刚体运动,从而导致弯矩预报不准确,而非线性模拟的中拱/垂弯矩均与实验吻合较好。     

An analysis of second-order wave forces on floating bodies by using a higher-order boundary element method

The wave diffraction and radiation problem is studied numerically by using a higher-order boundary element method. The convergence of the higher-order boundary element method is tested systematically for bodies of different shapes. For the second-order force, particular attention is paid to the contribution of the second-order potential, following the line of Molin's approach. For numerical evaluation, the free surface is divided into three subregions; inner, intermediate and outer ones. In the inner region, the integral is evaluated numerically by using higher-order boundary elements. In the intermediate region, semi-analytic form is constructed with the help of eigen functions. In the outer domain, the analytic solution is available. This subdivision scheme reduces the numerical burden remarkably.     

Analysis of the extreme wave elevation due to second-order diffraction around a vertical cylinder

Statistical analysis of nonlinear random waves is important in coastal and ocean engineering. One approach for modeling nonlinear waves is second-order random wave theory, which involves sum- and difference-frequency interactions between wave components. The probability distribution of the non-Gaussian surface elevation can be solved using a technique developed by Kac and Siegert [21]. The wave field can be significantly modified by wave diffraction due to a structure, and the nonlinear diffracted wave elevation can be of interest in certain applications, such as the airgap prediction for an offshore structure. This paper investigates the wave statistics due to second-order diffraction, motivated by the scarcity of prior research. The crossing rate approach is used to evaluate the extreme wave elevation over a specified duration. The application is a bottom-supported cylindrical structure, for which semi-analytical solutions for the second-order transfer functions are available. A new efficient statistical method is developed to allow the distribution of the diffracted wave elevation to be obtained exactly, accounting for the statistical dependency between the linear, sum-frequency and difference-frequency components. Moreover, refinements are proposed to improve the efficiency for computing the free surface integral. The case study yields insights into the problem. In particular, the second-order nonlinearity is found to significantly amplify the extreme wave elevation, especially in the upstream region; conversely, the extreme elevation at an oblique location downstream is attenuated due to sheltering effects. The statistical dependency between the linear and sum-frequency components is also shown to be important for the extreme wave statistics.     

Application of the state space model to the system of wave energy conversion Analytical method for wave forces on single oscillating rectangular buoy

A new analytical method is proposed to analyze the force acting on a rectangular oscillating buoy due to linear waves. In the method a new analytical expression for the diffraction velocity potential is obtained first by use of the eigenfunction expansion method and then the wave excitation force is calculated by use of the known incident wave potential and the diffraction potential. Compared with the classical analytical method, it can be seen that the present method is simpler for a two-dimensional problem due to the comparable effort needed for the computation of diffraction potential and for that of radiated potential. To verify the correctness of the method, a classical example in the reference is recomputed and the obtained results are in good accordance with those by use of other methods, which shows that the present method is correct.     

Higher Harmonics Induced by Waves Propagating over A Submerged Obstacle in the Presence of Uniform Current

To investigate higher harmonics induced by a submerged obstacle in the presence of uniform current, a 2D fully nonlinear numerical wave flume（NWF） is developed by use of a time-domain higher-order boundary element method（HOBEM） based on potential flow theory. A four-point method is developed to decompose higher bound and free harmonic waves propagating upstream and downstream around the obstacle. The model predictions are in good agreement with the experimental data for free harmonics induced by a submerged horizontal cylinder in the absence of currents. This serves as a benchmark to reveal the current effects on higher harmonic waves. The peak value of non-dimensional second free harmonic amplitude is shifted upstream for the opposing current relative to that for zero current with the variation of current-free incident wave amplitude, and it is vice versa for the following current. The second-order analysis shows a resonant behavior which is related to the ratio of the cylinder diameter to the second bound mode wavelength over the cylinder. The second-order resonant position slightly downshifted for the opposing current and upshifted for the following current.     

状态空间模型在海洋波能转换系统中的应用 Ⅱ.单个矩形振荡浮子装置运动受力分析的解析方法

对线性入射波作用下的矩形振荡浮子运动受力分析给出了一种新的解析方法.对该方法首先采用特征函数展开法得到绕射速度势的一种新解析式,然后直接由绕射势和入射势来计算波浪激励力.与经典方法相比可以看出,对于立面二维问题,由于绕射势和辐射势的计算量相当,因此该方法比过去常用的解析法更为简单.为了验证方法的正确性,采用该方法对过去文献介绍的经典算例进行了计算,得到了与过去方法完全一致的结果,从而证明该方法是正确的.     

Fully nonlinear analysis of near-trapping phenomenon around an array of cylinders

The wave diffraction around an array of fixed vertical circular cylinders is simulated in a numerical wave tank by using a fully nonlinear model in the time domain. The emphasis of the paper lies in the insightful investigation of the nonlinear properties of the near-trapping phenomenon associated with the multiple cylinders. The numerical model is validated by analytical solutions as well as experimental data for waves propagating past two and four vertical cylinders in certain arrangements. An array of four identical circular cylinders at the corners of a square with an incident wave along the diagonal of the square is the main focus here for investigating the near-trapping phenomenon. When near-trapping occurs, the present study shows that an extremely high wave elevation near the cylinders can be observed. At the same time, the hydrodynamic forces on different cylinders are found to be either in phase or out of phase, leading to some characteristic force patterns acting on the whole structure. Due to the nature of the numerical model adopted, nonlinearity at different orders can be captured using a harmonic analysis. In addition to first- and second-order near-trapping, the third-order (triple-frequency) nonlinear component is presented for the first time. For the configuration selected, it is found that at one specific incident wave frequency and direction one trapped mode is excited by second-order effects, while a different trapped mode (having similar symmetries) is excited by the third harmonic of the incident wave frequency.     

Experimental investigation of the load exerted by nonstationary internal solitary waves on a submerged slender body over a slope

Laboratory experiments are performed in a large stratified fluid flume to examine the characteristics of the load on a submerged slender body that is exerted by a nonstationary internal solitary wave (ISW) from its interaction with a gentle slope. The nonstationary ISW over the slope and its load on the body are measured by using multi-channel conductivity-probe arrays and a specially designed force measurement device, respectively, and the body’s vertical and horizontal positions on the load are determined by analyzing the effects of the incident ISW’s amplitude. The experimental results show that the load on the slender body increases as the incident ISW’s amplitude increases; additionally, the effect of oscillations is enhanced because of the ISW’s distortion, breaking and fission. The oscillating action from fission waves becomes dominant as the amplitude reaches a certain value. Additionally, the load is correlated with body’s vertical position relative to the pycnocline. The magnitudes of the vertical and horizontal forces reach a maximum and minimum in the pycnocline, respectively, and the horizontal force in this direction is the opposite above and below the pycnocline. Compared to a case without a slope, the load on the slender body increases because of the nonstationary ISW, and its effect on the maximum force is transferred to the pycnocline. When the body’s horizontal position is located close to the top of the slope, the direction of the horizontal and vertical forces remains consistent, but its acting time becomes longer. In addition, high-frequency actions on the slender body are impacted by nonstationary ISWs near the slope’s top.     

潜堤后高阶自由谐波的研究

基于高阶边界元方法的完全非线性数值水槽模型模拟潜堤地形上波浪的传播变形,通过与实验值进行比较,考察数学模型的正确性.采用两点法分离得到堤后高倍频自由波来研究入射波参数、水深对堤后高倍频自由波的影响.研究发现:基频波、二阶和三阶自由波幅值分别与入射波波幅成线性、二次和三次函数关系,基频波幅值基本不随波浪周期变化,而二阶和...