Procedure for Application-Oriented Optimisation of Marine Propellers

Florian Vesting 1,*, Rickard E. Bensow 1, Rikard Johansson 2, Robert Gustafsson 2 and Nicole Costa 1 1 Department of Shipping and Marine Technology, Chalmers University of Technology, 412 96 Gothenburg, Sweden; rickard.bensow@chalmers.se (R.E.B.); nicole.costa@chalmers.se (N.C.) 2 Rolls-Royce Hydrodynamic Research Centre, Rolls-Royce AB, PO Box 1010, 68129 Kristinehamn, Sweden; rikard.johansson@rolls-royce.com (R.J.); robert.gustafsson@rolls-royce.com (R.G.) * Correspondence: florian.vesting@chalmers.se; Tel.: +46-31-772-2694


Introduction-The Problem in Propeller Design
Suppliers in the competitive market of ship propellers and propulsion have experienced a change in customers' awareness for the selection criteria of an order.For merchant vessels it has predominantly been a question of investment cost to choose which supplier may deliver the propeller.The situation shifted during the last couple of years towards a situation where also the efficiency and fuel consumption became major selection criteria for the customer.With the expansion of routing in ecological sensitive areas, a recent concern is the propeller radiated noise level.However, in general, concern for operational and maintenance costs, environmental awareness and emission regulations have sharpened the perception in the merchant fleet, such that suppliers are more often forced into comparative performance races [1].
Competition between contestants is generally good since it ensures the customer to receive the product that fulfils the requirements best; and to choose from products that are obtained by different strategies and design philosophies of the competitors.However, in recent competitive situations, the designated winner is often the supplier that delivers the highest efficiency and thereby the lowest fuel consumption.The verification is commonly done by an impartial test facility, at model scale with self-propulsion tests and converted to full-scale according to standard procedures, e.g., the ITTC procedure.This competition and evaluation practice initiated a trend of sub-optimising the propeller designs towards efficiency performance and disregarding other performance characteristics [1].Such practice yield designs that are unsatisfying in operation and may ultimately end in failure of the equipment, e.g., due to cavitation erosion.
Propellers are highly adapted to the flow around the specific ship it is designed for and its operation conditions, to achieve high levels of efficiency.However, to avoid sub-optimal propellers that achieve the lowest power consumption yet ignores other performance characteristics, the designer has to consider requirements which are not only contradictory to efficiency but also have practical limitations.Since the propeller experiences varying inflow conditions while travelling through the circumferential wake, cavitation and propeller induced pressure pulses are the most evident propeller effects, contradictory to efficiency.To satisfy the shipowners' expectations and to deliver a practical design, the designer has to consider and control the cavitation, and constrain static and dynamic blade stresses, classification requirements and in case of controllable pitch propellers, hub strength and blade clearance.
Thus, propeller design is truly an art of compromising the performance and requirements in an iterative manner, with the final propeller geometry being a unique design.This is naturally a multi-objective optimisation task with constraints, intimately connected with multi-disciplinary numerical evaluation methods (e.g., for cavitation or stresses).The various requirements lead to a large number of alternatives to be studied and thereby put restrictions on the numerical analysis tools.The challenging factor for the designer is to develop a propeller that fulfils the requirements and expectations within a short time frame available in the competition races, before manufacturing of the supplier's model propeller candidate.Thus, although the evaluation of propeller performance demands highly unsteady and physically tremendously complex simulations, it is common practice to apply less accurate but faster potential methods (Mishima [2], Griffin and Kinnas [3], Lee et al. [4], Bertetta et al. [5]).
Automated propeller optimisation emerges more frequently to produce a large number of design alternatives and thereby support the designer during the design synthesis.Foeth [6] presents propeller geometry parametrisation based on database records and the multi-objective optimisation in behind condition, which applies the derived parameters.Thus, the optimisation is carried out with full interaction between the propeller changes and flow around the hull.However, this is computationally intensive and no constraints are considered.It is hypothesised that a preliminary optimisation is carried out to determine the main parameter ranges.Berger et al. [7] present such idea in a two-stage optimisation with propeller hull interaction applied in the second stage, after evaluation of a multi-objective optimisation in given wake condition.This is an efficient intermediate step, also indicated by Foeth [6], to assist the designer in finding the optimal design for a given hull geometry.However, in [7], only a constraint on thrust equality is applied and the objective function is biased by the applied weight factor.
Constrained propeller optimisation is , e.g., presented by Vesting and Bensow [8], with a multi-objective propeller optimisation in behind condition, considering sheet cavitation, or by Han et al. [9] with a multi-objective propeller blade optimisation in given wake and sheet cavitation as a constraint.Similarly, Kamarlouei et al. [10] apply an optimisation in fixed wake, with cavitation constraints on the blade loading (Keller criterion) and the two-dimensional cavitation behaviour of one blade section by the use of the cavitation bucket diagram, from which the cavitation-free area can be identified.Kamarlouei et al. [10] also include, unlike Han et al. [9], constraints for the blade strength, though with simplification to cantilever beam theory.Bertetta et al. [5] present the multi-objective optimisation of a controllable pitch propeller (CPP), with the main aim to reduce the face cavitation and thereby the resultant radiated noise.Common denominator of the enumerated examples is the applied numerical evaluation of the propeller performance, which is based on potential methods.However, applied cavitation constraints are based on general characteristics, e.g., the maximum cavitation volume or the cavity areas, and constraints on blade strength are only included by Kamarlouei et al. [10].
This shows a clear lack of optimisations applying practical limitations and thereby an efficient constraint handling and integration in the optimisation algorithms.Puisa and Streckwall [11] present a constrained propeller optimisation with an adapted constraint-handling technique.Results are compared with the standard non-dominated sorting genetic algorithm (NSGA-II) by Deb et al. [12], which show that the problem adapted constraint handling is superior.The applied constraints on cavitation are, however, conservative and based on the circumferential mean wake, and prohibit cavitation occurrence at all.Hence, they are not entirely practical for a broad range of propellers.
We have developed and applied new constraints for sheet cavitation, based on potential flow analysis and two novel optimisation algorithms based on the NSGA-II and the standard particle swarm optimisation (PSO) algorithm [13].We proposed, in addition, a constraint handling method for infeasibility handling, which is adapted to manual propeller design.In this paper, we present their application to a broad range of test propellers and evaluate their case dependent performance, in order to develop guidelines for efficient propeller optimisation.All optimisations additionally include cavitation constraints, limitations on blade strength, class regulations and geometrical blade limitations.We also propose a recommended optimisation procedure for marine propellers.
We will first introduce the objectives for the propeller design task and outline the computational methods, before a typical manual design process is explained for the design analysis and optimisation phase.Subsequently, in Section 4 we elaborate the automated design procedure with a summary of the applied algorithms, constraints and constraints handling routine.In Section 5 we provide information about the test cases and the different optimisation strategies which is followed by the results discussion.

Propeller Design Objective and Analysis Tools
What is the optimal propeller design?To answer this question we have to interrogate the purpose of the propeller and ship it needs to be designed for.The optimal propeller performs certainly different for a container vessel than for an icebreaker, a naval, or a cruise vessel.Depending on the perspective, strength or reliability are the critical characteristics (e.g., icebreaker or naval vessels) while in other applications low noise level is the most important factor (e.g., cruise vessels).However, on a certain level all applications share the same interest for low fuel consumption, less vibrations due to propeller induced pressure pulses, safety margin for cavitation erosion and classification conformity.Therefore we restrict the variation of the optimisation set-ups to the same objective of minimisation of pressure pulses and maximisation of propeller efficiency.Both are common characteristics in propeller design which can be measured in model and full scale, and are therefore often contractually guaranteed.Additionally, we apply constraints on blade stress, cavitation, classification regulations and manufacturing limitations.These objectives and constraints are interconnected and often contradictory which requires a simultaneous consideration during the optimisation.If the propeller is optimised disregarding some characteristics the result will be sub-optimal, and an adjustment of the optimal solution will only be patchwork.
We outline a typical design procedure for the analysis and optimisation phase to find the optimal blade geometry for a given wake, inspired by Rolls-Royce design philosophy and evaluation tools.The wake is in this procedure corrected to full scale and to include the effect of a propeller.As the initial step, a lifting line method, with lifting surface correction, is employed to deliver the initial pitch and camber distribution, which is determined from the optimal circulation distribution, according to Lerbs induction factor method [14].
For the evaluation of the propeller performance, we focus on the use of moderately low fidelity numerical tools to be utilised during the early design phase.For the propeller performance, including the prediction of sheet cavitation, we utilised the vortex lattice method "MPUF-3A" [15].This methods offer sufficient accuracy to compare the results in relation with each other and a benchmark, in a reasonable short evaluation time.The propeller blade geometry is represented by a lattice of discrete vortices distributed on the mean camber surface, which have to be determined from boundary conditions.First development of this lifting surface prediction method was presented by [16].The method was extended towards an unsteady analysis tool for cavitating flow [17].Further development, e.g., by introducing distinct set of source sheets for a blade thickness distribution and leading edge correction ( [3,18,19]), yield the current version of the tool, which also includes among other features, wake alignment, mid-chord cavitation, thickness-load coupling and hub influence.
The propeller induced pressure pulses are calculated in close relation to the MPUF-3A propeller performance computation.The applied tool is a boundary elements method to compute the pressure amplitudes of the first three blade harmonic frequencies.The method solves the diffraction potential on a flat hull-dummy surface, by applying Green's formula [20].However, the objective is the maximal amplitude of the first blade frequency which occurs on one of the surface panels.
Today's propeller geometry feature commonly high skew and curved rake distributions and evolve in a highly complex three-dimensional shape.Thus numerical method for the static stress level and the dynamic strength analysis require finite element methods (FEMs).Here we apply a code which uses shell elements of quadrilateral shape, generated on the mid-surface of the blade.Evaluated are the maximum von Mises stresses, which are evaluated at a generic strength section above the blade root filling.

Manual Propeller Design
A common practice in propeller design is to develop a preliminary design concept based on best knowledge and a design database before signing the contract with the customer, and improve this design successively to find the best compromise between the competing objectives and constraints.Hence it is naturally a multi-objective optimisation task in which various constraints, intimately connected with numerical methods (e.g., for cavitation or stresses), need to be accepted and considered.
Figure 1 demonstrates what the regular manual propeller design process involves.Input variables mainly provided by the customer are allocated in the database, resulting in a number of preliminary design outputs (principally regarding shaft line) which may be saved or discarded.From the database, the propeller designer can select the desired outputs and import them into the existing software, where propeller geometry is tweaked in two separate geometry modules.This depends not only on the customers' requirements but also on the designers' knowledge and experience.Once the geometry is fixed, the designer runs the analysis modules followed by the strength analysis modules.After obtaining the design results, if they are not entirely satisfactory, the designer may go back to the first or second geometry modules and iterate by modifying geometric parameters, until a fitting result is found.Although not usual, the designer may even go back to the preliminary design and modify shaft line characteristics.The propeller geometry is during the outlined procedure developed towards, what is believed, the optimal design.Running the iterations and making decisions closely depends on the experience and skill of the designer, where time available for the design iteration is the limiting factor.Thus, the number of tested design alternatives and the lead time of analysis tools is limited.Commonly around 30 variants in 3 to 8 iterations are investigated, which can readily be reduced to less variants in case of urgency.

Automated Propeller Design Procedure
Automated propeller optimisation emerges more frequently.However, it is often very ambitious regarding computational effort in case entirely viscous CFD simulations are applied, as described by Foeth [6].The ideal case for a propeller optimisation would include the optimisation of the hull geometry.A designer should iteratively account for interaction effects between the propeller and the hull.But this broadens the design space significantly and requires cooperation between the hull and the propeller designer as well as high fidelity simulation of the hull surrounding flow.Thus, we concentrate our optimisation on the most common situation in which a designer needs to find the best propeller design for a given hull form.
We adapted common optimisation algorithms (NSGA-II and PSO) to three algorithms for propeller design and implemented and applied six optimisation strategies [13] which are in this paper introduced to a broader range of propeller types.The algorithms are implemented in the design toolbox at Rolls-Royce Hydrodynamic Research Centre and thereby have access to all analysis tools that are part of the manual design process.Thus we have access to alternate objectives, geometrical limitation, e.g., for CPP, classification regulations and various analysis tools which differ in fidelity level.
Beside the standard NSGA-II, Deb et al. [12], we apply two modified versions of the NSGA-II, which are assisted by a meta-model (Kriging model [21]), implemented in the DACE MATLAB toolbox, Lophaven et al. [22].The first algorithm, the Adaptive Surrogate-assisted NSGA-II (AS NSGA-II) evaluates variants on the meta-model and updates the meta-model frequently every fourth generation with calculated results.The second algorithm, the Surrogate-Assisted Neighbourhood Assessment NSGA-II (SANA NSGA-II) evaluates all individuals in the population with calculated results.Yet before a new generation is evaluated, the SANA NSGA-II assesses the neighbourhood of each individual with the help of the meta-model and selects the most promising design as the individual to calculate.The third algorithm applied is a Particle Swarm Optimisation (PSO) algorithm which is a combination of the NSGA-II and a standard PSO algorithm by Kennedy and Eberhart [23] algorithm.Hence, our PSO is capable of multiple objectives and non-dominated sorting and utilises an archive of best solutions which continuously updates the rank and crowding distance, with respect to the latest swarm.A detailed explanation of the algorithms is presented by Vesting et al. [13].
The optimisation results in different designs when different strategies are applied regarding the selection of parameters modification, objective and constraint selection, initialisation and handling of infeasibility.In general we apply constraints on the required thrust, the blade strength, classification regulations and geometry limitations.Additionally constraints on the cavity shape and their location are applied to (i) enhance objectives sensitivity to the cavitation on the blade, predominantly propeller induced pressure pulses; and (ii) to monitor the formation of the cavity shape and thereby reduce the risk for potentially erosive cavitation.In particular, the constraints target the cavity sheet thickness at the blade tip, the maximal cavitation volume as well as its length in chord-wise direction and the change in cavity volume between two consecutive time steps.We have additionally developed new constraints which rate the shape of the cavity and assign a harm-factor according to the similarity of the cavity closure line to a convex closure line and the location of the maximum sectional cavity centroid (e.g., harmful at the tip, root and downstream in chord-wise direction) [13].We summarize the cavitation constraints with a brief description in Table 1.

Constraint Description
CMaxVol Maximum non-dimensional cavity volume at a single blade position.

CCent harmfactor
Calculates the chord-wise centroid from the finite set of cavity thickness at each span-wise section and assign a harm-factor for the furthest downstream position to rate that position.The blade is divided into harmless and harmful sectors and the harmfactor is accumulated.

Closure harmfactor
Cavity closure line compares the normalised (length between 0 and 1) cavity trailing edge with a given, faultily example for a closure line.Such example of distinct convex shape is given by a normalised polynomial.The pairwise difference between the given curve and the sections-wise cavity length are compared and yield a root mean square (RMS) error.A small RSM imply therefore a high similarity to a convex shape and vice verse.

Cavity change
Fluctuations of the cavity volume.Maximal change in non-dimensional volume between two time steps.

CLength
Maximal non-dimensional length of each cavitating blade section.
Unlike the standard optimisation algorithms, which commonly assign a constraint violation independent from the degree of violation, we utilise for each constraint a violation measure relative to the performance of the current population.Thus, a solution is assigned with a value of −1 if it features the worst constraint performance in the current population and 1, vice versa.The other individual are assigned relative to the best and worst performance respectively.The assessment is repeated with each constraint and eventually the constraint measure is accumulated in the Constraint Violation Measure (CVM) which determines the feasibility of a solution.A solution can thereby be feasible to the algorithm even if one or more constraints are marginal violated [13].A designer judges a design differently on the level of violation and rates constraints unequal and is possibly interested in designs that perform excellent although they might be marginally infeasible, to find the next design step.The procedure we utilise, is intended to feature the designer's perspective in a manual design process.
The strategies we pursue are directional classified in a 2-stage optimisation proposal to efficiently accomplish the propeller design, Figure 2. The complexity of the design procedure makes a systematic optimisation procedure valuable.We assume in this procedure that one starts from a preliminary design concept, based on project estimations.Thereby, we imply that the operation condition (e.g., RPM, ship speed etc.) and the main propeller dimensions, like the propeller diameter, are already determined; it is about to find the wake adapted blade geometry.A sectional shape variation is out of the scope of this paper to contain the design space.In the first primary optimisation the propeller geometry is modified using parametric geometry distribution curves.These geometry distributions are related to the propeller and its application type and are controlled by general input parameters, e.g., blade area ratio, skew at the blade tip or rake at the tip.Such standard distribution curves are based on experience and design philosophy and provide the first draft of the design, similar as described by Foeth [6].Here we utilize the Rolls-Royce standard distribution curve parameters which may be different for competitors, yet the optimisation strategies will still be the same.Table 2 provides an overview of the parameter distribution curves and their corresponding control parameters.The primary optimisation makes use of these global parameters to provide an insight into the propeller response and thereby uncover tendencies and parameter influences in a subsequent assessment step done by the designer.The designer selects also a suitable optimisation algorithm, the objectives, constraints and whether constraints are handled with amplified penalties or not.The combination of constraints is crucial for the CVM and requires the designer's experience to rate the achievable performance compared to the baseline design.In this paper we go beyond the scope of a primary optimisation by applying several strategies (as outlined in Section 5) in the primary optimisation to develop guidelines for their application.The intention for a practical application is, however, to utilise only one of the strategies in the primary optimisation.The selection of which strategy to use depends consequently on the performance of the baseline, e.g., cavitation extent or constraints which are challenging.The groups in Table 4 introduce a categorisation of the strategies: Group 1 differs in application of optimisation algorithms, Group 2 differs in constraints setting and Group 3 essentially in geometry modification method.Groups 1 and 2 correspond to a primary optimisation and Group 3 to support optimisation.The settings in Table 4 are elaborated in Section 5.
Support optimisation is applied to use the gathered knowledge on trends and parameter influences for the detailed design.Parameter curves with high impact on constraints or objectives are used for local blade geometry modification to further explore and exploit obtained trends or to improve the performance of constraints, e.g., together with cavity shape constraints locally manipulate the cavity.Local modifications are applied indirectly to the B-spline parameter distribution curves by modifying two 4th order Bézier curves which supply an offset to the B-spline curve control points.The desired parameter value is iteratively achieved until the parameter curve passes the desired point.A width factor determines in the process the span wise perturbation of the parameter curve and applies thereby smooth modifications to the original curve [13].

Test Cases
The Rolls-Royce Hydrodynamic Research Centre (RRHRC), Kristinehamn, provided in total seven propeller designs for this study.Table 3 provides an overview of the propeller outline, their design purpose and main characteristics.The designs represent common applications and a wide range of geometry, e.g., blade area (EAR 0.4-0.7),diameter (D 3-7 m), skew (6°-45°) and design pitch (P/D 1.05-1.4).Apart from Opt-6 all propellers are delivered designs and already passed through a manual optimisation.Thus for the optimisation assignment we "reset" each design to an early project calculation status, considering however, final settings like hub specification and possibly changed operation conditions according to the delivered design.Hence we have for each propeller one benchmark design (the delivered design) and an initial design (the baseline) to compare the final results of the optimisations with.Elaboration on design details and absolute values will be omitted in this paper for protection of commercial interests of Rolls-Royce AB.The optimisation cases are divided into three groups in which the first two groups emphasise on the primary optimisation according to Figure 2.Here we investigate general qualities of optimisation algorithms by altering the optimisation algorithms and utilise the parametrised geometry curve for variation, which are introduced in Table 2. Table 4 provides an overview of the different settings for several optimisation cases.Each column (Opt-X-1 to Opt-X-15 where X represents each propeller respectively) represents one strategy set-up from which the majority is applied to each of the test propellers.The selection of set-ups to apply in Group 2 and 3 depends on the outcome of Group 1 optimisations.Hence, optimisations in Group 2 and Group 3 are conducted subsequent to Group 1.
For optimisations in Group 2 (Opt-X-6 to Opt-X-9), we reduce the number of input parameters by observations of Group 1 and the results of a sensitivity analysis on the parameter influence.Opt-X-6 provides the basis with fewest modification compared to Group 1.For Opt-X-7, Opt-X-8 and Opt-X-9 we, additionally, introduce cavity shape constraints which are alternated among the cases.Optimisations Opt-X-13 and Opt-X-9 are conducted subsequent and utilise amplification on constraints which are frequently violated in earlier cases.Group 2 still belongs to the primary optimisation, yet it is composed of different optimisation strategies which carefully have to be selected by the designer, in accordance to the observations on the baseline design.
The principle difference for Group 3 is the application of local modification of blade parameters and the aim to particular improve either objectives or cavity shape.Hence, in Group 3 we apply either Opt-X-10 for objective improvement or Opt-X-11 for cavitation improvement.Opt-X-14 and Opt-X-15 are conducted supplementary according to either Opt-X-10 or Opt-X-11 to judge the capability of meta-model assisted optimisations for the local blade modifications.
For all optimisation cases we specify the same settings for the optimisation algorithms with settings for population size (N I = 24), number of generations (N G = 40), mutation probability (p mut = 0.75), crossover probability (p cross = 0.15) and, for the PSO, inertia weight (w = 0.4).Solely the population size and the number of generations is revised in case of Group 2 and 3 where the number of parameters is reduced.Parameter limits are generally set to a maximal variation of ±10% to the initial value of the baseline design.
However, special consideration is given to parameters like EAR for cases Opt 3 and Opt 5.The initial designs provide already a high blade area ratio.Hence, an increased value up to +10% yields, as expected, blade collision and other hub related problems.Thus, the upper limit is revised to meet the geometric requirement.The same holds for the rake parameter since this value is given in mm which results in too few variation within the 10% range.Consequently, we apply a wider range of 1000% compared to the baseline.The constraint limits are similar related to the baseline design with a limit of 5% improvement on cavitation and blade stress constraints.Only Opt-7 requires special consideration of the blade stress constraint since the baseline yields unacceptable high stresses and consequently the limit is set according to the benchmark design.The required thrust constraint (K T ) is as well set according to the benchmark design, to ensure a fair comparison with the benchmark design.

Results and Discussion
The optimisation survey involves the results of up to 15 optimisation set-ups for each of the seven propellers, which adds up to more than 80 optimisations and thereby a comprehensive number of potential optimal solutions that require assessment.Consequently, there is a need for qualitative and quantitative criteria to determine the convergence and quality of the optimisation itself, before appraising the optimal solution.In general we classify the assessment into three assessment ranks: • Potential of the optimal solutions, to rate if the optimisation was actually able to improve the design.• Convergence of the selected set-up, to rate the confidence of the optimisation to actually end up in a specific location, the global optimum, and thereby provide information about the computation effort of the optimisation approach.• Constraints feasibility, to judge if the approached optimal solutions comply with the constraints.
We divide the result presentation further according to the categories in Table 3 for either performance or constraint optimisation.Each propeller will be assessed separately, starting with a summary table of observations according to the assessment ranks of potential improvement and convergence followed by a discussion of constraints feasibility and comparison with the benchmark design.To examine the central tendency of input and output quantities, we evaluate predominantly the medians within the population of one generation.The median is a more robust measure of the location parameter in case of skewed or non-normal distributions.It reflects the central development of a parameter or objective over the progress of the optimisation remarkably well since it reduces the importance to outliers.

Performance Optimisation
The propeller Opt-1, Opt-2 and Opt-4 require thoughtfulness regarding the objectives due to e.g., contract penalties.These propeller are optimised to enhance the objective performance (η and p).In Groups 2 and 3 we select parameters with expected importance for the objectives (e.g., pitch and chord) and according to the parameter sensitivity.

Opt-1
Common geometry developments, among all optimisations on Opt-1, are exclusively obtained for the rake, which increased, and EAR, which decrease.Opt-1-3 decreases the blade area with 10% down to the parameter limitation, which is misled by an optimistic estimation of pressure pulse, cavitation and blade strength constraints (Figure A1a and Figure 3a).Both blade stresses and cavitation are up to 10% better estimated by the meta-model than simulated.Opt-1-7 and Opt-1-8 decreases its EAR parameter to a similar low level and thereby improved the blade weight objective in Group 2. Both fail, however, as a consequence in adherence of the blade strength constraint (Figure 3b).
Constraints development is given by the CVM in Figure 4 and the median of the constraint values in Figure 3.These figures reveal that the results of Opt-1-1 and Opt-1-6, with fair improvements according to objectives medians (Figures A1 and A2) are virtually infeasible after generation 20, although all constraints are feasible.The only constraint that gently declines from 5% to 4% after 10 generations is K T for Opt-1-1 (Figure 3a).A detailed analysis of the constraints of Opt-1-6 show that this is due to a steep increase in blade stress within the first 10 generation up to 10% higher than the baseline design (Figure 3b).
The algorithm tries to counteract the decline which leads to a sudden increase of cavitation volume by 5% at generation 13 and high blade stresses.However, this should still result in a positive CVM but is prevented by an increasing number of failed simulations, with the consequence of unreliably small values of K T and annulled cavitation.Since the CVM is a relative measure, depending on each the maximum and minimum constraint values within a generation, the consideration of such outliers yield small or negative violation measures for the other solutions in that generation.Our CVM has thereby the additional effect of penalising failed simulations.
To compare the influence of different cavitation constraints, we have to select particular designs from the Pareto optimal solutions for each case in Group 2. This is done manually to find the best trade-off among the objectives and to ensure that the K T constraint is to the greatest possible extent fulfilled, i.e., we except a higher K T of up to 4% compared to the required K T for equilibrium.Both selected solutions of Opt-1-6 and Opt-1-7 violate, however, the blade strength constraints with 9.2% and 3.4%, as indicated by the median assessment (Figure 3).The advantage of Opt-1-7 is otherwise a significant cavitation volume reduction of about 20%.From Figure 5 it is visible that this case is the only result in Group 2, where the cavitation thickness at the blade tip is reduced.Opt-1-7 utilises cavitation constraints particularly suited for the sheet thickness at the three outermost blade radii (r/R ≥ 0.95).Both Opt-1-6 and Opt-1-7 introduce geometrical changes that yield mid-chord cavitation (Figure 5a,b).This cavitation is on the other hand prevented in case of Opt-1-8 using the centroid harm-factor constraint.For the final comparison with the benchmark design, we select one solution from all Pareto fronts that satisfies the constraints and provide the best trade-off among the objectives to compare with.For Opt-1, the solution is obtained by the PSO algorithm in Group 1, with 0.83% higher efficiency, 14.14% lower pressure pulses and 32.56% smaller cavitation volume than the benchmark design.The maximal blade stresses are maintained on the same acceptable level as the benchmark propeller (+0.18%).In Figure 6 the simulation results are compared of both designs.The optimisation utilising only parametrised geometrical curves is successfully conducted and yields a design solution comparable with the manual work.Only the blade tip is different compared to the benchmark propeller which features a particular tip design; the sectional shape is modified to decrease the tip leading at the design point and suppress tip vortex inception.This requires manual interaction in the design process and is based on the designer's experience.Due to the absence of an analysis tool for the tip vortex inception, the optimisation yields a higher loading at the tip. Figure 6a indicates, however, a smaller sheet cavitation at the outermost radius.The main conclusions regarding Opt-1 are summarised in Table 5.Table 5. Potential improvement and convergence observations of Opt-1.
• Opt-1-2 and Opt-1-4 vary stronger in pressure pulse objective, which coincides with a variation of the blade tip loading (Figure A1).• Opt-1-4 varies parameters for root unloading and EAR (Figure A1a).Both EAR and unloading factors correlate: a reduction of EAR coincides with a higher loaded blade and thereby balances the efficiency.

II
• NSGA-II (Opt-1-6) superior objective capabilities (Figure A1b) but fails regarding constraint feasibility.• Opt-1-7 is similar capable in finding good objectives, however, not better with regard to the constraint feasibility.• Opt-1-8 (CCent constraint) is the only set up that manages to improve all three objectives and constraint feasibility, Figures A1 and 4. • No enhanced convergence, although less input parameters reduce the design space.
• Clearest approximation of final objective values by Opt-1-8 with virtually no improvement in η and no reduction blade weight (Figure A1).

Opt-2
The geometry of Opt-2 (both in Groups 1 and 2) evolves similarly for the majority of cases towards a higher EAR, Skew, Rake and root unloading, Figure 7.Only two parameters vary during the continuation of Opt-2-2: thickness and tip unloading, Figure 7a, which strengthen the assumption of low importance of these parameters.Hence, in Group 2 we reduce the number of decision variables down to three by omitting thickness reduction, blade unloading and skew line starting location on the hub (SkLFhub).The sensitivity analysis provide the fourth lowest importance for the root unloading and is therefore omitted in cases Opt-2-6 to Opt-2-9.This consequent reduction of parameters chokes off the optimisation, and the only true improvement in objectives and constraints is obtained by Opt-2-13.Which is conducted subsequent to Opt-2-6 to Opt-2-9 and uses again two more parameters (TMAX and root unloading).
The strongest effect can be assigned to the unloading root factor, by analysing the development of the additional parameter towards a root unloading of up to 50%, Figure 7b.Neither the tip unloading nor the blade thickness evolves to a significant geometry change.The positive efficiency median (around 0.3% in Figure A3b), despite larger blade area, can be seen as the result of less rotational losses due to the root unloading, which has not been captured by the sensitivity analysis.
The consequent violation of the K T constraint of Opt-2-1 is evident from Figure 8a.Yet it is not seen by the algorithm since the CVM is positive (>1%).The CVM balances infeasible K T with significantly feasible blade stresses and caviation (−10% in cavitation volume and −30% in blade stress).In case of the SANA NSGA-II (Opt-2-4) a sudden change of EAR at generation 16 and the consistent increase in cavitation abrogates the balance and CVM drops consequently, Figure 8a.The development of CVM reveals also that the meta-model of the AS NSGA-II (Opt-2-3) underestimates the constraints which spikes at every fours generation; particularly the stress prediction of the meta-model, deviates from the calculated values by 5% (Figure 8a).In Group 2 the K T constraint violation is improved by two of the cases, Opt-2-8 and Opt-2-13 (Figure 8b).Both Opt-2-8 and Opt-2-13 reduce the number of variants violating K T , or K T in combination with other constraints, from about 60% of variants (Opt-2-6 = 60.3%,Opt-2-7 = 69.2%) to 17.2% (Figure 9b) and 34.8% (Figure 9a) respectively.Opt-2-8 uses cavity shape constraints, which are hardly met or affected by applied geometrical changes, hence K T is the only characteristic that can be modified.Opt-2-13 amplifies the K T constraints and thereby increases its influence on the CVM.The additional parameters, which are utilised in Opt-2-13, are contradictory to an increase in thrust, since the optimisation decreased the blade loading.The improved K T is consequently a result of the amplified constraint handling.We select a Pareto optimal solution from Opt-2-13 for the the detailed analysis of cavitation on propeller Opt-2 (Figure 10).The sheet cavity volume of the baseline is substantial at the θ = 0°position and increases at the blade tip in thickness when the blade passes the θ = 30°position.Hence, two constraints are expected to improve the cavitation, the constraints on (i) maximal cavitation volume; and (ii) on cavitation sheet thickness at the tip.Opt-2-6 achieves with the CMaxVol constraint the best median improvement of −5% (Figure 8b).The Pareto optimal solution is, however, less improved (−0.12%, Table 6).The cavity thickness constraint fails in Opt-2-7 primarily due to a lack of sufficient parameters to modify the blade geometry locally at the tip and therefore only the cavity change constraint affects the design during the optimisation by reducing p (Table 6), which is the only case that achieves a reduction in pressure pulse together with increased cavity volume.Opt-2-8 enables cavity centroid constraints which, again, occurs most harmful towards the tip region where the parameters are only marginally effective.Table 6.Selected optimal solutions relative to baseline design.

Optimisation Variant η [%] p [%] CavVol [%] FEM [%] K T [%]
Opt In Group 3 we focus on the performance improvement which is achieved according to Figure A3c by local modification of pitch at r/R = 0.6 and r/R = 1.0, chord length at r/R = 0.5, camber at r/R = 0.8 and skew at r/R = 1.0.The analysis of constraints reveal that the CVM is in fact positive for all three cases in Group 3, yet the efficiency gain is only achieved by violating the required K T (Figure 11b).The CVM balances with considerable improvements in blade stress (Figure 11d) and cavity volume (Figure 11c).The performance of selected solutions in Group 3 which provide improvements in objectives and constraints and primarily fulfil the required K T are given in Table 7, relative to the performance of the benchmark design.All three algorithms in Group 3 yield the same level of improvement in objectives, blade stress and cavitation reduction and can be considered to achieve the same optimisation quality.This highlights an advantage of Opt-2-14, which utilises the adaptive meta-model, and reduces the lead time by a factor of five.The variations in constraints due to the difference between estimations and calculations are high at the beginning (Figure A3c), but shrink with improving meta-model, and yielding the same geometry as employing only calculated results.The main conclusions regarding Opt-2 are summarised in Table 8.
Table 7. Optimal solutions of Group 3, relative to performance of manual design.

II
• Less improvement possible compared with optimisation of Group 1.

Opt-4
The underestimation of pressure pulses in both cases that use the AS NSGA-II (Opt-4-3 and Opt-4-14) yields a thickness reduction which is opposed to the other algorithms.From Figure 12b it becomes particularly obvious that, at generations 4 and 8, and from Figure 12a at generation 4, a thickness reduction yields a positive effect on the objective (Figure A5a,c) and in both cases initiate the lasting trend for the thickness reduction.The input parameters are reduced in Group 2 by all three pitch parameters (unload factor root, unload factor tip and unload width).Their effect on constraints, objectives and convergence is marginal which is consistent with the result of a sensitivity analysis, employing the extendedFAST method by Saltelli and Bolado [24] (Figure 13).The importance measure is given as an aggregated sensitivity measure of the sensitivity indices of each objective and constraint.The unloading importance is assigned low which results from a varying unloading importance among the objective and constraints; both unloading parameters affect the cavitation rather than the efficiency.However, the span wise distribution of unloading (Unload-Width) is more important for the effect on efficiency than for the cavitation.This explains why the method assigned a higher importance to the unload width when the unloading itself is assigned less important.There are two more parameters, the skew line position on the hub and the rake, which are assigned with low importance, however, we included both in Group 2, since the skew at the tip is assigned with high importance and an effect of rake on the efficiency is noted in other cases.The results of Group 2 show, by strict convergence of rake, that this parameter is of importance for the obtained performance.However, for the performance increase in Group 3, we initially applied modifications of pitch and camber at r/R = 0.7 and r/R = 1.0 and of chord length at r/R = 0.6, which resulted in incapable performance (Figure 14a).The adoption of an additional chord length and the thickness parameters improved the performance (Figure 14b), yet only the inclusion of rake and skew at the tip modifications yield extended improvements (Figure 14c).This is consistent with the result of the sensitivity analysis.The margin for constraints violation is sufficiently large for Opt-4.The K T constraint is feasible according to the definition as an inequality constraint and the optimisations adapt from a too high K T towards the desired value of the baseline (Figure 15a).The stress constraint is similarly met and improved during the optimisation by −10% to −20% Figure 15a.
Comparing the baseline design with the manually designed benchmark propeller reveals 2% higher efficiency for the benchmark design, which is a significant difference between the baseline and benchmark design.The manual work on the design increased the performance which is similarly achieved only by one of the Pareto optimal solutions (Opt-4-10), featuring an increase in cavitation volume.However, even if the optimisations achieve pressure pulse reduction of up to 18%, compared to the baseline, they are still 3% higher than the benchmark design.The solution discovered by the optimisation is not quite as optimal as the manually designed blade.This is possibly the consequence of too narrow box limits on the selected parameters, since they converged in all cases rapidly to the limits and a baseline design which is far too off to provide the required performance.However, the design of the support optimisation is close to the benchmark design with 0.23% efficiency increase and 21.6% reduction in blade stress, compared to the benchmark.The main conclusions regarding Opt-4 are summarised in Table 9.Table 9. Potential improvement and convergence observations of Opt-4.

I, II, III
• The potential objective improvement is yet exceptionally good, compared to the baseline design.
All optimisations of the three groups in Figure A5 yield improvements in both η and p objectives.• In Group 2, the consideration of weight as an objectives, is successfully fulfilled, with a median improvement of 1.5% in efficiency and more than 12% reduction in pressure pulse (relative to baseline design).• Opt-4-2 (PSO) focuses on η objective.
• Variations of Opt-4-3 and Opt-4-14 results from underestimation of objectives by the meta-model (Figure A5a,c).• Large discrepancy between baseline and benchmark; gain is small compared to benchmark design.• Opt-4-4 balance EAR and blade loading.

Opt-3
Assessment of optimisation results have extensively been discussed in [13] and therefore only a brief summary is provided here.Several propeller alternatives originate from the baseline design and were tested in model tests as part of an internal design study.The results confirm that this design is all about to balance the blade loading such that tip vortex strength decreases, but erosive mid-chord cavitation is still contained.
The cavitation constraints and handling of constraints by amplification in Group 2 show that optimisations with enabled cavity shape constraints (Opt-3-7 to Opt-3-9), result in an increased number of solutions with lower pressure pulses.Additionally, it is possible to control the cavity shape and the objectives with the proposed constraints such that the optimisation yields different geometries.For instance, an amplification of cavitation length yields a reduction of the cavitation length median of 2% which on the other hand increases pressure pulses due to a constant cavity volume.Opt-3-8, with constraint on cavitation centroid, yields feasible solutions with a centroid relocated upstream by increasing the cavitation at the leading edge (Figure 16b).Optimisation case Opt-3-9 also utilises an amplified K T constraint and thereby enables a drop in K T at generation 6 (Figure 17b), down to the baseline value which eventually allows for raising efficiency between generation 1 and 10 (Figure 17a).Local modification are applied primarily to improve the cavitation.The results show that the selection of parameters, based on the location of cavities, and the results of a sensitivity analysis, yield the desired cavitation improvement and maintenance of the objectives.The optimal solutions are remarkably similar among Group 3, except for a difference in thickness.This highlights an advantage of the AS NSGA-II (opt-3-14), which reduced the lead time by a factor of 5, again.The SANA NSGA-II requires approximately the same computational time as the NSGA-II, yet achieves the best cavitation prediction (Figure 16c).
The Pareto optimal solutions of Group 3 are comparable with the manually designed benchmark propeller.The best trade-off design is achieved by Opt-3-15 which provides an increase in efficiency, reduction of pressure pulses but still a higher maximal cavitation volume, Table 10.The cavitation of this design (Figure 16c) is in fact distributed over a broader blade section, yet it is thinner compared to the manual design (Figure 18a), hence, contains less focusing energy and is therefore considered to offer a higher margin for cavitation erosion.Beneficial is also the stress distribution over the blade for Opt-3-15 design since it avoids higher peaks at the root (Figure 18).The main conclusions regarding Opt-3 are summarised in Table 11.
Table 10.Pareto optimal solutions of Group 3, relative to performance of manual design.

I
• Potential improvement between the algorithms diverge, with an absolute difference of 1% in η and 8% in p. • Opt-3-1 (NSGA-II) and Opt-3-3 yield no increase in η and only 2% reduction of p.

II
• Performance potential improved compared to NSGA-II in Group 1 (Opt-3-1): 0.1%-0.2%higher η.• The reduction of decision variables does not yield a positive effect on the number of variants needed to convergence.
• Median p increases up to 0.2%.
• All algorithms perform similarly.

Opt-5
Figure 19a provides that feasible CVM medians are achieved in Group 1 by PSO optimisation (Opt-5-2) and to a large extent for the NSGA-II (Opt-5-1) optimisation.The PSO convergences rapidly within 6 generations and provides a CVM which is annulled and thereby feasible.In case of SANA NSGA-II (Opt-5-4), it is obvious that the large increase in efficiency is traded for feasibility.A violation of the thrust coefficient K T is, however, present in almost all optimisations and where CVM balance mainly with feasible blade stress constraint but partly by a feasible cavitation volume (e.g., Opt-5-2).
No constraint amplification and no cavity shape constraints (Opt-5-6) results in an entirely infeasible set of solutions, mainly violating the required thrust constraint in combination with other constraints on the cavity (94.1%).The CVM is yet still positive (Figure 19b) due to improvements in blade stresses (Figure 20b).Cavitation, however, is violated which can be seen from the constraints development in Figure 20a towards a median increase of about 140%.In Figure 21a a solution is presented that improves the objectives and violates the K T with only −0.36%.Yet a cavity volume increase is evident at the root and over the blade surface as mid-chord cavitation.The cavitation behaviour improves significantly even without amplification solely by enabling the cavitation centroid constraint (CCent).Opt-5-8 utilises this constraint and yields in general a CCent median of the baseline level, but suffers significant violation of K T with 88.1% of the variants.The improved mid-chord cavitation is shown in Figure 21c.This variant is selected from the Pareto front, satisfies the K T constraints and yields an improvement in both objectives (η =1.5%, p = −2.7%).
The only optimisation case that distinctly reduces the violation of K T is Opt-5-9, which enables the amplification of the K T constraint.The number of variants that violates K T exclusively (31.8%) and K T in combination with other constraints (46.4%) is reduced to in total 78.2%.
Optimisations in Group 3 aim to improve the cavitation behaviour by reducing the root cavitation and removing the entire mid-chord sheet cavitation.The baseline design features no mid-chord cavitation, but the increase of efficiency eases the inception speed and mid-chord cavitation is present for all results in Group 1 and 2. To still improve or maintain the objectives, we introduce parameters according to Table 12.Cavitation constraints are set to a threshold of a 5% smaller maximal cavitation volume and cavitation length and a positive centroid harm-factor, which corresponds to an improvement of −133%.In Opt-5-11, the demand for a extreme centroid harm-factor reduction (CCent constraint) yields significant reduction of pressure pulses (median up to -11%, Figure A6c) and consequently no efficiency improvement.The cavitation volume and the centroid are, however, reduced by 80% and 150% respectively, while the K T is sacrificed and Opt-5-11 finalises the optimisation with a median K T of −15% (Figure 22b).
Consequently, for Opt-5-14 and Opt-5-15 the high demand for the centroid harm-factor is reduced to a moderate enhancement of 33%.This results in a better performance regarding efficiency and K T , Figure A6c and Figure 22b.The geometry development is still the same, except for an opposite rake development for both rake modifications (Figure 22a RAKE1).
In Group 3, all Pareto solutions are violating the K T constraint.We select therefore designs that improve the objectives and predominantly satisfy the required K T .Consequently, the gain in objectives is virtually not present (Table 13).The performance in Table 13 is given relative to the benchmark design.Although the cavitation volume is still larger than the benchmark, there is an improvement compared to the designs in Group 2. We compare the cavitation prediction of the benchmark design (Figure 23a) and the selected final designs in Figure 23 which provides quantitatively that the root cavitation is reduced in case of Opt-5-14 to the same level as the manual design and that the mid-chord cavitation vanishes.In terms of performance, the optimal design (Opt-5-15) is comparable with the benchmark, yet with 1.65% higher pressure pulses.The blade stress is reduced by 27.97% compared to the benchmark.The main conclusions regarding Opt-5 are summarised in Table 14.Table 13.Performance of optimal solutions, relative to manual design performance.
• Opt-5-1 yields a better trade-off between both objectives.• Fast convergence in the parameters EAR, Skew tip, Rake tip and unload factor root; converge within the first 15 generations.

II
• The reduction of decision variables does not yield a positive effect on the number of variants needed to convergence.

III
• Less improvement potential due to strict constraints on the mid-chord cavitation together with local modifications tailored for the cavitation improvement.• Underestimation of p by meta-model of Opt-5-14.

Opt-6
The convex closure line of the sheet cavitation can be predicted by MPUF-3A (Figure 24a) sufficiently well.Model tests and simulations with a boundary elements potential flow code confirm an even more pronounced convex closure line and the gap in radial directions of a cavitation section around r/R = 0.4.This case is therefore a favoured example to test the introduced cavity closure line constraints (Closure) and we are predominantly interested in the effect of that cavity shape constraint.Optimisation is hence only conducted within Group 3, aiming for a maintained efficiency and a sheet cavity that reduces the convexity of the cavity closure line and thereby reduce the risk for re-entrant jets and break-off.The harm-factor is intended to be improved by −60% by applying in total six parameter changes according to Table 15.The potential increase in efficiency is evident from Figure 24c, which shows a continuous rise until it converges to 0.5% increase at generations 15 to 20.The high threshold for the closure line harm-factor is, however, not achieved.The parameter development is similarly blunt as the convergence of efficiency and most constraints indicate.Within about 20 generations all parameters gravitate towards a final value, which is in general a reduction in both chord length parameters of about 4%, a mild reduction in camber, a significant reduction of thickness by 8% and in the pitch at the tip.All together, the changes indicate a dominance of the objective as driving force of the optimisation.Merely the thickness and camber reduction can be accounted to reduce the low pressure on the suction side and thereby the cavitation.A modification of camber and thickness distribution is preferable to locally change such cavitation.
Figure 24b indicates that the closure line constraint is violated in 97.3% of the variants.The specified centroid target value is missed in 83.7% of the solutions.A detailed analysis of the constraints in Figure 25 reveals that mainly the cavity centroid (CCent) is improved.The maximal cavity volume (not active) still increases, compared with the benchmark design.The same holds for the cavity closure line constraint which increases up to 10%.
However, the optimisation yields variants with improved cavitation behaviour, e.g., Opt-6-11-104 (Figure 26a) reduces the cavitation and more importantly changes the cavitation closure line such that the convexity vanishes with constant efficiency (η = +0.007%)relative to manual benchmark design.The maximal volume reduces with 30%.The variants Opt-6-11-334 and Opt-6-11-339 (Figure 26b,c) provide a different interpretation of the optimisation set-up, such that the cavitation satisfies the constraints in vicinity of the blade with the maximal volume and increases the efficiency, yet deteriorates the cavitation outside the range of considered blade positions.From these two variants, Opt-6-11-334 distinguishes itself by the best performance with regard to cavity centroid harm-factor, with an improvement of about 20%.The variant Opt-6-11-532, with the best closure line constraint, provides long cavitating sections which reach the trailing edge of the blade (Figure 26d).This behaviour explains the low harm-factor for the closure line constraint since it deviates significantly from the given convex closure line and discloses thereby the weakness of that constraint.

Opt-7
In all optimisations, variants are created along a suboptimal Pareto front which is evident from the convergence of the optimisations in Figure 27 where the colour scale indicates the progress through the generation number.The optimisations find Pareto optimal solutions similarly located within the first generations, yet all converge to different location when the algorithms progresses.Opt-7-1, for instance converges to a broad band parallel to the Pareto front while Opt-7-2 gravitates to a region of (−1%η, −20%p) and Opt-7-6 to (0.6%η, 5%p) which is on the opposite end of the Pareto front.Opt-7-4 splits and converges to both ends of the Pareto front.This behaviour results, however, in a difference between Pareto optimal solutions and solutions of the converged optimisation, which is particularly obvious from the blade stresses.
Pareto optimal solutions yield designs with short chord length and suffer a severe violation of blade stresses (Figure 28).The designs are similar and the blade stresses of the Pareto optimal solutions are alike the baseline design, which is the consequence of lack of exploration since all designs are created at the beginning of the optimisation.Selecting, on the other hand solutions from the final generations reveals a notably improvement in blade stresses by all optimisation cases, Figure 29, which are more similar to the benchmark design (Figure 29a).Yet Opt-7-2 provides a design which is similar to the benchmark design in terms of blade loads (increase of only 6%) while the pressure pulse objective is reduced by 25% (Table 16).The gain through optimisation is in case of Opt-7 small since the parametrised standard geometry curves offer not sufficient variation of the blade geometry.The chord length of the optimised solutions is, e.g., much shorter at the root than the manually designed blade (benchmark).The main conclusions regarding Opt-7 are summarised in Table 17.

I
• Potential gain in efficiency is the worst among the test cases (Figure A8).
• Opt-7-1(NSGA-II) achieves a median value of a level similar to the baseline design.
• The pressure pulse reductions is yet substantial with a possible reduction of about 20% (Opt-7-1 in Figure A8a).• Already the baseline design yields a severe violation of approved maximal stress (Figure 28a).• Opt-7-2 most distinct and rapid convergence.
• Opt-7-2 algorithm is the only that achieves the required reduction of stress (FEM in Figure 28c).• Opt-7-4 algorithm balances the blade loading parameters and the EAR parameter.• Opt-7-4 yields the second best stress reduction of around 20%.
• Opt-7-3 and Opt-7-1 achieve in fact a reduction of blade stress, yet only a median of about −10%.

II
• Reduction of input parameters, by omitting rake and skew line starting location at the hub, results in better convergence.• Opt-7-6 convergences strictly towards an improvement of η and weight, disregarding the p objective (Figure 30b) and the blade stress constraint (Figure 30b).

Summary of Results
Frequently, the NSGA-II provided the best compromise between two objectives and converged in both objectives towards an approximation with fewest variation.Yet compliance with the constraints was insufficient.In many cases, either the CVM was violated or the NSGA-II was outperformed by the other algorithms regarding compliance of individual constraints, e.g., the required thrust or maximal blade stress.The blade stress reduction is achieved with larger improvement margin (≈−20%), which can thereby balance the CVM.However, the parameters that are likely to yield a stress reduction, i.e., chord length and blade thickness, are contradictory to the objectives.Thus, the optimisation should find a natural balance.
The proposed PSO and SANA NSGA-II algorithms perform similar in most of the cases regarding objectives and particular on adherence of constraints, which were, in contrast to the NSGA-II, more often fulfilled.The PSO is characterised, beyond the constraints compliance, also through a rapid convergence, but focuses on one of the objectives.The SANA NSGA-II (Opt-X-4) yields frequently to high fluctuations of EAR and unloading which balance each other: designs with a small EAR provide the advantage of less wetted surface (given that the blade radius is kept constant) and thereby less viscous drag which increases the efficiency.However, the smaller blade needs to compensate the generated thrust with a higher loading which again results in higher pressure pulse.The intense variation is related to the neighbourhood size in the SANA algorithm, which defines the variation by the standard deviation of a normal distribution.The AS NSGA-II performed, with the setting in Group 1, not satisfactory.The estimations by the meta-model deviated widely from the calculated results which was obvious from the objectives and constraints.However, the meta-model adapted towards the end of the optimisation.In case of optimisations in Group 3, with localised modifications on the blade geometry, the AS NSGA-II performed significantly better and comparable with the NSGA-II, with the only exception of cruise ship propeller Opt-4.
In Group 2 the number of input parameters was reduced based on observations from Group 1 and supplementary with the results of a sensitivity analysis.Additionally, cavity shape constraints and amplification and constraints were introduced.The reduction of parameters, and thereby the design space, yielded, however, only improved convergence in one of the cases, independent whether additional shape constraints were applied.The cavity shape constraint on the cavity centroid harm-factor and the sheet thickness constraint performed well and managed to control the optimisation in the expected manner.Particularly during the support optimisation (Group 3) of the constraints category (Opt-3, Opt-5 and Opt-6), designs could be improved regarding cavitation utilising the shape constraints.
In comparison with the manually designed propellers, according to the established design procedure, the automated optimisation provide only in case of Opt-4 insufficient performance.In that case, the parameters converged exceptionally rapidly towards the parameter limits.Opt-4 baseline design deviated significantly from the benchmark design, which indicates that the parameter limits were too narrow to achieve the same performance.In the majority of the cases the optimisation yielded designs comparable with the manual work with advantages in constraints, e.g., lower blade loads or less cavitation or cavity shapes associated with higher margin to cavitation erosion, like removed mid-chord cavitation.The possible improvement in η is only limited since the preliminary design philosophy yields designs that are already optimised to find the circulation distribution in a lifting-line and lifting-surface design method.

Conclusions
Manual propeller design is supreme to our current method.Designing a propeller incorporates many considerations regarding constraints and limitations, truly a multi-objective and multi-disciplinary optimisation task.An experienced designer considers regulations, limitations due to the type of propulsion systems (e.g., podded propulsion) and expected performance output with experience and makes reasonable decisions.The challenging factor in the manual design process is, however, the time which limits the possibilities to find the optimal design.This paper presents optimisation algorithms and strategies applied to several propellers with different application purpose.In three cases the optimised design improved both, objectives and constraints, in three other cases the constraints are improved with maintained objectives, compared to the manual design.Only one case provided insufficient amendment due to deficient scope for development of the parameters.In all cases, optimal designs render, however, propeller designs that are comparable to the manually designed propeller blade geometries.
The optimisation handles multi-disciplinary evaluations and multiple constraints without user interaction.We recommend however a two-stage optimisation procedure with an assessment step in between to apply gathered knowledge and further exploit trends.The support optimisation, with local modification of the blade geometry, features particularly good performance to modify cavitation on the blade, together with the proposed cavity shape constraints.The AS NSGA-II is in this second optimisation phase comparable in performance with the NSGA-II and the SANA NSGA-II.The meta-model yields more accurate estimations due to reduced variation of the inputs.The AS NSGA-II is therefore beneficial in this stage as performance is better.The PSO or the SANA NSGA-II are advantageous in the primary optimisation due to superior constraint compliance.The aggregated constraint violation measure increases the flexibility and thereby the convergence of the optimisation.With amplification of a constraint, it is able to feature the designer's judgement of constraints but the CVM is greatly sensitive to the selection of constraints and their limits.
The purpose of this paper was to apply an automated optimisation procedure to a wide range of test propellers and evaluate the optimal results to determine propeller depending optimisation settings.The optimised propellers are, with respect to the applied settings, reasonably good; the designer may, however, consider design limitations that are beyond the applied evaluation methods.Yet, the generated designs guide the designer to find better solutions faster.To determine reliable guidelines for optimisation settings, further optimisations need to be assessed and categorised according to propeller design demands, e.g., a focus on cavitation, propeller efficiency or pressure pulse.In that context, the optimisation routine can be improved regarding stricter consideration of the K T constraint and in complexity, e.g., by consideration of several operation conditions or interchangeable settings of the objectives and constraints can be beneficial to achieve better and possibly novel designs faster.

Figure 1 .
Figure 1.Manual design process and iteration.
load distribution at the root by changing how unloaded it should be.(0.0 = No unloading, 1.0 = Full unloading) Pi/D Unload factor tip Modification of the load distribution at the tip by changing how unloaded it should be.(0.0 = No unloading,1.0= Full unloading) Pi/D Unload width factor Influence zone for adjustment of load distribution.

Propeller Characteristics Category 1 • 2 • 3 • 5 • 6 •
twin-screw cruise vessel • strict requirements on high η and low p with contract penalty • risk for tip vortex inception • CPP Performance FPP propeller for container vessel • integrate rudder and propeller (hub cap and propulsion bulb) to optimise hydrodynamic efficiency • high requirements on η Performance naval vessel • highly loaded • limitation on draft and propeller-hull clearance • high risk for erosive mid-chord cavitation Constraints 4 • cruise vessel • Mermaid pod, i.e., electric rotatable thruster • high requirements on η • low stresses at normal straight inflow Performance Ro-Pax vessel • high ice class requirement (LRS 1A Super) • high demand on low p • risk for mid-chord and root cavitation Constraints naval vessel • highly loaded • design alternative of Opt-3 • risk for travelling cavities (convex closure line of the sheet cavitation) Constraints 7 • cruise vessel • upgrade to an existing vessel to reduce fuel consumption • high requirements on η • risk for high blade stresses • no cavitation simulation Constraints

8 (b) Group 2 Figure 3 .
Figure 3. Development of the constraints, given as the generation median of K T , maximal cavity volume and blade strength, relative to baseline design; Opt-1.(a) Group 1; (b) Group 2.

8 (b) Group 2 Figure 4 .
Figure 4. Development of the feasibility, given as the generation wise validity ratio and the generation median of Constraint Violation Measure (CVM).(a) Group 1; (b) Group 2.

Figure 7 .
Figure 7. Parameter development, given as median per generation, relative to baseline design.(a) Selected geometry parameters of optimisations in Group 1; (b) Geometry parameters of optimisations in Group 2.

Figure 13 .
Figure 13.Importance of input parameter according to sensitivity analysis.

Figure 14 .
Figure 14.Pareto front in objective space for initial set-ups of opt-11-10 (a,b) and final set-up (c).

Figure 15 .
Figure 15.Development of K T (a) and FEM constraints (b) in Group 3.

9 (b) Group 2 Figure 19 .
Figure 19.Development of the feasibility, given as the generation wise validity ratio and the generation median of CVM.(a) Group 1; (b) Group 2.

Figure 20 .
Figure 20.Development of the constraints, given as the generation median relative to baseline design.(a) Cavitation volume constraint of Opt-4 (Group 2); (b) FEM constraint of Opt-4 (Group 2).

3 Figure A6 .
Figure A6.Analysis of the constraint violation in Group 2, with respect to the total number of solutions in each optimisation case.(a) Opt-5 Group 1; (b) Opt-5 Group 2; (c) Opt-5 Group 3.

2 Figure A8 .
Figure A8.Development of objectives, given as median of relative improvements per generation.(a) Group 1; (b) Group 2.

Table 2 .
Standard distribution curves and their corresponding controlling parameters.

Table 3 .
Benchmark propeller and their characteristics.

Table 4 .
Optimisation test matrix with alternating algorithms, geometry variation and constraint settings.(C: Applied in case of cavitation optimisation; E: Applied in case of efficiency performance optimisation.)

Table 8 .
Potential improvement and convergence observations of Opt-2.

Table 14 .
Potential improvement and convergence observations of Opt-5.

Table 16 .
Converged optimal solutions, relative to benchmark design.

Table 17 .
Potential improvement and convergence observations of Opt-7.