Optimization Approach for a Catamaran Hull Using CAESES and STAR-CCM+

Article information

J. Ocean Eng. Technol. 2020;34(4):272-276
Publication date (electronic) : 2020 August 31
doi : https://doi.org/10.26748/KSOE.2019.058
1Graduate Student, Interdisciplinary Program of Marine Convergence Design, Pukyong National University, Busan, Korea
2Professor, Department of Naval Architecture and Marine Systems Engineering, Pukyong National University, Busan, Korea
Corresponding author Dong-Joon Kim: +82-51-629-6614, djkim@pknu.ac.kr
Received 2019 July 12; Revised 2020 June 22; Accepted 2020 July 16.

Abstract

This paper presents an optimization process for a catamaran hull form. The entire optimization process was managed using the CAD-CFD integration platform CAESES. The resistance of the demi-hull was simulated in calm water using the CFD solver STAR-CCM+, and an inviscid fluid model was used to reduce the computing time. The Free-Form Deformation (FFD) method was used to make local changes in the bulbous bow. For the optimization of the bulbous bow, the Non-dominated Sorting Genetic Algorithm (NSGA)-II was applied, and the optimization variables were the length, breadth, and angle between the bulbous bow and the base line. The Lackenby method was used for global variation of the bow of the hull. Nine hull forms were generated by moving the center of buoyancy while keeping the displacement constant. The optimum bow part was selected by comparing the resistance of the forms. After obtaining the optimum demi-hull, the distance between two demi-hulls was optimized. The results show that the proposed optimization sequence can be used to reduce the resistance of a catamaran in calm water.

1. Introduction

Simulation-driven design is one of the most important methods for ship optimization. With the development of computer-aided design technology, research on hull shape optimization through computer simulation has gradually been applied to design more energy-efficient and environmentally friendly ships. In the preliminary design stage, it is very important to optimize the hydrodynamic characteristics of the hull form. As computers speed up and memory grows, researchers are experimenting more with Computer Aided Design (CAD) and simulation (CFD) methods.

Due to the complex geometry shape of a ship's hull, it is difficult to use numerical methods to describe it. Therefore, researchers often choose to modify the hull form by making changes to a basic design. Lackenby (1950) developed a method to modify a hull by controlling the position of the center of buoyancy and shifting the section curves. Since then, the Lackenby method has been widely used in hull modification.

In recent years, there have mainly been two popular ways of modifying a bulbous bow geometry: parametric modeling and the Free-Form Deformation (FFD) method. Chrismianto and Kim (2014) used a Cubic Bezier curve and curve-plan intersection methods to generate a parametric bulbous bow. Luo and Lan (2017) used a B-Spline curve and NURBS curve to generate a parametric bulbous bow in the CAD-CFD integration platform CAESES. Plug-in software called Grasshopper was used to generate a parametric bulbous bow from a few vertexes and NURBS curves.

The wave-making resistance of a ship hull depends largely on the bow part (the area between the stem and mid-ship). It is efficient to optimize the bow part of a ship to reduce the wave-making resistance. The bulbous bow and the hull between the bulbous bow and mid-ship are two main parameters to optimize.

Different governing equations are used in a CFD solver to predict a ship hull’s hydrodynamic performance. One of the most popular methods is the Reynolds averaged Navier-Stokes (RANS) method. Zhang et al. (2018) used the RANS method to calculate the total resistance in an optimization framework (Park et al., 2019; Kim et al., 2019). Usually, hundreds of simulations are carried out in an optimization process, and the hydrodynamic performance prediction can be quite a time consuming. Researchers have tried different ways to reduce the computing time. Han et al. (2012) selected a non-linear potential flow using the Rankine panel method to predict trim and sink during a simulation. Kostas et al. (2015) used the Neumann-Kelvin formulation and the boundary element method (BEM) to simulate the wave-making resistance.

Many kinds of optimization equations have been used to find an optimal solution quickly and accurately. Zhang et al. (2018) used the Particle Swarm Optimization (PSO) algorithm to help find an optimal bulbous bow. Huang et al. (2016) implemented a “new improved Artificial Bee Colony” (NIABC) algorithm in KCS hull optimization. Gammon implemented a Multi-Objective Genetic Algorithm (MOGA) in the optimization of a fishing vessel, and the optimization equations showed high accuracy and efficiency.

In this study, the wave-making resistance of a catamaran in calm water was selected as the objective function for optimization. The length, breadth, and angle of the bulbous bow were modified by the FFD method and then optimized by the Non-dominated Sorting Genetic Algorithm (NSGA)-II. The bow part was optimized by simulating nine different hull forms that were modified by the Lackenby method and then comparing them. Finally, after obtaining the optimal demi-hull, the distance between two demi-hulls was optimized.

2. Demi-hull Geometry Modification Method

2.1 Bulbous Bow Modification

The FFD and Lackenby methods were applied to modify the hull geometry. The bulbous bow area begins from station No. 20, and it was modified in three dimensions: the length, breadth, and its angle with the baseline. Fig. 1 shows how the FFD method is applied to the bulbous bow. The bulbous bow shape is modified by a control box with certain control points on it (Tomas and Scott, 1986). The bulbous bow angle α is defined in Fig. 1c as the angle between the base line and the rotated control box. The clockwise direction is defined as positive, and the anticlockwise direction is defined as negative.

Fig. 1.

The FFD method for the length (a), breadth (b) and angle (c) of the bulbous bow

Fig. 2 shows the length of the bulbous bow modified under the constraint of −0.02LOA ≤ ΔL ≤ 0.02LOA, where ΔL is the change of length of the bulbous bow, and LOA is the overall length. Fig. 3 shows the breadth of the bow modified under the constraint of 0.8B ≤ B’ ≤ 1.3B, where B is the original breadth of the bow, and B’ is the breadth of the new bow. Fig. 4 shows the angle of the bow modified under the constraint of −9° ≤ α ≤ 3°, where α is the difference from the original angle.

Fig. 2.

Right view of the bulbous bow at (a) ΔL of −0.02LOA, (b) original ΔL, and (c) ΔL of 0.02LOA

Fig. 3.

Front view (a) and Top view (b) of the bulbous bow at a breadth of 0.8B and front view (c) and top view (d) at a breadth of 1.3B

Fig. 4.

Right view of the bulbous bow at an angle of −9° (left) and an angle of 3° (right)

2.2 Bow Part Modification

The bow part is the part of the hull between the bulbous bow and mid-ship. The Lackenby method was used for the bow part modification. The longitudinal center of buoyancy (LCB) was changed slightly while the displacement was kept constant. The change in LCB (ΔLCB) was set as a design variable. Fig. 5 shows the bow part modified by moving the LCB by ±0.4%. The inflection point of the original sectional area curve was located around station 16, so this station was set to be the modification center.

Fig. 5.

Sectional area curve of bow part with LCB moved forward by 0.4% (blue) and backward by 0.4% (red)

3. Numerical Calculation by STAR-CCM+

A high-speed catamaran was selected as an example for optimization, and its main data are shown in Table 1. The high-speed catamaran travels with a high Froude number (i.e., Fn > 0.3). In this case, the viscous resistance is usually a smaller proportion of the total resistance, and the contribution of the bow part to the viscous resistance is negligible. Therefore, optimizing the resistance in an inviscid fluid is an effective way to optimize the performance of the high-speed ship. STAR-CCM+ was selected to perform the numerical simulation and evaluation, and the fluid model was set as an inviscid fluid model.

Main data of the demi-hull

To reduce the computing time and ensure the accuracy of the wave pattern of the free surface, the mesh was set to be relatively thin in the Z direction and wide in the X and Y directions, as in Fig. 6. A total of 460,000 mesh cells were generated, and the mesh is shown in Fig. 6.

Fig. 6.

Mesh distribution for the simulation

4. Software Integration and Demi-hull Optimization

4.1 Bulbous Bow Optimization

The resistance was simulated in STAR-CCM+ and fed back to CAESES by coupling the two programs. The FFD method was used to modify the bulbous bow shape with the design variables: the length, breadth, and angle. NSGA-II was then used to obtain the optimal bulbous bow.

Fig. 7 shows a flow chart of the NSGA-II process for the design variables of the bulbous bow. Firstly, 27 initial design plans were generated by changing the bulbous bow length, breadth, and angle defined by the FFD method as the first generation. 27 simulations were then carried out in STAR-CCM+ to obtain the resistance results. The plans that performed well were selected to mutate and crossover to obtain the next generation.. New generations of hull plans were simulated and selected again and again until the resistance results converged.

Fig. 7.

NSGA-II algorithm process in CASES

Fig. 8 shows the resistance results of the inviscid fluid after 5 generations of different bulbous bow plans. The results are from 135 simulations carried out by STAR-CCM+. The noted point represents the optimal solution. Fig. 9 shows the evolution of the length, breadth, and angle. It can be seen that the length converges to 1.018, the breadth converges to 1.14, and the angle converges to 2.64°.

Fig. 8.

Resistance of inviscid fluid after 5 generations

Fig. 9.

Design variables evolution in 5 generations

4.2 Bow Part Optimization

The bow part was modified by using the Lackenby method to shift the section curve, and the design variable was ΔLCB. A total of 9 different plans from the Lackenby method were computed, and the optimized demi-hull form was obtained. Fig. 10 shows a flow chart of choosing the optimal demi-hull.

Fig. 10.

Process of choosing the optimal demi-hull

Fig. 11 shows the simulation results of the 9 different hull forms from the Lackenby method. The results show that the resistance was reduced by about 1% when the longitudinal center of buoyancy was moved backward by 0.3%.

Fig. 11.

Resistance for different bow part plans at Fn 0.73

The wave pattern comparison is shown in Fig. 12. The 3 main differences are pointed out with 3 black squares. The wave pattern at the bow of the optimal demi-hull is smaller, and the wave pattern in the far field of the optimal demi-hull is lower than the original one.

Fig. 12.

Wave pattern comparison between original and optimal hull at design speed.

5. Optimization of Distance Between Demi-hulls

According to Millward (1992), the effect on the resistance coefficient is separated into three parts. The effect of the distance has no rules to follow when Fn is in the range of around 0.2 to 0.4. In the Fn range of 0.4 to 0.7, the distance between the demi-hulls has a positive correlation with the total wave-making resistance. When the distance between the demi-hulls is larger, the negative interference of the two demi-hulls is weaker. When Fn is beyond 0.7, the 2 lines of the resistance coefficient begin to cross, which means the effect of the distance also has no rules to follow. However, we do know that there is a potential optimal distance.

With the design speed of the catamaran (Fn=0.73), 19 different cases of separation between the two demi-hulls were generated and simulated. The results are shown in Fig. 13. The resistance of the catamaran is smallest when S/L is 0.42. This means the optimal distance of the two demi-hulls is 8.4 m. S is the distance between the two demi-hulls, and L is the length of the catamaran. The results show that the distance between the demi-hulls of the catamaran does have an optimal value in this case.

Fig. 13.

Resistance of the catamaran for different separation plans

6. Conclusion

The hull form of a high-speed catamaran was optimized by coupling the software CAESES and STAR-CCM+. The resistance simulation was carried out using the RANS method, and the fluid model was set to inviscid to reduce the computing time. The bulbous bow of the demi-hull was optimized by a genetic algorithm, and the forms were generated by the FFD method for various lengths, breadths, and angles. The bow part between the bulbous bow and the mid-ship was then optimized by simulating 9 different hull forms that were modified by the Lackenby method while keeping the displacement of the demi-hull constant. The design variable was ΔLCB. The distance between the two demi-hulls was then optimized by simulating 19 different separation cases.

The wave-making resistance of the optimal demi-hull was reduced by 6.2% compared to the original demi-hull. The total resistance of the catamaran had optimal performance when the distance between the two demi-hulls was 8.4 m (S/L=0.42). The results showed that this optimization loop is feasible and efficient. The NSGA-II algorithm was used for only the bulbous bow optimization, and the bow part was optimized independently. However, the bulbous bow and the bow part affect each other. Future work will focus on optimizing the bulbous bow and bow part of the hull together to find the best combination of the two parts.

Acknowledgements

This work was supported by a Research Grant of Pukyong National University (2020Year).

Notes

Author ORCIDs

References

Chrismianto D, Kim DJ. 2014;Parametric Bulbous Bow Design Using the Cubic Bezier Curve and Curve-plane Intersection Method for the Minimization of Ship Resistance in CFD. Journal of Marine Science and Technology 19(4):479–492. https://doi.org/10.1007/s00773-014-0278-x.
Han SH, Lee YS, Choi YB. 2012;Hydrodynamic Hull form Optimization Using Parametric Models. Journal of Marine Science and Technology 17:1–17. https://doi.org/10.1007/s00773-011-0148-8.
Huang FX, Wang LJ, Yang C. 2016;A New Improved Artificial Bee Colony Algorithm for Ship Hull form Optimization. Engineering Optimization 48(4):672–686. https://doi.org/10.1080/0305215X.2015.1031660.
Kim YC, Kim Y, Kim J, Kim KS. 2019;Application of the Overset Grid Scheme (Suggar++) for Flow Analysis around a Ship. Journal of the Society of Naval Architects of Korea 56(1):47–57. https://doi.org/10.3744/SNAK.2019.56.1.047.
Kostas KV, Ginnis AI, Politis CG, Kaklis PD. 2015;Ship-hull Shape Optimization with a T-spline Based BEM Iso-geometric Solver. Computer Methods in Applied Mechanics and Engineering 284:611–622. https://doi.org/10.1016/j.cma.2014.10.030.
Luo W, Lan L. 2017;Design Optimization of the Lines of the Bulbous Bow of a Hull Based on Parametric Modeling and Computational Fluid Dynamics Calculation. Mathematical and Computational Applications 22(4):1–12. https://doi.org/10.3390/mca22010004.
Lackenby H. 1950;On the Systematic Geometrical Variation of Ship Forms. RINA Transactions 92:289–309.
Millward A. 1992;The Effect of Hull Separation and Restricted Water Depth on Catamaran Resistance. Transactions of the Royal Institute of Naval Architects 134:341–349.
Park K, Kim DJ, Kim SY, Rhee SH. 2019;A Study on the Resistance Performance and Flow Pattern of High Speed Planing Hull using CFD. Journal of the Society of Naval Architects of Korea 56(1):23–33. https://doi.org/10.3744/SNAK.2019.56.1.023.
Thomas WS, Scott RP. 1986;Free-form Deformation of Solid Geometric Models. Proceedings of SIGGRAPH - Special Interest Group on GRAPHics and Interactive Techniques 20:151–159.
Zhang SL, Zhang BJ, Tezdogan T, Xu L, Lai YY. 2018;Computational Fluid Dynamics-based Hull form Optimization Using Approximation Method. Engineering Applications of Computational Fluid Mechanics 12(1):74–88. https://doi.org/10.1080/19942060.2017.1343751.

Article information Continued

Author ORCIDs

Author name ORCID
Zhang Yongxing 0000-0003-2610-6594
Kim, Dong-Joon 0000-0003-0730-8721

Fig. 1.

The FFD method for the length (a), breadth (b) and angle (c) of the bulbous bow

Fig. 2.

Right view of the bulbous bow at (a) ΔL of −0.02LOA, (b) original ΔL, and (c) ΔL of 0.02LOA

Fig. 3.

Front view (a) and Top view (b) of the bulbous bow at a breadth of 0.8B and front view (c) and top view (d) at a breadth of 1.3B

Fig. 4.

Right view of the bulbous bow at an angle of −9° (left) and an angle of 3° (right)

Fig. 5.

Sectional area curve of bow part with LCB moved forward by 0.4% (blue) and backward by 0.4% (red)

Fig. 6.

Mesh distribution for the simulation

Fig. 7.

NSGA-II algorithm process in CASES

Fig. 8.

Resistance of inviscid fluid after 5 generations

Fig. 9.

Design variables evolution in 5 generations

Fig. 10.

Process of choosing the optimal demi-hull

Fig. 11.

Resistance for different bow part plans at Fn 0.73

Fig. 12.

Wave pattern comparison between original and optimal hull at design speed.

Fig. 13.

Resistance of the catamaran for different separation plans

Table 1.

Main data of the demi-hull

Item Value
LOA 21.7 (m)
LPP 20.0 (m)
B 2.5 (m)
D 3.2 (m)
d 1.6 (m)
V 10.28 (m/s)
Fn 0.73