Numerical Analysis of Hydrodynamic Loads on Passing and Moored Ships in Shallow Water

In this study, hydrodynamic interactions between passing and moored ships were studied by applying a time-domain numerical simulation method. The boundary value problem for a fluid domain was formulated based on a potential flow theory. A numerical method was developed based on a finite element method with an efficient re-mesh algorithm. Regarding the free-surface boundary conditions, both double-body and free-surface models were considered for examining the free-surface effect on the hydrodynamic forces due to the passing ship. First, numerical results were validated by comparison with the model test results of Kriebel et al. (2005), where generic Series 60 hulls were considered as the target model for the passing and moored ships. In addition, hydrodynamic pressure fields and force time-series were investigated to understand the passing ship problem. Second, a series of numerical simulations were performed to study the effects of the passing ship speed, separation distance, water depth, and relative vessel size, which were used to compare the peak values of hydrodynamic forces. The applicability and limitations of the double-body and free-surface models are discussed for predicting passing ship loads.


Introduction
In ports or coastal areas, an advancing ship is often observed to pass near moored or anchored ships. In particular, as the speed and size of ships have increased recently, the effect of large and high-speed passing ships has become an important consideration for the safety of moored ships. It is known that significant changes in the hydrodynamic pressure field around a passing ship can cause an excessive horizontal motion in moored ships, which, in turn, may result in collisions with other ships or the quay wall. For instance, in the worst-case scenario, damage to a berthing or offloading ship's mooring and mechanical equipment could result in a massive explosion and injuries [1]. Consequently, for the safe design of the mooring system of anchored vessels or the establishment of reasonable operating regulations-such as the passing ship speed and separation distancecomputation and prediction methods for the effect of passing ship loads on a moored ship are necessary.
Researchers have conducted many model tests for the passing ship problem in order to understand the fundamental physics and directly measure the passing ship loads acting on a moored ship. Thus, various empirical formulas have been proposed based on model test results. Remery [2] conducted a 1:60 scale model test to examine a large tanker moored at an oil loading terminal in shallow water. Similarly, Muga [3] performed a 1:68 scale model test for supertankers and found that passing ship loads tend to increase with reduced separation distance or increased vessel speed or vessel size. Vantorre [4] conducted a 1:75 scale model test for a bulk carrier, container ship, and tanker and analyzed ship-ship interaction loads for encountering and overtaking cases. While previous experimental Processes 2021, 9,558 2 of 17 studies on the passing ship problem were performed on specific ship models with limited variables, Kriebel [5] performed a 1:135 scale model test using a general Series 60 ship model with several variables, such as passing vessel speed, vessel displacement, water depth, and separation distance. Moreover, by analyzing these experimental results, Kriebel developed empirical equations and evaluated the applicability of the two existing models of Flory [6] and Seelig [7] on the passing ship problem. Swiegers [8] also conducted a 1:100 scale model test for a moored bulk carrier and a passing container ship using several variables such as the passing speed, separation distance, water depth, and presence of walls and channels.
In addition to the model test approach, researchers have also attempted to solve the passing ship problem by introducing semi-empirical formulas based on slender-body theory [7,9]. However, the application of two-dimensional (2D) numerical methods for three-dimensional (3D) vessels has limitations, and performing numerous model tests to identify empirical coefficients using various ship shapes and sizes is difficult. Accordingly, in the 2000s, 3D potential flow models considering detailed shape and the surrounding terrain features of ship were developed to establish a generalized prediction method for passing ship loads. Initially, Pinkster [10] adopted potential flow methods based on the boundary element method for low-speed passing ships and performed a comparative study by applying the double-body and free-surface models. The analysis of free-surface effects demonstrated the applicability of the double-body model for the passing ship problem compared to the free-surface model. However, in this study, the presence of the moored ship was not considered in the free-surface model, and the free-surface effect was not analyzed rigorously enough. Molen [11] also simulated the passing ship problem using the potential flow model based on the double-body flow assumption applied by Pinkster [10] and validated their numerical results against Swiegers' model test [8]. In his study, simulations of the passing ship problem in open water or with the presence of a quay wall or channel wall provided a good approximation of the model test results. Recently, unlike previous studies where numerous 3D potential flow models were applied based on the boundary element method, Nam [12] and Park [13] solved the passing ship problem using the classical finite element method (FEM). Nam [12] analyzed complex hydrodynamic interaction forces, including the presence of passing and moored ships and the influence of surrounding terrain features, such as the seabed and sidewall. He studied a relatively high-speed passing ship and a moored barge alongside a quay by varying the simulation parameters, such as the passing ship speed, wall distances, separation distance, and waves, including the free-surface effect.
Direct numerical simulations based on computational fluid dynamics (CFD) have also been applied for the passing ship problem. Huang [14] analyzed the influence of the passing ship layout, speed, and separation distance on vessels moored at piers using a Reynolds-averaged Navier-Stokes (RANS)-based CFD model. Wang [15] also used the unsteady RANS-based CFD model to analyze the hydrodynamic interaction that occurs in a single-way lock. These studies were conducted by applying an undisturbed free-surface boundary condition using the double-body model, but the free-surface effects were not considered. Recently, with the rapid development of the CFD technique and computational hardware, CFD simulations have been increasingly used to analyze the passing ship problem, including the effects of vorticity and viscosity in limited passing ship problem conditions. However, thus far, this method still requires enormous computational time and hardware costs compared to existing potential flow models, and additional verification procedures are needed for accurate and reliable CFD simulation results.
In this study, we investigated the passing ship problem in order to predict the hydrodynamic forces acting on a moored ship resulting from a passing ship. The computations were performed using the 3D potential flow model based on the classical FEM. Two types of free-surface boundary conditions were applied: the double-body and free-surface models. The computational results were compared with those of Kriebel's model test [5] for validation. Simulation parameters, such as the passing ship speed, separation distance, Processes 2021, 9, 558 3 of 17 and water depth, were analyzed using generic Series 60 hulls as the passing and moored ships. Moreover, the relative length, breadth, and draft of the moored ship were varied to investigate the influence of the vessel size.

Boundary Value Problem
A numerical method based on 3D potential flow theory was developed for the computations of the passing ship problem. Here, an inviscid and incompressible fluid and an irrotational flow were assumed. Two numerical models, double-body and free-surface models, were applied to solve boundary value problems for the fluid domain. In the double-body model, the free surface is assumed as a rigid wall boundary under a low-speed condition of a passing ship, which was adopted by Korsmeyer [16] and Pinkster [10,17] for the passing ship problem. The boundary value problem of the double-body model is as follows (Figure 1): 2) becomes zero. Equation (3) represents the normal velocity boundary condition at the mean water level (z = 0) and on the seabed (z = −h). Previous studies confirmed that the suction effect, a primary factor of low-speed passing ships, can be calculated through the double-body model [18]. As this model neglects free-surface fluctuations, the time-dependent terms disappear in the boundary conditions. Therefore, the time integration method is not required, which reduces the total calculation time significantly. As a result, predicting real-time passing ship loads can be made possible using the double-body model [19].
Unlike the double-body model, the free-surface model applied Equations (4) and (5) as the free-surface boundary conditions instead of Equation (3). Here, ζ is the wave elevation on the free surface, and g is the gravitational constant.

Numerical Method: Finite Element Method
In this study, a classical FEM was introduced to solve the Laplace equation in the fluid domain. The weak formulation of the Laplace equation was obtained by multiplying the test function ψ and applying the integration by parts, which results in Equation (6).
In the FEM, the fluid domain is discretized with a finite number of elements, and the velocity potential and wave elevation are approximated with 3D and 2D elements, respectively, as shown in Equations (7) and (8). In this study, the velocity potential was approximated using the linear summation of the product of the nodal velocity potential and the 3D basis function of eight-node hexahedral elements. The wave elevation was approximated using the linear summation of the product of the nodal wave elevation and the 2D basis function of the four-node quadrilateral elements.
By integrating Equation (6) with Equations (7) and (8), a linear algebraic equation can be derived. In the free-surface model, the wave elevation, as well as the velocity potential, is updated at each time step by applying the time integration of the free-surface boundary conditions. However, in the double-body model, only the velocity potentials are calculated, without any time integration of the free surface. In this study, the conjugate gradient method and the fourth-order Adams-Bashforth-Moulton method were used for matrix solving and the time integration on the free surface, respectively. The numerical method applied here can also be validated by solving a conventional ship resistance problem (Nam [12]).
The computation procedure for solving the passing ship problem is depicted in Figure 2. It can be divided into two parts according to the free-surface boundary conditions, i.e., the double-body and free-surface models. First, for the simulation, the program reads inputs of computation variables and ship geometry. The computation variables include the analysis domain, grid information, and simulation parameters. The ship geometry may be defined as x-and z-coordinate information for each section of the ship. Then, according to the horizontal position of the passing ship, the computational mesh is modified based on the re-mesh algorithm proposed by Nam [12]. In this study, the mesh generation technique was supplemented so that the local mesh is generated by extending it to the seabed or sidewall, which can be applied to the simulation case where the passing ship moves at a shallow water depth or has a low separation distance from the moored ship.
In the re-mesh algorithm, the mesh is divided into background, local and global meshes, as shown in Figure 3. First, the background mesh is generated for the entire computational domain, while local meshes are generated in a smaller area as body-fixed meshes containing the passing ship or moored ship hulls. As the simulation time progresses, the local mesh of the passing ship is moved according to the forward speed, whereas that of the moored ship is fixed at a specific location. Both local meshes are placed at their corresponding locations on the background mesh. Subsequently, the background meshes, which overlap with local meshes, are extracted and replaced with local meshes. Finally, by connecting their interfaces, an earth-fixed global mesh is generated. For a more generalized connection process, a background mesh is generated based on a rectangular 3D rectilinear grid, and local meshes are generated based on a rectangular 3D curvilinear grid considering that the passing ship moves parallel to the moored ship. In this way, both mesh interfaces Processes 2021, 9, 558 5 of 17 become a 2D rectilinear grid, and the global mesh based on the 3D rectilinear grid is completed through a straightforward connecting process.

Convergence Study
Numerical studies were conducted on mesh convergence to determine the appropriate mesh size for accurate computations of the hydrodynamic forces due to the passing ship. Here, the mesh size was represented by the number of elements per mean ship length (L/dx), which varies from 30 to 120. Figure 4a shows the convergence test results in terms of the time histories of sway force according to different mesh sizes. It is evident that all the computational results were almost the same regardless of the mesh size, which indicates the convergence of numerical solutions. For a more precise comparison, the positive peak values of the surge and sway forces and yaw moment can be compared, as shown in Figure 4b. The computational results were represented by the difference in each load compared to that of the finest mesh (L/dx = 120). It is confirmed that the largest difference was less than 2.5% of the yawing moment for the coarsest grid. Thus, all hydrodynamic forces converged quickly as the mesh size decreased. Hereafter, numerical computations were performed by applying L/dx = 60, in which all load coefficients converge within an error of less than 0.5%.

Simulation Description
The numerical simulations were performed by varying the Froude number, the separation distance (S) between the passing and moored ships, and the water depth (h), as shown in Figure 1. Regarding the Froude numbers, the length-based Froude number is defined as F L = V p / gL, while the depth-based Froude number defined as F h = V p / gh. The computational results were validated by comparison with the data from Kriebel's model tests [5]. Generic Series 60 hull forms were used as the passing and moored ship hulls ( Figure 5, [20]), and the variable ranges in this study, as listed in Table 1, were the same as those of Kriebel's model test [5]. Further analysis was also conducted by varying the moored ship size (L m × B m × D m ) to examine its effect. In each part of the analysis, the target variable was changed within a given range (Table 2), whereas the other variables were fixed using the representative values summarized in Table 1.

Double-Body vs. Free-Surface Models
First, the hydrodynamic forces acting on a moored ship were calculated when a passing ship travels far from the moored ship and moves parallel to it. In this case, the distance between the moored and passing ships initially becomes closer as time passes. Then, the passing ship moves next to the moored ship, where the hydrodynamic interaction becomes most significant. Finally, the passing ship moves away. The horizontal hydrodynamic forces and moment acting on the moored ship due to the passing ship are shown in Figure 6. Here, the x-axis represents the non-dimensional horizontal distance, which is the ratio of the horizontal distance between two ships to the mean ship length. The y-axis represents the non-dimensional force and moment coefficients defined by Equation (9). In these plots, the computational results of the double-body (DB) and free-surface (FS) models are directly compared and the peak values from Kriebel's model test [5] are plotted together as dash-dot horizontal lines. As shown in Figure 6, the two numerical results show similar trends, including the positions at which the force peaked. However, it can be observed that the amplitudes of the surge, sway forces, and yaw moment of the FS model were slightly larger than those of the DB model. In addition, small wave-frequency force components are also observed due to the small waves caused by the passing ship in the numerical simulation of the FS model. When comparing the experimental data, the results of the FS model are relatively closer to the experimental data than that of the DB model. The speed of the passing ship is known to be the most important factor in the hydrodynamic problem of a passing ship [3,5]. In this study, numerical simulations were conducted according to various Froude numbers, as listed in Table 2. The present computational results were compared with the experimental data in Figure 7. The x-axis represents the Froude number, whereas the y-axis shows the peak values of the passing ship forces and moment. As expected, the hydrodynamic forces and moment significantly increase with the speed of the passing ship. Regarding the surge and sway forces, it can be seen that the computational results of the FS model are closer to the model test data than those of the DB model. In the case of the DB model, the hydrodynamic forces increase in proportion to the square of the passing ship speed. However, the model test and the FS model show that the passing ship forces increase more than the square of the passing ship speed due to the free-surface effect. Kriebel et al. [21] also reported that, if the depth-based Froude number is greater than 1.0 and the passing ship speed exceeds the wave celerity, the ship wave can be generated excessively, causing a significant free-surface effect. In addition, the free-surface effect also becomes important when the length-based Froude number is 0.4, and the transverse wavelength equals the ship length. As our simulation results indicate, the DB model could estimate the passing ship loads with similar accuracy to the FS model until the length-based Froude number became 0.12. However, if the passing ship speed exceeded this speed range, the results of the FS model and model test resulted in a difference of more than 10% from those of the DB model. When looking at the yaw moment results, the increasing tendency of the yawing moment from two present numerical calculations was similar to that of the model test even though some differences were observed.
For a more detailed analysis, the potential and pressure fields on the free surface were computed using the DB model, as shown in Figure 8. In each column, x/L of the passing ship (x: defined in Figure 1) corresponds to −0.5, 0, and +0.5, hereinafter referred to as P1, P2, and P3, respectively. When compared with the time series of the surge force in Figure  6a, the surge force showed negative and positive peak values at P1 and P3, respectively. At these locations, the center of the moored ship is close to the positive and negative pressure field regions. When the passing ship arrives at P2, the midships of the passing and moored ships are abreast, and it can be observed that an almost symmetrical pressure field also forms around the moored ship. As a result, the surge force acting on the moored ship becomes nearly zero at P2. Therefore, the time series of the surge force shows a sine-wave shape symmetrical at x/L = 0, as shown in Figure 6a.   Regarding the sway force, the maximum positive value is found at P2, where the moored ship is located at the center of the negative pressure field. The positive sign of the sway force indicates that the moored ship is pulled toward the passing ship. When x/L is −0.75 and +0.75, maximum negative sway forces are exerted on the moored ship in a direction away from the passing ship. At these positions, the moored vessel is directly affected by the positive pressure fields around the bow or stern of the passing ship, which results in a strong repulsive force between the passing and moored ships. The time series of sway force shows the shape of the cosine wave to be symmetrical about the y-axis, and the positive peak value is approximately twice the negative peak value. This is because, as shown in Figure 8b, when looking at the pressure field that affects the moored ship, the negative pressure field near the middle of the passing ship is almost twice as strong as the positive pressure field around the bow or stern. The yaw moment also shows a sine-wave shape, symmetric at x/L = 0. The negative peak, zero value, and positive peak of the yaw moment were observed at P1, P2, and P3, respectively. It can be seen that the yaw moment acting on the moored ship is determined by the lateral force and its application point.
The potential, pressure, and wave elevation fields obtained by the FS model are shown in Figure 9. The potential and pressure show a similar spatial distribution to the results of the DB model. However, the magnitude of the pressure in the FS model is slightly greater. In addition, a slightly asymmetric pressure field can be found around the bow and stern sections of the moored ship. The wave elevation field of the FS model also shows a similar spatial distribution to that of the pressure field. For a more rigorous comparison of the amplitude of the pressure for both models, the cross-sectional pressure fields at x/L = 0 are calculated in Figure 10. It can be seen that a high pressure field forms around the moored ship because of the wave elevation, which results in the higher forces of the passing ship.

Effects of Separation Distance and Water Depth
For analyzing the effect of separation distance, numerical simulations were performed under various separation distances, which were defined as the horizontal distances between the centerlines of the passing ship and moored ship. The range of the ratio of the separation distance to ship length (S/L) is 0.3-1.0. The minimum separation distance ratio S/L = 0.3 corresponds to the ship length-to-breadth ratio (S/B) of 2.0, which means that the distance between the two ships' sidewalls is equal to the ship breadth. In the present analyses, the Froude number was fixed at 0.142, and the depth-to-draft ratio at 1.324. Figure 11 shows the computational results of the DB model under various separation distance conditions. Here, the time series of the hydrodynamic forces and moment by the passing ship are directly compared. As the separation distance increases, the amplitudes of the passing ship loads significantly reduce, and the phase angles also slightly change. The peak sway forces are observed at the position where x/L is 0, regardless of the separation distance. However, the positions where the maximum surge forces or the maximum yaw moments are observed are far from the intersection position. Here, x/L changes from 0.3 to 0.5 with an increase of separation distance. This is because the maximum surge forces occur when the moored ship is located at the boundary between the positive and negative pressure, as shown in Figure 9b. As the separation distance increases, the boundary also moves farther from the intersection position.  Figure 12 shows the effect of the separation distance on the passing ship forces by comparing the results of the computations and experiments. Here, the x-axis shows two types of non-dimensional separation distance ratios, S/L and S/B. The y-axis denotes the non-dimensional force and moment coefficients divided by the square of the passing ship speed, as given in Equation (9). In general, as the separation distance decreases, the passing ship forces can be expected to increase. It can be clearly seen that the forces and moment sharply increase when the separation distance decreases, especially to less than four times the breadth. As shown in Figure 9b, when the two ships get closer, a strong hydrodynamic pressure field develops between them. Overall, the experimental result and numerical simulations showed similar trends for the passing ship forces and moment according to the various separation distances. For the surge and sway force coefficients, the FS model results showed slightly higher load coefficients than those of the DB model, and the model test results were closer to the former. For all separation distances, the sway forces were about four times larger than the surge force. Regarding the yaw moment coefficient (Cm), the model test data were slightly lower than the computational results. However, the increasing tendency of the yawing moment in both the model test data and numerical results was similar. It should be noted that the differences in computational results of the DB and FS models were still observed to be approximately 10-20% with respect to the values of the DB model. As the passing ship sometimes needs to pass a location in shallow waterways or coasts on a given route, the water depth is also a critical simulation parameter for the passing ship problem. Additional numerical simulations were performed to study its effect on the passing ship loads. Here, the depth-to-draft ratio (h/D) was varied from 1.108 to 3.919, and the other parameters were fixed at their representative values with a length-based Froude number of 0.142 and S/L of 0.4. The lowest water depth ratio, i.e., 1.108, represents a distance between the bottom of the ship and the seabed of only 10% of the ship's draft. Numerical computations were performed by applying both DB and FS models, and the load coefficients were compared with the model test results, as shown in Figure 13. It can be seen that the passing ship forces increase as the water depth decreases. In particular, as the water depth gets shallower, h/D < 1.5, the load increases exponentially due to the seabed effect. Under the deep water-depth conditions of h/D = 1.5 or more, the loads for the computation and model test results were similar, with the effect of the water depth on the ship-generated wave no longer being prominent. However, at shallow depths with h/D < 1.5, the computation and model test results showed significant discrepancies (in the order of the double-body model, free-surface model, and model test results). At shallower depths where h/D <1.2, the difference between the results of the FS model and the experiment significantly increased because, as the water depth decreases, nonlinear effects, such as those of viscosity and nonlinear waves, increase. Calculations for deeper water depth conditions with h/D > 4 were also performed as a representative of a water depth at which the seabed effect would not affect the passing ship load anymore. As a result, the variation of the passing ship load according to the increase of the water depth was less than 5% when h/D > 10. This result means that the pressure field near the water surface is formed without being affected by the seabed effect, and therefore the water-depth effect should be considered if the water depth is shallower than ten times the ship's draft.

Effect of Ship Size
In this section, the size effect on the passing ship problem is discussed. The passing ship problem is a phenomenon caused by the hydrodynamic interaction between passing and moored ships, so their relative sizes are a critical factor for the passing ship problem. The passing ship was selected as a Series 60 hull with a block coefficient of 0.80, while the moored ship was considered as a rectangular barge so that the length, breadth, and draft could be easily varied ( Figure 14 and Table 3). For analyzing the size effect, the sizes of the barge were varied in each direction independently. Numerical computations were performed by applying the DB model. Figure 15 compares the time series of the hydrodynamic forces and moment acting on the Series 60 and barge. The pressure fields on the water surface are also compared in Figure 16. It can be observed that Cy increases slightly, and Cm increases significantly owing to the shape of the moored ship, although Cx remains constant. The reasons for the increases in Cy and Cm were due to increases in the surface area at the bow and stern of the barge and in the pressure difference between the port and starboard sections of the barge. The peak loads were observed at similar horizontal distances for the results of both the Series 60 and barge.   For the analysis of the size effect, the breadth ratio (B/B 0 ) was varied within 75-100%, and the length ratio (L/L 0 ) within 80-120%. Figure 17a shows the changes in hydrodynamic forces according to variations of the barge breadth. The x-axis represents the change in the breadth ratio as a percentage. The y-axis represents the ratio of the force or moment coefficients acting on the re-sized barge compared to the original Series 60 with a block coefficient of 0.75. It can be seen that all the forces and moment linearly increase as the breadth of the moored ship increases. Specifically, as the ship breadth increases by 10%, the surge and sway forces and yaw moment coefficients also increase by approximately 10%.  According to the change in length of the barge, the load coefficients for the surge and sway forces and yaw moment were also calculated and compared to the basic loads acting on the original Series 60, as shown in Figure 17b. As the moored ship length increases, the surge force decreases linearly; however, the yaw moment increases linearly. As the moored ship's length increases, the bow or stern section of the moored ship move to the low-pressure field region, which results in the decrease of the surge force. However, the resultant yaw moment increases because the surface area at the port side on which the lateral force is exerted increases. It can be observed that the maximum sway force occurs when the passing and moored ships have the same length. This is because, if the moored ship length is lower than the passing ship length, the lateral area on which the positive sway force is exerted also reduces. On the contrary, if the length of the moored ship is greater than that of the passing ship, then the area under the negative pressure field remains the same, but the area of the positive pressure field around the bow and stern sections increases, as shown in Figure 16b. As a result, as the ship length ratio (L/L 0 ) increases from 90% to 100%, the surge force is reduced by 2%, sway force increases by 3%, and the yaw moment increases by 15%.
The computations were also conducted by changing the draft ratio (D/D 0 ) from 10% to 120%, as shown in Figure 17c. It can be seen that the surge force linearly increases due to the linear increase in the area at the bow and stern sections from which the pressure fields are formed; however, the sway force and yaw moment exponentially increase with D/h. As the draft increases, the exponential indices of the sway force and yaw moment increase from 2 and 1.4 to 3 and 2, respectively. Compared to the analysis of the breadth and length change, the passing ship load increases more rapidly when the draft varies. This is because the seabed effect is added to the effect of the increase in the surface area due to the reduction in the gap between the bottom of the ship and the seabed under the fixed water-depth condition.

Conclusions
In this study, the hydrodynamic forces acting on a moored ship caused by a passing ship were analyzed through a series of FEM-based numerical computations. The results are summarized as follows. - The hydrodynamic forces and moment acting on the moored ship increase significantly as the speed of the passing ship increases. In particular, the passing ship load increases in proportion to the square of the passing ship speed under low Froude numbers less than 0.12, for which both the DB and FS models can accurately predict the passing ship loads. However, when the Froude number exceeds this limit, the passing ship load increases at a rate faster than the square of the passing ship speed, mainly due to the free-surface effect. Under this condition, it was found that the numerical results based on the FS model were closer to the model test data than those based on the DB model. - The passing ship loads increase gradually with decreasing separation distance between the passing and moored ships. From the present numerical computations, it can be observed that the passing ship loads sharply increase when the separation distance ratio (S/B) becomes less than 4. Under low-speed conditions, the present numerical methods could predict the peak values of the passing ship loads with an error less than 10%, regardless of the separation distance. -As the water depth decreases, the passing ship load also increases. Especially under shallow water-depth conditions where the depth ratio (h/D) < 1.2, it could be clearly observed that the passing ship loads significantly increased. - The relative sizes of the moored and passing ships are a critical factor in the passing ship loads due to changing pressure fields formed around the moored ship. All the forces and moment increase linearly with increasing moored ship breadth. Furthermore, with an increase in the moored ship length, the surge force linearly decreases, and the yaw moment linearly increases. Compared to the changes in the breadth and length, passing ship loads are more sensitive to the change of the ship draft. - In the future, new validation data from model tests are required for a close examination of the free-surface effect in the passing ship problem. In addition, numerical studies on various ships should be performed by applying the developed numerical methods.