Hindawi Publishing Corporation Advances in Mechanical Engineering Volume 2014, Article ID 457542, 10 pages http://dx.doi.org/10.1155/2014/457542
Research Article Design Optimization of Centrifugal Pump Using Radial Basis Function Metamodels Yu Zhang, Jinglai Wu, Yunqing Zhang, and Liping Chen School of Mechanical Science & Engineering, Center for Computer-Aided Design, Huazhong University of Science & Technology, Wuhan, Hubei 430074, China Correspondence should be addressed to Yunqing Zhang;
[email protected] Received 2 January 2014; Revised 13 March 2014; Accepted 13 March 2014; Published 19 May 2014 Academic Editor: Marco Ceccarelli Copyright © 2014 Yu Zhang 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. Optimization design of centrifugal pump is a typical multiobjective optimization (MOO) problem. This paper presents an MOO design of centrifugal pump with five decision variables and three objective functions, and a set of centrifugal pumps with various impeller shroud shapes are studied by CFD numerical simulations. The important performance indexes for centrifugal pump such as head, efficiency, and required net positive suction head (NPSHr) are investigated, and the results indicate that the geometry shape of impeller shroud has strong effect on the pump’s performance indexes. Based on these, radial basis function (RBF) metamodels are constructed to approximate the functional relationship between the shape parameters of impeller shroud and the performance indexes of pump. To achieve the objectives of maximizing head and efficiency and minimizing NPSHr simultaneously, multiobjective evolutionary algorithm based on decomposition (MOEA/D) is applied to solve the triobjective optimization problem, and a final design point is selected from the Pareto solution set by means of robust design. Compared with the values of prototype test and CFD simulation, the solution of the final design point exhibits a good consistency.
1. Introduction Centrifugal pump, the most frequently used type of pump up to now, is an important member of turbomachine family for transporting liquids from a low level to a high level. As a typical kind of centrifugal pump, double suction centrifugal pumps with large flux and high lift are widely used in industry. However, traditional methods for pump design and optimization mainly adopt the trial-and-error strategy by manufacturing and testing prototype pumps, which are expensive and time consuming [1]. Since the huge market demand for centrifugal pumps is around 20 billion US dollar per year [2], even one percent increment of pump efficiency is of important significance. As the “heart” of centrifugal pump, the geometry shape of pump impeller has strong effect on the pump performance, so that the optimization of impeller shape should be paid much attention in order to enhance pump efficiency. Pierret and van den Braembussche [3] optimized the blade shapes by means of artificial neural networks and
simulated annealing to improve pump performance. Anagnostopoulos [4] studied the characteristic performance curves of centrifugal pump using CFD software and found the optimal geometry shape of impeller blade that could maximize the pump efficiency among a set of blade angles. Using incomplete sensitivity method and genetic algorithms, Derakhshan et al. [5] redesigned the shape of impeller blades to obtain a higher efficiency. Besides, Zhou et al. [6] investigated the best shape parameter combination of impeller to obtain higher pump efficiency by means of orthogonal experiment and CFD simulation. Obviously, these researches only viewed the design and optimization of centrifugal pump as a single objective optimization for pump efficiency but ignored other performance indexes such as head and NPSHr. Although higher efficiency seemed to be the most important performance for pumps which could produce considerable economic benefit, other performance indexes were also significant. Therefore, the optimization of double suction centrifugal pumps was actually a multiobjective optimization (MOO) problem with multiple design requirements.
Downloaded from ade.sagepub.com by guest on August 26, 2015
2
Advances in Mechanical Engineering
MOO, also called as vector optimization, is a numerical process of optimizing a collection of individual objective functions simultaneously and systematically [7]. Many practical engineering problems usually involve multiple design requirements conflicting with each other, and the optimization of those problems essentially belongs to MOO problems [8]. There are some works related to MOO and design for centrifugal pumps. For example, Oyama and Liou [9] used evolutionary algorithm to find the Pareto front of total head and input power of centrifugal pump, which were two conflictive objectives; Safikhani et al. [10] presented an MOO progress on centrifugal pumps and proposed the Pareto optimal solutions using genetic algorithm based on neural network metamodels; Papierski and Błaszczyk [11] decomposed the MOO and design problem of centrifugal pump to maximize the efficiency and simultaneously minimize NPSHr. These researches had good contributions to the MOO of centrifugal pump, which could provide meaningful references for the design of centrifugal pump. However, the optimization and design of centrifugal pump are an integrated process, and more than two performance indexes should be taken into account simultaneously. In practical optimization and design of centrifugal pump, there are several important performance indexes that are usually considered as the optimization objectives, such as head, efficiency, and NPSHr. Among these objectives, head and efficiency are required to be maximized, while NPSHr needs to be minimized. As a matter of fact, some of the optimization objectives conflict with each other and it is practically impossible to optimize the three objectives at the same time. Moreover, these three performance indexes are often obtained by means of experiment tests or CFD simulations and cannot be used for the optimization of pump performance iteratively, unless appropriate metamodels are established between the decision variables and the concerned objective functions. Radial basis function (RBF) was originally developed by Hardy to fit the irregular topographic contours of geographical data [12]. Compared with other metamodels such as Kriging and neural networks, RBF is relatively simple to implement [13]. Jin et al. [14] deeply compared four popular types of metamodels and concluded that RBF was the most reliable method in most situations in terms of accuracy and robustness. As aforementioned, head, efficiency, and NPSHr are usually the most concerned objectives in pump design. Therefore, an effective and powerful MOO algorithm is needed to obtain the Pareto optimal solutions. Since 1990s, many researchers have proposed a large amount of multiobjective evolutionary algorithms [15–18], which are considered as the most successful approaches in identifying the Pareto optimal set of MOO problems [19]. Thereinto, MOEA/D proposed by Zhang and Li is a well-known multiobjective evolutionary algorithm [20], which is characterized by the scalarizing functions with uniformly distributed weight vectors. Due to the good performance of global searching, MOEA/D is usually a good choice for dealing with multiobjective problems involving multiple conflictive objectives.
The objective of this work is to improve head, efficiency, and cavitation performance of centrifugal pump simultaneously, and MOEA/D is employed to solve the triobjective optimization problem of centrifugal pump. Although previous researches have optimized some indexes, they do not consider all these performance indexes at the same time. Moreover, most previous works about design optimization of centrifugal pump can be summarized into two categories: the first category of optimization methods couples CFD simulation with experiment, which is of high cost and time consuming; the other one combines metamodels with CFD simulation, which cannot guarantee the accuracy of the results without the verification of experiment. In this paper, the proposed method reasonably takes the advantages of CFD simulation, RBF metamodels, and experiment, and thus the above faults in the practical optimization of centrifugal pump can be well overcome.
2. Basic Concepts 2.1. Multiobjective Optimization. A general multiobjective optimization problem mainly consists of decision variables, objective functions, and constraint conditions. Multiobjective optimization is to find the decision variables x that minimize (or maximize) the vector of objective functions 𝐹(x) to a number of constraints or bounds. Multiobjective optimization problem can be formulated as 𝑇
Find x = [𝑥1 , 𝑥2 , . . . , 𝑟𝑛 ] 𝑇 Minimize 𝐹 (x) = min [𝑓1 (x) , 𝑓2 (x) , . . . , 𝑓𝑚 (x)] Subject to 𝑔𝑗 (x) ≤ 0 (𝑗 = 1, 2, . . . , 𝑝) ℎ𝑘 (x) ≤ 0 (𝑘 = 1, 2, . . . , 𝑞) 𝑥𝑖𝐿 ≤ 𝑥𝑖 ≤ 𝑥𝑖𝑈 (𝑖 = 1, 2, . . . , 𝑛) ,
(1)
where the vector function 𝐹(x) consists of 𝑓1 (x), 𝑓2 (x), . . . , 𝑓𝑚 (x) to be optimized at the same time and 𝑥𝑖𝐿 and 𝑥𝑖𝑈 are the lower boundary and upper boundary of 𝑥𝑖 , respectively. In practical engineering MOO problem, the objective functions and constraint conditions are often expensive computational models, such as computational fluid dynamics (CFD), multibody dynamics (MBD), and finite element analysis (FEA). 2.2. Pareto Optimal Set. The concept of Pareto optimal set is generally used to characterize the compromising solutions to an MOO problem, and the Pareto optimality is defined as follows [21]. A vector of x is Pareto optimal solution if and only if there exists no feasible vector x that would improve some objective functions without causing a simultaneous worsening of at least one other objective function. All Pareto optimal solutions constitute the Pareto optimal set.
3. Description and CFD Numerical Simulation of Centrifugal Pump 3.1. Decision Variables. The working process of centrifugal pump is actually accompanied with energy conversion and
Downloaded from ade.sagepub.com by guest on August 26, 2015
Advances in Mechanical Engineering
3
(a) Pump impeller
(b) Pump volute casing
Figure 1: The double suction centrifugal pump. P4
60
A
l P3
r2
333
𝛼2 r1 𝛼1
P2
B
59
155
P0
P1
(a) Main dimensions
(b) Meridional section
Figure 2: The main dimensions and meridional section of impeller (unit: mm).
loss. Impeller is the core component of a centrifugal pump, which converts the mechanical rotation to the kinetic energy of the fluid. As aforementioned, the geometry shape of impeller shroud has strong effect on pump performance. Therefore, the impeller design is significant for obtaining high performance centrifugal pump. Notably, too many blades may cause low blade loading due to higher friction losses; while fewer blades may result in a higher blade loading, and the turbulent dissipation losses will rise because of increased secondary flow and stronger deviation between blade and flow direction. For most centrifugal pumps, the number of impeller blade is usually recommended between 5 and 7 while the range of specific speed is 10∼120 [2]. Thus,
six blades are adopted here. Figure 1 shows the geometric models of the double suction centrifugal pump, including (a) pump impeller and (b) pump volute casing. Due to the symmetry character, the design of double suction impeller can be translated into the design of a single suction type. The impeller’s main dimensions in this paper are shown in Figure 2(a), and the meridional section of the impeller is also given in Figure 2(b), in which the solid line is controlled by quartic B´ezier curve with five control points. It is clear that the meridional section of pump impeller can be parameterized by five decision variables of 𝛼1 , 𝑟1 , 𝛼2 , 𝑟2 , and 𝑙. Here, 𝑟 is the relative position in the line segment. Taking line segment 𝑃0 𝑃2 , for example, 𝑟1 is the ratio of the line length 𝑃0 𝑃1 to 𝑃0 𝑃2 .
Downloaded from ade.sagepub.com by guest on August 26, 2015
4
Advances in Mechanical Engineering Table 1: Decision variables and their boundaries.
Design variable 𝛼1 (deg) 𝑟1 𝛼2 (deg) 𝑟2 𝑙 (mm)
Lower boundary 0 0.02 70 0.02 145
Table 2: Basic parameters for numerical simulations.
Upper boundary 30 0.98 90 0.98 195
𝑙 is the length of line segment 𝐴𝐵. The boundaries of the five decision variables are listed in Table 1. 3.2. Objective Functions. Head, efficiency, and NPSHr are significant performance indexes of centrifugal pump, which are usually the most concerned objectives in pump design. Therefore, they are chosen as the objective functions in this paper. 3.2.1. Head. Head is also termed as total dynamic head. The monometric head represents the total useful mechanical energy transferred by the pump to the fluid per unit mass. In practice, head can be expressed as the energy difference between pump inlet and outlet in unit of length, which can be formulated as the following equation: 𝐻 = 𝐸2 − 𝐸1 =
𝑢2 − 𝑢12 𝑃2 − 𝑃1 + 𝑧2 − 𝑧1 + 2 , 𝜌𝑔 2𝑔
(2)
where 𝐸1 and 𝐸2 denote total energy of per unit mass of fluid at pump inlet and outlet, respectively; 𝑃1 and 𝑃2 are pump inlet pressure and outlet pressure, respectively; 𝑧1 and 𝑧2 are used to denote pump inlet and outlet position head, respectively; 𝑢1 and 𝑢2 denote the fluid velocity at pump inlet and outlet, respectively. 3.2.2. Efficiency. The efficiency of centrifugal pump represents the ratio of useful power to shaft power, which can be defined as follows: 𝜂=
𝑃𝑢 𝜌𝑔𝐻𝑄 = , 𝑃𝑠 𝑃𝑠
(3)
where 𝑃𝑢 denotes the useful power; 𝑃𝑠 is the shaft power which is also called input power. Owing to various losses in the working process such as disc friction loss, mechanical loss, and frictional loss, the kinetic energy cannot be completely converted to pressure energy, so that 𝑃𝑠 is larger than 𝑃𝑢 . 3.2.3. NPSHr. Cavitation is a very dangerous phenomenon that must be avoided for centrifugal pump, and the damages caused by cavitation may range from minor pitting to catastrophic failure. Impeller is usually the component where most damages occur, since cavitation may erode the pump impeller after some period of time. NPSHr is a significant performance index of pump cavitation, which relates to the minimum pressure required at the suction port of the pump to prevent the pump from
Parameter Flow rate 𝑄 Rotational speed Number of blades Inlet operating pressure Flowing medium
Value 2000 (m3 /h) 1400 r/min 6 1 (atm) Water (20∘ C)
cavitation. NPSHr can be formulated with the following equation: NPSHr =
𝑃1 − 𝑃min 𝑢12 + , 𝜌𝑔 2𝑔
(4)
where 𝑃min is the minimum pressure of the impeller blade. The pump with lower value of NPSHr performs better in cavitation resistance. 3.3. CFD Numerical Simulation. CFD numerical simulation has been proved as a very effective method to investigate the flow inside centrifugal pumps at the stages of both pump design and performance analysis. Moreover, CFD numerical simulation can reduce the construction cost in the design and optimization of pump components [22]. Latin hypercube design (LHD) is a space-filling design technique, which can ensure that each decision variable has all portions of its represented range. Using the LHD method, 119 sample points and 30 test points can be obtained. The sample points are used as the input data of RBF metamodels, while the test points are used to validate the accuracy of the RBF metamodels. Based on the sample points and test points, CFD models of pumps can be built. Figure 3 shows one CFD model of the double suction centrifugal pumps, including a detailed view of the pump impeller. In this paper, the basic parameters used in CFD numerical simulations are shown in Table 2. The flow through the modeled centrifugal pump is simulated using the commercial code Fluent, and the flow models are established with the standard k-𝜀 turbulence models and logarithmic law functions, which are consistent with the no-slip condition. Static boundary condition is imposed on the boundary of volute flow domain, while rotary boundary condition is imposed on the boundary of impeller flow domain. Furthermore, the interaction between these two boundaries is taken into account through the multiple reference frame (MRF) model. The RANS equations are solved by finite volume method (FVM), and the pressure-velocity coupling is calculated by means of the SIMPLEC algorithm. Second order upwind discretizations are used for convective and diffusive terms of the turbulence model equations. The residual error is set as 1 × 10−5 to judge whether the calculation is convergent. Figure 4 shows the pressure and velocity distribution from one case of the CFD simulations, which provides an intuitive reflection on the characteristics of the pump’s internal flow field. As shown in the figure, after the fluid flows into the pump impeller, the pressure and velocity of each flow passage mainly present a uniform distribution along
Downloaded from ade.sagepub.com by guest on August 26, 2015
Advances in Mechanical Engineering
5
(a)
(b)
50 45 40 35 30 25 20 15 10 5
Pressure
1.3E + 06 1.2E + 06 1.1E + 06 1E + 06 900000 800000 700000 600000 500000 400000 300000 200000 100000 0 −100000 −200000
(a) Pressure distribution (unit: Pa)
Velocity magnitude
Figure 3: Three-dimensional mesh model of double suction centrifugal pump.
(b) Velocity distribution (unit: m/s)
Figure 4: One case of CFD simulation results.
the flow direction. The pressure first drops to the minimum from the impeller inlet to outlet and then rises gradually to the maximum at the outlet. During this process, the lowest pressure appears at the inlet, closing to the blade inlet side, where cavitation is easy to occur. As for the velocity, it rises gradually from the impeller inlet to outlet, and no velocity decline occurs all the way. Because of the volute, the pressure and velocity distributions are not symmetric, whereas both change periodically with the continuous rotation of the impeller.
4. RBF Modeling and Verification RBF model uses linear combinations of some radial basis functions based on Euclidean distance, and it is mathematically defined as follows: 𝑁
𝐹 (x) = ∑𝜆 𝑖 Φ (x − x𝑖 ) . 𝑖=1
(5)
Equation (5) can also be rewritten in the form of matrix: A𝜆 = f,
(6)
where f = [𝐹(x1 ), 𝐹(x2 ), . . . , 𝐹(x𝑁)]𝑇 , 𝜆 = [𝜆1 , 𝜆2 , . . . , 𝜆𝑁]𝑇 , and Φ (x1 − x1 ) ⋅ ⋅ ⋅ Φ (x1 − x𝑁) ] [ .. .. A=[ ]. . d . [Φ (x𝑁 − x1 ) ⋅ ⋅ ⋅ Φ (x𝑁 − x𝑁)]
(7)
There are several classes of radial basis functions that can be chosen for Φ, such as cubic, thin plate spline, Gaussian, multiquadric, and inverse multiquadric functions. However, the most popular radial function is Gaussian radial basis function in the form of (8). The value of user-defined constant 𝑐 is usually recommended between 0.002 and 0.03 [23], and the most commonly used is 0.01 [24]. In this paper,
Downloaded from ade.sagepub.com by guest on August 26, 2015
6
Advances in Mechanical Engineering Table 3: The sample points and corresponding results of CFD simulations.
No. 1 2 3 4 5 6 7 8 9 10 11 12 .. . 117 118 119
𝛼1 (deg) 0.00 0.25 0.51 0.76 1.02 1.27 1.53 1.78 2.03 2.29 2.54 2.80 .. . 29.49 29.75 30.00
Decision variable 𝛼2 (deg) 78.81 79.49 76.44 70.68 87.97 76.27 75.93 82.20 80.68 84.75 86.44 72.37 .. . 77.63 87.80 82.54
𝑟1 0.6546 0.7359 0.1746 0.8824 0.1664 0.6058 0.7603 0.1176 0.2397 0.4186 0.6220 0.2559 .. . 0.5000 0.9800 0.6953
𝑟2 0.3048 0.5163 0.9556 0.3861 0.0363 0.7685 0.4512 0.8824 0.7848 0.9068 0.2641 0.8254 .. . 0.7766 0.0281 0.5976
the Gaussian basis function with 𝑐 = 0.01 is used to develop and test the metamodels for centrifugal pump: −𝑐𝑟2
Φ (𝑟) = 𝑒
.
(8)
The sample points obtained by LHD and corresponding results of CFD simulation are shown in Table 3. According to the CFD numerical simulation results, pump head, efficiency 𝜂, and NPSHr can be calculated as the output data of the metamodels using (2), (3), and (4), respectively. Based on the input-output data in Table 3, the RBF metamodels can finally be established after calculating the interpolation coefficient vector 𝜆 in (6). To be specific, the geometric parameters of pump impeller, 𝛼1 , 𝑟1 , 𝛼2 , 𝑟2 , and 𝑙, are the input variables, while the corresponding sets of head, efficiency 𝜂, and NPSHr are the output variables. The metamodels can be applied to multiobjective optimization only if the prediction accuracy of the RBF models is higher enough; otherwise, the metamodels should be rebuilt. In order to validate the accuracy of the RBF predictor, an additional set of 30 points obtained by LHD are used as test points. The root mean square error (RMSE) between the estimated values obtained by the RBF metamodels and the simulated values from the CFD simulation is defined as follows: 2
∑𝑀 (𝑌𝑖 (RBF) − 𝑌𝑖 (CFD) ) . RMSE = √ 𝑖 𝑁
(9)
RSME is generally used to measure the average error between values predicted by a model and the actual values, which is an indicator of whether the metamodels satisfy the requirement or not. The RSME values of head, efficiency, and NPSHr calculated by (9) are shown in Table 4, and the results indicate that the prediction accuracy of the RBF metamodels is pretty reliable.
𝑙 (mm) 171.27 149.66 151.78 157.29 159.83 170.42 155.17 178.05 170.00 165.34 164.49 179.75 .. . 154.75 183.14 164.92
Head (m) 131.06 130.47 136.11 140.94 120.02 135.21 134.75 128.93 130.43 125.28 121.60 140.82 .. . 134.29 119.36 127.50
Objective parameter 𝜂 (%) NPSHr (m) 89.02 4.52 90.11 4.53 87.12 3.85 86.64 4.17 90.03 4.82 88.85 4.13 88.13 4.34 89.44 4.31 89.67 4.30 91.42 4.55 92.32 4.98 84.37 3.71 .. .. . . 79.35 3.72 85.46 4.67 81.77 4.13
Table 4: The RSME values of head, efficiency, and NPSHr. Head 1.0786
RSME
Efficiency 0.9441
NPSHr 0.0422
Figures 5, 6, and 7 show the comparison between the RBF predication values and the CFD simulation results of head, efficiency, and NPSHr, respectively. It can be found that the estimated values of RBF metamodels coincide well with the CFD simulation values at the test points. Therefore, the triobjective optimization design of centrifugal pump can be executed based on the RBF metamodels.
5. Multiobjective Optimization and Design of Centrifugal Pump 5.1. Problem Description. As aforementioned, the optimization design of double suction centrifugal pump is an MOO problem with more than two objectives, rather than a single objective optimization problem. This MOO problem can be described as follows: Find Maximize Maximize Minimize Subject to
𝑇
x = [𝛼1 , 𝑟1 , 𝛼2 , 𝑟2 , 𝑙] Head = 𝑓1 (𝛼1 , 𝑟1 , 𝛼2 , 𝑟2 , 𝑙) 𝜂 = 𝑓2 (𝛼1 , 𝑟1 , 𝛼2 , 𝑟2 , 𝑙) NPSHr = 𝑓3 (𝛼1 , 𝑟1 , 𝛼2 , 𝑟2 , 𝑙) 0∘ ≤ 𝛼1 ≤ 30∘ 0.02 ≤ 𝑟1 ≤ 0.98 70∘ ≤ 𝛼2 ≤ 90∘ 0.02 ≤ 𝑟2 ≤ 0.98 145 mm ≤ 𝑙 ≤ 195 mm,
(10)
where 𝑓1 , 𝑓2 , and 𝑓3 are the RBF approximation of head, efficiency 𝜂, and NPSHr, respectively.
Downloaded from ade.sagepub.com by guest on August 26, 2015
Advances in Mechanical Engineering
7
145
5 4.5
135
NPSHr (m)
Head (m)
140
130 125
3.5
120 115
4
0
5
10
15 20 Serial number
25
30
RBF CFD
3
0
5
10
15 20 Serial number
25
30
RBF CFD
Figure 5: The results of pump head at test points.
Figure 7: The results of NPSHr at test points.
95
Efficiency (%)
90 85 80 75 70
0
5
10
15 20 Serial number
25
30
RBF CFD
Figure 6: The results of efficiency at test points.
MOEA/D algorithm is employed to solve the MOO problem. Parameters of the MOEA/D algorithm are set as follows: SBX and polynomial mutation are used as operators for crossover and mutation, respectively, and the distribution indexes for both operators are 𝜅𝑐 = 20 and 𝜅𝑚 = 20, respectively; the crossover probability is set as 0.9 and the mutation probability as 1/𝑛, where 𝑛 is the number of decision variables; the weight vectors consist of 100 uniformly distributed vectors; 𝑇 is set as 25, which is the number of the weight vectors in the neighborhood of each weight vector. The MOEA/D algorithm runs with the population size of 150 and the maximum generation number of 300. In order to provide results with statistical confidence, the MOEA/D algorithm runs 30 times independently for the problem.
objective, the others will go worse, and vice versa. These nondominated optimal points obviously present tradeoffs between the objective function efficiency and head, as well as NPSHr and efficiency. For a centrifugal pump designer, each point in the Pareto optimal set is a feasible choice. As a matter of fact, it is impossible and not realistic to select all of the Pareto optimal solutions as the design schemes, especially for the practical engineering problem. Instead, a tradeoff point should be chosen among the Pareto optimal solutions as the final pump design point. However, the variation of design parameters is inevitable in the practical manufacturing process of centrifugal pumps, which can influence the characteristic parameters of centrifugal pumps. Thus, the final pump design point obtained by traditional methods cannot smoothly settle this problem. In order to overcome the defect of the traditional methods and to guide the tradeoff selection, a robust design method based on Monte Carlo simulation is developed to determine the best point. By evaluating a set of random design alternatives, Monte Carlo simulation can obtain the stochastic characteristics of an objective variable, including the mean value and variance. Besides, robust design can be used to improve the quality and reliability by reducing the functional variation of a system without removing the causes of variation [25]. More details of robust design can be found in the literature [26]. The formulation for implementing the proposed method is given as follows: 𝑛 𝑖
5.2. Optimization Results. Figure 8 shows the Pareto optimal front of the three objectives: head, efficiency, and NPSHr, among which head and efficiency are required to be maximized in order to improve the economic benefits while NPSHr should be minimized for the sake of enhancing cavitation resistance. Both Figures 8(a) and 8(b) present multilinearity more or less, while Figure 8(c) exhibits nonlinearity. In addition, some objectives in this triobjective optimization problem conflict with each other. As shown in Figures 8(a) and 8(c), when a beneficial choice leans to one
2
Minimize 𝐺 = ∑𝑎𝑖 [𝑏𝑖 (𝜇𝑖 − 𝑓𝑖 ) + 𝑑𝑖 𝜎𝑖2 ] Subject to 𝑓1 ≥ 𝐿 1 𝑓2 ≥ 𝐿 2 𝑓3 ≤ 𝐿 3 ,
(11)
where 𝜇𝑖 and 𝜎𝑖2 are the mean values and variance of 𝑓𝑖 in each point of the Pareto optimal solutions, respectively, and 𝐿 1 , 𝐿 2 , and 𝐿 3 are initial values of 𝑓1 , 𝑓2 , and 𝑓3 , respectively. In order to obtain 𝜇𝑖 and 𝜎𝑖2 of 𝑓𝑖 in each point of the Pareto optimal front, the decision variables 𝛼1 , 𝑟1 , 𝛼2 , 𝑟2 , and 𝑙 are assumed to obey normal distribution and Monte Carlo
Downloaded from ade.sagepub.com by guest on August 26, 2015
8
Advances in Mechanical Engineering 95
145 140 Head (m)
Efficiency (%)
90 85 80
130
75 70
135
126 128 130 132 134 136 138 140 142 144 146 Head (m)
125
3
3.5
4
4.5
NPSHr (m)
(a) Efficiency head
(b) Head-NPSHr
NPSHr (m)
4.5
4
3.5
3 70
75
80 85 Efficiency (%)
90
95
(c) NPSHr efficiency
Figure 8: Pareto fronts for the triobjective optimization problem.
(a) Optimized impeller
(b) Common impeller
Figure 9: The optimized impeller and common impeller.
method is carried out with 1000 simulations. Exemplifying by point 𝐴(𝑥1 , 𝑥2 , 𝑥3 , 𝑥4 , 𝑥5 ), it is one point of the Pareto set; 𝑦1 , 𝑦2 , and 𝑦3 are the values of the objective functions 𝑓1 , 𝑓2 , and 𝑓3 with respect to point 𝐴(𝑥1 , 𝑥2 , 𝑥3 , 𝑥4 , 𝑥5 ), respectively. In the Monte Carlo simulation of point 𝐴(𝑥1 , 𝑥2 , 𝑥3 , 𝑥4 , 𝑥5 ), the mean values are 𝑥1 , 𝑥2 , 𝑥3 , 𝑥4 , 𝑥5 , respectively, and the standard deviations of decision variables are all set as 5% of the corresponding mean values. The tradeoff optimal point is the point that makes the value of index function 𝐹 minimum, which indicates the lowest functional variation that can guarantee the robust of the selected point. Table 5 lists the values of decision variables at the optimal point and summarizes the rate of change between the robust results and the nondominated points with highest head and best efficiency, but smallest NPSHr.
The results indicate that the change rate of NPSHr is larger than that of head and efficiency. Figure 9 shows the optimized impeller and a common impeller, which are both from the Pareto solution set. It is clear that the impeller configurations are different and that the optimized impeller has a better comprehensive performance owing to the robust design. 5.3. Verification of Prototype Test. Prototype test plays an indispensable role in the design of centrifugal pump. However, both prototype pump and experiment equipment are expensive. Hence, it is essential to develop centrifugal pumps with best comprehensive performance at a minimum cost. Figure 10 shows the prototype of centrifugal pump corresponding to the tradeoff point, and Figure 11 shows
Downloaded from ade.sagepub.com by guest on August 26, 2015
Advances in Mechanical Engineering
9 Table 5: Results of optimization.
𝛼1 (deg) 10.36
𝑟1
𝛼2 (deg)
𝑟2
𝑙 (mm)
ΔHeadmax (%)
Δ𝜂max (%)
ΔNPSHrmin (%)
0.9051
77.16
0.7957
174.33
1.22
8.35
26.33
Table 6: Results of CFD simulation, RBF predictor, and test experiment.
CFD Predictor Experiment
Figure 10: The prototype of centrifugal pump corresponding to the tradeoff point. Electric machinery
Discharge outlet
Centrifugal pump
Test bench
Pool
Suction inlet
Head (m) 145.26 143.30 144.78
Efficiency (%) 85.71 86.97 87.84
effective on improving pump performance at a minimum cost. The RBF metamodels of pump performance indexes were established, and the Pareto optimal solutions were obtained using MOEA/D, which revealed the conflict among the objective functions of head, efficiency, and NPSHr. The RBF metamodels agreed with the CFD numerical simulation well at the test points, and the CFD simulation results showed the significant influence of decision variables on the behavior of centrifugal pump. A tradeoff point was selected from the Pareto optimal set using a robust design method based on Monte Carlo simulation, and the prototype pump was produced according to a set of decision variables provided by the tradeoff point. Based on the prototype pump, the results of CFD simulation and RBF predictor were in accord with those of experiment test.
Figure 11: Schematic layout of the prototype pump test.
Nomenclature the schematic layout of the prototype pump test. As shown in the figures, the prototype pump is fixed on the test bench, and the electric machinery drives the impeller rotating continuously. The water in the test pool is transferred from the suction inlet to discharge outlet circularly, and the basic test parameters are kept the same with those of CFD simulation in Table 2. Table 6 shows the CFD simulation, RBF predictor, and test values at the tradeoff optimal point, which indicates that the values at the tradeoff point obtained by the CFD simulation, RBF metamodels, and experimental test agree with each other well. The decision variables of the optimized pump can be obtained only by one experiment test. Considering the expensive cost (time and money) of preparing the prototype pump and experiment equipment, the proposed optimization design method in this paper can reduce the consumption to minimum.
6. Conclusions This paper presented a new optimization design method for double suction centrifugal pump, and it was proved to be
NPSHr (m) 3.85 3.79 3.71
𝑎, 𝑏, 𝑑: 𝑐: 𝐸: 𝑓𝑖 : 𝐹(x): 𝑔: 𝐺: 𝐻: 𝐿: 𝑀: 𝑛: 𝑁: 𝑃: 𝑄: 𝑇: 𝑢: 𝑥: x: 𝑧: 𝜌: 𝜂: 𝜆:
Weight coefficient User-defined constant Total energy of per unit mass (J/kg) 𝑖th objective function Vector consisting of objective functions Standard acceleration of gravity (m/s2 ) Index function Head (m) Initial value of objective function Number of test points Number of decision variables Number of sample points Pressure (Pa) Flow rate (m3 /h) Number of the weight vectors Velocity (m/s) Decision variable Vector consisting of decision variables Position head (m) Fluid density (kg/m3 ) Efficiency Unknown interpolation coefficient
Downloaded from ade.sagepub.com by guest on August 26, 2015
10 Φ: 𝜅: 𝜇: 𝜎2 : ‖ ‖:
Advances in Mechanical Engineering Radial function Distribution index Mean value Variance Euclidean norm.
Conflict of Interests The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment This research was supported by the National Natural Science Foundation of China (no. 11172108). This financial support is gratefully acknowledged.
References [1] B. Jafarzadeh, A. Hajari, M. M. Alishahi, and M. H. Akbari, “The flow simulation of a low-specific-speed high-speed centrifugal pump,” Applied Mathematical Modelling, vol. 35, no. 1, pp. 242– 249, 2011. [2] J. F. G¨ulich, Centrifugal Pumps, Springer, New York, NY, USA, 2nd edition, 2010. [3] S. Pierret and R. A. van den Braembussche, “Turbomachinery blade design using a Navier-Stokes solver and artificial neural network,” Journal of Turbomachinery, vol. 121, no. 2, pp. 326– 332, 1999. [4] J. S. Anagnostopoulos, “A fast numerical method for flow analysis and blade design in centrifugal pump impellers,” Computers & Fluids, vol. 38, no. 2, pp. 284–289, 2009. [5] S. Derakhshan, B. Mohammadi, and A. Nourbakhsh, “The comparison of incomplete sensitivities and Genetic algorithms applications in 3D radial turbomachinery blade optimization,” Computers & Fluids, vol. 39, no. 10, pp. 2022–2029, 2010. [6] L. Zhou, W. Shi, and S. Wu, “Performance optimization in a centrifugal pump impeller by orthogonal experiment and numerical simulation,” Advances in Mechanical Engineering, vol. 2013, Article ID 385809, 7 pages, 2013. [7] K. Deb, Multi-Objective Optimization Using Evolutionary Algorithms, John Wiley & Sons, 2001. [8] A. Konak, D. W. Coit, and A. E. Smith, “Multi-objective optimization using genetic algorithms: a tutorial,” Reliability Engineering & System Safety, vol. 91, no. 9, pp. 992–1007, 2006. [9] A. Oyama and M.-S. Liou, “Multiobjective optimization of rocket engine pumps using evolutionary algorithm,” Journal of Propulsion and Power, vol. 18, no. 3, pp. 528–535, 2002. [10] H. Safikhani, A. Khalkhali, and M. Farajpoor, “Pareto based multi-objective optimization of centrifugal pumps using CFD, neural networks and genetic algorithms,” Engineering Applications of Computational Fluid Mechanics, vol. 5, no. 1, pp. 37–48, 2011. [11] A. Papierski and A. Błaszczyk, “Multiobjective optimization of the semi-open impeller in a centrifugal pump by a multilevel method,” Journal of Theoretical and Applied Mechanics, vol. 49, no. 2, pp. 327–341, 2011. [12] R. L. Hardy, “Multiquadric equations of topography and other irregular surfaces,” Journal of Geophysical Research, vol. 76, no. 8, pp. 1905–1915, 1971.
[13] A. A. Mullur and A. Messac, “Extended radial basis functions: more flexible and effective metamodeling,” AIAA Journal, vol. 43, no. 6, pp. 1306–1315, 2005. [14] R. Jin, W. Chen, and T. W. Simpson, “Comparative studies of metamodelling techniques under multiple modelling criteria,” Structural and Multidisciplinary Optimization, vol. 23, no. 1, pp. 1–13, 2001. [15] J. Horn, N. Nafpliotis, and D. E. Goldberg, “Multiobjective optimization using the niched pareto genetic algorithm,” IlliGAL Report 93005, 1993. [16] N. Srinivas and K. Deb, “Muiltiobjective optimization using nondominated sorting in genetic algorithms,” Evolutionary Computation, vol. 2, no. 3, pp. 221–248, 1994. [17] E. Zitzler and L. Thiele, “Multiobjective evolutionary algorithms: a comparative case study and the strength Pareto approach,” IEEE Transactions on Evolutionary Computation, vol. 3, no. 4, pp. 257–271, 1999. [18] K. Deb, A. Pratap, S. Agarwal, and T. Meyarivan, “A fast and elitist multiobjective genetic algorithm: NSGA-II,” IEEE Transactions on Evolutionary Computation, vol. 6, no. 2, pp. 182– 197, 2002. [19] S. Shan and G. G. Wang, “An efficient Pareto set identification approach for multiobjective optimization on black-box functions,” Journal of Mechanical Design, vol. 127, no. 5, pp. 866–874, 2005. [20] Q. Zhang and H. Li, “MOEA/D: a multiobjective evolutionary algorithm based on decomposition,” IEEE Transactions on Evolutionary Computation, vol. 11, no. 6, pp. 712–731, 2007. [21] R. V. Tappeta and J. E. Renaud, “Interactive multiobjective optimization procedure,” AIAA journal, vol. 37, no. 7, pp. 881– 889, 1999. [22] K. A. Kaupert and T. Staubli, “The unsteady pressure field in a high specific speed centrifugal pump impeller—part I: influence of the volute,” Journal of Fluids Engineering, vol. 121, no. 3, pp. 621–626, 1999. [23] J. G. Wang and G. R. Liu, “A point interpolation meshless method based on radial basis functions,” International Journal for Numerical Methods in Engineering, vol. 54, no. 11, pp. 1623– 1648, 2002. [24] P. Rajaguru, S. Stoyanov, H. Lu, and C. Bailey, “Application of Kriging and radial basis function for reliability optimization in power modules,” Journal of Electronic Packaging, vol. 135, no. 2, Article ID 021009, 13 pages, 2013. [25] S. Chen, W. Chen, and S. Lee, “Level set based robust shape and topology optimization under random field uncertainties,” Structural and Multidisciplinary Optimization, vol. 41, no. 4, pp. 507–524, 2010. [26] W. Chen, J. K. Allen, K.-L. Tsui, and F. Mistree, “A procedure for robust design: minimizing variations caused by noise factors and control factors,” Journal of Mechanical Design, vol. 118, no. 4, pp. 478–485, 1996.
Downloaded from ade.sagepub.com by guest on August 26, 2015