Next Article in Journal
Effect of Seawater and Fly Ash Contaminants on Insulator Surfaces Made of Polymer Based on Finite Element Method
Previous Article in Journal
Suppression Method of Partial Discharge Interferences Based on Singular Value Decomposition and Improved Empirical Mode Decomposition
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
School of Civil Engineering and Architecture, Nanchang University, Nanchang 330031, China
2
State Key Laboratory of Simulation and Regulation of Water Cycle in River Basin, China Institute of Water Resources and Hydropower Research, Beijing 100038, China
3
Department of Hydraulic Engineering, Jiangxi Water Resources Institute, Nanchang 330013, China
4
State Key Laboratory of Hydroscience and Engineering, Tsinghua University, Beijing 100084, China
5
Shanxi Dongshan Water Supply Project Construction Management Co., Ltd., Taiyuan 030000, China
*
Author to whom correspondence should be addressed.
Energies 2021, 14(24), 8580; https://doi.org/10.3390/en14248580
Submission received: 22 November 2021 / Revised: 14 December 2021 / Accepted: 16 December 2021 / Published: 20 December 2021

Abstract

:
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 difficult 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 specific speed is established. Finally, by combining the mathematical model with the nonlinear relationship, the CPCs for a given specific 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.

1. 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 three-dimensional 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.

2. 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]:
[ M H e ] = [ M H η ] = [ ρ Q ω g ] [ V 2 u r 2 V 1 u r 1 ]
where subscripts 1 and 2 represent the impeller inlet and outlet, respectively; M is the shaft torque (N·m); He is the theoretical head (m); H is the practical head (m); η = H/(H + KQ2) is the hydraulic efficiency; K is the hydraulic loss coefficient; Q is the discharge (m3/s); ω is the angular velocity of rotation (rad/s); Vu is the peripheral component of the absolute velocity (m/s); r is the impeller radius (m); ρ is the flow density (kg/m3); and g is the gravitational acceleration (m/s2).
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]:
Q = V 1 m A 1 = 2 π b 1 r 1 V 1 m
where V1m is the meridional flow velocity at the impeller inlet (m/s), A1 is the section area at the impeller inlet (m2), and b1 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:
Q = V 2 m A 2 = 2 π b 2 r 2 V 2 m
where V2m is the meridional flow velocity at the impeller outlet (m/s), A2 is the section area at the impeller outlet (m2), 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]:
V 1 u = V 1 m cot α 1 = Q cot α 1 2 π b 1 r 1
V 2 u = U 2 V 2 m cot β 2 = ω r 2 Q cot β 2 2 π b 2 r 2
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 (°).
Substituting Equations (4) and (5) into Equation (1) yields:
[ M H e ] = [ M H η ] = [ ρ Q ω g ] [ ω r 2 2 ( cot α 1 2 π b 1 + cot β 2 2 π b 2 ) Q ]
According to the similarity principle, the following dimensionless similarity relationship exists when a pump has similar working conditions [34]:
q n = C o n s t .       n h = C o n s t .       q h = C o n s t .     m h = C o n s t .
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.
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/Qr, n = N/Nr, h = H/Hr, m = M/Mr, ω = 2πN/60, η = H/(H + KQ2) into Equation (6) and eliminating Q, ω, H and M, the following equations can be deduced:
m h = λ 1 q h n h + λ 2 ( q h ) 2
ξ 1 ( n h ) 2 + ξ 2 q h n h + ξ 3 ( q h ) 2 = 1
where λ 1 = π ρ r 2 2 Q r N r / 30 M r , λ 2 = ρ ( c o t α 1 / 2 π b 1 + c o t β 2 / 2 π b 2 ) Q r 2 / M r , ξ 1 = π 2 r 2 2 N r 2 / ( 900 g H r ) , ξ 2 = π ( c o t α 1 / 2 π b 1 + c o t β 2 / 2 π b 2 ) Q r N r / ( 30 g H r ) , and ξ 3 = K Q r 2 / 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:
m h = λ 1 q h n h + λ 2 ( q h ) 2 + λ 3
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.

3. Suter Transformation and COP Analysis of Complete Characteristic Curves

3.1. 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:
W H ( x ) = h q 2 + n 2
W B ( x ) = m q 2 + n 2
x = { arctan ( q n ) ,     n < 0 ,   q 0 π + arctan ( q n ) , n 0 2 π + arctan ( q n ) , n < 0 ,   q 0
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].
To use the CPCs for the numerical simulation, the curves in Figure 2 will be stored as data arrays (xi, WH(xi), WB(xi)) in a computer. Based on Equation (13), the value of q/n can be determined as follows:
q n ( x i ) = { t a n x i , x i [ 0 , π 2 ) tan ( x i π ) , x i ( π 2 , 3 π 2 )
Transforming the Equation (11) yields:
n 2 h ( x i ) = 1 W H ( x i ) [ 1 + ( q n ( x i ) ) 2 ]
q 2 h ( x i ) = n 2 h ( x i ) ( q n ( x i ) ) 2
Combining Equations (14)–(16), and according to the positive and negative relationship of n and q corresponding to the pump operation zones where xi is located, the values of n / h and q / h corresponding to any xi can be given as follows:
n h ( x i ) = { 1 W H ( x i ) [ [ 1 + ( q n ( x i ) ) ] 2 ,   x i [ 0 , π 2 ) 0 ,   x i = π 2 1 W H ( x i ) [ [ 1 + ( q n ( x i ) ) ] 2 , x i ( π 2 , x 1 ]
where x1 is the xi value corresponding to WH > 0 in the range of ( π 2 , 3 π 2 ] .
q h ( x i ) = { n h ( x i ) q n ( x i ) , x i [ 0 , π 2 ) ( π 2 , x 1 ] 1 W H ( x i ) , x i = π 2
Based on Equations (11) and (12), m / h can be expressed as:
m h ( x i ) = W B ( x i ) W H ( x i )
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.

3.2. 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 exists at ( q / h ) D = [ ( q / h ) C + ( q / h ) O ] / 2 = 0.5 , and point P in the dissipation region, whose included angle with the positive horizontal axis is always 135 degrees, i.e., ( n / h ) P = ( q / h ) P .
The position of the above COPs has a significant impact on the shape of the q / h n / h and m / h n / 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:
n s = N Q H 3 / 4
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].
Table 1. Eight centrifugal pumps with different specific speeds.
Table 1. Eight centrifugal pumps with different specific speeds.
Serial NumberSpecific SpeedData Source
US UnitsMetric Units (Used in This Paper)
1103220.00Ayder [40]
2108921.10Chen [41]
3132025.60Brown [42]
4149629.00Kittredge [43]
5155530.14Chen [41]
6193537.50Kittredge [43]
7367671.23Liu [44]
8419981.37Chen [41]
Since the dimensionless characteristic parameters are nonlinear functions of ns, the relationships can be described simply by the regression Model Y = axb based on the distribution of data points. To evaluate the fitting quality, the correlation coefficient is defined as follows [45]:
R = i = 1 n ( y i y ¯ ) ( y i y ¯ ) i = 1 n ( y i y ¯ ) 2 i = 1 n ( y i y ¯ ) 2
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.

3.2.1. 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 / h ) C = 1 / ξ 1 , 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 6 statistically analyzes the characteristic parameters ( n / h and m / h values) of points A and C according to eight sets of centrifugal pumps with different specific speeds in Table 1. The q / h value of point C varies within the range of 0.79 to 0.96 with a small variation, and the relationship between ( n / h ) C and ns obtained by fitting the power function Y = axb is shown in 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 d m / d n > 0 and d m / d n < 0 exist near point C while d m / d n > 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 ns of point C lower than that of point A, as shown in Figure 6b,d. The relationships between ns and the characteristic parameters of points A and C are as follows:
{ ( n h ) A = 0.9107 n s 0.1058 q h A = 0 ( m h ) A = 0.0719 n s 0.8316
{ ( n h ) C = 1.2759 n s 0.1101 q h C = 0 ( m h ) C = 0.0661 n s 0.4619

3.2.2. 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 / h ,   m / h ) of point B. The ( q / h ) B value tends to increase gradually with the growth of ns, while the ( m / h ) B value tends to decrease approximately with the growth of ns. The relationships between ( q / h ) B and ns,   ( m / h ) B 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 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 ( n / h ) M and ns. The q / h value of point M varies from −0.89 to −0.32 with relatively large variation, and the q / h value shows a gradual decreasing trend with the growth of ns. The fitted relationship between ( q / h ) M and ns is shown in Figure 7d, and the correlation coefficient R is 0.8655.
{ ( n h ) B = 0 q h B = 3.3672 n s 0.3084 ( m h ) B = 2.6528 n s 0.2148
{ ( n h ) M = 0.9500 n s 0.0581 q h M = 0.2339 n s 0.2996 ( m h ) M = 0

3.2.3. 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. The n / h value of point D is distributed between 0.87 and 0.96, and the range of variation is small. The relationship between ( n / h ) D and ns obtained by the power function fitting is shown in Figure 8a, and the correlation coefficient exceeds 0.86. Additionally, from the statistics of ( n / h ) P n s and ( q / h ) P n s shown in Figure 8c,d, respectively, the data points are concentrated above and below the regression model Y = axb, so the obtained correlation coefficient R is also high, close to 0.87. Figure 8b,e show that the m / h values of points D and P vary greatly, and the distribution of data points is relatively discrete. The correlation coefficients R obtained by the power function fitting are 0.8566 and 0.7622, respectively.
{ ( n h ) D = 1.1261 n s 0.0598 q h D = 0.5 ( m h ) D = 0.2335 n s 0.2582
{ ( n h ) P = 1.6204 n s 0.2519 ( q h ) P = 1.6204 n s 0.2519 ( m h ) P = 0.3115 n s 0.1962

4. Model Verification and CPCs Prediction

4.1. 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 q / h n / h and m / h n / h 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 q / h n / h 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 then the transformed m / h n / h curve with ns = 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.
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.

4.2. Prediction Comparison

When the rated parameters of the centrifugal pump (including the rated head, rated discharge and rated speed) are given, the specific speed ns can be yielded by Equation (20). Then, according to the relationship model between the characteristic parameters of COPs and ns 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, ns = 24.6 is selected for prediction and is compared with the measured CPCs. The results are shown in Figure 10 and Figure 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.

5. 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 Qr = 1.05 m3/s, Hr = 130 m, Nr = 993 rpm, ηr = 84.5%, and the specific speed of the pump is ns = 18.7, according to Equation (20). The combined pump-motor moment of inertia is 188.75 kg·m2. 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 pump-failure 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 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.
From Figure 13 and Table 3, for the condition one, although the specific speeds of ns = 24.6 and ns = 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 ns = 24.6 is 22.32 m higher than that obtained by using the measured CPCs with ns = 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 ns = 24.6 are in good agreement with the simulation results obtained by using the measured CPCs with ns = 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 ns = 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 ns = 24.6 and ns = 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 ns = 24.6 is 11.35 m higher than that simulated by the measured CPCs with ns = 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 ns = 24.6 under condition two is basically consistent with the simulation results of the measured CPCs with ns = 24.6. Moreover, for the condition two, the simulation results obtained by using the predicted CPCs and measured CPCs with ns = 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 ns = 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.

6. 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.

Author Contributions

Conceptualization, W.H. and H.L. (Huokun Li); methodology, H.L. (Huokun Li), W.H. and H.L. (Hongkang Lin); formal analysis, H.L. (Hongkang Lin), H.L. (Huokun Li) and J.L.; data curation, H.L. (Hongkang Lin), M.Z. and X.H.; writing original draft preparation, H.L. (Hongkang Lin); writing review and editing, W.H., H.L. (Huokun Li) and J.M.; supervision, J.M. and J.L.; project administration, W.H.; funding acquisition, W.H., M.Z. and X.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. 51909115, 52079061), Natural Science Foundation of Jiangxi, China (Grant No. 20192BAB216038), Open Research Fund of State Key Laboratory of Simulation and Regulation of Water Cycle in River Basin, China Institute of Water Resources and Hydropower Research (Grant No. IWHR-SKL-KF201904), Water Science and Technology Projects of Water Resources Department of Jiangxi Province (Grant No. 202123YBKT04) and 2021 Annual Scientific Research Program of Wanjiazhai Water Holding Group Co., Ltd. (Grant No. 2021−30).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

CPCcomplete pump characteristic
COPcharacteristic operating point
Mshaft torque, N∙m
Hetheoretical head, m
Hpractical head, m
Qdischarge, m3/s
Khydraulic loss coefficient
ωangular velocity of rotation, rad/s
Vuperipheral components of the absolute velocity, m/s
rimpeller radius, m
ρflow density, kg/m3
ggravitational acceleration, m/s2
Vmmeridional flow velocity, m/s
Asection area, m2
bsection width, m
Uperipheral velocity, m/s
α 1 absolute flow angle, °
β 2 relative flow angle, °
qrelative discharge
hrelative head
nrelative rotational speed
mrelative torque
Nrotational speed, rpm
nsspecific speed
WHhead characteristics of pump
WBtorque characteristics of pump
Subscripts
rrated condition
max.maximum value
min.maximum value
1impeller inlet
2impeller outlet

References

  1. Wang, F.J.; Tang, X.L.; Chen, X.; Xiao, R.F.; Yao, Z.F.; Yang, W. A review on flow analysis method for pumping stations. J. Hydraul. Eng. 2018, 49, 47–61. [Google Scholar]
  2. Dai, C.; Dong, L.; Lin, H.; Zhao, F. A Hydraulic Performance Comparison of Centrifugal Pump Operating in Pump and Turbine Modes. J. Therm. Sci. 2020, 29, 1594–1605. [Google Scholar] [CrossRef]
  3. Li, D.; Fu, X.; Zuo, Z.; Wang, H.; Li, Z.; Liu, S.; Wei, X. Investigation methods for analysis of transient phenomena concerning design and operation of hydraulic-machine systems—A review. Renew. Sustain. Energy Rev. 2019, 101, 26–46. [Google Scholar] [CrossRef]
  4. Huang, W.; Yang, K.; Guo, X.; Ma, J.; Wang, J.; Li, J. Prediction Method for the Complete Characteristic Curves of a Francis Pump-Turbine. Water 2018, 10, 205. [Google Scholar] [CrossRef] [Green Version]
  5. Knapp, R.T. Complete characteristics of centrifugal pumps and their use in the prediction of transient behavior. Trans. ASME 1937, 59, 683–689. [Google Scholar]
  6. Huang, B.; Wu, P.; Liu, X.S.; Feng, X.D.; Yang, S.; Wu, D.Z. Four-quadrant Full Characteristic Model Test of the ACP100 Reactor Coolant Pump. Fluid Mach. 2020, 48, 8–11. [Google Scholar]
  7. Shukla, S.N.; Kshirsagar, J.T. Numerical Experiments on a Centrifugal Pump. In Proceedings of the ASME 2002 Joint U.S.-European Fluids Engineering Division Conference, Montreal, QC, Canada, 14–18 July 2002; pp. 709–719. [Google Scholar]
  8. Donsky, B. Complete Pump Characteristics and the Effects of Specific Speeds on Hydraulic Transients. J. Fluids Eng. 1961, 83, 685. [Google Scholar] [CrossRef]
  9. Yang, Y.S.; Dong, R.; Jing, T. Influence of Full Characteristic Curve on Pump-off Water Hammer and Its Protection. China Water Wastewater 2010, 26, 63–66. [Google Scholar]
  10. Marchal, M.; Flesch, G.; Suter, P. The calculation of waterhammer problems by means of the digital computer. In Proceedings of the International Symposium on Waterhammer in Pumped Storage Projects ASME, Chicago, IL, USA, 7–11 November 1965; pp. 168–180. [Google Scholar]
  11. Zhang, C.; Peng, T.; Zhou, J.; Ji, J.; Wang, X. An improved autoencoder and partial least squares regression-based extreme learning machine model for pump turbine characteristics. Appl. Sci. 2019, 9, 3987. [Google Scholar] [CrossRef] [Green Version]
  12. Izquierdo, J.; Montalvo, I.; Pérez-García, R.; Ayala-Cabrera, D. Multi-agent simulation of hydraulic transient equations in pressurized systems. J. Comput. Civ. Eng. 2016, 30, 04015071. [Google Scholar] [CrossRef] [Green Version]
  13. Zhu, M.L.; Zhang, X.H.; Zhang, Y.H.; Wang, T. Study on prediction model of complete characteristic curves of centrifugal pumps. J. Northwest Sci-Tech Univ. Agric. For. Nat. Sci. Ed. 2006, 4, 143–146. [Google Scholar]
  14. Hu, X.Y.; Yu, B.; Guo, J.; Wang, S.K. Visualization for predictable model of complete characteristic curve of pump. Fluid Mach. 2012, 3, 37–39. [Google Scholar]
  15. Lima, G.M.; Luvizotto Júnior, E. Method to estimate complete curves of hydraulic pumps through the polymorphism of existing curves. J. Hydraul. Eng. 2017, 143, 04017017. [Google Scholar] [CrossRef]
  16. Wan, W.; Huang, W. Investigation on complete characteristics and hydraulic transient of centrifugal pump. J. Mech. Sci. Technol. 2011, 25, 2583–2590. [Google Scholar] [CrossRef]
  17. Liu, G.L.; Jiang, J.; Fu, X.Q. Predicting complete characteristics of pumps by using BP neural network. J. Wuhan Univ. Hydraul. Electr. Eng. 2000, 2, 37–39. [Google Scholar]
  18. Lu, W.; Li, L.G. Drawing the all-character curve of water pump on BP neural network. Ind. Control. Comput. 2001, 14, 18–20. [Google Scholar]
  19. Liu, Z.F. GARBF Network Method to Predict the Complete Characteristic Curve of Pump. Master’s Thesis, Chang’an University, Xi’an, China, 2011. [Google Scholar]
  20. Han, W.; Nan, L.; Su, M.; Chen, Y.; Li, R.; Zhang, X. Research on the prediction method of centrifugal pump performance based on a double hidden layer BP neural network. Energies 2019, 12, 2709. [Google Scholar] [CrossRef] [Green Version]
  21. Huang, S.; Qiu, G.; Su, X.; Chen, J.; Zou, W. Performance prediction of a centrifugal pump as turbine using rotor-volute matching principle. Renew. Energy 2017, 108, 64–71. [Google Scholar] [CrossRef]
  22. Wu, Q.; Wang, X.; Shen, Q. Research on dynamic modeling and simulation of axial-flow pumping system based on RBF neural network. Neurocomputing 2016, 186, 200–206. [Google Scholar] [CrossRef]
  23. Höller, S.; Benigni, H.; Jaberg, H. Investigation of the 4-Quadrant behaviour of a mixed flow diffuser pump with CFD-methods and test rig evaluation. IOP Conf. Ser. Earth Environ. Sci. 2016, 49, 032018. [Google Scholar] [CrossRef]
  24. Gros, L.; Couzinet, A.; Pierrat, D.; Landry, L. Complete pump characteristics and 4-quadrant representation investigated by experimental and numerical approaches. In Proceedings of the ASME-JSME-KSME 2011 Joint Fluids Engineering Conference, Hamamatsu, Japan, 24–29 July 2011; pp. 359–368. [Google Scholar]
  25. Wang, L.; Li, M.; Wang, F.J.; Wang, J.B.; Yao, C.G.; Yu, Y.S. Study on three-dimensional internal characteristics method of Suter curves for double-suction centrifugal pump. J. Hydraul. Eng. 2017, 48, 113–122. [Google Scholar]
  26. Muttalli, R.S.; Agrawal, S.; Warudkar, H. CFD simulation of centrifugal pump impeller using ANSYS-CFX. Intern. J. Innov. Res. Sci. Eng. Technol. 2014, 3, 15553–15561. [Google Scholar] [CrossRef]
  27. Frosina, E.; Buono, D.; Senatore, A. A performance prediction method for pumps as turbines (PAT) using a computational fluid dynamics (CFD) modeling approach. Energies 2017, 10, 103. [Google Scholar] [CrossRef] [Green Version]
  28. Chang, J.S. Transients of Hydraulic Machine Installations; Higher Education Press: Beijing, China, 2005. [Google Scholar]
  29. Zhang, Y.L.; Zhu, Z.C.; Dou, H.S.; Cui, B.L. A Generalized Euler Equation to Predict Theoretical Head of Turbomachinery. Int. J. Rotating Mach. 2015, 42, 26–38. [Google Scholar] [CrossRef]
  30. El-Naggar, M.A. A one-dimensional flow analysis for the prediction of centrifugal pump performance characteristics. Int. J. Rotating Mach. 2013, 2013, 473512. [Google Scholar] [CrossRef] [Green Version]
  31. Olimstad, G.; Nielsen, T.; Børresen, B. Dependency on runner geometry for reversible-pump turbine characteristics in turbine mode of operation. J. Fluids Eng. 2012, 134, 121102. [Google Scholar] [CrossRef]
  32. Tanaka, T.; Liu, C. An Investigation on Velocity Triangle Applied to Fluid Particle at Rotating Flow Passage of Impeller Blade in Centrifugal Pump. In Proceedings of the ASME 7th Biennial Conference on Engineering Systems Design and Analysis, Manchester, UK, 19–22 July 2004; pp. 433–442. [Google Scholar]
  33. Nicolet, C. Hydroacoustic Modelling and Numerical Simulation of Unsteady Operation of Hydroelectric Systems. Ph.D. Thesis, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland, 2007. [Google Scholar]
  34. Xu, H.Q.; Lu, L.; Tan, D.Q.; Li, T.Y. Specific speed of pump and unit parameter selection. In Proceedings of the 20th China Hydropower Equipment Symposium, Chengdu, China, 15–17 October 2015. [Google Scholar]
  35. Yang, K.L. Hydraulic Transients and Regulation in Power Station and Pump Station; China Water Conservancy and Hydropower Press: Beijing, China, 2000. [Google Scholar]
  36. Matteo, L.; Cerru, F.; Dazin, A.; Tauveron, N. Investigation of the pump, dissipation and inverse turbine operating modes using the CATHARE-3 one-dimensional rotodynamic pump model. In Proceedings of the ASME-JSME-KSME 2019 8th Joint Fluids Engineering Conference, San Francisco, CA, USA, 28 July–1 August 2019. [Google Scholar]
  37. Wu, D.Z.; Xu, B.J.; Li, Z.F.; Wang, L.Q. Numerical Simulation On Internal Flow Of Centrifugal Pump During Transient Operation. J. Eng. Thermophys. 2009, 781–783. [Google Scholar]
  38. Chang, J.S. Analysis of Whole Characteristic Curve and Their Special Condition Points in Mixed-flow Pump-turbine. J. Beijing Agric. Eng. Univ. 1995, 015, 77–83. [Google Scholar]
  39. Huang, W.; Kang, Q.; Li, S.S.; Zhu, Y.X.; Yan, F.; Li, J.Z. Analysis of the effect of valve characteristics on hydraulic transition process of pumping station. South-to-North Water Transf. Water Sci. Technol. 2019, 17, 6. [Google Scholar]
  40. Ayder, E.; Ilikan, A.N.; Sen, M.; Ozgur, C.; Kavurmacıoglu, L.; Kirkkopru, K. Experimental investigation of the complete characteristics of rotodynamic pumps. In Proceedings of the ASME 2009 Fluids Engineering Division Summer Meeting, Vail, CO, USA, 2–6 August 2009; pp. 35–40. [Google Scholar]
  41. Chen, N.X. Simulation and Control of Hydraulic Transients in Hydropower Projects; China Water Power Press: Beijing, China, 2005. [Google Scholar]
  42. Brown, R.J.; Rogers, D.C. Development of Pump Characteristics from Field Tests. J. Mech. Des. 1980, 102, 807–817. [Google Scholar] [CrossRef]
  43. Kittredge, C.; Princeton, N. Hydraulic transients in centrifugal pump systems. Trans. ASME 1956, 78, 1307–1321. [Google Scholar]
  44. Liu, G.L.; Guo, M.H. Computer simulation of universal complete pump characteristics. J. Hydraul. Eng. 1998, 8, 2–8. [Google Scholar] [CrossRef]
  45. Xie, A.; Li, D.H. Probability and Statistics; Tsinghua University Press: Beijing, China, 2012. [Google Scholar]
  46. Wylie, E.B.; Streeter, V.L. Fluid Transients; McGraw-Hill International Book Co.: New York, NY, USA, 1978. [Google Scholar]
Figure 1. The velocity triangles at the impeller of a centrifugal pump.
Figure 1. The velocity triangles at the impeller of a centrifugal pump.
Energies 14 08580 g001
Figure 2. Flow charts of two methods: (a) direct method; (b) indirect method.
Figure 2. Flow charts of two methods: (a) direct method; (b) indirect method.
Energies 14 08580 g002
Figure 3. Complete characteristic curves and operation zones of the pump.
Figure 3. Complete characteristic curves and operation zones of the pump.
Energies 14 08580 g003
Figure 4. Basic principle of mathematical transformation and inverse transformation of the CPCs.
Figure 4. Basic principle of mathematical transformation and inverse transformation of the CPCs.
Energies 14 08580 g004
Figure 5. The distribution of the COPs on the characteristic curves after mathematical transformation: (a) q / h n / h curve; (b) m / h n / h curve.
Figure 5. The distribution of the COPs on the characteristic curves after mathematical transformation: (a) q / h n / h curve; (b) m / h n / h curve.
Energies 14 08580 g005
Figure 6. The relationships between ns and the characteristic parameters of point A and point C: (a) ( n / h ) A n s ; (b) ( m / h ) A n s ; (c) ( n / h ) C n s ; (d) ( m / h ) C n s .
Figure 6. The relationships between ns and the characteristic parameters of point A and point C: (a) ( n / h ) A n s ; (b) ( m / h ) A n s ; (c) ( n / h ) C n s ; (d) ( m / h ) C n s .
Energies 14 08580 g006
Figure 7. The relationships between ns and the characteristic parameters of point B and point M: (a) ( q / h ) B n s ; (b) ( m / h ) B n s ; (c) ( n / h ) M n s ; (d) ( q / h ) M n s .
Figure 7. The relationships between ns and the characteristic parameters of point B and point M: (a) ( q / h ) B n s ; (b) ( m / h ) B n s ; (c) ( n / h ) M n s ; (d) ( q / h ) M n s .
Energies 14 08580 g007
Figure 8. The relationships between ns and the characteristic parameters of point D and point P: (a)   ( n / h ) D n s ; (b) ( m / h ) D n s ; (c) ( n / h ) P n s ; (d)   ( q / h ) P n s ; (e)   ( m / h ) P n s .
Figure 8. The relationships between ns and the characteristic parameters of point D and point P: (a)   ( n / h ) D n s ; (b) ( m / h ) D n s ; (c) ( n / h ) P n s ; (d)   ( q / h ) P n s ; (e)   ( m / h ) P n s .
Energies 14 08580 g008aEnergies 14 08580 g008b
Figure 9. Comparisons of the measured curve and verified curve: (a) q / h n / h curve; (b) m / h n / h curve.
Figure 9. Comparisons of the measured curve and verified curve: (a) q / h n / h curve; (b) m / h n / h curve.
Energies 14 08580 g009
Figure 10. Comparisons of measured curve and predicted curve: (a) comparison of predicted q / h n / h curve and measured data; (b) comparison of predicted m / h n / h curve and measured data.
Figure 10. Comparisons of measured curve and predicted curve: (a) comparison of predicted q / h n / h curve and measured data; (b) comparison of predicted m / h n / h curve and measured data.
Energies 14 08580 g010
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.
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.
Energies 14 08580 g011
Figure 12. Schematic diagram of pump station layout [39].
Figure 12. Schematic diagram of pump station layout [39].
Energies 14 08580 g012
Figure 13. Simulation results of hydraulic transients under valve rejection and valve two-stage closing conditions: (a) variation process of control valve outlet pressure under condition one; (b) variation process of control valve outlet pressure under condition two; (c) variation process of relative speed, discharge and head of the pump under condition one; (d) variation process of relative speed, discharge and head of the pump under condition two; (e) maximum and minimum pressure envelopes along the pipeline under condition one; (f) maximum and minimum pressure envelopes along the pipeline under condition two.
Figure 13. Simulation results of hydraulic transients under valve rejection and valve two-stage closing conditions: (a) variation process of control valve outlet pressure under condition one; (b) variation process of control valve outlet pressure under condition two; (c) variation process of relative speed, discharge and head of the pump under condition one; (d) variation process of relative speed, discharge and head of the pump under condition two; (e) maximum and minimum pressure envelopes along the pipeline under condition one; (f) maximum and minimum pressure envelopes along the pipeline under condition two.
Energies 14 08580 g013aEnergies 14 08580 g013b
Table 2. The characteristic parameters of COPs on the q / h n / h and m / h n / h curves corresponding to ns = 24.6.
Table 2. The characteristic parameters of COPs on the q / h n / h and m / h n / h curves corresponding to ns = 24.6.
COPs q h n h m h
A0.0000−1.2559−1.0789
M−0.6526−1.15200.0000
B−1.20210.00001.2500
P−0.70850.70850.5221
C0.00000.88110.3494
D0.50000.90510.6592
O1.00001.00001.0000
Table 3. Extreme values under valve rejection and valve two-stage closing conditions.
Table 3. Extreme values under valve rejection and valve two-stage closing conditions.
Adopted CPCsCondition OneCondition Two
Max./Min. Head at Control Valve Outlet (m)Max. Relative Reverse SpeedMax./Min. Head
Along Pipeline (m)
Max./Min. Head at Control Valve Outlet (m)Max. Relative Reverse SpeedMax./Min. Head
Along Pipeline (m)
Predicted
ns = 18.7
180.87/54.42−1.29180.87/6.66169.33/54.42−1.10169.33/6.64
Measured
ns = 24.6
172.26/54.95−1.32172.26/7.54167.59/54.95−1.12167.59/7.67
Predicted
ns = 24.6
173.89/56.81−1.28173.89/6.80167.10/56.81−1.09167.10/6.84
Interpolated
ns = 24.6
194.58/56.42−1.34194.58/7.31178.94/56.43−1.07178.94/7.64
Measured
ns = 25.6
186.40/58.82−1.33186.40/6.70173.70/58.86−1.23173.70/6.45
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, H.; Lin, H.; Huang, W.; Li, J.; Zeng, M.; Ma, J.; Hu, X. A New Prediction Method for the Complete Characteristic Curves of Centrifugal Pumps. Energies 2021, 14, 8580. https://doi.org/10.3390/en14248580

AMA Style

Li H, Lin H, Huang W, Li J, Zeng M, Ma J, Hu X. A New Prediction Method for the Complete Characteristic Curves of Centrifugal Pumps. Energies. 2021; 14(24):8580. https://doi.org/10.3390/en14248580

Chicago/Turabian Style

Li, Huokun, Hongkang Lin, Wei Huang, Jiazhen Li, Min Zeng, Jiming Ma, and Xin Hu. 2021. "A New Prediction Method for the Complete Characteristic Curves of Centrifugal Pumps" Energies 14, no. 24: 8580. https://doi.org/10.3390/en14248580

APA Style

Li, H., Lin, H., Huang, W., Li, J., Zeng, M., Ma, J., & Hu, X. (2021). A New Prediction Method for the Complete Characteristic Curves of Centrifugal Pumps. Energies, 14(24), 8580. https://doi.org/10.3390/en14248580

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop