The results are presented in four different subsections. The first subsection discusses the aerodynamic performance of the tandem fan stage under different rotating speeds. In the second subsection, the change rule of the choke mass flow of the tandem fan stage with the rotating speed is determined, and the internal flow mechanism is investigated by establishing a physical model. In the third subsection, the reason for the difference of ηisen in the overspeed condition is clarified by analyzing the flow features near the mid-span and tip region in the tandem rotor domain. Finally, in the fourth subsection, the reasons for the difference in the spanwise distribution of ηisen under the partial speed condition are explored, and the radial migration trend of TV under the partial speed condition is clarified.
4.1. Aerodynamic Performance at the Design and Off-Design Rotating Speeds
At partial speeds, the inflow at the middle-higher span of the fan rotor gradually changes from supersonic to transonic and that at the lower span gradually changes from transonic to subsonic. When the rotating speed of the tandem rotor exceeds N4, the intensity of the shock wave in the blade passage increases and the inflow at the lower span is transformed into supersonic. Therefore, the aerodynamic performance of the high-throughflow fan greatly changes.
Figure 7 displays the performance curves under typical rotating speeds (N1, N2, N3, N4, and N5). From the change of
π, the growth rate is relatively stable, which is positively correlated with the rotating speed. In terms of
ηisen, the fan stage reaches a high level at N2. At this point,
ηisen decreases when the speed is either increased or decreased. As for
, it mainly involves two aspects: one is the compression intensity of the airflow in front of the aerodynamic throat, which can be characterized by
π and the second is the area of the aerodynamic throat. The aerodynamic throat of the fan stage is usually located in the stator. It is extremely difficult to obtain the area in 3D flow. Thus, only qualitative analysis can be conducted. As the rotating speed increases,
is positively related to the rotating speed due to the increased intensity of airflow compression. Since the analysis of aerodynamic performance from the characteristic curve is only qualitative, PEP and choke point at each speed were further selected for comparison.
To describe the relative change trend of different aerodynamic parameters under various speed conditions, the variation of parameters corresponding to each speed relative to N4 is calculated, and the expression is as follows:
where
x represents different aerodynamic parameters, including
π and
ηisen under PEP and
. The relative change obtained using the above expression is shown in a columnar diagram in
Figure 8. With increasing rotating speed,
π and
increase and the growth rate of
π is promoted; in contrast, the growth rate of
decreases, especially when N4 rises to N5, and the relative increment decreases to 3.1%, exhibiting a relatively abnormal change. As for
ηisen, it is consistent with the characteristic curve analysis, but
ηisen sharply decreases in the overspeed condition. However, under the partial speed conditions,
ηisen obviously improves, but the internal mechanism is not yet clear. Generally, different rotating speeds lead to massive changes in the typical flow structure inside the tandem rotor.
4.2. Over-Speed Operating Condition
Under the overspeed condition, the peak efficiency of the high-throughflow fan stage sharply decreases, which is about 2.3%. Therefore, the aerodynamic performance of each component needs to be separately analyzed.
Figure 9 shows the characteristic curves of
ηisen and the total pressure loss coefficient at N4 and N5 of the high-throughflow fan, and the PEP is denoted in the figure. The total pressure loss coefficient is expressed by
Cpt:
where
p represents the pressure, subscripts 2 and 3 represent the inlet and outlet of the stator, respectively, and superscript * represents the total parameter. The
ηisen of the tandem rotor considerably decreases under the overspeed condition, and the maximum difference is about 1.6%. The
Cpt of the stator slightly changes, with a difference of approximately 0.7%. Therefore, it can be preliminarily predicted that the main component causing the reduction of
ηisen is the tandem rotor domain, and the additional loss caused by the stator domain is not an essential factor. However, in-depth comparative analysis needs to be performed for the aerodynamic performance under PEP to determine the reason of the decline of
ηisen.
Figure 10 demonstrates the
ηisen and
Cpt histograms under PEP. As shown in the figure,
ηisen decreases by about 1.5%,
Cpt increases by about 0.7%, and the efficiency loss of the tandem rotor is more pronounced. Therefore, the subsequent analysis of the related flow structure will focus on the rotor domain.
First, by comparing the difference in the spanwise distribution of
ηisen under different rotating speeds, the regions to be analyzed were determined. The spanwise distribution of the tandem rotor at N4 and N5 is shown in
Figure 11, and the difference is obtained via subtraction between the two. The region near the tip exhibits substantial differences. The maximum difference is located at nearly 90% span, while the difference of other spans is relatively small, which is all within 3%. However, note that massive efficiency loss occurs near 60% span, which has not been observed for a traditional fan or compressor. Thus, this region needs to be analyzed.
The most direct change to the rotor with an increase in rotating speed is the increase in the tangential velocity, which subsequently increases the relative Mach number (
Mare) of the inflow. Therefore, the effects of the shock wave, wake, and leakage flow are all enhanced, which is the source of the sharp increase in losses near the tip region.
Figure 12 shows the
Mare and entropy distribution of 90% span. The flow features are the same under different speeds, including the two shock waves and the separation caused by the shock wave/boundary layer interaction. The contour shows that the entropy directly caused by the shock wave at N4 is broader than that at N5, but the shock wave intensity is higher under the overspeed condition. In contrast, the boundary layer loss caused by the shock wave is more pronounced under N5, and this loss is the main reason for the reduction in
ηisen. As the tandem rotor and stator are aerodynamically matched, the former domain is closer to the choke condition, and the shock wave position approaches the TE streamwise. Moreover, the effect caused by the stronger shock wave is accordingly enhanced, resulting in the thickening of the boundary layer and increased loss. Moreover,
Mare at the lower span will increase with the rotating speed, which will generate a shock wave and induce boundary layer thickening or flow separation. Thus, there must be more low-energy fluid, which will migrate radially under the effect of the centrifugal force, accumulate to the tip region, and diffuse downstream, thereby increasing the loss.
The entropy distribution along the streamwise under the overspeed condition is further demonstrated, and the 3D streamlines colored by Mare are emitted from the high-entropy region. It is confirmed that the streamlines through the high-entropy region of SS are from the radial migration of low-energy fluid, while those of PS are from the tip leakage flow of the adjacent blade and the boundary layer fluid on PS. The commonality of the above flow features is that they will deflect in the counter-rotation direction when entering the uncovered region, thereby forming a high-entropy region in almost the entire circumferential direction. The driving force of the deflection is the Coriolis force. As the rotating speed is determined, the circumferential Coriolis force is mainly determined by the radial velocity. Based on the above analysis, at 90% span, multiple effects of the stronger shock wave/boundary layer interaction, radial migration of a great quantity of low-energy fluid, and Coriolis force cause the reduction of ηisen under the overspeed conditions.
The source of the spanwise efficiency loss also needs to be further explored.
Figure 13 displays the entropy distribution along streamwise at the middle-lower span. After the booster rotor, except for the spanwise high-entropy region of the wake, a circumferential high-entropy region is present, but the scope is relatively small. As the airflow diffuses downstream, the intensity of the circumferential high-entropy region does not weaken, and its effect expands in both circumferential and spanwise directions. The streamlines through the high-entropy region are emitted from the booster rotor, which mainly comprises three sources. The first is the leakage vortex driven by the pressure difference between SS and PS at the tip, the second is the radial vortex formed by the radial migration, and the third is the low-energy fluid formed by the boundary layer on the tip of the booster rotor.
Since no shroud restriction is present, low-energy fluids are influenced by the leakage vortex, forming a stronger TV and diffusing downstream. With increasing rotating speed, the tip loading of the booster rotor will increase and the driving force will be strengthened, which will subsequently increase the vortex strength. Therefore, the loss of the shear effect with the surrounding mainstream will also increase. Furthermore, the wakes of the booster and fan rotors will broaden, thereby increasing the loss at the rotor outlet.
Figure 14 shows the
Mare and entropy contour of the tandem rotor at 55% span. The interaction between the booster and fan rotors increases under the overspeed condition. On the one hand, the intensity of the accelerated flow at the gap between the tandem rotors increases, but from the perspective of entropy, no additional loss is observed. In contrast, the enhancement of the shock wave intensity at the LE of the booster rotor increases the loss. On the other hand, the width of the booster rotor wake gradually extends after passing through the shock wave of the fan rotor, denoting the expansion of the high-entropy region. As shown in the figure, the wake widens in the overspeed condition. Therefore, the loss of 55% span is relatively high. In summary, the enhancement in the intensity of TV and wake under the overspeed condition will cause the reduction of
ηisen in the middle-lower span.
4.3. Partial Speed Operating Conditions
When the engine is in the ground idle condition, the fan/compressor operates at low speeds. Therefore, studying the flow mechanism of the high-throughflow fan at the partial rotating speed is extremely important. The overall stage aerodynamic performance curve of the fan shows that the peak efficiency of N2 reaches the maximum, and the peak efficiency will decrease when the speed is increased or decreased. Consistent with the overspeed condition analysis, the aerodynamic performance of each domain should be separately considered.
Figure 15 displays the characteristic curves of
ηisen and
Cpt under different rotating speeds, and the PEP of the overall stage is marked. As the rotating speed decreases, the
ηisen curve of the rotor gradually shifts to the upper left in the figure and the peak efficiency gradually increases; the minimum
Cpt does not apparently change at different speeds, and the curve approximately shifts to the left with decreasing speed. According to the PEP position, the matching point between the tandem rotor and stator gradually shifts away from the choke condition. Therefore, the
ηisen curve of the fan stage shows the change trend in
Figure 7, while the peak efficiency shows the change trend in
Figure 8.
To quantitatively verify this conclusion, the aerodynamic performance of the tandem rotor and stator corresponding to the overall stage PEP were selected and represented by a histogram (
Figure 16). The
ηisen of the tandem rotor gradually increases, but it reaches a steady state as the speed decreases. The
Cpt of the stator does not significantly change with the speed, but it suddenly increases by 0.17% when N2 decreases to N1. Therefore, even if the
ηisen of the tandem rotor increases at N1, the peak efficiency of the fan stage still decreases. However, note that the shifting of the matching point caused by the aerodynamic design is the fundamental reason for the decrease, and the loss caused by the stator is not the main reason for determining the aerodynamic performance. Therefore, subsequent analysis will mainly focus on the tandem rotor region.
As shown in
Figure 17, the
ηisen spanwise distribution of the stage interface at different rotating speeds is compared, and the
ηisen corresponding to the partial speeds is compared with N4. The region with apparent changes is still near the tip, which exhibits an increment of approximately 24% compared with N4, while the changes at other spans are all within 5%. However, at some of these spans, the change trend of
ηisen is difficult to interpret. For example, at 75% span,
ηisen suddenly decreases at N1, while other speeds do not exhibit such a phenomenon; for example, the reduction of
ηisen also occurs at 50% span. Therefore, the flow field needs to be further combined for judgment and analysis.
The
Mare and entropy contour at 95% span of the tandem rotor at N4 and N1 are demonstrated in
Figure 18. Clearly, as the rotating speed decreases,
Mare will transition from supersonic to transonic and gradually break away from the choke condition. Moreover, the strength of the LE shock wave weakens, and the shock wave and boundary layer loss are greatly reduced. On the other hand, the reduction in
Mare leads to decreased loading near the tip, thereby weakening the leakage strength and further reducing the loss. Therefore, the combined reasons for these two aspects result in the change trend of
ηisen.
Similarly, the change trend of
ηisen at 75% span from N4 to N2 is relatively regular, which is caused by the reduction of the inflow
Mare. However, if it is further reduced to N1, the drop in
Mare will increase the loss.
Figure 19 displays the
Mare and entropy distribution at 75% span of N1 and N2. From the flow field, the inflow
Mare decreases, which is consistent with the change trend of 95% span. However, at N1, the scope of the wake is broader than that at N2. At N1, no strong passage shock wave is present, and the scale of the boundary layer or separation caused by the shock wave is bound to be weaker. Nevertheless, the boundary layer rapidly thickens streamwise. Such a phenomenon occurs due to the excessively low inflow velocity. The ability to resist the adverse pressure gradient is restricted, and it is more likely to form the boundary layer accumulation. Additionally, the shock wave position is closer to LE, and the airflow is reduced to subsonic through the shock wave with insufficient acceleration. Therefore, a broader wake forms at N1 compared to that at N2.
To quantitatively describe the wake extent and the evolutionary process of the near-wall velocity, the entropy distribution in the circumferential direction at the rotor outlet and isentropic Mach number (
Maisen) distribution on the blade surface at different speeds are further obtained, as shown in
Figure 20. The figure shows that at different speeds, the beginning positions of the PS branch for the wake are basically the same, while those of the SS branch are quite different. As the rotating speed decreased, the wake width became broader. At N1, the wake occupied about 40% of the blade passage. Compared with N2 and N3, the wake and high loss region significantly expanded at N1. At partial rotating speeds, there is little difference in entropy in the region not affected by the wake. Therefore, the source of the
ηisen change trend at 75% span is identified as the change in the wake width. From the
Maisen distribution, after the airflow passes through the shock wave near the wall,
Maisen is subsonic at N1, and the deceleration and pressurization processes begin from about 20% chord length. Therefore, the growth rate of the boundary layer thickness increases, which can be confirmed by the above speculation.
The
ηisen at 50% span under N1 suddenly decreases. According to the previous analysis of the
ηisen loss region,
ηisen is bound to be affected by TV. Therefore, the spanwise distribution of entropy is plotted further along streamwise at N4 and N1 in
Figure 21. The beginning position of TV is nearly at the same span under different speeds. As the speed decreases, the span of TV at the rotor outlet gradually decreases, corresponding to the
ηisen distribution of the tandem rotor. This flow phenomenon is considered to be related to the leakage flow of the booster rotor. When the rotating speed decreases, the loading at the booster rotor tip inevitably decreases, leading to a reduction in the pressure difference and thereby reducing the tangential velocity of the leakage flow. Therefore, the inertial force generated by the tangential velocity of the airflow is insufficient for enabling the TV to migrate higher and needs to be balanced with the radial adverse pressure gradient, resulting in a gradual reduction in its radial position at the rotor outlet.
To prove the above speculation about the migration trend of TV, the cylindrical coordinate form of the Navier–Stokes equation is employed. This paper focuses on the radial migration of TV; thus, the radial component equation is finally selected for further analysis and discussion, and it is expressed as follows:
where
w represents the relative velocity, the subscripts
r,
θ, and
z represent the radial, circumferential, and axial directions, respectively,
ω represents the angular velocity, and
ρ represents the density. The three terms on the left side of the above expression are acceleration terms. On the right side, the first term is the adverse pressure gradient term, the second term is the centrifugal force term, the third term is the Coriolis force term, the fourth term stems from the cylindrical coordinate system, and the fifth term is the viscosity term. Furthermore, the second and third terms on the right side are inertial force. In terms of the expression, the fourth term can be equivalent to the centrifugal force. Therefore, the second to the fourth terms are named as the equivalent inertial term. Additionally, the equivalent inertial term can be rewritten as
where
u represents the local tangential velocity. Based on the expression, the equivalent inertia term must be greater than 0 to resist the radial pressure gradient and move the fluid radially upward.
Figure 22 displays the distribution of the adverse pressure gradient and equivalent inertial force streamwise of the tandem rotor at N1. Additionally, the entropy isoline in the figure denotes the TV profile. The upper half of TV presents a positive pressure gradient and a higher equivalent inertial force, while the lower half presents an adverse pressure gradient and a lower equivalent inertial force. Therefore, the pressure gradient and equivalent inertial force are highly correlated. Since TV is driven by the pressure difference, the upper half of TV exhibits a larger relative circumferential velocity, while the lower half exhibits a smaller one. Therefore, after subtracting and squaring the relative circumferential velocity with the tangential velocity near the TV region, the upper and lower half regions are clearly divided. The figure also shows that the vortex core is roughly located at the positive pressure gradient region. Moreover, the left side of the expression is basically determined by the positive pressure gradient and equivalent inertial force (compared to other terms on the right side, the order of the viscous term in the main flow region is lower; thus, it can be ignored). As the rotating speed decreases, the intensity of the equivalent inertial force at the vortex core gradually decreases; therefore, the radial position of TV changes to downward, verifying the previous speculation.
Figure 23 demonstrates the
Mare and entropy contour at 50% span of the tandem rotor for N4 and N1. As the rotating speed decreases to N1, no shock wave forms in the blade passage at 50% span, the accelerated flow in the nozzle weakens, and no additional loss occurs in the mainstream. The difference between the two cases is mainly in the wake width, and the wake loss at N4 is higher than that at N1. Furthermore, a large-scale high-entropy region (white dotted circle) near the pressure surface (PS) of the fan rotor exists at N1, representing TV, which confirms the analysis of the radial migration trend of TV in this paper.
4.4. Gain of Choke Mass Flow at Different Rotating Speeds
A booster rotor is introduced into the fan stage to mainly improve the flow capacity near the hub.
Table 3 shows that, with N3 as the dividing line, when the speed exceeds it, the tandem fan stage can improve
, which is the flow capacity. However, the growth rate under the overspeed condition slows down. When the rotating speed drops below N3, the introduction of the booster rotor will have a negative effect. Therefore, the aerodynamic advantages of the high-throughflow fan can be gradually demonstrated at high speeds. This subsection combines aerodynamic theory with flow field analysis to clarify the change trend of the growth rate for
, fundamentally explain the problem, and reveal its internal flow mechanism.
In aerodynamics, it has been indicated that under the condition of subsonic inflow with a certain total temperature and total pressure, once the critical section (
Ma = 1) of the Laval nozzle appears, it will reach the choke condition and
cannot be increased by further reducing the back pressure. Under the choke condition, the airflow generally has a low total pressure loss before passing through the throat in the stator. In other words, the
π of the rotor can approximately measure the inflow compression intensity before the throat.
Figure 24 displays the relative change in the
π and
of the tandem rotor compared to that of the fan rotor inside the baseline fan stage under the choke condition. The gain of
π is always positive, while that of
is negative below N3. Such a phenomenon occurs due to the introduction of the booster rotor, which causes excessively high loading near the stator hub. To suppress the flow deterioration, the hub profile of the baseline fan stage is modified, that is, the slope angle of the hub profile in the stator domain is lifted to 15° to obtain the hub profile of the new concept tandem fan stage. Therefore, the aerodynamic throat area of the overall stage is greatly reduced. The gain of
π shows a quadratic trend with increasing speed. At N4 and partial speed conditions, the gain of
is positively correlated to that of
π. When the speed is increases to the overspeed condition, the gain of
is abnormal. Therefore, it can be inferred that the throat area contributed to the generation of such a trend.
The booster rotor was introduced to increase the core pressure ratio. Since the tangential velocity near the hub of the rotor domain is insufficient, the tandem rotor implants energy into the airflow by increasing the flow turning, thereby enhancing the compression intensity of the airflow before the aerodynamic throat. During the aerodynamic design, to improve the compression capability near the hub, the total camber angle of the tandem rotor (53°) is much larger than that of the baseline fan rotor (33°). Therefore, the reaction near the tandem rotor hub is reduced compared to that at the baseline fan rotor, and the airflow velocity at the stator inlet accordingly increases. The stator inflow can even reach the supersonic condition at N4 and N5. Therefore, a great negative incidence under choke condition was designed to avoid premature corner stall. In contrast to the subsonic inflow in the Laval nozzle, the determination of the aerodynamic throat position under such inflow condition is rarely studied; thus, physical modeling analysis is required.
Figure 25 presents a schematic to analyze the aerodynamic throat position. The green lines in the figure represent the detached shock wave, and the red and blue dashed lines represent the expansion waves emitted by the suction surface (SS) and PS, respectively. An expansion wave from PS will first intersect with the LE of the stator. Although the disturbance after the expansion wave cannot affect the flow field in front of the blade, the disturbance from SS can still be transmitted upstream. Therefore, when an expansion wave emitted from SS intersects the LE of the adjacent blade, the flow field in the channel will no longer affect the upstream. Then, the position of the aerodynamic throat can be determined. However, the two expansion waves from SS and PS will intersect with a high probability. Thus, the lower half of the throat profile is determined by the expansion wave from SS, while the upper half is determined by the expansion wave from PS.
Figure 26 displays the aerodynamic throat distribution near the hub and tip under the choke condition for the high-throughflow fan and baseline conventional fan, where the throat positions at different rotating speeds are denoted with dashed lines. When the rotating speed under the partial design speed condition is increased, the throat position of the high-throughflow fan stator and the baseline stator move upstream, signifying an increase in the throat area. Note that when N3 rises to N4, the throat position at the hub does not apparently change. However, since the Mach number of the supersonic inflow at the throat is greater than one, the gain of
at N4 rapidly increases. However, as N4 increases to N5, the throat position of the high-throughflow fan does not significantly change either at the hub or tip. Compared to the increase of the throat area for the baseline, the throat area of the high-throughflow fan nearly does not change under the overspeed condition, denoting a relative reduction of the throat area. Although the change of the rotor pressure ratio still increases at this time, the growth rate of the gain of
will generally be abnormal.