J Navig Port Res > Volume 46(2); 2022 > Article Vo, Mai, Jeon, and Yoon: Experimental Investigation of the Hydrodynamic Characteristics of a Ship due to Bank Effect

### ABSTRACT

When a ship moves in the proximity of the lateral bank, bank suction forces are generated due to bank effects. Thus, hydrodynamic forces can significantly impact the ship’s maneuverability and navigation safety. In this study, model tests were performed to investigate the hydrodynamic forces exerted on a ship, especially suction forces caused by bank effects, using captive model and bank effect tests. A low-speed condition was selected in this study, because of the perilous situation as the ship moves close to the bank. The accuracy of the hydrodynamic forces exerted on the hull was verified, by comparing the results of the static drift test with the results obtained from other institutes at design speed. The straight simulation caused by bank effects was then implemented using estimated hydrodynamic coefficients.

### 1. Introduction

Nowadays, waterway transport is growing, so the demand for transportation of ships to harbors, canals, rivers, and restricted water areas is increasing. One of the parameters in restricted water is the influence of the lateral bank which significantly affects ship maneuvering. As the ship approaches the lateral bank, displaced water gets squeezed in between the ship’s hull and the lateral bank, since it has nowhere to go. Increase in the water flow causes a drop in pressure, pulling the hull down to fill the vacuum, thereby increasing the hydrodynamic force. Therefore, it is crucial to estimate the hydrodynamic force acting on a ship which is close to the lateral bank.
In the past, several researchers investigated ship maneuvering due to bank effects by means of numerical, experimental techniques. Vantorre et al. (2003) developed an empirical formula for predicting ship-bank interaction force by performing several model ships based on various water depths, lateral distances, forward speeds, and propeller rpm. Sano et al. (2014) and Liu et al. (2018) presented a mathematical model of hull hydrodynamic forces. The hydrodynamic force characteristics in shallow water and bank effects were investigated by conducting an experimental study using the captive model test. Besides, Sano et al. (2014) reported the directional stability in the channel based on the eigenvalue analysis and dynamic maneuvering motion simulation. Moreover, Lataire et al. (2008) evaluated the bank effects induced by irregular bank geometries and formulated a mathematical model. Hoydonck et al. (2019) reported a viscous-flow solver on two different CFD codes to predict the hydrodynamic force on the ship by applying a numerical method to the bank effects. The comparison of CFD results with those of experimental results illustrated the reliable prediction of the ship-bank interaction.
In this study, a model test was performed to investigate the hydrodynamic characteristics acting on a ship, considering the bank effects using the captive model and bank effect tests in deep water. Low speed was selected in this study due to risk conditions. In fact, the ship moved slowly to ensure safe navigation. The accuracy of the experimental results was proved by comparing the hydrodynamic force in the static drift test with other institutes at design speed. On the other hand, the mathematical models of twin-rudders and twin-propellers were presented for this ship. The ship maneuvering characteristics close to the lateral bank were evaluated by straight simulation using estimated hydrodynamic coefficients.

### 2.1 Objective

A scale ratio of 1/11 was used as the target model ship for the captive model and bank effect tests. In fact, this ship was equipped with twin-propellers and twin-rudders. However, within the scope of this study, only the bare hull was used for the experiment. Besides, the mathematical models of twin-propellers and twin-rudders were introduced and applied for maneuvering simulation. The model ship is depicted in Fig. 1, the parameters of the ship are described in Table 1.

### 2.2 Equation of motion

#### 2.2.1 Coordinate system

Fig. 2 describes the two right-hand coordinate systems in the ship. The earth-fixed coordinate system Oxy defines the ship trajectory, and the body-fixed coordinate system oxbyb defines the hydrodynamic force acting on the hull. Three Degrees of Freedom (3-DoF) equation of motion was applied in this study, based on Newton’s second law of motion. Additionally, the MMG (Maneuvering Modeling Group) model introduced by Yasukawa et al (2015) was employed for maneuvering simulation of the ship. Thus, the equation of motion under the coordinate system was expressed as in Eq. (1).
##### (1)
m(u˙-vr-xGr2)=XH+XR+XPm(v˙-ur-xGr˙)=YH+YR+YPIzzr˙+mxG(v˙+ur)=NH+NR+NP
where, m is ship mass; u, v and r denote surge, sway, and yaw velocity, respectively; “.” above the variables denotes the rate of change with respect to time; xG is the longitudinal center of gravity; Izz represents the yaw mass moment of inertial; subscripts H, R, and P mean hull, rudder, and propeller, respectively. The external forces on the right side were divided as the hydrodynamic force, control force, and thrust.

#### 2.2.2 Hull hydrodynamic forces

The hydrodynamic force acting on the ship hull considering bank effect was investigated in this study. Static drift circular motion and combined circular motion with drift test were performed to estimate the linear and non-linear hydrodynamic forces with respect to surge velocity (u), sway velocity (v), and yaw rate (r), respectively. For bank effect, the hydrodynamic forces were expressed in terms of lateral distance from bank to the center of the ship (yBank) as shown in Fig. 3. Vantorre et al. (2003) proposed the variable ηb = 1/yBank for bank effect, which easily described the hydrodynamic force as a third-order function. Sano et al. (2014) and Liu et al. (2018) reported mathematical models of the hydrodynamic forces acting on the hull considering the bank effects. Yηbηb+Yηbηbηbηb3 and Nηbηb+Nηbηbηbηb3 was obtained by bank effect test, Yvvηbv2ηb+Yvηbηbvηb2 and Nvvηbv2ηb+Nvηbηbvηb2 combined bank effect with drift test. The mathematical model of the hydrodynamic force acting on the hull was written as Eq. (2).
##### (2)
XH=Xu˙u˙+Xuuu2+Xvvv2+Xrrr2+XvrvrYH=Yv˙v˙+Yr˙r˙+Yvv+Yv|v|v|v|+Yrr+Yr|r|r|r|+Yvvrv2r+Yvrrvr2+Yηbηb+Yηbηbηbηb3+Yvvηbv2ηb+Yvηbηbvη2bNH=Nv˙v˙+Nr˙r˙+Nvv+Nv|v|v|v|+Nrr+Nr|r|r|r|+Nvvrv2r+Nvrrvr2+Nηbηb+Nηbηbηbηb3+Nvvηbv2ηb+Nvηbηbvη2b
Non-dimensional of the hydrodynamic forces and variables were expressed in Eq. (3), where ρ is water density.
##### (3)
u=uU;v=vU;r=rLU;ηb=ηbB2X=FX0.5ρL2U2;Y=FY0.5ρL2U2;N=N0.5ρL3U2

#### 2.2.3 Propeller model

This ship was equipped with twin-propellers and twin-rudders. Kim, et al.(2021) proposed a mathematical model of twin-propellers and twin-rudder for Korea Autonomous Surface Ship (KASS). The propeller model is presents as follow,
##### (4)
XP=(1-t)ρD4[(nP)2KTPJPP+(nS)2KTSJPS]YP=0NP=yP(1-t)ρD4[(nP)2KTPJPP-(nS)2KTSJPS]
where, D is the propeller diameter; nP,S is the propeller revolution. The propeller rotational directions are opposite to each other based on the shape of propeller; t is the thrust deduction factor; yP is the lateral propeller position from the centerline; subscripts P and S mean port and starboard side, respectively. The thrust coefficient KT defined the function of advance coefficient JP and it can obtain from the POW (Propeller Open Water) test. Fig. 4 presents the experimental results of the propeller open water test performed by KRISO. Then, the function of thrust coefficients is displayed in Eq. (5).
##### (5)
KT=-0.01673JP2-0.47468JP+0.46066
JPP,S=uPP,SnP,SD
uPP,S=(1-wPP,S)u
##### (8)
wPP,S=wP0exp(-CPvP2)
##### (9)
CPP=CP-andCPS=CP+whenβP>0CPP=CP+andCPS=CP-whenβP<0vP=v+xPr
where, wP is the wake coefficient through the maneuvering motion; wP0 is the wake coefficient at propeller position in straight moving; βP is the geometrical inflow angle to propeller in maneuvering motion; CP is wake changing coefficients due to lateral motion, and xP is the longitudinal propeller position from midship.

#### 2.2.4 Rudder model

Effective rudder forces for twin-rudders were expressed as follows:
##### (10)
XR=-(1-tR)(FNP+FNS)sinδYR=(1+aH)(FNP+FNS)cosδNR=(xR+aHxH)(FNP+FNS)cosδ-yR(1-tR)(FNP-FNS)sinδ
The parameters were defined as follows:
##### (11)
FNP,S=12ρAR(URP,S)2fasinαRP,S
##### (12)
URP,S=(uRP,S)2+(vRP,S)2
fa=6.13ΛΛ+2.25
##### (14)
αRP,S=δ-tan-1(vRP,SuRP,S)
##### (15)
uRP,S=ɛuPP,Sη[1-κ(1+8KTP,Sπ(JPP,S)2-1)]2+(1-η)
##### (16)
vRP,S=γRP,S(v+lRr)γRP=γR-andγRS=γR+whenβR>0γRP=γR+andγRS=γR-whenβR<0
where FN is the normal rudder force; tR, aH, and xH mean the steering resistance deduction factor, the rudder increase factor, and the position of an additional lateral force component, respectively; ρ is water density; UR denotes the profile area of the rudder; UR represents the resultant rudder inflow velocity; fa is the rudder lift gradient coefficient; Λ denotes the rudder aspect ratio; αR is the effective inflow angle to rudder; uR and vR are longitudinal and lateral inflow velocity components to the rudder; ε is a ratio of a wake fraction at propeller and rudder position; γ refers to the flow straightening coefficients; and βR represents the effective inflow angle to rudder in maneuvering motion.

### 3.1 Experimental setup

The experiment was performed in the model basin in Changwon National University, Korea. The dimension of the model basin is 20.0 m in length, 14.0 m in breadth, and 1.95 m in depth. The ship was installed in the carriage where the ship was towed together with the support structure. The load-cell with two components for the force in the x-direction and y-direction was utilized to measure the hydrodynamic force acting on the ship hull. Two load-cells were located near the bow (point A) and stern (point B) as shown in Fig. 5. All motions were restrained throughout the experimental process.
Because the sidewall of the tank was far from the limitation of the carriage movement, the lateral bank, which was a straight line, was manufactured using wood plates which were fixed by wires and clamps to ensure sufficient strength (Fig. 6). The standard deviation of the lateral bank was 3 mm. Fig. 7 shows the experimental setup of the ship close to the bank in deep water.

### 3.2 Test condition

As mentioned in the mathematical model of the hull, the hydrodynamic force was obtained from the captive model test and bank effect in deep water. On the other hand, because of the influence of the bank effect, the ship was moved at low speed to ensure safe navigation. Therefore, the captive model test and bank effect were carried out at low speed. Besides, the static drift test was conducted at design speed to allow comparison with other results and thus ensure the experiment’s accuracy. Table 2 describes the test conditions in this study.

### 3.3 Experimental results

The hydrodynamic forces acting on the ship are the main concerns for evaluating hydrodynamic coefficients. Fig. 8 shows results of resistance test for various speed. It is the drag acting in the x direction when the ship goes straight and it was performed to obtain the hydrodynamic coefficient with respect to surge velocity as Xuu. Fig. 9 displays results of static drift test for various drift angle to achieve the hydrodynamic force with respect to different sway velocities, namely, Xvv, Yv, Yv|v|, Nv, and Nv|v|. The sway velocity is the function of the drift angle (β) defined as v =−U sinβ. In addition, the measured results of hydrodynamic forces of static drift test were compared with those performed by Kim et al.(2021) at design speed. A good agreement, which verified the quality of present experiment method, was observed in yaw moment and linear term of sway force. However, the non-linear term of sway force was slightly different, Kim et al.(2021) results were greater than the present results, it may be the difference in the experimental setup, ship’s scale ratio, and appendages attached on the ship. The present experiment fixed the roll and pitch motions throughout the performance. Moreover, the twin-propeller and twin-rudder were included in the Kim et al. (2021) model, this could also account for the different results.
Fig. 10 shows the results of circular motion test conducted at various angular velocities. This test was carried out to obtained the hydrodynamic forces with respect to different angular velocities, such as, Xrr, Yr, Yr|r|, Nr, and Nr|r|. Combined circular motion with drift test is shown in Fig. 11. Coupled hydrodynamic forces relative to sway velocity and angular velocity as Xvr, Yvvr, Yvrr, Nvvr, and Nvrr were considered by subtracting the hydrodynamic forces with respect to sway velocity in static drift test and angular velocity in circular motion test. A good fit was observed for all the tests.
On the other hand, the ship-bank interaction was conducted at various lateral distances. Figs. 12 and 13 show the suction force of the bank effect test and combined bank effect with drift, respectively. As explained above, when the water has nowhere to go, it gets squeezed in the space between ship and lateral bank, causing an increase in the velocity between ship and bank, creating a larger hydrodynamic force as the ship moves closer to the lateral bank. To easily obtain the third-order function of hydrodynamic force with respect to lateral distance to bank, the non-dimensional distance (ηb) was calculated by dividing half of the ship’s breadth by the lateral distance.

### 4. Maneuvering simulation

The hydrodynamic coefficients were estimated by fitting the hydrodynamic force using the least squares method. Table 3 presents the hydrodynamic coefficients used for maneuvering simulation in which damping hydrodynamic coefficients obtained from the present experiments, added mass coefficients as X, Y, Y, N, and N obtained from pure sway test and pure yaw test by Kim, et al(2021), and interaction coefficients as (1 − tR), (1 + aH), ε, κ, γR, CP obtained from the static rudder, combined static rudder with drift test by Kim, et al.(2021).
In order to validate the hydrodynamic coefficients obtained by the experiments, the maneuvering simulation was carried out and compared with the simulation results by Kim et al. (2021). The turning and zig-zag maneuvers did not include the terms of hydrodynamic coefficients generated by the bank effects. Fig. 14 shows the comparison of the trajectory, ratio of resultant speed to approach ship speed, and yaw rate in turning maneuver. Fig 15 shows the comparison of the trajectory, heading angle, rudder deflection angle, and yaw rate in zig-zag 10°/10° maneuver. The present simulation results give a good agreement with the results obtained by Kim et al. (2021). It confirmed the present hydrodynamic coefficients were reasonable to evaluate the ship’s maneuverability.
Fig. 16 shows a straight simulation when the ship moved in a straight direction close to the lateral bank by applying estimated hydrodynamic coefficients to the equation of motion. As shown in the figure, the ship did not move in a straight direction and tended to move far away from the bank. This was caused by a suction force when the velocity increases at the space between bank and ship which leads to a difference in the pressure on the each of the ship sides. Moreover, the increased pressure around the ship’s bow created a yaw moment pushed the ship away from the bank.
Additionally, the turning rate was reduced as the ship moved far from the bank and the ship could keep a straight path at a lateral distance far enough from the lateral bank. This is demonstrated in Figs. 17-18 which presents the yaw rate and heading angle throughout the straight simulation, it is easy to understand when hydrodynamic force becomes negligible as the ship is far from the bank enough. As the ship is quite close to the bank (yBank = 1.0B), the yaw rate increases dramatically and rapidly, of course, it also quickly decreases even smaller than other cases at x/L=12 where the ship reaches a large enough distance to be unaffected by the lateral bank. Because the ship position of yBank = 1.0B reaches y/L=1.5 at x/L = 12, while the ship position of yBank = 1.5B is being y/L=1.0. The yaw rate at yBank = 1.5B and yBank = 2.0B also smaller than yaw rate at yBank = 2.5B at a position, however, this process is slower, it can occur at x/L=15 even father.

### 5. Conclusion

In this paper, an experimental study was performed to investigate the hydrodynamic characteristics with respect to the bank effect. The hull mathematical model considering the bank effect was presented. Moreover, the mathematical models of twin-rudders and twin-propellers were given for this ship. The captive model and bank effect tests were implemented to obtain the hydrodynamic force with respect to variables relative to 3-DoF motion. The accuracy of the experiment method was verified by comparing the results of the static drift test with results from other institutes at design speed. A good agreement was observed without the non-linear term of sway force. On the other hand, the bank effect test illustrated that the hydrodynamic force of the ship increased as the ship moved close to the lateral bank, even though the ship moved in a straight direction. The hydrodynamic coefficients were obtained for 3-DoF motion by fitting the hydrodynamic force using the least square method. The straight simulation was presented when the ship moved close to the lateral bank by applying estimated hydrodynamic coefficients to the equation of motion. The ship did not move in a straight direction and tended to move far away from the bank as the ship moved close to the bank and kept a straight path at a lateral distance far enough from the lateral bank. The straight maneuver was performed to assess the influence of the lateral bank on maneuverability. Thereby, a control system by rudder deflection, as well as, the lateral distance will be calculated and reasonably designed in future work to keep straight maneuver when the ship moves close to the bank.

### Acknowledgments

This research was 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).
##### Fig. 1.
Model ship ##### Fig. 2.
Coordinate system ##### Fig. 3.
Vertical bank configuration ##### Fig. 4.
Thrust coefficient of POW test (KRISO) ##### Fig. 5.
Experiment installation ##### Fig. 6.
Lateral bank installation ##### Fig. 7.
Ship close to the bank ##### Fig. 8.
Results of resistance test ##### Fig. 9.
Results of static drift test ##### Fig. 10.
Results of circular motion test ##### Fig. 11.
Results of combined circular motion with dirft test ##### Fig. 12.
Results of bank effect test ##### Fig. 13.
Results of combined bank effect with drift test ##### Fig. 14.
Comparison of 35° starboard turning ##### Fig. 15.
Comparison of 10-10 zigzag maneuver ##### Fig. 16.
Straight running simulation nearby bank ##### Fig. 17.
Yaw rate at several Bank distances ##### Fig. 18.
Heading angle at several Bank distances ##### Table 1.
Parameters of ship
Item Full Model
Scale ratio 1 1:11
Length perpendicular, LPP (m) 22.000 2.000
Depth, D (m) 4.140 0.376
Draft, T (m) 1.250 0.114
Displacement volume, ∇ (m3) 86.681 0.064
Rudder area, AR (m2) 0.518 0.004
Rudder span length, HR (m) 0.900 0.082
Propeller diameter, Dp (m) 0.950 0.086
##### Table 2.
Test condition
Test Motion variable
Resistance U = 1, 2, 3, 4, 5, 6 (knots)
Static drift U = 3, 6 (knots)
β = 0, ±2, ±4, ±6, ±8, ±10, ±12, ±15, ±18 (deg.)
Circular motion U = 3 (knots)
r′ = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6 (-)
Combined circular motion with drift U = 3 (knots)
r′ = 0.2, 0.3, 0.4, 0.5 (-)
β = 0, ±3, ±6, 9, 12 (deg.)
Bank effect U = 3 (knots)
yBank = 0.9B, 1.1B, 1.3B, 1.5B, 1.7B, 2.0B, 2.5B (m)
Combined bank effect with drift U = 3 (knots)
yBank = 1.0B, 1.5B, 2.0B, 2.5B (m)
β = 0, ±2, ±4, ±6, ±8, ±10 (deg.)
##### Table 3.
Hydrodynamic coefficient
 Hull (×105) Xu̇ −81 Yv̇ −1034 Nv̇ 64 Xuu −627 Yṙ −126 Nṙ −33 Xvv −407 Yv −1444 Nv −278 Xrr 675 Yv|v| −1688 Nv|v| −473 Xvr 226 Yr 558 Nr −376 Yr|r| −586 Nr|r| −296 Yvvr 2273 Nvvr −1584 Yvrr −6897 Nvrr −676 Yηb −51 Nηb 11 Yηbηbηb −228 Nηbηbηb 21 Yvvηb −51 Nvvηb 46 Yvηbηb −45 Nvηb ηb −19 Propeller & rudder 1−tR 0.934 1+aH 0.702 γR+ 0.342 ε 0.960 CP+ −2.713 γR- 0.634 κ 0.695 CP- 11.211

### References

1. Kim, DJ, Choi, H, Kim, YG, Yeo, DJ.(2021), “Mathematical model for harbour manoeuvres of Korea Autonomous Surface Ship (KASS) based on captive model tests”, Conference of KAOSTS.
2. Khanfir, S, Hasegawa, K, Nagarajan, V, Shouji, K, Lee, SK.(2011), “Manoeuvring characteristics of twin-rudder systems:” Rudder-hull interaction effect on the manoeuvrability of twin-rudder ships”, Journal of Marine Science and Technology, Vol. 16, No. 4, pp. 472-490.  3. Liu, H, Ma, N, Gu, X.(2018), “Experimental Study on Ship-Bank Interaction of Very Large Crude Carrier in Shallow Water”, Journal of Shanghai Jiaotong University (Science), Vol. 23, No. 6, pp. 730-739.  4. Lataire, E, Vantorre, M.(2008), “Ship-bank interaction induced by irregular bank geometries”, Proceedings of the 27th Symposium on Naval Hydrodynamics, October 5-10 (2008), pp. 1-13.
5. Sukas, OF, Kinaci, OK, Bal, S.(2019), “System -based prediction of maneuvering performance of twin-propeller and twin-rudder ship using a modular mathematical model”, Applied Ocean Research, Vol. 84:pp. 145-162. 6. Sano, M, Yasukawa, H, Hata, H.(2014), “Directional stability of a ship in close proximity to channel wall”, Journal of Marine Science and Technology, Vol. 19, No. 4, pp. 376-393.  7. Vantorre, M, Delefortrie, G, Eloot, K, Laforce, E. (2003), Experimental investigation of ship-bank interaction forces, MARSIM 2003, 25th-28th. August;(2003), Kanazawa.
8. Van Hoydonck, W, Toxopeus, S, Eloot, K, Bhawsinka, K, Queutey, P, Visonneau, M.(2019), “Bank effects for KVLCC2”, Journal of Marine Science and Technology, Vol. 24, No. 1, pp. 174-199.  9. Yasukawa, H, Yoshimura, Y.(2015), “Introduction of MMG standard method for ship maneuvering predictions”, Journal of Marine Science and Technology, Vol. 20, No. 1, pp. 37-52.  Editorial Office
C1-327 Korea Maritime and Ocean University
727 Taejong-ro, Youngdo-gu, Busan 49112, Korea
Tel: +82-51-410-4127    Fax: +82-51-404-5993    E-mail: jkinpr@kmou.ac.kr     