Two-Objective Optimization of a Kaplan Turbine Draft Tube Using a Response Surface Methodology

: The overall cost of a hydropower plant is mainly due to the expenses of civil works, mechanical equipment (turbine and control units) and electrical components. The goal of a new draft tube design is to obtain a geometry that reduces investment costs, especially the excavation ones, but the primary driver is to increase overall machine e ﬃ ciency, allowing for a reduced payback time. In the present study, an optimization study of the elbow-draft tube assembly of a Kaplan turbine was conducted. First, a CFD model for the complete turbine was developed and validated. Next, an optimization of the draft tube alone was performed using a design of experiments technique. Finally, several optimum solutions for the draft tube were obtained using a response surface technique aiming at maximizing pressure recovery and minimizing ﬂow losses. A selection of optimized geometries was subsequently post-checked using the validated model of the entire turbine, and a detailed ﬂow analysis on the obtained results made it possible to provide insights into the improved designs. It was observed that e ﬃ ciency could be improved by 1% (in relative terms), and the mechanical power increased by 1.8% (in relative terms) with respect to the baseline turbine. order to assure improved numerical stability, 50 iterations were run at first using first-order schemes. These results were adopted as initial values for the high-resolution scheme. During the convergence runs, several variables of interest were monitored, such as flow rate, efficiency, static pressure at draft tube inlet, and mechanical power output. Final results were available after 3000 iterations on average and almost 72 h of CPU calculations. All calculations were performed on a server equipped with Intel ® Xeon ® CPU x5650 processors using parallel solutions on 10 multiple cores. Convergence was assessed when all variables of interest showed a variation lower than 0.08%. a multi-objective optimization based on the construction of a response surface. The shape


Introduction
In recent decades, environmental policies have been oriented towards reducing energetic dependence on fossil fuels, leading to renewed and increased investments in the hydroelectric sector; such investments have been primarily directed to optimizing both new designs and existing installations, thus allowing for more efficient plants and reduced payback times. Kaplan hydro turbines have been adopted for a long time to deal with efficient energy production in the range of high specific speeds [1]. The use of movable vaned distributors and runner blades, in fact, allows to limit incidence flow losses and achieve high conversion efficiency over a wide range of flow rates, from 40% to 120% of the nominal value. The main drawback of Kaplan turbines is that their size increases as the flow rate gets higher, and the cost attributed to civil works becomes more and more substantial.
In these machines, the draft tube is an essential component as it makes it possible to recover a significant fraction of the kinetic energy leaving the runner by conversion into static pressure. This, in turn, leads to creating an additional suction head downstream of the runner, which allows to increase the effective head that the runner is able to deliver [2]. In addition, the draft tube plays a fundamental role in determining turbine efficiency since the height at which the tube is installed contributes to a large percentage of the total net head that can be recovered [3].
The datum, or baseline, draft tube was connected to a Kaplan turbine runner installed in a 576 kW hydraulic power plant located in the Adda river (Italy), as depicted in Figure 1a. The original geometry adopted for the elbow draft tube is schematically represented in Figure 1b.
As mentioned above, the draft tube behaves as a diffuser acting to recover part of the kinetic energy from Section 2 to 3, according to its prescribed area distribution along the equivalent centerline. Two coefficients are used here to quantify the performance of a draft tube, namely a pressure recovery factor (Cp) and a drag coefficient (Cd), respectively defined as: where p is the static pressure and ptot the total pressure; subscripts 2 and 3 refer to the station downstream of the runner and on the draft tube outlet, respectively; ρ is the water density; c is the absolute velocity. The values of these coefficients are strictly related to the area distribution and are notoriously conflicting in such a way that the more intense the diffusion is (high Cp), the bigger the head loss is (high Cd) that might occur as a result of flow separation, secondary losses, etc. Therefore, a good draft tube features the highest possible Cp along with the minimum achievable Cd. These two coefficients were used as objective functions in the multi-objective optimization study described later.

CFD Model Setup and Validation
The steady-state RANS CFD model of the complete turbine adopted included complete spiral casing with stay vanes, guide vanes, and runner and draft tubes ( Figure 2). The final grids were chosen based on a mesh sensitivity analysis and their details and statistics are provided in Table 1. The influence of mesh size on the draft tube performance parameters is shown in Figure 3, where a limited sensitivity to the grid size can be observed, starting from the third level, having 1.8 M cells. As for boundary conditions, a total pressure was defined at the inlet of the spiral casing and a static pressure equal to 0 Pa (relative to the atmospheric pressure value of 1 atm) was set at the outlet of the draft tube. As mentioned above, the draft tube behaves as a diffuser acting to recover part of the kinetic energy from Section 2 to 3, according to its prescribed area distribution along the equivalent centerline.
Two coefficients are used here to quantify the performance of a draft tube, namely a pressure recovery factor (Cp) and a drag coefficient (Cd), respectively defined as: where p is the static pressure and p tot the total pressure; subscripts 2 and 3 refer to the station downstream of the runner and on the draft tube outlet, respectively; ρ is the water density; c is the absolute velocity. The values of these coefficients are strictly related to the area distribution and are notoriously conflicting in such a way that the more intense the diffusion is (high Cp), the bigger the head loss is (high Cd) that might occur as a result of flow separation, secondary losses, etc. Therefore, a good draft tube features the highest possible Cp along with the minimum achievable Cd. These two coefficients were used as objective functions in the multi-objective optimization study described later.

CFD Model Setup and Validation
The steady-state RANS CFD model of the complete turbine adopted included complete spiral casing with stay vanes, guide vanes, and runner and draft tubes ( Figure 2). The final grids were chosen based on a mesh sensitivity analysis and their details and statistics are provided in Table 1. The influence of mesh size on the draft tube performance parameters is shown in Figure 3, where a limited sensitivity to the grid size can be observed, starting from the third level, having 1.8 M cells. As for boundary conditions, a total pressure was defined at the inlet of the spiral casing and a static pressure equal to 0 Pa (relative to the atmospheric pressure value of 1 atm) was set at the outlet of the draft tube.
The total pressure at the inlet boundary corresponded to the total energy that the plant could process and was calculated from the net head and the flow rate that the turbine was subjected to. The flow solver was Ansys CFX 19.2. Since the runner was not modelled as a rotating mesh, Energies 2020, 13, 4899 4 of 13 a multiple-reference-frame approach was chosen, and a mixing-plane condition was set both at the distributor-runner and runner-draft tube interfaces. Among all the turbulence models (all the variants of both k-ε and k-ω available in Ansys CFX 19.2 were actually tested), it was observed that the k-ω SST [15] provided the most accurate predictions, but it was also the one that required the highest computational effort; for this reason, the k-ω SST was used in off-loop analyses on the complete machine model, while the standard k-ε was used in the optimization due to its cheaper requirement. A physical timescale corresponding to 1 • of runner rotation and a high-resolution scheme were set.
In order to assure improved numerical stability, 50 iterations were run at first using first-order schemes. These results were adopted as initial values for the high-resolution scheme.     The total pressure at the inlet boundary corresponded to the total energy that the plant could process and was calculated from the net head and the flow rate that the turbine was subjected to. The flow solver was Ansys CFX 19.2. Since the runner was not modelled as a rotating mesh, a multiplereference-frame approach was chosen, and a mixing-plane condition was set both at the distributorrunner and runner-draft tube interfaces. Among all the turbulence models (all the variants of both kε and k-ω available in Ansys CFX 19.2 were actually tested), it was observed that the k-ω SST [15] provided the most accurate predictions, but it was also the one that required the highest computational effort; for this reason, the k-ω SST was used in off-loop analyses on the complete machine model, while the standard k-ε was used in the optimization due to its cheaper requirement. A physical timescale corresponding to 1° of runner rotation and a high-resolution scheme were set.
In order to assure improved numerical stability, 50 iterations were run at first using first-order schemes. These results were adopted as initial values for the high-resolution scheme.
During the convergence runs, several variables of interest were monitored, such as flow rate, efficiency, static pressure at draft tube inlet, and mechanical power output. Final results were available after 3000 iterations on average and almost 72 h of CPU calculations. All calculations were performed on a server equipped with Intel ® Xeon ® CPU x5650 processors using parallel solutions on 10 multiple cores. Convergence was assessed when all variables of interest showed a variation lower than 0.08%. During the convergence runs, several variables of interest were monitored, such as flow rate, efficiency, static pressure at draft tube inlet, and mechanical power output. Final results were available after 3000 iterations on average and almost 72 h of CPU calculations. All calculations were performed on a server equipped with Intel ® Xeon ® CPU x5650 processors using parallel solutions on 10 multiple cores. Convergence was assessed when all variables of interest showed a variation lower than 0.08%.
Results from both experiments and simulations are presented in Figure 4. All the experimental data (net head, flow rate, efficiency, and mechanical power) were taken in compliance with the European IEC EN 60041:1991 standard [16] and included the appropriate uncertainty bandwidth. Power and efficiency data from CFD were calculated as follows: where M is the total torque on the runner blades, ω is the rotational speed, Q is the flow rate, and H is the net head. As can be observed in Figure 4, the power and efficiency data were within the experimental measurement uncertainty range. The steady RANS CFD model set up appeared, therefore, to be able to reproduce the machine performance with sufficient accuracy and acceptable computational time. From this validated model, the runner-draft tube interface conditions, that is, the distributions of the u, v, w (respectively tangential, absolute, and relative velocity components), were extracted to be used as inlet boundary conditions during the optimization process.
Energies 2020, 13, x FOR PEER REVIEW 6 of 14 Results from both experiments and simulations are presented in Figure 4. All the experimental data (net head, flow rate, efficiency, and mechanical power) were taken in compliance with the European IEC EN 60041:1991 standard [16] and included the appropriate uncertainty bandwidth. Power and efficiency data from CFD were calculated as follows: where M is the total torque on the runner blades, ω is the rotational speed, Q is the flow rate, and H is the net head. As can be observed in Figure 4, the power and efficiency data were within the experimental measurement uncertainty range. The steady RANS CFD model set up appeared, therefore, to be able to reproduce the machine performance with sufficient accuracy and acceptable computational time. From this validated model, the runner-draft tube interface conditions, that is, the distributions of the u, v, w (respectively tangential, absolute, and relative velocity components), were extracted to be used as inlet boundary conditions during the optimization process.

Draft Tube Geometry and Parameterization
The draft tube geometry of Kapan turbines is complicated by the different shapes of its inlet and outlet sections. The first is circular due to the interface with the runner, while the second is often rectangular. Moreover, the draft tube features a 90° elbow to minimize excavation costs and to

Draft Tube Geometry and Parameterization
The draft tube geometry of Kapan turbines is complicated by the different shapes of its inlet and outlet sections. The first is circular due to the interface with the runner, while the second is often rectangular. Moreover, the draft tube features a 90 • elbow to minimize excavation costs and to improve powerhouse compactness. As a result, a large number of design parameters is necessary to provide successful geometry parameterization. In the present paper, the draft tube was parameterized using a mean line and a number of cross-sections stacked along it.
The mean line was designed to be composed by a first straight section, made up of two segments, one related to the divergent cone, and one to the cylindrical section; a second curvilinear line was related to the elbow shape and controlled by a three-point Bézier polygon; finally, a segment was connected to the last section of the draft tube.
The first straight part of the mean line was described by two parameters-the two segments' lengths. For the Bézier polygon, four additional variables were introduced-the coordinates of the polygon-points (the coordinates of the first Bézier polygon point were derived from the first straight part). The last parameter was the length of the last straight segment. In total, 7 variables were involved in the parameterization of the mean line. The area distribution was managed using a 4th-degree polynomial as a function of the mean line coordinate x: In Figure 5, the area distribution of the original geometry is represented using dots along with the polynomial interpolation. From the original geometry, it was possible to obtain the polynomial coefficients (a, b, c, d, e) and use P1 and P2 as decision variables.
Energies 2020, 13, x FOR PEER REVIEW 7 of 14 In Figure 5, the area distribution of the original geometry is represented using dots along with the polynomial interpolation. From the original geometry, it was possible to obtain the polynomial coefficients (a, b, c, d, e) and use P1 and P2 as decision variables. Overall, the draft tube was parameterized using 9 decision variables. In order to avoid the need for further parameters to account for the shape change from inlet to exit sections, the geometric similarity with the baseline tube was preserved, keeping the geometric ratios R1/L, R2/L, and H/L for each section along the mean line equal to those of the datum.

Optimization: Problem Formulation and Tools
As previously mentioned, the optimization study was conducted on the draft tube alone by taking the inlet boundary conditions from the full machine model at the runner-draft tube interface. As a result, the flow domain included the runner discharge cone along with the runner tube interface (which form the inlet boundary), the first divergent part, the elbow, the last divergent part, and the draft tube prolongation leading to the outlet section.
The optimization was multi-objective, with the purpose of maximizing Cp and minimizing Cd at the same time using a Pareto approach. As for the constraints, the optimization had to lead to solutions that were interchangeable with the baseline case currently installed without any modifications regarding the runner. As a result, during the entire procedure, the baseline inlet and runner cone (both belonging to the runner) remained fixed. A further constraint was related to the Overall, the draft tube was parameterized using 9 decision variables. In order to avoid the need for further parameters to account for the shape change from inlet to exit sections, the geometric similarity with the baseline tube was preserved, keeping the geometric ratios R1/L, R2/L, and H/L for each section along the mean line equal to those of the datum.

Optimization: Problem Formulation and Tools
As previously mentioned, the optimization study was conducted on the draft tube alone by taking the inlet boundary conditions from the full machine model at the runner-draft tube interface. As a result, the flow domain included the runner discharge cone along with the runner tube interface (which form the inlet boundary), the first divergent part, the elbow, the last divergent part, and the draft tube prolongation leading to the outlet section.
The optimization was multi-objective, with the purpose of maximizing Cp and minimizing Cd at the same time using a Pareto approach. As for the constraints, the optimization had to lead to solutions that were interchangeable with the baseline case currently installed without any modifications regarding the runner. As a result, during the entire procedure, the baseline inlet and runner cone (both belonging to the runner) remained fixed. A further constraint was related to the global draft tube depth which affects the cost of civil works. The outlet depth could vary in all directions but with the limit of always remaining with a certain margin under the free water surface of the tailrace.
A response surface methodology for the optimization tool was used [17]. The first step was a design space exploration using design of experiments (DOE) in order to understand how the design parameters were related to each other and which were the most significant ones. A "central composite design" (CCD) DOE technique was adopted, which provided a screening set of 150 samples to determine the overall trends of the metamodel. The latter was built using the genetic aggregation method [18] and the Pareto optimal solutions were finally obtained (Figure 6a). A check-in-the-loop of optimal solutions was necessary to evaluate the response surface goodness using CFD (Figure 6b); the best-so-far optimum individuals from the response surface were taken and verified using a dedicated CFD simulation, then inserted back in the DOE, which, in turn, updated the response surface from which a new Pareto front was obtained. The loop was kept running until a stable and invariant Pareto optimal set was determined.
Energies 2020, 13, x FOR PEER REVIEW 8 of 14 It is worth mentioning that during the loop iterations, the response surface method allowed for an adaptive refinement in the search regions around maximum CP and minimum CD so as to make it possible to have a wider and more uniform Pareto front.

Results and Discussion
The final Pareto front is depicted in Figure 7, from which a subset of relevant configurations could be extracted: the baseline solution is marked with a black diamond, while the optimum CD, the optimum CP, and a good compromise between the CD and CP are marked using a red triangle, a yellow circle, and a green square, respectively. Dominated individuals are indicated as other points in the figure. The performance data of optimized individuals are given in Table 2. The solution corresponding to the minimum CD (Candidate 1) featured a reduction count ΔCD = −0.0254 (i.e., −2.54 percentage point reduction), while the one corresponding to the maximum CP (Candidate 3) showed ΔCP = +0.0574 (i.e., +5.7 percentage point improvement). It is worth mentioning that during the loop iterations, the response surface method allowed for an adaptive refinement in the search regions around maximum C P and minimum C D so as to make it possible to have a wider and more uniform Pareto front.

Results and Discussion
The final Pareto front is depicted in Figure 7, from which a subset of relevant configurations could be extracted: the baseline solution is marked with a black diamond, while the optimum C D , the optimum C P , and a good compromise between the C D and C P are marked using a red triangle, a yellow circle, and a green square, respectively. Dominated individuals are indicated as other points in the figure.
The performance data of optimized individuals are given in Table 2. The solution corresponding to the minimum C D (Candidate 1) featured a reduction count ∆C D = −0.0254 (i.e., −2.54 percentage point reduction), while the one corresponding to the maximum C P (Candidate 3) showed ∆C P = +0.0574 (i.e., +5.7 percentage point improvement).
Energies 2020, 13, 4899 8 of 13 could be extracted: the baseline solution is marked with a black diamond, while the optimum CD, the optimum CP, and a good compromise between the CD and CP are marked using a red triangle, a yellow circle, and a green square, respectively. Dominated individuals are indicated as other points in the figure. The performance data of optimized individuals are given in Table 2. The solution corresponding to the minimum CD (Candidate 1) featured a reduction count ΔCD = −0.0254 (i.e., −2.54 percentage point reduction), while the one corresponding to the maximum CP (Candidate 3) showed ΔCP = +0.0574 (i.e., +5.7 percentage point improvement).

Optimal Draft Tube Configurations
In Figure 8, the area distribution along the mean line coordinate for the Pareto optimal subset of solutions compared to the baseline one is shown. It can be seen that optimal solutions exhibited a much lower increase in area distribution from the beginning of the draft the tube (mean line coordinate > 0.05, see Figure 8b); this is beneficial to limit adverse phenomena during flow diffusion in this zone as confirmed, also, by the distribution of Cd versus the mean line coordinate (Figure 10a), particularly visible in the case of Candidate 1.

Optimal Draft Tube Configurations
In Figure 8, the area distribution along the mean line coordinate for the Pareto optimal subset of solutions compared to the baseline one is shown. It can be seen that optimal solutions exhibited a much lower increase in area distribution from the beginning of the draft the tube (mean line coordinate > 0.05, see Figure 8b); this is beneficial to limit adverse phenomena during flow diffusion in this zone as confirmed, also, by the distribution of Cd versus the mean line coordinate (Figure 10a), particularly visible in the case of Candidate 1. In Figure 9, a detailed comparison between the optimal solutions and the baseline is presented. For the first draft tube part (until the 20% of the mean line), very similar characteristics can be observed, while the most remarkable differences are visible in the elbow and in the last part.
Candidate 1 had a very similar geometry compared to the baseline, while Candidate 2 featured an elbow with a much higher radius of curvature. Candidate 3, on the other hand, had much larger dimensions than the baseline and this has a major impact on the excavation costs. As it will be seen In Figure 9, a detailed comparison between the optimal solutions and the baseline is presented. For the first draft tube part (until the 20% of the mean line), very similar characteristics can be observed, while the most remarkable differences are visible in the elbow and in the last part. In Figure 9, a detailed comparison between the optimal solutions and the baseline is presented. For the first draft tube part (until the 20% of the mean line), very similar characteristics can be observed, while the most remarkable differences are visible in the elbow and in the last part. Candidate 1 had a very similar geometry compared to the baseline, while Candidate 2 featured an elbow with a much higher radius of curvature. Candidate 3, on the other hand, had much larger dimensions than the baseline and this has a major impact on the excavation costs. As it will be seen later, Candidate 2 was the best in performance and this is due to the shape of the elbow, a fact that emphasizes the importance of producing a good design in the region where a big part of the diffusion is realized by the draft tube.  Candidate 1 had a very similar geometry compared to the baseline, while Candidate 2 featured an elbow with a much higher radius of curvature. Candidate 3, on the other hand, had much larger dimensions than the baseline and this has a major impact on the excavation costs. As it will be seen later, Candidate 2 was the best in performance and this is due to the shape of the elbow, a fact that emphasizes the importance of producing a good design in the region where a big part of the diffusion is realized by the draft tube.

Post-Check Validation on Full Machine and Result Transposition
A number of Pareto optimal draft tube geometries were obtained from the optimization and were subsequently analyzed, a posteriori, by connecting them into the CFD model of the entire machine in order to assess their influence on critical plant characteristics, such as mechanical power and efficiency. In Figure 10, both the C d and C p of the optimal draft tube geometries computed using the full CFD model of the turbine are plotted as functions of the mean line coordinate.

Post-Check Validation on Full Machine and Result Transposition
A number of Pareto optimal draft tube geometries were obtained from the optimization and were subsequently analyzed, a posteriori, by connecting them into the CFD model of the entire machine in order to assess their influence on critical plant characteristics, such as mechanical power and efficiency. In Figure 10, both the Cd and Cp of the optimal draft tube geometries computed using the full CFD model of the turbine are plotted as functions of the mean line coordinate. The CFD analyses performed on the full machine returned different values for the turbine performance figures, especially for the net head and mass flow rate (Table 3). For a better comparison, a result transposition was carried out according to the European IEC EN 60041:1991 standard [16]. In particular, the mass flow rate and the power output data were transposed using the baseline net head. This was acceptable since the condition 0.99 < √H/√H' < 1.01 on head H' compared to baseline H was met for each candidate. Optimal candidate solutions led to greater power and efficiency values compared to the baseline. In particular, Candidate 2 featured an increment of 1.8% on the produced power and 1% in efficiency ratio.  The CFD analyses performed on the full machine returned different values for the turbine performance figures, especially for the net head and mass flow rate (Table 3). For a better comparison, a result transposition was carried out according to the European IEC EN 60041:1991 standard [16]. In particular, the mass flow rate and the power output data were transposed using the baseline net head. This was acceptable since the condition 0.99 < √ H/ √ H' < 1.01 on head H' compared to baseline H was met for each candidate. Optimal candidate solutions led to greater power and efficiency values compared to the baseline. In particular, Candidate 2 featured an increment of 1.8% on the produced power and 1% in efficiency ratio.  From a fluid dynamic standpoint, it is apparent that the largest part of the pressure recovery took place in the first part of the diffuser (Figure 10b), which confirms the importance of having a very accurate design of the region upstream of the elbow. From a geometrical point of view (Figure 9), Candidates 1 and 2 were not very different from the baseline, although they clearly exhibited an elbow shape having a larger radius of curvature, which eventually caused the exit stations of both the elbows and the draft tubes to be placed further downstream with respect to the baseline. Additionally, Candidate 3 featured the same global tendencies but with a much longer draft tube, which, eventually, would lead to higher installation costs. In Figure 11 a 3D view of the several solutions is presented, with a comparison of the several streamlines. It can be seen that candidate 1 has the most uniform flow, a fact confirmed by his C d value (Figure 9).
From a fluid dynamic standpoint, it is apparent that the largest part of the pressure recovery took place in the first part of the diffuser (Figure 10b), which confirms the importance of having a very accurate design of the region upstream of the elbow.
As a matter of fact, approximately 70% of the pressure recovery is realized in the first 20% length of the diffuser. This also held true for the total pressure losses which were primarily generated in this region (see C d behavior in Figure 10, where the largest slope is evidenced in the C d distribution in the first 20% length of the baseline geometry); localized losses accumulated all along the draft tube mean line, finally leading to an overall large C d value in the baseline. Such an observation confirms what was previously pointed out in [5] and [13]. Furthermore, by looking at the total pressure distribution on consecutive station cuts along the draft tube ( Figure 12, right column), it is evident that the rotational region at the lowest total pressure in the core of the discharge flow was of a much lower intensity in the optimized solutions compared to the baseline, despite the latter having a separation wall in the middle of the diffuser. In fact, the analyses carried out showed that the vortex rope greatly influences the draft tube's performance, creating instabilities and turbulent zones in the flow field.
Energies 2020, 13, x FOR PEER REVIEW 12 of 14 As a matter of fact, approximately 70% of the pressure recovery is realized in the first 20% length of the diffuser. This also held true for the total pressure losses which were primarily generated in this region (see Cd behavior in Figure 10, where the largest slope is evidenced in the Cd distribution in the first 20% length of the baseline geometry); localized losses accumulated all along the draft tube mean line, finally leading to an overall large Cd value in the baseline. Such an observation confirms what was previously pointed out in [5] and [13]. Furthermore, by looking at the total pressure distribution on consecutive station cuts along the draft tube ( Figure 12, right column), it is evident that the rotational region at the lowest total pressure in the core of the discharge flow was of a much lower intensity in the optimized solutions compared to the baseline, despite the latter having a separation wall in the middle of the diffuser. In fact, the analyses carried out showed that the vortex rope greatly influences the draft tube's performance, creating instabilities and turbulent zones in the flow field.

Conclusions
A validated CFD model of a draft tube of a Kaplan turbine was successfully implemented and used for a multi-objective optimization based on the construction of a response surface. The shape optimization involved only the draft tube domain, where a one-way coupling with the runner was included in the form of a prescribed velocity distribution extracted from the full-component baseline model.
Two Pareto optimal designs have been extracted from the final set which outperformed the baseline. Regarding the performance of the draft tube alone, the optimum CD featured a reduction count of -0.0254 (-2.54 percentage point reduction), while the optimum CP showed an incremental

Conclusions
A validated CFD model of a draft tube of a Kaplan turbine was successfully implemented and used for a multi-objective optimization based on the construction of a response surface. The shape optimization involved only the draft tube domain, where a one-way coupling with the runner was included in the form of a prescribed velocity distribution extracted from the full-component baseline model.
Two Pareto optimal designs have been extracted from the final set which outperformed the baseline. Regarding the performance of the draft tube alone, the optimum C D featured a reduction count of -0.0254 (-2.54 percentage point reduction), while the optimum C P showed an incremental count of +0.0574 (i.e., +5.7 percentage point improvement) with respect to the baseline geometry.
A post-check carried out including the optimized draft tubes in the overall turbine confirmed the optimization trends, although the best improvements in terms of delivered hydraulic power-in the case where the overall turbine was considered-have been registered, including a third type of diffuser which exhibited an increment of 1.8% in the power produced.
The three best candidates overall were slightly deeper than the baseline in terms of excavation dimensions, so a cost/benefit assessment will be needed to determine which is the most appropriate draft tube for a given plant. In addition, future work will be related to extending the investigation to more operating points by performing a multi-objective multi-point optimization, where the Pareto solutions will account for the different conditions at the tube inflow along the turbine operating line.