Design Optimization of the Lines of the Bulbous Bow of a Hull Based on Parametric Modeling and Computational Fluid Dynamics Calculation

To reduce the ship wave-making resistance, the lines of the bulbous bow of a hull are optimized by an automatic optimization platform at the ship design stage. Parametric modeling was applied to the hull by using non-uniform rational basis spline (NURBS). The Rankine-source panel method was used to calculate the wave-making resistance. A hybrid optimization strategy was applied to achieve the optimization goal. A Ro-Ro ship was taken as an example to illustrate the optimization method adopted, with the objective to minimize the wave-making resistance. The optimization results show that wave-making resistance obviously reduces and the wave-shape of the near bow becomes gentle after the lines of the bulbous bow of the hull are optimized, which demonstrates the validity of the proposed optimization design strategy.


Introduction
Design of ship hull forms is one of the most important contents of the general ship design.How to reduce the ship resistance by means of well-designed forms is always the concern for ship designers and/or researchers.The parent ship based design method predominates in the field of ship design, by which the designed hull can inherit some excellent performance from the parent ship.For a long time, designer's experience played an important role in obtaining the desired hull forms.With the development of computer science and technology, the quality of ship design has been improved considerably by introducing some computational numerical calculation methods.For example, computational fluid dynamics (CFD) based numerical simulation can be used to evaluate the hydrodynamic performances of a ship after the hull forms are designed.By modifying the hull forms and repeatedly performing CFD simulation, the satisfactory hull forms can be obtained.To improve the design efficiency, computer aided design (CAD) was proposed to combine with CFD in the optimization design of hull forms [1].During the past decade, CAD/CFD integrated hull-form design systems have been developed and successfully applied to the drag reduction of ships [2].For example, Percival et al. combined a simple CFD tool and B-Spline method to the optimization design of a Wigley ship [3].Yusuke et al. presented the bulbous bow optimization of a container ship based on the naval architectural package (NAPA) system and a Rankine-source panel method [1].Tahara et al. applied NAPA and a RANS (Reynolds-averaged Navier-Stokes) solver to the optimization design of a container ship [4].Abt and Harries presented the generic, multi-objective ship design optimization based on software of the Ship Design Laboratory of the National Technical University of Athens (NTUA-SDL) in which well-established naval architectural and optimization software packages are integrated [5].Harries and Abt distinguished three types of coupling of CAD and CFD [6].Couser et al. explored the multi-dimensional design-space using first-principle methods within the FRIENDSHIP framework [7].Brizzolara and Vernengo proposed an integration of a parametric generation module, a multi-objective optimization algorithm and a CFD solver in the optimization of a new unmanned surface vehicle (USV) and small waterplane area twin hull (SWATH) [8].Biliotti et al. presented a fully automatic optimization chain to the optimization of fast naval vessels by adopting the ModeFrontier optimization environment to the FRIENDSHIP framework, the CFD codes and a multi-objective genetic algorithm [9].Ginnis et al. employed the NURBS (non-uniform rational basis spline) method and Kelvin-source method to the optimization of a container ship [10].Brenner et al. used the FRIENDSHIP-SHIPFLOW optimization scheme to retrofit the ships in operation [11].Zakerdoost et al. used the NURBS method and potential theory to the optimization design of a Wigley ship [12].Vernengo et al. presented the optimization of a Semi-SWATH by a fully parametric model of the unconventional hull form, a 3D linear Rankine-sources panel method, and a multiobjective global convergence genetic algorithm [13].Combined with the optimization of a high speed round bilge monohull, Brizzolara et al. addressed the significance of parametric hull form definition by comparing it with Free Form Deformation technique [14].Ang et al. proposed the Free Form Deformation and potential theory to the optimization design of a heavy-lift and pipe-lay vessel [15].Wang proposed the NURBS method and CFD calculation based on Neumann-Michell theory to the optimization design of several life ships [16].
Besides the parametric modeling of hull forms and calculation of hydrodynamic performances, optimization is another important factor in the hull-form design system.To guarantee the optimum design and the solving efficiency, an appropriate optimization algorithm is vital to the optimizer.During the last several decades, many optimization algorithms have been proposed such as successive quadratic programming (SQP) [17], genetic algorithm (GA) [18], simulated annealing (SA) [19], particle swarm optimization (PSO) [20], infeasibility driven evolutionary algorithm (IDEA) [21], and quasi-Newton method [10], among others.Generally, these algorithms can be categorized into gradient-based and derivative-free methods [22].For the gradient-based optimization algorithms (e.g., SQP and quasi-Newton approaches), the interest is the fast convergence; however, the limitation is the locally optimal solution.For the derivative-free optimization algorithms (e.g., GA, SA, PSO, and IDEA approaches), the advantage is the global optimization; however, the convergence rate cannot be guaranteed.Since each optimization algorithm has its advantages and disadvantages, a hybrid optimization scheme should be preferable.In this paper, we perform the combination of a global optimization method (Ensemble Investigation) and a local optimization method (T-Search algorithm) to the optimization design of hull forms of a surface ship.The design optimization is carried out in the SHIPFLOW-CAESES (former FRIENDSHIP) DESIGN environment.Based on the automatic optimization platform, the parametric modeling and the CFD calculation are integrated.The NURBS method is employed to the parametric modeling while the Rankine-source panel method is applied to calculate the wave-making resistance.The minimization of the wave-making resistance is taken as the optimization objective while the design variables are form parameters of the lines of the bulbous bow.
The rest of the paper is organized as follows.In Section 2, the parametric modeling method is introduced; in Section 3, the CFD method to calculate the wave-making resistance is described; in Section 4, the hybrid optimization algorithm is addressed; in Section 5, the entire optimization scheme is described; in Section 6, an example is studied to confirm the validity of the optimization strategy proposed in the study; and the final section is the concluding remarks.

Parametric Modeling of Hull Forms
In the geometry modeling of a hull, it is often required to obtain the forms rapidly.In addition, a large number of forms should be analyzed so that the optimization can be efficiently performed.
In this paper, the NURBS method is employed to meet the mentioned requirements.Based on this method, the high-quality of the lines can be guaranteed with fairly less design variables than conventional design methodology.This method provides the best way to establish the relationship between parameters and hull forms and makes sure that the design variation among the entire optimization process is effective and feasible.The parametric model is generated by CAESES-Modeler.The process of curves generation starts with a set of given data elements that are used to define the complete parametric curve.Within the modeler, one can define a set of longitudinal lines and express them as basic curves to describe the ship hull by means of differential, integral and topological information.Transversal basic curves are able to be formed as soon as the longitudinal basic curves are defined (by fair B-Spline parametric curves adopted in the study).Then, a series of surfaces can be set up by corresponding transversal sections that stem from the parametric hull model.Generally, the design of the fully parametric hull surfaces makes it possible to transform the hull form efficiently and effectively.Moreover, all of the hull topologies, the hull sections and the whole hull surface can be generated using the basic curves specifically and precisely.Usually, the hull forms are described by feature curves [23], e.g., the design waterline, flat of side, flat of bottom, etc., which are listed in Table 1 and some of the curves are shown in Figure 1.conventional design methodology.This method provides the best way to establish the relationship between parameters and hull forms and makes sure that the design variation among the entire optimization process is effective and feasible.The parametric model is generated by CAESES-Modeler.The process of curves generation starts with a set of given data elements that are used to define the complete parametric curve.Within the modeler, one can define a set of longitudinal lines and express them as basic curves to describe the ship hull by means of differential, integral and topological information.Transversal basic curves are able to be formed as soon as the longitudinal basic curves are defined (by fair B-Spline parametric curves adopted in the study).Then, a series of surfaces can be set up by corresponding transversal sections that stem from the parametric hull model.Generally, the design of the fully parametric hull surfaces makes it possible to transform the hull form efficiently and effectively.Moreover, all of the hull topologies, the hull sections and the whole hull surface can be generated using the basic curves specifically and precisely.Usually, the hull forms are described by feature curves [23], e.g., the design waterline, flat of side, flat of bottom, etc., which are listed in Table 1 and some of the curves are shown in Figure 1.For the bulbous bow, the basic feature curves include: the top longtitudinal section line, the low longtitudinal section line, the half beam line, and the height line of the maximal breadth, as shown in Figure 2. The shape of the bulbous bow can be modified by changing these basic feature curves that can be achieved by controlling design variables.Figure 3 presents an example of the variation of the shape of a bulbous bow by modifying these feature curves.For the bulbous bow, the basic feature curves include: the top longtitudinal section line, the low longtitudinal section line, the half beam line, and the height line of the maximal breadth, as shown in Figure 2. The shape of the bulbous bow can be modified by changing these basic feature curves that can be achieved by controlling design variables.Figure 3 presents an example of the variation of the shape of a bulbous bow by modifying these feature curves.Using the NURBS method, curves can be expressed as [24]: where n is the number vertex points, wi is the weight, di is the control point, and Ni,k(t) is the B-Spline basis function.For the hull forms, a NURBS surface is a bi-variate vector-valued piecewise rational function as [24]: where Pi,j forms a bidirectional control net [12].Ni,p(u) and Nj,q(v) are B-Spline basis functions in the direction of The hull forms can be expressed parametrically by using the form parameters to constitute the feature curves and hull section frames.The general procedure to generate fairing surfaces of the hull Using the NURBS method, curves can be expressed as [24]: where n is the number vertex points, wi is the weight, di is the control point, and Ni,k(t) is the B-Spline basis function.For the hull forms, a NURBS surface is a bi-variate vector-valued piecewise rational function as [24]: where Pi,j forms a bidirectional control net [12].Ni,p(u) and Nj,q(v) are B-Spline basis functions in the direction of The hull forms can be expressed parametrically by using the form parameters to constitute the feature curves and hull section frames.The general procedure to generate fairing surfaces of the hull Using the NURBS method, curves can be expressed as [24]: where n is the number vertex points, w i is the weight, d i is the control point, and N i,k (t) is the B-Spline basis function.For the hull forms, a NURBS surface is a bi-variate vector-valued piecewise rational function as [24]: where P i,j forms a bidirectional control net [12].N i,p (u) and N j,q (v) are B-Spline basis functions in the direction of The hull forms can be expressed parametrically by using the form parameters to constitute the feature curves and hull section frames.The general procedure to generate fairing surfaces of the hull form observes four steps, as presented in Figure 4. First, the form parameters should be suitably selected.Second, the parametric design of the longitudinal feature curves is performed.Third, transversal feature curves can be similarly generated and a set of surfaces can be derived by means of interpolation.Finally, fairing form variations can be obtained [23].It should be noted that the number and the range of the form parameters should be selected carefully to prevent any distortion of the ship.In the study, eight form parameters are selected for the generation of the bulbous bow.More details about these parameters are presented in the Example Study section.form observes four steps, as presented in Figure 4. First, the form parameters should be suitably selected.Second, the parametric design of the longitudinal feature curves is performed.Third, transversal feature curves can be similarly generated and a set of surfaces can be derived by means of interpolation.Finally, fairing form variations can be obtained [23].It should be noted that the number and the range of the form parameters should be selected carefully to prevent any distortion of the ship.In the study, eight form parameters are selected for the generation of the bulbous bow.More details about these parameters are presented in the Example Study section.

Calculation of Wave-Making Resistance
In the study, the optimization objective is to minimize the wave-making resistance.Suppose the fluid is inviscid and incompressible, and Rankine-source panel method is employed.
For a ship traveling with a steady forward speed U in the calm water of infinite depth, the wavemaking velocity potential ϕ(x,y,z) satisfies the following equations: in the fluid domain: on the free surface: on the free surface: ( ) on the hull surface: where , g is the acceleration of gravity, and ξ is the free surface elevation.
The velocity potential at the field point P can be calculated as [25] ( ) ( ) ( ) ( ) ( ) ( ) ( ) , , where Q and Q' are source points, σ(Q) is the source strength, and r and r' are the distances between the field point and source points.When the ship is sailing, the wave is steep and the floating state cannot be neglected.Therefore, the nonlinear free surface boundary condition should be used.The method of Taylor series expansion is applied to Equations ( 4) and ( 5) by making use of the approximate solution z = Z and ϕ = φ.Then, the boundary condition on the approximate free surface can be obtained as:

Calculation of Wave-Making Resistance
In the study, the optimization objective is to minimize the wave-making resistance.Suppose the fluid is inviscid and incompressible, and Rankine-source panel method is employed.
For a ship traveling with a steady forward speed U in the calm water of infinite depth, the wave-making velocity potential ϕ(x,y,z) satisfies the following equations: in the fluid domain : on the free surface : = on the free surface : on the hull surface : at infinity : ∇φ = (0, 0, 0) where U s = (U, 0, 0), g is the acceleration of gravity, and ξ is the free surface elevation.
The velocity potential at the field point P can be calculated as [25] φ where Q and Q' are source points, σ(Q) is the source strength, and r and r' are the distances between the field point and source points.When the ship is sailing, the wave is steep and the floating state cannot be neglected.Therefore, the nonlinear free surface boundary condition should be used.The method of Taylor series expansion is applied to Equations ( 4) and ( 5) by making use of the approximate solution z = Z and ϕ = φ.Then, the boundary condition on the approximate free surface can be obtained as: where The boundary condition on the ship surface is: Based on Bernoulli equation, the pressure can be calculated as where U is the ship velocity, and ρ is the fluid density.By integrating the pressure along the wetted boundary surface, one has the hydrodynamic force on the hull as Thus, the wave-making resistance and its nondimensional form can be obtained as and respectively, where S 0 is the wetted surface area in calm waters.C w is also named as the coefficient of wave-making resistance.

Optimization Algorithms
To guarantee the globally optimal solution and solving efficiency as well, a hybrid algorithm is applied in the study.Ensemble Investigation is employed to obtain the globally optimal solution.This algorithm allows permuting a set of design variables and generating DOE (design of experiments) tables, which is very effective to gather information about the optimization problem investigated in the whole design space.Generated DOE tables are useful to detect the trends of the optimization variables with regard to the optimization objectives, irrespective of constraints [26].This optimization approach is performed by creating several different versions of a proposed new design, and then selecting the option that best fits the design goals.The generation of a systematic series of design options can be completed soon.These options are tested against design goals-for example, the hydrodynamic performance.Computational techniques such as regression analysis can be used to determine the "best" option.
To improve the efficiency of solution-searching, the T-Search algorithm is incorporated.This algorithm is a method of feasible direction.It uses a direct search until some constraint gets violated.At the expense of a little extra effort, T-Search is characterized by efficient operation within the admissible region with satisfactory performance along the boundary of that region [27].The other two features of the algorithm are the allowance made for small perturbations of the base point and the flexible directions for exploration in the tangent hyperspace [28].It consists of exploratory moves that start from a so-called base point along the variable axes followed by global moves in the descent search direction found in successful exploratory moves.If a constraint bound is approached, then it returns to the last admissible point and performs a unidirectional optimization along the hyperplane tangent to the active constraint either to keep the search in the feasible domain or to bring it back to the feasible domain.T-Search is able to find the local minimum of an arbitrary, explicitly-stated function of more than one variable, subject to arbitrary, non-linear constraints.

Optimization Goal and Design Variables
The optimization objective in the study is to minimize the wave-making resistance under designed ship velocity, i.e., min{C w }.The constraints are set to be subject to (1) |∆ 0 − ∆ p |/∆ 0 ≤ 1% where ∆ 0 is the initial displacement, and ∆ p is the displacement of optimized ship; (2) |LCB 0 − LCB p |/LCB 0 ≤ 1%, where LCB 0 is the initial longitudal coordiante of the boyance center while LCB p is the longitudal coordiante of boyance center of the optimized ship; and (3) the ship length, breadth and draft are constant during the optimization process.
Since the lines of the bulbous bow are vital to the wave-making resistance, the design variables are selected as form control parameters of bulbous bow, including the distance between the center of the bulb bow and the front end of the bulb bow (D cf ); the height of the front end of the bulb bow (H f ); the contour fullness of the top longtitudinal section of bulb bow (F tl ) and its coefficient (C ftl ); the contour fullness of the low longtitudinal section of bulb bow (F ll ) and its coefficient (C fll ); the contour fullness at the half beam of the bulb bow (F hb ); and half of the maximal breadth of the bulb bow (B h ).These form control parameters are used to generate the basic feature curves of the bulbous bow, as shown in Figure 2.

Optimization Framework
In the optimization framework, there are mainly three modules: the parametric modeling, CFD calculation and lines modification.The optimization operation is performed in the SHIPFLOW-CAESES DESIGN environment.First, NURBS based parametric modeling is performed in the CAESES and the lines offsets of the given hull can be obtained.Then, the files of the lines of offsets are automatically imported into the SHIPFLOW so that the CFD calculation can be carried out.Afterwards, the CFD results are imported to CAESES as comparison with the optimization objective.If the optimization goal is achieved, the lines are exported as the optimal lines; otherwise, the design variables are modified based on the hybrid optimization algorithm, and the optimization operation is repeated until the optimization goal is achieved.For general procedure, the optimization design is presented in Figure 5.
The optimization objective in the study is to minimize the wave-making resistance under designed ship velocity, i.e., min{Cw}.The constraints are set to be subject to (1) |Δ0 − Δp|/Δ0 ≤ 1% where Δ0 is the initial displacement, and Δp is the displacement of optimized ship; (2) |LCB0 − LCBp|/LCB0 ≤ 1%, where LCB0 is the initial longitudal coordiante of the boyance center while LCBp is the longitudal coordiante of boyance center of the optimized ship; and (3) the ship length, breadth and draft are constant during the optimization process.
Since the lines of the bulbous bow are vital to the wave-making resistance, the design variables are selected as form control parameters of bulbous bow, including the distance between the center of the bulb bow and the front end of the bulb bow (Dcf); the height of the front end of the bulb bow (Hf); the contour fullness of the top longtitudinal section of bulb bow (Ftl) and its coefficient (Cftl); the contour fullness of the low longtitudinal section of bulb bow (Fll) and its coefficient (Cfll); the contour fullness at the half beam of the bulb bow (Fhb); and half of the maximal breadth of the bulb bow (Bh).These form control parameters are used to generate the basic feature curves of the bulbous bow, as shown in Figure 2.

Optimization Framework
In the optimization framework, there are mainly three modules: the parametric modeling, CFD calculation and lines modification.The optimization operation is performed in the SHIPFLOW-CAESES DESIGN environment.First, NURBS based parametric modeling is performed in the CAESES and the lines offsets of the given hull can be obtained.Then, the files of the lines of offsets are automatically imported into the SHIPFLOW so that the CFD calculation can be carried out.Afterwards, the CFD results are imported to CAESES as comparison with the optimization objective.If the optimization goal is achieved, the lines are exported as the optimal lines; otherwise, the design variables are modified based on the hybrid optimization algorithm, and the optimization operation is repeated until the optimization goal is achieved.For general procedure, the optimization design is presented in Figure 5.
It is noted that a batch file is necessary to connect SHIPFLOW and CAESES applications.SHIPFLOW provides the batch file in the installation package.The batch file includes echo commands, pause commands, call commands, start commands, go-to commands, and set commands.There is a boundary of surface between the surface of the bulbous bow and the surface of the front part of hull.This boundary is maintained between optimization circles so that the bulbous bow can connect to the front part of hull smoothly.It is noted that a batch file is necessary to connect SHIPFLOW and CAESES applications.SHIPFLOW provides the batch file in the installation package.The batch file includes echo commands, pause commands, call commands, start commands, go-to commands, and set commands.There is a boundary of surface between the surface of the bulbous bow and the surface of the front part of hull.This boundary is maintained between optimization circles so that the bulbous bow can connect to the front part of hull smoothly.

Example Study
A Ro-Ro ship is investigated to verify the optimization scheme proposed in the study.The main particulars of the ship are listed in Table 2.The ship velocity is 22 knots.In terms of Froude number, one has Fr = 0.28.In the CAESES framework, the form parameters are imported to the curve generator to obtain a set of feature curves.Then, fairing surfaces can be derived automatically through a surface generator.The parameterized hull surfaces are shown in Figure 6.

Example Study
A Ro-Ro ship is investigated to verify the optimization scheme proposed in the study.The main particulars of the ship are listed in Table 2.The ship velocity is 22 knots.In terms of Froude number, one has Fr = 0.28.In the CAESES framework, the form parameters are imported to the curve generator to obtain a set of feature curves.Then, fairing surfaces can be derived automatically through a surface generator.The parameterized hull surfaces are shown in Figure 6.Based on the optimization framework shown in Figure 5, the hybrid optimization algorithm is performed.In detail, the Ensemble Investigation is employed to obtain feasible solutions in the whole design space.A set of optimized hull forms that satisfy the constraints are derived.Then, the T-Search algorithm is used within the set to determine the solution with minimal wave-making resistance.The design variables and their variations are listed in Table 3.The results of optimization objective and constraints are presented in Table 4. Dcf, the distance between the center of the bulb bow and the front end of the bulb bow; Hf, the height Based on the optimization framework shown in Figure 5, the hybrid optimization algorithm is performed.In detail, the Ensemble Investigation is employed to obtain feasible solutions in the whole design space.A set of optimized hull forms that satisfy the constraints are derived.Then, the T-Search algorithm is used within the set to determine the solution with minimal wave-making resistance.The design variables and their variations are listed in Table 3.The results of optimization objective and constraints are presented in Table 4.As can be seen from the optimization results, the wave-making resistance is reduced by modifying the design variables while the constraints are satisfied.The initial values of design variables are one set of the feasible solutions family.The variation of the design variables indicates that smaller D cf and F tl , against larger H f , C ftl , F ll , C fll , F hb , and B h help reduce the wave-making resistance.The comparison of the wave contours on free surfaces before and after optimization is presented in Figure 7.The comparison of the longitudinal wave cuts before and after optimization is given as Figure 8.The right plot in Figure 7 presents the changes of wave contours on the forward part of the ship.Similarly, in Figure 8, the right plot highlights the changes of wave cuts on the area around the bulbous bow.As can be seen from the optimization results, the wave-making resistance is reduced by modifying the design variables while the constraints are satisfied.The initial values of design variables are one set of the feasible solutions family.The variation of the design variables indicates that smaller Dcf and Ftl, against larger Hf, Cftl, Fll, Cfll, Fhb, and Bh help reduce the wave-making resistance.The comparison of the wave contours on free surfaces before and after optimization is presented in Figure 7.The comparison of the longitudinal wave cuts before and after optimization is given as Figure 8.The right plot in Figure 7 presents the changes of wave contours on the forward part of the ship.Similarly, in Figure 8, the right plot highlights the changes of wave cuts on the area around the bulbous bow.As can be seen from the wave variation, the wave numbers and the wave amplitude decrease obviously around the bow after the lines of the bulb bow are optimized, which implies a decrease of wave-making resistance.Figures 9 and 10 present the variation of the lines of the bulb bow.
In Figure 9, the left plot is the original form of the bulb bow, while the right plot is the optimized form.In Figure 10, the yellow lines represent the optimized results while the black lines are with the original form.From the comparison results, it can be indicated that a stretch of the bulb bow upwards helps reduce the wave-making resistance.As can be seen from the wave variation, the wave numbers and the wave amplitude decrease obviously around the bow after the lines of the bulb bow are optimized, which implies a decrease of wave-making resistance.Figures 9 and 10 present the variation of the lines of the bulb bow.In Figure 9, the left plot is the original form of the bulb bow, while the right plot is the optimized form.In Figure 10, the yellow lines represent the optimized results while the black lines are with the original form.From the comparison results, it can be indicated that a stretch of the bulb bow upwards helps reduce the wave-making resistance.

Conclusions
The combination of CFD and CAD improves the quality of ship hull form design.In this paper, an optimization design strategy is proposed based on SHIPFLOW-CAESES framework.The optimization results demonstrate the feasibility of using successfully different software systems in the optimization process.A fully automatic optimization integration framework is achieved by using CAESES coupled with CFD solver SHIPFLOW.Through the framework, the full parametric modeling of the hull, hydrodynamic calculation, design evaluation and shape variation are implemented effectively, which improves the quality and efficiency of ship design.The NURBS method is employed in the parametric modeling of the hull forms, while potential theory is used to calculate the wave-making resistance of the ship hull.To improve the optimization efficiency, a As can be seen from the wave variation, the wave numbers and the wave amplitude decrease obviously around the bow after the lines of the bulb bow are optimized, which implies a decrease of wave-making resistance.Figures 9 and 10 present the variation of the lines of the bulb bow.In Figure 9, the left plot is the original form of the bulb bow, while the right plot is the optimized form.In Figure 10, the yellow lines represent the optimized results while the black lines are with the original form.From the comparison results, it can be indicated that a stretch of the bulb bow upwards helps reduce the wave-making resistance.

Conclusions
The combination of CFD and CAD improves the quality of ship hull form design.In this paper, an optimization design strategy is proposed based on SHIPFLOW-CAESES framework.The optimization results demonstrate the feasibility of using successfully different software systems in the optimization process.A fully automatic optimization integration framework is achieved by using CAESES coupled with CFD solver SHIPFLOW.Through the framework, the full parametric modeling of the hull, hydrodynamic calculation, design evaluation and shape variation are implemented effectively, which improves the quality and efficiency of ship design.The NURBS method is employed in the parametric modeling of the hull forms, while potential theory is used to

Conclusions
The combination of CFD and CAD improves the quality of ship hull form design.In this paper, an optimization design strategy is proposed based on SHIPFLOW-CAESES framework.The optimization results demonstrate the feasibility of using successfully different software systems in the optimization process.A fully automatic optimization integration framework is achieved by using CAESES coupled with CFD solver SHIPFLOW.Through the framework, the full parametric modeling of the hull, hydrodynamic calculation, design evaluation and shape variation are implemented effectively, which improves the quality and efficiency of ship design.The NURBS method is employed in the parametric modeling of the hull forms, while potential theory is used to calculate the wave-making resistance of the ship hull.To improve the optimization efficiency, a hybrid optimization scheme is taken.The validity of the proposed optimization design framework is confirmed by applying the optimization framework to a Ro-Ro ship.Based on the simulation results, some recommendations can be summarized as follows: (i) the smaller D cf and F tl , against larger H f , C ftl , F ll , C fll , F hb , and B h helps reduce the wave-making resistance; and (ii) a stretch of the bulb bow upwards also helps reduce the wave-making resistance.
In this study, only the bulbous bow optimization is concerned.In the next work, the optimization of the entire hull will be studied.Moreover, total resistance will be considered in the optimization design, instead of the only wave-making resistance.In addition, another calculation method for wave resistance, i.e., the wave cut analysis will be studied.As mentioned in literature, wave cut analysis provides a more accurate and robust solution.Moreover, for a planning hull, instead of the displacement hull that is investigated in the study, wave cut analysis might be preferable.

of 12 Figure 2 .Figure 3 .
Figure 2. The basic feature curves of a bulbous bow

Figure 4 .
Figure 4. Parametric modeling of a hull form.

Figure 4 .
Figure 4. Parametric modeling of a hull form.

Figure 6 .
Figure 6.Hull surfaces of the Ro-Ro ship.

Figure 6 .
Figure 6.Hull surfaces of the Ro-Ro ship.

Figure 7 .
Figure 7.Comparison of the wave contours on free surfaces before and after optimization; upper part,-before optimization; lower part,-after optimization.Figure 7. Comparison of the wave contours on free surfaces before and after optimization; upper part,-before optimization; lower part,-after optimization.

Figure 7 .
Figure 7.Comparison of the wave contours on free surfaces before and after optimization; upper part,-before optimization; lower part,-after optimization.Figure 7. Comparison of the wave contours on free surfaces before and after optimization; upper part,-before optimization; lower part,-after optimization.

Figure 7 .
Figure 7.Comparison of the wave contours on free surfaces before and after optimization; upper part,-before optimization; lower part,-after optimization.

Figure 8 .
Figure 8.Comparison of the longitudinal wave cuts before and after optimization (y/L = 0.1).L, ship length; units in m/m.Figure 8. Comparison of the longitudinal wave cuts before and after optimization (y/L = 0.1).L, ship length; units in m/m.

Figure 8 .
Figure 8.Comparison of the longitudinal wave cuts before and after optimization (y/L = 0.1).L, ship length; units in m/m.Figure 8. Comparison of the longitudinal wave cuts before and after optimization (y/L = 0.1).L, ship length; units in m/m.

Figure 9 .
Figure 9.Comparison of the bulbous bows before and after optimization: left,-original; right,optimized.

Figure 10 .
Figure 10.Comparison of the lines before and after optimization.

Figure 9 .
Figure 9.Comparison of the bulbous bows before and after optimization: left,-original; right,-optimized.

Figure 9 .
Figure 9.Comparison of the bulbous bows before and after optimization: left,-original; right,optimized.

Figure 10 .
Figure 10.Comparison of the lines before and after optimization.

Figure 10 .
Figure 10.Comparison of the lines before and after optimization.

Table 1 .
Basic feature curves describing hull forms.

Table 1 .
Basic feature curves describing hull forms

Table 2 .
Main particulars of a Ro-Ro ship.

Table 2 .
Main particulars of a Ro-Ro ship.

Table 3 .
Optimization results of the design variables.

Table 3 .
Optimization results of the design variables.cf , the distance between the center of the bulb bow and the front end of the bulb bow; H f , the height of the front end of the bulb bow; F tl , the contour fullness of the top longtitudinal section of the bulb bow; C ftl , the coefficient of the contour fullness of the top longtitudinal section; F ll , the contour fullness of the low longtitudinal section of the bulb bow; C fll , the coefficient of the contour fullness of the low longtitudinal section; F hb , the contour fullness at the half beam of the bulb bow; B h , half of the maximal breadth of the bulb bow. D

Table 4 .
Optimization results of the objective and constraints.

Table 4 .
Fhb, the contour fullness at the half beam of the bulb bow; Bh, half of the maximal breadth of the bulb bow.Optimization results of the objective and constraints.