On the parametric rolling of ships using a numerical simulation method

This paper has shown a numerical motion simulation method which can be employed to study on parametric rolling of ships in a seaway. The method takes account of the main nonlinear terms in the rolling equation which stabilize parametric rolling, including the nonlinear shape of the righting arm curve, nonlinear damping and cross coupling among all 6 degrees of freedom. For the heave, pitch, sway and yaw motions, the method uses response amplitude operators determined by means of the strip method, whereas the roll and surge motions of the ship are simulated, using nonlinear motion equations coupled with the other 4 degrees of freedom. For computing righting arms in seaways, Grim's effective wave concept is used. Using these transfer functions of effective wave together with the heave and pitch transfer functions, the mean ship immersion, its trim and the effective regular wave height are computed for every time step during the simulation. The righting arm is interpolated from tables, computed before starting the simulation, depending on these three quantities and the heel angle. The nonlinear damping moment and the effect of bilge keels are also taken into account. The numerical simulation tool has shown to be able to model the basic mechanism of parametric rolling motions. Some main characteristics of parametric rolling of ships in a seaway can be good reproduced by means of the method. Comprehensive parametric analyses on parametric rolling amplitude in regular waves have been carried out, with that the complicated parametric rolling phenomena can be understood better.     

基于非线性数学模型的半潜式平台共振运动特性研究

共振运动是深海浮式平台设计的关键考虑因素之一,对海洋平台的作业具有重要影响。采用半潜式平台运动的非线性耦合数学模型,考虑浮筒和横撑出入水以及垂荡、横摇和纵摇运动耦合对平台浮力和恢复力的影响,研究半潜式平台非线性共振运动特性,以及不规则波浪参数对运动的影响。研究表明:在非线性耦合运动和浮力变化的影响下,半潜式平台纵摇和垂荡运动的固有周期会随运动幅值的增大而逐渐减小,且最终趋于稳定,对纵摇运动周期的影响更为显著;非线性效应会使半潜式平台产生显著的低频纵摇共振响应,以及共振频率漂移的现象,且受随机种子和波浪周期的影响较小。     

Experimental and Numerical Investigations of Ship Parametric Rolling in Regular Head Waves

Parametric rolling is one of five types of the ship stability failure modes as proposed by IMO. The periodic change of the metacentric height is often considered as the internal cause of this phenomenon. Parametric rolling is a complex nonlinear hydrodynamic problem, often accompanied by large amplitude vertical motions of ships. In recent years,the Reynolds-averaged Navier–Stokes(RANS) equation simulations for viscous flows have made great progress in the field of ship seakeeping. In this paper, the parametric rolling for the C11 containership in regular waves is studied both experimentally and numerically. In the experiments, parametric rolling amplitudes at different drafts, forward speeds and wave steepnesses are analyzed. The differences in the steady amplitudes of parametric rolling are observed for two drafts. The effect of the incident wave steepness(or wave amplitude) is also studied, and this supports previous results obtained on limits of the stability for parametric rolling. In numerical simulations, the ship motions of parametric rolling are analyzed by use of the potential-flow and viscous-flow methods. In the viscousflow method, the Reynolds-averaged Navier–Stokes equations are solved using the overset grid method. The numerical accuracies of the two methods at different wave steepnesses are also discussed.     

Numerical investigation of the seakeeping behavior of a catamaran advancing in regular head waves

A numerical study was undertaken in order to assess the capability of an unsteady RANS code to predict the seakeeping characteristics of a high-speed multi-hull vessel in high sea states. Numerical analysis includes evaluation of ship motions, effects of wave steepness on ship response, catamaran natural frequency and added resistance in waves. Computations were performed for the DELFT 372 catamaran by the URANS solver CFDSHIP-Iowa V.4. The code was validated with encouraging results for high ship speeds (0.3≤Fn≤0.75) and high wave amplitudes (0.025≤Ak≤0.1). Comparison with strip theory solutions shows that the RANS method predicts ship motions with higher accuracy and allows the detection of nonlinear effects. Current computations evidence that heave peaks occur at resonance for all Fn, and reach the absolute maximum at Fn=0.75. Maximum pitch occurs at frequencies lower than resonance, for each speed, and absolute maximum occurs at medium Fn=0.6. Maximum added resistance, R_aw, was computed at Fn=0.45, which, interestingly, is near the catamaran Fn_coincidence. Overall, we found similar results as Simonsen et al. (2008) for KCS containership, though, herein, a multi-hull geometry and higher speeds were tested. Also, our results are useful to further evaluate the exciting forces and their correlation with f_e and λ/L_pp.     

Validation of a time-domain strip method to calculate the motions and loads on a fast monohull

The paper presents a comparison between experimental data and numerical results of the hydrodynamic coefficients and also of the wave induced motions and loads on a fast monohull model. The model with 4.52 m length was constructed in Fibre Reinforced Plastic (FRP), and made up of 4 segments connected by a backbone in order to measure sectional loads. The objective of the investigation was to assess the capability of a nonlinear time domain strip method to represent the nonlinear and also the forward speed effects on a displacement high speed vessel advancing in large amplitude waves. With this objective in mind the experimental program included forced oscillation tests in heaving and pitching, for a range of periods, three different amplitudes and several speeds of advance. In head regular waves comprehensive ranges of wave periods, wave steepness and speeds, were tested in order to measure heave, pitch and loads in three cross sections.

The numerical method assumes that the radiation and diffraction hydrodynamic forces are linear and the nonlinear contributions arise from the hydrostatics and Froude–Krilov forces and the effects of green water on deck. The assumption of linearity of the radiation forces is validated by comparing calculated hydrodynamic coefficients with experimental data for three different amplitudes of the forced oscillations. Both global coefficients and sectional coefficients are compared. The motions and loads in waves are compared in terms of first and higher harmonic amplitudes and also in terms of sagging and hogging peaks.     

An Irregular Wave Generating Approach Based on naoe-FOAM-SJTU Solver

In this paper, a wave generating approach for long-crest irregular waves in a numerical tank by our in-house solver naoe-FOAM-SJTU is presented. The naoe-FOAM-SJTU solver is developed using an open source tool kit, OpenFOAM. Reynolds-averaged Navier?Stokes (RANS) equations are chosen as governing equations and the volume of fluid (VOF) is employed to capture the two phases interface. Incoming wave group is generated by imposing the boundary conditions of the tank inlet. A spectrum based correction procedure is developed to make the measured spectrum approaching to the target spectrum. This procedure can automatically adjust the wave generation signal based on the measured wave elevation by wave height probe in numerical wave tank. After 3 to 4 iterations, the measured spectrum agrees well with the target one. In order to validate this method, several wave spectra are chosen and validated in the numerical wave tank, with comparison between the final measured and target spectra. In order to investigate a practical situation, a modified Wigley hull is placed in the wave tank with incoming irregular waves. The wave-induced heave and pitch motions are treated by Fourier analysis to obtain motion responses, showing good agreements with the measurements.     

An irregular wave generating approach based on naoe-FOAM-SJTU solver

In this paper, a wave generating approach for long-crest irregular waves in a numerical tank by our in-house solver naoe-FOAM-SJTU is presented. The naoe-FOAM-SJTU solver is developed using an open source tool kit, OpenFOAM. Reynolds-averaged Navier-Stokes (RANS) equations are chosen as governing equations and the volume of fluid (VOF) is employed to capture the two phases interface. Incoming wave group is generated by imposing the boundary conditions of the tank inlet. A spectrum based correction procedure is developed to make the measured spectrum approaching to the target spectrum. This procedure can automatically adjust the wave generation signal based on the measured wave elevation by wave height probe in numerical wave tank. After 3 to 4 iterations, the measured spectrum agrees well with the target one. In order to validate this method, several wave spectra are chosen and validated in the numerical wave tank, with comparison between the final measured and target spectra. In order to investigate a practical situation, a modified Wigley hull is placed in the wave tank with incoming irregular waves. The wave-induced heave and pitch motions are treated by Fourier analysis to obtain motion responses, showing good agreements with the measurements.     

Experimental measurements of the complex motion of a suspended axisymmetric floating body in regular and near-focused waves

Numerical models which account for the multiple response modes of floating wave energy converters (WECs) in operating conditions require experimental data for validation. Measurement and observation of complex hydrodynamic mechanisms are also required to inform the development of modelling tools suitable for the simulation of response to extreme waves. Experimental measurements are reported of the motion of an axisymmetric float to regular and near-focused waves. The mechanical system, incident wave conditions and response in a 2D vertical plane are detailed to facilitate comparison to numerical simulations. The system comprises a heaving float connected to a counterweight by an inextensible cable over two pulleys to provide a simplified representation of the slowly varying surge constraint of a mooring system. Translation of the float is measured using an optical encoder. Motion in heave, surge and pitch are also determined by a position identification method based on analysis of video footage. For low frequency regular waves, the float prescribes an elliptical trajectory and the variation of response amplitude with wave amplitude is linear. At higher frequencies, drift of up to one-third of the float radius is observed and the float oscillates along an arc. More complex motions are observed due to the three large amplitude waves of a near-focused wave group. During these waves the upper surfaces of the float are partly immersed and motion occurs in heave, surge and pitch.     

Coupled analysis of nonlinear sloshing and ship motions

A coupled numerical model considering nonlinear sloshing flows and the linear ship motions has been developed based on a boundary element method. Hydrodynamic performances of a tank containing internal fluid under regular wave excitations in sway are investigated by the present time-domain simulation model and comparative model tests. The numerical model features well the hydrodynamic performance of a tank and its internal sloshing flows obtained from the experiments. In particular, the numerical simulations of the strong nonlinear sloshing flows at the natural frequency have been validated. The influence of the excitation wave height and wave frequency on ship motions and internal sloshing has been investigated. The magnitude of the internal sloshing increases nonlinearly as the wave excitation increases. It is observed that the asymmetry of the internal sloshing relative to still water surface becomes more pronounced at higher wave excitation. The internal sloshing-induced wave elevation is found to be amplitude-modulated. The frequency of the amplitude modulation envelope is determined by the difference between the incident wave frequency and the natural frequency of the internal sloshing. Furthermore, the coupling mechanism between ship motions and internal sloshing is discussed.     

Investigation of motions of catamarans in regular waves—II

By extending the linear frequency domain theory, a quasi-non-linear time-domain technique has been developed to investigate the large amplitude motions of catamarans in regular waves. The non-linearity of hydrodynamic forces included in this practical method comes from variations of a ship's submerged portion. These forces are obtained from a database generated by the linear frequency domain method at each time step. The coupled equations, heave and pitch, are solved in the time domain by using the Runge-Kutta method with proper initial values. In order to investigate the non-linear effects of large amplitude motions of the V-1 catamaran in the head-sea condition, numerical results obtained from the linear and non-linear strip methods have been compared with those obtained from a series of experiments carried out in the towing tank of the Hydrodynamics Laboratory at the University of Glasgow. Based on the comparative studies, the numerical results obtained from the time-domain program can provide better predictions for the large amplitude motions of catamarans than the linear frequency domain method. It is concluded that the non-linear effects are significant when the model speeds and wave amplitudes increase. The peak values of large amplitude motions around the resonance frequencies, as obtained from the non-linear time-domain predictions as well as from measurements, are smaller than those obtained from the linear theory.     

Experimental investigation of a fast monohull in forced harmonic motions

The paper presents the results of an experimental investigation of added masses and damping coefficients of a model of a fast monohull. A model of 4.5 m length between perpendiculars was constructed of fiber glass reinforced plastic (FRP) with four segments connected by a backbone. The backbone was instrumented with load cells at the positions of the cuts. This configuration, combined with load cells measuring the force exerted by the forced motion actuators, made it possible to obtain the hydrodynamic coefficients for each of the four hull segments.

The investigation focused on the vertical motions. Thus, the experimental program included forced harmonic heave and pitch motions in calm water (no incident waves). Subtracting inertial and restoring forces from total measured forces, one obtained the hydrodynamic component, which then resulted in the hydrodynamic coefficients. The effects of steady forward speed on the radiation forces were investigated by conducting model tests at four forward speeds. Finally, nonlinear effects were assessed by conducting model tests for three amplitudes of forced heave and forced pitch motions.     

港口内靠码头系泊船运动的计算

本文以英国南海岸处Ｓｈｏｒｅｈａｍ港内系泊船为例，研究了港口内靠码头系泊船运动的数值计算问题。该船的实船实验和模型实验已经完成，为验证本文计算结果提供了依据。运动方程在时域内求解。在计算船体脉冲响应函数时，引入了船体阻尼系数在频率很低时的渐近表达式和一迭代算法。分析了港内共振波浪和其对船体的作用力。讨论了靠码头系泊船运动的非线性特征。计算结果与实验结果符合很好。     

Investigation of motions of catamarans in regular waves—I

In this paper an analytical technique based on the two-dimensional Green function method associated with a cross-flow approach for taking viscous effects into account to estimate the motion response of catamarans in the frequency domain is presented. In order to validate this method, the numerical results are compared with experimental values obtained for two different catamarans (ASR5061 [Wahab, R., Pritchett, C. and Ruth, L.C. 1971. On the behaviour of the ASR catamaran in waves. Marine Technology, 8, 334–360] and Marintek [Faltinsen, O., Hoff, J.R., Kvalsvold, J. and Zhao, R. 1992. Global loads on high speed catamarans. 5th Int. Symp. on Practical Design of Ships and Mobile Units, University of Newcastle-upon-Tyne, 1.360–1.373]).In the second part of the paper the tests carried out with a third catamaran configuration at the Hydrodynamics Laboratory of the University of Glasgow are presented to evaluate the non-linear effects. These test results cover different speeds and wave heights at a wide range of wave frequencies. The paper concludes that the two-dimensional method correlates very well with measurements of small amplitude motions. For large amplitude motion tests, the non-linear effects become significant when the model speed and wave amplitudes increase. The peak values of heave and pitch motions measured around the resonance frequency are smaller than those obtained from the linear theory.     

Investigation of ship hydroelasticity

The importance of hydroelastic analysis of large and flexible container ships of today is pointed out. A methodology for investigation of this challenging phenomenon is drawn up and a mathematical model is worked out. It includes definition of ship geometry, mass distribution, structure stiffness, and combines ship hydrostatics, hydrodynamics, wave load, ship motion and vibrations. Based on the presented theory, a computer program is developed and applied for hydroelastic analysis of a flexible segmented barge for which model test results of motion and distortion in waves have been available. A correlation analysis of numerical simulation and measured response shows quite good agreement of the transfer functions for heave, pitch, roll, vertical and horizontal bending and torsion. The tool checked in such a way can be further used for reliable hydroelastic analysis of ship-like structures.     

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.     

The dynamics of wave induced ship motions in shallow water

The analytical method developed by Svendsen (1968) for a forced heave motion is extended to the general problem of wave induced heave, roll and sway motions of a long ship at a depth of water which is only slightly larger than the draught of the ship. This corresponds, for example, to the situation of a fully loaded ship in a harbour area.After linearization of the problem, the water motion is considered for each of the three individual motions and for the wave reflection-transmission problem for a fixed ship. The ensuing results for the forces on the ship are then synthesized to form the equations of motion, which are presented with all coefficients given, including mooring forces.Analytical and numerical results are given for the three components of motion, for the associated resonance frequencies, and for the hydrodynamic masses and moments of inertia. Finally, the assumptions used are analyzed and evaluated by comparison with measurements and with other results for a special case.     

规则波中船舶复原力和参数横摇研究

为研究规则波中船舶复原力变化规律及其对参数横摇的影响,首先,基于切片理论求解出船舶无横倾时在波浪中时间序列垂荡和纵摇运动,确定出波面与船体的相对位置;其次,利用三个坐标系之间的转换关系进而确定规则波中船体各横剖面左右舷与波面瞬时交点,求得各浸水剖面面积;然后对波浪压力沿船长湿表面积分,得出规则波中船舶复原力的Froude-Krylov部分。同时,利用作用在横倾船舶上的辐射力和绕射力,求出规则波中船舶复原力辐射力和绕射力部分。在复原力计算的基础上,确定一个参数横摇模型,实现波浪中参数横摇计算。以一艘集装箱船为例,研究了规则波中复原力变化以及参数横摇规律,复原力变化幅度是影响参数横摇的一个重要因素。     

Analysis of hydrodynamic behaviors of gravity net cage in irregular waves

Based on the lumped-mass method and rigid-body kinematics theory, a mathematical model of a gravity cage system attacked by irregular waves is developed to simulate the hydrodynamic response of cage system, including the maximum tension of mooring lines and the motion of float collar. The normalized response amplitudes (response amplitude operators) are calculated for the cage motion response in heave and surge, and the mooring line tension response, in regular waves. In addition, a statistical approach is taken to determine the motion and tension transfer functions in irregular waves. In order to validate the numerical model of a gravity cage attacked by irregular waves, numerical predictions have been compared with the experimental observations in the time and frequency domain. The effect of wave incident angle on the float collar motion, mooring line tension and net volume reduction of the gravity cage system in irregular waves is also investigated. The results show that at high frequencies, the cage system has no significant heave motion. It tends to contour itself to longer waves. The variation amplitude of mooring line forces decreases as the wave frequency increases. With the increasing of wave incident angle, the horizontal displacement of the float collar increases.     

Parametric model identification of high-speed craft dynamics

Parametric models of heave, pitch and roll dynamics of a high-speed craft have been estimated for different wave incidence angles in the frequency domain. Several issues that make the identification problem interesting are the following: type of parameterization, starting values, non-quadratic functions, excitation signals and short data record. The method employed guarantees a fine linear approximation of the nonlinear dynamics of a fast ship for the ultimate goal of stabilization control to reduce motion sickness associated with heave, pitch and roll accelerations. In addition, the approach achieves high-quality starting values and avoids non-quadratic terms in the cost function, which results in less computational load and significantly more accurate models when compared with a previous method employed for the same problem.     

Computations of self-propulsion free to sink and trim and of motions in head waves of the KRISO Container Ship (KCS) model

Two computations of the KCS model with motions are presented. Self-propulsion in model scale free to sink and trim are studied with the rotating discretized propeller from the Hamburg Model Basin (HSVA) at Fr = 0.26. This case is particularly complex to simulate due to the close proximity of the propeller to the rudder. The second case involves pitch and heave in regular head waves. Computations were performed with CFDShip-Iowa version 4.5, a RANS/DES CFD code designed for ship hydrodynamics. The self-propulsion computations were carried out following the procedure described in Carrica et al. [1], in which a speed controller is used to find the propeller rotational speed that results in the specified ship velocity. The rate of revolutions n, sinkage, trim, thrust and torque coefficients K_T, K_Q and resistance coefficient CT_(SP) are thus obtained. Comparisons between CFD and EFD show that the rate of revolutions n, thrust and torque coefficients K_T and K_Q have higher prediction accuracies than sinkage and trim. For the simulation of pitch and heave in head waves, the geometry includes KCS hull and rudder under three conditions with two Froude numbers and three wave length and amplitude combinations. 0th and 1st harmonic amplitudes and 1st harmonic phase are computed for total resistance coefficient C_T, heave motion z and pitch angle θ. Comparisons between CFD and EFD show that pitch and heave are much better predicted than the resistance. In both cases comparisons with simulations by other authors presented at the G2010 CFD Workshop [2] using different CFD methodologies are included.