A NURBS-based high-order panel method for three-dimensional radiation and diffraction problems with forward speed

A high-order Rankine panel method based on Non-Uniform Rational B-Spline (NURBS) is developed for solving the three-dimensional radiation and diffraction problems with forward speed. A NURBS surface is used to precisely represent the body geometry. Velocity potential on the body surface is described by B-spline after the source density distribution on the boundary surface is determined. A collocation approach is applied to numerical computation and the integral equations are evaluated by applying Gauss–Legendre quadrature. The m_j-terms are evaluated by a desingularized method which utilizes NURBS technique. In order to verify the method proposed, it is firstly applied to the unbounded flow problem of a sphere and spheroids. The numerical results are found to be in good agreement with analytical solutions. Then the method is used to solve the radiation and diffraction problems of a sphere and the diffraction problem of a spheroid moving with a forward speed beneath the free surface in frequency domain. The numerical results are satisfactory in comparison with the published analytical results and experimental results.     

A B-spline-based method for radiation and diffraction problems

An open uniform B-spline-based panel method is developed for solution of potential flow problems. In this method, both geometry as well as the field variables are represented by the same open uniform B-spline basis function. The method is initially applied for the radiation problem in unbounded fluid. Computed results for a spheroid of different aspect ratio are found to be in excellent agreement with analytical results. The method is then applied for diffraction problem formulated based on the transient (time-domain) Green's function. Computed results for a hemisphere and Wigley hull are compared with published results and the comparison shows good agreement.     

On the evaluation of time-domain Green function

An analytical method has been developed to evaluate the wave part of the time-domain Green function and its derivatives. Based on Taylor series expansion, the Green function is obtained by solving a fourth-order ordinary differential equation. The method accelerates the convergence of the summation of an infinite series in the numerical computation. The accuracy of this method was demonstrated by its comparison with other method and its application to solve the radiation problem of a floating hemisphere using a panel-free method. The computed hydrodynamic coefficients agree well with the analytical solutions.     

Numerical study on added resistance of ships by using a time-domain Rankine panel method

This paper considers the numerical computation of added resistance on ships. As a method of solution, a higher-order Rankine panel method is applied in time domain. Furthermore, two different formulations, Neumann-Kelvin and double-body linearizations, are applied to investigate the difference between computational results. The quantities of added resistance are obtained by integrating the second-order pressure on the body surface. Computational results are validated by comparing with experimental data and other computational results on several models, and fair agreements are found for all the models. The components of added resistance are observed with regard to integral terms, radiation and diffraction components. This study continues to the computation of added resistance in irregular waves. This kind of computational approach can provide some important information to understand the uncertainty level of towing-tank measurement of added resistance. In this study, a series of long numerical simulations are carried out, and their statistical variations are observed for different time windows. Based on this observation, an appropriate criterion and averaging method are proposed for the prediction of added resistance in irregular waves.     

Second-order Taylor expansion boundary element method for the second-order wave radiation problem

A novel Boundary Element Method (BEM) named the second-order Taylor Expansion Boundary Element Method (the 2nd order TEBEM) is developed for the solution of the second-order wave radiation velocity potential and sum-frequency wave loads for floating bodies. The radiation condition is enforced by a hybrid method of the multi-transmitting formula and damping zone. For the interior domain problem of a cube and a sphere, numerical results demonstrate that the 2nd order TEBEM can accurately solve the first and second-order gradients of velocity potential on the no-smoothed and smoothed boundary compared to the low-order BEM. The double frequency forces acting on the truncated cylinder are calculated under finite water depth. The agreement between the 2nd order TEBEM and others' numerical results is good. Moreover, all of the singular integrals in the 2nd order TEBEM can be solved analytically, so its implementation is much easier compared to the high-order BEM.     

预修正快速傅里叶变换高阶边界元方法在波浪对物体作用问题中的应用

边界元方法被广泛应用于波浪对海上婕筑物作用领域,但由于传统边界元方法的存储量和计算量均为未知量的平方量级,很难满足大范围多未知量问题的计算需要.本文基于高阶边界元方法,应用预修正快速傅里叶变换方法,使计算量与存储量分别降低至O(NlogN)和O(N)量级,并得到一个连续的压强分布以适应结构设计的要求,同时可以通过使用满...     

海洋缆体系统的统一凝集参数时域分析法

提出了海洋缆一体系统的三维动态性能时域分析法，并给出了计算机程序。采用凝集参数法模拟了缆一体结构系统，缆为圆形截面，可伸长，不考虑弯曲刚度，采用了四阶龙格一库塔法积分系统运动方程，计算了二个实例以验证程序的可靠性，程序可用于系泊，拖曳，潜器脐带缆等具有组合成分的海洋缆体系统的初步设计。     

Diffraction–radiation of multiple floating structures in directional waves

The dynamics of multiple floating structures have been studied using the finite element method. The emphasis is on the hydrodynamic behaviour of multiple bodies under a multi-directional wave field. A two-dimensional numerical model has been adopted to evaluate hydrodynamic coefficients and forces in an oblique wave field. The responses in sway, heave and roll modes are reported. The linear filter technique is then used to extrapolate the responses under directional waves. The effect of mean wave direction and directional homogeneity on the hydrodynamic behaviour of the structure is studied. Based on the present study, it is found that the two-dimensional model is applicable to investigate the wave-structure interaction problems of the type herein considered.     

Statistical analysis of wave-induced extreme nonlinear load effects using time-domain simulations

A new hybrid method for the time-domain nonlinear simulation of the hydroelastic load effects and the peak-over-threshold (POT) method for the calculation of the short-term extreme responses are briefly described and applied to a flexible containership of the latest design. Statistical analysis has been carried out to study the sensitivity of the predicted extreme vertical bending moments and vertical shear forces to the changes in the threshold of the POT method, as well as the statistical uncertainty in the prediction due to the limited duration of the nonlinear simulation. It is recommended that 90%–95% quantile should be used as the threshold in the POT method and more than 100 h of time-domain simulation should be carried out in order to obtain satisfactory predictions of the short-term extreme nonlinear load effects.     

A wavelet for predicting the time-domain response of vertically tethered systems

The response wavelets for tether tension and the displacement of a package on the end of a tether are derived from a linear, frequency-domain model of vertically tethered systems using inverse fourier transform techniques. These wavelets are convolved with a record of ship displacement to predict the instantaneous tension and displacement to an accuracy of better than 20% for a typical system. Although snap loads are non-linear and the response cannot be calculated while the tether is slack, the onset of zero tension — a precursor to snap loads — is predicted by the wavelet. Thus, extreme and potentially hazardous conditions are foreseeable with the wavelet and this information can be used in real-time to guide operations at sea.     

浮基多体系统自激运动响应的时域分析方法

浮基多体系统上部机构作业时不仅存在多体系统内部各物体间的耦合作用 ,还存在浮基与流场间的耦合作用。本文将浮基在流场中引起的辐射势表示为浮基加速度与时域规范化速度势乘积的卷积形式 ,导出流体动压力及静水系泊恢复力与浮基运动的关系。使用多刚体力学的凯恩方法得到系统的动力学方程。最后给出了在时域中数值求解系统运动响应的具体步骤     

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.     

Numerical analysis and validation of weakly nonlinear ship motions and structural loads on a modern containership

This paper aims to validate a numerical seakeeping code based on a 3D Rankine panel method by comparing its results with experimental data. Particularly, the motion response and hull-girder loads on a real modern ship, a 6500 TEU containership, are considered in this validation study. The method of solution is a 3D Rankine panel method which adopts B-spline basis function in the time domain. The numerical code is based on the weakly nonlinear scheme which considers nonlinear Froude-Krylov and restoring forces. The main focus of this study is given to investigate the nonlinear characteristics of wave-induced loads, and to validate this present scheme for industrial use in the range of low Froude number. The comparisons show that the nonlinear motions and hull-girder loads, computed by the present numerical code, have good overall agreements with experimental results. It is found that, for the better accuracy of computational results, particularly at extreme waves in oblique seas, the careful treatment of soft-spring (or compatible) system is recommended to the control of non-restoring motions such as surge, sway, and yaw.     

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.     

Solution of the nonlinear roll model by a generalized asymptotic method

Mathematical modeling of the nonlinear roll motion of ships is one subject widely dealt with in nonlinear ship dynamics. This paper investigates setting up a form of nonlinear roll motion model and developing its periodic solution by the generalized Krylov–Bogoliubov asymptotic method in the time domain. In this model, nonlinearities are introduced through damping and restoring terms. The restoring term is approximated as a third-order odd polynomial whereas the quadratic term is favored to represent the nonlinear damping. The ship is assumed to be under the influence of a sinusoidal exciting force. Although the method is expressible to contain any order of the perturbing term, a single degree is chosen to avoid cumbersome mathematical complexity. In order to improve the solution a first-order correction term is also included. Moreover, a numerical example is carried out for a small vessel in order to validate the solution scheme.     

非线性波浪时域计算的三维耦合模型

将计算区域Ω划分为内域Ω₁和外域Ω₂(Ω₂=Ω-Ω₁),外域控制方程采用改进线性频散特性的二维Boussinesq方程,用预报一校正法数值求解;结构物附近的内域控制方程为三维Navier-Stokes方程,由VOF方法数值求解。通过在外域和内域相匹配的交界面上设置合适的速度和波面边界条件,建立了三维非线性波浪时域计算的耦合模型。模拟试验表明:(1)耦合模型数值波浪水池可以产生稳定的、重复性较好的波动过程;(2)用耦合模型数值波浪水池求解较大浅水区域上的非线性波浪数值计算问题可以取得较高的计算效率,同时又能得出结构物附近的复杂流场。     

Hybrid frequency–time domain models for dynamic response analysis of marine structures

Time-domain models of marine structures based on frequency domain data are usually built upon the Cummins equation. This type of model is a vector integro-differential equation which involves convolution terms. These convolution terms are not convenient for analysis and design of motion control systems. In addition, these models are not efficient with respect to simulation time, and ease of implementation in standard simulation packages. For these reasons, different methods have been proposed in the literature as approximate alternative representations of the convolutions. Because the convolution is a linear operation, different approaches can be followed to obtain an approximately equivalent linear system in the form of either transfer function or state-space models. This process involves the use of system identification, and several options are available depending on how the identification problem is posed. This raises the question whether one method is better than the others. This paper therefore has three objectives. The first objective is to revisit some of the methods for replacing the convolutions, which have been reported in different areas of analysis of marine systems: hydrodynamics, wave energy conversion, and motion control systems. The second objective is to compare the different methods in terms of complexity and performance. For this purpose, a model for the response in the vertical plane of a modern containership is considered. The third objective is to describe the implementation of the resulting model in the standard simulation environment Matlab/Simulink.     

Parametric interpolation method for probabilistic fatigue analysis of steel risers

Slender marine structures such as mooring lines and risers are susceptible to failures due to stress variations coming from environmental actions. Wave, wind and current are random phenomena, and consequently the most adequate methodology to estimate the fatigue damage accumulation on these structures is the probabilistic fatigue analysis. In practice, the estimation of fatigue life requires the numerical simulation of a huge number of loading cases to compute the multi-dimensional integral of the total fatigue damage.This paper presents an efficient approach to compute the long-term fatigue damage of marine structures. The proposed method needs only a few number of numerical simulations to estimate the structure fatigue life. It uses a parametric interpolation procedure to evaluate the fatigue damage for any individual short-term condition (sea state) required in the calculation of the multi-dimensional integral. In this way, the total number of short-term structural analyzes is considerably reduced.The effectiveness and accuracy of the proposed method is compared to the full direct integration by means of two comprehensive examples. The first studied case is an idealized theoretical model of a single-degree-of-freedom (SDOF) system under wave loading, and the second one is a Steel Catenary Riser (SCR) connected to a FPSO (Floating, Production, Storage and Offloading floating unit).     

Domain decomposition and matching for time-domain analysis of motions of ships advancing in head sea

A domain decomposition and matching method in the time-domain is outlined for simulating the motions of ships advancing in waves. The flow field is decomposed into inner and outer domains by an imaginary control surface, and the Rankine source method is applied to the inner domain while the transient Green function method is used in the outer domain. Two initial boundary value problems are matched on the control surface. The corresponding numerical codes are developed, and the added masses, wave exciting forces and ship motions advancing in head sea for Series 60 ship and S175 containership, are presented and verified. A good agreement has been obtained when the numerical results are compared with the experimental data and other references. It shows that the present method is more efficient because of the panel discretization only in the inner domain during the numerical calculation, and good numerical stability is proved to avoid divergence problem regarding ships with flare.     

考虑波-管-土耦合作用的海底管道在位稳定性分析方法

以水动力加载试验得到的描述管道在位稳定性的波浪环境参数、管道参数和海床土性参数之间的无量纲耦合关系为基础,提出了一种考虑波-管-土动力耦合作用的海底管道在位稳定性分析方法,并通过算例与挪威船级社推荐的DnV设计规范进行了比较分析。分析表明,考虑波-管-土耦合作用的管道稳定性分析方法与DnV管道设计规范有很好的可比性．而且物理机理更清晰,可为管道稳定性设计提供有益的参考。