1. Introduction
The barge, which is useful in various fields, has no self-propelled ability compared to other types of ships, and performs maritime transportation while being exposed to difficulties in operation. In the case of a barge drawn from the stern of a tugboat, the stability of the course is very unstable, which is accompanied by a large yaw movement during towing.
The hydrodynamic forces applied to ships on narrow waters, compared to those on open waters, cause the following differences that result in significant impact on the ship’s maneuvering performance. On narrow waters, if the ship deviates from the center of the channel, even though it moves straightforward, the flow between port and starboard becomes unbalanced, causing the ship to be under the lateral force sucked into the closest bank, and the yaw moment where the bow is expelled from the closest bank. Because of the relative location between the ship and the bank, the fluid movement becomes complicated such that unstable hydrodynamic forces result in the yaw movement of the ship in an unintended direction.
Note that, in terms of safety issues, there have been various studies on the bank effect of the hydrodynamic forces on the ship. Norrbin [
1] proposed an empirical equation to predict lateral force and yaw moment, while Li et al. [
2] and Vantorre et al. [
3] examined the ship–bank interaction via model tests. The studies based on empirical equations make it difficult to realize the free surface and viscosity effects; however, analytical methods based on water tank tests lead to considerable time and cost.
Recently, there have been studies that use the computational fluid dynamics (CFD) method to analyze the effects of the free surface, viscous, and rotating propeller on the bank effect. Van Hoydonck et al. [
4] conducted a wide range of water tank tests for KVLCC2. By comparing the results to the CFD solver and potential flow solver, they analyzed various factors that affect the bank effect. Note that Zou and Larsson [
5] analyzed the characteristics of hydrodynamic forces on KVLCC2 because of the difference in the distance from the ship to two different bank forms, i.e., the vertical bank and the slope bank, based on CFD. Ma et al. [
6] performed a simulation to examine the hydrodynamic interaction between the hull, rudder, and bank, because of the change in ship–bank distance and rudder angles. Lo et al. [
7] studied the bank effect and speed effect based on CFD, and confirmed that the closer the ships get to the bank and the faster the speed of the ship becomes, the greater the additional increase in the yaw angle and sway force.
There have been other studies to analyze the course stability of ships close to the bank. Liu et al. [
8] determined the hydrodynamic derivatives near the vertical bank, based on CFD, and used those values to evaluate the directional stability based on the bank distance. Yasukawa [
9] conducted an analysis of pure car carrier (PCC) course stability, by performing the captive model test under different water depths, as well as the distance between ship and bank, drift angle, and heel angle.
Lataire et al. [
10] performed a water tank test on a wide range of bank effects. Based on the test results, they proposed a mathematical model by which the bank effect can be predicted. In particular, they separated the lateral force into the force that works on the forward perpendicular and aft perpendicular, and analyzed the characteristics of the forces separated, based on the difference in ship speed, ship–bank distance, blockage ratio, and propeller action.
In the references Zou and Larsson [
5], Ma et al. [
6], and Liu et al. [
8] that used CFD for bank effect analysis, the effect of free surface was not considered; however, Hoydonck et al. [
4], which did consider the effect of free surface, compared the calculation results at
y′ = 3/4
B, the relatively large distance at intermediate depth such that they reported that the effect of free surface modelling was not significant. As ships will be influenced more by the reflected waves from the bank as they become closer to the bank, it is necessary to determine how this will affect the hydrodynamic forces on the ships.
Because the characteristics of the yaw movement differ depending on the bow shape of the barges [
11], this study will determine the characteristics of hydrodynamic forces that act on barges with different bow shapes, using the simulation results and analyzing the bank effects on the maneuvering stability of barges. To consider the nonlinear effect, this study calculates the surge force, sway force, and yaw moment on barges with different shapes in confined waters, using CFD. Moreover, it analyzes the effects of the flow field near the barges, particularly the effect of the free surface, such as the reflected waves from the bank.
3. Numerical Results and Discussion
3.1. Verification
The Grid Convergence Index (GCI) method, based on the Richardson extrapolation, was used to estimate the grid convergence uncertainty of the numerical simulations in this study.
The verification process using the GCI method was followed according to the method described by Celik et al. [
18]. The expression of the use of GCI method for verification is described with reference to [
19].
The numerical convergence ratio
is calculated as follows:
Here, it is obtained by and , each of which represents the difference in the solution value between medium-fine and coarse-medium, respectively. , , or each represent the solution calculated for the fine, medium, or coarse mesh.
The order of accuracy
can be obtained as follows, using a constant refinement ratio
,
:
The extrapolated values are calculated as shown in Equation (5):
In addition, the approximate relative error and extrapolated relative error are obtained by the following equations, respectively:
Finally, the medium-grid convergence index is calculated as follows:
The number of meshes for the grid convergence study is shown in
Table 2. The surge force (
XH) was calculated using the number of coarse, medium, and fine grids, and the results are shown in
Table 3. As shown in the calculation results, the numerical uncertainty for the medium mesh using the GCI method was estimated to be 1.40%. According to the results of the grid convergence uncertainty study, it was concluded that efficient numerical simulation is possible, using the medium mesh. Therefore, the medium mesh was selected as the base size in this study.
The authors conducted the research in the same way as the CFD method used in this paper for KVLCC2 and DTC ships in 2017, to take into account the shallow water effect [
20]. As a result, it was found that the numerical value was similar to the experimental value. In the medium mesh, KVLCC2 showed an error of 2.99%, and DTC 4.82%. Although this study did not directly compare with the experimental data, it is believed that there would be no significant errors in this study, because the authors confirmed good results from previous research.
3.2. Bank Effects on Hydrodynamic Forces and Moment
Figure 4 and
Figure 5 show the hydrodynamic force changes, based on the difference in the distance from KNU-020 and KNU-030 barge to the bank, respectively. Surge force (
XH), sway force (
YH), and yaw moment (
NH) used the nondimensionalized value shown below, which shows cases with, and without, skeg.
First, as for the results of the surge force, both barges demonstrated that the surge force value corresponding to the total resistance increased, because of the installation of the skeg. KNU-030 with a square-shaped bow demonstrated a higher surge force compared to KNU-020 with the round-shaped bow. KNU-020 showed the largest
X′ at
y′ = 0.2
B when it was closest to the bank. However, the
X′ of KNU-030 decreased, because it was closer to the bank. Zou and Larsson [
5] and Van Hoydonck et al. [
4] studied the bank effect of KVLCC2, showing that
X′ increased as the ships became closer to the bank. The bow of KNU-030 used in this study is different from that of KVLCC2 or KNU-020; thus, the distribution of surge force was shown differently.
Note that the sway force does not seem to be significantly affected by the installation of skeg. For KNU-030, there is a very small difference in sway force if the ship with skeg is closer to the bank. Regardless of the shape of the bow, both barges demonstrated that the suction force, the pulling force toward the bank, increased as they become closer to the bank.
The yaw moment gradually increased in both barges until y′ = 0.5B; however, at y′ = 0.2B, when the barges became closest to the bank, it tended to drastically increase; the bank effect had a significant effect from the distance that corresponded to 1/5 of the barge width. In particular, KNU-030 with skeg demonstrated that yaw moment rapidly increased; thus, it was determined that the bank effect at this condition might have been considerable.
3.3. Analysis of Flow Field and Hydrodynamic Pressure
Next, the study compared the characteristics of the flow field, based on the barge-bank distance and the subsequent hydrodynamic pressure distribution on the barge because of such fluid flow.
Figure 6 shows the (a) wave pattern, (b) velocity field with line integral convolution and hydrodynamic pressure distribution on the bottom, and (c) hydrodynamic pressure on the starboard and port side of KNU-020 barge at
y′ = 0.2
B, when it is closest to the bank. The hydrodynamic pressure is attributed to the subtraction of the hydrostatic pressure from the total pressure. Hydrodynamic pressure in
Figure 6c indicates the pressure distribution of 3 m up from the bottom of the barge. Note that the starboard side is closer to the bank.
Figure 6b,c show the increase in flow speed, and the pressure decrease because of the bank effect in the shoulder of the starboard side of the bow. The flow speed in the port side relatively decreases such that the pressure becomes higher than that of the starboard. However, the area toward the end of the bow from the shoulder demonstrated that the pressure on the starboard side becomes higher as the flow speed in the starboard becomes slower compared to the port side. It is estimated that because of this pressure difference from the shoulder to the end of the bow, the bow-out moment that pushes the bow outward from the bank occurs. The analysis of the pressure distribution from the shoulder to the stern demonstrates that the port side pressure gradually increases more than the starboard side pressure. Because the port side pressure tends to be higher than the starboard side pressure, the suction force that pushes the barge from port to starboard side occurs. The starboard side pressure shows oscillating because of the reflective waves from the bank, which can be shown in
Figure 6a. It is believed that the diverging waves due to the progress of the barge do not move away from the barge, and return to the barge after being reflected from the bank, causing the change in the hydrodynamic force.
Figure 7 shows KNU-020 at
y′ = 0.5
B. The port–starboard pressure difference in the shoulder was not significant, and the pressure in the area toward the end of the bow shows that the starboard side was somewhat higher than the port side. The port side pressure was somewhat higher at the end of the stern; thus, this and the higher starboard pressure from the shoulder to the end of the bow may have caused the bow-out moment.
Overall, the port side pressure is higher from the shoulder to the stern, which causes the suction force on the barge. A careful look at the reflected waves in
Figure 8a shows that the part where the diverging waves hit the quay wall and affects the ship again is from the center of the hull toward the stern. In other words, it is possible to understand the phenomenon in which the hydrodynamic force acting on the hull is changed by the bank effect, even at the current distance from the quay wall (50% of the ship width).
Figure 8 shows the simulation result at
y′ = 1.0
B, i.e., when the barge is distant from the bank by the barge width. This shows that the pressure distribution between the port and starboard sides is almost identical from the bow to about ¼
L of the stern; however, from ¼
L of the stern toward the stern, the port pressure is somewhat higher than the starboard pressure, which may have caused the bow-out moment and suction force.
The analysis of KNU-030 is shown in
Figure 9 at
y′ = 0.2
B, the closest to the bank. As shown in
Figure 9b, the flow speed at the bow shoulder is somewhat faster in the port side, further from the bank. Thus, the starboard pressure is higher than the port side, causing the bow-out moment. This flow speed in the bow is different from that shown in KNU-020 with the round-shaped bow. The stern pressure from the shoulder is higher in the port side than in the starboard side, causing the suction force that pushes the barge to the bank. Similar to KNU-020, the wave pattern shows that the diverging waves from the shoulder do not disperse, and are blocked by the bank and returned to the barge.
Figure 10 shows the calculation result at
y′ = 0.5
B, when the barge-bank distance is half of the barge width. Note that starboard pressure is slightly higher than in the port bow shoulder. The port pressure is slightly higher than the starboard side from the shoulder to the near center of the barge. However, from the center to the stern, the pressure changes between the port and starboard sides, such that the port side pressure becomes a bit higher than the starboard side pressure. As shown in
Figure 10a, such phenomena result in changing the hydrodynamic force on the barge, because the reflective waves from the bank interfere with the hydrodynamic force from the center of the barge.
Finally,
Figure 11 shows the result at
y′ = 1.0
B, where the pressure by fluid flow is very similar between the port and starboard sides. When the barge-bank distance is about the same as the barge width, the bank effect is negligible, such that the sway force and yaw moment on the barge are small.
Figure 12 shows the reflected waves from the bank when KNU-020 and KNU-030 are closest to the bank (
y′ = 0.2
B). The diverging waves from KNU-020 are reflected from the bank after the shoulder, and those from KNU-030 are reflected from the bank from a bit further toward the bow, such that they are reflected toward the barge along the bank. Such reflective waves from the bank do not push the close barge away; however, they cause pressure imbalance because of the difference in the height of free surface at the opposite of the bank, causing the suction force and bow-out moment that pull the ship toward the bank. Indeed, it is clear that the pressure change caused by the difference in the speed of flow coming into the narrow gap between the ship and the bank is an important factor to generating the bank effect. However, it is determined that the effect of free surface must be considered when conducting a comprehensive analysis of the bank effect when the ship is close to the bank.
4. Conclusions
In this study, we conducted a simulation on the change of hydrodynamic forces because of the bank effect and the generation and characteristics of reflected waves due to the bank, by targeting two barges with different bow shapes based on CFD; the results of the analysis are as follows.
X′ value, which represents the surge force, tended to differ, according to the bow shape of the two barges. In other words, KNU-020 with the round-shaped bow showed a higher X′ as the barge became closer to the bank, whereas KNU-030 with the square-shaped bow showed a smaller X′. Y′ value, which represents the sway force, tended to increase as both barges became closer to the bank, because it works as a suction force that pulls them toward the closest bank. N′ value, which represents the yaw moment, increased as the barges became closer to the bank, regardless of the shape of the bow. At y′ = 0.2B, when the barges were at the closest to the bank, it rapidly soared. The moment direction works as the bow-out moment that pushes the bow away from the bank.
As both vessels got closer to the bank, the suction force and bow-out moment were increasing, so caution must be taken to avoid contact with the bank when navigating near the quay wall. In particular, in the case of the KNU-030 barge with the skeg installed, the bow-out moment grew rapidly, so it seems to be necessary to operate with special attention.
The flow speed distribution by the bow shape showed different tendencies between the barges. For KNU-020, the flow speed near the shoulder toward the bank was faster than that of the port side, which is toward the opposite side. It was somewhat slow from the shoulder to the end of the bow; thus, these caused pressure differences in these areas, which resulted in the bow-out moment. However, the flow speed in the port side shoulder became faster, such that the pressure decreased compared to the pressure in the starboard side for KNU-030. This may cause the bow-out moment. Note that, in the port side, the pressure distribution from the shoulder toward the stern was higher than in the starboard side, such that the suction force that pushes the ship towards the bank was identical in both barges. The waves reflected from the bank were imbalanced with the diverging waves from the port side of the barges; thus, it is estimated that such a difference in free surface height may change the pressure on barges.