Research Article  Open Access
Multipoint and Multiobjective Optimization of a Centrifugal Compressor Impeller Based on Genetic Algorithm
Abstract
The design of high efficiency, high pressure ratio, and wide flow range centrifugal impellers is a challenging task. The paper describes the application of a multiobjective, multipoint optimization methodology to the redesign of a transonic compressor impeller for this purpose. The aerodynamic optimization method integrates an improved nondominated sorting genetic algorithm II (NSGAII), blade geometry parameterization based on NURBS, a 3D RANS solver, a selforganization map (SOM) based data mining technique, and a time series based surge detection method. The optimization results indicate a considerable improvement to the total pressure ratio and isentropic efficiency of the compressor over the whole design speed line and by 5.3% and 1.9% at design point, respectively. Meanwhile, surge margin and choke mass flow increase by 6.8% and 1.4%, respectively. The mechanism behind the performance improvement is further extracted by combining the geometry changes with detailed flow analysis.
1. Introduction
The reliance on numerical methods in the aerodynamic design process of turbomachinery components has considerably increased in the last decades. Nowadays, Computational Fluid Dynamics (CFD) codes have matured to a level where they are capable of not only providing a substantial insight into the threedimensional flow field in turbomachines, but also calculating aerodynamic performances of the machines [1–3]. Meanwhile, higher pressure ratio demand, impeller efficiency, and compressor map width are in severe tradeoff relations for transonic impellers [4]. Essentially, the aerodynamic design of transonic impellers is a multiobjective problem. It is a challenging task to design a high pressure ratio centrifugal impeller for high efficiency and wide operating range at the same time [5].
In order to obtain better designs and reduce design cost, automated design optimization of centrifugal impeller has received a widespread attention in recent years. Guo et al. [6] conducted an automated design optimization of a high pressure ratio centrifugal impeller by integrating an evolution algorithm, 3D blade parameterization method, CFD solver technique, and data mining technique. Verstraete et al. [7] combined a genetic algorithm with an artificial neural network (ANN) to optimize a centrifugal compressor. HyunSu et al. [8] carried out the optimal design of impeller for a centrifugal compressor under the influence of flowinduced vibration using fluidstructure interaction and response surface method (RSM). Although much progress has been made in this area, the research on the multiobjective design optimization of a blade is still insufficient especially for very high pressure centrifugal impellers, which has to take into account the tradeoff among objective functions. In addition, most published optimization design techniques were performed at one single operating point, usually design point, with the danger of serious deterioration of the performance at offdesign conditions such as poorer surge margin and smaller swallowing capacity.
The optimization for the whole speed line is even more challenging, especially near surge condition (certain impellers require good efficiency near the surge and most compressors need good surge margin), as CFD convergence is often not guaranteed under or near such a condition. Demeulenaere et al. [9] performed a multipoint optimization of a turbocharger compressor wheel. They simulated entire compressor stage including diffuser and housing with significantly increased computational time and effort. Their criterion of surge improvement is higher pressure ratios and the method is open to questions. Pini et al. [10] did a shape optimization of a supersonic turbine cascade at offdesign conditions. These two studies propose a pseudoobjective function by summing up all the penalty terms and the original performance objectives with weighting factors, which is not strictly multiobjectives at multipoints. The design variables are so strong and complex to affect the performances that the more variables must be considered in the optimization design. However, there are less analytical expressions available to directly correlate the design variables with the performance at present. Data mining techniques are considered to be able to provide a possible way to extract some useful information from the design space and make the optimization problems in an accessible way. By detecting the features of the data set such as parameter correlations, data mining method can help to gain the mechanism over the performance improvement by optimized designs. Guo et al. [6, 11] applied SOMbased data mining to optimization results of turbomachinery. Jeong et al. [12] conducted a data mining for aerodynamic design space. However, limited research on data mining in turbomachinery design means that the process of data mining and data mining results are still not very clear.
The purpose of this paper is to develop a multipoint and multiobjective design optimization method for high pressure ratio impellers to achieve better aerodynamic performances at both design and offdesign conditions and with a wider operating range. The remaining of this paper is organized as follows: first, a time series autoregressive (AR) model is developed to predict surge point from the CFD simulation of a single impeller flow passage and validated by experimental results. A selforganization map (SOM) is then carried out on samples of CFD results to explore the relation between objective functions and design parameters. A total of 27 key design variables selected by the SOM are then employed in the followup multiobjective optimization. Finally, the results from the optimization are shown and discussed, and some conclusions and remarks are drawn.
2. Time Series Based Surge Detection
Surge is instability of the centrifugal compressors and is associated with strong unsteadiness of the inlet and outlet pressures and temperatures of the compressor. Therefore, monitoring the time based signal by adopting a Fast Fourier Transform analysis makes it possible to state the instant when the compressor starts to surge and to highlight the typical frequency peak related to surge occurrence [13–15]. However, there is no universal standard for the upper limit of pressure pulsation amplitude in surge detection, making it inconvenient to use. In CFD simulation, though a very accurate CFD model will be the perfect candidate near surge for the rise of numerical instabilities due to the large temporal and spatial gradients related to the actual flow physics. In real design environment, often only a single impeller flow passage is employed in numerical simulation in order to reduce computational time and effort. As a result, the numerical calculation may become unstable and not convergent before or after real surge because of the simplification made to the real compression system in CFD model. This makes surge detection of a new impeller design or judging the surge margin relative to the baseline impeller difficult.
2.1. Autoregressive Model
Here we propose a new method for surge detection in CFD; it is based on autoregressive (AR) statistical pattern recognition algorithms [16] by monitoring model residual variances. As this paper was first written, we found a similar approach which was employed for early surge warning in axial compressors [17].
The slope of pressure rise is a reliable indicator of compressor stability. If the slope is positive the compressor will be unstable. The maximum or peak compressor pressure ratio thus defines the stability limit. Numerical studies of axial compressors [18, 19] showed that this criterion gives good results in predicting compressor stall. Centrifugal compressors may be able to operate into the left of the peak pressure ratio, but compressors will be in minor surge in this condition [20], and compressor outlet and inlet flow conditions will be pulsating and unsteady. Just like in experimental detection of surge, the unsteadiness of numerical simulation results could be utilized to find the peak compressor pressure ratio and thus help to find stability limit of a compressor.
In the surge detection, one may use the frequencies or amplitude of the unsteady signal. These two parameters are however dimensional so their values depend directly on compressor size and speed, making them unsuitable as universal thresholds in surge detection. By contrast, the time series model methods, such as AR and ARMA, are based on monitoring model residual variances, which is independent of specified compressor impellers and speed.
Autoregression is a data processing technique that is commonly used in constructing a model from time sequence data for extracting underline trends. For a stationary time series (), autoregressive process produces the following results:where is model coefficients, the random component of the model, and the variance. If is not zeromean, then the process can be similarly carried out by introducing a new time series where is the mean of . The model coefficients, , and the variance, , can be calculated by the following method: Autocovariance function,
Coefficients are calculated by where The variance is given by By monitoring the variance of compressor unsteady signal, one may judge whether the compressor surges or becomes unstable.
2.2. Surge Detection Method
For the outlet boundary condition, mass flow rate condition is not appropriate as the convergence of CFD cannot be guaranteed near surge condition, resulting in premature sinusoidal waves of flow parameters. Thus, a static pressure condition is adopted and the signal of mass flow rate at LE is monitored. It tends to be a periodic wave but not sinusoidal.
The rise of static pressure at outlet boundary is flattened out near surge condition; a dichotomy or sectional method is applied to the determination of the maximum outlet static pressure near surge condition. This method will be demonstrated in the next section with an example. The flow chart of surge detection is presented in Figure 1. An initial static pressure is specified at impeller outlet before a steady calculation. If this calculation is not convergent, the CFD simulation is switched to unsteady calculation, and the temporal mass flow rate signal at LE is monitored. The AR model is then applied to this time series signal to determine whether the compressor surges at the set pressure level and this is followed by a dichotomy method to find a new static pressure between previous converged steady and present unsteady pressures at the outlet. The detection process continues until a preset minimum pressure rise is reached, and the last steady point is regarded as the last stable point before surge.
2.3. Surge Detection of SRV2AB Impeller
This surge detection approach is applied to a higher pressure centrifugal impeller, SRV2AB impeller [21, 22].
The surge detection validation is carried out by the comparison of surge line between CFD prediction and experimental results of the SRV2AB impeller with vaneless diffuser, which serves as the baseline for the subsequent optimization. The basic design parameters of this compressor are given in Table 1.

A mesh independence investigation was firstly conducted by coarse, medium and fine meshes with total mesh number of 0.3, 0.9, and 1.7 million, respectively (for one blade channel). The mesh has a structured HO topology and the minimum value of is less than 10. The size of the first cell to wall is 5 × 10^{−6} m. As shown in Figure 2, the medium mesh seems to provide a good compromise between accuracy and computational efficiency and hence this mesh is used for the rest of the study (see Figure 3), wherein the gird numbers in , direction is 73, 53, and 201, respectively. The grid number in impeller clearance is 13. The 3D Reynoldsaveraged Navier–Stokes (RANS) equations are applied using the commercial software of NUMECA Fine. The SpalartAllmaras turbulence model, which is highly efficient and suitable for the 3D flow with strong pressure gradient, moderate curvature, and separating flows such as transonic compressor internal flow [11, 23], is adopted to take the turbulence effects into account. The spatial discretization used for the computation is central scheme and the time integration scheme is an explicit fourstep RungeKutta algorithm.
The sketch of dichotomy to reach the highest outlet static pressure at surge point is shown in Figure 4. For the impellers generated in the optimization process, the outlet static pressure near surge were first set to be 511,300 pa, 521,300 pa, 531,300 pa, and 541,300 pa according to Figure 1. A new value of static pressure was then obtained by dichotomy in between 531,300 pa and 541,300 pa and so on. The smallest step of pressure rise was set to be 200 pa. This value is large but was used to speed up the detection. In order to compensate the effect of a relatively large value of the minimum pressure increment, the surge point was actually located by decreasing the last detected stable mass flow rate value by 2.5% of the difference between the last two mass flow rates. The corresponding pressure ratio is then evaluated through a thirdorder polynomial interpolation. In the detection, CFD steady iteration was 2000, the upper border of variance was set to be . The spatial discretization and time integration schemes applied to the unsteady computation are the same as above in this section. CourantFriedrichLevy (CFL) number globally scales the timestep sizes used for the timemarching scheme of the flow solver. A higher value of the CFL number results in a faster convergence but will lead to divergence if the stability limit is exceeded [24]. Typical values of CFL are 1~10 and the maximum CFL number reached in unsteady calculations is set to be 3. The physical timestep size is recommended to be for turbomachinery, where denotes shaft speed (rpm), Nb is blade number [25], and here it is s (equals 0.001 second). The CFL and physical time step seem to provide a good tradeoff between computational stability and efficiency after trial.
The time based mass flow rate signals at LE, monitored at the last stable point and the first unstable point for design speed (50000 rpm), are presented in Figure 5. It shows a periodic but not sinusoidal wave for the first unstable point. The autocovariance and partial autocovariance versus order of AR model based on the time signal is given in Figure 6, where for the last stable point and for the first unstable point, which indicates that the AR model is suitable for the time based signals and their order, are and , respectively. As shown in Figure 7, variance of the time based signal fluctuates periodically for the first unstable point. The distance between two adjacent main peaks shows the surge period (frequency), which is corresponding to the main frequency (largest pulsation amplitude) in the Fast Fourier Transform analysis. The variance of a stable point is about zero, while that for an unstable point fluctuates with large pulsation amplitude, providing a convenient approach to detect instability.
Figure 8 compares predicted and measured surge lines of the SRV2AB impeller with vaneless diffuser at 50,000, 40,000, and 30,000 rpm, corresponding to 100%, 80%, and 60% of design speed. Generally, the prediction matches the experimental results quite well at all operating speeds. This gives confidence in the surge detection method.
3. Data Mining Based on SelfOrganizing Map for Optimization of Design Space of SRV2AB
In the turbomachinery design using CFDbased optimization, it is important to determine a small number of key design variables from design space to simplify design problem. The information about the design space, such as tradeoff between objective functions, the relations between design variables and objective functions, and why performance of the optimized designs has been improved will be useful for this purpose by eliminating the design variables which do not have a large influence on the objective functions, thereby the efficiency as well as the reliability of optimization process may be greatly improved. Furthermore, it is preferable for a designer to provide some alternative or suboptimum solutions for the decision making of the final design. The process to extract information from the design space by detecting the features of the optimization results is called “data mining.” This paper deals with the data mining technique based on SOM.
3.1. Initial Design Variables and Objective Functions
In terms of 3D parameterization of impeller geometry, refer to authors’ previous work [26]. Figure 9 shows the control points of endwalls (, two variables for one control point) and blade camber curves (, one variable for one control point). For the endwalls, 13 control points are applied to the hub and 8 for the shroud. The blade profile is parameterized by the root and the tip sections of both full and splitter blades (4 sections in total), and other camber sections are determined by linear interpolation from root and tip sections. 8 control points are selected for each full blade section while 5 for each splitter blade section.
(a) Endwalls
(b) Root section
(c) Tip section
Table 2 shows variable names and corresponding numbers of control points. Not all the variables of control points are active in the optimization. Table 3 shows active variable numbers, range of variations, and constraints. In total, 45 design variables are selected for the initial design space according to the constraints of mechanical design of SRV2AB.

 
Constraints. Inactive variables are unchanged in the optimization. Diffuser inclination against radial equals 13 deg. Shaft speed is 50000 rpm. 
The impeller will be optimized at 50000 rpm for maximizing isentropic efficiency and total pressure ratio at design operating condition (2.55 kg/s), while striving for smaller surge and larger choke mass flows. The corresponding mathematical expression of aerodynamic optimization of SRV2AB is as follows:
3.2. SelfOrganization Map
SOM expresses the information in a qualitative way and visualizes not only the relation between design variables and objective functions but also the tradeoff between the objective functions. It employs a nonlinear projection algorithm from high to lowdimensions and a clustering technique. This projection is based on selforganization of a lowdimensional array of neurons. In the projection algorithm, the weights between the input vector and the array of neurons are adjusted to represent features of the highdimensional data on the lowdimensional map. Figure 10 shows the schematic map of SOM, where is the number of neurons and is equal to the dimension of input vector. Each neuron is connected to adjacent neurons by a neighborhood relation and usually forms twodimensional hexagonal (see Figure 10(b)) topology.
The learning algorithm of SOM starts with finding the bestmatching unit (winning neuron) which is closest to the input vector , and the th winning neuron is selected as the one having minimal distance value:Once the bestmatching unit is determined, the weight vectors are adjusted not only for the bestmatching unit but also for its neighbors. As shown in Figure 11, based on the distance, the bestmatching unit and its neighboring neurons (situated on the cross of the solid lines, the weight vectors of neurons represent their locations) become closer to the input vector . The adjusted topology is represented with dashed lines. The adjustment of weight vector near the bestmatching unit can be formulated as follows:where is the weight vector of the th neuron at iteration. The amount of adjustment depends on the degree of similarity between a neuron and the input, represented by , and is scaled by the function that plays the role of a “learning rate.” This “learning rate” is called the neighborhood kernel [27]. It is a function of both time (iteration step) and the winning neuron spatial neighborhood . This spatial neighborhood is a timedependent function that defines the set of neurons that are topographically close to the winning neuron. The neurons in the spatial neighborhood adjust their weights according to the same learning rule but with amounts depending on their position with respect to the winner.
Repeating this learning algorithm, the weight vectors become smooth not only locally but also globally. Thus, the sequence of the vectors in the original space results in a sequence of the corresponding neighboring neurons in the twodimensional map.
Once the highdimensional data projected on the twodimensional regular grid, the map can be used for visualization and data mining. The same location on each component map corresponds to the same SOM neuron, which is colored according to its related neutron component values. By comparing the behavior of the color pattern in the same region, one can analyze the correlations among parameters. Parameters are correlated if there exist similar color patterns in the same region of the corresponding component maps.
3.3. SOMBased Data Mining for Choosing Design Space
The SOM was carried out by an inhouse code on initial design space to select key design parameters in the design space for further optimization. 100 impellers generated in random were simulated, in which 49 parameters (4 objectives and 45 active design variables) were analyzed. All the employed parameters are normalized (variance is normalized to one). 25 SOM neutrons were used in component map. Figure 12 presents SOMs colored by 4 objective functions, showing a nonlinear relation; no region exists with good performance for all the objectives. This means that it is difficult to improve all the objectives without any compromise. A tradeoff needs to make among the 4 objectives. For example, the clusters in bottom righthand corner are a better choice for total pressure ratio , surge mass flow , and choke mass flow , but bad in isentropic efficiency . SOMs colored by and show an inverse distribution of colors, so they are in a severe tradeoff relation. One physical example is a high trim impeller which has a large choke flow but tends to have low peak efficiency. Another potential conflict of objectives involves and . SOM colored by on the other hand shows diagonal behavior.
(a)
(b)
(c)
(d)
Figure 13 shows SOMs for some active design variables in the initial design space. Variable 33 (Figure 13(c)) represents and controls the exit blade height of impeller. Larger value of variable 33 means smaller blade tip width and vice versa (refer to Figure 9). The cluster with small values of variable 33 is located in the bottom righthand corner and brings good performance of , , and but deteriorates . The largest value of variable 33 situates on middle left and benefits and , while penalizing and . Figure 13(f), combined with Figure 12, shows the influence of variable 52 or which is the wrap angle of full blade at tip section. Figure 13(e) (variable 44, ) indicates the influence of the wrap angle of full blade at root section. It shows an approximate inverse distribution of color pattern with , indicating a positive correlation between increasing value of variable 44 and decreasing . Variable 58 (Figure 13(g)) represents which governs the back sweep of full blade at tip section (large value of variable 58 represents large sweepback). SOMs colored by variable 58 and show similar color patterns, implying an approximate linear relation between the two. Both smaller and larger values of variable 58 (located in bottom righthand and upper lefthand corner, resp.) bring good performance of ; this is due to smaller sweepback tending to increase the tangential component of absolute velocity to raise total pressure while larger sweepback tries to achieve the same by improving efficiency. Figure 13(h) (variable 68, ) shows the similar effects. Variable 21 (Figure 13(a)), variable 22 (Figure 13(b)), and variable 40 (Figure 13(d)) show severe nonlinear relations with the four objective functions, while variable 22 and variable 40 are more sensitive as the ranges of legends in Figures 13(b) and 13(d) are small. Thus, variable 22 and variable 40 are kept while variable 21 is ignored in the final design space.
(a) Variable 21
(b) Variable 22
(c) Variable 33
(d) Variable 40
(e) Variable 44
(f) Variable 52
(g) Variable 58
(h) Variable 68
Based on data mining results, one may derive the following conclusions. Total pressure ratio, isentropic efficiency, surge mass flow, and choke mass flow are in severe tradeoff relations and compromises are needed in optimization. The performance of high pressure ratio impellers is sensitive to the tip and root sections parameters located near the leading edge of full blade, especially for surge margin. The parameters located near the trailing edge of full and splitter blades also have large effects on the performance. The ranges of the control points of these parameters in optimization need careful planning. Some insensitive variables may be ignored in the design space to simplify the optimization and to improve the efficiency as well as the robustness of optimization.
Table 4 shows the active variables and their ranges in the final design space that were used in the optimization. The number of active variables reduces from 45 to 27. Figure 14 presents active variables of endwalls and blade camber curves in final design space, in which the arrows indicate the moving directions of the active variable in the optimization.
 
Constraints. Same as in Table 3. 
(a) Endwalls
(b) Root section
(c) Tip section
4. Aerodynamic Optimization of SRV2AB
4.1. Optimization Method
A multipoint and multiobjective design optimization method based on an improved NSGAII genetic algorithm [26], the time series surge detection, and data mining technique discussed earlier is applied to the optimization of SRV2AB impeller at design shaft speed (50000 rpm). The four objectives are higher efficiency, pressure ratio, larger choke flow, and smaller surge flow. A Paretobased ranking was adopted for the 4 objectives.
The flow chart of this multipoint and multiobjective design optimization method is shown in Figure 15. The righthand part of the chart presents the core optimization algorithm based on the genetic algorithm. The grid generation and CFD solution were automatically carried out using Autogrid and FineTurbo from NUMECA by templates, respectively. The prediction of total pressure ratio and isentropic efficiency at design operating condition was conducted by imposing design mass flow (2.55 kg/s) boundary condition at vaneless diffuser outlet. Choke mass flow rate was obtained by specifying a low outlet static pressure of 151,300 pa. The determination of surge mass flow is carried out by the surge detection method introduced earlier. The initial design space is presented in Table 3, and the final design space in Table 4. Table 5 shows some essential parameters in the genetic algorithm. Parallel computing technique was applied to the optimization process. It included parallel CFD tool and the multithreaded genetic algorithm in which the parallel CFD code runs. The optimization was carried out on a workstation with 48core Xeon(R) E52670 processor, and the CPU utilization was about 86 percent. The total computational time required to run the optimization was about 1400 hours. After the optimization, the optimum geometry is selected from Paretooptimal front solutions of the last iteration with the smallest surge mass flow.

4.2. Results and Discussions
The overall performances of baseline impeller and optimal impeller are presented in Table 6. It shows an improvement of both total pressure ratio (by 5.3%) and isentropic efficiency (by 1.9%). Meanwhile, surge margin shows a considerable improvement by 6.8% while choke mass flow rises by 1.4% as well. The whole design speed line at 50000 rpm shows improvements for all the 4 objectives.

Figure 16 compares the performances of the optimal impeller with the baselines at three speeds of 50000 rpm, 40000 rpm, and 30000 rpm. At all speeds, compressor surge is improved. Combined with a higher choke flow at the highest speed, this makes the new design particularly suitable for turbocharger application where a wide flow range is required and choke flow at low pressure ratios is relatively unimportant. Compressor peak efficiency and peak pressure ratio both increase at the offdesign speeds.
(a) Total pressure ratio
(b) Isentropic totaltototal efficiency
These results demonstrate the power of this multipoint and multiobjective optimization method. A singlepoint optimization is less likely to get similar results. Note that this multipoint and multiobjective optimization is only conducted at design shaft speed (50000 rpm), and it could further improve offdesign performance if the same technique is applied to both design and offdesign speeds. This will significantly increase computational time and resources (about 3 times of this study). A tradeoff is needed between the effort and the gains.
Figure 17 compares the geometry of the baseline and the optimal impeller. Figure 17(a) shows that the meridional passage (variable 22) is reduced after the inducer throat in the optimal impeller. This explains why the choke flow of the impeller is reduced at the two lower speeds but not at the highest speed. At the low pressure ratios, exducer passage area also plays a rule in impeller choking while inducer throat alone dictates the impeller choke flow at high speeds. Figures 17(b) and 17(c) show the profile or camber changes of the full and splitter blades at root and tip sections. At the leading edge of the full blade (variables 44 and 52), both sections display a larger blade angle and higher blade turning, contributing to a higher efficiency and a smaller inducer area that is unfavourable for surge. The reduction of full blade throat is compensated by a more open splitter (variables 59 and 64). The compressor stability is improved by the reduced impeller tip width (variable 33) and a smaller diffuser gap.
(a) Endwalls
(b) Blade camber curves at root section
(c) Blade camber curves at tip section
The backsweep angles (variables 50 and 63) are slightly reduced in the optimal impeller at its hub for both full blade and splitter blade, and this increases impeller work and pressure ratios. Because the same angle is not reduced at the shroud, the reduction at the hub decreases the diffusion imbalance between the shroud and the hub, resulting in more uniform impeller outflow. This compensates the higher diffuser inlet velocity caused by the reduced hub backsweep angle.
It is concluded that the design variables’ optimization for a better performance corresponds to the analysis in SOM. It proves that SOM analysis is effective in extracting the key design variables and their effects, and the proposed optimization method was able to find a proper tradeoff between all the objectives.
The flow field at design operating condition is analyzed. Figure 18 presents the comparison of entropy distributions at different streamwise sections and streamlines of tip leakage between the baseline and optimal impeller. In section I, a large high entropy region is found at midspan near suction side (circled in red) of the baseline impeller, and the corresponding area in the optimal impeller is considerably smaller. This is due to the weaker and detached LE shockwave moving back toward LE; see Figure 19. While the higher blade turning of the new impeller at LE area strengthens the tip leakage flow and thus higher loss is found in the tip region at Section I, the loss stops at Section III.
(a) Baseline
(b) Optimal
(a) Baseline 60% span
(b) Optimal 60% span
(c) Baseline 95% span
(d) Optimal 95% span
From Figure 19 it can be seen that impeller exit flow is now more uniform circumferentially between the two flow channels.
Figure 20 presents the spanwise performance of the impeller 10 mm downstream from the impeller exit. It shows that both the total pressure ratio and the isentropic efficiency have increased in the regions ranging from the hub to 80% span but reduced near the shroud.
(a) Total pressure ratio
(b) Isentropic efficiency
5. Conclusions and Remarks
A multipoint and multiobjective design optimization strategy of centrifugal impeller is proposed by integrating a genetic algorithm, 3D geometry parameterization, CFD tools, time series based surge detection method, and SOMs based data mining technique. This approach was successfully applied to high pressure ratio centrifugal impeller SRV2AB for higher total pressure ratio and better efficiency at design operating condition and smaller surge mass flow and larger choke mass flow at design speed. The main conclusions are drawn as follows:
(1) A time series based surge detection method was introduced. It uses autoregression to detect compressor instability and is independent of the impeller geometry, shaft speed, and boundary condition. This method was successfully applied to SRV2AB impeller surge prediction.
(2) By SOMbased data mining on initial design space, tradeoff relations between objective functions and correlations among design variables and objective functions were visualized and analyzed. The key design variables were then identified and kept in the final design space.
(3) The optimization improves the overall performance of the impeller at whole design speed line and widens the compressor flow range. At offdesign speeds and compressor efficiency, surge margin is also enhanced. The mechanism behind the performance improvement is further explained by combining geometry changes with detailed flow analysis.
The surge detection method still needs to be checked by experiment and by more applications. It is currently computational intensive. It essentially detects surge by the macroscopic time based signal instead of local flow features such as the development of stall cells, which may be more efficient. The method needs further refinement.
Another future work concerns with the method of selecting the data for data mining. Results of data mining techniques may depend on the used data. For the consistency of information obtained from data mining, a robust data selection method is necessary.
Nomenclature
AR:  Autoregressive processes 
ARMA:  Autoregressive moving average processes 
:  Amplitude of time based signal 
:  White noise series 
deg:  Degree 
GA:  Genetic algorithm 
LE:  Leading edge 
MA:  Moving average processes 
:  Mass flow rate 
:  Choke mass flow 
:  Weight vector 
:  Surge mass flow (last stable mass flow near surge condition) 
NID:  Normally and independently distributed 
:  Shaft speed (rpm) 
:  Probability 
:  Upper border of crossover probability 
:  Lower border of crossover probability 
:  Upper border of mutation probability 
:  Lower border of mutation probability 
:  Pressure 
:  Order of AR model 
:  Radial coordination (mm) 
:  Autocovariance function 
SOM:  Selforganization map 
TE:  Trailing edge 
:  Original time series signal 
:  Zeromean time series 
:  Axial coordination (mm). 
:  Upper border of variance in AR model 
:  Total pressure ratio 
:  Isentropic efficiency 
:  Azimuthal angle (rad) 
:  Mean value 
:  Autocorrelation function 
:  Variance 
:  Coefficients of AR model 
:  Partial autocovariance function. 
:  Crossover 
choke:  Choke condition 
design:  Design condition 
full:  Full blade 
is:  Isentropic 
:  Mutation 
root:  Root section 
splitter:  Splitter blade 
surge:  Surge condition 
tot:  Total conditions 
tip:  Tip section. 
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This research work is supported by National Natural Science Foundation of China (Grant no. 11672206) and State Key Laboratory of Internal Combustion Engine Burning (Grant no. K201607).
References
 P. G. Tucker, “Computation of unsteady turbomachinery flows: part 1—progress and challenges,” Progress in Aerospace Sciences, vol. 47, no. 7, pp. 522–545, 2011. View at: Publisher Site  Google Scholar
 L. He and J. Yi, “Twoscale methodology for URANS/large eddy simulation solutions of unsteady turbomachinery flows,” Journal of Turbomachinery, vol. 139, no. 10, p. 101012, 2017. View at: Publisher Site  Google Scholar
 D. Pasquale, G. Persico, and S. Rebay, “Optimization of turbomachinery flow surfaces applying a CFDbased throughflow method,” Journal of Turbomachinery, vol. 136, no. 3, p. 031013, 2014. View at: Publisher Site  Google Scholar
 R. Hunziker, H. Dickmann, and R. Emmrich, “Numerical and experimental investigation of a centrifugal compressor with an inducer casing bleed system,” Journal of Power and Energy, vol. 215, no. 6, pp. 783–791, 2005. View at: Publisher Site  Google Scholar
 H. Krain, “Review of centrifugal compressor's application and development,” Journal of Turbomachinery, vol. 127, no. 1, pp. 25–34, 2005. View at: Google Scholar
 Z. Guo, Z. Zhou, L. Song, J. Li, and Z. Feng, “Aerodynamic Analysis and MultiObjective Optimization Design of a High Pressure Ratio Centrifugal Impeller,” in Proceedings of the ASME Turbo Expo 2014: Turbine Technical Conference and Exposition, p. V02DT42A013, Düsseldorf, Germany, 2014. View at: Publisher Site  Google Scholar
 T. Verstraete, Z. Alsalihi, and R. A. Van den Braembussche, “Multidisciplinary optimization of a radial compressor for microgas turbine applications,” Journal of Turbomachinery, vol. 132, no. 3, pp. 1291–1299, 2007. View at: Publisher Site  Google Scholar
 K. HyunSu, S. YooJune, and K. YounJea, “Optimal design of impeller for centrifugal compressor under the influence of fluidstructure interaction, ASME/JSME/KSME joint fluids engineering conference,” in Proceedings of the Optimal Design of Impeller for Centrifugal Compressor Under the Influence of FluidStructure Interaction, ASME/JSME/KSME Joint Fluids Engineering Conference, p. 02248, 2015. View at: Google Scholar
 A. Demeulenaere, J. Bonaccorsi, D. Gutzwiller, L. Hu, and H. Sun, “Multidisciplinary multipoint optimization of a turbocharger compressor wheel,” in Proceedings of the ASME Turbo Expo 2015: Turbine Technical Conference and Exposition, p. V02CT45A020, Montreal, Quebec, Canada, 2015. View at: Publisher Site  Google Scholar
 M. Pini, G. Persico, and V. Dossena, “Robust adjointbased shape optimization of supersonic turbomachinery cascades,” in Proceedings of the ASME Turbo Expo 2014: Turbine Technical Conference and Exposition, p. V02BT39A043, Düsseldorf, Germany, 2014. View at: Publisher Site  Google Scholar
 Z. Guo, L. Song, J. Li, G. Li, and Z. Feng, “Research on metamodel based global design optimization and data mining methods,” in Proceedings of the ASME Turbo Expo 2015: Turbine Technical Conference and Exposition, p. V02CT45A007, Montreal, Quebec, Canada, 2015. View at: Publisher Site  Google Scholar
 S. Jeong, K. Chiba, and S. Obayashi, “Data mining for aerodynamic design space,” Journal of Aerospace Computing Information & Communication, vol. 2, no. 11, pp. 452–469, 2005. View at: Google Scholar
 N. González Díez, J. P. Smeulers, L. Tapinassi, A. S. Del Greco, and L. Toni, “Predictability of rotating stall and surge in a centrifugal compressor stage with dynamic simulations,” in Proceedings of the ASME Turbo Expo 2014: Turbine Technical Conference and Exposition, p. V02DT44A032, Düsseldorf, Germany, 2014. View at: Publisher Site  Google Scholar
 J. Galindo, J. Serrano R, C. Guardiola, and C. Cervelló, “Surge limit definition in a specific test bench for the characterization of automotive turbochargers,” Experimental Thermal & Fluid Science, vol. 30, no. 5, pp. 449–462, 2006. View at: Publisher Site  Google Scholar
 S. Marelli, C. Carraro, G. Marmorato, G. Zamboni, and M. Capobianco, “Experimental analysis on the performance of a turbocharger compressor in the unstable operating region and close to the surge limit,” Experimental Thermal & Fluid Science, vol. 53, no. 2, pp. 154–160, 2014. View at: Publisher Site  Google Scholar
 R. Yao and N. Pakzad S, “Autoregressive statistical pattern recognition algorithms for damage detection in civil structures,” Mechanical Systems & Signal Processing, vol. 31, no. 8, pp. 355–368, 2012. View at: Publisher Site  Google Scholar
 L. Fanyu, L. Jun, D. Xu, S. Dakun, and S. Xiaofeng, “Stall warning approach with application to stall precursorsuppressed casing treatment,” in Proceedings of the ASME Turbo Expo 2016: Turbomachinery Technical Conference and Exposition, p. V02DT44A038, Seoul, South Korea, 2016. View at: Publisher Site  Google Scholar
 C. Hah, D. C. Rabe, and A. R. Wadia, “Role of tipleakage vortices and passage shock in stall inception in a swept transonic compressor rotor,” in Proceedings of the ASME Turbo Expo 2004: Power for Land, Sea, and Air, pp. 545–555, Vienna, Austria, 2004. View at: Publisher Site  Google Scholar
 C. Hah, J. Bergner, and H. Schiffer, “Short lengthscale rotating stall inception in a transonic axial compressor: criteria and mechanisms,” in Proceedings of the ASME Turbo Expo 2006: Power for Land, Sea, and Air, pp. 61–70, Barcelona, Spain, 2006. View at: Publisher Site  Google Scholar
 X. Zheng, Z. Sun, T. Kawakubo, and H. Tamaki, “Experimental investigation of surge and stall in a turbocharger centrifugal compressor with a vaned diffuser,” Experimental Thermal and Fluid Science, vol. 82, pp. 493–506, 2017. View at: Publisher Site  Google Scholar
 G. Eisenlohr, H. Krain, F. Richter, and V. Tiede, “Investigations of the flow through a high pressure ratio centrifugal impeller,” in Proceedings of the ASME Turbo Expo 2002: Power for Land, Sea, and Air, pp. 649–657, Amsterdam, The Netherlands, 2002. View at: Publisher Site  Google Scholar
 M. Zangeneh, N. Amarel, K. Daneshkhah, and H. Krain, “Optimization of 6.2:1 pressure ratio centrifugal compressor impeller by 3D inverse design,” in Proceedings of the ASME 2011 Turbo Expo: Turbine Technical Conference and Exposition, pp. 2167–2177, Vancouver, British Columbia, Canada, 2011. View at: Publisher Site  Google Scholar
 F. Ning and L. Xu, “Numerical investigation of transonic compressor rotor flow using an implicit 3d flow solver with oneequation spalartallmaras turbulence model,” in Proceedings of the ASME Turbo Expo 2001: Power for Land, Sea, and Air, p. V001T03A054, New Orleans, LA, USA, 2001. View at: Publisher Site  Google Scholar
 D. X. Wang and X. Huang, “Solution stabilization and convergence acceleration for the harmonic balance equation system,” Journal of Engineering for Gas Turbines & Power, 2017. View at: Google Scholar  MathSciNet
 Ansys Fluent 16.0 User’s Guide, “Ansys Inc,” 2015. View at: Google Scholar
 X. Li, Y. Zhao, Z. Liu, and H. Chen, “The optimization of a centrifugal impeller based on a new multiobjective evolutionary strategy,” in Proceedings of the ASME Turbo Expo 2016: Turbomachinery Technical Conference and Exposition, p. V02CT39A022, Seoul, South Korea, 2016. View at: Publisher Site  Google Scholar
 K. J. Cios, W. Pedrycz, and R. W. Swiniarski, Data Mining Methods for Knowledge Discovery, Kluwer Academic Publisher, 1998. View at: Publisher Site
Copyright
Copyright © 2017 Xiaojian Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.