A New Prediction Method for the Complete Characteristic Curves of Centrifugal Pumps

: Complete pump characteristics (CPCs) are the key for establishing pump boundary conditions and simulating hydraulic transients. However, they are not normally available from manufacturers, making pump station design difﬁcult to carry out. To solve this issue, a novel method considering the inherent operating characteristics of the centrifugal pump is therefore proposed to predict the CPCs. First, depending on the Euler equations and the velocity triangles at the pump impeller, a mathematical model describing the complete characteristics of a centrifugal pump is deduced. Then, based on multiple measured CPCs, the nonlinear functional relationship between the characteristic parameters of the characteristic operating points (COPs) and the speciﬁc speed is established. Finally, by combining the mathematical model with the nonlinear relationship, the CPCs for a given speciﬁc speed are successfully predicted. A case study shows that the predicted CPCs are basically consistent with the measured data, showing a high prediction accuracy. For a pump-failure water hammer, the simulated results using the predicted CPCs are close to that using the measured data with a small deviation. This method is easy to program and the prediction accuracy meets the requirements for hydraulic transient simulations, providing important data support for engineering design.


Introduction
Due to its wide range of flow and head, a centrifugal pump is commonly used in the fields of regional water diversion, agricultural irrigation, industrial cooling cycles and urban water supply and drainage [1]. Its hydraulic performance is directly related to the operational safety and benefits of the project [2]. Therefore, obtaining various pump performance parameters is necessary. As the curves showing the relationships among the discharge, head, speed and torque under all possible operating conditions of the pump, the CPCs are a key factor for establishing pump boundary conditions, which provides important data support for the hydraulic transient analysis and the optimum design of the water conveyance system in the pump station [3][4][5]. However, because the impeller model test is time-consuming and laborious, pump manufacturers generally provide the pump characteristic data in the normal pumping zone but rarely provide the CPCs for all operating zones [6,7]. In addition, engineering practice has shown that in the absence of measured CPCs, simply applying the CPCs with an approximate specific speed may sometimes produce large deviations, which will adversely affect the operational safety of the pumping station project [8,9]. Therefore, obtaining a scientific and reasonable method to predict the CPCs for an arbitrary specific speed is necessary to accurately simulate the hydraulic transients and guide the design of pump stations.
In 1965, Suter and Marchal [10] first proposed the 'Suter transformation', which included the dimensionless head characteristic curve (WH) and torque characteristic curve (WB) under normal and abnormal operating zones. Since then, accurately predicting the WH and WB curves of a centrifugal pump has become a research hotspot. However, due to its complex shape and high nonlinearity, it is difficult to predict the WH and WB curves [11]. To expand the existing data of CPCs, scholars have conducted a series of studies and have achieved fruitful results in the prediction of CPCs. Currently, the main methods used to predict the CPCs can be divided into curve fitting, neural network prediction, and computational fluid dynamics (CFD) prediction.
The curve fitting method constructs a general formula to describe the CPCs at different specific speeds by setting parameters, such as the specific speed and the relative flow angle as independent variables based on existing CPCs data. Izquierdo [12] recommended the use of cubic spline interpolation to fit the CPCs, pointing out that this method can be better utilized for solving the boundary conditions of the pump while ensuring the good fitting effect. Zhu [13] takes the specific speed, relative flow angle and characteristic parameters corresponding to the zero flow condition as independent variables and establishes a new fitting formula for the CPCs by the least squares method. Hu [14] established a threedimensional Cartesian coordinate system with the relative flow angle, specific speed and WH value (or WB value) as independent variables and interpolated them with a bicubic polynomial to develop a three-dimensional surface visual model for predicting the CPCs. The predicted curves basically match the measured data. Lima [15] constructed an equation for fitting the CPCs in a similar triangular series form while using a particle swarm algorithm to optimize the fitting coefficients, and the predicted CPCs have shown a certain degree of reliability. Wan [16] analyzed the internal relationship between the normal performance curves and the CPCs, successfully inverting the 'pumping zone' data in the CPCs, and proved that this inversion method can provide more accurate data for the calculation of the hydraulic transients.
The inherent nonlinearity of the CPCs makes them difficult to predict by traditional system identification methods. A neural network, which is not limited by the nonlinear model and is easy to implement in engineering, shows good adaptability in the prediction. Liu [17] and Lu [18] adopted a three-layer backpropagation (BP) neural network to predict the CPCs for any specific speed based on the three to four measured characteristic data of pumps. Liu [19] used a genetic algorithm (GA) to optimize the number of hidden nodes and the center of the basis function in a radial basis function (RBF) neural network. She proposed a GARBF neural network model for predicting the CPCs, which can take the topology, connection weight, hidden node center and width of the network as a whole for global optimization. To overcome the problems of traditional BP neural networks, which exhibit low learning efficiency and can easily fall into local minima, Han [20] used a combination of the Levenberg-Marquardt (LM) algorithm and a double-hidden layer BP neural network to predict the normal performance curve of pumps, and the results show that the method can reduce the prediction error while shortening the neural network training time. Huang [21] introduced the loss in the inlet pipe, the incidence loss in the pump impeller and the disk friction loss into the BP neural network to further improve the prediction accuracy of the energy performance (i.e., the head, power and efficiency) of centrifugal pumps. Wu [22] built a prediction model of the normal performance curve of an axial flow pump based on an improved RBF neural network, and the prediction results were in good agreement with the experiments.
With the rapid development of modern computer technology, CFD technology has gradually become an effective means to predict the CPCs. CFD technology is used to model the pump in three dimensions, set the operation conditions and calculate the flow field to obtain the simulation results, such as the pump inlet and outlet head, discharge, speed and torque under different operation conditions (including turbine, dissipation, normal pump and reverse dissipation operation conditions). Finally, the mathematical transformation of these simulation results can realize the prediction of the CPCs. Höller [23] numerically simulated the steady-state and transient conditions of a mixed flow pump with CFD software, proposed a method to obtain the four quadrant characteristic curve of a mixed flow pump, and experimentally verified the accuracy of the method. Gros [24] analyzed and compared the four quadrant characteristic curves of a centrifugal pump obtained by experiments and CFD technology. Wang [25] established a three-dimensional simulation model for a double-suction centrifugal pump by using CFD software and proposed a novel three-dimensional internal characteristics method for determining the WH and WB curves of double-suction centrifugal pumps, which greatly improved the accuracy of the water hammer analysis. Muttalli [26] evaluated the normal performance curve of a centrifugal pump using the ANSYS CFX-14 and discussed the effect of cavitation flow phenomenon on the simulation results. Frosina [27] predicted the inverse characteristic of centrifugal pumps using CFD technology and evaluated the deviation between the performance of pumps-as-turbines obtained by different prediction methods and the experimental results.
Since the existing prediction methods for the CPCs rely mostly on sufficient pump design parameters and the number of training samples or derive the prediction from the perspectives of function fitting and trend prediction, the inherent operating characteristics of the pump are ignored, and the mathematical relationships among parameters, such as head, speed and torque, are not revealed. Therefore, this paper takes several sets of the measured CPCs with different specific speeds as the research object. First, based on the Euler equations and the velocity triangles at the pump impeller when the pump is running at a steady state, a mathematical model describing the complete characteristics of a centrifugal pump is deduced in Section 2. Second, multiple measured CPCs at various specific speeds are transformed to obtain the characteristic parameters corresponding to each COP and the functional relationship between the characteristic parameters and the specific speed is established using a nonlinear regression model in Section 3. Then, this relational model is substituted into the mathematical model describing the CPCs, the undetermined coefficients in the model are solved, and the CPCs for a given specific speed are successfully predicted in Section 4. Finally, a comparison with the measured CPCs and a case study are conducted to validate the proposed prediction method in Section 5. Conclusions are presented in Section 6.

Materials and Methods
According to the theorem of moment of momentum, the torque and the theoretical pumping head of a centrifugal pump under steady conditions can be determined by Euler's equation [28,29]: where subscripts 1 and 2 represent the impeller inlet and outlet, respectively; M is the shaft torque (N·m); H e is the theoretical head (m); H is the practical head (m); η = H/(H + KQ 2 ) is the hydraulic efficiency; K is the hydraulic loss coefficient; Q is the discharge (m 3 /s); ω is the angular velocity of rotation (rad/s); V u is the peripheral component of the absolute velocity (m/s); r is the impeller radius (m); ρ is the flow density (kg/m 3 ); and g is the gravitational acceleration (m/s 2 ). Without considering the flow loss at the impeller and assuming that the flow distribution at the impeller inlet or outlet is uniform, the discharge at the impeller inlet can be expressed as [30,31]: where V 1m is the meridional flow velocity at the impeller inlet (m/s), A 1 is the section area at the impeller inlet (m 2 ), and b 1 is the section width at the impeller inlet (m). According to the continuity of water flow, the discharge at the impeller outlet is equal to the flow at the inlet, which is: where V 2m is the meridional flow velocity at the impeller outlet (m/s), A 2 is the section area at the impeller outlet (m 2 ), and b 2 is the section width at the impeller outlet (m). According to the velocity triangle of water flow inside the centrifugal pump impeller shown in Figure 1, V 1u and V 2u in Equation (1) can be written as [32,33]: where U 2 is the peripheral velocity at the impeller outlet (m/s), α 1 is the absolute flow angle at the impeller inlet ( • ), and β 2 is the relative flow angle at the impeller outlet ( • ).
to the flow at the inlet, which is: where V2m is the meridional flow velocity at the impeller outlet (m/s), A2 is the section area at the impeller outlet (m 2 ), and b2 is the section width at the impeller outlet (m). According to the velocity triangle of water flow inside the centrifugal pump impeller shown in Figure 1, V1u and V2u in Equation (1) can be written as [32,33]: where U2 is the peripheral velocity at the impeller outlet (m/s), 1 is the absolute flow angle at the impeller inlet (°), and 2 is the relative flow angle at the impeller outlet (°).
According to the similarity principle, the following dimensionless similarity relationship exists when a pump has similar working conditions [34]: where q = Q/Qr is the relative discharge, h = H/Hr is the relative head, n = N/Nr is the relative rotational speed, m = M/Mr is the relative torque, and subscript r represents the rated condition of the pump. Substituting Equations (4) and (5) into Equation (1) yields: According to the similarity principle, the following dimensionless similarity relationship exists when a pump has similar working conditions [34]: where q = Q/Q r is the relative discharge, h = H/H r is the relative head, n = N/N r is the relative rotational speed, m = M/M r is the relative torque, and subscript r represents the rated condition of the pump. If q/ √ h and m/h are the abscissas and n/ √ h is the ordinate, q/ √ h − n/ √ h and m/h − n/ √ h curves can be drawn, respectively, which can reflect the relationship between q, n, h and m of the pump at different speeds. Substituting q = Q/Q r , n = N/N r , h = H/H r , m = M/M r , ω = 2πN/60, η = H/(H + KQ 2 ) into Equation (6) and eliminating Q, ω, H and M, the following equations can be deduced: where λ 1 = πρr 2 2 Q r N r /30M r , λ 2 = −ρ(cotα 1 /2πb 1 + cotβ 2 /2πb 2 )Q 2 r /M r , ξ 1 = π 2 r 2 2 N 2 r / (900gH r ), ξ 2 = −π(cotα 1 /2πb 1 + cotβ 2 /2πb 2 )Q r N r /(30gH r ), and ξ 3 = −KQ 2 r /H r . The above characterization of the relationship curves among q, h and m under different values of n during steady operating conditions provides the complete characteristics of a centrifugal pump. When a hydraulic transient occurs, the internal flow pattern of the pump is complex, which causes parameters, such as flow rate and torque, to change significantly over time. Analysis of the existing measured CPCs indicates that the corresponding torque is not equal to zero when the discharge is zero. Thus, Equation (9) can be corrected as follows: where λ 3 is the correction factor, which is related to the specific speed of the pump. Equations (9) and (10) constitute a mathematical model describing the complete characteristics of a centrifugal pump, where Equations (9) and (10) are used to determine the flow characteristics (q/ √ h − n/ √ h) and the torque characteristics (m/h − n/ √ h), respectively. The shapes of the complete characteristic curves depend on the parameters (such as λ 1 , λ 2 , λ 3 , λ 1 , ξ 1 , ξ 2 , and ξ 3 ) that are related to the geometric parameters of the impeller and the working condition parameters. However, the above coefficients are unknown before the model test of the pump impeller is carried out [4].
There are two solutions: the direct method and indirect method, and flow charts for these two methods are shown in Figure 2. The direct method refers to the statistics of the geometric parameters and working condition parameters of pumps with different specific speeds, including the impeller diameter, the section width at the impeller inlet and outlet, the absolute flow angle and the relative flow angle. Then, the relationship between these parameters and the specific speeds can be determined by regression analysis; thus, the values of λ 1 , λ 2 , λ 3 , λ 1 , ξ 1 , ξ 2 , and ξ 3 of the pump at a given specific speed can be directly determined according to the definitions of these parameters. Generally, the statistics of working condition parameters in the direct method are more difficult. Therefore, the indirect method can be used.
Since the CPCs are comprised of a series of curves reflecting the relationship between q, n, h and m, the q/ √ h − n/ √ h relation curve and m/h − n/ √ h relation curve can be obtained by mathematical transformation. Then, several COPs are collected on the complete characteristic curves of centrifugal pumps with different specific speeds and the regression model is developed for their dimensionless characteristic parameters (q/ √ h, n/ √ h and m/h) as functions of specific speed. When the specific speed is given, the characteristic parameters corresponding to these COPs on the CPCs can be approximately determined according to the regression model. Thus, the values of λ 1 , λ 2 , λ 3 , λ 1 , ξ 1 , ξ 2 , and ξ 3 can be indirectly determined by substituting the values of q/ √ h, n/ √ h and m/h of COPs into Equations (9) and (10). Considering that both Equations (9) and (10) have only three undetermined coefficients, the values of these coefficients in theory can be determined if the dimensionless characteristic parameters (q/ √ h, n/ √ h and m/h) of three COPs in the corresponding operating region are known. Since the characteristic parameters of the measured COPs are easier to obtain, the indirect method is applied to predict the complete characteristics of a centrifugal pump. COPs in the corresponding operating region are known. Since the characteristic parameters of the measured COPs are easier to obtain, the indirect method is applied to predict the complete characteristics of a centrifugal pump.

Suter Transformation
Because q, n, h and m may change signs and pass through zero during a transient process, it is not convenient to use the q/ √ h − n/ √ h relation curve and m/h − n/ √ h relation curve for the hydraulic transient simulation of a pump station directly. Therefore, Suter and Marchal [10] proposed the following mathematical transformations based on the dimensionless similar characteristics of pumps: where x represents the instantaneous position of pump operation and WH and WB represent the head characteristics and torque characteristics of the pump, respectively. Equations (11) and (12) are valid for any real numbers of n and q except that n and q are equal to zero simultaneously when the pump is shut down. As shown in Figure 3, according to the positive and negative relationship between n and q, the entire operating conditions of the pump can be divided into the turbine zone (q ≤ 0, n < 0, first quadrant), dissipation zone (q < 0, n ≥ 0, second quadrant), normal pump zone (q > 0, n ≥ 0, third quadrant) and reverse dissipation zone (q > 0, n < 0, fourth quadrant) [35]. When the pump is stopped in the case of an accident, water does not usually flow out of the pump while the pump is reversed. In other words, the operating condition trajectory of the pump generally does not enter the reverse dissipation zone (shadow area in Figure 3), so the data for the first to third quadrants in the CPCs can meet the calculation requirements of the hydraulic transients [36,37].
according to the positive and negative relationship between n and q, the entire oper conditions of the pump can be divided into the turbine zone (q ≤ 0, n < 0, first quad dissipation zone (q < 0, n ≥ 0, second quadrant), normal pump zone (q > 0, n ≥ 0, quadrant) and reverse dissipation zone (q > 0, n < 0, fourth quadrant) [35]. When the p is stopped in the case of an accident, water does not usually flow out of the pump the pump is reversed. In other words, the operating condition trajectory of the pump erally does not enter the reverse dissipation zone (shadow area in Figure 3), so the for the first to third quadrants in the CPCs can meet the calculation requirements hydraulic transients [36,37]. To use the CPCs for the numerical simulation, the curves in Figure 2 will be s as data arrays (xi, WH(xi), WB(xi)) in a computer. Based on Equation (13), the value can be determined as follows: Transforming the Equation (11)  To use the CPCs for the numerical simulation, the curves in Figure 2 will be stored as data arrays (x i , WH(x i ), WB(x i )) in a computer. Based on Equation (13), the value of q/n can be determined as follows: Transforming the Equation (11) yields: Combining Equations (14)- (16), and according to the positive and negative relationship of n and q corresponding to the pump operation zones where x i is located, the values of n/ √ h and q/ √ h corresponding to any x i can be given as follows: where x 1 is the x i value corresponding to WH > 0 in the range of π 2 , 3π 2 ].
Based on Equations (11) and (12), m/h can be expressed as: When the CPCs data after the Suter transformation are given, the values of q/ √ h, n/ √ h and m/h can be determined through Equations (17) to (19). The basic principle of mathematical transformation and inverse transformation of the CPCs is shown in Figure 4.
Based on Equations (11) and (12), ℎ ⁄ can be expressed as: When the CPCs data after the Suter transformation are given, the values of √ℎ ⁄ , √ℎ ⁄ and ℎ ⁄ can be determined through Equations (17) to (19). The basic principle of mathematical transformation and inverse transformation of the CPCs is shown in Figure  4.

COPs Analysis
The generic q/ √ h − n/ √ h and m/h − n/ √ h curves of the pump are presented in Figure 5, showing the distribution of the COPs on the curves, in which four operating regimes are distinguished: turbine, dissipation, normal pump and reverse dissipation [38]. In Figure 5, point A represents the demarcation point of the turbine region and reverse dissipation region, while point C is the demarcation point of the pump region and dissipation region. Their discharges are equal to zero. Point B is the demarcation point of the turbine region and dissipation region, and its speed is equal to zero. Point M represents the runaway point of the pump, whose torque is equal to zero. Point O represents the rated working point of the pump, i.e., the best efficiency point (BEP), and its coordinate value is always (1,1). Moreover, to improve the accuracy of the prediction curves, the COPs would be increased to ensure that there are at least three COPs in each quadrant. Therefore, it is necessary to add some other COPs, e.g., point D in the pump region, which , and point P in the dissipation region, whose included angle with the positive horizontal axis is always 135 degrees, i.e., exists at ( √ℎ ⁄ ) = [( √ℎ ⁄ ) + ( √ℎ ⁄ ) ] /2 = 0.5, and point P in the dissipation region, whose included angle with the positive horizontal axis is always 135 degrees, i.e., ( √ℎ ⁄ ) = −( √ℎ ⁄ ) .
The position of the above COPs has a significant impact on the shape of the √ℎ ⁄ − √ℎ ⁄ and ℎ ⁄ − √ℎ ⁄ curves. Indeed, as long as the number of COPs is sufficient and the position coordinates are accurate, the √ℎ ⁄ − √ℎ ⁄ and ℎ ⁄ − √ℎ ⁄ curves of the pump can be determined by combining Equations (9) and (10). On this basis, the CPCs in the form of the Suter transformation (x-WH and x-WB curves) can be drawn by the inverse transformation.
(a) (b) Figure 5. The distribution of the COPs on the characteristic curves after mathematical transformation: Since the locations of these COPs are related mainly to the specific speed of the pump, for constructing the regression model between the characteristic parameters of each COP and specific speed, i.e., √ℎ ⁄ = ( ), √ℎ ⁄ = ( ) and ℎ ⁄ = ( ), the measured characteristic data of eight centrifugal pumps were collected in Table 1, and the specific speeds of these pumps range from 20 to 90. The specific speed of the pump is defined as: where ns represents the specific speed. For a two-stage pump, the specific speed of the pump shall be calculated with H/2 instead of H in Equation (20). For a double-suction pump, the specific speed shall be calculated with Q/2 instead of Q [39].
The position of the above COPs has a significant impact on the shape of the q/ h curves. Indeed, as long as the number of COPs is sufficient and the position coordinates are accurate, the q/ √ h − n/ √ h and m/h − n/ √ h curves of the pump can be determined by combining Equations (9) and (10). On this basis, the CPCs in the form of the Suter transformation (x-WH and x-WB curves) can be drawn by the inverse transformation.
Since the locations of these COPs are related mainly to the specific speed of the pump, for constructing the regression model between the characteristic parameters of each COP and specific speed, i.e., q/ √ h = f (n s ), n/ √ h = f (n s ) and m/h = f (n s ), the measured characteristic data of eight centrifugal pumps were collected in Table 1, and the specific speeds of these pumps range from 20 to 90. The specific speed of the pump is defined as: where n s represents the specific speed. For a two-stage pump, the specific speed of the pump shall be calculated with H/2 instead of H in Equation (20). For a double-suction pump, the specific speed shall be calculated with Q/2 instead of Q [39]. Since the dimensionless characteristic parameters are nonlinear functions of n s , the relationships can be described simply by the regression Model Y = ax b based on the distribution of data points. To evaluate the fitting quality, the correlation coefficient is defined as follows [45]: where y and y are the original data and the fitting results, respectively; and y and y are the mean values of y and y , respectively.

Characteristic Parameters of Points A and C
The q/ √ h values of both points A and C are zero (i.e., the pump discharge is zero).
According to Equation (9), n/ where ξ 1 is related to the diameter of the impeller outlet, the rated speed and rated head of the pump, i.e., related to the specific speed at the rated operating condition.  Figure 6c, with the correlation coefficient R close to 0.87. The n/ √ h value of point A varies in the range of −1.43 to −1.16 with relatively large variation; therefore, the correlation coefficient R is relatively low and does not exceed 0.60, as shown in Figure 6a. In addition, both dm/dn > 0 and dm/dn < 0 exist near point C while dm/dn > 0 near point A, when the pump is operating near point C, a small oscillation of pump speed will also cause a large change in the pump power and torque, resulting in the correlation coefficient between m/h and n s of point C lower than that of point A, as shown in Figure 6b,d. The relationships between n s and the characteristic parameters of points A and C are as follows:

Characteristic Parameters of Points B and M
The n/ √ h value of point B is zero (i.e., the speed of the pump is zero). At this point, the pump is about to reverse, and small changes in the discharge and shaft torque will also cause large changes in speed. Figure 7 shows the data distribution of the characteristic parameters (q/  Figure 7a,b, and the correlation coefficients R are 0.9454 and 0.7078, respectively. Point M is the runaway point, and its m/h value is always zero. At this point, the reverse rotational speed of the pump is the largest, and the centrifugal effect of the water flow through the impeller is significant, so the small change in the reverse rotational speed will also lead to a large change in the reverse flow. From Figure 7c, the n/ √ h value of point M presents less variation, ranging from −1.24 to −1.08, with a small variation. Therefore, although the correlation coefficient R is less than 0.62, power function fitting can still be used to approximately determine the functional relationship between  Figure 7d, and the

Characteristic Parameters of Points B and M
The √ℎ ⁄ value of point B is zero (i.e., the speed of the pump is zero). At this point, the pump is about to reverse, and small changes in the discharge and shaft torque will also cause large changes in speed. Figure 7 shows the data distribution of the characteristic parameters ( √ℎ ⁄ , ℎ ⁄ ) of point B. The ( √ℎ ⁄ ) value tends to increase gradually with the growth of ns, while the ( ℎ ⁄ ) value tends to decrease approximately with the growth of ns. The relationships between ( √ℎ ⁄ ) and ns, ( ℎ ⁄ ) and ns obtained by the power function fitting are shown in Figure 7a,b, and the correlation coefficients R are 0.9454 and 0.7078, respectively. Point M is the runaway point, and its ℎ ⁄ value is al-

Characteristic Parameters of Points D and P
Points D and P are set as 'custom' COPs to satisfy at least three COPs per quadrant, and the statistical data of their corresponding characteristic parameters are shown in Figure 8.

Model Verification
To verify the mathematical model proposed in Section 2, the CPC data (x, WH, WB) with ns = 24.6 in the reference [46] are selected as the validation object. The characteristic parameters of COPs on the √ℎ ⁄ − √ℎ ⁄ and ℎ ⁄ − √ℎ ⁄ curves corresponding to ns = 24.6 are obtained by Equations (14) to (19), and their values are shown in Table 2. Then, the characteristic parameters of COPs are substituted into Equation (9), and the values of undetermined coefficients (ξ1, ξ2, and ξ3) in the corresponding quadrant are solved by successively selecting each adjacent three COPs. For example, using the adjacent three points A, M, and B to draw the A-M segment. In this way, the transformed √ℎ ⁄ − √ℎ ⁄ curve with ns = 24.6 is obtained. Similarly, substituting the characteristic parameters into Equation (10), the values of undetermined coefficients (λ1, λ2 and λ3) are determined, and

Model Verification
To verify the mathematical model proposed in Section 2, the CPC data (x, WH, WB) with n s = 24.6 in the reference [46] are selected as the validation object. The characteristic parameters of COPs on the q/ √ h − n/ √ h and m/h − n/ √ h curves corresponding to n s = 24.6 are obtained by Equations (14) to (19), and their values are shown in Table 2. Then, the characteristic parameters of COPs are substituted into Equation (9), and the values of undetermined coefficients (ξ 1 , ξ 2 , and ξ 3 ) in the corresponding quadrant are solved by successively selecting each adjacent three COPs. For example, using the adjacent three points A, M, and B to draw the A-M segment. In this way, the transformed q/ √ h − n/ √ h curve with n s = 24.6 is obtained. Similarly, substituting the characteristic parameters into Equation (10), the values of undetermined coefficients (λ 1 , λ 2 and λ 3 ) are determined, and then the transformed m/h − n/ √ h curve with n s = 24.6 can be easily obtained. The q/ √ h − n/ √ h and m/h − n/ √ h curves theoretically constructed from the mathematical model in Section 2 are compared with the q/ √ h − n/ √ h and m/h − n/ √ h curves obtained from the measured CPCs after mathematical transformation in Section 3.1, as shown in Figure 9. (a) (b) Figure 9. Comparisons of the measured curve and verified curve:

Prediction Comparison
When the rated parameters of the centrifugal pump (including the rated head, rated m h Figure 9. Comparisons of the measured curve and verified curve: As expected, the theoretically constructed q/ √ h − n/ √ h and m/h − n/ √ h curves are highly consistent with the measured q/ √ h − n/ √ h and m/h − n/ √ h curves in each quadrant, also proving that the mathematical model derived in Section 2 to describe the complete characteristic curves of the centrifugal pump is reliable and valid. As long as the characteristic parameters of COPs on the q/ √ h − n/ √ h and m/h − n/ √ h curves are yielded, the CPCs for an arbitrary specific speed can be accurately predicted by combining the established mathematical model. In other words, the engineer needs to test only each characteristic operating condition when conducting the costly and time-consuming model test for CPCs, instead of tediously performing tests on all complex operating conditions such as turbine zone, dissipation zone and reverse dissipation zone, so as to save costs and accelerate the engineering design.

Prediction Comparison
When the rated parameters of the centrifugal pump (including the rated head, rated discharge and rated speed) are given, the specific speed n s can be yielded by Equation (20).
Then, according to the relationship model between the characteristic parameters of COPs and n s established in Section 3.2, the characteristic parameters of COPs can be determined. Substituting these characteristic parameters into the mathematical model in Section 2 (Equations (9) and (10)) and obtaining each of the coefficients to be determined in the model, the q/ √ h − n/ √ h and m/h − n/ √ h curve data are available. Ultimately, the q/ √ h − n/ √ h and m/h − n/ √ h curves can be inverted using the inverse transformation shown in Figure 4 to obtain the CPCs (WH and WB curves). In this paper, n s = 24.6 is selected for prediction and is compared with the measured CPCs. The results are shown in Figures 10 and 11. The CPCs predicted by the theoretical construction are basically in good agreement with the CPCs measured by the pump impeller model test. Nevertheless, due to the complex internal flow regime when the pump operates in the dissipation zone, the statistics of the characteristic parameters of COPs such as points B and C are prone to deviation. However, these deviations are not large. The maximum absolute deviations between the predicted values of WH and WB and their measured values are only 0.06 and 0.15, respectively. Indeed, the accuracy of the CPCs constructed by the prediction method proposed in this paper depends on the statistical accuracy of the characteristic parameters corresponding to the COPs. The higher the statistical accuracy of the characteristic parameters, the higher the accuracy of the CPCs constructed theoretically.

Case Study and Discussion
A water delivery pump station in northern China is equipped with four single-stage double-suction centrifugal pumps of the same type, with three working pumps and one

Case Study and Discussion
A water delivery pump station in northern China is equipped with four single-stage double-suction centrifugal pumps of the same type, with three working pumps and one reserve pump, and its layout is shown in Figure 12. The reserve pump unit is considered Figure 11. Comparisons of measured complete characteristic curves of pump and predicted complete characteristic data of pump (0 ≤ x ≤ 5π/4): (a) comparison of predicted WH curve and measured WH data; (b) comparison of predicted WB curve and measured WB data.

Case Study and Discussion
A water delivery pump station in northern China is equipped with four single-stage double-suction centrifugal pumps of the same type, with three working pumps and one reserve pump, and its layout is shown in Figure 12. The reserve pump unit is considered activated with the other three working pumps when the water intake increases the discharge. The rated parameters of a single pump are Q r = 1.05 m 3 /s, H r = 130 m, N r = 993 rpm, η r = 84.5%, and the specific speed of the pump is n s = 18.7, according to Equation (20). The combined pump-motor moment of inertia is 188.75 kg·m 2 . The installation position of the pump is close to the intake pool. Four pumps share one outlet pipe with a total length of 366 m, a diameter of 1400 mm, and an average roughness of 0.014. The outlet control valve adopts type HQ7x40H-10 slow-closing butterfly valve with a diameter of 1000 mm. The water levels in the upstream intake pool and downstream outlet pool are 1264.65 m and 1395.00 m, respectively. Since measured CPCs with ns = 18.7 were lacking, to quantitatively evaluate and verify the reliability of the CPCs prediction method proposed in this paper, numerical simulations of the hydraulic transient process are carried out under the same accidental pumpfailure transient conditions using the following CPCs: the predicted CPCs with ns = 18.7 and ns = 24.6 obtained by the prediction method in this paper, the measured CPCs with ns = 24.6, the CPCs with ns = 24.6 obtained by linear interpolation of ns = 21.1 and ns = 29.0 measured data and the measured CPCs with ns = 25.6. Among the transient conditions, the following are included: condition one, in which three working pumps suddenly encounter an accidental power failure and the control valve at the pump outlet fails to close; and condition two, in which three working pumps suddenly encounter an accidental power failure and the control valve at the pump outlet is subjected to a optimized twostage closing strategy (3 s fast shutdown to 20% opening and 10 s to completely closed). The simulation results of condition one and condition two are shown in Figure 13. The maximum and minimum pressures at the outlet of the control valve and along the pipeline and the maximum relative reverse rotational speed of the pump under these two transient conditions are collected in Table 3. Since measured CPCs with n s = 18.7 were lacking, to quantitatively evaluate and verify the reliability of the CPCs prediction method proposed in this paper, numerical simulations of the hydraulic transient process are carried out under the same accidental pump-failure transient conditions using the following CPCs: the predicted CPCs with n s = 18.7 and n s = 24.6 obtained by the prediction method in this paper, the measured CPCs with n s = 24.6, the CPCs with n s = 24.6 obtained by linear interpolation of n s = 21.1 and n s = 29.0 measured data and the measured CPCs with n s = 25.6. Among the transient conditions, the following are included: condition one, in which three working pumps suddenly encounter an accidental power failure and the control valve at the pump outlet fails to close; and condition two, in which three working pumps suddenly encounter an accidental power failure and the control valve at the pump outlet is subjected to a optimized two-stage closing strategy (3 s fast shutdown to 20% opening and 10 s to completely closed). The simulation results of condition one and condition two are shown in Figure 13. The maximum and minimum pressures at the outlet of the control valve and along the pipeline and the maximum relative reverse rotational speed of the pump under these two transient conditions are collected in Table 3. and condition two, in which three working pumps suddenly encounter an accidental power failure and the control valve at the pump outlet is subjected to a optimized twostage closing strategy (3 s fast shutdown to 20% opening and 10 s to completely closed). The simulation results of condition one and condition two are shown in Figure 13. The maximum and minimum pressures at the outlet of the control valve and along the pipeline and the maximum relative reverse rotational speed of the pump under these two transient conditions are collected in Table 3.    From Figure 13 and Table 3, for the condition one, although the specific speeds of n s = 24.6 and n s = 25.6 are close, the simulation results obtained by using the measured CPCs corresponding to the two specific speeds are quite different. The maximum pressure difference at the control valve outlet reaches 14.14 m, and the minimum pressure difference is close to 4 m. The relative deviation of the minimum pressure along the pipeline is more than 11%. The maximum pressure at the control valve outlet obtained by using the linear interpolated CPCs with n s = 24.6 is 22.32 m higher than that obtained by using the measured CPCs with n s = 24.6, and the relative deviation is close to 13%. Moreover, the fluctuation amplitude of the pump head and discharge simulated by interpolated CPCs is larger, which is not conducive to the safe operation of the project. In addition, different from the others, the simulation results obtained by using the predicted CPCs with n s = 24.6 are in good agreement with the simulation results obtained by using the measured CPCs with n s = 24.6, in which the relative deviations of the maximum and minimum pressures at the outlet of the control valve are 0.9% and 3.4%, respectively, and the relative deviation of the maximum reverse rotational speed of the pump is 3.0%. The relative deviation of the minimum pressure along the pipeline is 9.8%, but if the relative deviation of 9.8% is converted into the absolute value, it is only 0.74 m, which is considered acceptable. When the predicted CPCs with n s = 18.7 are used, the maximum and minimum pressures along the pipeline are 180.87 m and 6.66 m, respectively. The maximum relative reverse rotational speed is −1.29, which exceeds 1.2 times of the rated speed of the pump and does not meet the specification requirements.
Similar to condition one, the maximum and minimum pressures at the outlet of the control valve simulated by using the measured CPCs with n s = 24.6 and n s = 25.6 differ greatly, and the differences are 6.11 m and 3.91 m respectively. The relative deviation of the minimum pressure along the pipeline is close to 16%. The maximum pressure at the outlet of the control valve simulated by the interpolated CPCs with n s = 24.6 is 11.35 m higher than that simulated by the measured CPCs with n s = 24.6. However, compared with condition one, since the two-stage closing strategy can alleviate the fluctuation amplitude of the pump head and discharge, the variation process of pump head and discharge simulated by the interpolated CPCs with n s = 24.6 under condition two is basically consistent with the simulation results of the measured CPCs with n s = 24.6. Moreover, for the condition two, the simulation results obtained by using the predicted CPCs and measured CPCs with n s = 24.6 are also in good agreement, in which the relative deviations of the maximum and minimum pressures at the control valve outlet do not exceed 3.4%, and the relative deviation of the maximum reverse rotational speed of the pump does not exceed 2.7%. The deviation of maximum and minimum pressures along the pipeline does not exceed 0.83 m. In addition, the maximum and minimum pressures along the pipeline simulated by using the predicted CPCs with n s = 18.7 under condition two are smaller than that under condition one. The maximum pressure along the pipeline shows the largest difference, which is 11.54 m. The maximum relative reverse speed is −1.10, which is less than 1.2 times that of the rated speed of the pump, and meets the specification requirements.
The above analysis shows the following. (1) For some projects, when the measured CPCs are not available, applying the CPCs with an approximate specific speed to simulate the hydraulic transients of the pump station may cause large errors, introducing adverse effects to the safe operation of the project. (2) In addition to applying the CPCs with an approximate specific speed, linear interpolation is often used to obtain the complete characteristic curves of pumps. However, this approach may also produce large errors in the simulation results, causing safety hazards. Furthermore, when the actual specific speed is not within the range of collected pump specific speeds, the corresponding CPCs data cannot be obtained by linear interpolation. (3) The CPCs prediction method proposed in this paper can simulate the common transient processes of pump stations. The simulation results are basically consistent with the simulation results using the measured CPCs, and the error is smaller than the traditional method of using linear interpolation and using the CPCs with an approximate specific speed. It can fully meet the accuracy requirements for the hydraulic transient simulation in the preliminary design of the pump stations, and provide important data support for the design and safe operation of the project.

Conclusions
Considering that the traditional prediction method for the CPCs is relatively low in terms of prediction accuracy and ignores the actual operation characteristics of the pump, a novel method considering the inherent operating characteristics of the centrifugal pump is proposed to predict the CPCs, which serves for the design of pump stations. First, a mathematical model describing the complete characteristics of a centrifugal pump is constructed based on the Euler equation and the velocity triangles at the pump impeller. Then, the nonlinear functional relationship between the characteristic parameters corresponding to the COPs on the q/ √ h − n/ √ h and m/h − n/ √ h curves and the specific speed are established based on multiple measured CPCs at various specific speeds. Ultimately, through model verification and a case study, the following conclusions are drawn: (1) Under the condition that the characteristic parameters corresponding to all COPs on the q/ √ h − n/ √ h and m/h − n/ √ h curves are known, the CPCs constructed by the mathematical model derived in this paper are in good agreement with the measured CPCs. This proves that the prediction model of the CPCs proposed in this paper is effective. When the characteristic parameters are unknown, by substituting the regression model of the characteristic parameters into the mathematical model describing the CPCs, the CPCs for a given specific speed can be predicted successfully, and the error of WH and WB values is within the acceptable range. Theoretically, the higher the statistical accuracy of characteristic parameters of COPs, the higher the prediction accuracy of the method proposed in this paper.
(2) The CPCs prediction method proposed in this paper can simulate the common transient processes of pump stations, and the error is smaller than the traditional method of using linear interpolation and using the CPCs with an approximate specific speed. The simulation accuracy of the prediction method can meet the requirements for the hydraulic transient simulation in the preliminary design, and provide important data support for the design and safe operation of the project.