Journal of Applied Mathematics

Journal of Applied Mathematics / 2013 / Article
Special Issue

Mathematical and Numerical Modeling of Flow and Transport 2013

View this Special Issue

Research Article | Open Access

Volume 2013 |Article ID 136195 |

Qiaorui Si, Shouqi Yuan, Jianping Yuan, Chuan Wang, Weigang Lu, "Multiobjective Optimization of Low-Specific-Speed Multistage Pumps by Using Matrix Analysis and CFD Method", Journal of Applied Mathematics, vol. 2013, Article ID 136195, 10 pages, 2013.

Multiobjective Optimization of Low-Specific-Speed Multistage Pumps by Using Matrix Analysis and CFD Method

Academic Editor: Mohamed Fathy El-Amin
Received16 Jun 2013
Revised22 Aug 2013
Accepted22 Aug 2013
Published10 Dec 2013


The implementation of energy-saving and emission-reduction techniques has become a worldwide consensus. Thus, special attention should be provided to the field of pump optimization. With the objective of focusing on multiobjective optimization problems in low-specific-speed pumps, 10 parameters were carefully selected in this study for an (310) orthogonal experiment. The parameters include the outlet width of the impeller blade, blade number, and inlet setting angle of the guide vane. The numerical calculation appropriate for forecasting the performance of multistage pumps, such as the head, efficiency, and shaft power, was analyzed. Results were obtained after calculating the two-stage flow field of the pump through computational fluid dynamics (CFD) methods. A matrix method was proposed to optimize the results of the orthographic experiment. The optimal plan was selected according to the weight of each factor. Calculated results indicate that the inlet setting angle of the guide vane influences efficiency significantly and that the outlet angle of blades has an effect on the head and shaft power. A prototype was produced with the optimal plan for testing. The efficiency rating of the prototype reached 58.61%; maximum shaft power was within the design requirements, which verifies that the proposed method is feasible for pump optimization.

1. Introduction

Pumps are widely utilized in various fields of the national economy [1]. Where there is water, there are pumps. Low-specific-speed centrifugal pumps are commonly employed in sewage treatment and well pumping because of their small flow and high head. They consume large amounts of energy and have a potential for energy saving. Low efficiency and overload tendency in long-term large-flow operations are the two obstacles in designing such pumps; the impeller diameter of low-specific-speed centrifugal pumps is relatively large, and the flow channel is long and narrow, resulting in significant disc friction and hydraulic loss. Thus, the efficiency of low-specific-speed centrifugal pumps is relatively low; the shaft power curves of such pumps increase sharply. The greater flow method is usually adopted in the design of low-specific-speed pumps to enhance flow efficiency. Thus, the ratio of macroshaft power to designed power becomes much higher than that of common centrifugal pumps. A motor whose export pipeline has no valve burns out easily [2, 3]. Pumps are usually designed in multistages to increase the amount of pressure because of limitations imposed by costs or diameter, such as the case of well pumps. Increasing the single-stage head may decrease the number of stages to save costs, energy, and materials. High efficiency in maintaining the head and low power are the two goals of multistage pump designs. Thus, the study on multitarget optimal pump designs is important in implementing energy saving and emission reduction.

Mathematical methods, such as neural networks, orthographic experiments, genetic algorithms, and grey theory, are often adopted in multitarget optimization design [4, 5]. An orthographic experiment design can optimize the experimental conditions to achieve the target with fewer experiments [6, 7]. It has been widely adopted in pump design. Shouqi et al. [8] employed an orthogonal tableto thoroughly study the effect of impeller geometric parameters and the throat area on pump performance. He proposed a practical framework for a centrifugal pump without overload and its design method. Wang et al. [9] focused on the effect of main impeller geometric parameters on deep-geometric performance via orthographic experiments. Shen et al. [10] designed nine models of complex impeller centrifugal pumps with orthographic experiments and discovered the effect order of geometric parameters on pump performance. Zhou et al. [11] identified the key factors that affect the performance of guiding blades by studying conduit guiding blades and by conducting orthographic experiments and designing high-performance guiding blades. The orthogonal experiment method is therefore a multivariate and multilevel optimization method suitable for pump design. The studies mentioned above examined the influence of parameters on the impeller or guiding blades separately; however, the factors were not enough to completely reflect the effect of geometric parameters on the experiment. Furthermore, the multiobjective test was usually transformed into a single-objective test in the experimental analysis. Plans were then comprehensively selected through single-index analyses, including intuitive analysis and variance analysis. These methods neglect the significance and differences of indexes.

QS10-68, a typical multistage electric submersible pump (ESP), was regarded as an example in the present study. A method that combines orthogonal experimentation and numerical simulation was applied for optimization. A simple matrix method was also introduced to calculate the effect of factors on each index and to directly decide the factor order according to weight, which would perfectly solve the issue of selecting an optimal plan in an orthographic experiment design. The effect of 10 geometric parameters of the impeller and guiding vane on the multistage head, efficiency, and shaft power was also studied. Lastly, a set of optimum geometric parameters was obtained based on weight calculation. A prototype that employs the optimal plan was then tested for verification.

2. Pump Model

The structure of the multistage ESP is displayed in Figure 1. The figure shows a multistage centrifugal pump operating in a vertical position [12]. The pump shaft is connected to the protector by a mechanical coupling at the bottom of the pump. Fluids from the well enter the pump through an intake screen and are lifted by the pump stages. Produced liquids, after being subjected to great centrifugal forces caused by the high rotational speed of the impeller, lose their kinetic energy in the diffuser. Kinetic energy is converted to pressure energy in the diffuser. The design parameters for this pump are as follows: flow rate  m3/h, total head  m, speed  r/min, four stages, and the efficiency of the necessary electric motor is 4 kW.

3. Matrix Analysis Model and Experimental Scheme

In a multiobjective optimization problem, engineers adopt an orthographic experimental method as a solvent. The orthographic experimental method is a scientific method utilized to arrange and analyze multiple factor experiments by means of a table based on orthogonal principles. Through scientific arrangement and analysis of the result, a study can discover the ideal production conditions and techniques. However, problems, including large calculated amounts and confirming the weights unreasonably, exist in the multitarget orthogonal test method. A matrix analysis model was presented to solve this problem. In this model, a three-layer structure and the layer structure matrix of the orthogonal test are established first. The weight matrix of the test index is then calculated by multiplying the matrix of each layer. The weighs of the factors and levels that affect the tests results are calculated. Finally, the optimal plan and the importance order of the factors that influence the test index values are determined according to the weights.

3.1. Matrix Analysis Model

An orthographic experiment was designed (ignoring the interaction for the first). The three-layer structural model based on data structure was established. The first layer in Table 1 is the investigation index, the second layer is the factor layer, and the last one is the level layer. The data of each layer determine the definition of the matrix as follows.

First levelExperimental index

Second levelFactor A 1Factor A 2Factor A l
Third level

Definition 1. The following is the matrix of the experimental investigation index. The study order is . The matrix is established on the conditions that factors exist in the experiment, each factor has levels, is the average index of factor at level, and the investigation index is good when it is high. A low investigation index is considered good when the study establishes a matrix supposing .

Definition 2. The matrix of the factor layer: to build (2) matrix.

Definition 3. The matrix of the level layer: the range of in the orthographic experiment is set as . The study order is . Equation (3) matrix is then established.

Definition 4. The weight matrix that may affect the test index:.
In the matrix above mentioned, , refers to the ratio of ’s first level index to the sum of all ’s level indexes. is the ratio of’s range to the sum of all ranges. The product of these two data reflects the effect of the first level of on the index and range of . The other factors and levels are identical. The weight of each factor and level was obtained after calculation. The optimal plan and the factor order in the index were established based on the weight value.

3.2. Experimental Scheme

According to Euler equations, theoretical head is where (m/s) is peripheral speed at the inlet; is the peripheral velocity component of the blade inlet; in the straight cone suction chamber; (m/s) refers to the circular velocity at the impeller outlet; is acceleration because of gravity; is the stodala slip coefficient; represents blade numbers; is the theoretical flow; is the cross-section area of the impeller outer; is the blade expelling coefficient;is the width of the blade outlet; is the blade outlet angle. , the hydraulic power input, can be deduced by the following equation:

In (6), (kg/m3) refers to water density. The equation shows that both theoretical head and hydraulic input power are determined by the impeller. However, rolling flow causes , additional shaft power, as it enters the guiding blade channel. Hence, shaft power and pump efficiency . Pump shaft power and efficiency are determined by both the impeller and the guide vane in hydraulics.

The following are the main geometric parameters of the model: width of the blade outlet , suction diameter of impeller , inclination of the back cover board , inlet attack angle , inlet angle in axial plane , blade outlet angle , blade number, outer diameter of the back cover , blade thickness , inlet angle of the guide vane , outlet angle of the guide vane , import width of the guide vane , blade number of the guide vane , axial length of the guide vane , inlet diameter of the guide vane , and outlet diameter of the guide vane . is equal to 90° to eliminate the rotation component. The last five factors were determined by the basic geometric parameters of the impeller. The first 10 geometric parameters were selected as testing factors.

According to the specific speed equation, where refers to the flow and refers to the head of the pump. The design-specific speed of the pump is low, namely, 65.6, because it is lower than 80. Relative parameters can be obtained with the velocity-coefficient method under structure demand; these factors are shown in Table 2. Three levels were selected in each factor. The experimental scheme, shown in Table 3, is based on the orthogonal table. In the 27 group tests, different numerals appear at the same frequency in any column, and any two-digit number appears equivalent in any row. After deciding the experimental plan, tests were implemented in numerical order without any arbitrary changes.

Level Factor


Ad 2min/mm, Bb 2/mm, C /(°), Dz, ES 2/mm, F /(°), Gd 1/mm, H /(°), I /(°), J /(°).

Test numberA/mmB/mmC/(°)D/E/mmF/(°)G/mmH/(°)I/(°)J/(°)


4. Numerical Simulation

In traditional orthogonal experiments, research procedures involve manufacturing prototypes, performance tests, and results analyses. However, creating prototypes from the 27 groups of impellers and guide vanes would be a waste of time and money; moreover, a large number of prototype tests would inevitably cause significant manufacturing and test errors. The prediction of pump performance has become possible in engineering applications owing to the rapid development of computational fluid dynamics (CFD) [13]. Reasonable CFD calculations reflect the actual internal flow of the pump and accurately predict the pump head performance, efficiency, shaft power, and so forth at specific conditions [14, 15]. Therefore, the use of CFD technology is more feasible than the use of a prototype in pump design optimization and in the establishment of a preliminary forecast of pump performance.

4.1. Governing Equations and Boundary Conditions

The flow field information of pumps can be described by Navier-Stoke equations. A numerical simulation is performed to solve the following governing equations [16]:

Fluent 6.2 was utilized in this study. The 3D unsteady flow of centrifugal pumps was calculated based on Reynolds-averaged equations that resemble the standard renormalization group turbulence model. Velocity, turbulent kinetic energy, and eddy viscosity coefficients were provided as a first-order upwind scheme. The velocity at the inlet and free flow at the outlet were selected with boundary conditions. A solid wall was set as the no-slip condition, and a smooth wall condition was employed as the near-wall function. The convergence precision was set to 10−5.

4.2. Computational Domain

The flow domain of the submersible pump is composed of the inlet, a multistage distorting impeller, a multistage space diffuser, and an outlet. Its flow pattern is more complex than that of a single-stage pump. A rotational flow exists at the impeller inlet, but the first stage and flow in the channel are similar. Demands on computer performance would be made if all stages are considered in the calculation. Therefore, selecting only the appropriate stages is important. Four kinds of computational domains in different stages are simulated for scheme. The relevant results are presented in Tables 4, 5, and 6.

Stage1st head 2nd head 3rd head 4th head Total head


Stage1st Shaft power2nd Shaft power3rd Shaft power4th Shaft powerTotal shaft power


Stage1st efficiency2nd efficiency3rd efficiency4th efficiencyTotal efficiency


Tables 4 and 5 indicate that the head and shaft power of each stage in the pump are very similar that the two-stage values can be regarded as the corresponding values of the pump. Table 6 shows that only a few differences exist in terms of efficiency; however, many differences were noted between the first stage and other stages after the second. No changes were observed after stage two. In summary, it is appropriate to select two full-flow models to analyze pump performance. Similarly, the use of the head and shaft power of two stages as well as the efficiency of stage two to predict performance is also feasible. The computational domain composed of the inlet, the distorting impeller of two stages, the space diffusers, and the outlet is shown in Figure 2. The first pump water chamber is not shown in the figure.

4.3. Model Meshing

The quantity and quality of the grid are two important factors that affect computation accuracy and duration. The unstructured tetrahedral mesh provided by GAMBIT, the professional software for grid generation, was utilized on the entire basin. The mesh has strong adaptability. Encryption and nonequidistance were processed in the near-wall region. Five different mesh sizes were selected to determine the appropriate number of grids and to conduct a grid-independent analysis. Schemeat was regarded as an example. The calculated results are shown in Table 7.

ParameterMesh size/mm

Efficiency /%61.8261.8161.7761.4660.53
Shaft Power P/kW1.7161.7061.7111.6941.686

Table 7 shows that the difference between pump head does not exceed 0.5%, efficiency does not exceed 0.06%, and shaft power does not exceed 0.3% when the mesh size is smaller than 2.0 mm. , which was utilized to examine the closest node to the surface, is not greater than 150. Thus, the criteria of standard renormalization group turbulence model calculation were met. Therefore, numerical simulation grid size of 2.0 mm was selected after considering computation accuracy and duration. The grid view of computational domain is shown in Figure 3.

One of the convergence plots is shown in Figure 4; this plot can be used to justify the ultimate choice of grid for the analysis. According to Figure 4, all the convergence precisions reached 10−5, which satisfied the requirement of the calculation.

5. Results Analysis

Given that the shaft power of low-specific-speed centrifugal pumps increases rapidly with the enlargement of flow, the shaft power at was selected as the monitoring value for overload judging. An electric pump works safely at the range of to . If the shaft power at is less than the motor power, then the pump will not experience an overload. A total of 27 numerical simulations were performed. The results of two-stage head under the rated condition, shaft power at , and efficiency of stage two under the rated condition are shown in Table 8.

Test numberH/mη/%P/kW


H (m) is two-stage head under Q; η (%) is the efficiency of stage two under Q; P (kW) is two-stage shaft power at 1.5 .
5.1. Visual Analysis

Orthographic experimental data were analyzed. With the first column of rated point head index as an example, Where , and in the equations stand for the average efficiency values at the rated point of factor ., the range, is introduced to reflect the effect of various factors on the target value. A large range signifies an important factor whose change greatly affects the objective. If the range is small, then the factor has minimal influence and is usually not important.was obtained by the subtraction of the minimum from the maximum among , ,  and in each column. The calculation methods for the other nine columns were similar to the first one. The efficiency under the rated condition and shaft power at were obtained. , , ,  and of each factor are listed in Table 9.

IndexA B C D E F G H I J

 k 134.97734.40133.86234.62535.78435.47335.28135.36635.41235.087
 k 235.62735.35535.83136.04535.29535.25235.44735.65635.43436.535
 k 335.51336.36236.42435.44735.03835.39335.38935.09635.27134.495
 s 0.6501.9612.5621.4200.7460.2210.1670.5600.1632.040
 k 160.79560.93561.05961.11660.63060.73960.28560.53060.67860.600
 k 260.66160.83860.75460.96560.74760.75560.92060.89760.71361.795
 k 360.57160.25460.21459.94660.65060.53260.82260.60160.63659.632
 s 0.2240.6810.8441.1700.1170.2230.6350.3670.0772.163
 k 11.9631.9271.8551.9362.0332.0272.0092.0152.0172.018
 k 22.0332.0072.0452.0532.0172.0002.0082.0222.0312.025
 k 32.0382.1012.1342.0451.9852.0072.0171.9981.9871.992
 s 0.0750.1740.2790.1160.0480.0270.0090.0240.0430.033

Visual analysis indicates that the order of the importance of the index for the head is , and its optimal plan is . The order of importance of the index for efficiency is , and its optimal plan is . Similarly, the order of importance for shaft power is, and its optimal plan is . The integrated balance method even drawing was utilized to select the optimal program, resulting in large amounts of calculation and difficulty in establishing choices for multitarget optimization. If the matrix method introduced in this study is selected, an optimal plan can be immediately established because only three investigation indexes are calculated. The optimal program is determined by the weight value.

5.2. Weight Matrix Method

Weight matrix analysis is necessary to calculate the weight matrix of the investigation index. The weight calculation of the head was regarded as an example; high values were considered good. Thus, , , , and .

In general, higher efficiency is also considered good. Hence, , , . Furthermore, low shaft power is desirable. Thus, , , and . The general weight matrix in (14) is obtained after substituting the above calculations into (4) as follows:

The calculations provide the weight of the result of each factor with three levels. Hence, the optimal plan for orthographic experiments is determined quickly. The optimal plan is , and the importance order of the index is .

6. Experimental Verification

The optimal model, , was prototyped for the experiment to test the optimized programs. The geometric parameters of the hydraulic model were as follows: width of the blade outlet , inclination of back cover board , inlet attack angel , inlet angle in axial plane , blade outlet angle , blades number , outer diameter of back cover , blade thickness , inlet angle of guide vane , and axial length of the guide vane . The test was performed in accordance with the national standards GB/T 12785-2002 of China. Figure 5 presents the results of the experiment and the entire flow field simulation.

Comparison of the simulation and experimental results indicates that the simulated head and efficiency are slightly higher than those of the experiment; however, the simulated power is slightly lower. Error analysis indicates the following:the numerical simulation did not consider the loss of leakage at the oral ring and between stages;small differences exist because of limitations imposed by casting accuracy, casting errors, casting model (particularly the impeller), and simulated model; these differences can cause disparity between experimental and simulated values. However, the error is less than 5%, and both experimental and simulated values exhibit similar trends when the flow changes. Hence, the simulation methods employed in this study can meet the demand of practical usage in programs for design optimization.

The efficiency of QS10-68/4 based on the national standard of China is 51%. After optimization, the pump’s rated head is 68.9 m and rated efficiency is 58.61%, which is approximately 7% higher than the national standard. The maximum shaft power is 3.83 kW. The performance of nonoverload also meets the demand, which is 4 kW. All of the above data indicate that optimization was successful.

7. Conclusions

The influence of 10 factors and three levels on the head, efficiency, and shaft power was analyzed by conducting 27 orthographic experiments through numerical simulation. Weight matrix analysis was then performed to determine the optimal model. The principal conclusions drawn from this research are as follows.(1)The simulation results indicate that the inlet setting angle of the guide vane significantly affects the efficiency and that the outlet angle of blades has a significant effect on the head and shaft power.(2)The importance order of the index for the head is , and the geometric parameters of optimal plan are as follows: , , , , , , , , , and . The importance order of the index for efficiency is , and the geometric parameters of optimal plan are as follows: , , , , , , , , , and . For shaft power, the importance order of the index is , and the geometric parameters of optimal plan are as follows: , , , , , , , , , and .(3)By performing weight matrix analysis, optimal model was obtained based on the influence of factors and levels on the head, efficiency, and shaft power. The geometric parameters of the hydraulic model were as follows: , , , , , , , , , and . The performance of the prototype in the test indicates that efficiency at the rated point is 7% higher than the national standard of China. This result verifies the potential use of the orthographic design combined matrix analysis with numerical simulation in pump optimization.(4)The changes in the head, efficiency, and shaft power determined from the simulation are in accordance with the test results. The error between numerical simulation and test result is less than 5% in the rated flow, which further verifies the possibility of forecasting the performance of multistage pumps through numerical simulation.


This work was financially supported by the State Key Program of National Natural Science of China (Grant no. 51239005) and the National Science & Technology Pillar Program (Grant no. 2011BAF14B04) of China. This research was also sponsored by the Research and Innovation Project for College Graduates of Jiangsu Province (CXZZ12_0679).


  1. F. G. Johann, Centrifugal Pumps, Springer, New York, NY, USA, 2008.
  2. Y. Shouqi, C. Wuling, and L. Shiying, “Theory and design method of non-overload centrifugal pumps,” Chinese Journal of Mechanical Engineering, vol. 5, no. 4, pp. 252–260, 1992. View at: Google Scholar
  3. Z. peifang, Y. Shouqi, and X. Jianqiang, “Present status and development trends of low-specific-speed centrifugal pumps,” Drainge and Irrigation Machinery, vol. 21, no. 6, pp. 42–45, 2003. View at: Google Scholar
  4. M. A. Moreno, J. I. Córcoles, D. A. Moraleda, A. Martinez, and J. M. Tarjuelo, “Optimization of underground water pumping,” Journal of Irrigation and Drainage Engineering, vol. 136, no. 6, pp. 414–420, 2010. View at: Publisher Site | Google Scholar
  5. M. Moradi-Jalal, S. I. Rodin, and M. A. Marino, “Use of genetic algorithm in optimization of irrigation pumping stations,” Journal of Irrigation and Drainage Engineering, vol. 130, no. 5, pp. 357–365, 2004. View at: Publisher Site | Google Scholar
  6. A. S. Hedayat, N. J. A. Sloane, and J. Stufken, Orthogonal Arrays, Theory and Applications, Springer, New York, NY, USA, 1999. View at: Publisher Site | MathSciNet
  7. Z. Lixue, S. Yaozhen, L. Qingsheng et al., “Mathematic model of seed peeling efficiency with critical structure parameters based on orthogonal multinomial tests,” Transactions of the CSAE, vol. 19, no. 5, pp. 22–26, 2003. View at: Google Scholar
  8. Y. Shouqi, Z. Jinfeng, Y. Jianping et al., “Orthogonal experimental study effect of main geometry factors of splitter blades on pump performance,” And Irrigation Machinery, vol. 26, no. 2, pp. 1–5, 2008. View at: Google Scholar
  9. H. Wang, W. Shi, W. Lu, L. Zhou, and C. Wang, “Optimization design of deep well pump based on latin square test,” Transactions of the Chinese Society of Agricultural Machinery, vol. 41, no. 5, pp. 56–63, 2010 (Chinese). View at: Publisher Site | Google Scholar
  10. Y. Shen, S. Yuan, W. Lu, J. Zhang, and J. Yuan, “Orthogonal test design method based on numerical simulation for non-overload centrifugal pump with complex impeller,” Transactions of the Chinese Society of Agricultural Machinery, vol. 41, no. 9, pp. 22–26, 2010 (Chinese). View at: Publisher Site | Google Scholar
  11. L. Zhou, W. Shi, W. Lu, R. Xu, and C. Wang, “Orthogonal test and optimization design of submersible pump guide vanes,” Journal of Drainage and Irrigation Machinery Engineering, vol. 29, no. 4, pp. 312–342, 2011. View at: Publisher Site | Google Scholar
  12. W. Shi, T. Jiang, W. Cao, W. Lu, D. Zhang, and B. Chen, “Optimal design and experiment on a high-head non-overload submersible sewage pump,” Transactions of the Chinese Society of Agricultural Engineering, vol. 27, no. 5, pp. 151–155, 2011. View at: Publisher Site | Google Scholar
  13. S. N. Shukla and J. T. Kshirsagar, “Numerical experiments on a centrifugal pump,” in Proceedings of the 2002 ASME Joint U.S.-European Fluids Engineering Division Conference (FEDSM '02), pp. 709–720, July 2002. View at: Publisher Site | Google Scholar
  14. R. K. Byskov, C. B. Jacobsen, and N. Pedersen, “Flow in a centrifugal pump impeller at design and off-design conditions. Part II: large eddy simulations,” Journal of Fluids Engineering, vol. 125, no. 1, pp. 73–83, 2003. View at: Publisher Site | Google Scholar
  15. R. Spence and J. Amaral-Teixeira, “A CFD parametric study of geometrical variations on the pressure pulsations and performance characteristics of a centrifugal pump,” Computers and Fluids, vol. 38, no. 6, pp. 1243–1257, 2009. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  16. J. H. Ferziger, Computational Methods for Fluid Dynamics, Springer, Berlin, Germany, 2002.

Copyright © 2013 Qiaorui Si 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.

Related articles

No related content is available yet for this article.
 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

No related content is available yet for this article.

Article of the Year Award: Outstanding research contributions of 2021, as selected by our Chief Editors. Read the winning articles.