基于三维波浪水池方法的浮体气隙估计方法研究（英文）

Computations for air gap response of a semisubmersible platform based on a 3D numerical wave tank approach are presented.The developed method is in time domain and can consider nonlinearities associated with incident wave and hydrostatic forces exactly in determining the body response, but the interaction hydrodynamics of radiation and diffraction are based on simplified linearization assumptions. The incident wave can be defined by any suitable wave theory and here defined by a fully nonlinear numerical wave model. After verifying the present computations results in its degenerated linearized version against the usual linear 3D Green function–based frequency-domain results for air gap predictions, systematic comparative studies are undertaken between linear and the approximate nonlinear solutions. It is found that nonlinear computations can yield considerably conservative predictions as compared to fully linear calculations, amounting to a difference of up to 30%–40% in the minimum air gap in steep ambient incident waves at high and moderate frequencies.     

3-D computational method of wave loads on turret moored FPSO tankers

A three-dimensional method of calculating wave loads of turret moored FPSO (Floating Production Storage and Offloading) tankers is presented. The linearized restoring forces acting on the ship hull by the mooring system are calculated according to the catenary theory, which are expressed as the function of linear stiffness coefficients and the displacements of the upper ends of mooring chains. The hydrodynamic coefficients of the ship are calculated by the three-dimensional potential flow theory of the linear hydrodynamic problem for ships with a low forward speed. The equations of ship motions are established with the effect of the restoring forces from the mooring system included as linear stiffness coefficients. The equations of motions are solved in frequency domain, and the responses of wave-induced motions and loads on the ship can be obtained. A computer program based on this method has been developed, and some calculation examples are illustrated. Analysis results show that the method can give satisfying prediction of wave loads.     

船体周围Stokes波浪流耦合VOF模型的数值模拟（英文）

The surface wave generated by flow around a ship hull moving near free surface of water is simulated numerically in this study. The three-dimensional implicit finite volume method(FVM) is applied to solve Reynolds averaged Navier-Stokes(RANS) equation. The realizable k-ε turbulence model has been implemented to capture turbulent flow around the ship hull in the free surface zone. The volume of fluid(VOF) method coupled with the Stokes wave theory has been used to determine the free surface effect of water. By using is a six degrees of freedom model, the ship hull's movement is numerically solved with the Stokes wave together. Under the action of Stokes waves on the sea, the interface between the air and water waves at the same regular pattern and so does the pressure and the vertical velocity. The ship hull moves in the same way as the wave. The amplitude of the ship hull's heave is less than the wave height because of the viscosity damping. This method could provide an important reference for the study of ships' movement, wave and hydrodynamics.     

高速滑行艇在规则波中运动的RANSE数值模拟(英文)

This paper presents a study on the numerical simulation of planing crafts sailing in regular waves. This allows an accurate estimate of the seas keeping performance of the high speed craft. The simulation set in six-degree of freedom motions is based on the Reynolds averaged Navier Stokes equations volume of fluid (RANSE VOF) solver. The trimming mesh technique and integral dynamic mesh method are used to guarantee the good accuracy of the hydrodynamic force and high efficiency of the numerical simulation. Incident head waves, oblique waves and beam waves are generated in the simulation with three different velocities (Fn =1.0, 1.5, 2.0). The motions and sea keeping performance of the planing craft with waves coming from different directions are indicated in the flow solver. The ship designer placed an emphasis on the effects of waves on sailing amplitude and pressure distribution of planing craft in the configuration of building high speed crafts.     

采用解析和数值方法分析浮动式矩形浮箱的波浪诱导载荷研究（英文）

A three-dimensional mathematical hydrodynamic model associated with surface wave radiation by a floating rectangular box-type structure due to heave, sway, and roll motions in finite water depth is investigated based on small amplitude water wave theory and linear structural response. The analytical expressions for the radiation potentials, wave forces, and hydrodynamic coefficients are presented based on matched eigenfunction expansion method(MEFEM). The correctness of the analytical results of...     

Numerical Simulation of the Stokes Wave for the Flow around a Ship Hull Coupled with the VOF Model

The surface wave generated by flow around a ship hull moving near free surface of water is simulated numerically in this study. The three-dimensional implicit finite volume method(FVM) is applied to solve Reynolds averaged Navier-Stokes(RANS) equation. The realizable k-ε turbulence model has been implemented to capture turbulent flow around the ship hull in the free surface zone. The volume of fluid(VOF) method coupled with the Stokes wave theory has been used to determine the free surface effect of water. By using is a six degrees of freedom model, the ship hull’s movement is numerically solved with the Stokes wave together. Under the action of Stokes waves on the sea, the interface between the air and water waves at the same regular pattern and so does the pressure and the vertical velocity. The ship hull moves in the same way as the wave. The amplitude of the ship hull’s heave is less than the wave height because of the viscosity damping. This method could provide an important reference for the study of ships’ movement, wave and hydrodynamics.     

冰盖下两层流体水波与可渗透起伏海床的作用(英文)

The scattering problem involving water waves by small undulation on the porous ocean-bed in a two-layer fluid,is investigated within the framework of the two-dimensional linear water wave theory where the upper layer is covered by a thin uniform sheet of ice modeled as a thin elastic plate.In such a two-layer fluid there exist waves with two different modes,one with a lower wave number propagate along the ice-cover whilst those with a higher wave number propagate along the interface.An incident wave of a particular wave number gets reflected and transmitted over the bottom undulation into waves of both modes.Perturbation analysis in conjunction with the Fourier transform technique is used to derive the first-order corrections of reflection and transmission coefficients for both the modes due to incident waves of two different modes.One special type of bottom topography is considered as an example to evaluate the related coefficients in detail.These coefficients are depicted in graphical forms to demonstrate the transformation of wave energy between the two modes and also to illustrate the effects of the ice sheet and the porosity of the undulating bed.     

The choice of speed and clearance for RAS on 3D method

In this paper, a 3D source distribution technique is used to calculate the coupled motions between two ships which advance in the wave with the same speed. The numerical results of coupled motions for a frigate and a supply ship have a good agreement with the experimental results. Based on the 3D coupled motions of two ships, a spectral analysis is employed to clearly observe the effect of speed, clearance and wave heading on the significant relative motion amplitude (SRMA) of two ships.The method presented in this paper will be helpful to select suitable clearance, speed and wave heading for underway replenishment at sea(RAS).     

Coupled motions of two ships in irregular waves in time domain

A three-dimensional time domain approach is used to study the coupled motions of two ships with forward speed in waves. In this approach, the boundary condition is satisfied on the mean wetted hull surface of the moving bodies and the free surface condition is linearized. The problem is solved by using a transient free-surface Green function source distribution on the submerged hulls. After solving the response amplitude operator, the method of spectral analysis is employed to clearly express the motion energy spectrum and significant amplitude of two ships. For verifying the code, two same circular cylinders at beam wave are selected to calculate coupled motions by comparison with the results obtained by 3D frequency method which has been proved to be efficient for solving such problems.Two Wigley ships of different sizes with the same forward speed are chosen for numerical calculation of the interaction effect, and some useful suggestions are obtained for underway replenishment at sea.     

可渗透海床上具有小底部变形的斜弯曲重力波的相互作用:格林函数法（英文）

The interaction of oblique incident water waves with a small bottom deformation on a porous ocean-bed is examined analytically here within the framework of linear water wave theory. The upper surface of the ocean is assumed to be covered by an infinitely extended thin uniform elastic plate, while the lower surface is bounded by a porous bottom surface having a small deformation. By employing a simplified perturbation analysis, involving a small parameter δ(1), which measures the smallness of the deformation, the governing Boundary Value Problem(BVP) is reduced to a simpler BVP for the first-order correction of the potential function. This BVP is solved using a method based on Green's integral theorem with the introduction of suitable Green's function to obtain the first-order potential, and this potential function is then utilized to calculate the first-order reflection and transmission coefficients in terms of integrals involving the shape function c(x) representing the bottom deformation. Consideration of a patch of sinusoidal ripples shows that when the quotient of twice the component of the incident field wave number propagating just below the elastic plate and the ripple wave number approaches one, the theory predicts a resonant interaction between the bed and the surface below the elastic plate. Again, for small angles of incidence, the reflected wave energy is more as compared to the other angles of incidence. It is also observed that the reflected wave energy is somewhat sensitive to the changes in the flexural rigidity of the elastic plate, the porosity of the bed and the ripple wave numbers. The main advantage of the present study is that the results for the values of reflection and transmission coefficients obtained are found to satisfy the energy-balance relation almost accurately.     

CFD技术在船舶与海洋工程复杂粘性流动中的应用进展（英文）

Complex flow around floating structures is a highly nonlinear problem, and it is a typical feature in ship and ocean engineering.Traditional experimental methods and potential flow theory have limitations in predicting complex viscous flows. With the improvement of high-performance computing and the development of numerical techniques, computational fluid dynamics(CFD) has become increasingly powerful in predicting the complex viscous flow around floating structures. This paper reviews the recent progress in CFD techniques for numerical solutions of typical complex viscous flows in ship and ocean engineering.Applications to free-surface flows, breaking bow waves of high-speed ship, ship hull–propeller–rudder interaction, vortexinduced vibration of risers, vortex-induced motions of deep-draft platforms, and floating offshore wind turbines are discussed.Typical techniques, including volume of fluid for sharp interface, dynamic overset grid, detached eddy simulation, and fluid–structure coupling, are reviewed along with their applications. Some novel techniques, such as high-efficiency Cartesian grid method and GPU acceleration technique, are discussed in the last part as the future perspective for further enhancement of accuracy and efficiency for CFD simulations of complex flow in ship and ocean engineering.     

Wave load computation in direct strength analysis of semi-submersible platform structures

A wave load computation approach in direct strength analysis of semi-submersible platform structures was presented in this paper. Considering the differences in shape of pontoon, column and beam, the combination of accumulative chord length cubic parameter spline theory and analytic method was adopted for generating the wet surface mesh of platform. The hydrodynamic coefficients of platform were calculated by the three-dimensional potential flow theory of the linear hydrodynamic problem for platform with low forward speed. The equation of platform motions was established and solved in frequency domain, and the responses of wave-induced loads on the platform can be obtained. With the interpolation method being utilized, the pressure loads on shell elements for finite element analysis(FEA) were converted from those on the hydrodynamic computation mesh, which pave the basis for FEA with commercial software. A computer program based on this method has been developed ,and a calculation example of semi-submersible platform was illustrated. Analysis results show that this method is a satisfying approach of wave loads computation for this kind of platform.     

液舱晃荡与浮体非线性耦合运动分析(英文)

Nonlinear interactions among incident wave, tank-sloshing and floating body coupling motion are investigated. The fully nonlinear sloshing and body-surface nonlinear free surface hydrodynamics is simulated using a Non-Uniform Rational B-Spline (NURBS) higher-order panel method in time domain based on the potential theory. A robust and stable improved iterative procedure (Yan and Ma, 2007) for floating bodies is used for calculating the time derivative of velocity potential and floating body motion. An energy dissipation condition based on linear theory adopted by Huang (2011) is developed to consider flow viscosity effects of sloshing flow in nonlinear model. A two-dimensional tank model test was performed to identify its validity. The present nonlinear coupling sway motion results are subsequently compared with the corresponding Rognebakke and Faltinsen (2003)’s experimental results, showing fair agreement. Thus, the numerical approach presented in this paper is expected to be very efficient and realistic in evaluating the coupling effects of nonlinear sloshing and body motion.     

船首形状对自由液面水下航行体阻力影响的试验研究（英文）

In this paper, towing tank experiments are conducted to study the behavior of flow on a model of the underwater vehicle with various shapes of bows, i.e. tango and standard bows in free surface motion tests. The total resistances for different Froude numbers are considered experimentally. The towing tank is equipped with a trolley that can operate in through 0.05–6 m/s speed with ±0.02 m/s accuracy. Furthermore, the study is done on hydrodynamic coefficients i.e. total, residual and friction resistance coefficients, and the results are compared. Finally, the study on flow of wave fields around bows is done and wave filed around two bows are compared. The Froude number interval is between 0.099 and 0.349. Blockage fraction for the model is fixed to 0.005 3. The results showed that the residual resistance of the standard bow in 0.19 to 0.3 Froude number is more than the tango bow in surface motion which causes more total resistance for the submarine. Finally, details of wave generated by the bow are depicted and the effects of flow pattern on resistance drag are discussed.     

基于B样条应用的船舶波浪增阻计算（英文）

Making an exact computation of added resistance in sea waves is of high interest due to the economic effects relating to ship design and operation. In this paper, a B-spline based method is developed for computation of added resistance. Based on the potential flow assumption, the velocity potential is computed using Green's formula. The Kochin function is applied to compute added resistance using Maruo's far-field method, the body surface is described by a B-spline curve and potentials and normal derivation of potentials are also described by B-spline basis functions and B-spline derivations. A collocation approach is applied for numerical computation, and integral equations are then evaluated by applying Gauss–Legendre quadrature. Computations are performed for a spheroid and different hull forms; results are validated by a comparison with experimental results. All results obtained with the present method show good agreement with experimental results.     

完全非线性方法的两个三维自由浮体瞬态行为数值模拟研究(英文)

Two floating structures in close proximity are very commonly seen in offshore engineering. They are often subjected to steep waves and, therefore, the transient effects on their hydrodynamic features are of great concern. This paper uses the quasi arbitrary Lagrangian-Eulerian finite element method (QALE-FEM), based on the fully nonlinear potential theory (FNPT), to numerically investigate the interaction between two 3-D floating structures, which undergo motions with 6 degrees of freedom (DOFs), and are subjected to waves with different incident angles. The transient behaviours of floating structures, the effect of the accompanied structures, and the nonlinearity on the motion of and the wave loads on the structures are the main focuses of the study. The investigation reveals an important transient effects causing considerably larger structure motion than that in steady state. The results also indicate that the accompanied structure in close proximity enhances the interaction between different motion modes and results in stronger nonlinearity causing 2nd-order component to be of similar significance to the fundamental one.     

Mathematical modelling of response amplitude operator for roll motion of a floating body: Analysis in frequency domain with numerical validation

This paper investigates mathematical modelling of response amplitude operator （RAO） or transfer function using the frequency-based analysis for uncoupled roll motion of a floating body under the influence of small amplitude regular waves. The hydrodynamic coefficients are computed using strip theory formulation by integrating over the length of the floating body. Considering sinusoidal wave with frequency （ω ） varying between 0.3 rad/s and 1.2 rad/s acts on beam to the floating body for zero forward speed, analytical expressions of RAO in frequency domain is obtained. Using the normalization procedure and frequency based analysis, group based classifications are obtained and accordingly governing equations are formulated for each case. After applying the fourth order Runge-Kutta method numerical solutions are obtained and relative importance of the hydrodynamic coefficients is analyzed. To illustrate the roll amplitude effects numerical experiments have been carried out for a Panamax container ship under the action of sinusoidal wave with a fixed wave height. The effect of viscous damping on RAO is evaluated and the model is validated using convergence, consistency and stability analysis. This modelling approach could be useful to model floating body dynamics for higher degrees of freedom and to validate the result.     

斜浪下FLNG与LNGC旁靠时相对运动的时域分析（英文）

Strong hydrodynamic interactions during the side-by-side offloading operation between floating liquefied natural gas(FLNG) and liquefied natural gas carrier(LNGC) can induce high risks of collision. The weather vane effect of a single-point mooring system normally results in the satisfactory hydrodynamic performance of the side-by-side configuration in head seas. Nevertheless, the changes in wave directions in real sea conditions can significantly influence the relative motions. This article studies the relative motions of the side-by-side system by using the theoretical analysis method and the numerical calculation method. Based on the three-dimensional potential theory modified by artificial damping-lid method, the frequency-domain hydrodynamic coefficients can be improved to calculate the retardation functions for the multi-body problem. An in-house code is then developed to perform the time-domain simulation of two vessels, through which the relative motions are subsequently obtained. A range of oblique waves are chosen for the extensive calculation of relative motions between the two vessels, which are further analyzed in terms of the phase shift of motion responses induced by specific resonant wave patterns. Investigation results show that wave directions have a significant influence on the relative sway, roll, and yaw motions. Under the circumstance that the absolute phase shift between the roll motions of two vessels approaches 180°, stronger relative motions are induced when LNGC is on the weather side.Moreover, the gap water resonances at high frequencies tend to cause the dangerous opposed oscillation of two vessels in the sway and yaw modes, whereas FLNG reduces the gap water resonances and relative motions when located on the weather side.     

高速滑行船体的气水界面效应模拟（英文）

High-speed planing crafts have successfully evolved through developments in the last several decades. Classical approaches such as inviscid potential flow–based methods and the empirically based Savitsky method provide general understanding for practical design. However, sometimes such analyses suffer inaccuracies since the air–water interface effects, especially in the transition phase, are not fully accounted for. Hence, understanding the behaviour at the transition speed is of fundamental importance for the designer. The fluid forces in planing hulls are dominated by phenomena such as flow separation at various discontinuities viz.,knuckles, chines and transom, with resultant spray generation. In such cases, the application of potential theory at high speeds introduces limitations. This paper investigates the simulation of modelling of the pre-planing behaviour with a view to capturing the air–water interface effects, with validations through experiments to compare the drag, dynamic trim and wetted surface area.The paper also brings out the merits of gridding strategies to obtain reliable results especially with regard to spray generation due to the air–water interface effects. The verification and validation studies serve to authenticate the use of the multi-gridding strategies on the basis of comparisons with simulations using model tests. It emerges from the study that overset/chimera grids give better results compared with single unstructured hexahedral grids. Two overset methods are investigated to obtain reliable estimation of the dynamic trim and drag, and their ability to capture the spray resulting from the air–water interaction. The results demonstrate very close simulation of the actual flow kinematics at steady-speed conditions in terms of spray at the air–water interface, drag at the pre-planing and full planing range and dynamic trim angles.