Next Article in Journal
An Invasive Mussel (Arcuatula senhousia, Benson 1842) Interacts with Resident Biota in Controlling Benthic Ecosystem Functioning
Next Article in Special Issue
Stability and Seakeeping of Marine Vessels
Previous Article in Journal
Multidisciplinary Analysis of Cystoseira sensu lato (SE Spain) Suggest a Complex Colonization of the Mediterranean
Previous Article in Special Issue
Evaluation of the Effect of Container Ship Characteristics on Added Resistance in Waves
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Comparative Seakeeping Analysis of the Full Scale KCS by Several Hydrodynamic Approaches

Department of Naval Architecture, “Dunarea de Jos” University of Galati, 800008 Galati, Romania
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2020, 8(12), 962; https://doi.org/10.3390/jmse8120962
Submission received: 17 October 2020 / Revised: 15 November 2020 / Accepted: 19 November 2020 / Published: 25 November 2020
(This article belongs to the Special Issue Stability and Seakeeping of Marine Vessels)

Abstract

:
The main transport channel of the global economy is represented by shipping. Engineers and hull designers are more preoccupied in ensuring fleet safety, the proper operation of the ships, and, more recently, compliance with International Maritime Organization (IMO) regulatory incentives. Considerable efforts have been devoted to in-depth understanding of the hydrodynamics mechanism and prediction of ship behavior in waves. Prediction of seakeeping performances with a certain degree of accuracy is a demanding task for naval architects and researchers. In this paper, a fully numerical approach of the seakeeping performance of a KRISO (Korea Research Institute of Ships and Ocean Engineering, Daejeon, South Korea) container ship (KCS) container vessel is presented. Several hydrodynamic methods have been employed in order to obtain accurate results of ship hydrodynamic response in regular waves. First, an in-house code DYN (Dynamic Ship Analysis, “Dunarea de Jos” University of Galati, Romania), based on linear strip theory (ST) was used. Then, a 3D fully nonlinear time-domain Boundary Element Method (BEM) was implemented, using the commercial code SHIPFLOW (FLOWTECH International AB, Gothenburg, Sweden). Finally, the commercial software NUMECA (NUMECA International, Brussels, Belgium) was used in order to solve the incompressible unsteady Reynolds-averaged Navier–Stokes equation (RANSE) flow at ship motions in head waves. The results obtained using these methods are represented and discussed, in order to establish a methodology for estimating the ship response in regular waves with accurate results and the sensitivity of hydrodynamical models.

1. Introduction

Shipping, often considered as the main transport channel of the global economy, is responsible for approximately 80 percent of world trade. An increasing focus on more environmentally friendly shipping pushes hull designers to further ensure fleet safety, proper operation of the ships, and, more recently, compliance with International Maritime Organization (IMO) regulatory incentives regarding the energy efficiency operational indicator (EEOI) and emissions reduction. Considerable efforts have been devoted to an in-depth understanding of the hydrodynamics mechanism and the prediction of ship behavior in waves. In addition to safety, efficiency and operability, the waves’ loads may determine structural failure. Prediction of seakeeping performances with a certain degree of accuracy is a demanding task for naval architects and of great practical interest for shipbuilders, owners, operators, as it affects both the ships’ design and operation [1]. More recently, the shipping industry has embraced more and more digitalization. Digital twin concept couples physical–numerical modelling and simulation to evolve solutions that improve vessel predictability, behavior control and response.
There are several aspects concerning seakeeping that make it one of the most challenging problems in ship hydrodynamics. The nonlinearities induced by the fluid viscosity and hydrodynamic pressure, boundary conditions at free-surface formulation, interference between external waves and ship’s body, specific geometry of the hull shape, all require advanced iterative time-domain procedures for the ship’s oscillations response analysis, involving significant computational resources [2].
The problem of a moving body interacting with waves has been pursued since the very first Froude and Michel studies. Since then, different approaches, such as experimental fluid dynamics, potential flow and recently computational fluid dynamics, have been developed to estimate the seakeeping performances. Traditionally, the hydrodynamic performances of a full-scale ship are determined by extrapolating the results of model-scale towing tank tests. Another approach commonly adopted in ship hydrodynamics is the employment of numerical techniques to solve the equations’ system of ship motions. The prediction of non-linear phenomena specific to ship motions in waves is difficult to be solved numerically, not only due to the complexity of the problem, but also the results are highly dependent on the details of the hull form and the incident wave condition. Due to the development of numerical modelling methods and computing power increase, direct prediction of full-scale ship performance became a practical approach. Most of the available techniques used to predict ship motions rely on assumptions from potential flow theory. Even if the natural trend in ship hydrodynamics is to move from frequency-domain to time-domain approach, from linear 2D strip theory type to fully 3D nonlinear techniques, and from potential to viscous computation, potential flow methods are still highly used in ship design, since they provide robust and quite accurate results in low to moderate sea states [3].
The frequency-domain approach is carried out mainly for linear or weakly nonlinear wave theories, where time dependence can be removed by assuming that solution is harmonic in time, in consequence only steady solutions are solved. Assuming that the ship body is slender, strip theory simplifies the 3D flow problem into a 2D formulation, modelling the ship hull as a set of multiple 2D ship stations [4]. The independent boundary value problem for each station can be solved analytically (Lewis form method) or by boundary element method. For the 2D hydrodynamic formulation, since the 1950s, various strip theories have been developed for the seakeeping problem [5]. The pioneer work of Korvin-Kroukovsky [6] set the principal feature of the strip theory for calculating ship motions based on the slender body assumption. This was the first suitable theory for numerical computations of ship motions that had adequate accuracy for engineering applications [7]. A modified strip theory approach of Gerritsma and Beukelman [8] was shown to obtain a good agreement with experimental tests for head wave case. Ogilvie and Tuck [9] developed a mathematically consistent approach based on short wave-length approximation, conducting a systematic analysis for the slender body problem to determine the added mass and damping for heave and pitch motions. Most of today’s strip methods are variations of the approach proposed by Salvesen et al. [7], which is one of the most complete versions, solving five degrees of freedom in ship’s motion equations, as the surge component was neglected, being not relevant for a slender body oscillation. Since then, more comprehensive strip theories have been developed for ship design, such as [10,11,12,13,14]. Few comprehensive reviews of strip theory variations have been reported by [4,15,16]. Most recently, a combination of two-dimensional strip theory with the two-dimensional Green function based on the potential theory to solve boundary values and motion responses of a semi-planing craft have been reported by [17]. Despite theoretical shortcomings, strip theory has the advantage of being fast, cheap and sufficiently accurate for a range of hull forms and moderate speeds. However, for higher speed vessels, highly flared hull forms, wave loads or extreme motions, other approaches have to be used.
The effort devoted to overcoming the weakness of the strip theory methods and to modelling the non-linear phenomena leads to the development of other potential flow alternatives based on numerical 3D methods that enhance the modelling technique of the physical domain of interest. By the late 1970s, a new treatment of the boundary conditions led to the Neumann–Kelvin approach [2], which supposes that the body boundary condition is applied to the mean position of the exact body surface and linearized free-surface boundary condition [2]. The most efficient way to solve the Neumann–Kelvin problem is to employ boundary integral methods where the solution is formulated in terms of singularities, as sources or dipoles, over the hull and free surface [2]. Two different categories of methods have been developed based on the type of singularities used in integral equations, one relies on Green function [18], which automatically satisfies both the radiation and the linearized free surface condition, and the other one relies on the distribution of the simple Rankine sources [18], not only on the hull surface, but also on the free surface. Some of the weaknesses of strip theory and Green function approach can be overcome by using a 3D time-domain panel method based on Rankine sources, including the fully 3D effects of the flow and forward speeds effects.
Hess and Smith [19] proposed in 1964 a method for evaluating velocity and pressure fields around a fully immersed arbitrary three-dimensional body. The method, also called the panel method or the boundary element method allowed, for the first time, to calculate the flow around an arbitrary three-dimensional body imposing the boundary condition exactly on the body surface. Firstly Gadd [20], then Dawson [21], applied the panel method for steady free surface flow distributing Rankine sources on hull and undisturbed free surface. A wide variety of different approaches for solving the nonlinear steady forward motion problem, based on Rankine sources distribution, has been reported by [22,23,24,25]. Due to the development of computing capabilities, the panel method has been continuously improved to solve the problems of ship motions in waves considering different modelling of nonlinearities of the hull and free surface boundary conditions. Nakos et al. [26] used a Rankine source method to solve transient wave–body interactions. This method has linearized the solutions about the double body flow. Large ship motion analysis performed by [27] was based on the application of the desingularised source method. Söding et al. [28] applied patch method instead of the collocation method to satisfy boundary conditions on the solid body surface, considering for the free surface a nonlinear condition. Recently, an improved desingularised Rankine method has been proposed by Mei et al. [29] to solve 3D diffraction and radiation problems. Dai and Wu [30] combined a method based on a Rankine panel method for the near-field and a transient Green function method for the far-field to predict large amplitude ship motions.
Potential flow solvers for seakeeping prediction have been continuously improved in terms of efficiency and accuracy, but still, the effect of viscosity, turbulence, wave dispersion, wave breaking, green water and slamming impact loading, deck green sea cannot be captured properly using the method based on potential flow assumption [31,32]. Considering the steadily increasing computational power and parallel computing, computational fluid dynamics methods based on the Reynolds-averaged Navier–Stokes (RANS) approach is more frequently applied to solve unsteady seakeeping problems. Most of the commercial solvers adopt a combination of unsteady RANS method, based on finite volume method and multi-phase flow approach, using the volume of fluid method for the treatment of free surface, which is rapidly gaining popularity for ship’s motions applications. The advantages of the nonlinear computation techniques without using analytical formulas for added resistance or empirical values for viscous effect are significant for the seakeeping analysis accuracy, but they are considerably more time consuming than potential flow approaches. Therefore, a good balance between accuracy and computational speed is required, especially for the ship design process.
The majority of RANS seakeeping simulations have been performed at model scale. The first attempt to solve ship motions in waves was presented in [33]. The results show some problems regarding the accuracy of free-surface due to the limited grid quality. Later, Simonsen et al. [34] carried out motions analysis in heave and pitch motions in regular head waves for KRISO container ship (KCS) hull using the CFDSHIP-IOWA code. Recently, Lungu and Bekhit [35,36] have assessed the seakeeping performances of the KCS and KVLCC (KRISO Very Large Crude Carrier) ship models for regular head wave conditions using NUMECA/FineMarine code. Their results show good agreement with the experiment for model scale. On the other hand, the scale effect seems to be important, as Hochkirch and Mallol [37] presented significant scale effects due to the differences between model-scale flows and full-scale flows. Tezdogan et al. [38] investigate seakeeping behavior and performance of the KCS model at a slow forward speed using Star-CCM+ package (Siemens PLM Software, Plano, Texas, United States).
Considering the context described above, the present work is focused on a systematic comparison study concerning the analysis of ship’s pitch, heave and roll motions in regular waves, for a full-scale ship model. The study was made using three different numerical approaches: the linear strip theory based on the Lewis form shape parameterization developed by Domnisoru [39], as in-house code DYN (module OSC - Oscillations), the non-linear potential flow based on Rankine source distribution using the commercial code SHIPFLOW and, for a limited number of cases, due to it being a highly time consuming solution, the viscous RANS method using NUMECA/FineMarine commercial code. Several series of computations have been performed for the regular wave height of 2 m, heading angle in the range of 0 to 180 deg, with a step of 45 deg. For the BEM model also, a wave height of 8 m has been considered. The computations have been carried out for the full-scale KCS hull numerical model, considering two forward speeds, 12 and, respectively, 24 knots.
The authors of this article aim to integrate all the presented methods in the same study in order to achieve a methodology that serves the rapid decision-making process in the case of further studies concerning the choice of the most accurate estimation method for ship motions on waves, depending on the specific requirements of the study and taking into account costs and calculation time. The proposed methodology meant approaching the mentioned methods, globally, with in-house code and commercial software, and the major advantage consisted of a better management of the input data, working conditions, calculation hypotheses, analyzed cases. It is worth mentioning that in the present article, full-scaled ship computations are considered, while most of the literature references reported analyses on model-scaled ship.

2. Methods for Ship Motions Prediction

As most of the available techniques to predict ship motions for design purposes rely on assumptions of potential flow theory, two different methods from this category have been chosen for comparison. The third approach considered for the present study is based on RANS solver, but, despite the parallel computing, this method does not prove yet to be able to solve the seakeeping problem in reasonable time to be consider in the practical design procedures. Taking into account that the methods involved in the present study rely on well-known mathematical formulation presented in seakeeping references [4,7,22,40,41,42], in this chapter only a brief presentation for each method is included.
The first method is based on a potential flow hydrodynamic linear 2D strip theory, with Lewis ship’s stations parameterization and regular wave Airy model excitation, with a frequency-domain solution of the coupled heave–pitch motion equations and uncoupled roll motion equation. This method is implemented by Domnisoru [39]—as in-house code DYN, module OSC—for the linear oscillations response amplitude operator (RAO) response amplitude operator computation. The mathematical formulation of the method used is described in [39]. The DYN code, module OSC, has been validated by experimental tests on scaled models of a fishing vessel [43] and a survey vessel [44], on follow, head, beam and quartering regular wave conditions.
For the second set of computations, a 3D fully nonlinear time-domain boundary element method, implemented into the commercial code SHIPFLOW [22,40] has been employed for ship motion estimation, considering only potential flow formulation being suitable for practical use. The method solves the missed boundary value problem for Laplace equation by distributing sources on the hull and free surface. Integrating the kinematic and dynamic boundary conditions with a fourth order Adam–Bashford–Moulton method, the elevation of free surface is obtained at each time step based on a Mixed Euler–Lagrange method. A blending zone based on an analytical solution is introduced in order to avoid the reflection from domain boundaries and in order to generate a certain incident wave. Once the velocity potential is obtained on the hull, using Unsteady Bernoulli formulation, the pressure distribution on the hull can be obtained. Finally, after integrating pressure on the hull, forces, moment and then motions are determined. A full description of the mathematical model and the numerical method may be found in [40,45]. Validation of the above-described method has been reported by Larsson et al. [46] for added resistance, pitch and heave, and by Coslovich et al. [47] mainly for roll.
The commercial software NUMECA/FineMarine [41,42] has been used in this study to solve the incompressible unsteady RANSE flow at ship motions in head waves. The solver relies on the finite volume method to build the spatial discretization for the governing equation. Closure to the turbulence is achieved by making use of the k-ω SST model with wall function formulation [41]. Pressure–velocity coupling is enforced through a SIMPLE (Semi-Implicit Method for Pressure Linked Equations)-like approach, where the velocity updates come from the momentum equation and the pressure is extracted from the mass conservation constraints transformed into pressure equation. Convection and diffusion terms in the RANSE are discretized using a second-order up-wind scheme and a central difference scheme, respectively. Free-surface capturing strategy is based on multi-phase flow approach using a volume of fluid method with high-resolution interface schemes [42]. Ship’s 6 DOF (Degrees Of Freedom) can be solved by the solver, but also some degree of freedom can be restrained. Mathematical model and the numerical framework are nicely described in [41]. Several validation studies of the solver have been reported by [48] for ship resistance and self-propulsion, and by [31,35,36,49] for ship motions in head waves. Their results show good agreement with the experiment for the scaled model.
In the preliminary ship design process, time is an important aspect to be considered due to the fact that naval architect usually cannot afford too long computational time. Considering the present methods comparison for one case (v = 12 Kn, ω = 0.4, hw = 2 m and head wave), by far strip theory is the less time-consuming method, the time spend for a single calculation case was 3 min using a single processor. The nonlinear BEM computations have been performed on 4 processor machines and the physical computational time for the same case was 19 h and 5 min for a number of about 68,000 panels. The RANS computations have been carried out on a high-performance computing (HPC) machine using 120 processors for a grid of about 20 million cells, which led to a 41 h and 47 min computational time. On the other hand, time consuming should be also correlated with the accuracy and with the amount of information received from a specific flow post-processing. The RAO functions for heave, pitch and roll based on the linear strip theory implemented in DYN code have been validated in references [43,44] for heading angles, 0, 90, 180 degrees, proving good results for practical investigations, revealing average differences from 12% to 18%. A validation BEM method in terms of heave and pitch was performed by the authors in [50] and the computed results showed good agreement of 3% to 10% with towing tank measurements. As mentioned before, extensive validation studies may be found in [45,46]. In addition, grid convergence and validations studies based on FineMarine RANS solver have been reported by Lungu and Bekhit [35,36]. Their results revealed good agreement with the towing tank test measurements for model scale, 1–7%. Taking all the above mentioned into consideration, a trade-off between accuracy and computational costs must be counted in order to obtain feasible results for the ship design practical use.
In order to compare directly the computational fluid dynamics (CFD) results with those obtained by using linear strip theory, for all numerical analyses, the regular Airy model is considered as excitation source. Both CFD methods are hydrodynamic nonlinear, so that even if the excitation has one harmonic, the time-domain response has several harmonics. For the two CFD methods, heave, pitch and roll time records are spectrally analyzed by using Direct Fourier Transformation Method. Regarding the response amplitude spectra (see Section 6), only the motions amplitudes with encountering wave frequency harmonic are considered. The response amplitude operators (RAO s) for ship motions obtained by the two CFD methods represent the ratio between heave, pitch and roll amplitudes on wave’s harmonic and the regular wave amplitude.

3. Hull Geometry and Conditions

The objective of the present study is to investigate the seakeeping performances of the KRISO container ship (KCS) while sailing in regular wave at two speeds, 12 and 24 Kn, corresponding to Froude numbers of 0.13, and 0.26, respectively. The KCS hull is a benchmark test case for ship hydrodynamics widely used in the marine hydrodynamic scientific community [34]. The hull is a modern commercial container vessel with bulbous bow and flare fore above the waterline. The stern is a typical pram type stern with transom. Bare hull, without rudder, has been considered for the seakeeping calculations. KCS hull geometry is presented in Figure 1 and the main dimensions and conditions used for numerical simulations are presented in Table 1.
In Table 2, the test cases for seakeeping analysis are presented. The ship motions are approached considering the ship on regular Airy wave, having the wave height of 2 and 8 m, the wave frequency between 0 and 2 rad/s and heading angle in range of 0 to 180 deg, with a step of 45 deg. Even though the ship response in a regular wave at zero speed was computed, in the following paragraphs the ship motions for only two speed values of 12 and 24 Kn, respectively, are presented. Although, the study was carried out for three values of vertical center of gravity, KG = 7.3, 10.8 and 14.3 m, in the paper, the results for only KG = 10.8 m, which corresponds to a medium loading case of the ship, are presented.

4. Numerical Setup

Several details about the computational conditions will be briefly described in the following. For the discretization of the hull used for the linear 2D strip theory calculations, a model of 526 stations has been used.
For the BEM method, the boundary value problem modelling the field equation and boundary conditions have been solved, making use of the commercial code SHIPFLOW. The computations carried out with BEM is a fully nonlinear unsteady three-dimensional potential flow method, that takes into account the geometric non-linearities of ship shape and the induced hydrodynamic non-linearities, although the external excitation is kept linear. Wave reflection from domain’s boundaries is a pure numerical problem that affects all the simulations with a restricted domain. In order to avoid wave reflection from the free surface truncation boundaries, the flow on the intersection between the outer and inner domain has to be the same. To achieve this condition, a numerical damping zone is located close to the domain’s boundaries and such a zone only has to dampen the difference between the velocity potential in the inner and outer domain. This difference in the solution between the inner and the outer domain is mainly due to the presence of the body with radiation and diffraction effects but can also derive from a numerical error. This damping is achieved adding a term in the free surface boundary condition for all those panels that belong to such a zone [51,52]. The fluid domain assumes that the flow is known a priori in the outer domain and the current method uses analytically described waves to represent the flow field in the outer domain. The hull and free surface have been discretized by panels. The calculations have been performed for a surface domain of one ship length upstream, two ship length downstream, the width of the free surface being two ship length. For BEM method, the ship hull is discretized both sides without the request of symmetry condition, so that any heading angle condition can be computed. The number of panels distributed on hull and free surface varied between 60,000 and 90,000 panels function of the wave frequency case studied, which assures a minimum number of 40 panels per wavelength. The computations have been performed on local desktop machines with four cores of 3.1 GHz.
For the unsteady RANS simulation of incompressible flow of ship motions in head waves, the commercial solver NUMECA/FineMarine was used. A reference length has been considered Lref = max (LBP, λ). The dimensions of the computation domain have been chosen according to the ship and wavelength, considering 2.0 Lref upstream, 4.0 Lref downstream, 2.0 LBP on the side, 4.0 LBP underneath, and 2.0 LBP above the undisturbed free-surface level. The boundary conditions imposed on the solid wall and domain boundaries, but also the dimensions considered for the computational domain, are depicted in Figure 2. A mirror condition is applied on the centerline plane of the ship and on the side boundary of the domain, to avoid wave reflections at lateral boundary. Considering the head wave case computed by the RANS method, the corresponding wave is generated at the upstream boundary and the inlet boundary is chosen to ensure at least two full waves to be generated before encountering the ship. The free surface is refined for the entire domain using three different refinement boxes, as can be observed in Figure 2. The first starts from the inlet boundary and is extended until it reaches 1.0 Lref behind the ship. The cell size is chosen for this refinement zone in x- and y-direction to provide minimum 60 cells per wavelength, such that Δx1 = Δy1 = λ/60, while in the z-direction the refinement criterion is selected as Δz = hw/16, where hw represents the wave height. The refinement depth in z-direction is extended for 3.0 hw equally distributed above and beneath the non-disturbed free-surface level, which is set at the design draft T = 10.8 m from the ship base line. The second and third boxes are coarsened gradually in x- and y-directions to generate sufficient numerical damping zones to prevent any reflections from the exit boundary. The second box is extended for 1.0 Lref and it is used as an initial damping relaxation zone, in which the cell size is coarsened in x- and y-directions by a factor of 4 (i.e., Δx2 = Δy2 = 4.0 × Δx1). The third refinement zone is used as a final damping zone, having the cell size coarsened in x- and y-directions by a factor of 8.0 (i.e., Δx3 = Δy3 = 8.0 Δx1), while the refinement in the z-direction for the entire domain is maintained unchanged. The effect of damping zone on the RANS solution is clearly seen in Figure 3.
Body-fitted full hexahedral unstructured meshes of about 22 million cells have been generated (Figure 4). The time step Δt has been chosen in order to fulfil the condition of 300-time steps per wave period, with fourth order convergence criteria. All the computations have been performed on a HPC machine with 120 cores of 3.3 GHz.
In general, methodology used in the present paper for the seakeeping performance prediction is based on the practical guidelines recommended by the International Towing Tank Conference (ITTC), by software providers and different sensitivity and validation studies [35,36,49,50,51,52], but for some particular cases the recommendations could not be accomplished due to the practical reasons.

5. Grid Convergence Test

A grid study for SHIPFLOW calculations was conducted for three grids based on the Richardson extrapolation [53] in order to investigate the numerical simulation error and uncertainty. Two calculation cases have been considered for the grid test: ship on regular waves, head waves (μ = 180 deg) and beam sea (μ = 180 deg). For both cases ship speed was v = 12 Kn. The results of grid convergence study for RAO heave, pitch and roll are summarized in Table 3. The convergence ratio is defined as RG = ε2132, where ε21 = S2S1 and ε32 = S3S2 stands for simulation error, S1 represents the solution calculated for fine grid (68,716 panels), S2 for medium grid (32,692 panels) and S3 for coarse grid (16,346 panels), according to Richardson extrapolation approach [53], with the expressions (1) ÷ (4). The grid corresponding to each level of refinement is depicted in Figure 5 for hull and free surface.
Considering the computed values of RG for all variables of grid convergence study, one can observe that monotonic convergence is achieved for all the simulation variable. The order of accuracy pG can be calculated based on the refinement ration rG, as follows:
p G = ln ε 32 ε 21 ln r G .
The calculated order of accuracy is used further to determine the correction factor CG and the error δG, as follows:
c G = r G p G 1 r G p t h 1 ,
δ G = ε 21 r G p G ,
where pth is the theoretical grid convergence order usually considered pth = 2. Finally, the uncertainty is computed by the expression:
U G = C G × δ G + 1 C G × δ G
Based on this study, bearing in mind that the percentage estimated error with respect to the fine grid is very small, i.e., ε21%S1 << 1, which concludes that the solution computed by SHIPFLOW is grid independent.

6. Results and Discussions

In this section, the ship hydrodynamic response in regular waves is presented and discussed. The results are represented in the frequency domain, in terms of response amplitude operators (RAO) for heave, pitch and roll motion, respectively. The charts contain the results obtained by linear strip theory method (DYN code) and the boundary element method (SHIPFLOW software), also including a few results obtained using the RANS method (NUMECA/FineMarine).
In Figure 6, the ship on regular waves, for the case of top speed v = 12 Kn, μ = 90 deg, ω = 0.4 rad/s and hw = 2 m, is presented. The images represent the hull position at five different fractions of ship motion periods (ΔT) related to free surface topologies. Hydrodynamic forces and moments are determined by integrating pressure over the ship hull based on the quadratic pressure distribution. Pictures on the left, Figure 6a1–e1, present the pressure distribution over the ship hull. One may see that the flow solution for ω = 0.4 rad/s reveals that the wave is fully developed, and artificial damping condition imposed on the boundaries works fine. Moreover, the interference of the wave system generated by the ship and the incident wave is well captured by the boundary element method and some nonlinear effects are expected (see the pictures on the right, Figure 6a2–e2).
In Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11, the response amplitude operators for heave, pitch and roll motions of KCS ship advancing in regular waves are represented. Heave RAO (m/m), Pitch RAO (rad/m) and Roll RAO (rad/m), respectively, for the speed case of ship speed 12 Kn and the wave height 2 m (amplitude aw = 1 m). In Figure 7, the heave response amplitude operator, Heave RAO (m/m), for the speed case of 12 Kn, obtained using a 2D strip and BEM methods is represented. As a general remark, for both methods, the ship response is significant in the lower frequency domain, while, as the frequency increases, the ship response decreases considerably. The maximum values are recorded for the heading angle μ = 90 deg, at ω = 0.756 rad/s. Once the circular frequency of the regular wave increases, one may observe that both methods are in good agreement.
In Figure 8, the response amplitude operator for pitch motion, Pitch RAO (rad/m) is depicted, for ship speed 12 Kn, heading angles from 0 to 180 deg, and wave height of 2 m. For pitch motion at beam sea condition, the response amplitude is the smallest one, recording scattered values between the 2D strip and the BEM hydrodynamic models. In particular, at 90 deg, the variation of BEM results has two peaks, the second one at similar frequency as the DYN results; however, the values obtained by the boundary element method are higher than the strip model results at beam sea. The maximum value is recorded at μ = 135 deg for ω = 0.57 rad/s. A good agreement between both hydrodynamic models is found at a heading angle of 45 and 135 deg.
In Figure 9 and Figure 10, the following are presented: numerical results for Heave RAO (m/m) and Pitch RAO (rad/m) obtained with linear 2D strip method (DYN), BEM boundary element method (SHIPFLOW) and RANS method (NUMECA/FineMarine), considering the ship on head regular wave, µ = 180 deg, with 2 m height and a ship speed of 12 Kn.
In Figure 11, the following are given: the results for the roll motion in terms of response amplitude operator, Roll RAO (rad/m). The results are represented for the heading angle of 45, 90 and 135 deg, respectively. The maximum value for Roll RAO is recorded for μ = 90 deg and, as long as it moves away from the beam waves, the roll motion of the ship attenuates.
Next, a comparison between the results computed with linear 2D strip and BEM hydrodynamic models, for a regular wave height hw = 8 m (aw = 4 m) for BEM model, is depicted in Figure 12 and Figure 13. First, one may see the Heave RAO (m/m) values at five different heading angles (Figure 12). The trend of the differences between the two methods remains the same as for the case of wave height hw = 2 m. Very reduced differences occur in the case of the BEM method results at 2 and 8 m wave height, due to the reduced hydrodynamic nonlinearities. Higher differences appear at beam wave condition between 0.4 and 0.6 rad/s.
In Figure 13, the variation curves of Pitch RAO (rad/m) for wave height hw = 8 m for BEM model, with the peak value for the heading angle of µ = 135 deg around ω = 0.5 rad/s, are depicted. At beam waves, µ = 90 deg, differences between the hydrodynamic approaches are scattered, where the values of Pitch RAO are very reduced.
The next results are computed for the second ship speed case, v = 24 Kn, for both wave heights of 2 and 8 m, respectively. There were also considered the same medium loading conditions of the ship. In Figure 14, the RAO heave (m/m) is presented. The differences between the linear 2D strip and BEM hydrodynamic approaches have slightly increased compared to 12 Kn ship speed, especially for wave circular frequency lower than ω = 0.4 rad/s.
In Figure 15, the Pitch RAO (rad/m) for 24 Kn ship speed and 2 m wave height are presented. The Pitch RAO has the maximum peak values for the head wave case, around ω = 0.46 rad/s. At the same time, one may observe that the hydrodynamic models provide better agreement for the beam wave case µ = 90 deg.
In Figure 16, the Roll RAO (rad/m) for 24 Kn ship speed and 2 m wave height are presented. The trend of results is maintained for the roll motion as for the 12 Kn ship speed case, the boundary element method being in agreement with the strip theory, for all the heading waves considered. The maximum values appear at beam wave, 0.16 rad/m around ω = 0.47 rad/s. At the same time, as long as the ship is in oblique waves, the values of roll motion decrease.
At 24 Kn ship speed and 8 m wave height for the BEM model, the Heave RAO (Figure 17) records visible differences to the BEM solution at the same speed, but with 2 m wave height, due to hydrodynamic nonlinearities for the wave’s circular frequency up to 0.6 rad/s. Figure 18 presents the pitch RAOs for 24 Kn ship speed and 8 m wave height. The differences between the BEM model, with 8 m wave height, and the linear 2D strip model become more significant in comparison to the case with wave height 2 m, due to the recorded hydrodynamic nonlinearities.
The Heave RAO peak is 1.071 m/m around ω = 0.5 rad/s for μ = 135 deg (as can be seen in Figure 17) while the maximum value of the Pitch RAO is 0.02 rad/m around ω = 0.45 rad/s (shown in Figure 18), when the ship is in head waves μ = 180 deg.
In Figure 19, a comparison between heave amplitude spectrum, by DFT (Direct Fourier Transformation) method implemented in DYN code, for the two speed cases (12 and 24 Kn) and wave height hw = 2 m, heading angle µ = 180 deg, wave circular frequency ω = 0.4 rad/s (with ship–wave fe encountering frequencies 0.080 and 0.096 Hz), is made. The same comparison for pitch amplitude spectrum, by DFT method is made. Both diagrams represent the results obtained using the BEM boundary element method by SHIPFLOW. In Figure 20, the same comparison is made for wave height of hw = 8 m.
Both heave and pitch amplitude spectra by DFT method, in the case of higher ship speed of 24 Kn, have narrowed frequency bands. Additionally, at hw = 8 m wave height, the heave and pitch amplitude spectra present a series of higher harmonics, with small amplitudes, due to the hydrodynamic nonlinearities.
In Figure 21 and Figure 22, a comparison between BEM and RANS hydrodynamic models is performed. A comparison of Heave Time Record and Pitch Time Record computed for the KCS vessel at v = 12 Kn at heading angle μ = 180 deg and circular wave frequency ω = 0.4 rad/s (fe = 0.080 Hz) is presented in Figure 21.
A comparison of Heave Amplitude Spectrum and Pitch Amplitude Spectrum computed for the KCS vessel at v = 12 Kn, at heading angle μ = 180 deg and circular wave frequency ω = 0.4 rad/s (fe = 0.080 Hz) is presented in Figure 22. While for the Heave Amplitude Spectrum, SHIPFLOW and NUMECA show good agreement, for the Pitch Amplitude Spectrum some difference occurs between the two models.

7. Conclusions

In this paper, for the KCS full-scale model (Section 3), several numerical seakeeping analyses have been performed, in order to compare three hydrodynamic methods for the prediction of the dynamic response on heave, pitch and roll oscillations, in terms of response amplitude operators. The selected hydrodynamic methods are the most used approaches for ships motions, linear 2D strip theory, BEM and RANS methods, involving in-house codes and commercial codes.
A grid convergence study has been conducted to determine the accuracy of the grid used to compute the flow based on the potential BEM approach (Section 5). The study revealed that solution computed by SHIPFLOW is grid independent and the grid error for the grid used is less than 1%.
The BEM and RANS methods, with time-domain solutions, were compared to the linear strip (ST) method, with frequency-domain solution; although the external excitation is a regular wave, the ship’s response shows higher harmonics (Figure 19 and Figure 20), due to the hydrodynamic nonlinearities. Therefore, in order to ensure the same reference for all three methods, the nonlinear response amplitude operators (RAOs), obtained using BEM and RANS methods, were computed considering only the main spectral amplitude for the ship-wave encountering frequency, being the only result delivered by the linear 2D strip method. Nevertheless, the BEM and RANS studies have revealed, for the KCS model, reduced nonlinear dynamic response components, even if the regular wave height was up to 8 m and the speed of the ship was 24 Kn (Figure 19, Figure 20 and Figure 22).
For a detailed benchmark between the three hydrodynamic models (Section 2), this parametric study developed for the KCS model (Section 3 and Section 4), at one loading case, has been focused on the influence of ships speed (12 and 24 Kn), wave heading angle (0, 45, 90, 135, 180 deg) and height (2 and 8 m), on the ship’s dynamic response in regular wave, using deterministic analyses (Section 6).
Comparing the Heave RAO functions between the three hydrodynamic models’ results (Figure 7, Figure 9, Figure 12, Figure 14 and Figure 17) (Section 6), for ship speed of 24 Kn versus 12 Kn and 8 m wave height versus 2 m, the BEM method delivers visible differences due to the hydrodynamic nonlinearities, especially for the wave circular frequencies up to 0.6 rad/s. In the case of beam waves, for wave circular frequency within the range of 0.4 to 0.6 rad/s, the BEM values are quite scattered around the 2D strip model results. In addition, visible differences, for head waves, have been seen between the RANS method and the other two (Figure 9).
Comparing the Pitch RAO functions between the three hydrodynamic models (Figure 8, Figure 10, Figure 13, Figure 15 and Figure 18) (Section 6), it turned out that the maximum values are obtained for head wave conditions and are very small for beam wave conditions. As the ship’s speed and wave height increase, the hydrodynamic nonlinearities induce visible differences between the results, especially close to the peak values of Pitch RAO. Except for 0.3 rad/s, on head waves, the differences between the RANS and the other two methods are quite scattered for the analyzed wave’s frequency domain (Figure 10).
Comparing the Roll RAO functions between the BEM and 2D strip hydrodynamic results (Figure 11 and Figure 16) (Section 6), it turned out that the maximum values are obtained for beam waves condition, being a narrow frequency band response amplitude operator. The differences are noticeable at a beam wave condition close to the RAO peak value, especially due to the different hydrodynamic damping formulation of the two methods.
As an overall analysis, we can conclude that the RAO’s values predicted by the linear 2D strip theory can be considered as average values compared to the BEM and RANS results, being suitable for practical seakeeping analyses of mono-hull slender body ships.
Further studies will extend this parametric analysis to more numerical results based on the RANS method considering other ship types, being focused on the sensitivity of the hydrodynamic methods used for ship’s motions prediction. Additionally, the next step consists of short-term seakeeping analysis in irregular waves, where, besides the BEM and RANS methods, also a nonlinear strip theory with time-domain solution is to be applied.

Author Contributions

Conceptualization, F.P. and L.D.; methodology, F.P., L.D. and S.P.; software, F.P., L.D. and S.P.; validation, F.P. and L.D.; formal analysis, S.P.; investigation, F.P., L.D. and S.P.; resources, F.P., L.D. and S.P.; data curation, F.P.; writing—original draft preparation, F.P., L.D. and S.P.; writing—review and editing, S.P., F.P. and L.D.; visualization, S.P. and F.P.; supervision, L.D.; project administration, S.P.; funding acquisition, S.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The research study presented in this article was developed in the frame of Dunarea de Jos University of Galati, Naval Architecture Research Centre.

Conflicts of Interest

The author declares no conflict of interest.

Abbreviations

KCSKRISO container ship
KRISOKorea Research Institute of Ships and Ocean Engineering
IMOInternational Maritime Organization
DYNDynamic Ship Analysis
RANSEReynolds-averaged Navier–Stokes equation
BEMBoundary element method
STStrip theory
EEOIEnergy efficiency operational indicator
CFDComputational fluid dynamics
RAOResponse amplitude operator
KVLCCKRISO Very Large Crude Carrier
OSCOscillations
DOFDegrees of freedom
DFTDirect Fourier Transformation
Nomenclature
k-ω SSTk-ω turbulence model, SST shear stress transport
LrefReference length
LBPLength between the perpendiculars
LWLLength of waterline
BWLBeam at waterline
DDepth
TDesign draft
ΔDisplacement
CBBlock coefficient
vShip’s speed
LCBLongitudinal center of buoyancy (%LBP), fwd+
LCGLongitudinal center of gravity from the aft peak
KGVertical center of gravity from keel
Kxx/BMoment of inertia
Kyy/LBP, Kzz/LBPMoment of inertia
SShip wetted area
ω Wave circular frequency
feShip-wave encountering frequency
h w Wave height
a w Wave amplitude
λWavelength
μ Heading angle
εSimulation error
RGConvergence ratio
SiSolution calculated for fine (S1), medium (S2) and coarse grid (S3)
rGRefinement ration
pGOrder of accuracy
CGCorrection factor
δGError
pthTheoretical grid convergence order

References

  1. Shukui, L.; Papanikolau, A. On Nonlinear Simulation Methods and Tools for Evaluating the performance of Ships and Offshore Structures in Waves. J. Appl. Math. 2012, 2012, 1–21. [Google Scholar]
  2. Beck, R. Modern Seakeeping Computations for Ships. In Twenty-Third Symposium on Naval Hydrodynamics; The National Academies Press: Washington, DC, USA, 2001. [Google Scholar] [CrossRef]
  3. ITTC, The Seakeeping Committee. Final report and recommendations to the 25th ITTC. In Proceedings of the 25th ITTC, Fukuoka, Japan, 14–20 September 2008; Volume I. [Google Scholar]
  4. Bertram, V. Practical Ship Hydrodynamics, 2nd ed.; Butterworth-Heinemann: Oxford, UK, 2012. [Google Scholar]
  5. Sclavounos, P.D.; Lee, C. Topics on Boundary Element Solutions of Wave Radiation-Diffraction Problems. In Proceedings of the 4th International Conference on Numerical Ship Hydrodynamics, Washington, DC, USA, 24–27 September 1988. [Google Scholar]
  6. Korvin-Kroukovsky, B. Investigation of ship motions in regular waves. Trans. SNAME 1955, 63, 386–435. [Google Scholar]
  7. Salvesen, N.; Tuck, E.O.; Faltinsen, O.M. Ship motions and sealoads. Trans. SNAME 1970, 78, 250–287. [Google Scholar]
  8. Gerritsma, J.; Beukelman, W. Analysis of the Modified Strip Theory for Calculation of Ship Motions and Wave Bending Moments. Int. Shipbuilding Prog. 1967, 14, 319–337. [Google Scholar] [CrossRef]
  9. Ogilvie, T.F.; Tuck, E.O. Rational Strip Theory of Ship Motions; University of Michigan: Ann Arbor, MI, USA, 1969. [Google Scholar]
  10. Holloway, D.; Davis, M. Ship Motion Computations Using a High Froude Number Time Domain Strip Theory. J. Ship Res. 2006, 50, 15–30. [Google Scholar]
  11. Jensen, J.J.; Pedersen, P.T.; Shi, B.; Wang, S.; Petricic, M.; Mansour, A.E. Wave Induced Extreme Hull Girder Loads on Containership; The Society of Naval Architects and Marine Engineers Transactions: Jersey City, NJ, USA, 2008; Volume 116. [Google Scholar]
  12. Sutulo, S.; Guedes Soares, C. A Generalized Strip Theory for Curvilinear Motion in Waves. In Proceedings of the 27th International Conference on Offshore Mechanics and Arctic Engineering, Estoril, Portugal, 15–20 June 2008; pp. 359–368. [Google Scholar]
  13. Prpic-Orsic, J.; Faltinsen, O.M. Estimation of ship speed loss and associated CO2 emissions in a seaway. Ocean Eng. 2012, 44, 1–10. [Google Scholar] [CrossRef]
  14. Bandyk, P.J.; Hazen, G.S. A forward-speed body-exact strip theory. In Proceedings of the 34th International Conference on Ocean, Offshore and Arctic Engineering, St. John’s, NL, Canada, 31 May–5 June 2015. [Google Scholar]
  15. Newman, J.N. Marine Hydrodynamics, 1st ed.; MIT Press: Cambridge, MA, USA, 1977. [Google Scholar]
  16. Faltinsen, O.M. Hydrodynamics of High-Speed Marine Vehicles, 1st ed.; Cambridge University Press: Cambridge, UK, 2005. [Google Scholar]
  17. Lin, H.; Lin, C.W. Numerical Simulation of Seakeeping Performance on the Preliminary Design of a Semi-Planing Craft. J. Mar. Sci. Eng. 2019, 7, 199. [Google Scholar] [CrossRef] [Green Version]
  18. Rorvik, J. Application of Inviscid Flow CFD for Prediction of Motions and Added Resistance of Ships. Master’s Thesis, NTNU Trondheim, Trondheim, Norway, 2016. [Google Scholar]
  19. Hess, J.L.; Smith, M.O. Calculation of Nonlifting Potential Flow about Arbitrary Three-dimensional body. J. Ship Res. 1964, 8, 22–44. [Google Scholar]
  20. Gadd, G.E. A Method of Computing the Flow and Surface Wave Pattern around Full Forms; The Royal Institution of Naval Architects: London, UK, 1976; pp. 207–219. [Google Scholar]
  21. Dawson, C.W. A Practical Computer Method for Solving Ship-Wave Problem. In Proceedings of the 2nd International Conference on Numerical Ship Hydrodynamics, Berkeley, CA, USA, 19–21 September 1977; pp. 30–38. [Google Scholar]
  22. Jensen, G.; Mi, Z.X.; Söding, H. Rankine Source Methods for Numerical Solutions of the Steady Wave Resistance Problem. In Proceedings of the 16th Symposium on Naval Hydrodynamics, Berkeley, CA, USA, 13–18 July 1986. [Google Scholar]
  23. Raven, H.C. Variations on a Theme by Dawson. In Proceedings of the 17th Symposium on Naval Hydrodynamics, The Hague, The Netherlands, 29 August–2 September 1988; pp. 151–172. [Google Scholar]
  24. Sclavounos, P.D.; Nakos, D.E. Stability Analysis of Panel Methods for Free Surface Flow with forward Speed. In Proceedings of the 17th Symposium on Naval Hydrodynamics, The Hague, The Netherlands, 29 August–2 September 1988. [Google Scholar]
  25. Janson, C.E. Potential Flow Panel Methods for the Calculation of Free-Surface Flows with Lift. Ph.D. Thesis, Chalmers University of Technology, Gothenburg, Sweden, 1997. [Google Scholar]
  26. Nakos, D.E.; Kring, D.E.; Sclavounos, P.D. Rankine Panel Method for Time Domain Free Surface Flows. In Proceedings of the 6th International Conference on Numerical Ship Hydrodynamics, Iowa City, IA, USA, 2–5 August 1993. [Google Scholar]
  27. Zhang, X.; Bandyk, P.; Beck, R.F. Large Amplitude Body Motion Computations in the Time-Domain. In Proceedings of the 9th International Conference in Numerical Ship Hydrodynamics, Ann Arbor, MI, USA, 5–8 August 2007. [Google Scholar]
  28. Söding, H.; von Gräfe, A.; el Moctar, O.; Shigunov, V. Rankine Source Method for Seakeeping Predictions. In Proceedings of the 31st International Conference on Ocean, Offshore and Arctic Engineering, Rio de Janeiro, Brazil, 1–6 July 2012. [Google Scholar]
  29. Mei, T.; Candries, M.; Lataire, E.; Zou, Z. Numerical Study on Hydrodynamics of Ships with Forward Speed Based on Nonlinear Steady Wave. J. Mar. Sci. Eng. 2020, 8, 106. [Google Scholar] [CrossRef] [Green Version]
  30. Dai, Y.Z.; Wu, G.X. Time Domain Computation of Large Amplitude Body Motion with the Mixed Source Formulation. In Proceedings of the 8th ICHD, Nantes, France, 30 September–3 October 2008; pp. 441–452. [Google Scholar]
  31. Deng, D.G.; Queutey, P.; Visonneau, M. RANS prediction of KVLCC2 tanker in head waves. In Proceedings of the 9th International Conference on Hydrodynamics, Shanghai, China, 11–15 October 2010. [Google Scholar]
  32. Choi, J.; Yoon, S.B. Numerical simulations using momentum source wave-maker applied to RANS equation model. Coast. Eng. 2009, 56, 1043–1060. [Google Scholar] [CrossRef]
  33. Sato, Y.; Miyata, H.; Sato, T.J. CFD simulation of 3-dimensional motion of a ship in waves: Application to an advancing ship in regular heading waves. J. Mar. Sci. Technol. 1999, 4, 108–116. [Google Scholar] [CrossRef]
  34. Simonsen, C.D.; Otzen, J.F.; Stern, F. EFD and CFD for KCS Heaving and Pitching in Regular Head Waves. J. Mar. Sc. Tech. 2013, 18, 435–459. [Google Scholar] [CrossRef]
  35. Lungu, A. Unsteady Numerical Simulation of the Behaviour of a Ship Moving in Head Sea. In Proceedings of the 38th International Conference on Ocean, Offshore and Arctic Engineering, Glasgow, UK, 9–14 June 2019. [Google Scholar]
  36. Beckhit, A.; Lungu, A. URANSE simulation for the Seakeeping of the KVLCC2 Ship Model in Short and Long Regular Head Waves. IOP Conf. Ser. Mater. Sci. Eng. 2019, 591, 012102. [Google Scholar] [CrossRef] [Green Version]
  37. Hochkirch, K.; Mallol, B. On the importance of full-scale CFD simulations for ships. In Proceedings of the 12th International Conference on Computer Applications and Information Technology in the Maritime Industries (COMPIT 2013), Cortona, Italy, 15–17 April 2013; pp. 85–95. [Google Scholar]
  38. Tezdogan, T.; Demirel, Y.K.; Kellett, P.; Khorasanchi, M.; Incecik, A.; Turan, O. Full-scale unsteady RANS CFD simulations of ship behaviour and performance in head seas due to slow steaming. Ocean Eng. 2015, 97, 186–206. [Google Scholar] [CrossRef] [Green Version]
  39. Domnisoru, L. Ship Dynamics. Oscillations and Vibrations; Technical Publishing House: Bucharest, Romania, 2001. [Google Scholar]
  40. Kjellberg, M.; Janson, C.E.; Contento, G. Fully Nonlinear Potential Flow Method for Three-Dimensional Body Motion Proceedings from NAV 2012. In Proceedings of the 17th International Conference on Ships and Shipping Research, Naples, Italy, 17–19 October 2012. [Google Scholar]
  41. Duvigneau, R.; Visonneau, M.; Deng, G.B. On the role played by turbulence closures in hull shape optimization at model and full-scale. J. Mar. Sci. Technol. 2003, 8, 1–25. [Google Scholar] [CrossRef]
  42. Queutey, P.; Visonneau, M. An interface capturing method for free-surface hydrodynamic flows. Comput. Fluids 2007, 36, 1481–1510. [Google Scholar] [CrossRef]
  43. Obreja, D.; Nabergoj, R.; Crudu, L.; Domnisoru, L. Seakeeping performance of a Mediterranean fishing vessel. In Proceedings of the IMAM Maritime Transportation and Harvesting of Sea, Lisbon, Portugal, 9–11 October 2017; pp. 483–491. [Google Scholar]
  44. Burlacu, E.; Domnisoru, L.; Obreja, D. Seakeeping Prediction of a Survey Vessel Operating in the Caspian Sea. In Proceedings of the 37th International Conference on Offshore Mechanics and Arctic Engineering, Madrid, Spain, 17–22 June 2018. [Google Scholar]
  45. Zhang, B.J.; Ning, X. The research of added resistance in waves based on nonlinear time-domain potential flow theory. J. Mar. Sci. Technol. 2018, 26, 343–351. [Google Scholar]
  46. Larsson, L.; Stern, F.; Visonneau, M. Numerical Ship Hydrodynamics—An Assessment of the Gothenburg 2010 Workshop; Springer: Dordrecht, The Netherlands, 2013. [Google Scholar]
  47. Coslovich, F.; Contento, G.; Kjellberg, M.; Janson, C.E. Computations of Roll Motions in Waves Using a Fully Nonlinear Potential Flow Method, Technologyand Science for the Ship of the Future; Marino, A., Bucci, V., Eds.; IOS press: Amsterdam, The Netherlands, 2018; pp. 186–193. [Google Scholar]
  48. Lungu, A. Numerical simulation of the resistance and self-propulsion model tests. J. Offshore Mech. Arctic Eng. 2020, 142, 021905. [Google Scholar] [CrossRef]
  49. Bekhit, A.; Lungu, A. Numerical Simulation for Predicting Ship Resistance and Vertical Motions in Regular Head Waves, ASME 2019. In Proceedings of the 38th International Conference on Ocean, Offshore and Arctic Engineering, Glasgow, Scotland, 9–14 June 2019. [Google Scholar]
  50. Pacuraru, S.; Domnisoru, L.; Pacuraru, F. Numerical study on motions of a containership on head waves. In Proceedings of the 6th International Conference—Modern Technologies in Industrial Engineering, Constanta, Romania, 13–16 June 2018; Volume 400. [Google Scholar]
  51. Kjellberg, M. Fully Nonlinear Unsteady Three-Dimensional Boundary Element Method for Ship Motions in Waves. Ph.D. Thesis, Chalmers University of Technology, Gothenburg, Sweden, 2013. [Google Scholar]
  52. Coslovich, F. Computations of Ship Motions in Waves Using a Fully Nonlinear Time Domain Potential Flow Method. Master’s Thesis, University of Trieste, Trieste, Italy, 2015. [Google Scholar]
  53. International Towing Tank Conference (ITTC), Uncertainty Analysis in CFD Verification and Validation Methodology and Procedures 7.5-03-01-01. 2017, pp. 1–13. Available online: https://www.ittc.info/media/8153/75-03-01-01.pdf (accessed on 20 August 2020).
Figure 1. KRISO container ship (KCS) hull geometry.
Figure 1. KRISO container ship (KCS) hull geometry.
Jmse 08 00962 g001
Figure 2. RANS Computational domain, dimensions and boundary conditions.
Figure 2. RANS Computational domain, dimensions and boundary conditions.
Jmse 08 00962 g002
Figure 3. RANS damping zone effect.
Figure 3. RANS damping zone effect.
Jmse 08 00962 g003
Figure 4. Computational grid. (a) Domain discretization; (b) mesh details in bow area.
Figure 4. Computational grid. (a) Domain discretization; (b) mesh details in bow area.
Jmse 08 00962 g004
Figure 5. Grid refinement levels used for the grid convergence test. (a1) Coarse grid on hull surface; (a2) Coarse grid on free-surface; (b1) Medium grid on hull surface; (b2) Medium grid on free-surface; (c1) Fine grid on hull surface; (c2) Fine grid on free-surface.
Figure 5. Grid refinement levels used for the grid convergence test. (a1) Coarse grid on hull surface; (a2) Coarse grid on free-surface; (b1) Medium grid on hull surface; (b2) Medium grid on free-surface; (c1) Fine grid on hull surface; (c2) Fine grid on free-surface.
Jmse 08 00962 g005
Figure 6. Pressure distribution and free surface topology computed for the case μ = 90 deg. (a1) Pressure distribution on the hull, ΔT = 0; (a2) Free-surface topology, ΔT = 0; (b1) Pressure distribution on the hull, ΔT = 0.125; (b2) Free-surface topology, ΔT = 0.125; (c1) Pressure distribution on the hull, ΔT = 0.250; (c2) Free-surface topology, ΔT = 0.250; (d1) Pressure distribution on the hull, ΔT = 0.375; (d2) Fre-surface topology, ΔT = 0.375; (e1) Pressure distribution on the hull, ΔT = 0.500; (e2) Free-surface topology, ΔT = 0.500.
Figure 6. Pressure distribution and free surface topology computed for the case μ = 90 deg. (a1) Pressure distribution on the hull, ΔT = 0; (a2) Free-surface topology, ΔT = 0; (b1) Pressure distribution on the hull, ΔT = 0.125; (b2) Free-surface topology, ΔT = 0.125; (c1) Pressure distribution on the hull, ΔT = 0.250; (c2) Free-surface topology, ΔT = 0.250; (d1) Pressure distribution on the hull, ΔT = 0.375; (d2) Fre-surface topology, ΔT = 0.375; (e1) Pressure distribution on the hull, ΔT = 0.500; (e2) Free-surface topology, ΔT = 0.500.
Jmse 08 00962 g006aJmse 08 00962 g006b
Figure 7. Heave response amplitude operator (RAO) (m/m) for v = 12 Kn, hw = 2 m. (a) μ = 0 deg; (b) μ = 45 deg; (c) μ = 90 deg; (d) μ = 135 deg; (e) μ = 180 deg.
Figure 7. Heave response amplitude operator (RAO) (m/m) for v = 12 Kn, hw = 2 m. (a) μ = 0 deg; (b) μ = 45 deg; (c) μ = 90 deg; (d) μ = 135 deg; (e) μ = 180 deg.
Jmse 08 00962 g007
Figure 8. Pitch RAO (rad/m), for v = 12 Kn, hw = 2 m. (a) μ = 0 deg; (b) μ = 45 deg; (c) μ = 90 deg; (d) μ = 135 deg; (e) μ = 180 deg.
Figure 8. Pitch RAO (rad/m), for v = 12 Kn, hw = 2 m. (a) μ = 0 deg; (b) μ = 45 deg; (c) μ = 90 deg; (d) μ = 135 deg; (e) μ = 180 deg.
Jmse 08 00962 g008
Figure 9. Heave RAO (m/m), DYN–SHIPFLOW–NUMECA comparison, v = 12 Kn, μ = 180 deg, hw = 2 m.
Figure 9. Heave RAO (m/m), DYN–SHIPFLOW–NUMECA comparison, v = 12 Kn, μ = 180 deg, hw = 2 m.
Jmse 08 00962 g009
Figure 10. Pitch RAO (rad/m), DYN–SHIPFLOW–NUMECA comparison, v = 12 Kn, μ = 180 deg, hw = 2 m.
Figure 10. Pitch RAO (rad/m), DYN–SHIPFLOW–NUMECA comparison, v = 12 Kn, μ = 180 deg, hw = 2 m.
Jmse 08 00962 g010
Figure 11. Roll RAO (rad/m), v = 12 Kn, KG = 10.8 m, hw = 2 m. (a) μ = 45 deg; (b) μ = 90 deg; (c) μ = 135 deg.
Figure 11. Roll RAO (rad/m), v = 12 Kn, KG = 10.8 m, hw = 2 m. (a) μ = 45 deg; (b) μ = 90 deg; (c) μ = 135 deg.
Jmse 08 00962 g011
Figure 12. Heave RAO (m/m), v = 12 Kn, hw = 8 m for BEM model. (a) μ = 0 deg; (b) μ = 45 deg; (c) μ = 90 deg; (d) μ = 135 deg; (e) μ = 180 deg.
Figure 12. Heave RAO (m/m), v = 12 Kn, hw = 8 m for BEM model. (a) μ = 0 deg; (b) μ = 45 deg; (c) μ = 90 deg; (d) μ = 135 deg; (e) μ = 180 deg.
Jmse 08 00962 g012
Figure 13. Pitch RAO (rad/m), v = 12 Kn, hw = 8 m for BEM model. (a) μ = 0 deg; (b) μ = 45 deg; (c) μ = 90 deg; (d) μ = 135 deg; (e) μ = 180 deg.
Figure 13. Pitch RAO (rad/m), v = 12 Kn, hw = 8 m for BEM model. (a) μ = 0 deg; (b) μ = 45 deg; (c) μ = 90 deg; (d) μ = 135 deg; (e) μ = 180 deg.
Jmse 08 00962 g013
Figure 14. Heave RAO (m/m), v = 24 Kn, hw = 2 m for BEM model. (a) μ = 0 deg; (b) μ = 45 deg; (c) μ = 90 deg; (d) μ = 135 deg; (e) μ = 180 deg.
Figure 14. Heave RAO (m/m), v = 24 Kn, hw = 2 m for BEM model. (a) μ = 0 deg; (b) μ = 45 deg; (c) μ = 90 deg; (d) μ = 135 deg; (e) μ = 180 deg.
Jmse 08 00962 g014
Figure 15. Pitch RAO (rad/m), v = 24 Kn, hw = 2 m. (a) μ = 0 deg; (b) μ = 45 deg; (c) μ = 90 deg; (d) μ = 135 deg; (e) μ = 180 deg.
Figure 15. Pitch RAO (rad/m), v = 24 Kn, hw = 2 m. (a) μ = 0 deg; (b) μ = 45 deg; (c) μ = 90 deg; (d) μ = 135 deg; (e) μ = 180 deg.
Jmse 08 00962 g015
Figure 16. Roll RAO (rad/m), v = 24 Kn, KG = 10.8 m, hw = 2 m. (a) μ = 45 deg; (b) μ = 90 deg; (c) μ = 135 deg.
Figure 16. Roll RAO (rad/m), v = 24 Kn, KG = 10.8 m, hw = 2 m. (a) μ = 45 deg; (b) μ = 90 deg; (c) μ = 135 deg.
Jmse 08 00962 g016
Figure 17. Heave RAO (m/m), v = 24 Kn, hw = 8 m for the BEM model. (a) μ = 0 deg; (b) μ = 45 deg; (c) μ = 90 deg; (d) μ = 135 deg; (e) μ = 180 deg.
Figure 17. Heave RAO (m/m), v = 24 Kn, hw = 8 m for the BEM model. (a) μ = 0 deg; (b) μ = 45 deg; (c) μ = 90 deg; (d) μ = 135 deg; (e) μ = 180 deg.
Jmse 08 00962 g017aJmse 08 00962 g017b
Figure 18. Pitch RAO (rad/m), v = 24 Kn, hw = 8 m for the BEM model. (a) μ = 0 deg; (b) μ = 45 deg; (c) μ = 90 deg; (d) μ = 135 deg; (e) μ = 180 deg.
Figure 18. Pitch RAO (rad/m), v = 24 Kn, hw = 8 m for the BEM model. (a) μ = 0 deg; (b) μ = 45 deg; (c) μ = 90 deg; (d) μ = 135 deg; (e) μ = 180 deg.
Jmse 08 00962 g018aJmse 08 00962 g018b
Figure 19. SHIPFLOW—Amplitude Spectrum by DFT method, v = 12 and 24 Kn, hw = 2 m, μ = 180 deg, ω = 0.4 rad/s (fe = 0.080 and 0.096 Hz). (a) Heave Amplitude Spectrum (m); (b) Pitch Amplitude Spectrum (rad).
Figure 19. SHIPFLOW—Amplitude Spectrum by DFT method, v = 12 and 24 Kn, hw = 2 m, μ = 180 deg, ω = 0.4 rad/s (fe = 0.080 and 0.096 Hz). (a) Heave Amplitude Spectrum (m); (b) Pitch Amplitude Spectrum (rad).
Jmse 08 00962 g019
Figure 20. SHIPFLOW— Amplitude Spectrum by DFT method, v = 12 and 24 Kn, hw = 8 m, μ = 180 deg, ω = 0.4 rad/s (fe = 0.080 and 0.096 Hz). (a) Heave Amplitude Spectrum (m); (b) Pitch Amplitude Spectrum (rad).
Figure 20. SHIPFLOW— Amplitude Spectrum by DFT method, v = 12 and 24 Kn, hw = 8 m, μ = 180 deg, ω = 0.4 rad/s (fe = 0.080 and 0.096 Hz). (a) Heave Amplitude Spectrum (m); (b) Pitch Amplitude Spectrum (rad).
Jmse 08 00962 g020
Figure 21. NUMECA–SHIPFLOW Time Records comparison v = 12 Kn, hw = 2 m, μ = 180 deg, ω = 0.4 rad/s (fe = 0.080). (a) Heave Time Record (m) (b) Pitch Time Record (rad).
Figure 21. NUMECA–SHIPFLOW Time Records comparison v = 12 Kn, hw = 2 m, μ = 180 deg, ω = 0.4 rad/s (fe = 0.080). (a) Heave Time Record (m) (b) Pitch Time Record (rad).
Jmse 08 00962 g021
Figure 22. Comparison of NUMECA–SHIPFLOW Amplitude Spectrum by DFT method, v = 12 Kn, hw = 2 m, μ = 180 deg, ω = 0.4 rad/s (fe = 0.080). (a) Heave Amplitude Spectrum [m]; (b) Pitch Amplitude Spectrum (rad).
Figure 22. Comparison of NUMECA–SHIPFLOW Amplitude Spectrum by DFT method, v = 12 Kn, hw = 2 m, μ = 180 deg, ω = 0.4 rad/s (fe = 0.080). (a) Heave Amplitude Spectrum [m]; (b) Pitch Amplitude Spectrum (rad).
Jmse 08 00962 g022
Table 1. KCS ship main dimensions.
Table 1. KCS ship main dimensions.
DimensionValueUnits
Length between the perpendiculars (LBP)230 m
Length of waterline (LWL)232.5 m
Beam at waterline (BWL)32.2m
Depth (D)19.0m
Design draft (T)10.8 m
Displacement (Δ)52,030m3
Block coefficient (CB)0.6505-
Ship wetted area (S)9422m2
Longitudinal center of buoyancy (LCB) (%LBP), fwd+−1.48
Longitudinal center of gravity (LCG) from the aft peak111.603 m
Vertical center of gravity (KG) from keel10.8 m
Moment of inertia (Kxx/B)0.35-
Moment of inertia (Kyy/LBP, Kzz/LBP)0.25-
Speeds (v)12.24Kn
Table 2. Test cases.
Table 2. Test cases.
hw
(m)
v
(Kn)
Methodµ (deg)ω
(rad/s)
212ST0, 45, 90, 135, 1800Step 0.012
BEM0.30.350.40.450.50.550.60.650.70.80.91
RANS1800.3 0.4-0.5-0.6-0.7---
24ST0, 45, 90, 135, 1800Step 0.012
BEM0.30.350.40.450.50.550.60.650.70.80.91
812ST0, 45, 90, 135, 1800Step 0.012
BEM0.30.350.40.450.50.550.60.650.70.80.91
24ST0, 45, 90, 135, 1800Step 0.012
BEM0.30.350.40.450.50.550.60.650.70.80.91
Table 3. Grid convergence test results.
Table 3. Grid convergence test results.
µ = 90 degµ = 180 deg
VariableHeaveRollHeavePitch
S1 (fine)1.1011880.0313960.7131470.013927
S2 (medium)1.1013660.0314280.7149340.013929
S3 (coarse)1.1042840.0316830.7197300.014070
rG2.12.12.12.1
ε210.0001790.0000320.0017880.000001
ε320.0029170.0002540.0047960.000141
RG0.0611950.1271750.3727260.010014
pG3.76542.7794631.3301846.109572
δG0.0000120.0000050.0010621.432 × 10−8
UG0.0000930.0000140.0010627.913 × 10−7
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pacuraru, F.; Domnisoru, L.; Pacuraru, S. On the Comparative Seakeeping Analysis of the Full Scale KCS by Several Hydrodynamic Approaches. J. Mar. Sci. Eng. 2020, 8, 962. https://doi.org/10.3390/jmse8120962

AMA Style

Pacuraru F, Domnisoru L, Pacuraru S. On the Comparative Seakeeping Analysis of the Full Scale KCS by Several Hydrodynamic Approaches. Journal of Marine Science and Engineering. 2020; 8(12):962. https://doi.org/10.3390/jmse8120962

Chicago/Turabian Style

Pacuraru, Florin, Leonard Domnisoru, and Sandita Pacuraru. 2020. "On the Comparative Seakeeping Analysis of the Full Scale KCS by Several Hydrodynamic Approaches" Journal of Marine Science and Engineering 8, no. 12: 962. https://doi.org/10.3390/jmse8120962

APA Style

Pacuraru, F., Domnisoru, L., & Pacuraru, S. (2020). On the Comparative Seakeeping Analysis of the Full Scale KCS by Several Hydrodynamic Approaches. Journal of Marine Science and Engineering, 8(12), 962. https://doi.org/10.3390/jmse8120962

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop