Next Article in Journal
Node Load and Location-Based Clustering Protocol for Underwater Acoustic Sensor Networks
Previous Article in Journal
A Self-Adaptive Compression Method for Ship Trajectories without Threshold Setting
Previous Article in Special Issue
Variable Natural Frequency Damper for Minimizing Response of Offshore Wind Turbine: Effect on Dynamic Response According to Inner Water Level
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Four-DOF Maneuvering Motion of a Container Ship in Shallow Water Based on CFD Approach

1
Faculty of Transportation Mechanical Engineering, The University of Danang—University of Science and Technology, Da Nang 50000, Vietnam
2
Faculty of Transportation Engineering, Ho Chi Minh City—University of Science and Technology—Vietnam National University (VNU-HCM), Ho Chi Minh City 700000, Vietnam
3
Changwon National University, Changwon 51140, Republic of Korea
*
Authors to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2024, 12(6), 981; https://doi.org/10.3390/jmse12060981
Submission received: 19 April 2024 / Revised: 30 May 2024 / Accepted: 7 June 2024 / Published: 11 June 2024
(This article belongs to the Special Issue Offshore Structures and Hydrodynamic Modeling)

Abstract

:
With the continuous increase in ship size combined with the generally slower increase in the sizes of waterways, the need for the prediction of ship maneuvering in shallow waterways continues to attract attention from the international scientific community. Ship behavior in shallow water is relevant in seabed effects that result in changing the hydrodynamic forces acting on a ship. In this study, the maneuvering characteristics of a container ship with four degrees of freedom in shallow water are analyzed. The Reynolds-Averaged Navier Stokes approach in Ansys Fluent code is used to produce the maneuvering coefficients through the simulations of forward running, static drift, static heel, circular motion, the combined motions, and the pure roll motion of the KRISO container ship. The maneuvering characteristics of the ship are estimated for evaluating the ship behaviors in shallow-water conditions. The obtained results show that the roll has a significant decrease and the ship’s turning diameter has a significant increase when the ship operates in a shallow waterway. The predicted maneuvering characteristics of the ship are in good agreement with those of free-running model tests, indicating that the numerical simulation based on the Computational Fluid Dynamics method has good capability to predict the maneuvering derivatives and the four-DOF ship maneuvering motion in shallow water as well.

Graphical Abstract

1. Introduction

A ship is designed for not only oceangoing transport but also moving in a restricted area. Maneuverability is one of the most important characteristics of the ship’s performance that allows the ship to operate in safe and efficient conditions. Ship maneuverability in the open sea focuses on the navigation of ship journeys in calm water or in a random sea [1,2,3,4]. However, ship behavior in shallow water is associated with interactions between the ship and the seabed that affect the dynamic behavior of the ship, changing its attitude and increasing risks [5]. Therefore, the International Maritime Organization (IMO) recommends that the maneuverability of a ship in shallow water be included in the maneuvering booklet (IMO Resolution A.601, 1987) [6]. Previously, most studies addressed ship resistance and three-degree-of-freedom (3-DOF) maneuvering behavior in shallow-water conditions [2,3,7,8,9,10]. One of the first examples of a study on ship maneuvering in shallow water can be found in the work of Gronarz [11], where the author performed an extended study of the hydrodynamics of the Panamax container ship at different water depths using a CFD-based simulation. The maneuvering derivatives of the ship in the shallow-water case were represented as a function of the derivatives of the ship in deep water by using both the numeral results and experimental data. The results of the 3-DOF maneuvering simulation showed that the decrease in water depth resulted in an increase in the turning radius and a decrease in the drift angle. Milanov and Chotukova carried out experiments on a KRISO Container Ship (KCS) in shallow water in the BSHC (Bulgarian Ship Hydrodynamics Centre) tank for three values of metacentric height [12]. It could be seen that the ship roll response to the rudder action in shallow water was much more moderate compared with the deep-water case. Toxopeus performed a numerical study on the motion of the KVLCC2 model to predict variations in the surge force, sway force, and yaw moment for different water depths [13]. The calculated results showed that the hydrodynamic forces acting on the hull in very shallow water have nonlinear forms of the sway velocity and yaw rate. Form factors also increase as the ship moves in shallow waters.
Delefortrie et.al. formulated the six-DOF maneuvering model of KVLCC2 at water depth ratios of 1.2, 1.3, and 1.8 [14]. The hydrodynamic force and moment acting on the hull were mainly obtained from the captive model tests, where the ship was forced to move with three horizontal degrees of freedom accompanying free heave and pitch. Some numerical results on the roll were introduced to the maneuvering model of the ship. It could be seen that the prediction of KVLCC2’s maneuvering characteristics were in good agreement with those of the free-running tests. However, the prediction of roll was underestimated in comparison with the roll obtained from the free-running model tests (FRMTs). Additionally, Ruiz et al. [15] studied the wave effects on the motion of the KVLCC2 model with six DOFs in shallow water using both numerical and experimental methods. The experiment on a four-DOF KCS with a rudder and propeller was conducted by Kim et al. [16]. The captive model tests with a Computerized Planar Motion Carriage were performed to estimate hydrodynamic coefficients, especially those related to roll. Hamid et al. investigated a four-DOF ONR tumblehome based on the CFD-generated maneuvering derivatives [17]. Then, maneuvering simulations were performed to evaluate the ship maneuvering behavior, measuring variables such as the surge velocity, sway velocity, yaw rate, roll angle, heading angle, and ship trajectory. The CFD method was confirmed to have sufficient ability to predict the ship’s maneuvering characteristics in calm water. The results also showed zigzag overshoot and the under-prediction of the ship’s turning diameter. Yasukawa and Hirata experimentally investigated a ship sailing in the heeled state [18]. Oblique towing and turning circle tests were performed on the container ship, pure car carrier, and ferry to measure the surge force, sway force, roll moment, and yaw moment. Hydrodynamic derivatives obtained from forces and moments were used to investigate the course stability of the ship. It was proven that the course stability was deeply related to nonlinear derivatives of the heel angle. Fukui [19] performed a circular motion test by setting a constant roll angle for a container ship and a passenger ferry in deep water. The effects of the roll angle on the hydrodynamic derivatives were extracted from the measured forces and moments. The suitable four-DOF mathematical model was proposed to simulate the turning trajectories and the overshoot of the zigzag movement at several metacentric height values. It was concluded that maneuvering variables were well predicted by the mathematical model. Gu et al. [20] also predicted hydrodynamic derivatives via the Reynolds-Averaged Navier Stokes (RANS) method to investigate the maneuvering of a four-DOF tumblehome ship in deep water. The static heel angle, static drift, and circular motion tests were performed to estimate the surge force, sway force, roll moment, and yaw moment. These forces and moments were approximated through Taylor series expansion to determine the maneuvering derivatives. The obtained results shows that the derivatives were in good agreement with the experiment values except for the cross-coupled derivatives Nvvr, Yϕ, and Kϕ. The obtained derivatives were introduced into the four-DOF MMG model to simulate the turning circle and zigzag tests. The outcome variables of the ship motion were well predicted in comparison with the FRMT results.
Although many studies have been conducted to investigate the effects of shallow water on ship maneuvering, little effort has been made when it comes to a four-DOF maneuvering model of a container ship in shallow-water conditions. Therefore, this study addresses the prediction of the maneuvering characteristics of a four-DOF KCS model based on maneuvering derivatives that are generated through a CFD-based simulation.

2. Methodology

2.1. Maneuvering of the Container Ship with Four Degrees of Freedom

Normally, the 3-DOF model is considered in the maneuvering study of a surface ship; however, the roll of a hull such as in a container ship needs to be considered. To describe the ship maneuvering motion, two coordinate systems are considered: an earth-fixed coordinate system (Ox0y0z0) and a body-fixed coordinate system (Oxyz). The coordinate systems follow the right-hand rule with the z-axes pointing downward and the x-axes pointing to the bow. The origin of the body-fixed coordinate system is set at midship. The body-fixed coordinate system is moving with the ship, and it is used to express the linear velocities and angular velocities of the vehicle. The earth-fixed reference is the inertial reference frame adopted to represent the position and orientation of the ship. According to SNAME (The Society of Naval Architects and Marine Engineers), the position, orientation, and velocities of the ship motion in 4 DOFs (Degrees Of Freedom) are described in the two coordinate systems as shown in Figure 1.
The Son and Nomoto model for maneuvering motion in the surge, sway, roll, and yaw [21] is chosen to be the baseline of a KCS maneuvering study. The general form of the model is shown in Equation (1). The hydrodynamic force model is shown in the following equations.
m u v r x G r 2 + y G r = X
m v + u r y G r 2 + p 2 + x G r = Y
I x p I x z r + I y z r 2 + I x y p r + m y G v p z G v + u r = K
I z r I y z r p I x y p 2 I z x p + m x G v + u r y G u v r = N
Regarding the ship’s motion in calm water, the external force and moment on the right sides of the equations consist of hydrodynamic (FH) and hydrostatic (FS) forces acting on the hull, propeller forces (FP), and rudder forces (FR). The hydrodynamic forces (X, Y) and moments (K, N) of the hull, rudder, and the propeller are modeled as follows:
X = X u u u 2 + 1 t T J + X v r v r + X v v v 2 + C v δ v δ + X r r r 2 + X ϕ ϕ ϕ 2 + X δ δ δ 2
Y = Y v v + Y r r + Y p p + Y ϕ ϕ + Y v v v v + Y r r r r 3 + Y δ δ + Y δ δ δ δ 3 + Y v v r v 2 + Y v r r v r 2 + Y v v ϕ v 2 ϕ + Y r r ϕ r 2 ϕ + Y ϕ ϕ r r ϕ 2
K = K v v + K r r + K p p + K ϕ ϕ + K v v v v 3 + K r r r r 3 + K v v r v 2 r + K r r v r 2 + K v v ϕ v 2 ϕ + K r r ϕ r 2 + K ϕ ϕ r r ϕ 2 ( 1 + a H ) z R Y R cos δ
N = N v v + N r r + N p p + N ϕ ϕ + N v v v v + N r r r r 3 + N δ δ + N δ δ δ δ 3 + N v v r v 2 r + N r r v r 2 v + N v v ϕ v 2 ϕ + N r r ϕ r 2 ϕ + N ϕ ϕ r r ϕ 2 + N v v δ v 2 δ + N δ δ v v δ 2
Here, (u′, v′) and (p′, r′) stand for the dimensionless linear velocities along x and y axes and rotational velocities about the x and z axes, respectively. The hydrodynamic coefficients are derived from dimensionless values of the hydrodynamic force and moment by (½ρL2U2) and (½ ρL3U2).

2.2. Numerical Method

The fluid flow around the ship will be modeled using the Reynolds-Averaged Navier Stokes (RANS) method. The governing equations [22,23] of the method for the incompressible flow over the ship can be expressed in the following equations:
u ¯ i x i = 0
u ¯ i t + u ¯ j u ¯ i x j = p ¯ x i + ν 2 u ¯ i x j x j + g τ i j x j
The Reynolds stress, τ i j = u i u j ¯ , is modeled with the function of eddy viscosity and kinetic energy based on the Boussinesq hypothesis. The turbulent model, k-ω SST (Shear Stress Transport) [22,24,25,26,27], for solving turbulent flow around the ship is presented as follows:
k t + u ¯ i k x i = P k β * ω k + x j ( v + σ k v t ) k x i
ω t + u ¯ i ω x i = γ μ t P k β ω 2 + x j v + σ ω v t ω x i + 2 ( 1 F 1 ) σ ω 2 ω k x i ω x i
μ t = ρ a 1 k m a x a 1 ω , Ω F 2
The variables are defined as follows:
P k = τ i j u i x j
τ i j = μ t 2 S i j 2 3 u k x k δ i j 2 3 ρ k δ i j
S i j = 1 2 u i x j + u j x i
β i = F 1 β i , 1 + 1 F 1 β i , 2
σ k = 1 F 1 / σ k , 1 + 1 F 1 σ k , 2
σ ω = 1 F 1 / σ ω , 1 + 1 F 1 σ ω , 2
F 1 = t a n h a r g 1 4
F 2 = t a n h a r g 2 2
a r g 1 = m i n m a x k β * ω y ; 500 ν y 2 ω , 4 ρ σ ω 2 k C D k ω y 2
a r g 2 = m a x 2 k 0.09 ω y ; 500 ν y 2 ω
C D k ω = m a x 2 ρ σ ω , 2 1 ω k x i ω x i , 10 20
Here, the model constants are as follows: σk,1 = 1.176, σω,1 = 2.0, σk,2 = 1.0, σω,2 = 1.168, a1 = 0.31, βi,1 = 0.075, and βi,2 = 0.0828.
The variations in free water surface [28] are calculated from the formulation of VOF (Volume of Fluid) in which the volume fraction (αq) will be 0 or 1 for the cell empty of water or full of water, correspondingly. The cell contains both water and air phases if the volume fraction is between 0 and 1. Continuity equation for the volume fraction of the phases, volume fraction of the primary phase, and implicit scheme for calculating the face fluxes for all cells are presented in the following equations.
1 ρ q t α q ρ q + α q ρ q v q = S α q + p = 1 n m ˙ p q m ˙ q p
q = 1 n α q = 1
α q n + 1 ρ q n + 1 α q n ρ q n Δ t V + f ρ q n + 1 U f n + 1 α q , f n + 1 = S α q + p = 1 n m ˙ p q m ˙ q p V
Here, indexes (n, n+1) stand for previous time step and current time step. αq,f is the face value of the qth volume fraction. V is the volume of cell and Uf is volume flux through the face. m ˙ p q   a n d   m ˙ q p are the mass transfers from phase q to phase p and vice versa. S α q is the source term.

3. Numerical Study

3.1. KCS Ship Model and Simulation Setup

The candidate ship presented in this study is KRISO Container Ship (KCS) with a bulbous bow and transom. The ship is equipped with a fixed-pitch five-blade propeller and a semi-balance horn rudder as shown in Figure 2. As mentioned in SIMMAN (2020) [29], the full-scaled KCS is assumed to be operated at a speed of 8.75 knots in shallow-water area without considering the bank effects. Thus, a corresponding scaled model with a length between perpendiculars of 1.5 m is used and a flat sea bottom is specified in this numerical study. The stationary straight-line motions, circular motions, the combined motions, and the pure roll are simulated for predicting the maneuvering derivatives of the 4-DOF model of KCS ship at water depth ratios of 1.5 and 2.0. The main particulars of the hull, propeller, and rudder are tabulated in Table 1. The constraints of the tests are shown in Table 2 for forward, drift, heel, circular motion, and combined motion tests, and Table 3 for harmonic roll test.
In this study, the hydrodynamic forces acting on the ship constitute the main concern for estimating the maneuvering derivatives. The results of the CFD-based simulation are verified with the experimental results from KRISO [30] and Gronaz [11]. Straight-line tests, circular tests, and combined tests are numerically simulated using the RANS-based method for incompressible flow. The fluid domains are chosen to avoid the influence of boundary conditions such as inlet, outlet, and side walls. For the stationary straight-line tests and the harmonic tests, the rectangular domain has opted. In the stationary straight-line tests, the domain is divided into two parts, the inner part and outer part, and the fluid flow over the ship is simulated in the stationary reference frame. On the other hand, the circular domain is used to simulate the steady circular tests in the rotating reference frame. According to International Towing Tank Conference (ITTC) Recommendation on Practical Guidelines for Ship CFD Application [31], the boundaries should be located 1.0~2.0 Lpp away from the hull. In addition, the domain size must be increased in cases of ship simulation with the free surface condition. Therefore, the dimensions of the rectangular domain are chosen to be 6.5 Lpp in length, 4.0 Lpp in width, and 1.0 Lpp in height. The inner radius, outer radius, and thickness of the circular domains are 0.2 Lpp, 5.0 Lpp, and 0.6 Lpp. In addition, physical conditions will be applied to the boundary of the domain. The front and back sides of the rectangular domain are assigned as corresponding pressure inlet and pressure outlet. The slip condition is set for the side walls and the no-slip condition is specified for hull surfaces and sea bottom. The top face of the domain is considered as representing the symmetry conditions. For the circular domain, the front face and outer surface are assigned to the velocity inlet condition. The pressure outlet condition is set for the back face. The ship hull and bottom face are considered as representing the no-slip conditions. Figure 3 and Figure 4 show the typical domains and their dimensions for the straight-line and oblique tests, 2D pure roll, and circular tests, respectively.

3.2. Mesh Independent Study

Grid generation is the following step for discretizing the fluid domain into a number of control volumes. The hybrid mesh is adopted for the stationary straight-line tests oblique towing tests. The unstructured mesh is more flexible to facilitate the mesh generation of complex geometry and it, therefore, is used to discretize the inner domain including the hull. Specifically, the boundary layer domain is discretized using a prism mesh grown out of triangle mesh on the hull surface for resolving the viscous flow surrounding the hull. The tetrahedral mesh is used elsewhere away from the hull. The height of the first prism layer next to the ship wall is calculated for satisfying the y+ value of 10, corresponding to the transition flow regime at the Reynolds number of 5.4 × 105. On the other hand, the outer part is discretized using the structured mesh that can help reduce the grid size. Hexahedral elements are applied to the outer domain. The typical meshes of the oblique motion and the circular motion are illustrated in Figure 3 and Figure 4, respectively. The numbers of elements in the corresponding meshes are approximate 5.6 million and 3.4 million. The maximum skewness of the meshes is kept smaller than 0.84. The RANS-based solver in ANSYS Fluent code is chosen as a tool for simulating the fluid flow around the hull. According to the ITTC Practical Guidelines for Ship CFD Application [31], the two-equation turbulence models have shown to be able to provide an accurate prediction of the ship hydrodynamics. Thus, the k-ω SST turbulence model is opted for modeling turbulent flow around the hull. SIMPLE (semi-implicit method for pressure link equations) algorithm is employed to solve the governing equations iteratively. The least-square cell-based method is used for evaluating the gradient of flow variables. The quantities at cell faces are calculated from the cell-centered values by using the second-order upwind method.
For the concerned geometry, a series of grids are generated and the CFD computations are performed. The variations in the results of the CFD computation are then analyzed. Grid Convergence Index (GCI) method [32] is applied to evaluate the variation and the mesh discretization error. According to this method, at least three levels of grids—coarse, medium, and fine are generated. The grid convergence parameters for the fine mesh (1), medium mesh (2), and coarse mesh (3) are defined as follows:
  • Grid convergence rate:
p = ln ε 32 ε 21 ln r G
  • Grid convergence ratio:
R G = ε 21 ε 32
  • Grid convergence index:
G C I i j = F S e i j r p 1
Here, εij = ϕjϕi, eij = (ϕjϕi)/ϕj, and ϕi denotes the solution on the ith grid. Fs is 1.5 for the study with unstructured mesh.
GCI (Grid Convergence Index) method is applied to investigate the discretization error for the three meshes: fine mesh (1), medium mesh (2), and coarse mesh (3). The cases of static drift test at β = 12° and combined drift-heel test at β = 12° and ϕ = 8° are chosen to examine the discretization errors. The generated mesh sizes for the static drift test are 1978999 nodes, 2638389 nodes, and 2794898 nodes, corresponding to the coarse grid, medium grid, and fine grid while these are 2433339 nodes, 2561101 nodes, and 2873365 nodes for the combined test case. The slices of the three meshes for these cases are shown in Figure 5. Table 4 shows the hydrodynamic forces and moments with respect to the three meshes. It is visible that the errors of the solutions are decreased when the meshes are generated more finely. According to Table 5, the numerical uncertainties of the X, Y, and N obtained from the static drift test for the coarse mesh, medium mesh, and fine mesh are 7.76%, 4.62%, and 2.22%, respectively. Similarly, the numerical uncertainties for the mesh cases of the combined drift-heel test are 6.56%, 1.72%, 2.18%, and 0.85% for X, Y, K, and N, respectively. These quantities indicate that the solutions of Y, K, and N of the fine mesh are independent of grid generation while the X force has a small change when the mesh quality is varied. From the grid-independent study, the fine mesh is chosen for all the CFD-based simulations in this study.

4. Results

4.1. Hydrodynamic Forces and Moments Acting on the KCS Model

Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12, Figure 13, Figure 14 and Figure 15 show the computation results of hydrodynamic forces and moments of the straight-line tests at water depth ratios of 1.5 and 2.0. In general, the overall trend shows that the hydrodynamic forces and moments are augmented when the water depth level is decreased. The roll moment linearly increases with the increase in the roll angle while the other forces and moments have nonlinear variation with respect to the motion constraints. The similar trends for the results of the steady circular tests and the combined heel-CMT and combined drift-heel tests are shown in Figure 8, Figure 10, Figure 11, Figure 12, Figure 13, Figure 14 and Figure 15. The relationship between the surge force, sway force, and yaw moment with respect to the yaw rate is nonlinear. These forces and moments increase with the increase in the yaw rate.
For the harmonic motions, pure roll tests of the midship section are carried out. Fourier’s analysis is applied to analyze the resulting roll moment to derive the in-phase and out-phase components. The roll moment of the ship is obtained by integrating the roll moment of the section over the length of the ship. Figure 16 shows the in-phase and out-phase components of the roll moment for the midship section of the KCS. The added roll moment and roll damping are estimated through the linear approximation of the in-phase component and out-phase component, respectively.

4.2. Verification Study

The maneuvering derivatives determined from the current CFD analysis of the KCS in shallow water are verified with the derivatives obtained from the experimental results on KRISO from SIMMAN (2020) [30], the experimental work by Gronarz [11]. From the comparison shown in Figure 17, Figure 18, Figure 19 and Figure 20, it is clear to note that the sway force and yaw moment obtained from the static drift tests are in good agreement with those experimental results from the abovementioned researchers. On the other hand, there is a large discrepancy in surge force between the CFD results and the experimental results on KRISO and some discrepancy in the sway-velocity dependency of the surge force between the CFD result and the experimental data in [11] can be found. Thus, the differences between the current CFD results and the KRISO experimental results might be stem from the interaction between the hull, propeller, and rudder. The sway velocity-dependent derivatives of sway forces and yaw moments are in good agreement with the experimental results on KRISO. However, the nonlinear derivatives of the yaw moment obtained from the KRISO experiment are greater than the derivatives determined from the current CFD computation. On the contrary, the rotary derivatives obtained from Gronaz’s experiments are rather small in comparison with the current CFD results. It is curious that the current CFD study yields a rotary derivative of sway force that is greater than the rotary derivative of the yaw moment but the opposite trend is found in the results of Gronaz’s study. The comparison of rotary derivatives and the cross-coupled derivatives is depicted in Figure 19.

4.3. Maneuvering Simulation

The obtained maneuvering derivatives of the surge force, sway force, roll moment, and yaw moment are tabulated in Table 6 and Table 7 for the water depth ratios of 1.5 and 2.0, respectively. It should be noticed that the surge, sway, and yaw inertia terms have been taken from Mucha [33]. The hydrodynamic forces and moments of the propeller and rudder are referred from SIMMAN (2020) [30]. The vertical center of gravity of the model has been taken from the research of Tezdogan [34]. For evaluating the KCS maneuvering in shallow water, standard maneuvering simulations are performed at a speed of 8.75 knots for the cases of water depth ratios of 1.5 and 2.0. The results from this maneuvering study are then compared with the results of free-running model tests involving the hard rudder turning maneuvers at δR = −35° [12].
The standard maneuvers are simulated for evaluating the ship maneuvering characteristics in shallow-water conditions. The equation of motion is established after introducing the hydrodynamic forces, the rudder forces, and the propeller forces. The equation is then numerically solved by using Microsoft Visual C++ 10 Express. The fourth-order Runge–Kutta method is employed to find the approximate values of the ship motion variables for a given rudder execute and propeller RPM. Figure 21 shows the simulation results of the turning circle maneuvers with a rudder deflection of −35 degrees by using the maneuvering derivatives in Table 7. The starboard turning trajectory, roll angle, and relationship between rudder and heading angle obtained from the FRMT [12] are also presented for the comparison. It is seen that the numerical turning trajectory has been adapted to the experimental result. The rate of turn is decreased, resulting in a larger turning circle in comparison with the turning circle in the deep-water case. In the case of the roll, the roll angle obtained from the current CFD is in good agreement with the estimated roll in the FRMT. It has been demonstrated that the ship advance and tactical diameter have been well predicted in the current CFD-based maneuvering study of a container ship. Figure 22 shows the simulation results of zig-zag 20/20 maneuver. The heading angle obtained from the current CFD-based maneuvering study is fairly good agreement with FRMT.

5. Conclusions

The need for the maneuverability of a ship in shallow water increases due to the increase in ship size in recent years that calls for the evaluation of maneuvering prediction methods. CFD-based computation is one of the efficient methods, and this can compromise between accuracy and high-cost performance in the preliminary design of a ship. Therefore, the dependence of the CFD-based results on the mesh density has been analyzed and the four-DOF maneuverability of the container ship in shallow-water conditions has been analyzed through the CFD-based estimation of the maneuvering derivatives. All the CFD-based analyses have been implemented for the incompressible flow around the bare hull by using the RANS solver in Ansys Fluent code as shown below.
-
The CFD RANS-based simulations of the static drift tests, static heel tests, steady circular tests, and the combined tests were carried out to estimate the hydrodynamic forces and moments acting on the KCS model. The roll damping and added roll moment of inertia were obtained from the virtual pure roll test.
-
The grid independence study was implemented for checking the convergence of the numerical results. The obtained results showed that the numerical uncertainties of the combined drift-heel test were 6.56%, 1.72%, 2.18%, and 0.85% for X, Y, K, and N, correspondingly, and the process of discretization in a fine mesh was chosen to generate the grids for simulating the virtual captive model tests in this study.
-
The maneuvering derivatives are derived from the results of the virtual tests by using the Taylor series. These were introduced to the maneuvering mathematical model of the ship with four DOFs for evaluating the maneuvering characteristics in shallow water and medium-shallow-water conditions. In general, the predicted derivatives were in good agreement with the experimental data. The KCS maneuvering characteristics obtained from the current CFD-based maneuvering study complied with the results of free-running model tests.
Although the maneuvering characteristics of the container ship using the hydrodynamic derivatives obtained from the CFD-based simulation were verified with the experimental results, the rotary derivatives of surge force, sway force, and yaw moments were not totally in agreement with the corresponding derivatives obtained from the experimental study. The difference may have been due to estimation errors, so uncertainty analysis in the CFD validation methodology should be considered in further studies.

Author Contributions

Conceptualization, H.K.Y.; methodology, T.L.P. and T.T.N.; software, T.T.N. and T.L.P.; validation, T.T.N. and H.K.Y.; formal analysis, T.L.M.; investigation, T.L.M. and T.-H.L.; resources, T.T.N. and T.-H.L.; data curation, T.T.N.; writing—original draft preparation, T.T.N. and T.L.P.; writing—review and editing, T.T.N. and T.L.P.; supervision, H.K.Y.; project administration, T.T.N.; funding acquisition, T.T.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research is funded by the Funds for Science and Technology Development of the University of Danang under project number B2020-DN02-85.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in this article (Tables and Figures).

Acknowledgments

We acknowledge the support of time and facilities from Ho Chi Minh City University of Technology (HCMUT) and Vietnam National University Ho Chi Minh City (VNU-HCM) for this study.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Chen, J.; Zou, Z.J.; Xi, C.; Li, X.; Lu, Z. CFD-based simulation of the flow around a ship in turning motion at low speed. J. Mar. Sci. Technol. 2017, 22, 784–796. [Google Scholar] [CrossRef]
  2. He, S.; Kellett, P.; Yuan, Z.; Incecik, A.; Turan, O.; Boulougouris, E. Manoeuvring prediction based on CFD generated derivatives. J. Hydrodyn. 2016, 28, 284–292. [Google Scholar] [CrossRef]
  3. Jachowski, J. Assessment of Ship Squat in Shallow Water using CFD. Arch. Civ. Mech. Eng. 2008, 8, 27–35. [Google Scholar] [CrossRef]
  4. Sung, Y.J.; Park, S.H. Prediction of ship manoeuvring performance based on virtual captive model tests. J. Soc. Nav. Archit. Korea 2015, 52, 407–417. [Google Scholar] [CrossRef]
  5. Vantorre, M. Review of practical methods for assessing shallow and restricted water effects. In Proceedings of the International Conference on Marine Simulation and Ship Maneuverability, Kanazawa, Japan, 25–28 August 2003. [Google Scholar]
  6. IMO Standard for Ship Manoeuvrability: ANNEX 6. Available online: https://wwwcdn.imo.org/localresources/en/KnowledgeCentre/IndexofIMOResolutions/MSCResolutions/MSC.137(76).pdf (accessed on 29 May 2024).
  7. Patel, P.K.; Premchand, M. Numerical Investigation of the Influence of Water Depth on Ship Resistance. J. Comput. Appl. 2015, 116, 11–17. [Google Scholar]
  8. Terziev, M.; Tezdogan, T.; Oguz, E.; Gourlay, T.; Demirel, Y.K.; Incecik, A. Numerical investigation of the behaviour and performance of ships advancing through restricted shallow waters. J. Fluids Struct. 2018, 76, 185–215. [Google Scholar] [CrossRef]
  9. Eloot, K.; Vantorre, M. Ship behaviour in shallow water and confined water: An overview of hydrodynamic effects through EFD. In Proceedings of the 2nd International Conference on Ship Manoeuvring in Shallow and Confined Water, Trondheim, Norway, 18–20 May 2011. [Google Scholar]
  10. Liu, Y.; Zou, Z.J.; Zou, L. RANS-based numerical simulation of captive model test in shallow water for the DTC container carrier. In Proceedings of the Fourth International Conference on Ship Manoeuvring in Shallow and Confined Water, Hamburg, Germany, 23–25 May 2016. [Google Scholar]
  11. Gronarz, A. Rechnerische Simulation der Schiffsbewegung beim Manövrieren unter Besonderer Berücksichtigung der Abhängigkeit von der Wassertiefe. Ph.D. Thesis, Unversity of Duisburg-Essen, Duisburg/Essen, Germany, 11 September 1997. [Google Scholar]
  12. Milanov, E.; Chotukova, V.Y. Roll motion of container ship in shallow water. In Proceedings of the International Conference on Ship Manoeuvring in Shallow and Confined Water, Antwerp, Belgium, 13–15 May 2009. [Google Scholar]
  13. Toxopeus, S.L. Viscous-flow Calculations for KVLCC2 in Deep and Shallow Water. In Proceedings of the International Conference on Computational Method in Marine Engineering, Lisbon, Portugal, 28–30 September 2011. [Google Scholar]
  14. Delefortrie, G.; Eloot, K.; Lataire, E.; VanHoydonck, W.; Vantorre, M. Captive model tests bases 6DOF shallow water model. In Proceedings of the 4th MASHCON Conference, Hamburg, Germany, 23–25 May 2016. [Google Scholar]
  15. Ruiz, M.T.; Caluwé, S.D.; Zwijnsvooorde, T.V.; Delefortrie, G.; Vantorre, M. Wave effects in 6DOF on a ship in shallow water. In Proceedings of the MARSIM Conference, Newcastle upon Tyne, UK, 8–11 September 2015. [Google Scholar]
  16. Kim, Y.G.; Yeo, D.; Namsun, S.; Kim, S.; Yun, K.; Oh, B. Prediction maneuverability of KCS with 4 degrees of freedom. J. Soc. Nav. Archit. Korea 2011, 48, 267–274. [Google Scholar] [CrossRef]
  17. Hamid, S.H.; Araki, M.; Umeda, N.; Sano, M.; Yeo, D.J.; Toda, Y.; Carrica, P.M.; Stern, F. CFD, system–based, and EFD preliminary investigation of ONR tumblehome Instability and capsize with evaluation of the mathematical model. In Proceedings of the 12th International Ship Stability Workshop, Washington, DC, USA, 12–15 June 2011. [Google Scholar]
  18. Yasukawa, H.; Hirata, N. Maneuverability and hydrodynamic derivatives of ships traveling in heeled condition. J. Jpn. Soc. Nav. Archit. Ocean Eng. 2013, 17, 19–29. [Google Scholar]
  19. Fukui, Y.; Yokota, H.; Yano, H.; Kondo, M.; Nakano, T.; Yoshimura, Y. 4-DOF Mathematical model for manoeuvring simulation including roll motion. J. Soc. Nav. Archit. Korea 2016, 24, 167–179. [Google Scholar]
  20. Gu, H.; Zou, Z. System-based investigation on 4-DOF ship maneuvering with hydrodynamic derivatives determined by RANS simulation of captive model tests. Appl. Ocean Res 2018, 68, 11–25. [Google Scholar] [CrossRef]
  21. Fossen, T.I.B. Guidance and Control of Ocean Vehicles. Chichester; Wiley: New York, NY, USA, 1994. [Google Scholar]
  22. Piquet, J. Turbulent Flows: Models and Physics; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2013. [Google Scholar]
  23. Bakker, A. Applied Computational Fluid Dynamics. Available online: http://www.bakker.org/dartmouth06/engs150/11-bl.pdf (accessed on 29 May 2024).
  24. Menter, F.R. Two-Equation Eddy-Viscosity Turbulence Models for Engineering Applications. AIAA J. 1994, 32, 1598–1605. [Google Scholar] [CrossRef]
  25. Shih, T.H.; Liou, W.W.; Shabbir, A.; Yang, Y.; Zhu, Y. A New k – Eddy Viscosity Model for High Reynolds Number Turbulent Flows - Model Development and Validation. Comput. Fluids 1995, 24, 227–238. [Google Scholar] [CrossRef]
  26. Wilcox, D.C. Re-assessment of the scale-determining equation for advanced turbulence models. AIAA J. 1988, 26, 1299–1310. [Google Scholar] [CrossRef]
  27. CFD-Online. Turbulence Intensity. CFD-Online. Available online: https://www.cfd-online.com/Wiki/Turbulence_intensity (accessed on 29 May 2024).
  28. PIANC. Capability of Ship Manoeuvering Simulation Models for Approach Channels and Fairways in Harbours; PIANC Report No 20; PIANC: Brussels, Belgium, 1992. [Google Scholar]
  29. ANSYS Fluent Theory Guide. Available online: https://www.afs.enea.it/project/neptunius/docs/fluent/html/th/node297.htm (accessed on 29 May 2024).
  30. SIMMAN. Workshop on Verification and Validation of Ship Manoeuvring Simulation Methods, Copenhagen, Denmark. Available online: https://www.simman2020.kr (accessed on 29 May 2024).
  31. International Towing Tank Conference (ITTC). ITTC-Recommended Procedures and Guidelines: Practical Guidelines for Ship CFD Applications; International Towing Tank Conference (ITTC): St. John S, NL, Canada, 2011. [Google Scholar]
  32. Blancas, F.E.; Roache, P.J. An Evaluation of the GCI for Unstructured Grids. Available online: https://cstools.asme.org (accessed on 29 May 2024).
  33. Mucha, P. On simulation-based ship maneuvering prediction in deep and shallow water. Ph.D. Thesis, University of Duisburg-Essen, Duisburg/Essen, Germany, 9 January 2017. [Google Scholar]
  34. 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]
Figure 1. Coordinate systems and symbols.
Figure 1. Coordinate systems and symbols.
Jmse 12 00981 g001
Figure 2. KCS model and appendages: (Black line) Ship hull; (Green line) Rudder; (Red line) Propeller; (Blue line) draft level; (magenta line) base line.
Figure 2. KCS model and appendages: (Black line) Ship hull; (Green line) Rudder; (Red line) Propeller; (Blue line) draft level; (magenta line) base line.
Jmse 12 00981 g002
Figure 3. The fluid domain for straight-line motion and oblique motion: (Orange line) Ship hull; (Green line) Refinement region around the ship; (Blue line) Boundary domain.
Figure 3. The fluid domain for straight-line motion and oblique motion: (Orange line) Ship hull; (Green line) Refinement region around the ship; (Blue line) Boundary domain.
Jmse 12 00981 g003
Figure 4. Fluid domain and boundary conditions for (a) 2D pure roll and (b) circular motion.
Figure 4. Fluid domain and boundary conditions for (a) 2D pure roll and (b) circular motion.
Jmse 12 00981 g004
Figure 5. Mesh sizes for static drift test (left) and combined drift-heel test (right), h/T = 2.0.
Figure 5. Mesh sizes for static drift test (left) and combined drift-heel test (right), h/T = 2.0.
Jmse 12 00981 g005
Figure 6. Forward running test.
Figure 6. Forward running test.
Jmse 12 00981 g006
Figure 7. Static heel test.
Figure 7. Static heel test.
Jmse 12 00981 g007
Figure 8. Circular motion test.
Figure 8. Circular motion test.
Jmse 12 00981 g008
Figure 9. Static drift test.
Figure 9. Static drift test.
Jmse 12 00981 g009
Figure 10. Combined heel-drift test, h/T = 1.5.
Figure 10. Combined heel-drift test, h/T = 1.5.
Jmse 12 00981 g010
Figure 11. Combined heel-drift test, h/T = 2.0.
Figure 11. Combined heel-drift test, h/T = 2.0.
Jmse 12 00981 g011
Figure 12. Combined drift-CMT test, h/T = 1.5.
Figure 12. Combined drift-CMT test, h/T = 1.5.
Jmse 12 00981 g012
Figure 13. Combined drift-CMT test, h/T = 2.0.
Figure 13. Combined drift-CMT test, h/T = 2.0.
Jmse 12 00981 g013
Figure 14. Combined heel-CMT test, h/T = 1.5.
Figure 14. Combined heel-CMT test, h/T = 1.5.
Jmse 12 00981 g014
Figure 15. Combined heel-CMT test, h/T = 2.0.
Figure 15. Combined heel-CMT test, h/T = 2.0.
Jmse 12 00981 g015
Figure 16. Results of 2D pure roll test, h/T = 2.0.
Figure 16. Results of 2D pure roll test, h/T = 2.0.
Jmse 12 00981 g016
Figure 17. Comparison of velocity-dependent derivatives of surge force, sway force, and yaw moment, h/T = 1.5.
Figure 17. Comparison of velocity-dependent derivatives of surge force, sway force, and yaw moment, h/T = 1.5.
Jmse 12 00981 g017
Figure 18. Comparison of velocity-dependent derivatives of surge force, sway force, and yaw moment, h/T = 2.0.
Figure 18. Comparison of velocity-dependent derivatives of surge force, sway force, and yaw moment, h/T = 2.0.
Jmse 12 00981 g018
Figure 19. Comparison of velocity-dependent derivatives of surge force, sway force, and yaw moment, h/T = 1.5.
Figure 19. Comparison of velocity-dependent derivatives of surge force, sway force, and yaw moment, h/T = 1.5.
Jmse 12 00981 g019
Figure 20. Comparison of velocity-dependent derivatives of surge force, sway force, and yaw moment, h/T = 2.0.
Figure 20. Comparison of velocity-dependent derivatives of surge force, sway force, and yaw moment, h/T = 2.0.
Jmse 12 00981 g020
Figure 21. Turning circle of KCS model circle towards the starboard in comparison with the results of FRMT [22], δR = −35°, h/T = 2.0.
Figure 21. Turning circle of KCS model circle towards the starboard in comparison with the results of FRMT [22], δR = −35°, h/T = 2.0.
Jmse 12 00981 g021
Figure 22. Results of 20/20 zigzag test in comparison with the results of FRMT [22], h/T = 2.0.
Figure 22. Results of 20/20 zigzag test in comparison with the results of FRMT [22], h/T = 2.0.
Jmse 12 00981 g022aJmse 12 00981 g022b
Table 1. Main particulars of KCS and appendages.
Table 1. Main particulars of KCS and appendages.
DesignationUnitFull ScaleModel Scale
Lppm230.01.500
Bwlm32.200.210
Hm19.00.124
Tm10.80.070
CB-0.6510.651
m352,0300.0144
LCB%, fwd+−1.48−1.48
LCGm111.60.728
KGm14.320.0934
kxx/B-0.4000.400
kzz/Lpp-0.2500.250
GMTm0.6000.0039
Table 2. Test conditions for forward, drift, heel, and circular motion and combined motion.
Table 2. Test conditions for forward, drift, heel, and circular motion and combined motion.
CaseV (Knots)β (Degree)ϕ (Degree)r′ [-]
Straight forward7, 8.75, 10, 12000
Static drift8.754, 8, 12, 1600
Static heel8.7500, 2, 4, 6, 8, 100
Circular motion8.75000.3, 0.35, 0.40, 0.50
Combined drift-heel8.754, 6, 8, 122, 4, 80
Combined heel-CMT8.7502, 4, 60.3, 0.4, 0.5
Combined drift-CMT8.754, 6, 8, 1200.3, 0.4, 0.5
Table 3. Test conditions of the harmonic roll.
Table 3. Test conditions of the harmonic roll.
ϕmax [Degree]ω (rad/s)p (rad/s) p ˙ (rad/s2)
50.65.236 × 10−2−3.141 × 10−2
50.86.981 × 10−2−5.585 × 10−2
51.08.726 × 10−2−8.726 × 10−2
Table 4. Hydrodynamic forces and moments in static drift test and drift-heel test, h/T = 2.0.
Table 4. Hydrodynamic forces and moments in static drift test and drift-heel test, h/T = 2.0.
Static Drift Test, β = 12°Combined Drift-Heel Test, β = 12°, ϕ = 8°
Mesh sizeX [N]Y [N]N [N.m]X [N]Y [N]K [N.m]N [N.m]
Coarse−0.199−0.760−0.320−1.39 × 10−3−3.31 × 10−31.71 × 10−31.83 × 10−3
Medium−0.202−0.778−0.327−1.44 × 10−3−3.46 × 10−31.75 × 10−31.85 × 10−3
Fine−0.204−0.789−0.330−1.47 × 10−3−3.5 × 10−31.77 × 10−31.85 × 10−3
Table 5. Grid convergence study.
Table 5. Grid convergence study.
Static Drift Test, β = 12°Combined Drift-Heel Test, β = 12°, ϕ = 8°
XYNXYKN
ε32−2.02 × 10−5−1.26 × 10−4−2.94 × 10−5−5.03 × 10−5−1.49 × 10−43.79 × 10−51.93 × 10−5
ε21−1.50 × 10−5−7.20 × 10−5−1.39 × 10−5−3.14 × 10−5−7.57 × 10−51.71 × 10−58.01 × 10−6
ε321.46 × 10−22.35 × 10−21.96 × 10−23.50 × 10−24.30 × 10−22.17 × 10−21.05 × 10−2
ε211.07 × 10−21.33 × 10−29.19 × 10−32.14 × 10−32.14 × 10−29.71 × 10−34.32 × 10−3
pG8.45 × 10−11.612.171.351.952.292.54
RG0.7460.5730.4720.6250.5090.4520.415
GCIfine7.76%4.62%2.22%6.56%1.72%2.18%0.85%
Table 6. Non-dimensional maneuvering derivatives of KCS model (×103), h/T = 1.5.
Table 6. Non-dimensional maneuvering derivatives of KCS model (×103), h/T = 1.5.
X-CoefficientsY-CoefficientK-CoefficientN-Coefficient
X’uu = −1.206Y′v = −14.425K′ϕ = −0.251N′v = −8.829
X’vv = −9.909Y′v|v| = −137.059K′v = −0.221N′v|v| = −6.325
X’ϕϕ = −4.348Y′r = 1.539K′vvv = 21.571N′r = −1.783
X’rr = −4.022Y′rrr = 25.05K′r = −0.0651N′rrr = −4.805
X’vr = 51.498Y′vrr = −90.860K′rrr = 0.4687N′vrr = 5.464
X’δδ = −2.223Y′vvr = 15.979 N′vvr = −66.939
X’ = −8.970Y′ϕ = −1.214 N′ϕ = −0.561
Y′ϕvv = 238.265 N′vvϕ = −108.258
Y′v|ϕ| = 107.489 N′v|ϕ| = 13.926
Y′ϕrr = −298.522 N′ϕrr = −5.052
Y′r|ϕ| = 157.118 N′r|ϕ| = 0.742
Y′δ = 5.652 N′δ = −2.304
Y′δδδ = −5.522 N′δδδ = 11.496
Y′vvδ = 44.616 N′vvδ = 1.496
Y′vδδ = 1759.4 N′vδδ = 787.52
Table 7. Non-dimensional maneuvering derivatives of KCS model (×103), h/T = 2.0.
Table 7. Non-dimensional maneuvering derivatives of KCS model (×103), h/T = 2.0.
X-CoefficientsY-CoefficientK-CoefficientN-Coefficient
X u ˙ = −0.417 Y v ˙ = −8.573 K p ˙ = −0.0193 N r ˙ = −1.242
X’uu = −1.107 Y r ˙ = 0.0414K′p = −0.0635 N v ˙ = 0.217
X’vv = −8.246Y′v = −8.685K′ϕ = −0.200N′v = −7.433
X’ϕϕ = −0.391Y′v|v| = −79.797K′v = 0.711N′v|v| = 1.334
X’rr = −2.257Y′r = 1.830K′vvv = 2.485N′r = −2.385
X’vr = 27.417Y′rrr = 7.673K′r = −0.071N′rrr = −1.026
X’δδ = −2.668Y′vrr = −64.647K′rrr = 0.312N′vrr = 0.372
X’ = −3.747Y′vvr = −28.413 N′vvr = −44.538
Y′ϕ = −1.180 N′ϕ = −0.619
Y′ϕvv = 87.101 N′vvϕ = −156.10
Y′v|ϕ| = 39.321 N′v|ϕ| = 36.530
Y′ϕrr = −125.100 N′ϕrr = −10.481
Y′r|ϕ| = 60.300 N′r|ϕ| = 4.318
Y′δ = −5.459 N′δ = 2.428
Y′δδδ = 5.258 N′δδδ = −1.916
Y′vvδ = −64.981 N′vvδ = 38.004
Y′vδδ = 762.950 N′vδδ = 322.200
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Nguyen, T.T.; Phan, T.L.; Le, T.-H.; Mai, T.L.; Yoon, H.K. Four-DOF Maneuvering Motion of a Container Ship in Shallow Water Based on CFD Approach. J. Mar. Sci. Eng. 2024, 12, 981. https://doi.org/10.3390/jmse12060981

AMA Style

Nguyen TT, Phan TL, Le T-H, Mai TL, Yoon HK. Four-DOF Maneuvering Motion of a Container Ship in Shallow Water Based on CFD Approach. Journal of Marine Science and Engineering. 2024; 12(6):981. https://doi.org/10.3390/jmse12060981

Chicago/Turabian Style

Nguyen, Tien Thua, Thanh Long Phan, Tat-Hien Le, Thi Loan Mai, and Hyeon Kyu Yoon. 2024. "Four-DOF Maneuvering Motion of a Container Ship in Shallow Water Based on CFD Approach" Journal of Marine Science and Engineering 12, no. 6: 981. https://doi.org/10.3390/jmse12060981

APA Style

Nguyen, T. T., Phan, T. L., Le, T. -H., Mai, T. L., & Yoon, H. K. (2024). Four-DOF Maneuvering Motion of a Container Ship in Shallow Water Based on CFD Approach. Journal of Marine Science and Engineering, 12(6), 981. https://doi.org/10.3390/jmse12060981

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