Next Article in Journal
Influence of Farming Intensity and Climate on Lowland Stream Nitrogen
Previous Article in Journal
Investigating Atrazine Concentrations in the Zwischenscholle Aquifer Using MODFLOW with the HYDRUS-1D Package and MT3DMS
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Seakeeping Analysis of Planing Craft under Large Wave Height

Science and Technology on Underwater Vehicle Laboratory, Harbin Engineering University, Harbin 150001, China
*
Author to whom correspondence should be addressed.
Water 2020, 12(4), 1020; https://doi.org/10.3390/w12041020
Submission received: 14 February 2020 / Revised: 26 March 2020 / Accepted: 31 March 2020 / Published: 2 April 2020
(This article belongs to the Section Hydraulics and Hydrodynamics)

Abstract

:
The purpose of this paper is to conduct a seakeeping analysis of planing craft under regular wave with large wave height. To obtain a reliable numerical method to simulate the sailing of planing craft in waves, Reynolds-averaged Navier–Stokes (RANS) solver and overset method are adopted. The motion response and resistance of the planing craft USV01 in regular wave were numerical predicted and compared with the corresponding seakeeping experimental tests. The results show that the numerical method has high accuracy. For further study, a new planing craft whose name is improved vessel is selected for simulation, the low steaming of the USV01 and improved vessel in regular wave with large wave height was simulated, and the seakeeping of the two vessels was studied. The analysis about the influence of wave length on the motion response and navigation configurations of the improved vessel under regular wave was carried out. Meanwhile, the influence of speed on different navigation configurations of the improved vessel was also analyzed. The improved vessel can provide better seakeeping, and a reduction in the speed of the vessel will benefit its seakeeping, irrespective of its navigation configuration.

1. Introduction

Planing craft have been widely used in various fields. The characteristic of the planing craft is that the bottom of the vessel will come into contact with the free surface when sailing at high speed. Since the bottom is relatively flat, the water pressure at the bottom is increased because the bottom of the boat squeezes water forward [1].
For a better hydrodynamic performance, numerical prediction is becoming increasingly important in the design of planing craft. The methods for predicting the resistance of planing craft include model test data, regression formulas, and semiempirical theory [2]. With the development of the computational fluid dynamics (CFD) method, the accuracy of numerical prediction in predicting the hydrodynamic performance of planing craft under calm water and wave conditions has been improved. In 2017, Diez et al. [3] selected the DTMB 5415 model (US Navy Combatant) as the parent hull of the destroyer, and the hull of the destroyer was improved based on the CFD method. In 2018, Campana et al. [4] used the numerical method to optimize the hull of a high-speed catamaran, and a real ocean environment was also considered. Their results indicated that larger computational power leads to an increase in computational efficiency and that the numerical method is becoming a necessary key to optimizing hull design.
In 2001, Azcueta et al. [5] conducted deep research into the free motion simulation of planing craft, based on the Comet, and the flow field around the hull were studied based on the k-ε turbulence model. Their results indicated that the numerical method can be adopted as an important complement to the study of the free surface flows around planing craft. In 2001, Caponnetto [6] conducted a hydrodynamic analysis of a planing hull with an unknown center of gravity with Reynolds-averaged Navier–Stokes (RANS) solver. They also compared their results with Savitsky’s method. In 2015, Weymouth et al. [7] conducted an analysis on the sailing of a ship under head wave conditions based on the unsteady RANS solver, and their work present accurate numerical results is uncertainties less than 2%. Their results indicate that the RANS solver has a better accuracy in simulating the sailing of a vessel and that the solver is suitable for both high speed and head wave conditions.
In 2015, Sun et al. [8] researched the grid factor in the numerical calculation of planing craft resistance, based on the prismatic glider. The effect of different grid parameters on the accuracy of the numerical results was analyzed. With this grid scheme, the resistance of planing craft in various navigation states is calculated, and the numerical results present good accuracy. Their results are useful for the grid setup and boundary condition definition in this paper. In 2013, Su et al. [9] put forward a numerical method to simulate the freedom motion of planing craft under head wave conditions, based on the six degrees of freedom (6-DOF) motion equation and volume of fluid (VOF) solver. The calculation results present a good agreement with the attitude of vessels under wave conditions. Their results show that the 6-DOF solver is suitable for simulating the motion of planing hull under head wave conditions. In 2017, De Marco et al. [10] conducted a seakeeping analysis of a stepped planing craft. A hydrodynamic experiment of the single-step hull model, as a new systematic series, was carried out, and the corresponding tests were based on numerical simulation, with an overset and morphing grid. Their high-accuracy numerical results show that the overset mesh can improve the efficiency, while ensuring the accuracy of simulation, and the setting of the overset grid is of great referential significance to our work. In 2017, Dashtimanesh et al. [11] provided a numerical prediction of the performance of two-stepped planing craft. based on the mathematical model and the obtained results present the good accuracy compared with the test results. Their results indicate that the CFD method has become a fundamental support for the hydrodynamics study and design of planing craft, and it can be highly beneficial with respect to the cost and duration of hydrodynamic tests.
With the increasing demand for planing craft, the ability to sail in rough sea condition is necessary, and this requires an adequate seakeeping performance. However, traditional design of planing craft makes it difficult to navigate safely. Due to the maturity of the numerical method, the complex motion problems associated with viscous flow, the nonlinear problem, and the transient response problem can all be directly solved by the numerical method. In 2000, Ikeda et al. [12] simulated the motion of planing craft under waves, and the porpoising oscillations and motion responses were numerically forecasted. This means that the unsteady motion of vessels can also be numerically simulated. This has great significance for the study of the extremely high-speed navigation of planing craft. In 2003, Azcueta et al. [13] conducted an analysis on the speed performance of a power boat, based on the free-surface RANS, and their method is suitable for predicting the severe responses of the vessel in waves, making it possible to predict the nonlinear motion of planing craft in rough sea conditions. In 2014, Begovic et al. [14,15] conducted research on the responses of hull forms of a vessel in waves, and the response amplitudes in both regular and irregular waves were obtained and compared with the corresponding test data. The Weibull distribution was adopted to analyze the motion response of the model in irregular waves, especially in relation to the acceleration of the center of gravity (CG). Irregular waves are closer to actual sea conditions, and the study of sailing in irregular waves is more reflective of actual sailing conditions. In 2016, Jiang et al. [16,17] conducted a study on the resistance performance of a planing trimaran with numerical and experimental methods. The CFD simulations were based on 2-DOF motion equations, and high speed sailing of the vessel was simulated. The results were validated with a good accuracy by comparing them with the test data. Their experimental design and setting is referential for our work, and their results indicate that numerical method can also be used to analyze the aerodynamic and hydrodynamic performance of trimaran.
In this paper, the STAR-CCM+ software (Germany) was adopted to simulate the sailing of a small-scale planing craft under regular wave. In regular wave with large wave height, due to the large wave height and wave disturbance, planing craft need to maintain low-speed navigation. This paper will focus on the research of the motion response and navigation configuration of planing craft in regular wave with large wave height at a low speed. The RANS solver and overset grid were adopted in the numerical simulation of the navigation of the vessel under head waves to obtain better accuracy, and validation was carried out by comparing the numerical results with the corresponding test results. Additionally, a new planing craft with better seakeeping is selected for further analysis, the low steaming of both vessels in regular wave with large wave height was numerically simulated and compared; and the influence of wave length and speed on the motion response and navigation configuration of the vessels is also analyzed.
Previous studies are mainly about the prediction of the motion response of planing craft, this paper focus on the seakeeping and changes of navigation configuration of planing under the regular wave with large wave height. In 2019, we have carried out a research on the seakeeping analysis of planing craft [18], mainly focus on the optimization effect of hydrofoil. This paper aims to illustrate a study on the navigation configuration of small planing craft under large wave height. And the numerical methods used in the two papers are similar.

2. The Numerical Method

2.1. RANS Equation

To simulate the viscous flow field around a sailing ship is to solve the Navier–Stokes equation. In response to this problem, the RANS Equation is commonly used in engineering. In the actual solving process, the time average value is used to replace the statistical average value, the momentum equations is as follows:
( ρ u i ) t + x i ( ρ u i u j ) = p x i + x j ( μ u i x j ρ u i u j ¯ ) + S i
Here u i and u j are the time mean of the velocity component, (i, j = 1, 2, 3), p is piezometric pressure coefficient, ρ is fluid density, μ is coefficient of dynamic viscosity, ρ u i u j ¯ is the Reynolds stress term, S i represents the generalized source term
Equation (1) and the continuity equation of incompressible fluid motion constitute the control equations for solving the viscous flow field around the hull. The continuity equation is as follows:
ρ t + ( ρ u i ) x i = 0
In the RANS equation, the Reynolds stress term is introduced in addition to the progressive equalization. In order to close the control equations, an appropriate turbulence model is introduced to calculate the Reynolds stress term.

2.2. Turbulence Model

In this paper, the Shear Stress Transport turbulence (SST) model is used to seal the control equations. The SST k–ω model [19] has been widely adopted in the solving of turbulence problems, especially the high Reynolds number flow problems, which is particularly suitable for the simulation of the high-speed navigation of planing craft. Similar to the k–ω model, the transport equation of k (turbulent kinetic energy) and ω (dissipation rating) are presented as follows:
ν t = a 1 k max ( a 1 ω ; Ω F 2 )
D ρ k D t = x j [ ( μ + σ k μ t k x j ) ] + τ i j u i x j β * ρ ω k
( τ i j = ρ u i u j ¯ )
D ρ ω D t = x j [ ( μ + σ ω μ t ω x j ) ] + γ v t τ i j u i x j β ρ ω 2 + 2 ( 1 F 1 ) ρ σ ω 2 1 ω k x j ω x j
where Ω is the absolute value of the vorticity, ρ is fluid density.
F 2 = tanh ( arg 2 2 )
arg 2 = max ( 2 k 0.09 ω y ; 500 ν y 2 ω )
where y is the distance to the wall.
The constants Φ are obtained from Φ 1 and   Φ 2   .
Φ = F 1 Φ 1 + ( 1 F 1 ) Φ 2
F 1 = tan h ( arg 4 1 )
arg 1 = min ( max ( k 0.09 ω y ; 500 v y 2 ω ) ; 4 ρ σ ω 2 k C D k ω y 2 )
C D k ω = max ( 2 ρ σ ω 2 1 ω k x j ω x j ; 10 20 )
The constants of Set Φ 1 are Wilcox:
a 1 = 0.31 , σ k 1 = 0.85 , σ ω 1 = 0.5 , β 1 = 0.075 , β * = 0.09 , κ = 0.41 , γ 1 = β 1 / β * σ ω 1 κ 2 / β *
The constants of Set Φ 2 are Jones–Launder:
a 1 = 0.31 , σ k 2 = 1 , σ ω 2 = 0.856 , β 2 = 0.0828 , β * = 0.09 , κ = 0.41 , γ 1 = β 2 / β * σ ω 2 κ 2 / β *

2.3. Boundary Conditions

The setting of the boundary conditions has a significant influence on solving physical problems of flow field. When the governing equation is given to solve different physical problems, the different boundary value conditions are adopted. The boundary conditions are presented as follows:
Velocity inlet:
The incoming boundary of the flow field, the boundary is usually set far from the ship. The incoming velocity is given, the pressure is obtained by using reconstructed gradient interpolation.
v = ( U , 0 , 0 )
d p d n = 0
The estimation formula of k and ω are:
k = 3 2 ( I k U ) 2
ω = ρ k μ ( μ t μ ) 1
α = { 1 z 0 0 z > 0
I k is turbulence intensity, μ t / μ is turbulent viscosity ratio, α is volume fraction.
Pressure outlet:
The outlet of the flow field in the computation domain, it is set far behind the stern. The variables of outlet boundary are generally unknown, but the boundary is usually far away from the hull, the physical variables change little. Except for the pressure, the normal gradient of other physical variables is zero:
u n = 0 , k n = 0 , ω n = 0 , α n = 0
Symmetry plane:
To reduce the total number of grids, only half computation domain is typically built, the symmetry plane has no physical meaning. The normal velocity of the boundary is zero, the gradient of the physical quantity of the flow is zero.
u n = n v = 0
u n = 0 , k n = 0 , ω n = 0 , α n = 0
Wall condition:
In the computation domain, only the hull needs wall conditions. The wall surface of the hull is non-slippable:
v = ( 0 , 0 , 0 )
d p d n = 0

2.4. Wall y+

In the numerical simulation of turbulence, especially with respect to the high Reynolds number problems, the hull is regarded as a non-slip wall, and because of the viscous damping of the wall, the velocity gradient near the hull is very large. In the numerical method, the wall function is used for the near wall treatment as a hybrid approach. In 2005, Wang [20] conducted research on the turbulence problem, indicating that the y+ value near the hull should be controlled between 30 and 300. The wall y+ value is dimensionless quantity, and the formula for calculating the y+ value is as follows:
y + = y υ 1 2 U 2 0.074 Re L 1 5
Here y is the distance from the wall to the boundary layer grid; U is the speed; v is the fluid viscosity coefficient; Re is the Reynolds number.
In the numerical simulation, the prism layer grid is generated to replace boundary layer. However, due to the intense motion of the planing craft, the waterline length will be sharply reduced as well. A low y+ value will greatly affect the accuracy, the value of y+ should be larger.

2.5. Volume of Fluid Method

In the numerical simulations of the sailing of vessel, the flow field is complex, tracking the free surface with a high accuracy method will ensure the accuracy of numerical simulation. To solve this problem, in 1981, Nichols and Hirt [21] proposed the volume of fluid method (VOF). The VOF method is to obtain and track the function F, which is defined as the volume ratio of different fluid occupancy in the grid. If the value of this function on each grid is obtained, the motion interface of the two phases of the flow will be tracked. The function F is defined as follows [22]:
F t + u F x + v F y + w F z = 0
Here u, v, w are the velocity components.

2.6. Overset Method

Overset grid is the method of region segmentation and grid combination, which consists of background and overset region. The overset method consists of two parts: (1) Divided grid cells in the whole flow field area into effective cells and invalid cells by hole cutting. Effective cell refers to the grid cell which participates in solving discrete governing equations, the grid cells that do not participate in computation are invalid cells. (2) Find the contributor cell, which configures the corresponding contributor cell for the recipient cell. The receiver cell layer is formed along the boundary of the overset region, and it is this layer of receiver cell layer that divides the background grid into effective cell region and invalid cell region. The overset region and the background region are closely coupled, and the governing equation is solved simultaneously in the effective cells of the two regions. The overset grid around the hull is presented in Figure 1.

3. Towing Tank Experiments

To validate the accuracy of the numerical method, a planing craft model (USV01) is selected as the study object. The seakeeping tests of USV01 were carried out in 2011 by Hailong Shen and colleagues [18]. The hull model and the specific parameters of the model are presented in Figure 2 and Table 1.
The setup of test is presented in Figure 3. All the experimental installations were installed under the carriage platform, and the guide rods are inserted into the guide plates, which were set at the bow and stern to ensure stability in the drag. In this way, the yaw and roll motions are prevented, and the heave and pitch motions of the model are free. The resistance is measured by a dynamometer in the front of the model, the pitch angle is measured by a gyroscope, which is set at the center of gravity, the heave is measured by the cable-extension displacement sensor, which is set on the above model, and the acceleration of the CG and bow of the model are measured by two acceleration sensors, which are set at the center of gravity and 0.275 m from bow. A snapshot of the model is presented in Figure 4.
In the regular wave test, the towed velocity (U) was 6 m/s (Fn = 1.16), the wave height is 0.21 m, and the wave length ranged from 1.5 L to 5 L. The regular wave test matrix is presented in Table 2, where the wave length is λ, the wave height is H, and the encounter frequency is ω . The motion parameters of the model were read when it crossed ten wave lengths in waves.
Under the same speed and wave height, based on the curves of the model’s motion, the peak–peak values of the motion response and the average resistance of the model in different conditions were obtained.

4. Numerical Simulation Validation

4.1. Computational Domains

For a better accuracy, the boundaries of the calculation domain were set far enough away to reduce the reflection of waves, and the entire computational domain should be no less than five times the wave length in the navigation direction. To improve the computing efficiency, only half of the model was built for the simulation. The definition of the boundary conditions and the size of domains are illustrated in Figure 5 and Table 3, respectively.
The point at the bottom of the stern is set as the origin. The distance between each boundary and the origin are presented in Table 3, and all the distances are dimensionalized in terms of the overall length of the model.
The definitions of boundary conditions in the domain are presented as follow: in the background region, the inlet, top and bottom are defined as the velocity inlet; the outlet is defined as the pressure outlet; and the sym and side are defined as the symmetry plane. In the overset region, all the boundaries of the region are defined as the overset grid, and the surface of the hull is defined as the nonslip wall. The definition of the boundary conditions and the setting of the calculation domains are illustrated in Figure 6.

4.2. The Verification of Overset Method

The heave and pitch motion of the model in waves is intense, the accuracy of the numerical method, based simply on the Dynamic Fluid Body Interaction, is reduced due to grid quality deterioration. Therefore, the overset grid is more appropriate in this case. The overset mesh method in this paper refers to the solution adopted in [8,18,23].

4.2.1. Grid Parameter

In the computational domains, the grids of two regions are refined: the region around the overset region and the free surface. The grids in the overset region are small, the size of the background region grid around the overset region needs to be consistent with that of the boundary of the overset region. And the grids around the free surface are also refined in order to track the free surface accurately, and the whole free surface is guaranteed to be included in the refinement region. The refinement regions are presented in Figure 7.
To verify the grid parameters, the sailing of USV01 in calm water was numerical simulated, the numerical results of resistance was compared with the test results, the model and the setting of the test is consistent with Section 3, the test condition is changed to calm water, towing speed is 8 m/s.
The grid parameter is presented in Table 4, four grids were adopted in the verification. Only in grid 1, the free surface is not refined, the refinement of the other three grids are different in the size. The grids on the hull and the free surface are presented in Figure 8.
The numerical results of the four grids are presented in Table 5, the resistance curves are presented in Figure 9.
Table 5 and Figure 9 indicate that, due to the lack of refinement grid around free surface, the resistance curve of grid 1 cannot converge, and the resistance curves of the other three grids show good convergence. Comparing the numerical results of grid 2 and 3, both the grids were refined, but the numerical results of the grid 2 show a large deviation, the size of the refinements of grid 2 is too large to ensure the accuracy of numerical results. It proves that the refinement grid around free surface is necessary in tracking the free surface. Moreover, the size of grid refinement should be small enough to ensure the accuracy of numerical simulation.
Comparing the numerical results of grid 3 and 4, the deviations and difference of the grids are both small. The two schemes of grids refinement can both guarantee the requirement of numerical simulation, but the size of grid 4 is too small, and the total number of the grids is too large, the computing efficiency will drop significantly.
Summing up the above, parameters of grid can satisfy the accuracy of numerical simulation and improve the computational efficiency at the same time. Subsequent grid settings are all referenced to that of grid 3.

4.2.2. Value of y+

As mentioned in Section 2.4, the value of y+ has influence on the accuracy of calculation, under normal conditions the value of y+ ranges from 30 to 300.
For verification, the influence of y+, the values of 50, 100, 200, 250, and 300 were adopted, the boundary layer grid generated according to different values of y+ are presented in Figure 10.
The Figure 10 indicates that the trim motion of the planing craft leads to a great change of the waterline of the vessel. The decrease of the waterline length makes the y+ value of hull less than the theoretical value during the sailing.
The numerical results and the resistance curves of different y+ are presented in Table 6 and Figure 11, respectively.
The Table 6 and Figure 11 indicate that resistance curves corresponding different y+ have similar convergence. The deviations of different values of y+ are all less than 10%, but the deviation which y+ = 250 is the smallest. Due to the intense motion of the planing craft, the waterline length sharply reduces. A low y+ value greatly affect the accuracy. In subsequent numerical simulations, the y+ value was set at 250.

4.2.3. Time Step of Iteration

The ∆t (time step) of iteration should be set at the value which allows the physical field to move no less than the distance of the minimum free surface mesh in each iteration, thus improving the calculation efficiency, with the aim of ensuring accuracy.
To verify the value of ∆t, four values: 0.001 s, 0.002 s, 0.006 s, 0.01 s, 0.02 s were adopted and the numerical results corresponding to different values were compared in Table 7 and the resistance curves of different y+ are presented Figure 12.
The comparison between results of different values of ∆t indicate that the convergence trends of the five resistance curves are consistent, but different values will affect the accuracy of numerical results. The deviations of the results become larger with the increasing of ∆t, this indicates a short time interval between two iterations will benefit the accuracy of the numerical results. But a too small value of ∆t can lead to a waste of computational efficiency. The results of ∆t = 0.001 s and 0.002 s are almost the same and the accuracy is high, ∆t = 0.002 s is adopted for its better accuracy and computational efficiency.
Summing up the above, the influence of the grid parameter, value of y+, and the time step of iteration on the accuracy and computational efficiency are verified. The results show that the numerical method has good convergence and high accuracy in simulating of the sailing of USV01 in calm water. Based on the overset method which has been verified, the numerical simulation validation for the seakeeping test of USV01 will be carried out.

4.3. Validation of Numerical Method

On the basis of the verification of overset method, the grid parameters refer to grid 3, the value of y+ is set at 250, the ∆t is set at 0.002 s. In this way both the accuracy and computational efficiency are guaranteed.
In the numerical simulation of seakeeping tests, the setting of calculation domains, grid parameters, and ∆t are referenced to the Section 4.1 and Section 4.2, meanwhile, in the wave conditions the flow field of the free surface is more complicated, the refinement region around the free surface needs to contain the entire wave surface and the vertical refinement grid should be not less than 10 layers. When the motion of the model in regular wave tends to stabilize and sails at least five wave lengths, the response amplitude operators of the model were obtained. The calculation grids of the numerical simulations of seakeeping tests are presented in Figure 13.
The 3D model in the numerical simulation is quite the same as the model of USV01, which is shown in Figure 1, and regular wave conditions (1–6; Table 4) were adopted in the simulation. The response amplitude operators (RAOs) of the model, including the amplitude values of the heave, pitch angle and acceleration of the CG and bow and the average resistance under head wave conditions, were numerically predicted. After the model sailed steadily in regular wave at least five wavelengths, the curves of the motion response and resistance were acquired. The peak peak values of the motion response can be accessed directly, the average resistance is obtained by averaging the convergence resistance curve. The comparisons between the numerical and test results are presented in Table 8 and Figure 14, respectively.
The RAOs of numerical results presented a similar variation tendency with those of the tested results: they reach the maximum values in the resonance wave length region and then decrease.
By comparing the RAOs of the model in different conditions, the resonance region of the pitch and heave is quite the same (the region is in condition 4 (λ/L = 3.5)). Before reaching the resonance region, the encounter frequency Is larger than the natural frequency of pitching and heave; the resonance phenomenon does not occur and the values of amplitude are far less than 1. In the resonance region, the encounter frequency is the same as the natural frequency of heave and pitch (the frequency is 6.446 rad/s (λ/L = 3.5)), and the amplitudes reach the max. After reaching the resonance region, the wave slope angle decrease as the wave length increase, causing the amplitudes to decrease, but they always remain around 1.
The resonance region of the acceleration of the CG and bow is different (the region is in condition 2 (λ/L = 2.25)). Under the short wave length conditions (λ/L ≤ 2.25), the vertical motion of the model is more intense due to the larger wave slope angle and wave disturbances, when the λ/L reaches 3.5, the pitch motion become more intense due to the resonance phenomenon. This indicates that the resonance region of the acceleration shifts to the shorter wave length direction compared with that of heave and pitch, the natural frequency of the vertical motion of the vessel is larger than that of pitch and heave, and the natural frequency of vertical motion is 9.261 rad/s (λ/L = 2.25). Under any condition, the amplitudes of the acceleration of the bow are significantly greater than that of CG. This indicates that the bow will slap on the free surface during navigation, and the excessive value of bow acceleration may threaten navigation safety.
The deviation of the numerical method is mainly due to the deformation of the grid and the complex flow fields around the model. The relative deviations of the numerical methods are greatly affected by the wave conditions. Before reaching the resonance region, the deviation between the calculated and test values is small. However, when the waves reach the resonance region, the amplitudes of the heave and pitch increased gradually, the motion of the model in the waves become more intense, and the deviation become larger.
The maximum deviations of the calculated heave, pitch, and acceleration of the CG and bow are 9.33%, 9.14%, 8.83%, and 8.48%, respectively. The average deviations are 4.95%, 4.79%, 4.81%, and 5.73%, respectively. All the deviations of the RAOs are less than 10%. The deviations of the numerical results of the RAOs all become larger in their resonance regions. In the short wave length condition (λ/L = 1.5 and 2.25), the short encountering period increases the frequency of the vertical motion, and the amplitudes of the acceleration increase sharply. This makes the nonlinear characteristics of the vertical motion increase sharply, and the deviation between the numerical and test results become larger in these conditions. As for the pitch and heave, the variation in the deviations is similar to that of the acceleration, but the wave length of their resonance region is larger.
In terms of average resistance, the values of average resistance are similar in any condition, the deviation between the numerical and test values is similar in different conditions, and the numerical results of the average resistance of the vessel in regular wave are all lower than those of the test results. The maximum value is 12.54%, the average value is 11.55%. Meanwhile, compared with the RAOs, the deviation of the resistance is greater. In wave conditions, the proportion of hull wave-making resistance in the total resistance is large as the flow field around the hull becomes more complicated by the wave disturbance. The grids around the hull cannot accurately track the free surface, the complex flow field around the hull brings deviation to the numerical prediction of wave-making resistance, this makes the numerical result of the resistance mean a certain deviation compared with the experimental value, but the variation tendency are similar.
Summing up the above, the numerical method accurately predicted the variation trends of the model’s motion response in regular wave. Compared with the corresponding test results, the relative deviations of the numerical results of the RAOs are less than 10%. This shows that the numerical method has a high accuracy in predicting the motion response of the model in regular wave. However, increasing nonlinear characteristics of the motion will make deviation of the numerical results become larger. Meanwhile, the complex flow field around the hull made it difficult to predict the wave-making resistance accurately and the deviation of the numerical results of the resistance are larger than that of the RAOs.

5. Numerical Simulation of Planing Craft in Regular Wave with Large Wave Height

Based on the numerical method which was been validated in Section 4, the low steaming of small planing craft under regular wave with large wave height is numerical analyzed.
The USV01 is designed for high-speed offshore navigation, its seakeeping is not enough to ensure safe navigation under regular wave with large wave height. The subsequent numerical simulation in large wave height conditions are not suitable to take USV01 as the simulation object.
A new small planing craft with better seakeeping is selected as the object of numerical simulation. The vessel’s length is 6 m, displacement is 1.8t, the deadrise angle at the stern is 24 deg, and the model scale is set at 1:3. The vessel has a larger deadrise angle and more drafts, in pitching and heave motion caused by wave disturbance, the hull can provide more damping.
To compare and validate the seakeeping performance of the new planing craft, the model of USV01 was scaled down to the same size and displacement as the new planing craft. Under the same wave conditions, the sailing of the two vessels under large wave height is predicted and compared.
The scaled model of USV01 is selected as prototype. The new vessel improved in seakeeping compared with the prototype, the new vessel is named as improved vessel.
The hull and specific parameters of the improved vessel are presented in Figure 15 and Table 9, respectively.
The model of USV01 is scaled down to the size shown in Table 10, but the shape remains the same.
The low steaming of the two vessels in the same conditions is numerically simulated. The regular wave conditions refer to moderate sea conditions. The wave height of the moderate sea conditions ranges from 1.25 m to 2.5 m. In the numerical simulation, the scale ratio is 1:3, the wave height of regular wave is set at 0.45 m, and the wave length range from 2 L to 6 L. The parameters of regular wave are shown in Table 11.

6. Results

6.1. The Influence of Wave Length on the Navigation Configuration

To study the seakeeping of the two vessels, the sailing of the prototype and improved vessel in conditions 7 to 11 were simulated, their motion responses were compared. The comparisons are presented in Table 12 and Figure 16, respectively.
Figure 16 indicates that the motion tendency of the prototype and improved vessel are similar: As the wave length increased, the amplitudes of the pitch and acceleration began to decrease gradually after reaching the resonance region, while the amplitudes of heave tend to remain stable.
The dimensionless amplitudes of the pitch and heave in the short wave length conditions (conditions 7 and 8) are less than 1. This indicates that, in these conditions, no resonance occurred during the sailing and the heave amplitudes are far lower than those of the wave amplitude. Meanwhile, the vessel’s pitching motion lag behind the sinusoidal motion of regular wave, the large wave height allows the prototype to plough through the wave before reaching the crest, the vessel is in the wave-piercing status. In these conditions, the heave and pitch amplitudes are small, while that of the acceleration are large, this indicates that the vertical motion is intense. The wave-piercing navigation is accompanied by a heavy green water, and the navigation safety will be threatened. The simulation of the wave-piercing navigation of the vessel is presented in Figure 17.
When the wave length reaches 4 L, the dimensionless amplitudes of the heave and pitch float around 1, indicating that in these conditions, the response of vessel to wave disturbance is better. The vessel starts to depart from the wave-piercing status. As the wave length increases, the encountering period increases, and the nonlinear characteristics of the motion begins to decrease. The vessel can normally pass through the crests and troughs of regular wave. At this time, the vessel is in a normal navigation configuration. However, the speed of the vessel in regular wave is too low to provide sufficient dynamic lift, and the vessel is in the hull-borne status. In this navigation configuration, the heave amplitude is greatly affected by the wave height and tends to be stable in the condition of an invariant wave height. The amplitudes of the pitch are greatly affected by the wave slope angle. Figure 18 shows the hull-borne status of the vessel. The model of the vessel can pass the crest smoothly, and the green water phenomenon did not occur.
The dimensions of the improved vessel are the same as those of the prototype, and the motion trends are consistent in the same condition, both of which have the wave-piercing status and hull-borne status.
In the wave-piercing status, the RAOs of the improved vessel are slightly larger than those of the prototype, and while the two vessels are similar in size, the large deadrise angle of the improved vessel makes its pitching period shorter, allowing the vessel to better respond to wave disturbances. Due to the large wave height and short encountering period, the improved vessel is also in the wave-piercing status, and because of the good response of the improved vessel to the waves, the phenomenon of green water on deck during navigation will be greatly reduced. Therefore, the seakeeping of the improved vessel is better than that of the prototype in wave-piercing status.
When the vessels are in the hull-borne status, the amplitudes of the heave and pitch become larger, compared with those in the wave-piercing status, and the amplitude of the acceleration decreases sharply. Due to the better seakeeping of the improved vessel, the hull is able to provide more pitch damping under the regular wave, compared with the prototype. After entering the hull-borne status, the RAOs of the improved vessel are all smaller than those of the prototype in the same wave conditions.
Summing up the above, the improved vessel provides a better seakeeping performance in both the wave-piercing and hull-borne statuses. In the wave-piercing status, the design of the improved vessel can effectively reduce the occurrence of wave-piercing and green water; in the hull-borne status, the RAOs of the vessel sailing in regular wave are reduced.

6.2. The Influence of Speed on the Navigation Configuration

To study the influence of speed on the RAOs of improved vessel, the sailing of the improved vessel at different speeds and with different wave lengths was simulated. Another speed (U = 3 m/s) is selected. The wave matrix is shown in Table 13.
The sailing of the improved vessel in conditions 7–11 and 12–16 were numerically simulated and compared. The amplitude values of the motion response of model are obtained, and a dimensionless treatment is also carried out. A comparison of the RAOs of the vessel at two speeds is presented in Table 14 and Figure 19, respectively.
As Figure 19 shows, there are still two kinds of navigation configurations: the wave-piercing and hull-borne statuses. Comparing the RAOs of the improved vessel in the same wave condition at different speeds, speed has an obvious effect on the RAOs of the vessel. In terms of heave, when λ/L > 3, the vessel is in the hull-borne status. Increasing the speed makes the hull generate more dynamic lift, and this greatly increase the heave amplitudes of the vessel. However, when the wave length decreases, the vessel goes into the wave-piercing status, and the effect of the speed on the heave is opposite to that in the hull-borne status. Meanwhile, the effect of the speed on pitch is consistent with that of the heave in any condition, but the effect is weaker. In terms of the acceleration, the acceleration amplitudes of t vessel at 3 m/s (Fn = 0.69) are larger than that of the vessel at 2 m/s (Fn = 0.46) in any condition, and the speed has a more obvious influence on the acceleration, compared with the heave and pitch.
For further study, the influence of speed on the RAOs of improved vessel in the wave-piercing and hull-borne statuses have been numerically predicted. The wave conditions where the wave height is 0.45 m, λ/L = 3 and 6 are selected, λ/L = 3, the improved vessel is in the wave-piercing status, and when λ/L = 6, the improved vessel is in the hull-borne status. The wave matrix is shown in Table 15. The RAOs of the improved vessel are obtained, based on the numerical method, and a comparison of RAOs at different speeds is presented in Figure 20 and Table 16.
As Figure 20 shows, the effect of increasing the speed is completely opposite to that of changing the navigation configuration. In the conditions where the wave length is 3 L, U = 2 m/s (Fn = 0.46), and the vessel is in the wave-piercing status, the amplitudes of the heave and pitch decrease sharply as the speed increases. When the wave length reaches 3 L, the period of the encounter decreases with the increasing of the speed, while the pitch period of the vessel remains unchanged. This means that, with the increasing of the speed, the gap between the encounter period and the resonance period increases, and the resonance phenomenon gradually becomes weaker, the amplitudes begin to decrease, the nonlinear characteristics of pitch increase, and the wave-piercing phenomenon of the vessel becomes more obvious, leading to more green water.
In the hull-borne status (λ/L = 6), the amplitudes of the heave and pitch both increase with the increasing of the speed. Compared with the wave-piercing status, amplitudes of the vessel in the hull-borne status are much higher than those in the wave-piercing status. With the increasing of the speed, the lift force received by the craft gradually increases, making the RAOs of the vessel in regular wave increase obviously.
In terms of the acceleration, the influence of the speed in different navigation configurations is quite the same, while in the wave-piercing status, the influence of the speed on the acceleration is more obvious. The amplitudes of the acceleration of the vessel increase with the increasing of the speed in any condition, and the increase will become sharp at a high speed. The amplitudes of the acceleration of the vessel in the wave-piercing status (λ/L = 3) are always larger than that in the hull-borne status (λ/L = 6). This shows that in the wave-piercing status, the pitching motion is more intense, and the nonlinear characteristics of the motion increase. The influence of the speed on the acceleration of CG and bow is consistent. However, the amplitude of the acceleration of the bow increases quite sharply as the speed continues to increase, indicating that the nonlinear characteristics of the motion of the bow at a high speed increase sharply, which will greatly endanger the navigation safety.
Summing up the above, increasing the speed will increase the nonlinear characteristics of the motion in any condition: in the wave-piercing status, the occurrence of wave-piercing and green water will become more serious, in the hull-borne status, the RAOs will increase sharply. This indicates that reducing the speed is the most effective means to ensure navigation safety under regular wave with large wave height.

7. Conclusions

This paper presents a seakeeping analysis of planing craft, based on numerical and experimental methods. In order to validate numerical method, the sailing of USV01 was numerically simulated, based on the RANS solver and overset method. The numerical results of the RAOs were compared with the corresponding test results, and the overall uncertainty was within 10%. The deviation of the numerical method became slightly larger in the resonance region, especially for the acceleration. These results indicate that the numerical method has a high accuracy in simulating the seakeeping test of planing craft. Additionally, the numerical prediction of wave-making resistance has a slightly larger deviation due to complex flow field around the hull.
For further study on the seakeeping performance of vessels in regular wave with large wave height, a new small size planing craft with better seakeeping performance is selected for numerical simulation, and the vessel is named as improved vessel, the model of USV01 shrunk to the same size as improved vessel and named as prototype. The low steaming of the prototype and improved vessel in the same conditions were numerically simulated and compared. Under the influence of the large wave height and different wave lengths, there are two navigation configurations in the sailing of the vessels: The wave-piercing and hull-borne status. In the wave-piercing status, the vessel’s pitching motion lags behind the sinusoidal motion of regular wave, and there is no resonance phenomenon. The pitch and heave amplitudes are small, and the vessel traverses the waves, accompanied by severely green water. In the hull-borne status, the amplitudes of the pitch and heave are affected by the wave conditions. Compared with the prototype, the improved vessel shows a better seakeeping performance in both navigation configurations: in the wave-piercing status, the phenomenon of green water will be greatly reduced with the improved vessel, and in the hull-borne status, the RAOs of the improved vessel are lower than those of the prototype in the same wave conditions. For further study, the influence of the speed on the RAOs of the improved vessel in different navigation configurations will be analyzed. The increasing of the speed increases the nonlinear characteristics of the vessel’s motion in regular wave and in both the wave-piercing and hull-borne statuses. However, in different navigation configurations, the effect of the speed on the RAOs will be different, and reducing speed is the most effective means to ensure the navigation safety of planing craft under regular wave with large wave height.

Author Contributions

Conceptualization, X.B., J.Z. and Y.S.; Formal analysis, J.Z. and Y.S.; Methodology, X.B.; Project administration, J.Z.; Software, X.B.; Validation, X.B.; Writing—original draft, X.B.; Writing—review & editing, J.Z. and Y.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financially funded by the National Natural Science Foundation of China (Grant No. 51409054), the Foundation Research Funds for the Central Universities (Grant No. HEUCF180101).

Acknowledgments

The author would like to express his gratitude to the reviewers for their suggestions and observations that helped in improving the present work.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Faltinsen, O.M. Hydrodynamics of High-Speed Marine Vehicles; National Defense Industry Press: Beijing, China, 2007. [Google Scholar]
  2. Dong, W.-C.; Guo, R.-X. State of the art of prediction on resistance of planing crafts. Ship Mech. 2000, 4, 68–81. [Google Scholar]
  3. Diez, M.; Serani, A.; Campana, E.F. CFD-based stochastic optimization of a destroyer hull form for realistic ocean operations. In Proceedings of the 14th International Conference on Fast Sea Transportation, Nantes, France, 27–29 September 2017; pp. 1–9. [Google Scholar]
  4. Campana, E.F.; Diez, M.; Liuzzi, G.; Lucidi, S.; Pellegrini, R.; Piccialli, V.; Rinaldi, F.; Serani, A. A multi-objective DIRECT algorithm for ship hull optimization. Comput. Optim. Appl. 2019, 71, 53–72. [Google Scholar] [CrossRef] [Green Version]
  5. AZCUETA, R. Computational of Turbulent Free Surface Flows Around Ships and Floating Bodies; Technical University Hamburg-Harburg: Hamburg, Germany, 2001. [Google Scholar]
  6. Caponnetto, M. Practical CFD simulations for planing hulls. In Proceedings of the Second International Euro Conference on High Performance Marine Vehicles, Hamburg, Germany, 2–5 May 2001; pp. 128–138. [Google Scholar]
  7. Weymouth, G.D.; Wilson, R.V.; Stern, F. RANS computational fluid dynamics predictions of pitch and heave ship motions in head seas. J. Ship Res. 2005, 49, 80–97. [Google Scholar]
  8. Sun, H.; Ma, W.; Zhu, J. Research on grid factor in numerical calculation of planing craft resistance. Ship Build. China 2015, 56, 170–178. [Google Scholar]
  9. Su, Y.-m.; Zhao, J.-x.; Chen, Q.-t.; Zhu, W.; Jin, X. Numerical simulation of the planing vessel in regular waves. Ship Mech. 2013, 17, 583–590. [Google Scholar]
  10. De Marco, A.; Mancini, S.; Miranda, S.; Scognamiglio, R.; Vitiello, L. Experimental and numerical hydrodynamic analysis of a stepped planing hull. Appl. Ocean Res. 2017, 64, 135–154. [Google Scholar] [CrossRef]
  11. Dashtimanesh, A.; Tavakoli, S.; Sahoo, P. A simplified method to calculate trim and resistance of a two-stepped planing hull. Ships Offshore Struct. 2017, 12 (Suppl. 1), S317–S329. [Google Scholar] [CrossRef]
  12. Ikeda, Y.; Katayama, T. Porpoising oscillations of very-high-speed marine craft. Phil. Trans. R. Soc. Lond. A 2000, 358, 1905–1915. [Google Scholar] [CrossRef]
  13. AZCUETA, R. Steady and unsteady ranse simulations for planing crafts. In Proceedings of the 7th International Conference on Fast Sea Transportation, Ischia, Italy, 7–10 October 2003; pp. 41–48. [Google Scholar]
  14. Begovic, E.; Bertorello, C.; Pennino, S. Experimental seakeeping assessment of a warped planing hull model series. Ocean Eng. 2014, 83, 1–15. [Google Scholar] [CrossRef]
  15. Begovic, E.; Bertorello, C.; Pennino, S.; Piscopo, V.; Scamardella, A. Statistical analysis of planing hull motions and accelerations in irregular head sea. Ocean Eng. 2016, 112, 253–264. [Google Scholar] [CrossRef]
  16. Jiang, Y.; Sun, H.; Zou, J.; Hu, A.; Yang, J. Experimental and numerical investigations on hydrodynamic and aerodynamic characteristics of the tunnel of planing trimaran. Appl. Ocean Res. 2017, 63, 1–10. [Google Scholar] [CrossRef]
  17. Jiang, Y.; Sun, H.; Zou, J.; Hu, A.; Yang, J. Analysis of tunnel hydrodynamic characteristics for planing trimaran by model tests and numerical simulations. Ocean Eng. 2016, 113, 101–110. [Google Scholar] [CrossRef]
  18. Bi, X.; Shen, H.; Zhou, J.; Su, Y. Numerical analysis of the influence of fixed hydrofoil installation position on seakeeping of the planing craft. Appl. Ocean Res. 2019, 90, 101863. [Google Scholar] [CrossRef]
  19. Menter, F.R. Improved Two-Equation Turbulence Models for Aerodynamic Flows; NASA Technical Memorandum: California, CA, USA, 1992.
  20. Wang, F. Computational Fluid Dynamics Analysis-Theory and Application of CFD Software; Tsinghua University Press: Beijing, China, 2005. [Google Scholar]
  21. Scardovelli, R.; Zaleski, S. Direct numerical simulation of free-surface and interfacial flow. Annu. Rev. Fluid Mech. 1999, 31, 567–603. [Google Scholar] [CrossRef] [Green Version]
  22. Hirt, C.W.; Nichols, B.D. Volume of Fluid (VOF) method for the dynamics of free boundary. J. Comput. Phys. 1981, 39, 201–225. [Google Scholar] [CrossRef]
  23. Bertorello, C.; Begovic, E.; Mancini, S. Hydrodynamic performances of small sizes swath craft. Brodogr. Shipbuild. 2015, 66, 1–22. [Google Scholar]
Figure 1. Background grid and overset grid.
Figure 1. Background grid and overset grid.
Water 12 01020 g001
Figure 2. Three-dimensional (3D) model of the USV01 model.
Figure 2. Three-dimensional (3D) model of the USV01 model.
Water 12 01020 g002
Figure 3. Experimental setup.
Figure 3. Experimental setup.
Water 12 01020 g003
Figure 4. Snapshot of the model.
Figure 4. Snapshot of the model.
Water 12 01020 g004
Figure 5. The size of the calculation domains.
Figure 5. The size of the calculation domains.
Water 12 01020 g005
Figure 6. Boundary conditions of the calculation domains.
Figure 6. Boundary conditions of the calculation domains.
Water 12 01020 g006
Figure 7. Refinement of grids.
Figure 7. Refinement of grids.
Water 12 01020 g007
Figure 8. Grids of the free surface: (a) Grid 1, (b) Grid 2, (c) Grid 3, (d) Grid 4.
Figure 8. Grids of the free surface: (a) Grid 1, (b) Grid 2, (c) Grid 3, (d) Grid 4.
Water 12 01020 g008
Figure 9. The resistance curve of different grids.
Figure 9. The resistance curve of different grids.
Water 12 01020 g009
Figure 10. Boundary layer grid: (a) y+ = 50, (b) y+ = 100, (c) y+ = 200, (d) y+ = 250, (e) y+ = 300.
Figure 10. Boundary layer grid: (a) y+ = 50, (b) y+ = 100, (c) y+ = 200, (d) y+ = 250, (e) y+ = 300.
Water 12 01020 g010
Figure 11. The resistance curve of different y+ s.
Figure 11. The resistance curve of different y+ s.
Water 12 01020 g011
Figure 12. The resistance curve of different ∆t.
Figure 12. The resistance curve of different ∆t.
Water 12 01020 g012
Figure 13. Calculation grids: (a) Domain grids; (b) Medium-profile local grids; (c) Surface grid on the hull; (d) Grid near free surface; (e) Free surface of regular wave.
Figure 13. Calculation grids: (a) Domain grids; (b) Medium-profile local grids; (c) Surface grid on the hull; (d) Grid near free surface; (e) Free surface of regular wave.
Water 12 01020 g013
Figure 14. Comparison of the numerical and tested results: (a) Amplitudes of the heave; (b) Amplitudes of the pitch; (c) Amplitudes of the acceleration of the CG; (d) Amplitudes of the acceleration of the bow; (e) Average values of the resistance.
Figure 14. Comparison of the numerical and tested results: (a) Amplitudes of the heave; (b) Amplitudes of the pitch; (c) Amplitudes of the acceleration of the CG; (d) Amplitudes of the acceleration of the bow; (e) Average values of the resistance.
Water 12 01020 g014
Figure 15. Threedimensional (3D) model and line plan of the improved vessel.
Figure 15. Threedimensional (3D) model and line plan of the improved vessel.
Water 12 01020 g015
Figure 16. The comparison of the RAOs of the prototype and improved vessel: (a) Amplitude of the acceleration of the CG; (b) Amplitudes of the acceleration of the bow; (c) Amplitudes of the heave; (d) Amplitudes of the pitch.
Figure 16. The comparison of the RAOs of the prototype and improved vessel: (a) Amplitude of the acceleration of the CG; (b) Amplitudes of the acceleration of the bow; (c) Amplitudes of the heave; (d) Amplitudes of the pitch.
Water 12 01020 g016
Figure 17. The simulation of the sailing of the improved vessel in the wave-piercing status: (a) Free surface; (b) Gas–liquid phase diagram; (c) Free surface representation of the hull.
Figure 17. The simulation of the sailing of the improved vessel in the wave-piercing status: (a) Free surface; (b) Gas–liquid phase diagram; (c) Free surface representation of the hull.
Water 12 01020 g017
Figure 18. The simulation of the sailing of an improved vessel in the hull-borne status: (a) Free surface; (b) Gas–liquid phase diagram; (c) Free surface representation of the hull.
Figure 18. The simulation of the sailing of an improved vessel in the hull-borne status: (a) Free surface; (b) Gas–liquid phase diagram; (c) Free surface representation of the hull.
Water 12 01020 g018
Figure 19. Comparison of the RAOs at two speeds: (a) Amplitudes of the acceleration of the CG; (b) Amplitudes of the acceleration of the bow; (c) Amplitudes of the heave; (d) Amplitudes of the pitch.
Figure 19. Comparison of the RAOs at two speeds: (a) Amplitudes of the acceleration of the CG; (b) Amplitudes of the acceleration of the bow; (c) Amplitudes of the heave; (d) Amplitudes of the pitch.
Water 12 01020 g019
Figure 20. Comparison of the RAOs under different navigation configurations: (a) Amplitudes of the acceleration of the CG; (b) Amplitudes of the acceleration of the bow; (c) Amplitudes of the heave; (d) Amplitudes of the pitch.
Figure 20. Comparison of the RAOs under different navigation configurations: (a) Amplitudes of the acceleration of the CG; (b) Amplitudes of the acceleration of the bow; (c) Amplitudes of the heave; (d) Amplitudes of the pitch.
Water 12 01020 g020
Table 1. Parameters of the model of USV01.
Table 1. Parameters of the model of USV01.
Main featureSymbolValue
Model scalek1:4
Length overallL2.75 m
Beam overallB0.78 m
Mouded depthh0.325 m
DisplacementΔ125.4 kg
Draftd0.1325 m
Longitudinal position of the centre of gravity L C G 1.048 m
Rotational inertiaJ58.6 (kg·m2)
Deadrise angleβ18 deg
Table 2. Regular wave test matrix of the seakeeping test.
Table 2. Regular wave test matrix of the seakeeping test.
No.U (m/s)FnH (m)λ/Lλ(m) ω   ( rad / s ) H/λ
161.160.211.54.12513.0310.051
261.160.212.256.18759.2610.034
361.160.2138.257.3010.025
461.160.213.59.6256.4460.022
561.160.214115.7870.019
661.160.21513.754.8160.015
Table 3. Sizes of the computational domains.
Table 3. Sizes of the computational domains.
BoundaryDistanceBoundary Conditions
Inlet2.5 LVelocity inlet
Outlet4 LPressure outlet
Top2 LVelocity inlet
Bottom1 LVelocity inlet
Side1.5 LSymmetry plane
Inlet of overset region1.25 LOverset grid
Outlet of overset region0.25 L
Top of overset region0.25 L
Bottom of overset region0.25 L
Side of overset region0.3 L
Table 4. Size of the grid.
Table 4. Size of the grid.
Grid parametersGrid 1Grid 2Grid 3Grid 4
Total150 k280 k650 k1400 k
Grid on hullRelative Size (% L)1110.5
Refinement grid on the
boundary of overset region
Relative Size (% L)42.1.51
Refinement Grid around overset regionRelative Size (% L)42.51.51
Refinement Grid around free surfaceRelative Size of grid in X Y Z direction (% L)N4
4
2
2.5
2.5
0.8
1.5
1.5
0.5
Table 5. Numerical results of four grids. CFD, computational fluid dynamics; EFD, experimental fluid dynamics.
Table 5. Numerical results of four grids. CFD, computational fluid dynamics; EFD, experimental fluid dynamics.
GridResistance/NDeviation (%)
CFDEFD
1N235.8N
2206.52235.812.41
3229.41235.82.71
4230.24235.82.35
Table 6. Numerical results of different y+.
Table 6. Numerical results of different y+.
y+Resistance/NDeviation (%)
CFDEFD
50223.06235.85.4
100220.82235.86.35
200221.78235.85.95
250229.41235.82.71
300228.56235.83.14
Table 7. Numerical results of different ∆t (time steps).
Table 7. Numerical results of different ∆t (time steps).
∆t (s)Resistance/NDeviation (%)
CFDEFD
0.001228.84235.83.11
0.002229.41235.82.71
0.006209.72235.811.06
0.01200.82235.814.83
0.02180.66235.823.38
Table 8. Numerical and test results of the response amplitude operators (RAOs) and resistance: (a) Heave and pitch; (b) Acceleration; (c) Resistance.
Table 8. Numerical and test results of the response amplitude operators (RAOs) and resistance: (a) Heave and pitch; (b) Acceleration; (c) Resistance.
(a)
Heave/mmPitch/deg
λ/LEFDCFDDeviation (%)EFDCFDDeviation (%)
1.53.8443.489.332.2612.212.25
2.258.5178.661.753.8444.168.22
312.96413.020.435.355.441.68
3.513.49214.195.245.5025.80455.49
413.34113.984.85.0515.1511.97
511.45612.398.153.9154.2739.14
(b)
Acceleration of Center of Gravity (CG)/gAcceleration of Bow/g
λ/LEFDCFDDeviation (%)EFDCFDDeviation (%)
1.50.9050.8258.832.111.9318.48
2.250.981.0032.342.3552.2056.36
30.8290.8492.411.8091.7055.74
3.50.6780.6376.041.3561.3510.36
40.5270.4946.260.9790.9255.51
50.3010.3102.990.5270.4857.97
(c)
Resistance/N
λ/LEFDCFDDeviation (%)
1.5240.806 213.993 11.13
2.25255.584 225.988 11.58
3246.725 223.342 9.47
3.5227.517 199.802 12.18
4221.607 193.844 12.52
5217.178 190.169 12.43
Table 9. Model parameters of the improved vessel.
Table 9. Model parameters of the improved vessel.
Main FeatureSymbolValue
Model scalek1:3
Length overallL2 m
Beam overallB0.54 m
Mouded depthh0.256 m
DisplacementΔ66.67 kg
Draftd0.147 m
Longitudinal position of the center of gravityLCG0.82 m
Deadrise angle β 24 deg
Table 10. Model parameters of the prototype.
Table 10. Model parameters of the prototype.
Main FeatureSymbolValue
Length overall L2 m
Beam overallB0.67 m
Mouded depth h0.23 m
Displacement Δ66.67 kg
Draft d0.118 m
Longitudinal position of the center of gravityLCG0.82 m
Deadrise angle β 18 deg
Table 11. The matrix of regular wave with large wave height.
Table 11. The matrix of regular wave with large wave height.
NoU (m/s)FnH (m)λ/Lλ (m) ω   ( rad / s ) H/λ
720.460.45247.316 0.1125
820.460.45365.391 0.075
920.460.45484.390 0.05625
1020.460.455103.749 0.045
1120.460.456123.324 0.0375
Table 12. RAOs of the prototype and improved vessel: (a) Heave and pitch; (b) Acceleration
Table 12. RAOs of the prototype and improved vessel: (a) Heave and pitch; (b) Acceleration
(a)
Amplitude of Heave (m)Amplitude of Pitch (deg)
λ/Lprototypeimproved vesselprototypeimproved vessel
20.1030.1186.8257.895
30.2130.2098.76410.685
40.2280.2099.6539.44
50.2330.2088.9427.835
60.2350.2077.4836.195
(b)
Amplitude of Acceleration of CG (g)Amplitude of Acceleration of Bow (g)
λ/Lprototypeimproved vesselprototypeimproved vessel
20.5740.6511.0171.278
30.6770.6591.4111.309
40.5310.5030.8710.783
50.3980.3670.6170.549
60.2720.2510.4510.347
Table 13. The parameters of regular wave, in which U = 3 m/s.
Table 13. The parameters of regular wave, in which U = 3 m/s.
NoU (m/s)FnH (m)λ/Lλ (m) ω   ( rad / s ) H/λ
1230.690.45248.897 0.1125
1330.690.45366.440 0.075
1430.690.45485.177 0.05625
1530.690.455104.392 0.045
1630.690.456123.845 0.0375
Table 14. RAOs of the improved vessel at different speeds and with different wave lengths: (a) Heave and pitch; (b) Acceleration.
Table 14. RAOs of the improved vessel at different speeds and with different wave lengths: (a) Heave and pitch; (b) Acceleration.
(a)
Amplitude of Heave (m)Amplitude of Pitch (deg)
λ/LFn = 0.46Fn = 0.69Fn = 0.46Fn = 0.69
20.1180.0977.8956.597
30.2090.23210.68510.434
40.2090.2469.4410.225
50.2080.2447.8358.673
60.2070.2426.1957.465
(b)
Amplitude of Acceleration of CG (g)Amplitude of Acceleration of Bow (g)
λ/LFn = 0.46Fn = 0.69Fn = 0.46Fn = 0.69
20.6510.7711.2782.059
30.6590.9181.3092.709
40.5030.7520.7831.403
50.3670.6120.5490.783
60.2510.4810.3470.574
Table 15. The regular wave matrix at different speeds.
Table 15. The regular wave matrix at different speeds.
NoU (m/s)FnH (m)λ/Lλ (m) ω   ( rad / s ) H/λ
1720.460.45365.3910.075
1830.690.45366.4400.075
1940.920.45367.4890.075
2051.150.45368.5390.075
2161.390.45369.5940.075
2220.460.456123.3240.0375
2330.690.456123.8450.0375
2440.920.456124.3670.0375
2551.150.456124.8880.0375
2661.390.456125.4110.0375
Table 16. The RAOs of the improved vessel at different speeds: (a) Heave and pitch; (b) Acceleration.
Table 16. The RAOs of the improved vessel at different speeds: (a) Heave and pitch; (b) Acceleration.
(a)
Amplitude of Heave (m)Amplitude of Pitch (deg)
Fnλ/L = 3λ/L = 6λ/L = 3λ/L = 6
0.460.2090.20710.6856.159
0.690.2320.24210.4347.465
0.920.1960.2558.7457.651
1.150.1710.27046.9518.408
1.390.1420.30356.1159.483
(b)
Amplitude of Acceleration of CG (g)Amplitude of Acceleration of Bow (g)
Fnλ/L = 3λ/L = 6λ/L = 3λ/L = 6
0.460.6590.2511.3090.347
0.690.9180.4812.7090.574
0.921.5830.5574.0220.766
1.152.1930.7124.7881.276
1.392.4361.0384.8712.381

Share and Cite

MDPI and ACS Style

Bi, X.; Zhuang, J.; Su, Y. Seakeeping Analysis of Planing Craft under Large Wave Height. Water 2020, 12, 1020. https://doi.org/10.3390/w12041020

AMA Style

Bi X, Zhuang J, Su Y. Seakeeping Analysis of Planing Craft under Large Wave Height. Water. 2020; 12(4):1020. https://doi.org/10.3390/w12041020

Chicago/Turabian Style

Bi, Xiaosheng, Jiayuan Zhuang, and Yumin Su. 2020. "Seakeeping Analysis of Planing Craft under Large Wave Height" Water 12, no. 4: 1020. https://doi.org/10.3390/w12041020

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