4.1. Energy Performance Analysis
Figure 7 presents the external characteristic curves for the seven design schemes.
Figure 7a shows that for a single blade deviation, the flow rate of the optimal hydraulic efficiency point of the mixed-flow pump gradually increased when blade No. 2 rotated from −4° to +4° (Design Schemes I–V). This result occurred because the flow capacity in the impeller passage generally tended toward larger flow rates with increases in the rotation angle. In all the design schemes, when the rotation angle of each blade was adjusted to 0° (Design Scheme III), the optimal hydraulic efficiency of the mixed-flow pump reached a maximum value. Therefore, it was concluded that BRADs changed the flow field and decreased the optimal hydraulic efficiency.
When the rotation angles of the other blades remained at 0° and the flow rate remained constant, the head of the mixed-flow pump gradually increased when blade No. 2 rotated counterclockwise from −4° to +4°, as shown in
Figure 7b.
Figure 8 and
Figure 9 show the pressure distribution of each span in the impeller domain under the design flow rate, which can explain the curve rules described previously. The meaning of span is as shown in Equation (4). When there were no BRADs, the pressure distributions near the three blades were relatively centrally symmetrical, and the pressure gradients were distributed evenly. As blade No. 2 gradually deviated counterclockwise, the area of low pressure on the suction surface of the impeller gradually increased, while the pressure distribution on the pressure surface changed slightly. As a result, when the rotation angle of the blade deviated in a counterclockwise direction, the pressure difference between the pressure and suction surfaces of the impeller gradually increased. Therefore, the working capability of the impeller was gradually enhanced, and the head of the mixed-flow pump increased. When two blades deviated in opposite directions simultaneously (Design Schemes VI–VII), the working capacities of the blades canceled each other out, so the head curves did not change significantly.
In Equation (4), r is the radius at the span, rhub is the hub radius of the impeller outlet, and rshroud is the shroud radius of the impeller outlet.
Figure 10 depicts comparisons of the energy performance of the mixed-flow pump for the different design schemes under the design flow rate.
Figure 10a shows that the hydraulic efficiency point for each deviation scheme was reduced from that of Design Scheme III. When the rotation angle of the blade deviated clockwise (Design Scheme I–II), the pumping capacity of the impeller decreased. The optimal hydraulic efficiency corresponded to the condition of a small flow rate, so the hydraulic efficiency decreased. When the rotation angle of the blade deviated counterclockwise (Design Scheme IV–V), the power of the impeller increased. However, hydraulic efficiency is inversely proportional to power. As the power increased, the hydraulic efficiency decreased. When the two blades had opposite deviations (Design Scheme VI–VII), the force on the impeller was uneven, which caused more hydraulic losses and reduced hydraulic efficiency. The most significant reduction occurred for Design Scheme I, which had a decrease of 2.43%, followed by that for Design Scheme VII. Greater counterclockwise and clockwise blade deviations led to lower efficiencies.
Figure 10b shows that when the blade deviation angle gradually increased from −4° to +4°, the head rose, and the gradient changed regularly. The head of Design Scheme V was 1.51 m higher than that of Design Scheme III. The head when two blades deviated in opposite directions was slightly higher than that when the blades had no deviations. Additionally, the head of Design Scheme VII was slightly lower than that of Design Scheme VI, perhaps because the fluid distribution in the pump was more uneven, thus causing more flow loss. The power was proportional to the head. When the head increased, so did the power. Therefore, the trend in
Figure 10c was essentially the same as that in
Figure 10b.
4.2. Pressure Fluctuation Characteristics
To study the variation trends of the pressure pulsation at different sections under the design flow rate, three monitoring points were set at the impeller inlet and outlet, named P11–P13 and P21–P23, respectively. Their locations are shown in
Figure 11. Pij represents the positions of the monitoring points, where i = 1–2 represents the axial position from the inlet to the outlet of the pump, and j = 1–3 represents the radial position from the hub to the shroud.
The pressure coefficient,
Cp, was selected to reflect the pressure fluctuation characteristics, and its relational expression is as follows [
26]:
In Equation (5), Pi represents the instantaneous pressure at the monitoring point, Pave is the average pressure, ρ is the density of water, and u is the circumferential velocity at the outlet of the impeller.
4.2.1. Impeller Inlet
Figure 12 shows the time-domain variations of the impeller inlet pressure fluctuations for each design scheme under the design flow rate. The last rotation period of the impeller was selected for analysis. In Design Scheme III, which had no BRADs, the impeller regularly passed through the monitoring points three times every rotation period. Therefore, the pressure distribution around the hub side, the middle position, and the shroud side presented three sets of very regular peaks and valleys, and the number of peaks was equal to the number of impeller blades. In other design schemes, although there were three peaks and valleys near the shroud side and the middle position, the periodic curves near the hub side were not obvious. These results show that when the blade rotation angle deviated, the pressure fluctuations near the hub side of the impeller inlet monitoring point were less affected by the rotation of the blade, and that the fluid flow was disordered. When there was a deviation of a single blade (Design Schemes I, II, IV, and V), the peak values for the three blades at P12 and P13 changed little compared to the changes for Design Scheme III. However, the minimum value for one set of blades decreased with increases in the blade rotation angle and increased with decreases in the blade rotation angle. This result indicates that the central symmetry of the pressure distribution at the impeller inlet was destroyed, which primarily affected the fluctuation values in the valleys. When two blades deviated, the pressure coefficient in one valley for Design Scheme VII was higher than that for Design Scheme VI, and another was lower.
Table 2 displays the pressure coefficient pulsation amplitudes for the last rotation period in each design scheme. The values for all the deviation design schemes were higher than those for the scheme without deviation, and the most obvious increase occurred for Design Scheme VII. These results indicate that when BRADs occurred, the pressure fluctuation range at the impeller inlet increased.
The calculated time-domain data were converted into frequency-domain data by a fast Fourier transform (FFT) to obtain the frequency distribution of the pressure signal. The shaft frequency was
fn =
n/60 = 24.17 Hz, and the horizontal axis of the frequency-domain diagram was set as multiples of the shaft frequency.
Figure 13 depicts the frequency-domain diagrams of the pressure pulsations at the impeller inlet for different design schemes. In all the schemes, the amplitude of the dominant pressure pulsation frequency increased gradually from the hub side to the shroud side. Unlike Design Scheme III, which had no BRADs, the dominant pressure pulsation frequency at P11 changed to the shaft frequency for the other schemes, and large secondary frequency amplitudes appeared at P12 and P13. This result may have occurred because BRADs caused the velocity distributions at the hub side to be uneven, resulting in increases in vorticity and changes in hydraulic excitation. This type of deviation had the greatest influence on Design Scheme VII because the dominant pressure pulsation frequency at the middle position became the shaft frequency, and a secondary frequency corresponding to the amplitude of the dominant frequency emerged on the side of the shroud.
Table 3 shows the amplitude of the dominant frequency at the impeller inlet for each design scheme. At P11, the amplitude changed significantly, and the amplitudes for Design Schemes I, II, IV, V, VI, and VII were 8.07, 3.61, 4.34, 7.91, 7.61, and 15.16 times that of Design Scheme III, respectively. For a single blade deviation, the amplitudes of the dominant frequency pulsations at P12 and P13 increased with a counterclockwise deviation. Deviations of two blades in opposite directions led to extremely uneven distributions of the impeller inlet pressure and significant increases in the amplitude.
4.2.2. Impeller Outlet
Figure 14 shows the time-domain variations in the impeller outlet pressure fluctuations for each design scheme. The impeller outlet results had more tortuous curves than did the results for the impeller inlet. This occurred because BRADs caused differences in speed and direction in each blade passage, leading to violent collisions between the high-speed water at the outlets. Due to the rotor-and-stator interference of the impeller and guide vane, there were three distinct peaks and valleys in the last rotation period.
Table 4 lists the pulsation amplitudes of the pressure coefficient at the impeller outlet for each design scheme. The occurrence of BRADs increased the pulsation amplitude of the pressure coefficient at the impeller outlet, and the range of the pressure change increased.
Figure 15 presents the frequency-domain diagrams of the pressure coefficient at the impeller outlet for all design schemes. The amplitude of the dominant pressure fluctuation frequency increased gradually from the hub side to the shroud side in all design schemes. When there were no BRADs, the dominant and secondary pressure pulsation frequencies were equal to the blade frequency and multiples of the blade frequency, respectively. When there were BRADs, the dominant pressure pulsation frequencies at each monitoring point for Design Schemes I, II, IV, V, and VI remained unchanged and equal to the blade frequency, but the secondary frequency changed to the shaft frequency. These results show that, despite the deviations in the blade rotation angle, the blade frequency was still dominant at the impeller outlet. However, the secondary frequency cannot be ignored, and greater BRAD led to a greater rise in the amplitude. Design Scheme VII had the most significant influence on the pressure pulsations at the impeller outlet because the dominant frequencies near the hub side and middle position were converted into the shaft frequency, while a secondary frequency with a larger amplitude appeared near the shroud side. In this case, the entire pump operation remained in a low-frequency pulsation state, which can cause resonance of the pump.
Table 5 shows that BRADs increased the amplitude of the dominant pressure pulsation frequency at each monitoring point, which was different from the situation at the inlet. These results occurred because the P2 monitoring points were located behind the impeller, and the asymmetric outflow in the three blade passages had high velocities. Therefore, the collision was greater and more severe than for Design Scheme III. The dominant pressure pulsation frequencies near the shroud side were not affected by BRADs (all of them were equal to blade frequencies). At P23, the dominant pressure pulsation frequencies of Design Schemes I, II, IV, V, VI, and VII were 1.021, 1.01, 1.037, 1.083, 1.056, and 1.061 times that of Design Scheme III, respectively.
4.3. Force on the Impeller
When the mixed-flow pump was running, the rotation of the impeller generated axial lift on the fluid, and the fluid produced an equal and opposite axial force on the impeller.
Figure 16 depicts the trajectory distribution of the axial force in the last rotation period under the condition of 1.0
Qdes. The curves are essentially closed, indicating that the axial force of the impeller reached a stable state. The fluctuation of axial force is defined as the difference between the maximum value and the minimum value of axial force in the last period. The axial force fluctuated little with the rotation of the pump, and the entire trajectory was approximately circular.
Figure 16b shows that when a BRAD occurred for a single blade, the total axial force on the impeller increased with increases in the blade rotation angle and decreased with decreases in the rotation angle. This result occurred because the axial force was determined by the pressure on the blade surface. In
Figure 9, the pressure difference on the blade surface increased with increases in the blade rotation angle, resulting in increases in the axial force, and vice versa. The fluctuation values of axial force for design schemes I, II, III, IV, V, VI, and VII are 75.05 N, 68.76 N, 16.59 N, 36.02 N, 66.68 N, 64.08 N, and 89.75 N, respectively. When two blades deviated in opposite directions, the fluctuations in the axial force increased sharply, indicating that the axial force stability was poor.
Table 6 presents the average axial force on the impeller at flow rates of 0.8
Qdes, 1.0
Qdes, and 1.2
Qdes. When the flow rate increased, the axial force decreased in each scheme. For the 1.0
Qdes flow rate, the average axial forces of Design Schemes I, II, IV, V, VI, and VII were 0.898, 0.972, 1.071, 1.112, 1.027, and 0.999 times that of Design Scheme III, respectively. When the deviation directions of two blades were opposite each other, the axial force of the impeller was essentially the same as that when there were no deviations (Design Scheme III). The results show that when two blades deviated symmetrically and reversely, the axial forces of the blades canceled each other out, leading to a situation similar to that in the scheme without deviations.
Figure 17 shows the trajectory distributions of the radial force on the impeller in the last rotation period under the condition of 1.0
Qdes. When there were no deviations, the radial force fluctuated periodically with six lobes, which was twice the number of impeller blades. The force center was essentially located at the center of the circle, indicating that the radial force of the impeller was stable and the torque of the pump shaft was relatively balanced. When BRADs existed, each design scheme presented seven obvious peaks and valleys. This result occurred because as the impeller passed through seven guide vane blades within one rotation period, the radial force changed seven times periodically.
Figure 18 shows that the dominant frequencies of the radial force on the impeller for the design schemes with deviations were different from those of the design scheme without BRADs and that they were related to the number of guide vane blades. When BRADs occurred, the amplitude of the dominant pulsation frequency was larger than that when there were no deviations, which could cause the torque of the pump shaft to be unstable.
Table 7 presents the average radial force on the impeller at flow rates of 0.8
Qdes, 1.0
Qdes, and 1.2
Qdes. In the absence of angle deviations, the radial force decreased with increases in the flow rate. However, the radial force increased with increases in the flow rate when there were angle deviations. For the condition of 1.0
Qdes, the average radial forces of Design Schemes I, II, IV, V, VI, and VII were 15.14, 7.34, 7.85, 14.39, 13.13, and 25.16 times that of Design Scheme III, respectively. In the radial velocity nephogram shown in
Figure 19, the area of radial high-speed zone under the condition of blade angle deviation is obviously larger than that under the condition of no deviation (design scheme III), which may cause an increase in the radial force of the impeller. These results show that the radial force increased dramatically when BRADs occurred. When the rotation angles of two blades were in opposite directions, the radial force generated on the impeller reached a maximum. Therefore, the occurrence of BRADs enhanced the resultant radial force on the impeller, which would lead to more serious swing and deformation of the pump shaft, thereby causing a certain degree of pump instability and shortening its service life.