An Experimental Study on Hydrodynamic Forces of Korea Autonomous Surface Ship in Various Loading Condition

Article information

J Navig Port Res. 2022;46(2):73-81
Publication date (electronic) : 2022 April 30
doi : https://doi.org/10.5394/KINPR.2022.46.2.73
*Ph.D. Candidate, Graduate School of Changwon National University, Gyeongsangnam-do, 51140, Korea
**,***Master’s Student, Graduate School of Changwon National University, Gyeongsangnam-do, 51140, Korea
Professor of Changwon National University, Gyeongsangnam-do, 51140, Korea
Corresponding author, hkyoon@changwon.ac.kr 055)213-368
Note) This paper was presented on the subject of "An Experimental Study on Effect of Loading Conditions on Hydrodynamic Forces Acting on a Ship" in Asian Navigation Conference (ANC) proceeding(Tianjin, China 18th-20th October, 2021).
Received 2021 November 19; Revised 2021 December 13; Accepted 2022 March 26.

Abstract

Currently, shipping by sea is becoming common because of the low price and the safety of goods. The ship is designed as a larger vessel to meet the need of this development. In the design stage, the investigation of hydrodynamic forces acting on the ship hull is very important in predicting the ship’s maneuverability. Given that the ship docks at various ports for loading or discharging goods, the ship usually operates in various loading conditions, depending on the site condition and other various factors. Hence, it is necessary to investigate the effect of the loading condition on the hydrodynamic forces acting on the ship, to most accurately determine the maneuverability of the ship. In this study, an experiment of Korea Autonomous Surface Ship (KASS) was conducted at the towing tank of Changwon National University to measure the hydrodynamic forces acting on the KASS. The loading condition considered in this experiment is determined based on the draft, which was decreased by 5% for each loading condition. The smallest draft is 85% of the design draft. The static test as Oblique Towing Test (OTT), Circular Motion Test (CMT), Circular Motion Test with Drift (CMTD) is performed in the various loading conditions. First, the hydrodynamic forces in the Oblique Towing test (OTT) are compared with the result of other institutes. Second, the hydrodynamic forces in various drift angle, yaw rate and loading conditions are measured. Finally, the influence of the loading conditions on the hydrodynamic coefficient is discussed.

1. Introduction

In recent years, shipping has increased due to the development of the industry and freight transport. To reduce the price of shipping, the ship is designed bigger to transport more goods at a time. However, the ship berths at various ports to loading and discharging goods during its voyage. Therefore, the ship usually operate in various loading conditions depending on site condition and other factors. To ensure that the ship has good maneuverability in given loading conditions such as the turning ability and course keeping ability, it is necessary to investigate the hydrodynamic forces acting on the ship hull.

In the past, researchers were researching the hydrodynamic forces acting on the ship hull by performing numerical simulation or captive model test. Lee (2003) predicted the hydrodynamic forces in laterally berthing maneuvers using Computational Fluid Dynamics (CFD). Park (2018) investigated the effect of the vertical center of gravity on the hydrodynamic force acting on the ship hull. The numerical method was performed in various drift angles to estimate the hydrodynamic forces. Nguyen (2019) studied the hydrodynamic forces acting on the ship hull with and without a propeller by performing a model test. In addition, the hydrodynamic forces in the regular waves were estimated. Dai (2019) estimated the hydrodynamic forces and hydrodynamic derivatives of the KVLCC2 by a virtual captive model test. The Oblique Towing Test (OTT) and Circular Motion Test (CMT) were carried out to estimate the hydrodynamic forces. Wang (2011) estimated the hydrodynamic forces in oblique motion through experiment and numerical method. The effect of the drift angle was investigated. Islam (2018) estimated the linear hydrodynamic derivatives using static drift simulation. The result of the simulation of hydrodynamic forces was compared with experimental results. RANS based solver and SHIP_Motion were used for the simulation. Yang (2019) investigated the hydrodynamic forces of the KVLCC2 model in various drift conditions by numerical method. In addition, the effect of the mesh conditions on the hydrodynamic forces was considered in that study.

In our study, the effect of loading conditions on hydrodynamic forces acting on the ship hull in various loading conditions is investigated. After that, the hydrodynamic derivatives are estimated based on the estimated hydrodynamic forces obtained by performing experiments for each condition. The Korea Autonomous Surface Ship (KASS) was selected to perform captive model test. This test is carried out at towing tank of Changwon National University to measure the hydrodynamic forces acting on the KASS. In actual ship, the loading condition was determined based on loading as full loading, half loading or ballasting. However, in this study, the loading condition considered in this study is determined based on the draft, which is decreased by 5 percent for each loading condition. The static test as Oblique Towing Test (OTT) and Circular Motion Test (CMT) and Circular Motion Test with Drift (CMTD) are performed in the various loading conditions. OTT and CMT are performed to measure the hydrodynamic forces acting on the KASS model in oblique moving and steady turning, respectively. The hydrodynamic derivatives related to sway velocity and yaw velocity are estimated based on OTT and CMT results. In order to determine the coupling coefficients of sway velocity and yaw velocity, the CMTD is performed. The hydrodynamic forces in the OTT are compared with the result of other institutes in the case of the full loading condition. Then, the effect of the drift angle and yaw rate on the hydrodynamic forces are evaluated. The changing of the hydrodynamic forces in various loading conditions is measured and the influence of the loading conditions on the hydrodynamic coefficients is discussed. The hydrodynamic coefficients obtained from this research can be used to perform maneuvers simulation of the same ship type as KASS. The standards of maneuvers can be checked in this state.

2. Test facility and test conditions

2.1 Test facility

The experiment of KASS was performed in a towing wave tank in Changwon National University. Two load cells with 2-components force transducer were used to measure the hydrodynamic force acting on the KASS’s hull. The model of the load cell is MCL-2A01-60N and the capacity of the load cell is ±60N for both the force in-x direction and y-direction. Load cells that were used in this test are shown in Fig. 1. Fig. 2 shows the towing tank used in this experiment. Towing tank dimensions are 14 m breadth, 20 m long and 1.8 m deep. The limit of the carriage speed is 1.0 m/s. In this experiment, the water deep was set at 1.0 m.

Fig. 1.

Two components load cell

Fig. 2.

Towing tank

2.2 Model ship and test conditions

The KASS model was selected to investigate the effect of the loading condition on hydrodynamic forces. The model test was performed in various drafts to obtain hydrodynamic coefficients of the KASS model. The model test has a scale ratio of 1/11. Fig. 3 shows the KASS model used in this experiment. The principal dimensions of the both full scale and model scale are listed in Table 1. The experiments were carried out at the condition of the model scale with ship speed set at 0.931 m/s (equivalent to 6 knots in full scale). The model test was performed in the OTT, CMT CMTD to investigate the effect of the drift angle and yaw rate on hydrodynamic forces. In addition, the effect of the loading condition on hydrodynamic forces is investigated. The loading conditions are listed in Table 2. Table 3 lists the test conditions in this experiment.

Fig. 3.

KASS model

Principal particular of initial hul form of KASS

Loading conditions

Test conditions

3. Experiment

3.1 Equation of motion

In this experiment, the ship’s maneuvering mathematical model considers three degrees of freedom. Therefore, the surge force, sway force and yaw moment are measured. Two right-handed coordinate systems were used to describe the hydrodynamic force acting on the ship hull. A body-fixed coordinate system and an earth-fixed coordinate system are shown in Fig. 4. The mathematical model can be described by the following surge motion, sway motion and yaw motion by Eq. (1).

(1) {X=m(u˙-vr-xGr2)Y=m(v˙+ur+xGr˙)N=Izz+mxG(v˙+ur)
where X , Y and N are denoted for hydrodynamic forces acting on the KASSs hull as surge force, sway force and yaw moment, respectively. m, Izz and xG are the mass of KASS, moment of inertia in yaw motion and the longitudinal center of gravity of the ship, respectively. u and v is the velocity in x-axis and y-axis, respectively. is the acceleration of angular.

Fig. 4.

Coordinate system

In order to estimate the effect of the loading condition on the hydrodynamic coefficients, the hull forces are modeled based on the Maneuvering Modeling Group (MMG) model for ship maneuvering (Yasukawa, 2014). In addition, the 3rd order polynomial function is replaced by the 2nd order polynomial function and absolute value of sway velocity and yaw rate to determine the hydrodynamic coefficients as Yv|v|, Yr|r|, Nv|v| and Nr|r|. The mathematical model of the hull forces is described in Eqs. (2)– (4). The prime represents for non-dimensional value.

(2) X=Xvvv2+Xrrr2+Xvrvr
(3) Y=Yvv+Yv|v|v|v|+Yrr+Yr|r|r|r|+Yvvrv2r+Yvrrvr2
(4) N=Nvv+Nv|v|v|v|+Nrr+Nr|r|r|r|+Nvvrv2r+Nvrrvr2
where X′vv, X′rr and X′vr are the surge derivatives. Y′v, Y′v|v|, Y′r, Y′r|r|, Y′vvr and Y′vrr are the sway derivatives. N′v, N′v|v|, N′r, N′r|r|, N′vvr and N′vrr are the yaw derivatives. v′ and r′ are non-dimensional sway velocity and angular velocity, respectively.

3.2 Ballasting test and inertia test

Before performing the experiment, the mass distribution of the ship was made identical to the real scale. Fig. 5 shows the ballasting test of the KASS model ship in various loading conditions. In addition, the inertia table was used to obtain the mass moment of inertia Izz. Fig. 6 shows the inertia test in full loading condition. The results of the inertia test are listed in Table 4. The differences of measured Izz and target Izz are smaller than 5%. This error value can be accepted in captive model test.

Fig. 5.

Ballasting in various loading conditions

Fig. 6.

Inertia test in case 1 (Full loading condition)

Result of inertia test in various loading conditions

3.3 Experimental setup

The KASS model test was carried out at the towing tank in Changwon National University. The load cells are installed on the strong back near the bow and stern of the model ship. The potentiometers are connected between load cells and the measuring frame which is then connected with the sub-carriage in the yaw table. During the test, the model test is towed by the towing carriage. The yaw table can change the drift angle by rotating around its center. In the case of the circular test, the yaw table rotates with the yaw rate that was set. The signal obtained from the load cell is then transferred from the electrical signal into the digital signal by an A/D converter NI USB-6212 was made by National Instruments Corporation Company. The surge force and sway forces are measured by the total surge force and sway force obtained from the load cells. In addition, the yaw moment is obtained by the sway force of the load cells and the distance from the center of these load cells to the center of the model ship. Fig. 7 shows the real experimental setup of the KASS.

Fig. 7.

Real experimental setup in the case of full loading condition

3.4 Data analysis

The measured data is recorded as time series when the carriage speed reaches the target speed of 0.931 m/s while stable. The time series are constantly checked after each run to avoid obvious errors. The hydrodynamic forces obtained from two load cells are estimated by the mean value. X , Y and N are the hydrodynamic forces acting on the KASSs hull as surge force, sway force and yaw moment, respectively. X′, Y′ and N′, are the non-dimensional of hydrodynamic forces acting on the KASSs hull as surge force, sway force and yaw moment, respectively. The definition of the non-dimensional hydrodynamic forces is shown in Eqs. (5) – (7).

(5) X=X0.5ρL2U2
(6) Y=Y0.5ρL2U2
(7) N=N0.5ρL3U2

4. Result

4.1 Validation result

The experimental results of the sway force acting on the ship hull in the case 1 (full loading condition) are compared with the results from the Korean Research Institute of Ship and Ocean Engineering (KRISO) (Kim, 2021) as shown in Fig. 8. Sway force results are similar to the results from KRISO. In the small drift angle, the error of the model’s hull is barely affected. When the drift angle increase, the effect of the error on the model’s hull rapidly increases. In the large drift angle, there is a slight difference in sway force due to the model’s hull. In experiment, many thing can be affect to the result. The length between perpendiculars of the KRISO’s model ship is 5.789 m. The length of model which used in this experiment is 2 m. In addition, the boundary layer of the KRISO’s model is larger than our model ship. Therefore, there may a little difference between KRISO’s result and result of the present study in large drift angle due to the scale effect.

Fig. 8.

Comparison sway force in OTT with KRISO’s results

4.2 Experimental results

In the oblique towing test, the effect of drift angle is investigated. Besides the effect of drift angle, the effect of the loading condition is the main investigated in this study. The hydrodynamic forces acting on the model ship hull in various drift angles and various loading conditions are shown in Fig. 9. The hydrodynamic forces increase slightly in various loading conditions and its effect is shown clearly in the large drift angle. The hydrodynamic forces are large in the case of full loading condition and reduce which decrease in the loading condition. The hydrodynamic forces are affected by the fluid’s pressure around the ship’s hull. When the loading condition decreases, the wetted area also decreases. Therefore, the fluid’s pressure also decreases due to the decreasing of the wetted area. The hydrodynamic forces acting on the ship hull increase significantly with increase in the drift angle. The results of the hydrodynamic forces in various yaw rates and various loading conditions are shown in Fig. 10. The loading condition also affects the hydrodynamic forces in CMT. In the case of full loading condition, the hydrodynamic forces increase dramatically and nearly 2 times of case 4. In addition, when the yaw rate increases, the hydrodynamic forces increase significantly. The changing of hydrodynamic forces in various loading conditions is affected by the wetted area and wetted volume. The wetted area and wetted volume of case 1 are 1.174 m2 and 0.0651 m3. The wetted area and wetted volume of case 4 are 0.993 m2 and 0.0509 m3. The wetted area and wetted volume decrease due to the decreasing of loading condition. In addition, the hydrodynamic forces acting on the ship hull is the water pressure acting on the wetted area of ship. Therefore, the hydrodynamic forces decrease in lower loading condition.

Fig. 9.

Results of OTT in various loading conditions

Fig. 10.

Results of CMT in various loading conditions

The effect of the drift angle, yaw rate and loading condition on hydrodynamic forces are shown in Fig. 11. The loading condition effect on small yaw rate is similar to trend of CMT’s results. However, the hydrodynamic forces increase in large yaw rate due to the effect of the loading conditions. Moreover, the hydrodynamic forces reduce and become a positive value when the drift angle increases.

Fig. 11.

Results of CMTD in various loading conditions

4.3 Hydrodynamic coefficients

The hydrodynamic coefficients on maneuvering are determined using the OTT, CMT and CMTD results. The centrifugal force is removed in the hydrodynamic forces measured in both CMT and CMTD. Then, the hydrodynamic coefficients are estimated using the expression of the 2nd order polynomial function. X′vv, Y′v, Y ′v|v|, N ′v and N ′v|v| coefficients are determined based on the OTT results. X ′rr, Y ′r, Y ′r|r|, N ′r and N ′r|r| coefficients are determined based on the CMT results. Y ′vvr, Y ′vrr, N ′vvr and N′vrr coefficients are determined based on the CMTD results and the hydrodynamic coefficients which are estimated from OTT and CMT. The hydrodynamic coefficients are obtained by Least Square Method. Fig. 12 shows the fitting value and measured value to obtain the Xvv, Yv, Yv|v|, Nv and Nv|v| coefficients in various loading conditions. In the case of OTT, the effect of drift angle on hydrodynamic forces are investigated. Fig. 13 shows the fitting value and measured value used to obtain the Xrr, Yr, Yr|r|, Nr and Nr|r| coefficients in various loading conditions. Fig. 14 shows the fitting value and measured value to obtain the Yvvr, Yvrr, Nvvr and Nvrr coefficients in various loading conditions. To confirm the accuracy of the expression, the fitting curves expressed as the solid line and symbol are the measured value are plotted in Figs. 1214. The fitting curves accuracy is a good agreement with the measured value. The hydrodynamic coefficients depend on the loading condition obtained from the Least Square Method are listed in Table 5. The hydrodynamic coefficients change dramatically due to the effect of the loading condition.

Fig. 12.

Fitting value of OTT in various loading conditions

Fig. 13.

Fitting value of CMT in various loading conditions

Fig.14

Fitting value of CMTD in various loading conditions

Hydrodynamic coefficients

Almost the hydrodynamic coefficients slightly reduce with decrease in the loading condition. Yv and Nv coefficient slightly decreases when the loading condition reduces. However, Yr and Nr coefficient changes dramatically in various loading conditions. In addition, the linear coefficients as Yvlinear, Nvlinear, Yrlinear and Nrlinear are derived from the drift angle −6 deg. to 6deg. and non-dimensional yaw rate smaller than 0.5 to checking the tendencies of the hydrodynamic coefficients. The linear coefficients reduce when the loading condition decreases as listed in Table 5.

5. Conclusion

In our study, the model test of KASS was carried out in towing tank at Changwon National University and the effect of the loading condition on hydrodynamic forces in various drift angles and various yaw rates were investigated. The concluding remarks are as follows:

First, the hydrodynamic forces slightly increase in various loading conditions. The effect of the loading condition is shown clearly in the large drift angle. The hydrodynamic forces are increase dramatically in the case of full loading condition and decreases when the loading condition decreases.

Second, the hydrodynamic forces change dramatically in various yaw rates due to the effect of loading condition. In the small yaw rate, the hydrodynamic forces change slightly in various loading conditions. However, the hydrodynamic forces increase significantly in various loading conditions especially in the largest yaw rate while in the case of full loading condition, the hydrodynamic forces increase dramatically and nearly 2 times of case 4.

Third, the effect of the loading condition on the hydrodynamic forces in various yaw rates and various drift angles was investigated. The loading condition effect in small yaw rate is identical to the trend in CMT’s results.. However, the hydrodynamic forces increase significantly in large yaw rate due to the effect of loading conditions. Moreover, the hydrodynamic forces reduce and become positive value when the drift angle increases.

Lastly, the hydrodynamic coefficients depending on the loading conditions are obtained based on the experiment results. The hydrodynamic coefficients change dramatically due to the effect of the loading conditions. The hydrodynamic coefficients slightly reduce when the loading condition decreases.

Acknowledgements

This research was partly supported by the ‘Development of Autonomous Ship Technology (PJT201313, Development of Autonomous Navigation System with Intelligent Route Planning Function)’ funded by the Ministry of Oceans and Fisheries (MOF, Korea) and by Institute for Information & communications Technology Promotion(IITP) grant funded by the Korea Government(MSIT) and Korea Institute for Advancement of Technology(KIAT) grant funded by the Korea Government(MOTIE) (N0009999, The Competency Development Program for Industry Specialist).

References

1. Dai K, Li Y. 2019;Manoeuvering prediction of KVLCC2 with hydrodynamic derivatives generated by a virtual captive model test. Journal of Polish Maritime Research 26(4):16–26.
2. Islam H, Rahaman MM, Afroz L, Akimoto H. 2018;Estimation of linear hydrodynamic derivatives of a VLCC using static drift simulation. Proceeding of AIP conference :1–7.
3. Kim DJ, Choi HJ, Kim YG, Yeo DJ. 2021;Mathematical model for harbor manoeuvers of Korea Autonomous Surface Ship (KASS) based on captive model test. Conference of the Korean Association of Ocean Science and Technology Societies :1–8.
4. Lee YS. 2003;The prediction of hydrodynamic forces acting on ship hull in laterally berthing maneuver using CFD. International Journal of Navigation and Port Research 27(3):253–258.
5. Nguyen VM, Seo JW, Yoon HK, Kim YG. 2019;Experimental study on hydrodynamic forces acting on ship hull and rudder behind the propeller in regular waves. Proceeding of 11th International Work shop on Ship and Marine Hydrodynamics :1–12.
6. Park TC, Lee SW, Paik KJ, Moon SH. 2018;Study on hydrodynamic forces acting on tanker hull with consideration of various vertical center of gravity in drift test. Journal of Ocean Engineering and Technology 32(6):433–439.
7. Wang HM, Xie YH, Liu JM, Zou ZJ, He W. 2011;Experimental and numerical study of hydrodynamic forces on ship in oblique motion. Proceeding of International Conference on Remote Sensing, Environment and Transportation Engineering :328–331.
8. Yang X, Yin Y, Lian JJ. 2019;A numerical study on flow field and maneuvering derivatives of KVLCC2 model at drift condition. Journal of Marine Engineering and Technology :1–12.
9. Yasukawa H, Yoshimura Y. 2014;Introduction of MMG model standard method for ship maneuvering predictions. Journal of Marine Science and Technology 20:37–52.

Article information Continued

Fig. 1.

Two components load cell

Fig. 2.

Towing tank

Fig. 3.

KASS model

Fig. 4.

Coordinate system

Fig. 5.

Ballasting in various loading conditions

Fig. 6.

Inertia test in case 1 (Full loading condition)

Fig. 7.

Real experimental setup in the case of full loading condition

Fig. 8.

Comparison sway force in OTT with KRISO’s results

Fig. 9.

Results of OTT in various loading conditions

Fig. 10.

Results of CMT in various loading conditions

Fig. 11.

Results of CMTD in various loading conditions

Fig. 12.

Fitting value of OTT in various loading conditions

Fig. 13.

Fitting value of CMT in various loading conditions

Fig.14

Fitting value of CMTD in various loading conditions

Table 1.

Principal particular of initial hul form of KASS

Item Unit Full scale Model scale
Scale ratio - 1 1/11
Length, L m 22.000 2.000
Breadth, B m 6.000 0.545
Draft, T m 1.250 0.114
Depth, D m 4.138 0.376
Displacement, Δ kgf. 88873.641 66.772
Froude number, Fr - 0.210 0.210
Ship speed, U knots 6.000 1.809
Ship speed, U m/s 3.087 0.931

Table 2.

Loading conditions

Item Unit Case of loading condition
Case 1 (100%T) Case 2 (95%T) Case 3 (90%T) Case 4 (85%T)
Draft, T m 0.114 0.108 0.102 0.097
Dis., Δ kgf. 66.772 61.682 56.805 52.184

Table 3.

Test conditions

Item Variable Value
Oblique Towing Test (OTT) β[°] 0, ±3, ±6, ±9, ±10 ±12, ±15
Circular Motion Test (CMT) r′[−] 0.2, 0.3, 0.4, 0.5, 0.6
CMT with Drift (CMTD) r′[−] 0.2, 0.3, 0.4, 0.5, 0.6
β[°] 3, 6, 9, 12, 15

Table 4.

Result of inertia test in various loading conditions

Item Unit Case of loading condition
Case 1 Case 2 Case 3 Case 4
Target Izz kgm2 16.693 15.420 14.201 13.046
Measured Izz kgm2 16.969 15.061 13.803 13.574
Difference % 1.626 2.387 2.807 4.050

Table 5.

Hydrodynamic coefficients

Item Case 1 Case 2 Case 3 Case 4
X′vv 1.63E-03 3.44E-04 8.79E-04 4.03E-04
X′rr 1.90E-02 1.82E-02 1.81E-02 1.59E-02
X′vr 1.44E-02 1.22E-02 1.25E-02 1.39E-02
Y′v −1.44E-02 −1.43E-02 −1.31E-02 −1.25E-02
Y′v|v| −7.92E-03 −4.00E-03 −5.21E-03 −2.32E-03
Y′r 1.28E-03 −3.87E-03 −3.77E-03 −1.80E-03
Y′r|r| −2.16E-02 −1.01E-02 −9.76E-03 −1.02E-02
Y′vvr −9.74E-02 −3.60E-02 −7.06E-02 −5.35E-02
Y′vrr −8.86E-02 −4.76E-02 −6.12E-02 −4.62E-02
N′v −4.02E-03 −3.72E-03 −3.30E-03 −3.05E-03
N′v|v| −2.22E-03 −1.98E-03 −2.37E-03 −1.44E-03
N′r −3.30E-03 −3.20E-03 −3.98E-03 −3.73E-03
N′r|r| −8.57E-03 −6.79E-03 −5.17E-03 −4.23E-03
N′vvr −2.78E-02 −1.64E-02 −2.49E-02 −2.03E-02
N′vrr −2.77E-02 −1.90E-02 −2.20E-02 −2.05E-02
Yv'linear −1.53E-02 −1.45E-02 −1.33E-02 −1.28E-02
Yr'linear −1.39E-02 −1.12E-02 −1.07E-02 −8.90E-03
Nv'linear −4.22E-03 −3.83E-03 −3.48E-03 −3.22E-03
Nr'linear −9.30E-03 −7.81E-03 −7.56E-03 −6.68E-03