Computational Analysis for Estimation of Mooring Force Acting on Various Ships in Busan New Port

: Recently, smart port systems connected to autonomous ships have attracted increasing interest. Thus, an intelligent port system can automatically berth to create an intelligent port system. To ensure the safety of large ships moored in open coastal terminals, which are often subjected to bad weather such as strong winds, it is necessary to calculate and evaluate their mooring security on a case-by-case basis. In this study, the mooring capacities of the large ships were estimated according to the port and ﬁshing port design criteria of the Ministry of Ocean and Fisheries. Under the wind speed of 14 m/s, the longitudinal and lateral forces acting on the JBC, KCS, and KVLCC ships are 41.2 and 340 kN, 38.7 and 837 kN, and 77.2 and 222 kN, while under the wind speed of 30 m/s, they are 43 and 1674 kN, 132.7 and 4118 kN, and 159.2 and 1091 kN, respectively, for the mooring forces.


Introduction
With the rapid development of the global shipping business, there are increasing ports and terminals, and ships tend to be larger; meanwhile, with the improvement of loading and unloading efficiency, frequent mooring operations have become the norm in the shipping industry. The operation of mooring cables by crew and shore workers is the riskiest task in port operations due to hydro-meteorological conditions and frequent changes in the ship draft caused by loading and unloading operations. The statistics regarding mooring accidents and injuries in the maritime industry are concerning and emphasize the need for improved mooring methods.
According to the European Harbor Masters' Committee (EHMC) [1], 95% of registered mooring injuries are caused by ropes and wires, with 60% occurring during mooring operations. Reports indicate that losses from major accidents exceeded USD 34 million between 1995 and 2016 [2]. These accidents primarily occur during the handling of ropes and wires. Rope and wire breakages account for 53% of the accident rate, while rope and wire slipping or jumping from the end of the drum or bite contribute to 42% of the accident rate. Equipment failure is responsible for 5.0% of these accidents (Figure 1) [3][4][5]. The UK P&I Club states that approximately 1/10 of personal injuries in ports are caused by mooring failures, with separate ropes and wire ropes being the primary cause [6][7][8][9]. The Australian Maritime Safety Authority (AMSA) received 227 mooring-related incident reports between 2011 and 2015, resulting in 51 injuries [10]. The International Maritime Health Association collected 331 mooring reports from 2006 to 2016 [11]. Data from the Transportation Safety Board of Canada (TSB) database reveal that from 2007 to 2017, 24 mooring incidents were reported involving both Canadian and foreign-flag ships, leading to twenty-four injuries, including two serious injuries. Over a seven-year period from 2001 to 2017, other member states of the International Maritime Organization (IMO) investigated 25 marine incidents and 14 very serious marine incidents related to mooring. These incidents resulted in a total of 27 injuries and 13 serious injuries, caused by the buckling or trapping of mooring lines, Figure 1. Types of incidents resulting in personal injury from 1995 to 2016 (adapted with permission from Ref. [15]).
These statistics underscore the importance of addressing the challenges associated with traditional mooring methods and developing improved techniques in the field of ocean engineering to enhance safety and efficiency in mooring operations. The research on mooring systems has evolved from single-cable monitoring to system monitoring, with that on automatic mooring systems in port terminals being widely emphasized and gradually applied in recent years, and automation technology can be used to ensure the safety of ship operations during mooring instead of the traditional cable method. This area is represented by the results developed by European companies, such as the Shore Tension ® and MoorMaster TM mooring systems, both of which are already in use as automated mooring solutions ( Figure 2). According to the mooring arrangement, automated mooring branched into two types, with and without the cable [16,17]. Jan et al. [18] proposed a shore-based, ship-independent mooring aid that supports master and pilot maneuvering in tight harbor pools and inclement weather by measuring distance and speed to determine precise mooring positions through reference points aligned with dock metrics, thus meeting the challenges of damage that can be caused by too-rapid approaches and steep angles of attack. Chen et al. [19] analyzed and calculated the ship's heading and normal vector, determined the bow and stern position using feature points, obtained the line segment between the ship and the pier with the help of the area growth method, and applied a visibility analysis to identify the bow and stern. Using 3D LiDAR data, the distance, speed, and approach angle between the bow and stern of the dynamic ship and the quay were analyzed qualitatively and quantitatively, and the berthing situation of the Ro-Ro ship "Ocean Island" in Lushun port was effectively monitored in field experiments, which fully demonstrates the feasibility and effectiveness of this method in the identification of dynamic ships' targets and safe berthing. Çağatay et al. [20] analyzed the operational strategies, technology applications, renewable energy and alternative fuels, and energy management systems through a systematic literature review to indicate the research directions and research gaps that exist in order to improve the energy efficiency and environmental performance of the ports and terminals, revealing the existence of a great potential for the ports in the field of energy efficiency and providing influential research opportunities for These statistics underscore the importance of addressing the challenges associated with traditional mooring methods and developing improved techniques in the field of ocean engineering to enhance safety and efficiency in mooring operations. The research on mooring systems has evolved from single-cable monitoring to system monitoring, with that on automatic mooring systems in port terminals being widely emphasized and gradually applied in recent years, and automation technology can be used to ensure the safety of ship operations during mooring instead of the traditional cable method. This area is represented by the results developed by European companies, such as the Shore Tension ® and MoorMaster TM mooring systems, both of which are already in use as automated mooring solutions ( Figure 2). According to the mooring arrangement, automated mooring branched into two types, with and without the cable [16,17]. Jan et al. [18] proposed a shore-based, ship-independent mooring aid that supports master and pilot maneuvering in tight harbor pools and inclement weather by measuring distance and speed to determine precise mooring positions through reference points aligned with dock metrics, thus meeting the challenges of damage that can be caused by too-rapid approaches and steep angles of attack. Chen et al. [19] analyzed and calculated the ship's heading and normal vector, determined the bow and stern position using feature points, obtained the line segment between the ship and the pier with the help of the area growth method, and applied a visibility analysis to identify the bow and stern. Using 3D LiDAR data, the distance, speed, and approach angle between the bow and stern of the dynamic ship and the quay were analyzed qualitatively and quantitatively, and the berthing situation of the Ro-Ro ship "Ocean Island" in Lushun port was effectively monitored in field experiments, which fully demonstrates the feasibility and effectiveness of this method in the identification of dynamic ships' targets and safe berthing. Çagatay et al. [20] analyzed the operational strategies, technology applications, renewable energy and alternative fuels, and energy management systems through a systematic literature review to indicate the research directions and research gaps that exist in order to improve the energy efficiency and environmental performance of the ports and terminals, revealing the existence of a great potential for the ports in the field of energy efficiency and providing influential research opportunities for the researchers. Çagatay et al. [21] showed that smart grids can reduce costs and that deploying energy storage systems in port microgrids will result in important cost savings. Overall, ports utilizing renewable energy can realize significant cost savings. Kuzu et al. [3] conducted a comparative analysis of vacuum-based automated, magnetic, and conventional mooring systems involving ropes and windlass, considering operation safety, operating cost, maintenance cost, environmental impact, ease of handling, and system limitations. Their findings suggest that the vacuum-based automated mooring system is more suitable for future applications. Çagatay et al. [22] developed a recoverable robust optimization approach to enhance the performance of container terminals by considering uncertainty parameters in berth and quay crane planning and minimizing costs due to resource efficiency. This study provides a feasible theoretical basis for the future installation of mooring equipment.
ploying energy storage systems in port microgrids will result in important cost savings. Overall, ports utilizing renewable energy can realize significant cost savings. Kuzu et al. [3] conducted a comparative analysis of vacuum-based automated, magnetic, and conventional mooring systems involving ropes and windlass, considering operation safety, operating cost, maintenance cost, environmental impact, ease of handling, and system limitations. Their findings suggest that the vacuum-based automated mooring system is more suitable for future applications. Çağatay et al. [22] developed a recoverable robust optimization approach to enhance the performance of container terminals by considering uncertainty parameters in berth and quay crane planning and minimizing costs due to resource efficiency. This study provides a feasible theoretical basis for the future installation of mooring equipment.
The automated mooring systems in ports hold a significant role regarding operational reliability and safety enhancement, the minimization of human involvement during berthing, and the elimination of potential errors and risks. Additionally, automated mooring systems improve port operation efficiency and reduce downtime. These systems shorten the mooring time and streamline berth operation procedures, leading to shorter transit times and enhanced terminal operation efficiency and throughput.
(a) Shore Tension ® mooring systems (b) MoorMaster TM mooring systems The intricate interactions among the ship, dock, and water pose significant challenges due to wind, wave, and current disturbances. Achieving precise control in automated mooring requires a comprehensive understanding of the forces acting on the ship during the process. This study aims to investigate the primary forces involved, including wind, wave, and current forces, considering both the transverse and longitudinal directions. By quantifying these forces under various environmental conditions, valuable insights can be gained to improve the design and optimization of automated mooring systems, ensuring their reliability and effectiveness in challenging operational scenarios.
Andersen [23] studied the influence of the cargo-hold configuration on wind forces on the deck of a 9000 TEU container ship using wind tunnel tests with a 1:450 scale model. Longitudinal and transverse forces, together with yawing moments, were measured. Wang et al. [24] investigated the wind and wave loads regarding the motion characteristics of a 1900 TEU container ship through wind tunnel tests, towing tank experiments, and computational fluid dynamics (CFD) simulations. The results indicated that the combined wind and wave loads resulted in higher resistance than the sum of independent wind and wave loads. Ricci et al. [25] conducted three-dimensional steady-state Reynolds-averaged Navier-Stokes (RANS) simulations to study the wind forces acting on a large cruise ship docked at the Rotterdam cruise terminal. The study revealed that tall buildings nearby may increase the local surface pressure on the ship, reducing horizontal forces but increasing vertical uplift forces. In smaller ships, tall buildings nearby may increase both horizontal and vertical forces. Zhou et al. [26] investigated the impact of wind and currents on the dynamic behavior of ships in ports and waterways by analyzing automatic identification system data, meteorological data, and hydrological data from the Port of Rotterdam. The results demonstrated variations in the influence of wind and currents on the ground The automated mooring systems in ports hold a significant role regarding operational reliability and safety enhancement, the minimization of human involvement during berthing, and the elimination of potential errors and risks. Additionally, automated mooring systems improve port operation efficiency and reduce downtime. These systems shorten the mooring time and streamline berth operation procedures, leading to shorter transit times and enhanced terminal operation efficiency and throughput.
The intricate interactions among the ship, dock, and water pose significant challenges due to wind, wave, and current disturbances. Achieving precise control in automated mooring requires a comprehensive understanding of the forces acting on the ship during the process. This study aims to investigate the primary forces involved, including wind, wave, and current forces, considering both the transverse and longitudinal directions. By quantifying these forces under various environmental conditions, valuable insights can be gained to improve the design and optimization of automated mooring systems, ensuring their reliability and effectiveness in challenging operational scenarios.
Andersen [23] studied the influence of the cargo-hold configuration on wind forces on the deck of a 9000 TEU container ship using wind tunnel tests with a 1:450 scale model. Longitudinal and transverse forces, together with yawing moments, were measured. Wang et al. [24] investigated the wind and wave loads regarding the motion characteristics of a 1900 TEU container ship through wind tunnel tests, towing tank experiments, and computational fluid dynamics (CFD) simulations. The results indicated that the combined wind and wave loads resulted in higher resistance than the sum of independent wind and wave loads. Ricci et al. [25] conducted three-dimensional steady-state Reynolds-averaged Navier-Stokes (RANS) simulations to study the wind forces acting on a large cruise ship docked at the Rotterdam cruise terminal. The study revealed that tall buildings nearby may increase the local surface pressure on the ship, reducing horizontal forces but increasing vertical uplift forces. In smaller ships, tall buildings nearby may increase both horizontal and vertical forces. Zhou et al. [26] investigated the impact of wind and currents on the dynamic behavior of ships in ports and waterways by analyzing automatic identification system data, meteorological data, and hydrological data from the Port of Rotterdam. The results demonstrated variations in the influence of wind and currents on the ground speed, yaw angle, and drift angle for different-sized ships. Kobayashi et al. [27] performed a grid sensitivity analysis on the superstructure of a ship using their in-house solver and overset grid technology with the Japan Bulk Carrier (JBC) model. The study revealed that steady-state CFD analysis reliably estimates the wind forces and moments on the superstructure of different ship models in wind tunnel tests, showing good agreement with measurement data and outperforming regression formulas in certain cases. Martić et al. [28] conducted a study on the influence of the prismatic coefficient, the longitudinal position of the center of buoyancy, the trim, the pitch radius of gyration, and the ship speed on the added resistance in waves for the KCS ship under regular head waves and different sea conditions using the 3D panel method based on the Kelvin-type Green function. The research findings provide insights into how variations in ship characteristics affect the added resistance in waves. Sasa et al. [29] proposed a novel algorithm to simulate the sliding anchor and grounding situations under wave and wind conditions, comparing the results of dynamically computed anchor chain forces with those from the catenary method. Their aimed to provide valuable insights for improving port planning related to anchorage. Lee et al. [30] conducted a failure analysis on a ship that experienced anchor failure in a port, revealing that wave drift forces could lead to the failure of mooring systems. They suggested new mooring system designs considering wave loads for port planning. Lin et al. [31] conducted a resistance and propulsion performance analysis on a 25 m ship model in actual wind, wave, and current flow environments. Additionally, corresponding correction methods were developed to ascertain the performance under actual sea conditions. The study also evaluated the impact of an energy-saving device (ESD) on the ship's performance in actual sea conditions. Farkas et al. [32] investigated the resistance and propulsion characteristics of a handy-size bulk carrier through full-scale numerical and experimental evaluations. They detailed the analysis result, verified the numerical results, and compared the different turbulence models provided, along with the pros and cons of the extrapolation methods discussion. The study emphasized the significance of viscous flow based CFD in enhancing the reliability of extrapolation methods and reducing the increment of the resistance coefficient. It also highlighted the advantages of assessing ship hydrodynamic characteristics using full-scale evaluations.
Many scholars have previously studied the mooring forces of ships; however, the combined effects of wind, waves, and currents have not been considered. This research aimed to calculate the mooring forces on three types of ships commonly operated in ports: cargo ships, container ships, and oil tankers. The results obtained from this study can be adopted to estimate the mooring forces for these three types of ships, thus enhancing the safety of ships under mooring conditions. Numerical simulations for wind resistance were performed using the unsteady Reynolds-averaged Navier-Stokes solver in STAR CCM+. The simulations investigated the variation of resistance on ships under different wind directions and uniform wind velocities while considering the design draught of the ship. In addition, under moored conditions, the port and fishing port design criteria were utilized to calculate the maximum current force acting on the ship within the harbor. Furthermore, ANSYS AQWA software was employed to calculate the maximum wave force on the ship during a three-hour period under irregular wave conditions. The combination of these calculations provided a more accurate estimation of the mooring forces experienced by ships in moored conditions.

Target Ship
The JBC (Japan Bulk Carrier) is a capsize bulk carrier equipped with a stern duct as an ESD. The National Maritime Research Institute (NMRI), Yokohama National University, and the Ship Building Research Centre of Japan (SRC) were jointly involved in the design of a ship hull, a duct, and a rudder. Towing tank experiments were planned at NMRI, SRC, and Osaka University, which included resistance tests, self-propulsion tests, and PIV measurements of stern flow fields. The KCS was conceived to provide data for the explication of both the flow physics and CFD validation for a modern container ship with a bulbous bow. The Korea Research Institute for Ships and Ocean Engineering (KRISO) performed towing tank experiments to obtain the resistance, mean flow data, and free surface waves. The KRISO 300 K Very Large Crude Carrier (KVLCC) is a specific type of VLCC ship developed by the KRISO designed for large volumes of crude oil across long distances ( Figure 3, Table 1).

Ocean Environmental Condition
The target ship will be moored at the Busan New Port, Republic of Korea, as shown in Figure 4. The mooring force acting on the ship in the port is affected by various environmental disturbances. In general, ports are designed to minimize disturbances to the environment. However, wind loads are inevitably generated. To reduce the loads caused by waves, tidal barriers and breakwaters (including offshore breakwaters) are installed around the harbor, resulting in the size of the waves being proportional to the height of the sea level. Thus, the wave size in the harbor decreased significantly. The JBC, KCS, and

Ocean Environmental Condition
The target ship will be moored at the Busan New Port, Republic of Korea, as shown in Figure 4. The mooring force acting on the ship in the port is affected by various environmental disturbances. In general, ports are designed to minimize disturbances to the environment. However, wind loads are inevitably generated. To reduce the loads caused by waves, tidal barriers and breakwaters (including offshore breakwaters) are installed around the harbor, resulting in the size of the waves being proportional to the height of the sea level. Thus, the wave size in the harbor decreased significantly. The JBC, KCS, and KVLCC ships are large ships that can also be moored in the harbor when a typhoon occurs. The wind speeds of 14 and 30 m/s were considered in the mooring forces calculation of the ship's moored condition. In Korea, a work alert is issued when the wind speed reaches 14 m/s at sea, and during a typhoon, the wind speed can reach 30 m/s. In addition, the water depth at Busan New Port is approximately 17 m, and the maximum tidal current of 1 m/s is listed in the available information [34]. As shown in Figure 4, on the left and right sides of the outer periphery of the harbor, breakwaters have been installed to ensure the smoothness of the water surface within the port. Due to the ship's angle with the port when it is moored, the current velocity and direction were set at 2 knots and 30 degrees, respectively ( Figure 5).
KVLCC ships are large ships that can also be moored in the harbor when a typhoon occurs. The wind speeds of 14 and 30 m/s were considered in the mooring forces calculation of the ship's moored condition. In Korea, a work alert is issued when the wind speed reaches 14 m/s at sea, and during a typhoon, the wind speed can reach 30 m/s. In addition, the water depth at Busan New Port is approximately 17 m, and the maximum tidal current of 1 m/s is listed in the available information [34]. As shown in Figure 4, on the left and right sides of the outer periphery of the harbor, breakwaters have been installed to ensure the smoothness of the water surface within the port. Due to the ship's angle with the port when it is moored, the current velocity and direction were set at 2 knots and 30 degrees, respectively ( Figure 5).  The wave height is 1 m where the wave direction is established at 30 degrees, according to the port and fishing port design criteria and the information on Busan New Port [34]. Due to the breakwater at the harbor entrance, there are no long waves within the port, and the wave height is not very high. Based on the data from the Busan New Port, the wave period is 3.5 s, and the wave height is 0.5 m.

Viscous Fluid Solver
In the commercial software package STAR-CCM+, a mathematical model based on the Reynolds-averaged Navier-Stokes (RANS) equations is used. These equations are derived from the laws of conservation of mass and momentum. For the incompressible, nonstationary case, the equations are as follows: KVLCC ships are large ships that can also be moored in the harbor when a typhoon occurs. The wind speeds of 14 and 30 m/s were considered in the mooring forces calculation of the ship's moored condition. In Korea, a work alert is issued when the wind speed reaches 14 m/s at sea, and during a typhoon, the wind speed can reach 30 m/s. In addition, the water depth at Busan New Port is approximately 17 m, and the maximum tidal current of 1 m/s is listed in the available information [34]. As shown in Figure 4, on the left and right sides of the outer periphery of the harbor, breakwaters have been installed to ensure the smoothness of the water surface within the port. Due to the ship's angle with the port when it is moored, the current velocity and direction were set at 2 knots and 30 degrees, respectively ( Figure 5).  The wave height is 1 m where the wave direction is established at 30 degrees, according to the port and fishing port design criteria and the information on Busan New Port [34]. Due to the breakwater at the harbor entrance, there are no long waves within the port, and the wave height is not very high. Based on the data from the Busan New Port, the wave period is 3.5 s, and the wave height is 0.5 m.

Viscous Fluid Solver
In the commercial software package STAR-CCM+, a mathematical model based on the Reynolds-averaged Navier-Stokes (RANS) equations is used. These equations are derived from the laws of conservation of mass and momentum. For the incompressible, nonstationary case, the equations are as follows: The wave height is 1 m where the wave direction is established at 30 degrees, according to the port and fishing port design criteria and the information on Busan New Port [34]. Due to the breakwater at the harbor entrance, there are no long waves within the port, and the wave height is not very high. Based on the data from the Busan New Port, the wave period is 3.5 s, and the wave height is 0.5 m.

Viscous Fluid Solver
In the commercial software package STAR-CCM+, a mathematical model based on the Reynolds-averaged Navier-Stokes (RANS) equations is used. These equations are derived from the laws of conservation of mass and momentum. For the incompressible, non-stationary case, the equations are as follows: where µ is the viscosity modulus of the fluid, P is the static pressure, −ρu i u j is the Reynold stress, and f i is the force of the object per unit volume. The RANS turbulence model is based on the modeling of the Reynolds stress transport equation (RST). The stress tensor is modeled as a function of the mean flow. The model used in this study belongs to the class of eddy-viscosity models, which use the concept of the turbulent eddy-viscosity coefficient µ t to model the stress tensor. According to the Boussinesq hypothesis, the RST is modeled as follows [35]: where u t and k represent turbulent viscosity and turbulent kinetic energy, respectively, while the symbol δ ij denotes the Kronecker delta. All the turbulence models employed in the present study are based on the Boussinesq hypothesis. The realizable k-ε (RKE) turbulence model is a two-equation model that solves the transport equations for the turbulent kinetic energy k and turbulent dissipation rate ε. The eddy viscosity in this model is described by the following equation [35]: where C µ is a model coefficient. The k and ε transport equations are as follows: where σ k and σ ε are turbulent Schmidt numbers; f c is the curvature correction factor, G k is the turbulence generation term; G b is the buoyancy generation term; ε 0 is the ambient turbulence value that counteracts the turbulence attenuation; γ M is the expansion dissipation; Cε 1 , Cε 2 , and Cε 3 are the model correlation coefficients; S is the modulus of the mean strain tensor; and S k and Sε are the user-defined source terms.

Calculation of the Current Force
According to fluid dynamics, fluid pressure consists of two parts: frictional resistance and pressure resistance. Frictional resistance is predominant in the bow direction of the ship. For lateral flow, pressure drag is predominant. However, it is very difficult to strictly distinguish between frictional and pressure resistance and perform calculations. According to [36] (7), the fluid pressure caused by the tidal flow can be calculated as Rf : fluid Pressure (kN). P w : specific gravity of seawater (1.03 ton/m 3 ). T: temperature (15 • C). S: flooded surface area (m 2 ). V: flow velocity (m/s). λ: coefficients (0.14).

Calculation of the Wave Force
For the hydrodynamic analysis of a ship under regular waves based on three-dimensional potential flow theory, the following conditions must be satisfied: the water is an ideal and incompressible fluid with no rotation, and the waves are small-amplitude waves. The velocity potential function should satisfy the following conditions [37]: Control equation: Free-surface boundary condition: Seabed boundary condition: Free-surface motion condition after linearization: Free-surface dynamic condition after linearization: In general, the above conditions are solved using the Green function method, which establishes the integral equation between the velocity potential and the Green function. When solving the integral equation, the boundary element method is employed to discretize the object surface into a certain number of units. Subsequently, the strength at the center of each unit (referred to as the source point) is calculated under the boundary conditions, and the velocity potential at each node is obtained based on the pulse source distribution.
After calculating the radiation and diffraction velocity potential of the structure, the first-order hydrodynamic pressure can be expressed as [37]: The total wave force is obtained by integrating over the wet surface of the object: where S is the wet surface of the object; in general, the first-order wave force consists of three components: where F is the first wave force; F HS is the rigid hydrostatic restoring force; F R is the radiation force; and F EX is the first wave excitation force and moment, generated in the interaction of wave incidence and diffraction velocity potential.
To consider the ship's sailing characteristics, it is necessary to select the wave spectrum suitable for the sea area. The common wave spectra used to evaluate the ship's sailing performance in irregular waves are the Pearson Moskowitz (P-M) spectrum and the Bretschneider JONSWAP spectrum. In this study, the standard spectrum representing a sufficiently open sea area, i.e., the P-M spectrum, was used according to the load analysis guidelines of the Korean Register of Shipping (KR) and Chapter 5 of Part 3 of the Ship Rules [38].
The P-M spectrum is given by Equation (16) [39] with the peak frequency. Figure 6 shows the relationship between the wave spectrum and wave frequency for a significant wave height of 0.5 m.
The following relationship exists between T Z , T 1 , and T 0 (17): where T 1 is the mean wave period and T 0 is the peak period.
sufficiently open sea area, i.e., the P-M spectrum, was used according to the load guidelines of the Korean Register of Shipping (KR) and Chapter 5 of Part 3 of Rules [38]. The P-M spectrum is given by Equation (16) [39] with the peak frequency. shows the relationship between the wave spectrum and wave frequency for a sig wave height of 0.5 m. where T1 is the mean wave period and T0 is the peak period.

Grid Validation and Verification
Considering the variability in the wind speed and the influence of the hull sectional shape on the aerodynamic coefficients, it is crucial to calculate the win acting on a streamlined ship. The wind forces acting on the ship branch into long forces (X direction) and lateral forces (Y direction) (Figure 7). The force exerted wind varies significantly based on the upper water surface area. Specifically, wi the bow direction affects the frontal projection area of the upper water surface, r in relatively small longitudinal forces. Conversely, wind blowing from the later tion affects the lateral projection area of the water surface, leading to considerab lateral forces. In this study, a uniform wind flow assumption was proposed t

Grid Validation and Verification
Considering the variability in the wind speed and the influence of the hull's crosssectional shape on the aerodynamic coefficients, it is crucial to calculate the wind loads acting on a streamlined ship. The wind forces acting on the ship branch into longitudinal forces (X direction) and lateral forces (Y direction) (Figure 7). The force exerted by the wind varies significantly based on the upper water surface area. Specifically, wind from the bow direction affects the frontal projection area of the upper water surface, resulting in relatively small longitudinal forces. Conversely, wind blowing from the lateral direction affects the lateral projection area of the water surface, leading to considerably larger lateral forces. In this study, a uniform wind flow assumption was proposed to obtain relatively conservative results, and calculations were accomplished using equations specified in the "Port and Fishing Port Design Criteria" [36].
C X : wind drag coefficient in the X direction.     In this study, a CFD approach was purposed to obtain the wind pressure coefficients, with STAR-CCM+ (Ver. 15.06) as the analysis software. When using STAR-CCM+, we utilized a computer cluster subsisting of six servers, each equipped with an Intel Xeon Gold 6150 CPU featuring 108 cores. Domain boundaries were placed sufficiently far away to avoid their effect on the airflow around the ship. The inlet and outlet boundaries were placed at 3.5 Lpp from the hull, while the top boundary was placed at 0.5 Lpp from the hull. The left-and right-side boundaries were set at 2 Lpp from the hull. The boundary conditions were as follows: Velocity distribution at the inlet boundary; Constant pressure is applied at the outlet boundary; Both side boundaries are symmetric surfaces; Symmetric boundary conditions are applied to the top boundary while the bottom boundary is set as a non-slip wall according to the type of numerical simulation. The reason for setting the non-slip wall boundary condition at the bottom is that a natural boundary layer is formed due to the wind. The velocity distribution is constant, independent of the height since the airflow generated by the ship's motion does not form a boundary layer. The ship's surface is defined as a non-slip wall boundary condition. Special attention was paid to the discretization around the wind-loaded region of the ship. The prism layer was carefully created to maintain the non-dimensional wall distance of Y + > 30. The distance between the center of the first cell near the wall and the wall is defined by the following equation [33]. Figure 8 shows the computational domain with the interface mesh enclosing the ship along with the described boundary conditions. center of the first cell near the wall and the wall is defined by the following equation [33].   The JBC ship was selected for the calculations during the grid validity verification. As shown in the table, to increase the accuracy of the numerical calculations and reduce the computation time, the convergence speed of the calculations is determined using the CFL (Courant-Friedrichs-Lewy) number. Five cases are selected to verify the influence of mesh partitioning on the simulation drag coefficient to identify the optimal mesh partitioning method. The mesh partitioning methods for the five cases are: high refinement, refinement, conventional, coarsening, and high coarsening (Figure 9). The number of generated grids varies with the change in CFL, as shown in Table 2. (20) δt: time step. δx: minimum grid size. |U|: absolute value of velocity.  The JBC ship was selected for the calculations during the grid validity verification. As shown in the table, to increase the accuracy of the numerical calculations and reduce the computation time, the convergence speed of the calculations is determined using the CFL (Courant-Friedrichs-Lewy) number. Five cases are selected to verify the influence of mesh partitioning on the simulation drag coefficient to identify the optimal mesh partitioning method. The mesh partitioning methods for the five cases are: high refinement, refinement, conventional, coarsening, and high coarsening (Figure 9). The number of generated grids varies with the change in CFL, as shown in Table 2.
δt: time step. δx: minimum grid size. |U|: absolute value of velocity. When the wind direction at the inlet is set to 180 degrees, the ship is only subjected to drag from the X direction. As shown in Figure 10, the calculated wind resistance coefficient decreases as the CFL value increases. However, when the CFL value is lower than or equals two, the wind resistance coefficient does not change significantly as the CFL value increases. When the CFL value is more than two, the wind resistance coefficient changes drastically as the CFL value increases, so the results cannot be calculated accurately when the CFL value is greater than two. However, at CFL values of less than two, the wind resistance coefficient remains stable while the number of meshes increases exponentially, resulting in a longer computation time. Therefore, in this paper, a CFL value of two was chosen for all grid divisions to ensure accurate and efficient calculation results.  When the wind direction at the inlet is set to 180 degrees, the ship is only subjected to drag from the X direction. As shown in Figure 10, the calculated wind resistance coefficient decreases as the CFL value increases. However, when the CFL value is lower than or equals two, the wind resistance coefficient does not change significantly as the CFL value increases. When the CFL value is more than two, the wind resistance coefficient changes drastically as the CFL value increases, so the results cannot be calculated accurately when the CFL value is greater than two. However, at CFL values of less than two, the wind resistance coefficient remains stable while the number of meshes increases exponentially, resulting in a longer computation time. Therefore, in this paper, a CFL value of two was chosen for all grid divisions to ensure accurate and efficient calculation results.

Estimation of the Wind Force for Various Ships
The detailed configuration of the fine mesh used in the numerical simulations to determine the wind resistance is shown in Figure 11.  When the wind direction at the inlet is set to 180 degrees, the ship is only subjected to drag from the X direction. As shown in Figure 10, the calculated wind resistance coefficient decreases as the CFL value increases. However, when the CFL value is lower than or equals two, the wind resistance coefficient does not change significantly as the CFL value increases. When the CFL value is more than two, the wind resistance coefficient changes drastically as the CFL value increases, so the results cannot be calculated accurately when the CFL value is greater than two. However, at CFL values of less than two, the wind resistance coefficient remains stable while the number of meshes increases exponentially, resulting in a longer computation time. Therefore, in this paper, a CFL value of two was chosen for all grid divisions to ensure accurate and efficient calculation results.

Estimation of the Wind Force for Various Ships
The detailed configuration of the fine mesh used in the numerical simulations to determine the wind resistance is shown in Figure 11.

Estimation of the Wind Force for Various Ships
The detailed configuration of the fine mesh used in the numerical simulations to determine the wind resistance is shown in Figure 11.

Wind Force of JBC Ship
The wind force was separated by the lengthwise and lateral directions of the ship and was dimensionless, as calculated by Equation (7). The X-direction indicates the direction of the wind, 0 degrees indicates that the wind is blowing from the bow, 90 degrees shows that the wind is blowing from the starboard side, and 120 degrees shows that the wind is blowing from the starboard side aft.
The wind resistance of the ship was calculated with a total of seven cases from 0 to 180 degrees with 30-degree intervals. Resulting in the calculations derived to determine the wind resistance coefficients in the X-and Y-directions and the moment coefficients in the Z-direction, respectively, with the same comparison trend as the experimental data of KUME Kenichi [40] (Figures 12 and 13). It was found that the longitudinal and lateral forces are maximum when the wind blows at 150 and 90 degrees, respectively.

Wind Force of JBC Ship
The wind force was separated by the lengthwise and lateral directions of the ship and was dimensionless, as calculated by Equation (7). The X-direction indicates the direction of the wind, 0 degrees indicates that the wind is blowing from the bow, 90 degrees shows that the wind is blowing from the starboard side, and 120 degrees shows that the wind is blowing from the starboard side aft.
The wind resistance of the ship was calculated with a total of seven cases from 0 to 180 degrees with 30-degree intervals. Resulting in the calculations derived to determine the wind resistance coefficients in the X-and Y-directions and the moment coefficients in the Z-direction, respectively, with the same comparison trend as the experimental data of KUME Kenichi [40] (Figures 12 and 13). It was found that the longitudinal and lateral forces are maximum when the wind blows at 150 and 90 degrees, respectively.

Wind Force of JBC Ship
The wind force was separated by the lengthwise and lateral directions of the ship and was dimensionless, as calculated by Equation (7). The X-direction indicates the direction of the wind, 0 degrees indicates that the wind is blowing from the bow, 90 degrees shows that the wind is blowing from the starboard side, and 120 degrees shows that the wind is blowing from the starboard side aft.
The wind resistance of the ship was calculated with a total of seven cases from 0 to 180 degrees with 30-degree intervals. Resulting in the calculations derived to determine the wind resistance coefficients in the X-and Y-directions and the moment coefficients in the Z-direction, respectively, with the same comparison trend as the experimental data of KUME Kenichi [40] (Figures 12 and 13). It was found that the longitudinal and lateral forces are maximum when the wind blows at 150 and 90 degrees, respectively.   Figure 14 shows the moment coefficients at the center of gravity for the JBC ship. The moment coefficients were amplified by a factor of ten due to their small values, of which a maximum value of about −1.3 was observed at approximately 120 degrees. The calculated moment coefficients will be used as input data for future control systems.  Figure 14 shows the moment coefficients at the center of gravity for the JBC ship. The moment coefficients were amplified by a factor of ten due to their small values, of which a maximum value of about −1.3 was observed at approximately 120 degrees. The calculated moment coefficients will be used as input data for future control systems.  Figure 14 shows the moment coefficients at the center of gravity for the JBC ship. Th moment coefficients were amplified by a factor of ten due to their small values, of whic a maximum value of about −1.3 was observed at approximately 120 degrees. The calcu lated moment coefficients will be used as input data for future control systems. In the validation of grid validity, CFL = 2 is used for grid division. The experimenta data of the JBC ship model obtained from KUME Kenichi and the calculation results o the JBC ship in actual dimensions are shown as red and green, respectively, as shown i Figures 12 and 13. Figure 12 shows that the error of the X-direction damping coefficient i larger at 60 and 120 degrees because the CFD model is the result of the actual ship siz calculation, and the difference in the area can be ignored because the X-direction forc does not have a significant effect on the ship's mooring state. Figure 13 shows that th CFD results are almost the same as the experimental results. Therefore, it can be verifie again that the chosen meshing method is appropriate for the calculations. Figure 15 shows how the velocity distribution around the ship is extracted at 60-de gree intervals, depending on the wind direction, at Z = 10 m. At 0 and 180 degrees, it i difficult to observe the main features of the flow, but at 60 and 120 degrees, a stagnan flow marked in blue can be seen on the face colliding with the wind, and vortex separatio in the opposite direction can also be observed. Therefore, it can be inferred that latera forces will be found to have an important role. In the validation of grid validity, CFL = 2 is used for grid division. The experimental data of the JBC ship model obtained from KUME Kenichi and the calculation results of the JBC ship in actual dimensions are shown as red and green, respectively, as shown in Figures 12 and 13. Figure 12 shows that the error of the X-direction damping coefficient is larger at 60 and 120 degrees because the CFD model is the result of the actual ship size calculation, and the difference in the area can be ignored because the X-direction force does not have a significant effect on the ship's mooring state. Figure 13 shows that the CFD results are almost the same as the experimental results. Therefore, it can be verified again that the chosen meshing method is appropriate for the calculations. Figure 15 shows how the velocity distribution around the ship is extracted at 60-degree intervals, depending on the wind direction, at Z = 10 m. At 0 and 180 degrees, it is difficult to observe the main features of the flow, but at 60 and 120 degrees, a stagnant flow marked in blue can be seen on the face colliding with the wind, and vortex separation in the opposite direction can also be observed. Therefore, it can be inferred that lateral forces will be found to have an important role.

Wind Force of the KCS Ship
The wind resistance of the ship was calculated for a total of seven cases from 0 to 180 degrees and 30-degree intervals, and the results of the calculations were obtained for the wind resistance coefficients in the X-and Y-directions and the moment coefficients in the Z-direction, respectively (Figures 16 and 17). It was found that the longitudinal and lateral forces are maximum when the wind direction is at 150 and 90 degrees, respectively. Figure 18 shows the moment coefficients at the center of gravity for the JBC ship. The moment coefficients were amplified by a factor of ten due to their small values. The maximum value occurred at approximately 120 degrees, measuring approximately −1.5. The calculated moment coefficients will be used as input data for future control systems. Figure 19 shows how the velocity distribution around the ship is extracted at 60-degree intervals, depending on the direction of the wind, at the Z = 10 m position. At 0 and 180 degrees, it is difficult to observe the main features of the flow, but at 60 and 120 degrees, a stagnant flow marked in blue can be seen on the face colliding with the wind, and vortex separation in the opposite direction can also be observed. Therefore, it can be inferred that lateral forces will be found to have an important role.

Wind Force of the KCS Ship
The wind resistance of the ship was calculated for a total of seven cases from 0 to 180 degrees and 30-degree intervals, and the results of the calculations were obtained for the wind resistance coefficients in the X-and Y-directions and the moment coefficients in the Z-direction, respectively (Figures 16 and 17). It was found that the longitudinal and forces are maximum when the wind direction is at 150 and 90 degrees, respectively   Figure 18 shows the moment coefficients at the center of gravity for the JBC shi moment coefficients were amplified by a factor of ten due to their small values. The imum value occurred at approximately 120 degrees, measuring approximately −1. calculated moment coefficients will be used as input data for future control system  Figure 19 shows how the velocity distribution around the ship is extracted at gree intervals, depending on the direction of the wind, at the Z = 10 m position. At Z-direction, respectively (Figures 16 and 17). It was found that the longitudinal and forces are maximum when the wind direction is at 150 and 90 degrees, respectively   Figure 18 shows the moment coefficients at the center of gravity for the JBC shi moment coefficients were amplified by a factor of ten due to their small values. The imum value occurred at approximately 120 degrees, measuring approximately −1. calculated moment coefficients will be used as input data for future control system   forces are maximum when the wind direction is at 150 and 90 degrees, respectively   Figure 18 shows the moment coefficients at the center of gravity for the JBC shi moment coefficients were amplified by a factor of ten due to their small values. The imum value occurred at approximately 120 degrees, measuring approximately −1. calculated moment coefficients will be used as input data for future control system   vortex separation in the opposite direction can also be observed. Therefore, it can be inferred that lateral forces will be found to have an important role.

Wind Force of the KVLCC Ship
The wind resistance of the ship was calculated for a total of seven cases from 0 to 180 degrees and 30-degree intervals, and the results of the calculations were obtained for the wind resistance coefficients in the X-and Y-directions and the moment coefficients in the z-direction, respectively (Figures 20 and 21). It was found that the longitudinal and lateral forces are maximum when the wind direction is at 150 and 90 degrees, respectively. The wind resistance of the ship was calculated for a total of seven cases from 0 to 180 degrees and 30-degree intervals, and the results of the calculations were obtained for the wind resistance coefficients in the X-and Y-directions and the moment coefficients in the z-direction, respectively (Figures 20 and 21). It was found that the longitudinal and lateral forces are maximum when the wind direction is at 150 and 90 degrees, respectively.   Figure 22 shows the moment coefficients at the center of gravity for the JBC ship. The moment coefficients were amplified by a factor of ten due to their small values. The maximum value occurred at approximately 120 degrees, measuring approximately −1.1. The calculated moment coefficients will be used as input data for future control systems.  The wind resistance of the ship was calculated for a total of seven cases from 0 to 180 degrees and 30-degree intervals, and the results of the calculations were obtained for the wind resistance coefficients in the X-and Y-directions and the moment coefficients in the z-direction, respectively (Figures 20 and 21). It was found that the longitudinal and lateral forces are maximum when the wind direction is at 150 and 90 degrees, respectively.   Figure 22 shows the moment coefficients at the center of gravity for the JBC ship. The moment coefficients were amplified by a factor of ten due to their small values. The maximum value occurred at approximately 120 degrees, measuring approximately −1.1. The calculated moment coefficients will be used as input data for future control systems.   Figure 22 shows the moment coefficients at the center of gravity for the JBC ship. The moment coefficients were amplified by a factor of ten due to their small values. The maximum value occurred at approximately 120 degrees, measuring approximately −1.1. The calculated moment coefficients will be used as input data for future control systems. The wind resistance of the ship was calculated for a total of seven cases from 0 to 180 degrees and 30-degree intervals, and the results of the calculations were obtained for the wind resistance coefficients in the X-and Y-directions and the moment coefficients in the z-direction, respectively (Figures 20 and 21). It was found that the longitudinal and lateral forces are maximum when the wind direction is at 150 and 90 degrees, respectively.   Figure 22 shows the moment coefficients at the center of gravity for the JBC ship. The moment coefficients were amplified by a factor of ten due to their small values. The maximum value occurred at approximately 120 degrees, measuring approximately −1.1. The calculated moment coefficients will be used as input data for future control systems.    Figure 23 shows how the velocity distribution around the ship is extracted at 60-degree intervals, depending on the direction of the wind, at the Z = 10 m position. At 0 and 180 degrees, it is difficult to observe the main features of the flow, but at 60 and 120 degrees, a stagnant flow marked in blue can be seen on the face colliding with the wind, and vortex separation in the opposite direction can also be observed. Therefore, it can be inferred that lateral forces will be found to have an important role.  Based on the calculations, the wind resistance coefficients in the X-and Y-directions of the JBC, KCS, and KVLCC ships were obtained. Considering these coefficients, the maximum wind force in the X-and Y-directions for each ship at wind speeds of 14 m/s and 30 m/s were further calculated as shown in Table 3.

Estimation of Current Force for Various Ships
Under the full-load condition, the submerged areas of the JBC, KCS, and KVLCC ships are 19,556 m 2 , 9330 m 2 , and 27,320 m 2 , respectively, with a flow speed of 2 knots (1 m/s). In addition, we assume that the fender and port support structure can provide sufficient support, so the lateral (Y-axis direction) force caused by the current is neglected. According to the formula, the current force is calculated as shown in Table 4.

Estimation of Wave Force for Various Ships
The calculation of wave forces on the ships is performed using the ship hydrodynamic calculation software ANSYS AQWA, and the mapped mesh method is employed within the AQWA Workbench to conduct more refined time-domain simulations for JBC, KCS, and KVLCC. Figure 24 shows the mesh samples of the model simulated by ANSYS AQWA.
of the JBC, KCS, and KVLCC ships were obtained. Considering these coefficients, the maximum wind force in the X-and Y-directions for each ship at wind speeds of 14 m/s and 30 m/s were further calculated as shown in Table 3.

Estimation of Current Force for Various Ships
Under the full-load condition, the submerged areas of the JBC, KCS, and KVLCC ships are 19,556 m 2 , 9330 m 2 , and 27,320 m 2 , respectively, with a flow speed of 2 knots (1 m/s). In addition, we assume that the fender and port support structure can provide sufficient support, so the lateral (Y-axis direction) force caused by the current is neglected. According to the formula, the current force is calculated as shown in Table 4.

Estimation of Wave Force for Various Ships
The calculation of wave forces on the ships is performed using the ship hydrodynamic calculation software ANSYS AQWA, and the mapped mesh method is employed within the AQWA Workbench to conduct more refined time-domain simulations for JBC, KCS, and KVLCC. Figure 24 shows the mesh samples of the model simulated by ANSYS AQWA. In ship wave force, the maximum force is generated in the pitch direction, followed by the heave and roll directions. Therefore, in this study, the surge, sway, and yaw directions are assumed to be fixed. Since no long-period waves were present when the ship In ship wave force, the maximum force is generated in the pitch direction, followed by the heave and roll directions. Therefore, in this study, the surge, sway, and yaw directions are assumed to be fixed. Since no long-period waves were present when the ship moored in Busan New Port, a significant wave height of 0.5 m and a wave period of 3.5 s are assumed. The longitudinal wave forces for the JBC, KCS, and KVLCC ships were determined using ANSYS AQWA, as shown in Figure 25. Additionally, the transverse forces induced by waves are neglected because they are assumed to be adequately supported by the fenders and port structures.
J. Mar. Sci. Eng. 2023, 11, x FOR PEER REVIEW 21 of moored in Busan New Port, a significant wave height of 0.5 m and a wave period of 3. are assumed. The longitudinal wave forces for the JBC, KCS, and KVLCC ships were d termined using ANSYS AQWA, as shown in Figure 25. Additionally, the transverse forc induced by waves are neglected because they are assumed to be adequately supported the fenders and port structures.

Estimation of Total Mooring Force for Various Ships
When JBC, KCS, and KVLCC ships are at their design draft with a wind speed of m/s and a tidal current of 2 knots, considering a significant wave height of 1 m and a wa period of 3.5 s, the estimated forces that the ships may experience under normal sea co

Estimation of Total Mooring Force for Various Ships
When JBC, KCS, and KVLCC ships are at their design draft with a wind speed of 14 m/s and a tidal current of 2 knots, considering a significant wave height of 1 m and a wave period of 3.5 s, the estimated forces that the ships may experience under normal sea conditions are calculated and presented in Table 5. When the JBC, KCS, and KVLCC ships are at their design drafts with a wind speed of 30 m/s, a current speed of 2 knots, and significant wave heights of 0.5 m; with wave periods 3.5 s, the estimated forces exerted on the ships under typical sea conditions are presented in Table 6. The results in Tables 5 and 6 show that the wind force is the primary external force acting on the ship in the moored state. Due to the small projected area in the longitudinal direction and the absence of strong currents or waves, the wind force becomes relatively small. However, in the transverse direction of the ship, the presence of upper structures and other factors results in a larger projected area, leading to a higher wind resistance. When the wind speed is 14 and 30 m/s, the wind resistance accounts for over 90% of the total resistance. Thus, the wind force has the maximum impact on the mooring system when the ship in the moored state. The calculated ships in this paper are the most operated ships in ports, including the JBC, KCS, and KVLCC ships which represent cargo ships, container ships, and oil tankers, respectively. By calculating the values of C X , C Y , and C M for these three types of ships, the forces F X and F Y can be estimated. Furthermore, they can be used to determine whether the ship is suitable for operations in the harbor and the number of mooring devices required since a ship in a moored state does not connect to a single mooring device.

Conclusions
In this study, an evaluation was conducted to estimate the mooring forces acting on JBC, KCS, and KVLCC for the basic design of an automatic mooring device suitable for large ships in the future. Most mooring forces result from the wind force, while the forces caused by waves and currents are relatively small.
Considering factors such as wind, waves, and currents, the maximum mooring forces required in the longitudinal direction for JBC, KCS, and KVLCC, with a wind speed of 14 m/s, are 41.2, 38.7, and 77.2 kN, respectively, which are relatively small. However, in the transverse direction, due to the influence of upper structures, the maximum mooring forces required are 340, 837, and 222 kN, respectively. Notably, under extremely severe weather conditions with a wind speed of 30 m/s, the maximum mooring forces required in the longitudinal direction for JBC, KCS, and KVLCC are 43, 132.7, and 159.2 kN, respectively, which are relatively small. However, in the transverse direction, due to the influence of upper structures, the maximum mooring forces required are 1674, 4118, and 1091 kN, respectively. Therefore, these factors should be considered in the design of an automatic mooring device. The calculation results provide information on the mooring forces required for the ship in the moored state and the number of mooring devices needed for the operation. By calculating the results for a wind speed of 30 m/s, it can be determined whether the ship can be moored in the port during a typhoon. If the forces acting on the ship exceed the capacity of the mooring system, separation between the mooring system and the ship may occur, leading to significant accidents at the port. In this study, the wind resistance coefficients for the three types of ships that are most operated in the port have been calculated, which could serve as a reference for estimating the forces for other types of ships. This can enhance the safety of ships under mooring conditions. Furthermore, the calculated mooring forces did not consider the maritime traffic conditions, such as additional force caused by the interaction among navigating ships. Additionally, if long-period waves exist within the port, the possibility of resonance with the ship cannot be excluded. Hence, future studies may require the use of tension calculation methods to analyze the ship's motion.

Data Availability Statement:
The study did not report any data.