Nonlinear heave radiation forces on two-dimensional single and twin hulls

A study of nonlinear heave radiation of two-dimensional single and double hulls has been carried out in the time domain. The problem is analyzed by means of a fully nonlinear mathematical model, referred to as the mixed Eulerian–Lagrangian (MEL) model, which is based on an integral relation formulation coupled with time-integration of the nonlinear free-surface boundary conditions. The integral equation solver is based on a cubic-spline boundary-element scheme in which both potential and velocity continuity conditions can be enforced through the intersection points. The body undergoes periodic forced heave oscillation. By implementing effective wave-absorbing beaches at the two ends of the rectangular numerical tank, long-term steady-state force-histories could be achieved consistently in all computations.Results in terms of radiation forces for rectangular and triangular single- and twin-hull geometries are presented and discussed. Linear hydrodynamic forces in terms of added-mass and damping are validated for the rectangular hull. The Fourier-analyzed results reveal the extent of nonlinear (higher-order) components in the force-signals over different parameters which include the amplitudes of oscillation, hull-spacing for the twin-hulls and water depth.     

Time-domain computation of large amplitude 3D ship motions with forward speed

Time-domain computations of 3D ship motions with forward speed are presented in this paper. The method of computation is based upon transient Green function. Both linear and nonlinear (large-amplitude) computations are performed where the included nonlinearities are those arising from the incident wave, but the diffraction and radiation forces are otherwise retained as linear. The incident wave can be described by any explicit nonlinear model. Computations over a variety of wave and speed parameters establish the robustness of the algorithm, which include high speed and following waves. Comparison of linear and nonlinear computations show that nonlinearities have a considerable influence on the results, particularly in predicting the instantaneous location of the hull in relation to the wave, which is crucial in determining forefoot emergence and deck wetness.     

Time-domain wave diffraction of two-dimensional single and twin hulls

The nonlinear diffraction of 2D single and twin hulls are studied by employing a mixed Eulerian–Lagrangian model based on a higher-order cubic-spline boundary element solver. Two types of simulations are considered. In the first, waves are generated by a piston-type wave-maker in a rectangular tank and in the second case a nonlinear incident wave is assumed to exist in the tank in which the body is introduced. For the application of this model, the full nonlinear diffraction problem is recast in terms of a perturbation wave-field. Computations are performed for rectangular and triangular hull geometries. Computed results show significant nonlinearities, particularly in the heave force. The twin hull results show the influence of wave interference on the diffraction forces. This interference influences the surge force considerably, but heave force is less affected.     

Time domain three-dimensional fully nonlinear computations of steady body–wave interaction problem

Three-dimensional fully nonlinear waves generated by moving disturbances with steady forward speed without motions are solved using a mixed Eulerian–Lagrangian method in terms of an indirect boundary integral method and a Runge–Kutta time marching approach which integrates the fully nonlinear free surface boundary conditions with respect to time.A moving computational window is used in the computations by truncating the fluid domain (the free surface) into a computational domain. The computational window maintains the computational domain and tracks the free surface profile by a node-shifting scheme applied within it. An implicit implement of far field condition is enforced automatically at the truncation boundary of the computational window.Numerical computations are applied to free surface waves generated by Wigley and Series 60 hulls for the steady problem. The present numerical results are presented and compared with existing linear theory, experimental measurements, and other numerical nonlinear computations. The comparisons show satisfactory agreements for these hydrodynamic problems.     

An efficient domain decomposition strategy for wave loads on surface piercing circular cylinders

A fully nonlinear domain decomposed solver is proposed for efficient computations of wave loads on surface piercing structures in the time domain. A fully nonlinear potential flow solver was combined with a fully nonlinear Navier–Stokes/VOF solver via generalized coupling zones of arbitrary shape. Sensitivity tests of the extent of the inner Navier–Stokes/VOF domain were carried out. Numerical computations of wave loads on surface piercing circular cylinders at intermediate water depths are presented. Four different test cases of increasing complexity were considered; 1) weakly nonlinear regular waves on a sloping bed, 2) phase-focused irregular waves on a flat bed, 3) irregular waves on a sloping bed and 4) multidirectional irregular waves on a sloping bed. For all cases, the free surface elevation and the inline force were successfully compared against experimental measurements.     

系泊船非线性波浪力时域计算:二维模型

为找到具有工程实用价值的港口系泊船波浪力的时域计算方法,建立了在港口中存在系泊船时非线性波浪力时域计算的垂直二维耦合模型:用Boussinesq方程计算船的两侧的外域,用欧拉方程计算船底面下的内域,两域在交界面处的连接条件是流量连续和压力相等.将复平面内的边界元方法应用于所研究问题,对耦合模型进行了验证.进行了相关模型实验,实验结果与数值计算结果比较表明这两种数值计算模型都具有满意的精度,但耦合模型的计算效率要远远高于边界元方法的计算效率.本耦合模型的数学处理简单,可适用于工程计算.     

Practical evaluation of sinkage and trim effects on the drag of a common generic freely floating monohull ship

A practical method to account for the influence of sinkage and trim on the drag of a freely floating (free to sink and trim) common monohull ship at a Froude number F ≤ 0.45 is considered. The sinkage and the trim are estimated via two alternative simple methods, considered previously. The drag is also estimated in a simple way, based on the classical Froude decomposition into viscous and wave components. Specifically, well-known semiempirical expressions for the friction drag, the viscous pressure drag and the drag due to hull roughness are used, and the wave drag is evaluated via a practical linear potential flow method. This simple approach can be used for ship models as well as full-scale ships with smooth or rough hull surfaces, and is well suited for early ship design and optimization. The method considered here to determine the sinkage and the trim, and their influence on the drag, yields theoretical predictions of the drag of the Wigley, S60 and DTMB5415 hulls that are much closer to experimental measurements than the corresponding predictions for the hull surfaces of the ships in equilibrium position at rest. These numerical results suggest that sinkage and trim effects, significant at Froude numbers 0.25 < F, on the drag of a typical freely floating monohull ship can be realistically accounted for in a practical manner that only requires simple potential flow computations without iterative computations for a sequence of hull positions.     

压力公式的修正及其在波浪荷载计算中的应用

计入前进速度流动与非定常流动之间的干扰,修正了切片法中计算压力的公式,从而建立了计算波浪荷载的方法,开发了相应的计算机程序。计算实例同国内外有关资料的比较表明结果吻合良好,为船舶与海洋结构物结构有限元分析程序提供了有效的前处理模块     

Efficient primary energy conversion in irregular waves

This paper studies reactive control of wave-energy devices in irregular waves. Estimates of future velocity of the primary energy converter are based on a time history of past velocity measurements, and prediction of incident waves is not attempted. Time–domain computations for a single-degree-of-freedom oscillating water column in irregular waves reveal that control is not optimal; but would enhance efficiency at low peak frequencies where energy-rich swells predominate, and allow appreciably smaller devices to fulfill a given power requirement.     

An inverse hull design approach in minimizing the ship wave

The Levenberg–Marquardt Method (LMM) and a panel code for solving the wave-making problem are utilized in an inverse hull design problem for minimizing the wave of ships. A typical catamaran is selected as the example ship for the present study. The hull form of the catamaran is described by the B-spline surface method so that the shape of the hull can be completely specified using only a small number of parameters (i.e. control points). The technique of parameter estimation for the inverse design problem is thus chosen. The LMM of parameter estimation, which is the combination of steepest descent and Newton’s methods, has been proven to be a powerful tool for the inverse shape design problem. For this reason it is adopted in the present study.In the present studies, the inverse hull design method can not only be applied to estimate the hull form based on the known wave data of the target ship but can also be applied to estimate the unknown hull form based on the reduced wave height. The optimal hull forms of minimizing wave for a typical catamaran in deep water at service speed and at the critical speed of shallow water are estimated, respectively. Moreover, a new hull form with the combining feature of the optimal hull forms for deep water and shallow water is performing well under both conditions. The numerical simulation indicates that the hull form designed by inverse hull design method can reduce the ship wave significantly in comparison with the original hull form.     

Numerical wave tank computations of nonlinear motions of two-dimensional arbitrarily shaped free floating bodies

Time domain simulations of nonlinear motions of two-dimensional floating bodies in waves are presented. The so called `body exact' approach is adopted in a numerical wave tank. A new scheme for pressure evaluation on the wetted hull is developed and systematically used with good results in terms of accuracy and stability. Strongly flared geometries are successfully handled even at very large amplitude motions. The validation of the code is carried out according to the 20th ITTC Seakeeping Committee recommendations through internal checks of consistency and through comparisons with available experimental data.     

A proposed adiabatic formulation of 3‐dimensional global atmospheric models based on potential vorticity

A 2‐time‐level finite difference atmospheric general circulation model based on the semi‐Lagrangian advection of pseudo potential vorticity (which becomes potential vorticity in that part of the domain where the hybrid vertical coordinate becomes isentropic) has been formulated. At low levels, the hybrid vertical coordinate is terrain following. The problem of isentropic potential vorticity possibly becoming ill‐defined in the regions of planetary boundary layer is thus circumvented. The divergence equation is a companion to the (pseudo) potential vorticity equation and the model is thus called a PV‐D model. Many features of a previously developed shallow water PV‐D model are carried over: a modification of the PV equation needed to give computational stability of long Rossby waves; a semi‐Lagrangian semi‐implicit treatment of both the linear and the nonlinear terms; the use of an unstaggered grid in the horizontal; the use of a nonlinear multigrid technique to solve the nonlinear implicit equations. A linear numerical stability analysis of the model's gravity–inertia waves indicates that the potential temperature needs to be separated into horizontal mean and perturbation parts. This allows an implicit treatment of the vertical advection associated with the mean in the thermodynamic equation. Numerical experiments with developing baroclinic waves have been carried out and give realistic results.     

Direct time domain analysis of floating structures with linear and nonlinear mooring stiffness in a 3D numerical wave tank

In this paper, motion response of a moored floating structure interacting with a large amplitude and steep incident wave field is studied using a coupled time domain solution scheme. Solution of the hydrodynamic boundary value problem is achieved using a three-dimensional numerical wave tank (3D NWT) approach based upon a form of Mixed-Eulerian–Lagrangian (MEL) scheme. In the developed method, nonlinearity arising due to incident wave as well as nonlinear hydrostatics is completely captured while the hydrodynamic interactions of radiation and diffraction are determined at every time step based on certain simplifying approximations. Mooring lines are modelled as linear as well as nonlinear springs. The horizontal tension for each individual mooring line is obtained from the nonlinear load-excursion plot of the lines computed using catenary theory, from which the linear and nonlinear line stiffness are determined. Motions of three realistic floating structures with different mooring systems are analyzed considering various combinations of linear and approximate nonlinear hydrodynamic load computations and linear/nonlinear mooring line stiffness. Results are discussed to bring out the influence and need for consideration of nonlinearities in the hydrodynamics and hydrostatics as well as the nonlinear modelling of the line stiffness.     

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.     

Methods for ultimate limit state assessment of ships and ship-shaped offshore structures: Part I—Unstiffened plates

The present paper is Part I of a series of three papers prepared by the authors on the methods useful for ultimate limit state assessment of marine structures, that have been developed in the literature during the last few decades. It is considered that such methods are now mature enough to enter day-by-day design and strength assessment practice. The aims of the three papers are to conduct some benchmark studies of such methods on ultimate limit state assessment of (unstiffened) plates, stiffened panels, and hull girders of ships and ship-shaped offshore structures, using some candidate methods such as ANSYS nonlinear finite element analysis (FEA), DNV PULS, ALPS/ULSAP, ALPS/HULL, and IACS common structural rules (CSR) methods. As an illustrative example, an AFRAMAX-class hypothetical double hull oil tanker structure designed by CSR method is studied. In the present paper (Part I), the ultimate limit state assessment of unstiffened plates under combined biaxial compression and lateral pressure loads is emphasized using ANSYS, DNV PULS, and ALPS/ULSAP methods, and their resulting computations are compared. Part II will deal with methods for the ultimate limit state assessment of stiffened panels under combined biaxial compression and lateral pressure using ANSYS, DNV PULS, and ALPS/ULSAP methods, and Part III will treat methods for the progressive collapse analysis of the hull structure using ANSYS, ALPS/HULL, and IACS CSR methods.     

Computations of forces on, and particle orbits around, horizontal cylinders under steep waves

A numerical method, based on a boundary integral equation combined with a non-linear time stepping procedure for the free water surface, is developed for simulations of the interaction between highly non-linear water waves and submerged horizontal cylinders. The method is based on potential theory, and the omission of viscous effects restricts the wave-structure interaction computations to low Keulegan-Carpenter numbers where inertia forces are dominant. The numerical scheme is verified by computations with a steep wave of exact form during several wave periods, and by computations of a breaking wave. A new method for tracing the orbits of water particles in the fluid domain is developed, and the influence from submerged structures on the orbits is visualized through several computational examples. The wave forces on submerged structures are computed and are found to correspond well with other computed results for low Keulegan-Carpenter numbers.     

Time stepping solutions of the two-dimensional nonlinear wave radiation problem

The two-dimensional nonlinear time domain free surface flow problem is analysed using potential flow theory. The problem is solved by a time marching method. At each time step two numerical approaches are used. One is based on the boundary element method in the complex plane. The complex potential is assumed to vary linearly within each element and the solution is obtained by imposing the boundary conditions at the nodes of the elements. The other approach is based on the finite element formulation. Triangular elements and linear shape functions are used. The solution is obtained by the Galerkin method. Numerical results are obtained for the wave elevation generated by a vertical wave maker. Results are also provided for a circular cylinder oscillating below the free surface. For these cases the finite element method is found to provide substantially more efficient computations than the boundary element method using equivalent discretizations.     

Simulation of free surface flow around a VLCC hull using viscous and potential flow methods

Numerical simulations have been carried out to determine the incompressible free surface flow around a VLCC hull form for which experimental results are available. A commercial viscous flow finite volume code using the two-phase Eulerian–Eulerian fluid approach and a potential flow code based on the Rankine source method have been used in this study. The simulation conditions are the ones for which experimental results exist. The shear stress transport (SST) turbulence model has been used in the viscous flow code. A tetrahedral unstructured grid was used with the viscous flow code for meshing the computational domain, while quadrilateral structural patches were used with the potential flow code for meshing the VLCC hull surface and the water surface around it. The results compare well with the available experimental data and they allow an understanding of the differences that can be expected from viscous and potential flow methods as a result of their different mathematical formulations, which make their complementary application useful for determining the total ship resistance.     

Coupled dynamic analysis for wave interaction with a truss spar and its mooring line/riser system in time domain

A full time-domain analysis program is developed for the coupled dynamic analysis of offshore structures. For the hydrodynamic loads, a time domain second order method is developed. In this approach, Taylor series expansions are applied to the body surface and free-surface boundary conditions, and the Stokes perturbation procedure is then used to establish the corresponding boundary value problems with time-independent boundaries. A higher-order boundary element method (HOBEM) is developed to calculate the velocity potential of the resulting flow field at each time step. The free-surface boundary condition is satisfied to the second order by fourth order Adams–Bashforth–Moultn method. An artificial damping layer is adopted on the free surface to avoid the wave reflection. The mooring-line/tendon/riser dynamics are based on the rod theory and the finite element method (FEM), with the governing equations described in a global coordinate system. In the coupled dynamic analysis, the motion equation for the hull and dynamic equations for mooring-lines/tendons/risers are solved simultaneously using the Newmark method. The coupled analysis program is applied for a truss Spar motion response simulation. Numerical results including motions and tensions at the top of mooring-lines/risers are presented, and some significant conclusions are derived.     

Simulating nonlinear waves on the free surface in surf zones with two-dimensional sloping beach

Unsteady nonlinear wave motions on the free surface in shallow water and over slopes of various geometries are numerically simulated using a finite difference method in rectangular grid system. Two-dimensional Navier–Stokes equations and the continuity equation are used for the computations. Irregular leg lengths and stars are employed near the boundaries of body and free surface to satisfy the boundary conditions. Also, the free surface which consists of markers or segments is determined every time step with the satisfaction of kinematic and dynamic free surface conditions. Moreover, marker-density method is also adopted to allow plunging jets impinging on the free surface. Either linear or Stokes wave theory is employed for the generation of waves on the inflow boundary. For the simulation of wave breaking phenomena, the computations are carried out with various wave periods and sea bottom slopes in surf zone. The results are compared with other existing computational and experimental results. Agreement between the experimental data and the computation results is good.