4.2. Internal Flow Parameters in the Increasing Flow Rate Section
The velocity and total pressure for case 1 under the design flow rate was selected as the base data for the differential analysis. The variation between case 11 and 12 was the result of subtracting cases 11 and 12 from case 1. Additionally, cases 2 and 3 were the result of the subtraction of case 1 from case 2 and 3. Cases 11, 12, 2 and 3 were the cases in the increasing flow rate section.
Large variation values indicated huge changes in velocity/total pressure. In
Figure 7a,b, the velocity variations for cases 11 and 12 show a similar distribution. Fields with high velocity increases were mainly located at the leading edge and the trailing edge of the suction surface, and also at the leading edge and the middle region of the pressure surface. It is important to note that the velocity was observed to change more rapidly in the areas near the pressure surface.
In
Figure 7c,d, the velocity variation for case 2 and 3 reached the maximum in the “swirling” region, marked in
Figure 7d near the leading edge. During the increasing of the flow rate from case 1 to case 3, there was still a large velocity variation in the leading edge of the suction surface. In addition, the region with large positive velocity variation was located in the middle suction surface and extended out to the pressure surface’s trailing edge.
The differential analysis of the total pressure on the center plane for cases 11, 12, 2 and 3 was carried out. In
Figure 8a,b, with the increasing of the flow rate from a low flow rate to the design flow rate, the distribution of the total pressure variation shows an overall increasing tendency. The main areas of increase were concentrated around the leading edge of the suction surface and the trailing edge of the pressure surface. With the increasing of the flow rate, the region with increasing total pressure and the region with decreasing total pressure were mainly located around the leading edge of the suction surface. This was mainly affected by the small separating vortices generated at the blade’s leading edge.
Figure 8c,d shows the total pressure variation from the design flow rate to the high flow rate. During this process, the areas of increasing total pressure were mainly concentrated in the impeller passage and regions near the pressure surface. Conversely, in case 11 and 12, for the regions near the leading edge of the suction surface, the areas of increasing total pressure transferred to the areas where pressure was reduced. Both the fields with increasing total pressure and the fields with decreasing total pressure existed simultaneously in the vicinity of the leading edge of the suction surface. This phenomenon at the leading edge was mainly caused by the effect of the small-scale separation vortices.
In
Figure 9a,b, regions with strong vortex intensity were mainly located in the leading edge and the middle area of the suction surface. For the miniature centrifugal pump, there was a narrow space between the impeller and the back chamber. As a result of the rotor–stator interaction between the high speed impeller and the static back case, as shown in
Figure 9a,b, the flow field inside this narrow space had high vortex intensity. It is worth nothing that the regions with separation vortices near the leading edge of the suction surface showed good agreement with the regions with high variations of velocity and total pressure, as shown in
Figure 7 and
Figure 8. This indicated that the separating vortices had a huge effect on the local velocity and total pressure.
In
Figure 9c, in the middle area near the suction surface, one vortex core shifted into two vortex cores. This indicates that with the flow rate entering into the large flow regions, the shedding of the small-scale vortices on the suction surface became more frequent. This further led to the increasing of the number of small-scale vortices in the impeller passage. Furthermore, as a result of the geometry curvature of the impeller’s inlet, “circular ring” regions with strong vortex intensities appeared near the impeller’s inlet. The emerging “circular ring” regions could be explained in terms of the strong collision impact between the fluid particles and the impeller’s solid wall.
According to the energy gradient theory, the velocity gradient and the pressure gradient are one of the most important parameters. In
Figure 10, it can be seen that with the increasing of the flow rate, the regions with the high values for the energy gradient function
K under the periodic inlet flow presented an expanding trend. Large values were observed for the energy gradient function
K in the areas where the velocity/total pressure had drastically changed, which further characterized those areas that had high
K.
As shown in
Figure 10, the regions with
Kmax were located along the middle and the trailing edge of the pressure surface and the impeller exit near the trailing edge of the pressure surface. In
Figure 10, the “circular ring” it can be seen that regions near the impeller entrance also presented high mechanical energy gradients. This was the result of the sudden change of velocity direction at the impeller’s inlet as well as the strong collision impact between the fluid particles and the impeller’s solid wall.
4.3. Internal Flow Parameters in the Decreasing Flow Rate Section
In this section, the velocity and total pressure for case 7 under the design flow rate are selected as the data to be used in the differential analysis. Cases 5, 6, 8 and 9 are the cases in the decreasing flow rate section.
As presented in
Figure 11, during the process in which the flow rate decreased from a large flow rate to the designed flow rate, the main regions on the center plane showed a decreasing velocity. Fields with the increased velocity were mainly located in the areas of the leading edge of the suction surface, the middle position of the pressure surface and the exit of the impeller. In
Figure 11a, two regions can be identified with large velocity variation around the suction surface’s leading edge. This was mainly because the shedding vortices influenced the local velocity field and the pressure field, which further intensified the velocity change.
Figure 11c,d shows that with the flow rate decreasing from the design flow rate to a low flow rate, the distribution of the positive velocity variation was different from that of case 11 and case 12 in
Figure 7. This indicates that at the same flow rate, different tendencies in the variation of the flow rate led to different degrees of velocity variation.
The results of the differential analysis of the total pressure on the center plane for cases 5, 6, 8 and 9 are shown in
Figure 12. The total pressure variation for case 5 and case 6 was mainly negative, and the decreasing areas were located in the leading edge and the trailing edge of the pressure surface. In the leading edge of the suction surface, both the fields with increasing total pressure and the fields with decreasing total pressure existed simultaneously. With the flow rate decreasing from a high flow rate to the design flow rate, the areas with increased total pressure, for case 5 and case 6, were mainly concentrated around the leading edge and the trailing edge of the pressure surface. With the decreasing of the flow rate from the design flow rate to the low flow rate, the total pressure, for case 8 and case 9, exhibited an increasing trend. Furthermore, the regions with increased total pressure were located at the leading edge and the trailing edge of the pressure surface.
Regions with strong vortex intensity in the impeller for the decreasing flow rate section are shown in
Figure 13. Fields with strong vortex intensity were mainly located in the leading edge, the middle area of the suction surface as well as the narrow clearance between the impeller and the back chamber. In the leading edge and middle area of the suction surface, there were regions with different vortex intensities and different vortex scales, that is, the flow structures in these areas were dominated by strong vortex intensity.
In the decreasing flow rate section, the vortex intensity of internal flow was more complicated than that of the increasing flow rate section. In
Figure 13a,b, it can be seen that there was stronger vortex intensity for the shedding vortices near the leading edge of the suction surface.
Figure 14 shows the distribution of the energy gradient function
K in impeller for the decreasing flow rate section. In
Figure 14a,b, it is shown that fields with the
Kmax were mainly located in the middle area and trailing edge of the pressure surface. With the decreasing of the flow rate, regions with
Kmax inside the impeller showed a gradually decreasing trend. It is worth noting that, considering the above in combination with
Figure 1c,d, it can be seen that the vortices in the impeller passage affected the local velocity field and the total pressure field, which further influenced the local mechanical energy field and local shear strength. Furthermore, the changed mechanical energy and shear strength were then captured by the energy gradient function
K. In addition, the “circular ring” regions near the impeller’s inlet with
Kmax decreased gradually as the flow rate decreased.
4.4. Internal Flow Parameters for Special Positions of the Flow Function
After the differential analysis of the flow characteristics of the increasing flow rate section and the decreasing flow rate section, the flow characteristics of the special position of the periodic variable flow function were further observed and analyzed.
In this section, the two extremum positions of the flow function and the two design flow rates (the same flow rates with the different flow rate trends) were investigated, respectively. The corresponding relationship among the four cases was: case 1 was the design flow rate with a positive slope of the flow function, case 4 was the maximum point of the flow function, case 7 was the design flow rate with a negative slope of the flow function and case 10 was the minimum for the flow function.
In
Figure 15a,c, it is shown that the fields with strong vortex intensities presented a similar distribution for both case 1 and case 7 under the design flow rate. Both of them had two vortex cores in their separation vortices at the leading edge of the suction surface. However, in the middle area of the suction surface, the vortex intensity for case 1 was lower than that of case 7. The contour diagram for case 1 shows the single vortex core distribution, and the contour diagram for case 7 shows two vortex cores in the same region. In
Figure 15, the vortex core distribution inside the impeller for case 4 is shown to be more extensive than those of the other three operating conditions.
The distribution of the fields with strong vortex intensities were different for the fixed constant flow rate and the periodic inlet pulsation flow. For the same flow rate of 1.0Qn, fields with strong vortex intensities were mainly located in the outer layer of the large-scale vortices under the fixed constant flow rate, as well as the leading edge and middle area of the suction surface.
To further study the peculiar features under this rate of periodic pulsation flow of the inlet, time-domain analysis and frequency-domain analysis were carried out. The pressures measured from the monitoring points around the impeller in the four different cases are shown in
Figure 16. The monitors are numbered clockwise, and the monitor that is located near the impeller tongue is named monitor 1.
In addition, pressure pulsation at monitors near to/far from the pump tongue presented similar time-domain and frequency-domain characteristics. Therefore, the following investigation on the pressure pulsation is divided into two parts: the near tongue monitors (monitor 1, monitor 2, monitor 3 and monitor 8) and the far tongue monitors (monitor 4, monitor 5, monitor 6 and monitor 7).
In the following investigation, the pressure pulsation for both the constant incoming flow and periodic incoming pulsation flow is carefully studied and compared. Considering that the incoming periodic pulsation flow was constantly changing with time, the time-domain and frequency-domain analysis of the pressure were conducted based on the monitor data of ±0.5 rotation cycles around the analyzed cases. For cycles with variable flow regimes containing 36 centrifugal pump cycles, the pressure changes in the selected ±0.5 rotation cycles were negligible.
The impeller diameter
D1 of this pump was 33 mm, and the circumference diameter of each monitoring point was 35 mm. As shown in
Figure 16, eight monitors were evenly distributed in the circumferential direction. The miniature centrifugal pump had a design speed of 8000×
g rpm, the impeller’s rotational frequency
F = 133.33 Hz, and the blade-passing frequency
Fb = 5
F ≈ 666.6 Hz.
A periodicity characteristic of the time-domain pressure was observed, and is presented
Figure 17. The pressure value at monitor 2 was the smallest, but it showed the largest change in pressure amplitude among all the detection points. Monitor 2 was located in the area near the pump outlet and the impeller tongue, where unstable structures such as the “jet-wake” structure and the separation vortex could not be neglected. These unstable structures aggravated the complexity of local flows and were reflected in the range of pressure amplitude variation at monitor 2. Similarly, the circumferential position of the monitor points 4, 5, 6 and 7 were far away from the tongue, which resulted in the corresponding pressure amplitude variation range being small and the overall regularity being relatively uniform.
Under the different incoming flow rates, the main frequency of the pressure pulsation was always located at the blade-passing frequency Fb and the multiple blade-passing frequency m × Fb (m = 1, 2, 3…), and the amplitude of the other frequency division was much lower. For the same monitor (the impeller’s circumferential position), the pressure frequency at the 0.6Qn flow rate was higher than that in the design condition 1.0Qn and the large flow condition 1.4Qn, and the main pressure frequency tended to decrease with the increasing of the flow rate. This might have been due to the fact that the flow structure in the impeller was more complex under the small flow rate, and the larger dominant frequency peak may have been related to the collision of the multi-scale vortex structures.
Furthermore, monitor 1 and monitor 2 were located near the tongue, and their dominant frequency amplitudes at 0.6Qn were much higher than those at other monitors, which further verified that the internal flow conditions near the tongue were poor and that there was a huge pressure instability. With the increasing of the flow rate, the amplitude of pressure frequency at monitor 1 and 2 decreased significantly.
Figure 18 visualizes the transient Energy gradient function
K (define
Kmax = 30) performance in the pump impeller. Compared with the vortex core distribution (visualized by
QT), it can be seen that the distribution of the energy gradient function
K in space had a certain overlap with the vortex core distribution, especially in the outer layer of the large-scale vortices. In addition, the flow instability, which was affected by the wake-jet regions, was also reflected by the
K function.
At the times of T, T + 9 × ∆t and T + 18 × ∆t, the regions with Kmax agreed fairly well with the jet-wake structure and the large-scale vortices near the impeller tongue, which implied that the local flow instability was caused by the strong jet-wake structure and the influence of the impeller. When the time increased, the relative location of the impeller and impeller tongue was changed, which led to a strong interaction occurring in the local field. The corresponding changes in the separation passage vortices and the local flow instability were also reflected by the energy gradient function K.
The variation of the impeller pressure under periodic pulse flow was very different from that observed at a fixed constant flow rate. The different manifestations were as follows: the time-domain analysis of the pressure pulsation under the fixed flow rate was highly periodic, while that under periodic pulse flow was not. There was also a great difference between the pressure fluctuation spectrum range and the pressure fluctuation period. The large difference in pressure pulsation indicated that the flow characteristics under periodic pulse flow had their own features, and could not be fitted with the fixed constant flow rate cases.
Figure 19 shows the time-domain and frequency-domain analysis of pressure pulsation spectrum for case 1 and case 7. The two cases had the same flow rate and different flow rate trends. The monitors 1, 2, 3 and 8 near the tongue were selected to study the circumferential pressure performance.
Figure 19 shows that, for the same monitoring points, the pressure fluctuation spectrum amplitudes for case 1 and case 7 were different. In the frequency-domain analysis of pressure pulsation, the main pressure fluctuation frequency of case 1 was smaller than
Fb ≈ 666.6 Hz, and was close to the impeller’s rotational frequency
F = 133.33 Hz. The frequency-domain for case 1 remained in the low-frequency composition category, and the overall frequency composition was mainly concentrated at a level that was twice that of the blade-passing frequency. The main pressure fluctuation frequency of case 7 was close to its impeller’s rotational frequency
F = 133.33 Hz; the fractional frequency of the pressure was mainly concentrated at a level that was one to two times the impeller’s rotational frequency, denoted as
F. Compared with case 7, the amplitude of the dominant pressure frequency for case 1 was much larger. The generation, separation and fragmentation of vortices and the collision between flow particles were the source of low-frequency pressure pulsation. The case with the complex low-frequency pressure pulsation indicated that the internal flow was highly complex and unstable. Case 1 and case 7 showed different distributions in terms of the fractional frequency of the pressure. This indicated that the flow instability in the impeller was more likely to occur in the cases under the rising section than those under the decreasing section.
In
Figure 20, case 4 is the case with the maximum flow rate, and case 10 is the case with the minimum flow rate. The pressure fluctuation spectrum for case 4 was mainly located within one or two times the blade-passing frequency
Fb. Furthermore, there were some fractional frequencies of pressure that exceeded two times the blade-passing frequency
Fb for case 10. For the same monitoring point, the pressure amplitude for case 4 was also lower than that for case 10.
It is seen from
Figure 20 that the main pressure fluctuation frequency for case 4 was equal to the blade-passing frequency
Fb. Comparing the pressure spectrum for case 4 with that of 1.4
Qn under a fixed constant flow rate, the amplitude of main pressure fluctuation frequency under periodic pulse flow showed a smaller amplitude than that obtained under a fixed constant flow rate. Subjected to the effect of periodic pulse flow, many low-frequency components for case 4 were less than twice the blade-passing frequency (2
Fb), and they were mainly located around the impeller’s rotational frequency
F (133.33 Hz).
In addition, different from case 4, with the minimum flow rate, the main pressure fluctuation frequency for case 10 was around the impeller’s rotational frequency F = 133.33 Hz and the pressure amplitude was much higher than that for case 4. The components’ pressure pulsation spectra were concentrated at a level of two times the blade-passing frequency Fb.
The far more complex and lower amplitude of the fractional frequency of the pressure further illustrated the complex internal flow structures around the impeller at low flow rates. In the impeller, the unstable flow structures included the “jet-wake” structure and various scales of the separation vortices. The collisions and static pressure exchange between these unstable flow structures led to the local pressure and local pressure pulsation spectrum reaching its low flow value for case 10.
Finally, based on the energy gradient method, the distribution of the energy gradient function
K of the four cases with special flow rates were analyzed. The parameters of the flow inside the impeller were captured and investigated from the perspective of the energy field. In general, the distributions of the flow fields with large
K values were roughly the same for the four cases. As shown in
Figure 21, the regions with the
Kmax were mainly located in the middle and the trailing edge of the pressure surface, the exit of the impeller passage and the “circular ring” regions near the impeller’s inlet.
There was a certain change in regions with
Kmax in the middle and the trailing edge of the pressure surface. For the same flow rate, case 1 and case 7 showed different distributions of regions with
Kmax. In
Figure 21a,b, the area with the
Kmax at the outlet of the impeller for case 1 had a smaller scale than that for case 7. It is noteworthy that areas with the
Kmax for the middle and trailing edge of the pressure surface had a large change. This indicated that the flow behavior at the impeller exit was sensitive to the flow rate. Furthermore, “circular ring” regions with the
Kmax near the impeller’s inlet are shown for the four cases.
Comparing
Figure 10,
Figure 14 and
Figure 21, it was found that the most unstable section appeared in the first half-period of the flow rate variation (large flow rate, cases 2, 3, 4, 5 and 6), according to the distributions of Q criteria of the vortex and the energy gradient function
K. In this section, it was shown that the motions of strong vortices led to large gradients of the mechanical energy.