3.1. Circumferential and Axial Analysis of the Steady Velocity and Pressure Fields
Figure 1b shows the location of the section in
Figure 4. The internal flow characteristics of the centrifugal compressor are analyzed from the axial x = 0 section based on the velocity streamline. At the high mass-flow condition (1.1 Qn = 30.56 kg/s), a distinct low-velocity region mainly occurs near the vaneless diffuser and tongue as illustrated in
Figure 4a. The streamline deviates from the arc of the wall in the volute and strikes the inner wall of the volute, causing a lot of impact losses. Compared with small mass-flow rates, the velocity of the impeller is obviously higher than that of other regions and the velocity of the vaneless diffuser is obviously lower than that of other regions.
Figure 4b exhibits the velocity stream line of the x = 0 mm section of a centrifugal compressor under the designed conditions (Qn = 28.27 kg/s). The streamline of the centrifugal compressor is very smooth in the designed conditions.
Figure 4c shows the velocity stream line of the x = 0 mm section of a centrifugal compressor inthe low mass-flow rate case of 0.9 Qn = 25.48 kg/s. It is clearly seen that a little separated flow appears at the suction surface of the impeller. The streamline at the blade tip gap is chaotic, and the incident angle of the fluid is inconsistent with the blade design angle, which reduces mainstream energy. The fluid is aggravated by viscous force, and boundary-layer separation occurs, which is confirmed by the lower impeller speed. At the junction of the impeller and the vaneless diffuser, large velocity gradients are exhibited, and the fluid mixing losses increase, which inevitably leads to acomplexity of flow in the vaneless diffuser. The axial cross-sectional schematic diagram shows that eddies exist in the vaneless diffuser, and the velocity is significantly lower than that of the design conditions in the volute.
Figure 4d plots the velocity streamline of the x = 0 mm section of a centrifugal compressor in the lower mass-flow rate case of 0.8 Qn = 22.62 kg/s (near-stall condition). It is clearly observed that with the decrease in mass-flow rate, flow separation occurs on the blade pressure surfaceof the impeller, and back-flow is even observed at the inlet of the impeller. Moreover, the flow instability is more serious in the vaneless diffuser, and the low-speed area becomes increasingly wide in the volute. To sum up, the blade tip clearance is a sensitive zone for unstable flow. If the mass-flow rate deviates from design conditions, the airflow angle will not match the blade design angle, which enhances the energy loss of the mainstream and causes flow instability inside the centrifugal compressor.
Figure 5 displays the distribution of the velocity and streamline in the Z = −12 mm section at different mass-flow rates.
Figure 5a shows the distribution of the velocity and streamline at the Z = −12 mm section and the high mass-flow rate (1.1 Qn = 30.56 kg/s). The high mass-flow rate represents excessive velocity in the impeller. Slight flow separation occurs at the end of the blade channels with the streamlines converging towards the pressure surface. Terrible velocity gradients are observed at the junction of the impeller outlet and the vaneless diffuser, which results in considerable mixing losses. The low-speed zones and eddies near the volute tongue are noted as causing blockage of the centrifugal compressor. The presence of vortices indicates that the fluid needs to overcome a large negative pressure gradient to work, which is not conducive to the development of turbulence.
Figure 5b shows the distribution of the velocity and streamline at the Z = −12 mm section in the design condition (Qn = 28.28 kg/s). The streamline of the volute and impeller is smooth and the permeability of blade passage is excellent. At the same time, the vaneless diffuser shows the positive effects of speed reduction and pressurization.
Figure 5c displays the streamline distribution with velocity contour as background at the Z = −12 mm section in the 25.48 kg/s condition (0.9 Qn). The overall velocity in the impeller decreases and secondary flow occurs in the two flow channels located at 0°. At a low mass-flow rate, flow separation occurs preferentially in one or several channels, which is confirmed by Dodds [
7]. The upper instability of the vaneless diffuser under centrifugal force was produced by the low-speed fluid at the impeller outlet.
Figure 5d exhibits the distribution of the velocity and streamline at the Z = −12 mm section at the lower mass-flow rate of 0.8 Qn = 22.62 kg/s (near-stall condition). When the mass-flow rate is further reduced, the density of streamlines changes at 60° and 210°, which is caused by the non-uniformity of the velocity gradient. This phenomenon develops from the outlet of the impeller to the vaneless diffuser. Additionally, slow flow is observed in more blade channels. The airflow angle changes at low mass-flow rates, while the back-flow of the tip clearance reduces the energy of the main flow, resulting in a more pronounced effect of viscous force and more severe boundary-layer separation. Furthermore, as the mass-flow rate decreases, the stall extent and area of the vaneless diffuser are further expanded, and large-scale, low-speed vortices exist, causing more mixing losses. Flow instability is exacerbated downstream of the vaneless diffuser, which is consistent with the work by Elias et al. [
29].
Our research reveals the stall mechanism in two different off-design conditions. For high mass-flow rate conditions, high-speed fluid impacts the volute tongue, resulting in a large amount of impact loss eddies beingformed at the volute tongue, which causes the blockage of the centrifugal compressor. As for low mass-flow rate conditions, the flow separation at the trailing edge of the blade is caused by the backflow at the tip clearance As the mass-flow rate decreases, flow separation intensifies, which gives rise to the main flow energy becoming too low.
Figure 6a–d show the contours corresponding to the static pressure value of four different mass-flow rates at the z = −12 mm section.
Figure 1b presents the cross-section.
Figure 6a exhibits the static pressure contour of the high mass-flow rate(30.56 kg/s). A considerable radial static pressure gradient mainly appears in the impeller. Compared with the low mass-flow rates, the case of 1.1 Qn = 30.56 kg/s shows globally low pressure and a slight static pressure gradient across the vaneless diffuser and volute.
Figure 6b shows the static pressure distribution under the design conditions (Qn = 28.27 kg/s). Naturally, the unique diverging runners of the centrifugal compressor produce a superb deceleration and supercharging effect. Furthermore, the pressure gradient of the vaneless diffuser and the volute is stable.
Figure 6c shows the static pressure contour at mass-flow rate of 0.9 Qn = 25.48 kg/s. It is clearly seen that the vaneless diffuser and volute exhibit a higher static pressure area; the radial static pressure gradient of the vaneless diffuser and volute are larger than the design condition. However, there is no such obvious phenomenon in the impeller. A zone of high pressure is observed at the volute tongue, which is a sensitive zone for pressure fluctuations. We found a low-pressure region appearing at 45°; it can be related to the jet-wake structure. Unstable flow is frequently found in vaneless diffusers, as previously shown by Hellstrom et al. [
30]. They argued that the low-speed area in the vaneless diffuser and volute is mainly caused by a large pressure gradient.
Figure 6d exhibits the static pressure distribution at near-stall conditions (0.8 Qn = 22.62 kg/s). As shown in
Figure 6d, one can see that the radial pressure gradient of the vaneless diffuser and volute obviously enhance with a decrease in the mass-flow rate. A dramatic pressure gradient appears circumferential of the diffuser under unsteady flow, as explained by Liu et al. [
31].
At a high mass-flow rate, more energy exists in the form of kinetic energy, which is why the static pressure of the centrifugal compressor is lower than that of the design conditions. The static pressure gradient is almost slight, while the radial static pressure gradient is enhanced in the vaneless diffuser and volute. Nevertheless, the influence of the pressure gradient in the impeller is milder, which is confirmed by Semlitsch et al. [
32]. For low mass-flow rates, the static pressure value is greater than that of the design conditions, which is a symbol of low-speed flow. In addition, the excessive static pressure gradient appears at the vaneless diffuser of centrifugal compressor, and the work required to overcome the adverse pressure gradient is unreasonable.
Figure 7 represents the distribution of turbulent kinetic energy at 50% blade height at four different mass-flow rates.
Figure 7a shows the distribution of turbulent kinetic energy at the high mass-flow rate (1.1 Qn = 30.56 kg/s). It is obvious that a large amount of TKE is generated at the blade’s outer diameter at the high mass-flow rate (1.1 Qn = 30.56 kg/s). Turbulent kinetic energy exists for overcoming large energy dissipation to continue the development of turbulence, which cause kinetic energy loss of the mainstream. As shown in
Figure 6a, the massive energy is lost due to the excessive back-pressure gradient at the impeller outlet.
Figure 7b shows that the status of turbulent kinetic energy underthe design conditions (28.27 kg/s). When the mass-flow rate reduces to 25.48 kg/s (see
Figure 6c), TKE is generated on the suction surface of several flow runners due to the flow separation. The development of turbulence is maintained by the resulting TKE in poor flow channels. Flow channels that generate TKE are often identified as appearing with flow separation, which is an unfavorable phenomenon coming at the expense of mainstream energy.
Figure 7d displays the statistic of turbulent kinetic energy at the lowest mass-flow rate (0.8 Qn = near-stall condition). It can be noted that more channels produce TKE with decreased mass-flow rates. In general, the presence of the TKE is an adverse phenomenon for the efficiency of the centrifugal compressor. Excessive TKE production means clogged centrifugal compressors at the high/low mass-flow rates. From the analysis of
Figure 6, the work provides an insight into the flow instability of impellers under different off-design conditions. The high flow losses in the impeller mainly occur in different suction surfaces of the blade passage at high and low mass-flow rates.
3.2. Unsteady Analysis of Internal Flow in the Centrifugal Compressor
For unsteady computational data acquisition, the time-step size is set to 3.16 × 10
−5 s in this paper. We obtain a set of data for every two degrees of rotation. The unsteady numerical results of the centrifugal compressor after 10 revolutions are adopted to ensure the accuracy of the numerical simulation.
Figure 8 shows the development of the velocity field and the distribution of streamlines at the z = −12 mm cross-section and the high mass-flow rate (1.1 Qn = 30.56 kg/s). The characteristic of the high mass-flow rate is blockage; as shown the red elliptical area in the
Figure 8, the velocity contours present the low-velocity areas and eddies near the volute tongue. From
Figure 8a–c, we can clearly see the development process of the eddies at the tongue, which is decomposed into two small eddies from a large eddy. The separation of the vortex is more conducive to the circulation of the centrifugal compressor. As the separation of the vortex progresses, the low-speed blocking area at the outer ring of the volute tongue disappears and the low-energy fluid is discharged from the outlet of the volute, as can be observed in
Figure 8c. This provides insight into turbulent energy transfer. The large eddies cause serious destruction of the average mainstream, the existence of eddies inevitably leading to a sudden change in the pressure gradient where a large amount of the energy conversion takes place. The mechanism of the small eddy is mainly to convert the energy of the large eddy into heat and dissipate it. Moreover, the large eddy current exhibits anisotropy with certain stagnation characteristics, while the isotropic performance of the small eddy current is indicated. These small vortices are better able to drain the compressor with the mainstream. Compared with the design conditions, a larger flow velocity in the impeller can be noted with a high mass-flow rate. The flow separation in the blade channel is mainly manifested in the outer diameter of the blade. Serious boundary-layer separation is observed at the suction surface of the trailing edge of the blade, resulting in a jet structure. Under the action of centrifugal force and Coriolis force, bundled streamlines are observed and the streamlines show inferior conformability with the inner wall of the volute.
Figure 9 exhibits the velocity field situation and streamline distribution at the z = −12 mm cross-section and the near-stall mass-flow rate (0.8 Qn = 22.62 kg/s). In general, acomplexity of internal flow is observed at a low mass-flow rate. The extreme inhomogeneity of the velocity field brings about an excessive velocity gradient. At a constant speed of 10,545 rpm, centrifugal and Coriolis forces intensify the effect on low-energy flow in the impeller, accompanied by the squeezing of high-speed groups in the vaneless diffuser—which reveals that most streamlines are distributed in bundles.
Figure 10 shows the distribution and development of pressures at the z = −12 mm cross-section and near-stall mass-flow rate (0.8 Qn = 22.62 kg/s). A widely low-pressure zone at 45° is seen and the appearance of the lowpressure leads to excessive energy exchange. With time evolution, the low-pressure zone gradually increases. From
Figure 10a–c, the high-pressure region at the vortex tongue gradually expands, resulting in a greater squeeze of the downstream fluid on the main flow. In addition, almost all blade channels in the impeller present low-speed zones and flow separation, which is caused by the stratification of the suction-surface boundary layer due to the low mainstream energy.
From
Figure 9a–c, one can see that the fluid is deflected circumferentially and the flow instability within the impeller is exacerbated, which is a precursor to the appearance of secondary flow. Overtime, flow deterioration is observed at 120° of the vaneless diffusers and at the volute outlet. Abrupt changes occur in the velocity gradient, resulting in more energy losses and lower efficiency than that of design conditions.
Figure 11 describes how three points are selected in a flow channel at 30%, 50%, and 90% blade height.
Figure 12 shows the velocity fluctuation of three detection points located at different axial heights (30%bh, 50%bh, 90%bh) in a blade channel, and the straight line is the average value of the fluctuation. According to the unsteady calculation, the temporal and spatial variation of the internal flow in the centrifugal compressor impeller is analyzed to provide some flow characteristics under off-design conditions.
Figure 12a exhibits the velocity fluctuations of different blade heights (30%bh, 50%bh, and 90%bh) under the high mass-flow rate conditions (1.1 Qn = 30.46 kg/s). Velocity fluctuation is not violent at 30% blade height, and the flow is relatively stable in the tip clearance, indicating that small impacts are caused by the internal-flow instability of the tip clearance at the high mass-flow rate. The velocity fluctuation of the monitoring point at 50% and 90% of the blade height is very strong, which is caused by the velocity field disturbance in the boundary-layer separation of the impeller suction surface. Furthermore, the velocity fluctuation at 50% blade height is obviously higher than that of 90% blade height. This is mainly due to the self-regulating ability of the expanding channel of the centrifugal compressor. Flow instability mainly occurs in impellers and vaneless diffusers, which is well-consistent with that studied by Lou et al. [
33]. The locations of these starting points of flow disorder are refined in the work. In general, the boundary layer separation of impeller suction surface appears at 30-to-50% blade height at a high mass-flow rate. The flow instability characteristics at 50% blade height are more serious than those at 90% due to the narrow diameter.
Figure 12b shows the velocity fluctuation of different blade heights (30%bh, 50%bh, and 90%bh) at the near-stall mass-flow rate(0.8 Qn = 22.62 kg/s). Compared with 50% blade height, the fluctuation of the monitoring points at 90% blade height becomes moderate, indicating that the starting point of the velocity disturbance is at 30–50% blade height. The velocity fluctuates drastically at the monitoring points of 50% blade height. In the same near-stall mass-flow rate case (0.8 Qn = 22.62 kg/s), the fluctuation of the trailing edge of the blade at 90% blade height is more stable than that at 50%blade height. By comparing velocity fluctuations, the fluctuation amplitude of the impeller at the near-stall mass-flow rate is more apparentthan that at the high mass-flow rate (30.46 kg/s). Flow separation occurs earlier at the near-stall mass-flow rate when comparing the velocity fluctuation of 30% blade height in
Figure 12a,b. At 50% blade height, the average flow velocity forthe high mass-flow rate is much larger than that ofthe near-stall mass-flow rate, indicating that more serious diversion losses and lower mainstream energy appear under near-stall conditions (0.8 Qn = 22.62 kg/s).
Fast Fourier transform (FFT) is implemented to analyze unsteady flow in rotating machineries, which can provide some insights to understanding the internal flow by transforming the time domain into the frequency domain.
Figure 13 plots the amplitude spectra obtained by FFT for pressure fluctuations at points A1, A2, and A3 at different mass-flow rates.
Figure 13a exhibits the amplitude spectra of monitors at the high mass-flow rate (30.46 kg/s). The pulsation at 50% and 90% of the blade height is dominated by frequency conversion. However, high amplitudes occur at90% of the blade height of the impeller outlet.
Figure 13b shows the amplitude spectra of monitors at near-stall conditions (0.8 Qn = 22.62 kg/s). The pulsation at 30%, 50%, and 90% of the blade height is dominated by frequency conversion. The rotational stall frequency and frequency multiplication show a difference in off-design conditions. The amplitude of pointA2 (50% BH) is lower than that of pointA3 (90% BH) at the same mass-flow rate, but the number of characteristic frequencies at point A2 is more than that at point A3. In general, the stronger instability occurs in the impeller of the centrifugal compressor at near-stall conditions (0.8 Qn = 22.62 kg/s).
Figure 14 displays the amplitude spectra of monitoring point A3 under the two off-design conditions. In
Figure 14a,b, the amplitude of point A3 at the impeller outlet is analyzed. The stall frequency is calculated to obtain a high mass-flow rate of 21.2 Hz, while the near-stall mass-flow rate is 31.8 Hz. The frequency spectrum of the impeller outlet at the high mass-flow rate (30.46 kg/s) is calculated in
Figure 14a. It is obtained that the first, second, and third characteristic frequencies are 3.014 fn, 1.025 fn, and 1.989 fn, respectively. As shown in
Figure 14b, the frequency spectrum of the impeller outlet at 20.62 kg/s (near-stall) is analyzed, and its characteristic frequencies are 1.989 fn, 3.014 fn, and 3.979 fn, respectively. These characteristic frequencies may be related to boundary-layer separation in the blade channels and the jet at the impeller outlet. Comparing to the two different mass-flow rate conditions, there is a characteristic frequency with a prominent amplitude at the near-stall mass-flow rate.
Figure 15a shows the velocity fluctuations at the monitoring points of the tongue and the four calculation cases. Compared with the design conditions (28.27 kg/s), the velocity fluctuation is clearly lower near the tongue at the high mass-flow rate. Moreover, several amplitudes deviating from the mean value are observed, which is closely related to low-velocity eddies at the tongue (see
Figure 8). At the near-stall mass-flow rate, periodic high-amplitude fluctuations are exhibited at the tongue, which aggravates with decreased flow.
Figure 15b shows the velocity fluctuation at monitoring points of the outlet. It is further found thata sudden, irregular velocity fluctuation appears in the outlet of the centrifugal compressor at all off-design conditions, but the peak value of the fluctuation gradually increases with the distance from the design working point, which quantitatively reveals that the stall cell is gradually increasing at the cochlear tongue and exit sites of the centrifugal compressor, and the intensity of the stall cell blockage is also enhanced with decrease inflow rate.
Figure 15 further verifies that blockage occurs near the outlet and the internal worm tongue of the compressor with decreasing flow rate in
Figure 5.
Figure 16 shows the pressure fluctuations at the monitoring points of tongue and outlet at different mass-flow rates. The same pressure-fluctuation trend is found at the tongue and outlet. However, the change influctuation frequency at the tongue is greater than that at the outlet. At the high mass-flow rate, the pressure fluctuation at the volute tongue and outlet is more intense than that at the near-stall mass-flow rate (0.8 Qn = 22.62 kg/s). The fluctuation of pressure is more intense and many spike signals appear at the near-stall mass-flow rate when compared with the design condition (Qn = 28.27 Kg/s). The unstable flow complexity of the centrifugal compressor volute is demonstratedunder off-design conditions. Pressure fluctuations at the volute deteriorate the flow field, which is well-consistent with the previous work by Guo et al. [
34].
Figure 17 shows the schematic diagram of the monitoring points on the vaneless diffuser. The data of these monitoring points are collected to analyze the spectrum of circumferential pressure fluctuations on the vaneless diffuser by FFT.
Figure 18 displays the amplitude spectrum of monitoring points of the vaneless diffuser for one circumferential turn at two different mass-flow rates.
Figure 18a illustrates the amplitude spectrum at the high mass-flow rate (1.1 Qn = 30.46 kg/s). Illustrated in
Figure 18a, it is noted that the pressure fluctuation period in the vaneless diffuser is mainly dominated by the blade passing frequency and its frequency multiplication. Meanwhile, the rotating stall frequency is almost negligible.
Figure 18b shows the spectral analysis of monitoring points regarding pressure fluctuations of the vaneless diffusers at the near-stall mass-flow rate, including blade passing frequencies and theirfrequency multiplication. Meanwhile, numerous low-frequency fluctuations are observed in the vaneless diffusers of the centrifugal compressor. The amplitude value of the blade frequency at different locations is distinguished at the near-stall mass-flow rate (0.8 Qn = 22.62 kg/s), which is related to the jet-wake structure in the vaneless diffuser. These things considered, the low-frequency amplitude is significantly smaller than that of the blade frequency.
The blade-frequency amplitude value at the high mass-flow rate is greater than that in near-stall mass-flow rate conditions. In addition, the fluctuation is dominated by the blade frequency and its frequency multiplication; substantial low-frequency fluctuations exist at the near-stall mass-flow rate, which is not observed at the high mass-flow rate.
Figure 19 describes six monitoring points position in the volute of the centrifugal compressor. A circle of monitoring points (D1, D2, D3, D4, D5, and D6) are measured inside the volute along the flow direction.
Figure 20 plots the amplitude spectrum of the monitoring points of the volutefor one circumferential turn at two different mass-flow rates.
In
Figure 20, the characteristic frequencies of monitoring points’ amplitude spectra are shown to be low at these two mass-flow rates, while the blades’ passing frequencies increase. Additionally, the low-frequency amplitude value is significantly larger than that of the blade-frequency amplitude. The low characteristic frequency amplitude at the high mass-flow rate is higher than that at the near-stall mass-flow rate. The amplitude at the near-stall condition is higher than that at the high mass-flow rate from points D1 to D3. It is worth nothing that the amplitude spectrum from points D1 to D3 is more complex, which is caused by the jet-wake in the vaneless diffuser.