1 Introduction
Increasing the size and speed of a vessel gives economic benefits, but at the same time the risk of accidents increases. Thus, understanding a ship's manoeuvrability is important to reduce risk and simulation technology is a representative measure to estimate manoeuvrability. Estimation of the hydrodynamic coefficients for the forces and moments acting on the hull is important in a ship modelling process for the simulation. The International Towing Tank Conference (ITTC) presents a classification for estimation methods and their accuracy in effort/cost characteristics, as shown in Figure 1 and Figure 2 (ITTC, 2008).
Unlike a captive model test and computational fluid dynamics (CFD), which are widely used for ships under construction (Lee, 2015), a system identification using sea trials is an efficient method for existing ships (Kim et al., 2017). This system identification estimates the hydrodynamic coefficients using a mathematical optimization algorithm. The algorithm conducts manoeuvre simulation and compares its result with benchmark data, such as sea trial data, at every iteration and it provides updated target variables to be optimized.
The Extended Kalman Filter (EKF) was used at the beginning of this method (Abkowitz, 1980; Hwang, 1980), and System Based (SB) free running tests were included in the EKF algorithm (Rhee and Kim, 1999; Zahng and Zou, 2011). Other mathematical algorithms were also introduced with the development of computer technology, such as the sequential quadratic programming (SQP) and the Broyden-Fletcher-Goldfarb-Shanno method (BFGS) and an interior point (Saha and Sarker, 2010; Tran et al., 2014).
This study compares the optimization algorithms with the same benchmark from the sea trial and optimization conditions. The interior point and the SQP were selected for comparison and three manoeuvres were used to conduct the optimization.
2 Optimization algorithms
2.1 Interior point
An interior point algorithm finds a solution for linear and nonlinear convex optimization problems by searching the interior of the feasible region. Since von Neumann suggested the concept of the algorithm as the classical simplex method, there has been great progress with the interior point algorithm (Dantzig and Thapa, 2003). Khachiyan provided an ellipsoid method, which is slower than the simplex method for the worst case (Khachiyan, 1980). Karmarkar discovered a faster one than Khachiyan’s and it could apply to convex problems (Karmarkar, 1984). Beyond Karmarkar, the interior point algorithm has several varied subclasses, and it is accepted that a primal-dual method is the most efficient (Gondzio, 2012).
An idea of the primal-dual method starts from that a linear program can be solved by using another linear program, which is the opposite of the original one. The primal problem is the following:(1)
subject to Ax = b, x ≥ 0
where and A is an m × n matrix. The dual problem can be formulated as follows.(2)
subject to
where . The optimality conditions are the following:(3)
-
where X is a form of diagonal matrix.
-
Benefits and disadvantages of the interior point algorithm are as follows (Bertsimas and Tsitsiklis, 1997):
-
Benefits
-
- Good for large, sparse problems
-
- Have polynomial time asymptotic complexity Disadvantages
-
- Slower calculation for small problems than the simplex method
-
- Easily stuck at a local minimum due to the great influence of the starting points
2.2 Sequential quadratic programming
Sequential quadratic programming (SQP) is a reliable method of solving optimization problems that have nonlinear constraints. Like most optimization algorithms, SQP is not a single algorithm, but a concept with which various sub-algorithms are associated (Boggs et al., 1995). The idea of SQP is to solve a nonlinear problem by quadratic programming subproblems to get a better optimization result. A standard problem with nonlinear constraints can be written as follows:(4)
where f is an objective function is an equality constraint and cI is an inequality constraint .
The SQP algorithm solves a problem locally (Gilbert, 2017). At each iteration k, the algorithm calculates a solution to the quadratic programming subproblem with its Lagrangian:(5)
where Mk is the Hessian of the Lagrangian ℒ or an approximation.(6)(7)
The benefit of the SQP algorithm is that each iteration does not have to be a feasible point. It is advantageous when a critical constraint is near the feasible points. However, this also means that an iteration that is not the final minimum is not guaranteed regardless of whether it is in the feasible area.
3 Estimation of hydrodynamic coefficients
3.1 Mathematical model
All simulations are carried out under the three degrees of freedom ship- and earth-fixed coordinate systems. As shown in Fig. 3, the Earth-fixed-coordinate plane and the Ship-fixed-coordinate O-xy place on the undisturbed free surface. The x0 axis points in the direction of a ship's heading, and the z0 and z axis point downward.
where,
In the mathematical model, a ship is considered as a massive and rigid body. Equation (8) presents the forces and moment acting on a vessel according to the Newtonian law of motion (Rheinmetall Defense, 2008).
All forces and moment in the mathematical model are composed of several modules: hull, propeller, rudder and external forces and moments, if necessary. Equation (9) shows the detailed composition of the hydrodynamic forces and moment, especially acting on the hull. This study applied the empirical estimation formulas of Norrbin and Clarke to calculate initial values for the optimization (Norrbin, 1971; Clarke et al., 1983). As shown in Equation (10), every coefficient can be formalized by the function of the ship's main dimensions: length, beam, draught, displacement and trim. Ynon and Nnon are non-linear components for sway force and yaw moment.
3.2 Benchmark data
Sea trials are carried out with the training ship 'Hanbada' to collect the benchmark data. A preparation for the sea trials is borrowed from the trials of the training ship ‘Sae Yu Dal’ (Im et al., 2011). For the purpose of environment calibration, the International Maritime Organization’s (IMO) correction method is applied to raw measurement data (IMO, 2002). The correction requires a record of a turning manoeuvre, of at least a 720° change of heading. From the data, two half circles can be obtained after the ship’s heading change of 180° from the beginning of the test. The local current velocity Vi can be defined by the two corresponding positions and from the half circles as equation (11):
From the set of local velocity, the estimated current velocity can be calculated by equation (12):
The magnitude of the current velocity can be calculated from the equation (13):
The final corrected trajectories of the measurement data can be obtained from equation (14):
where X(t) is the measured position vector and is the corrected position vector of the ship and . Fig. 4 shows the result of the environmental calibration for the measurement data.
The number of target variables is an important factor to determine calculation time in an optimization. Unimportant target variables may cause confusion in the decision of the optimization algorithm. Thus, is it important to conduct a sensitivity analysis for all hydrodynamic coefficients prior to optimization. Table 1 presents the results of the sensitivity analysis, from a previous study (Kim et al., 2017).
3.3 Optimization conditions
Table 2 presents conditions for the whole optimization process which is composed of three steps. Different from steps 2 and e, step 1 does an optimization without constraints. A manoeuvre characteristic of straight motion, distance by LBP, is chosen for the objective function of step 1. For steps 2 and 3, an objective function is the trajectory difference between the benchmark and a simulation result. A nonlinear equality function is the difference in the first overshoot angle of the zigzag manoeuvre and the tactical diameter of the turning manoeuvre, respectively. The conditions are applied for both the Interior Point and SQP algorithms.
4 Comparison of algorithms
4.1 Optimization results
Table 3 presents the optimized hydrodynamic coefficients for all steps. The coefficients optimized by both algorithms differ from each other, but the difference is not very large. This is also shown in Table 4, which compares the manoeuvre characteristics and number of iterations for the benchmark, with simulation results using Clark estimation and simulation results optimizations. The number of iterations for both algorithms is the same, and derived manoeuvre characteristics are similar to each other.
Fig. 5 through Fig. 8 show the comparison of trajectories and heading changes for the straight motion, zigzag manoeuvre and turning manoeuvre, respectively. The results for both algorithms are close to the benchmark, compared to the result of the Clarke estimation. In the case of the turning manoeuvre, the advance and tactical diameter for the results of both algorithms are to with the benchmark, but the steady-state diameter is bigger than the benchmark. This can be solved by applying the constraints. Fig. 6, 7
4.2 Review of results
The comparison of simulations shows that there is only a minor difference between the results using the two algorithms. The results of the straight motion and the zigzag manoeuvre are almost the same. For the case of the turning manoeuvre, the advance and tactical diameter of the interior point method are fit to the benchmark more than the ones for the SQP. However, as shown in Fig. 8, the diameter of the steady-state circle for both methods is greater than that of the benchmark. It is supposed that this result is caused by the optimization settings: improper design variables, insufficient constraints and so on. An enhanced sensitivity analysis for design variables and strict constraints can be helpful for satisfactory and accurate optimization. This should be considered in future studies with additional sea trial measurement data.
5 Conclusion
This paper compared two algorithms for efficient optimization of the hydrodynamic coefficients for ship modelling. The mathematical optimization and its algorithm is the main idea of the system identification method and its result is dependent on choosing a proper algorithm. For this, an interior point and an SQP algorithm are chosen for comparison. The study can be summarized as follows:
-
1) The optimization is carried out under three steps, and each step has a different benchmark manoeuvre. Prior to the conduct of the optimization, an environmental correction for the measurements and a sensitivity analysis for the coefficients is carried out.
-
2) Simulation results using optimized coefficients for both algorithms are almost the same as each other, and they are close to the benchmark compared to the results of the Clarke estimation.
Based on the results of this study, it is shown that both algorithms are suitable to apply in the optimization. However, comparison is not satisfactory because their results are similar despite the difference in finding the optimal solution. A stricter sensitivity analysis for design variables and constraints should be considered in future studies for acceptable optimization results. Also, preparation of additional measurements with various type of ships is required for the reliability of the optimization.