Numerical Analysis and Validation of an Optimized B-Series Marine Propeller Based on NSGA-II Constrained by Cavitation

<jats:p>Constantly growing environmental concerns focused on reducing pollution, in addition to rising fuel costs in recent years, have led the maritime industry to develop and implement fuel-saving solutions. Among them is the optimization of marine propeller efficiency, as marine propellers are a crucial part of ship’s propulsion system. During the initial design stage, selecting the optimal propeller is considered a multi-objective optimization process. This research focused on maximizing propeller open water efficiency, while minimizing engine brake power constrained by thrust and cavitation. Optimization was applied to Wageningen B-series propellers and conducted using the Non-Dominated Sorting Genetic Algorithm II (NSGA-II). The algorithm selected optimum parameters to create the propeller model, which was then evaluated numerically through computational fluid dynamics (CFD) with a multiple reference frame (MRF) and under the SST k-ω turbulence model, to obtain the open water hydrodynamic characteristics. In addition, the cavitation effect was evaluated using the Zwart–Gerber–Belamri cavitation model. The numerical model results were validated through comparison with open water experimental data from the Netherlands Ship Model Basin for the Wageningen B-series propellers. The results showed errors of 3.29% and 2.01% in efficiency under noncavitating and cavitating conditions, respectively. Correct performance of the functions was shown, based on neural networks trained to estimate thrust and torque coefficients instead of polynomials. The proposed optimization process and numerical model are suitable for solving multi-objective optimization problems in the preliminary design of fixed-pitch marine propellers.</jats:p>

section_content: The primary goal of marine propeller design is to achieve the highest level of propulsive efficiency. Higher propeller efficiency impacts positively on fuel consumption, decreasing operational costs. Constantly increasing oil prices, in addition to stringent regulations imposed on ships concerning air pollution and restrictions on NO X and SO X emissions, have led the shipping industry to further develop efficient propellers [1] [2] [3] [4]. On the other hand, in order to prevent vibration and sound issues in the hull, it is important to avoid cavitation, despite this being a common occurrence in almost all ship propellers [5, 6]. Thus, it is crucial to maintain low cavitation levels from the initial design phase. Designers often face a complex task in balancing efficiency and cavitation [7]. 

Cavitation occurs when the local pressure drops below the vapor pressure of the liquid at the corresponding temperature. Propeller blades are exposed to local lift forces that can lead to low-pressure peaks. When the surrounding water pressure drops below the vapor pressure, cavities or bubbles will form. Because of their specific geometry and the rotational conditions in which marine propellers operate, cavitation is likely to occur [8, 9]. Cavitation phenomena, when present, can lead to safety problems, major repair and maintenance costs, and additional hydrodynamic energy losses [7, 10]. 

Wageningen B-series propellers [11] have one of the largest and most common datasets available, as a result of experimental data. Propellers analyzed in this series provide parameterization, in order to be able to select a specific propeller among a continuous set of combinations that are valid when properly bounded [12]. In order to estimate a ship's required power and achieve the designed speed and thrust, certain parameters and hydrodynamic characteristics must be matched. These selected parameters include four geometric and one operational parameter: propeller diameter (D), pitch ratio (P/D), expanded area ratio (A E /A 0 ), propeller blade number (Z), and propeller rotational speed (n) [13, 14]. The hydrodynamic characteristics are the dimensionless coefficients: thrust coefficient (K T ), torque coefficient (K Q ), and open water efficiency (η O ). In the conceptual design stage, optimization of the propeller is complex, due to the number of variables and constraints involved, resulting in a multi-objective problem [15]. There is no unique solution, as one solution may be better for one objective but not for others [7]. 

Since the Wageningen B-series propeller dataset has existed, it has formed the basis for a large number of studies, given it provides the experimental data required to validate any theory or mathematical model related to phenomena occurring in a propeller. 

Genetic algorithms (GA) are part of the family of heuristic optimization methods useful for performing optimizations for either single or multi-objective problems. GAs have been applied in the marine industry for solving a variety of problems related to propellers. Lee and Lin [16] implemented a GA to determine the stacking sequence needed to fabricate a composite propeller useful at a wide range of speeds. Pluci ński et al. [17] studied the design of a self-twisting composite propeller with an efficiency optimized for a variety of inflow velocities, using a GA to determine the orientation angles for the fiber in each layer. Gypa et al. [18] used an interactive GA and machine learning approach to optimize the parameters of a propeller based on cavitation constraints; for this purpose, the authors used a vortex lattice method for hydrodynamic evaluation. Ebrahimi et al. [19] carried out a study by applying GA to optimize a baseline propeller, in order to reduce propeller noise, while maximizing hydrodynamic efficiency. They performed experimental tests and validated their results using the Ffowcs Williams-Hawkings (FW-H) model compared to noise measurements. 

GAs applied to the Wageningen B-series propellers have focused on addressing singleobjective or multi-objective optimizations, considering a variety of goals and constraints. Optimization, in most cases, is not shown as a validated process, as the purpose is only to show the behavior and convergence of optimization. Chen and Shih [20] used a GA to optimize the design of a B-series propeller, considering vibration and efficiency as objective functions. Xie [15] implemented the Non-Dominated Sorting Genetic Algorithm II (NSGA-II) to maximize the propeller efficiency and thrust coefficient at a single design speed for a B-series propeller. Esmailian et al. [21] proposed a methodology to optimize the hull-propeller system of a Series-60 ship equipped with a B-series propeller, using NSGA-II with two objective functions given by lifetime fuel consumption (LFC) and a cost function including several propeller parameters. Permadi and Sugianto [14] applied NSGA-II to obtain parameters of three-bladed and five-bladed B-series propellers, in order to optimize thrust and open water efficiency, subject to a thrust constraint. Huisman and Foeth [22] applied NSGA-II to optimize the parametric propeller geometries of two propellers under two case studies, to maximize efficiency under different constraints and minimize cavitation. This study is of great relevance, given it was conducted by the Maritime Research Institute Netherlands (MARIN), who created the original Wageningen B-Series propeller database at the Netherlands Ship Model Basin. 

The present study aimed to investigate the advantages of using the NSGA-II genetic algorithm to optimize the preliminary design of a Wageningen B-series propeller. The objectives were to maximize open water efficiency, while minimizing brake power and inherent fuel consumption, with cavitation and thrust coefficient as the constraints. Convergence of the GA determined the geometry and characteristics of the propeller as follows: J = 0. 73, P/D = 1. 1049, A E /A 0 = 0. 5823, η O = 0. 6363, Z= 5, n = 368 rpm, and P B = 2765 kW. The hydrodynamic characteristics in open water and in cavitation conditions were studied using CFD. Subsequently, validation of the results against experimental data from the Wageningen B-series was performed.
section_name: Introduction

section_name: Theory and Methods

section_content: The overall process for the multi-objective optimization design and validation of a Wageningen B-series propeller is described in Figure 1. 

Maritime Research Institute Netherlands (MARIN), who created the original Wageningen B-Series propeller database at the Netherlands Ship Model Basin. 

The present study aimed to investigate the advantages of using the NSGA-II genetic algorithm to optimize the preliminary design of a Wageningen B-series propeller. The objectives were to maximize open water efficiency, while minimizing brake power and inherent fuel consumption, with cavitation and thrust coefficient as the constraints. Convergence of the GA determined the geometry and characteristics of the propeller as follows: J = 0. 73, P/D = 1. 1049, AE/A0 = 0. 5823, ηO = 0. 6363, Z= 5, n = 368 rpm, and PB = 2765 kW. The hydrodynamic characteristics in open water and in cavitation conditions were studied using CFD. Subsequently, validation of the results against experimental data from the Wageningen B-series was performed.
section_name: Description of the Multi-Objective Optimization Design Process

section_content: The overall process for the multi-objective optimization design and validation of a Wageningen B-series propeller is described in Figure 1. The available data and boundary conditions used to perform this study are given by a case study of an offshore patrol vessel (OPV). A parametric description of the Wageningen B-series propellers is used, given their dataset is publicly available to perform validation of any numerical model. This was a multi-objective study as it was necessary to optimize two functions; open water efficiency was maximized, while power brake was minimized. Constraints for optimization were given as achieving the required thrust while avoiding cavitation. Optimization was based on the NSGA-II algorithm; however, for comparison, a variation of the goal attainment method called minimizing the maximum objective was used, which is also suitable for multi-objective optimization. Both algorithms were implemented, and their results were presented, in order to select the best suboptimal set of parameters to describe the Wageningen B-series propeller, which were then drawn and used to run numerical models in both cavitation and non-cavitation conditions. The final results were compared to experimental data to validate the optimization and numerical model.
section_name: Description of the Multi-Objective Optimization Design Process

section_content: In order to apply our research methodology, a twin-screw offshore patrol vessel (OPV) was used, as shown in Figure 2, and its main characteristics can be found in Table 1. For the purposes of optimization, the propeller diameter was kept constant. The total resistance was obtained through performing an analysis with the Holtrop method at cruising speed, i. e., R T = 163,040. 37 N. 

geningen B-series propellers is used, given their dataset is publicly available to perform validation of any numerical model. This was a multi-objective study as it was necessary to optimize two functions; open water efficiency was maximized, while power brake was minimized. Constraints for optimization were given as achieving the required thrust while avoiding cavitation. Optimization was based on the NSGA-II algorithm; however, for comparison, a variation of the goal attainment method called minimizing the maximum objective was used, which is also suitable for multi-objective optimization. Both algorithms were implemented, and their results were presented, in order to select the best suboptimal set of parameters to describe the Wageningen B-series propeller, which were then drawn and used to run numerical models in both cavitation and non-cavitation conditions. The final results were compared to experimental data to validate the optimization and numerical model.
section_name: Case Study

section_content: In order to apply our research methodology, a twin-screw offshore patrol vessel (OPV) was used, as shown in Figure 2, and its main characteristics can be found in Table 1. For the purposes of optimization, the propeller diameter was kept constant. The total resistance was obtained through performing an analysis with the Holtrop method at cruising speed, i. e., RT = 163,040. 37 N. Twin-screw
section_name: Case Study

section_content: Systematic propeller series based on experimental data from open water model tests are commonly used in the early design stage to determine the optimal geometric parameters. The systematic series used in this project was the common Wageningen B-series. The B-series was first developed by Troost [23] [24] [25], and later improvements were published by Lammeren et al. [11]. Subsequently, Oosterveld and Oossanen [26] fitted the KT-KQ data and presented them as regression polynomials, shown in Equations ( 1 ) and ( 2 ). An alternative approach to computing KT and KQ coefficients based on neural networks is presented in Section 5. 3.
section_name: Wageningen B-Series Propellers

section_content: Systematic propeller series based on experimental data from open water model tests are commonly used in the early design stage to determine the optimal geometric parameters. The systematic series used in this project was the common Wageningen B-series. The B-series was first developed by Troost [23] [24] [25], and later improvements were published by Lammeren et al. [11]. Subsequently, Oosterveld and Oossanen [26] fitted the K T -K Q data and presented them as regression polynomials, shown in Equations ( 1 ) and ( 2 ). An alternative approach to computing K T and K Q coefficients based on neural networks is presented in Section 5. 3. 

Usually, open water test characteristics are presented in dimensionless form and plotted in terms of K T, K Q, and η O on a J-basis [8]. Their respective equations are as follows:

(5)
section_name: Wageningen B-Series Propellers

section_content: Blade preliminary design can be performed using one of the basic procedures such as a Burrill cavitation table or Keller's equation. Keller formulated an equation for obtaining an early estimate of the required expanded area ratio (EAR). This equation provides the minimum required EAR, as shown in Equation ( 7 ) [13, 27]. 

In Equation ( 7 ), T represents the thrust required at the designed speed; p 0 is the absolute pressure located at the shaft's center given by the atmospheric pressure plus the hydrostatic pressure p 0 = p a + ρgh; p v is the vapor pressure given by 3. 170 × 10 -3 MPa at 25 • C according to ITTC; D is the propeller diameter; and K is a constant given by K = 0. 2 for single propeller ships. K = 0. 1 for twin propeller ships. Depth for propeller shaft's center h can be estimated from Equation (8). 

The start of cavitation can be estimated using CFD and analyzing cavitation number by means of the local pressure at the propeller surface. In this sense, it is necessary to model the phase transition [8]. One of the goals during the design stage of hydrodynamic systems is to avoid or at least delay cavitation phenomena. The cavitation number σ is used as an index to quantify the probability of cavitation occurring [13]. Helal et al. [9] described the beginning of cavitation phenomena as cavitation inception. Cavitation number σ and pressure coefficient c p are defined in Equations ( 9 ) and (10) ; and if inequality σ ≤ -c p occurs, then cavitation begins. The higher the cavitation number, the lower the probability cavitation forming [13].
section_name: Cavitarion Criteria

section_content: The genetic algorithm (GA) was created by Holland [28], and it was inspired by Darwin's biological evolution theory. In the beginning, a GA creates a set of chromosomes, which are random initial populations. Then, genetic operators such as selection, crossover, and mutation are used to produce a new generation [29]. As each generation passes, the population produces increasingly fit individuals. This process eventually stops after a certain number of generations. During this process, the fittest individual overall, which represents the global optimum, is identified [17]. 

In particular cases, it is necessary to address a problem with more than one objective function, which leads to a multi-objective optimization problem. NSGA-II is one of the most widely used and generalized methods for solving multi-objective problems. It was proposed by Deb et al. [30] and it introduces the concept of elitism, to improve its convergence properties [31]. 

NSGA-II is an elitist GA that prioritizes individuals with better fitness values and that can increase the population diversity, even if their fitness value is lower. 

Figure 3 presents a flowchart of the general process of NSGA-II. The algorithm begins by creating a random population of solutions named generation 1. The second step is evaluating the fitness functions of the population using corresponding optimization constraints; the results are classified into nondominated layers named Pareto fronts with an assigned rank. The upper layer contains the nondominated or nonsurpassed solutions. In order to enhance the diversity among the population and prevent early convergence, a crowding distance is computed and assigned to each solution. Prioritizing low nondominated ranks with higher crowding distances, selection is performed using tournaments to find the best adapted individuals. Further operations are applied to selected individuals, such as crossover and mutation, with the goal of creating new child solutions. Elitist selection is applied in order to select the best individuals from the current population and pass them unaltered to the next generation. The process of selection, crossover, mutation, and elitism keeps repeating until the convergence criteria or the maximum number of generations is reached. According to Akbari et al. [32], in a Pareto front set of solutions, at least one solution is inferior to another, even though the set as a whole is better than other solutions in the remaining search space. Nevertheless, according to Kramer [33], GAs are excellent for optimizing difficult problems where classical methods fail due to their complexity.
section_name: Non-Dominated Sorting Genetic Algorithm-II (NSGA-II)

section_name: Propeller Optimization

section_content: The first objective function was used to maximize the open water efficiency, as seen According to Akbari et al. [32], in a Pareto front set of solutions, at least one solution is inferior to another, even though the set as a whole is better than other solutions in the remaining search space. Nevertheless, according to Kramer [33], GAs are excellent for optimizing difficult problems where classical methods fail due to their complexity.
section_name: Multi-Objective Optimization Algorithms

section_content: The first objective function was used to maximize the open water efficiency, as seen in Equation (11), from polynomials K T and K Q, in Equations ( 1 ) and ( 2 ), respectively. The second objective function helps the ship attain cruise speed with the minimum amount of power delivered to the propeller (P D ) and brake power (P B ) supplied by the engines. Assuming that there is no friction or transmission losses (η S = 1), the second objective function to be minimized is defined in Equation (12). 

Two explicit constraints are used for optimization. The first is related to Keller's cavitation criterion [27, 34] and can be expressed as the nonlinear inequality shown in Equation (13). The second constraint is the required thrust (T), which is based on the definition of the thrust coefficient (KT); therefore, a nonlinear equality as expressed in Equation ( 14 ) was established. 

Both functions and constraints were implemented in MATLAB ® R2023b software for the NSGA-II, and the minimizing the maximum objective algorithms used the "gamultiobj" and the "fminimax" functions, respectively. The parameters used to run the gamultiobj function are presented in Table 2, while the parameters needed to initialize and run the fminimax function are shown in Table 3.
section_name: Multi-Objective Optimization Algorithms

section_content: Maximum number of evaluations 5000 J 0 0. 4 (P/D) 0 1. 1 (A E /A 0 ) 0 0. 5 Table 4 summarizes the limits for the variables and search space for optimization. The standard range characteristics of the B-series are represented by the limits of J, P/D, and A E /A 0, whereas n covers a reasonable design range and D remains a fixed variable. During execution of the code, the number of blades (Z) was varied between 4 and 5, to cover the most common commercial range. Table 5 presents a comparison of the results obtained from the algorithms for each number of blades. It was observed that the best results were obtained with the NSGA-II algorithm for Z = 5; therefore, this number was used to construct the 3D geometry of the propeller.
section_name: Parameter Value

section_content: The propeller geometry was obtained from the table of general geometric properties for the B-series, as well as equations from Oosterveld and Van Oossanen [26]. Using a standardized parametric propeller series, such as the B-series, has an enormous advantage, as foils to form the blade sections can be generated with coordinates, as depicted in Figure 4. After computing the blade section coordinates, they were exported to create the 3D geometry shown in Figure 5a. Based on the hydrodynamic terms of ITTC [35], the global Cartesian coordinate system (X, Y, Z) of the propeller was defined as outlined in Figure 5b. When looking from aft (-X), a positive rotational speed indicates clockwise rotation and a negative speed indicates counterclockwise rotation, following the right-hand rule.
section_name: Geometric Model

section_name: Numerical Model

section_content: The computational domain of the study used a multiple reference frame (MRF), also known as the frozen rotor approach, as shown in Figure 6. The static or frozen domain was 3 diameters wide and 10 diameters long, to allow open fluid development after it accelerates through the propeller. The rotor, which is situated near the propeller, had a width of 0. 38 diameters and a height of 1. 2 diameters. Considering the size of the rotor domain is crucial, as it impacts the computational cost. If the domain is too large, this will increase the cost; on the other hand, if it is too small, this will result in the formation of large eddies near the propeller, leading to imprecise results. Meanwhile, the "rotor" subdomain was assigned an angular rotational speed n, subject to the advance coefficient J. In the "frozen" domain, a static condition was assigned. According to the MRF approach methodology [36], the interface zone between subdomains was also specified.
section_name: Computational Domain

section_content: Research on turbo-machinery, aircraft propellers, and ship propellers is focused on studying rotational motion. A practical method of representation is using the multiple reference frame (MRF) approach [36], where the geometric model includes a rotational component with respect to a fixed reference. Transitory phenomena are explained through separating the fluid domain into subdomains with at least one fixed subdomain, one rotational subdomain, and interfaces among them.
section_name: Governing Equations

section_content: The computational domain of the study used a multiple reference frame (MRF), also known as the frozen rotor approach, as shown in Figure 6. The static or frozen domain was 3 diameters wide and 10 diameters long, to allow open fluid development after it accelerates through the propeller. The rotor, which is situated near the propeller, had a width of 0. 38 diameters and a height of 1. 2 diameters. Considering the size of the rotor domain is crucial, as it impacts the computational cost. If the domain is too large, this will increase the cost; on the other hand, if it is too small, this will result in the formation of large eddies near the propeller, leading to imprecise results. Meanwhile, the "rotor" subdomain was assigned an angular rotational speed n, subject to the advance coefficient J. In the "frozen" domain, a static condition was assigned. According to the MRF approach methodology [36], the interface zone between subdomains was also specified.
section_name: Numerical Model 4.1. Computational Domain

section_content: Research on turbo-machinery, aircraft propellers, and ship propellers is focused on studying rotational motion. A practical method of representation is using the multiple reference frame (MRF) approach [36], where the geometric model includes a rotational component with respect to a fixed reference. Transitory phenomena are explained through separating the fluid domain into subdomains with at least one fixed subdomain, one rotational subdomain, and interfaces among them. 

MRF uses stationary approximations of the transitory rotation at an instant of time. Flow rotation around moving parts is simulated, adding Coriolis source terms and centripetal accel-erations to the momentum equation. The flow in fixed subdomains is governed by continuity and momentum equations (Equations ( 15 ) and ( 16 )); meanwhile, the rotational subdomains are governed by continuity and momentum equations for the MRF, as seen in Equations ( 17 ) and (18). 

Equations for the MRF can be expressed using relative or absolute velocities as dependent variables [37]. Pressure-based solvers from ANSYS ® Fluent can be configured to use any of those two options. This numerical model was implemented using the absolute velocity.
section_name: Governing Equations

section_content: Reynolds-averaged Navier-Stokes (RANS) turbulence models are a set of two transport equation turbulence models. The kω turbulence model is an empiric model containing the turbulent kinetic energy (k) and the turbulent dissipation frequency (w) that was developed by Wilcox [38]. The kω turbulence model has been improved over time; this is the case of the shear stress transport (SST) kω turbulence model developed by Menter [39]. It adds a laminar boundary layer model that includes a mix function to be applied to the kω and the boundary layer models, to give different responses near the wall and in the freestream depending on the value of the function. 

Even though there are several turbulence models that can be applied to model phenomena happening at the boundary layer, this research used the SST kω turbulence model, as the aforementioned features make it one of the most widely used turbulence models in CFD simulations related to aerodynamics and hydrodynamics, as it provides good results for flows with strong pressure gradients and flow separation for surfaces with curvatures [37]. The SST kω turbulence model is expressed in Equations ( 19 ) and (20).
section_name: Turbulence Model

section_content: As described by Helal et al. [9], flow under cavitation is a two-phase mix made of vapor and liquid that is simplified and simulated as a mix of a single phase. Considering an homogeneous mix, multiphase flow can be solved using modified RANS equations. The continuity and momentum for the mix are described in Equations ( 21 ) and (22). 

In the multiphase approach, density and viscosity are represented as a volume fraction of the secondary phase, as expressed in Equations ( 23 ) and (24), respectively. The liquid-vapor mass transport is governed by the vapor transport in Equation (25). Likewise, the mass transfer between fluid and bubbles is modeled using a derivation of the generalized Rayleigh-Plesset equation, as shown in Equation (26). 

Bubble formation was applied to the numerical model in ANSYS ® Fluent, using the Zwart-Gerber-Belamri cavitation model from Zwart et al. [40]. This cavitation model assumes all bubbles have the same size and is expressed in Equations ( 27 ) and (28).
section_name: Multiphase Model

section_content: Table 6 summarizes the boundary conditions for each variation in J applied to the rotating domain. At the inlet, the flow was homogeneous at a constant advance velocity of 8. 276 m/s. The flow went downstream from bow to stern, with a turbulence intensity of 1%. The static domain of the outer walls were subjected to symmetry conditions, while the propeller walls were established using a no-slip condition. For CFD analysis, two configurations were used to perform different simulations: a noncavitating simulation, and a cavitation behavior simulation. The former was conducted at various advance coefficients J, as previously listed in Table 5. The latter was performed at the optimal operating conditions (Case II). Both configurations involved the use of a C-mesh. The freshwater properties were obtained from ITTC [41] data at 25 • C. Tables 7 and 8 provide a summary of the numerical configurations used.
section_name: CFD Configuration and Boundary Conditions

section_content: Three unstructured meshes (A, B, and C) were generated, varying the number of elements between 2. 21 × 10 6 and 9. 27 × 10 6. The tetrahedron method was applied to the complete domain, with a growth factor of 1. 2. In the critical region near the propeller walls, the inflation technique was applied to increase the resolution and better represent the boundary layer, as shown in Figures 7 and 8.
section_name: Domain Discretization

section_content: Three unstructured meshes (A, B, and C) were generated, varying the number of elements between 2. 21 × 10 6 and 9. 27 × 10 6. The tetrahedron method was applied to the complete domain, with a growth factor of 1. 2. In the critical region near the propeller walls, the inflation technique was applied to increase the resolution and better represent the boundary layer, as shown in Figures 7 and 8. A mesh sensitivity study was conducted to estimate the balance between computational cost and mesh density, obtaining results with a good correlation compared with the experimental data, thus ensuring mesh independence. Metrics for seven meshes were analyzed, in order to meet the criteria established by ANSYS ® Inc. (Canonsburg, PA, USA). [42] of having a good mesh quality, as observed in Table 9. The numerical model was run with every mesh using fresh water in noncavitating conditions (Case II) to calculate the thrust (T), torque (Q), and corresponding open water efficiency (ηO). The results were then compared against experimental B-series data, as presented in Table 10. A mesh sensitivity study was conducted to estimate the balance between computational cost and mesh density, obtaining results with a good correlation compared with the experimental data, thus ensuring mesh independence. Metrics for seven meshes were analyzed, in order to meet the criteria established by ANSYS ® Inc. (Canonsburg, PA, USA). [42] of having a good mesh quality, as observed in Table 9. The numerical model was run with every mesh using fresh water in noncavitating conditions (Case II) to calculate the thrust (T), torque (Q), and corresponding open water efficiency (η O ). The results were then compared against experimental B-series data, as presented in Table 10. The plots in Figures 9 10 11 show the convergence of the thrust, torque, and open water efficiency in the numerical model as the number of elements/nodes was increased to discretize fluid domains using the finite volume method (FMV). As the number of elements/nodes was increased, the convergence of the results in order to make them mesh independent became evident. The mesh selected for the study was mesh F, as it resulted in the first point of the convergence curve after the values had settled.
section_name: Domain Discretization

section_name: Results

section_content: Table 11 and Figures 12 and 13 show a comparison between the numerical model results and the B-series experimental data by means of computing the relative error. In general, a good correlation is considered having a maximum error of less than 10%, according to [43]. The highest error percentage was 8. 32% in the thrust coefficient, while the highest error in the torque coefficient was 6. 23%, both at J = 0. 90. In the case of the open water efficiency, this obtained the lowest percentage error of the set, at J = 0. 90. Figures 14 and 15 illustrate a comparison of the thrust and torques from the numerical models against the expected values from the experimental data, showing a good correlation, with errors not exceeding 10%.
section_name: Noncavitating Condition

section_content: This section shows the results obtained from the numerical model implemented considering cavitation; it used a multiphase mixing model and the cavitation model. After optimization, J = 0. 73 (Case II from Table 5 ) was the optimal point at which the two objective functions were optimized, as power brake was minimal, while
section_name: Cavitating Condition

section_content: This section shows the results obtained from the numerical model implemented considering cavitation; it used a multiphase mixing model and the Zwart-Gerber-Belamri cavitation model. After optimization, J = 0. 73 (Case II from Table 5 ) was the optimal point at which the two objective functions were optimized, as power brake was minimal, while the open water efficiency was the maximum. At the same time, for an optimization that minimizes cavitation as a constraint, the propeller operating point and its characteristics should show minimal or no presence of cavitation. Table 12 shows the results of the numerical model considering cavitation and experimental data from B-series. Likewise, with the noncavitating numerical model, it was observed that the error percentage did not exceed 5%. The best representation of the occurrence of cavitation is the vapor fraction as a result of the numerical model. The proportion of vapor present in the two-phase mixture is visualized with respect to the liquid phase. A value of 1 or close to this number will show that vapor bubbles are fully formed and cavitation is happening. Figure 16 depicts the different values of vapor volume fraction present on both sides of the optimized propeller. In Figure 16a, a low volume fraction at the propeller tips can be observed, while in Figure 16b a wide range of values for the vapor volume fraction can be observed; this means the suction side is the surface of the propeller with the higher probability of having cavitation. The areas with the cooler shades signify a null or lower vapor fraction, whereas warmer colors signify a higher presence of vapor fraction, or in other words cavitation. 

The development of cavitation was further analyzed. The vapor volume iso-surface of fractions 0. 1, 0. 5, and 0. 90 are presented in Figures 17 and 18. From these Figures, it can be observed that high values of vapor volume fraction, i. e., 0. 9 or higher, were present in a minimal portion of the propeller. These contours suggest that if cavitation occurred, it would be of the sheet cavitation type occurring at the suction face. The development of cavitation was further analyzed. The vapor volume iso-surface of fractions 0. 1, 0. 5, and 0. 90 are presented in Figures 17 and 18. From these Figures, it can be observed that high values of vapor volume fraction, i. e., 0. 9 or higher, were present in a minimal portion of the propeller. These contours suggest that if cavitation occurred, it would be of the sheet cavitation type occurring at the suction face. The development of cavitation was further analyzed. The vapor volume iso-surface of fractions 0. 1, 0. 5, and 0. 90 are presented in Figures 17 and 18. From these Figures, it can be observed that high values of vapor volume fraction, i. e., 0. 9 or higher, were present in a minimal portion of the propeller. These contours suggest that if cavitation occurred, it would be of the sheet cavitation type occurring at the suction face.
section_name: Cavitating Condition

section_content: Polynomials for retrieving KT and KQ coefficients for the Wageningen B-series propeller were first introduced by Oosterveld and Oossanen [26] ; since then they have been used as the primary means to obtain these values, together with corresponding plots. Comparatively, in a similar way to how polynomials were developed through multiple regressions in order to fit experimental data, nowadays neural networks (NN) can be trained with a set of data points in order to predict the same coefficients. Data from [26, 44] were used to train two independent neural networks using MATLAB ® with the parameters seen in Table 13. The performance in training neural networks for KT and KQ is shown in Figures 19 and 20, respectively; end of training is marked with the green circle.
section_name: Approximation of KT and KQ Coefficients Using Neural Networks

section_content: Polynomials for retrieving K T and K Q coefficients for the Wageningen B-series propeller were first introduced by Oosterveld and Oossanen [26] ; since then they have been used as the primary means to obtain these values, together with corresponding plots. Comparatively, in a similar way to how polynomials were developed through multiple regressions in order to fit experimental data, nowadays neural networks (NN) can be trained with a set of data points in order to predict the same coefficients. Data from [26, 44] were used to train two independent neural networks using MATLAB ® with the parameters seen in Table 13. The performance in training neural networks for K T and K Q is shown in Figures 19 and 20, respectively; end of training is marked with the green circle. In order to compare the results from using polynomials and neural networks, coefficients and efficiencies were plot corresponding to Wageningen B-series for Z, AE/A0, and P/D using the suboptimal parameters obtained from the NSGA-II algorithm. Figure 21 shows a plot using polynomials, Figure 22 shows a plot using neural networks, and Figure 23 depicts the absolute error between the polynomial and neural network estimations for KT and KQ. From Figure 23, it can be observed that errors fell within a range of 10 -4 in this specific case. In order to compare the results from using polynomials and neural networks, coefficients and efficiencies were plot corresponding to Wageningen B-series for Z, AE/A0, and P/D using the suboptimal parameters obtained from the NSGA-II algorithm. Figure 21 shows a plot using polynomials, Figure 22 shows a plot using neural networks, and Figure 23 depicts the absolute error between the polynomial and neural network estimations for KT and KQ. From Figure 23, it can be observed that errors fell within a range of 10 -4 in this specific case. In order to compare the results from using polynomials and neural networks, coefficients and efficiencies were plot corresponding to Wageningen B-series for Z, A E /A 0, and P/D using the suboptimal parameters obtained from the NSGA-II algorithm. Figure 21 shows a plot using polynomials, Figure 22 shows a plot using neural networks, and Figure 23 depicts the absolute error between the polynomial and neural network estimations for K T and K Q. From Figure 23, it can be observed that errors fell within a range of 10 -4 in this specific case.
section_name: Approximation of K T and K Q Coefficients Using Neural Networks

section_content: In this research, the benefit of using the NSGA-II genetic algorithm to optimize the preliminary design of a Wageningen B-series propeller was studied. Among the many possible configurations, the algorithm selected one that increased the open water efficiency and minimized the ship brake power within the thrust and cavitation constraints.
section_name: Conclusions

section_content: Acknowledgments: A. V. -S. extends her gratitude to the Consejo Nacional de Humanidades, Ciencias y Tecnologías for the scholarship granted during her master's studies.

section_content: Data Availability Statement: Data are contained within the article.

section_name: Conclusions

section_content: In this research, the benefit of using the NSGA-II genetic algorithm to optimize the preliminary design of a Wageningen B-series propeller was studied. Among the many possible configurations, the algorithm selected one that increased the open water efficiency and minimized the ship brake power within the thrust and cavitation constraints. CFD simulations using the SST k-ω turbulence model and MRF approach allowed the analysis of the hydrodynamic characteristics in noncavitating and cavitation conditions. Based on the numerical results obtained, the following conclusions were reached:

• The numerical model indicated a close correlation between the numerical predictions and experimental data in noncavitating conditions. With maximum errors of 8. 32%, 6. 23%, and 3. 29% for K T, K Q, and η O, respectively. • The use of a multiphase approach and the Zwart-Gerber-Belamri cavitation model enabled a detailed investigation into the propeller cavitation phenomenon. The numerical model revealed that the propeller underwent partial cavitation on its suction side, while also revealing a 2. 01% error between the efficiency η O of the numerical model's operating point and the experimental data.
section_name: Conclusions

section_content: Performing mesh refinement, particularly in the areas of interest, such as the tips and leading edge of the blades, could improve the capture of tip and sheet type cavitation effects in future research.
section_name: •

section_content: Functions derived from trained neural networks could be used for estimating K T and K Q coefficients instead of the original polynomials. 

The presented methodology was to be feasible for addressing complex and multivariable problems in the field of naval engineering by combining two advanced techniques: genetic algorithm optimization and CFD analysis.
section_name: •

section_content: The APC was funded by the Master in Applied Engineering of Universidad Veracruzana.
section_name: Funding:

section_content: The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.
section_name: Conflicts of Interest: