Next Article in Journal
Electro-Thermal Parameters of Graphene Nano-Platelets Films for De-Icing Applications
Previous Article in Journal
TDE-Based Adaptive Integral Sliding Mode Control of Space Manipulator for Space-Debris Active Removal
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Blade Fracture on the Flow of Rotor-Stator Systems with Centrifugal Superposed Flow

1
School of Energy and Power Engineering, Beihang University, Xueyuan Road No. 37, Haidian District, Beijing 100191, China
2
Research Institute of Aero-Engine, Beihang University, Xueyuan Road No. 37, Haidian District, Beijing 100191, China
*
Author to whom correspondence should be addressed.
Aerospace 2022, 9(2), 106; https://doi.org/10.3390/aerospace9020106
Submission received: 17 January 2022 / Revised: 10 February 2022 / Accepted: 13 February 2022 / Published: 16 February 2022

Abstract

:
Rotor-stator cavities are often found in turbomachinery; they supply cold air that is bled from the compressor to the turbine blades. The pressure of the outlet of a rotor-stator cavity is axisymmetric under normal circumstances. However, its pressure would be non-axisymmetric in the event of blade fracture. The impact of blade fracture on a rotor-stator cavity with centrifugal superposed flow is studied in this paper. The Euler number E, the rotational Reynolds number R e φ , and the low-pressure zone range θ are investigated and, for the first time, with the non-axisymmetrical boundary conditions employing numerical simulation. The results of the numerical calculations show that after turbine blade fracture, the velocity is more affected in the downstream region at a high radius, especially when the R e φ is large. As for the distribution of the mass flow rate, there may be a critical θ c at which the other blades are least affected. The θ c would increase as the R e φ or the E increase, and the θ c 0.2 when C w = 10 , 137 , R e φ = 5.12 × 10 5 , and 0.2 E 0.4 . In addition, the thrust coefficient increases as the E or the θ increases, and the increase in the thrust coefficient does not exceed 4% when the E = 0.2   and   the   θ = 0.1 in this paper. However, the moment coefficient on the rotating shaft is almost independent of the E and the θ. An increase in the R e φ will reduce the effect of turbine blade fracture on the thrust and moment coefficients, when the R e φ is small.

1. Introduction

By raising turbine-entry temperatures, some axial turbomachines may reach an efficiency of approximately 80% with well-designed impellers. While part of the temperature increase may be attributed to the discovery of new materials, the majority would be attributed to advancements in cooling technologies. For this purpose, a small amount of the air collected from the compressor is utilized to cool the nozzle guide vanes, turbine blades, and disks. This type of flow may be found in almost all turbomachinery cavities positioned between the rotating impeller and the stationary disk, which is called the rotor-stator system.
Concerns about the design of turbines are divided into two categories: blade design and rotor-stator cavity design. The design of a rotor-stator cavity is critical, since it is directly connected to numerous practical concerns, such as swirl ratio, axial thrust, and moment coefficients. To fulfill industrial expectations, the causes of axial thrust and disk frictional loss are explored.
The core swirl ratio β (the ratio of the fluid’s angular velocity to that of the disk at half the axial gap width) is used to demonstrate the fluid’s dominating tangential motion. The pressure distribution along the rotor may be approximated using the core swirl ratio. The axial thrust operating on the disks may be calculated using the pressure distribution. In addition, the moment coefficient is affected by the tangential velocity profile. Thus, it can be seen that the core swirl ratio is one of the most important parameters in the design of a turbine.
In the case of an enclosed rotor-stator cavity, the core swirl ratio depends exclusively on geometry. Daily and Nece [1] classify the flow regimes into four categories according to core swirl ratio, and turbulent flow regimes (regimes III and IV) are more likely to occur in turbines. Will [2] summarizes some of the usual β values reported in the literature and uses β = 1 / ( 1 + 1 + 5 G ) to forecast the effect of G on β for enclosed rotor-stator cavities.
As for the impact of superposed flow on rotor-stator cavity, Owen [3] provides a flow model to determine the core swirl ratio by solving the Ekman equations for flow regime IV. In his configuration, the flow direction is usually radially outward. The Ekman equations are linearized, simplified versions of the equations of motion. Finally, the functional relation for the core swirl ratio in laminar flow is derived by ignoring the outer cylindrical wall in the motion equation. Owen [3] used the integral technique to solve the nonlinear equations of motion for the disk boundary layer in a turbulent flow. The overall decelerating effect of an outer (stationary) shroud on core rotation may be seen in both measured and calculated data. Will [2] introduced a correction for the stator’s friction factor in order to reflect such effects in the flow model. Please see the review literature [2,3,4] for further information.
When evaluating the above works from a safety perspective, these efforts that can minimize the likelihood of failure might be referred to as “active safety designs.” Despite significant attempts to protect the blade from high strain, the service experience of aero engines suggests that turbine blade fracture is a common problem [5,6]. As a result, substantial attention is devoted to secondary failure induced by initial turbine blade fracture, with an emphasis on failure development. Then, to avoid and manage the spread of failure, many designs were developed that are examples of “passive safety designs.” This article focuses on the influence of turbine blade fracture on the engine, which is a concern of “passive safety designs.”
The outlet boundary of the rotor-stator cavity is axisymmetric in the case of common conditions (that is, the backpressure of the cavity is uniform). However, when turbine blades fracture in actual work, the channel inside the broken blades is exposed to the low-pressure mainstream, which results in discrepant pressure at the outlet of the rotor-stator cavity. Because of the discrepant pressure, a large volume of cold air will rush to the broken blades, reducing the cold air acquired by the normal blades, assuming that the total volume of cold air is constant. A normal blade may potentially fail, due to overlimiting the thermal stress, because the amount of cold air obtained by the normal blade is less than the needed quantity. A continuing development is that the distribution of cold air will become increasingly uneven, and eventually cause a cascade of turbine blades to break [7]. In addition, the probability of turbine blade fracture exceeds 10−7 times per flight hour [8], which is above the probability of hazardous occurrences [9] required by FAA AC 33.75. Therefore, the influence of blade fracture on gas turbines must be demonstrated. In previous publications, the impact of blade fracture on turbines has been investigated [8]. This paper focuses on the influence of the non-axisymmetric boundary caused by turbine blade fracture on a simple rotor-stator cavity.
The axisymmetric boundary and the axisymmetric configuration were used in previous studies of the rotor-stator cavity. Few papers focused on the flow with non-axisymmetric boundary conditions in a rotor-stator system. Below is a synopsis of the relevant research.
Bein et al. [10,11] explored a lubricating oil sealing problem in a narrow rotor-stator cavity. At the outlet of the cavity, there is a high-pressure zone and a low-pressure zone, and a certain pressure lubricating oil is provided at the disk center. The problem has the following characteristics: first, the rotational Reynolds number is very small, R e s = Ω s 2 / ν 1 ; second, the Euler number is sufficiently large, E = ( p 2 p 1 ) / 0.5 ρ Ω 2 b 2 R e s 1 ; third, the gap is narrow, G 1 . Based on these assumptions and characteristics, the velocity and pressure can be obtained by solving the simplified NS equation.
Another type of problem is gas ingress. There has been abundant research on this problem [12,13]. This paper only refers to one significant theoretical model for analyzing this issue: the orifice model. The orifice model was proposed by Owen [14]. It is based on two assumptions: (1) that the flow is from a large reservoir through a small nozzle; and (2) that there is a discontinuous flow across an imaginary surface of the actuator disk. For an axial clearance seal, the orifice model is built on an imaginary “orifice ring.” Through the tiny regions δ A e and δ A i , which add up to the clearance area A c of the seal, egress and ingress occurs across the separate regions of the orifice ring at the same time. For the inviscid equations, mass and energy are considered to be continuous inside the separate stream tubes for egress and ingress, but there is a pressure discontinuity across the orifice ring. The principal “orifice assumptions” are that ( r 2 r 1 ) / r 1 1 and V r , 1 2 V r , 2 2 for egress, and vice versa for ingress. When the pressure distribution inside and outside the orifice ring is known, the velocity distribution of ingress and egress can be solved via the orifice model.
For the rotor-stator cavity studied in this paper, the R e φ 1 and the inertial force term cannot be ignored, so the model of Bein et al. [10] cannot be used. Additionally, because the radial velocity difference between different radii in the cavity is so small, the orifice model cannot be utilized. This paper investigates the effect of turbine blade fracture on the rotor-stator cavity, employing numerical simulation in terms of swirl ratio, mass flow rate distribution, thrust, and moment coefficients.

2. Computational Setup

A typical diagram of a rotor-stator system is shown in paper [15], where the radius of the nozzles is lower than that of the receiver holes and the cold air flow is radially outward. Generally, when a turbine blade fractures, structures such as the mortise are intact. Therefore, this can lead to a low outlet pressure zone in the rotor-stator system. To investigate the effect of turbine blade fracture on a rotor-stator system in isolation, a simple disk cavity model, as shown in the Figure 1, was selected for this paper.
When turbine blades fracture, the backpressure of the outlet in the rotor-stator cavity is no longer axisymmetric, assuming that the backpressure profile of the outlet is a stepped distribution, as shown in Figure 1, P ( θ ) . The high-pressure zone has a pressure of P 2 , and the low-pressure zone has a pressure of P 1 . The low-pressure zone covers the area where the turbine blades are fractured and the high-pressure zone covers the area where the blades are intact. The upstream area is located in the counter-rotating direction of the low-pressure zone, while the downstream area is located in the cis-rotating direction of the low-pressure zone. The upstream border and the downstream border are the interfaces between the high- and low-pressure zones, where the former is located at φ = π (shown as the green line in Figure 1) and the latter is located at φ = π + θ (shown as the red line in Figure 1). The middle of the low-pressure zone means the position in the middle of the low-pressure zone, located at φ = π + θ / 2 (shown as the yellow line in Figure 1). A sector calculation domain cannot be adopted because the boundary conditions are not axisymmetric. The calculation uses a simple 360°domain without a peripheral wall, as shown in Figure 1. The flow enters radially at a low radius and exits at a high radius, and the coordinate is fixed at the rotor’s center. The direction of rotation is pointing to the stator and is perpendicular to the rotor. The gap between the rotor and the stator is 9 mm, and the radius of the rotor and stator is 250 mm.
Poncet et al. [16,17,18] found that when the flow rate is modest, the flow structure in a rotor-stator cavity with radial outflow is similar to that of a closed rotor-stator cavity, while the flow structure becomes Stewartson type when the flow rate is large. Combined with the research of Owen et al. [3,19], the turbulence parameter λ T affects the flow structure in the cavity without pre-swirl. A Stewartson-type structure is found at a low radius and a Batchelor-type structure is found at a high radius in the case of radial outflow. The Stewartson region at a low radius will expand as the λ T increases, while the Batchelor region will shrink until the Stewartson region is filled throughout the entire cavity. Therefore, the λ T is a vital parameter in the case of an intact turbine (in axisymmetric boundary conditions). We also wanted to know what role it plays in the case of a non-intact turbine (in non-axisymmetric boundary conditions). Therefore, the λ T covers the typical working condition of a rotor-stator cavity [3,20] in this paper. Bein et al. [10,11] employed the Euler number, E, to describe the axisymmetry of boundary conditions relative to centrifugal effect. This paper also adopts the E, assuming that the pressure at the mortise (outlet of the cavity) equals the mainstream pressure after the turbine blades fracture. The backflow margin of the film hole is roughly 15% to 20% under design conditions [21]; therefore, the pressure ratio between the high-pressure zone and the low-pressure zone is about 1 P 2 / P 1 1.3 (assuming that the total pressure relative to the turbine rotor is equal to the static pressure). Considering the range of the rotational Reynolds number, the Euler number is 0 E 0.4 in this paper.
Another important parameter is the range of low-pressure, denoted as θ, which describes the number of broken blades. Modern gas turbines typically have 60–70 turbine blades, so a single turbine blade accounts for approximately 0.1 rad. A single turbine blade fracture results in a low-pressure zone range of 0.1 rad, and a succession of N blade fractures results in a low-pressure zone range of 0.1 * N rad, using 0 θ 0.4 in this paper. The cases and parameters involved in this paper are shown in Table 1. When the E or the θ is zero, the turbine blades are intact and the boundary conditions are axisymmetric (corresponding to cases I, J, K, L in Table 1).
The grid used in this paper was meshed by rotating and copying a 2D planar grid to eliminate the effects of a non-axisymmetric grid on calculation (meshed by ICEM19.1). The local grid is shown in Figure 2. Grid setting refers to [22], with 200 radial grids, 140 axial grids, 360 tangential grids, a 0.001 mm thickness of the first layer in the wall boundary layer, and a 1.1 growth rate. The total number of grids surpasses ten million, and Y + < 1 in almost all the disks in the above cases.
The working medium was air with constant density, viscosity, and heat capacity, calculated by commercial software, CFX19.1. A nonslip adiabatic wall condition was defined for all walls. According to the operating conditions in the pre-swirl system, the wall was set to be stationary or rotating, and the rotational speed was set for rotor. The turbulence model of SST and a high-order discrete scheme were employed in this article, according to Poncet et al. [22] and Da [23]. To ignore the effect of heat transfer on the flow, the air was set as an isothermal fluid. Computational convergence was accepted when the residual of the continuity equation converged to R M S < 10 6 .

3. Results and Discussion

3.1. Swirl Ratio

To calculate the axial thrust, it was first necessary to calculate the distribution of the pressure across the disk surface. The distribution of pressure was, in turn, highly dependent on the flow structure in the cavity. Based on previous studies [1], the flow structures were classified by different distributions of the swirl ratio in the axial direction. Therefore, the effect of turbine blade fracture on the swirl ratio was analyzed initially, as set out below.
The distribution of the β in the radial direction for different circumferential positions is shown in Figure 3. As can be seen from that figure, β at ϕ = 0 when blade fracture is exactly the same as when the blades are intact. This means that the position away from the low-pressure zone was not affected by turbine blade fracture in terms of the β. Furthermore, in the low-radius region ( r * < r c * 0.7 ), the β values at different circumferential positions were equal, implying that the effects of turbine blade fracture did not propagate to the low-radius position. In the high-radius region ( r * > r c * 0.7 ), on the other hand, the β was affected to a greater extent. Specifically, in the upstream half of the low-pressure zone (ϕ < π + θ/2), there was a larger β when blades fractured, while in the downstream half of the low-pressure zone (ϕ > π + θ/2), the β was smaller. Clearly, r c * decreases as the degree of boundary asymmetry increases (i.e., the E or the θ increases). It is also worth noting that the β values at the upstream and downstream borders of the low-pressure zone (ϕ = π and ϕ = π + θ) were almost equal in magnitude and opposite in sign, with the absolute value of the β at the downstream border (ϕ = π + θ) being somewhat larger. The reason for the rapid zeroing of the β near the outlet is due to the stepped outlet pressure. Figure 4 shows the velocity vector and contour of the β when the R e φ = 5.12 × 10 5 , E = 0.2, θ = 0.1 at z/s = 0.5, from which it can be seen that the β is roughly symmetrically distributed. However, at locations very close to the outlet (r/b = 1), the small amount of fluid near the upstream and downstream borders is squeezed by the fluid at the lower-radius locations, causing it to be unable to flow out from the outlet of the low-pressure zone and only be able to flow radially at a lower velocity. As a result, the β tends to zero at these extreme locations and, again, these locations are more likely to experience backflow.
Figure 5 shows the distribution of the β in the circumferential direction for different radii. From preliminary calculations, it can be concluded that the rate of change of the β relative to that of an unfractured blade hardly varies with the radius after turbine blade fractures. Furthermore, as the R e φ increases (from case A/I to D/L), the downstream area of the low-pressure zone is progressively more affected than the upstream area.
Another point of interest is that as the R e φ increases, the rate of change of the β decreases (see Figure 6 and Figure 7). For example: at r/b = 0.92, the maximum β increases by 180%, 102%, 62%, and 50% when R e φ = 2.39 × 10 5 ,   5.12 × 10 5 , 1.02 × 10 6 ,     and   2.39 × 10 6 , respectively. This is because the influence of centrifugal forces becomes greater as the R e φ increases and the circumferential imbalance is gradually suppressed by the centrifugal forces, so that the variations in the β decrease. However, as the R e φ continues to increase, the centrifugal forces, which are still gradually increasing, slowly plateau. As a consequence, the rate of change in the β gradually tends to a constant.
The distribution of the β in the circumferential direction for different Euler numbers and different θ numbers is shown in Figure 8 and Figure 9, respectively. The dashed vertical lines in the figures represent the position of the downstream border of the low-pressure zone at different θ numbers. It can be seen from the two figures that the effect of the E and the θ on the β after turbine blade fracture is similar. As the E or the θ increases, the fluctuation of the β increases, especially at high-radius locations. The standard deviation of the β at r/b = 0.44, 0.68, and 0.92 is 0.00381, 0.00576, and 0.01612, respectively, when E = 0.4 (case E in Figure 8). It is worth noting that the change in the maximum and minimum values of E 0 or θ 0 , relative to that of E = 0 or θ = 0 , was not very different, and is slightly larger downstream, which was probably due to the coriolis force of the counter-rotational direction. It is known that when the rotor is stationary, the β will be distributed symmetrically by the middle of the low-pressure zone. However, when the rotor has a rotational speed, the fluid near the upstream and downstream borders is subjected to a reversing tangential coriolis force. As a result, the negative β at downstream borders is smaller, while the positive β at upstream borders is somewhat reduced. Another interesting point in Figure 9 is that at high radii (r/b = 0.92), the peaks of both negative and positive β appearing around the upstream and downstream borders are shifted towards the middle of the low-pressure zone as the θ increases, compared to a symmetrical case.

3.2. Radial Velocity and Mass Flow Rate Distribution

3.2.1. Radial Velocity

When the turbine blades are not fractured, the flow in the rotor-stator cavity is axisymmetric, and therefore the amount of cold air obtained by each turbine blade is equal. However, when the turbine blade fractures and fails, the flow in the rotor-stator cavity is no longer symmetrical, resulting in an uneven velocity as well as a mass flow rate distribution. This leads to problems with uneven cooling of the turbine disk and the turbine blades. This section will focus on the distribution of radial velocity and mass flow rate after a turbine blade fracture.
The distribution of the dimensionless radial velocity in the circumferential direction for different rotational Reynolds numbers is shown in Figure 10. It is clear that as the R e φ increases, the maximum value at the high-radius position (r/b = 0.92) gradually shifts from the middle of the low-pressure zone towards the downstream border. At the same time, the radial velocity distribution becomes more asymmetric and steep, indicating that radial velocity is more sensitive to turbine blade fracture when the R e φ is high. This is also evident from the tangential velocity distribution.
Figure 11 and Figure 12 show the distribution of the dimensionless radial velocity in the circumferential direction for different E and θ numbers. As can be seen from the figures, the E and the θ have a similar effect on the radial velocity distribution. The upstream areas are more affected at the low-radius locations, while the downstream areas are more affected at the high-radius locations. In addition, the high-radius position is more affected than the low-radius position. The standard deviation of the radial velocity at high-radius locations is about one order of magnitude larger than that at low radii. For E = 0.4, the standard deviation of the radial velocity at r/b = 0.44, 0.68, and 0.92 is 0.00447, 0.02628, and 0.14153, respectively. One phenomenon is worth noting in Figure 11. An increase in the E has little effect on the distribution pattern of the radial velocities, which are essentially symmetrically distributed. For r/b = 0.92 and E = 0.04, 0.2, and 0.4, the skewness (which can be used to measure the asymmetry of the distribution) is 5.39, 5.70, and 5.84, respectively, while the kurtosis (which can be used to measure the steepness of the distribution) is 32.07, 35.63, and 37.34, respectively.
As the θ increases, the distribution of radial velocities becomes increasingly asymmetrical, with the downstream gradually taking the lead in the emergence of recirculation zones (as shown in Figure 12). In addition, the minimum values in the downstream areas decrease by a larger amount relative to θ = 0 than those in the upstream areas. (This phenomenon becomes more pronounced as the θ increases, as θ = 0.4 and r/b = 0.92). This is mainly due to the different magnitudes of the radial coriolis force on the upstream and downstream areas. When θ ≠ 0 (turbine blade fracture), the tangential velocities on either side of the low-pressure zone are in opposite directions, and the fluid in the upstream areas is subject to a smaller radial inward coriolis force, while the fluid in the downstream areas is subject to a larger force, so the radial velocity in the downstream areas is smaller and more prone to backflow. This phenomenon is more pronounced at a larger θ, where the difference in tangential velocity between the two sides of the low-pressure zone is greater, and therefore the effect of the coriolis force is more pronounced. This means that when a turbine blade fractures, the downstream area at the high-radius location is more affected, in terms of radial velocity.

3.2.2. Mass Flow Rate

Figure 13 shows the distribution of the mass flow rate in the circumferential direction at different R e φ values with the vertical coordinates indicating the ratio of the actual mass flow rate to the average mass flow rate. It can be seen from the figure that as the R e φ increases, backflow gradually appears in the area near the upstream and downstream borders, and the outflow in the low-pressure area also increases significantly due to mass flow rate conservation. At R e φ = 2.39 × 10 6 , the maximum mass flow rate is 27 times larger than the average value. Meanwhile, as the R e φ increases, the location of the maximum mass flow rate gradually moves from the middle of the low-pressure zone towards the downstream border. Another point of interest is that at the current condition (θ = 0.1), the mass flow rate at the upstream border is always smaller than that at the downstream border, regardless of how the R e φ varies. (This phenomenon is of course more pronounced when the R e φ is larger). This is because the fluid at the low-pressure zone is subject to a higher counter-rotating tangential coriolis force, due to the higher radial velocity than that at other areas, which pushes the low-pressure fluid in the low-pressure zone upstream, eventually leading to a large inverse pressure at the upstream border (as shown in Figure 14, the pressure coefficient is defined by Equation (4)). Therefore, this pressure is more likely to produce backflow upstream.
In addition, as the R e φ increases, the radial inward coriolis force on the downstream region gradually increases (i.e., the tangential velocity relative to the rotor increases, as does the rotational speed, as shown in Figure 5), so backflow also begins to occur. From the above analysis, it is clear that when a particular blade fractures (θ = 0.1), the flow rate of its upstream blade will be reduced or even experience gas intrusion, and the blade is therefore more likely to fracture; continuing to increase the R e φ will also affect the downstream blade.
Figure 15 shows the distribution of the mass flow rate in the circumferential direction for different Euler numbers. As in the analysis above, the variation in the E does not affect the distribution pattern of the mass flow rate; the mass flow rate remains maximum in the middle of the low-pressure zone and minimum at the upstream and downstream borders. As the E increases, the maximum and minimum values increase and decrease, respectively, and the mass flow rate at the upstream border is always smaller than that at the downstream border. The reasons for this are analyzed above.
Figure 16 shows the distribution of the mass flow rate in the circumferential direction for different θ values It is clear from the figure that the maximum values of the mass flow rate occur in the middle of the low-pressure zone when R e φ = 5.12 × 10 5 , regardless of the variation of the θ, and that the maximum values are almost equal. When the θ is small (θ = 0.1), backflow occurs at both the upstream and downstream borders, and the mass flow rate at the upstream border is smaller than that at the downstream border; as the θ increases (θ = 0.2), the backflow zone disappears and the mass flow rate at the upstream and downstream borders is almost equal; upon continuing to increase the θ (θ = 0.4), the backflow zone reappears, but the mass flow rate at the downstream border is smaller than that at the upstream border. It seems that when θ = θ c ( θ c = 0.2 , when C w = 10 , 137 , R e φ = 5.12 × 10 5 , and E = 0.2), the backflow zone disappears and the mass flow rate of the upstream and downstream border are equal; when θ > θ c , the backflow zone gradually appears first at the downstream border; and when θ < θ c , the backflow zone gradually appears first at the upstream border.
The reason for this phenomenon is that when the θ is small, the tangential velocity is small, the radial coriolis forces have little influence, and the pressure dominates, while the inverse pressure gradient near the upstream border is larger (see Figure 14), so the mass flow rate near the upstream border is smaller. As the θ increases, the radial coriolis force gradually increases, but the radial inward coriolis force near the downstream border is larger than at the upstream border, superimposed on the influence of pressure, until θ = θ c , when the upstream and downstream borders are subject to the same radial combined force. In this case, the mass flow rate of upstream and downstream borders is equal. Upon continuing to increase the θ, the effect of the coriolis force is greater, while the Euler number remains the same (the effect of pressure remains the same), so the downstream border is subject to a greater radial inward force, resulting in a smaller mass flow rate at the downstream border than at the upstream border. For the cooling of turbine blades, it is vital to find the θ c . If a blade fractures at a point where θ = θ c , then the distribution of cold air will be as balanced as possible, thus potentially avoiding the serious consequences of gas intrusion and ensuring the safety of the engine.
From the above analysis, it is clear that the θ c will increase as the R e φ and the E increase. Figure 17 shows the distribution of the mass flow rate in the circumferential direction for different θ when C w = 10 , 137 , R e φ = 5.12 × 10 5 , and E = 0.4. The figure shows that for θ = 0.2 , the flow rate at the upstream border is still smaller than that at downstream border, so the value of the θ c should be slightly greater than 0.2.

3.3. Pressure and Thrust

3.3.1. Pressure Coefficient

The axial thrust of a gas turbine is critical to the performance and safety of the engine. However, to the authors’ knowledge, there are no studies to date on the change in turbine disk axial thrust after turbine blade fracture. This section investigates the distribution of pressure and then analyses the variation of axial thrust. In previous studies, Will et al. [2] used Equation (1) to evaluate the pressure distribution along the radius of the disk in a rotor-stator cavity with through-flow for an incompressible, steady flow.
p r = ρ β C q r 2 Ω 2 r + ρ Q 2 4 π 2 s 2 r 3
Based on Equation (1), the pressure along the radius of the disk can be calculated with Equation (2) by Hu et al. [24].
p ( r ) = p b + b r ρ β C q r 2 Ω 2 r d r + ρ Q 2 8 π 2 s 2 ( 1 b 2 1 r 2 )
Dimensionless pressure is defined as follows:
P * = p ρ Ω 2 b 2
The pressure coefficient is defined as:
C p = P * ( r b = 1 ) P * ( r b )
Figure 18 shows the pressure coefficient distribution for the different circumferential positions of cases B, D, F, and H, numbered a, b, c and d, respectively. As can be seen from Figure 18a, the pressure coefficients obtained by Equation (2) are in very good agreement with the SST model, indicating that the SST model is suitable for this problem. Throughout these four plots, it can be seen that the pressure coefficients at the low-radius locations are almost unaffected, in line with the previous analysis. In contrast, at high radii, the pressure coefficients are generally larger at the middle of the low-pressure zone (ϕ = π + θ/2) than at other locations (except in case H), and larger at the downstream border (ϕ = π + θ) than at the upstream border (ϕ = π). Comparing the plots of cases B and D, the pressure coefficients of upstream and downstream borders tend to be equal at high-radius locations as the R e φ increases. Comparing the two plots of cases B and F, a change in the Euler number does not significantly affect the distribution pattern of the pressure coefficients; only the magnitude of the values changes. Comparing the graphs of cases B and H, the pressure coefficient at the downstream border even exceeds that at the middle of the low-pressure zone as the θ increases. This corresponds to the fact that the downstream border is the first to experience backflow as the θ increases (see Figure 16).
Figure 19 shows the distribution of pressure coefficients at different R e φ values, where a, b, c and d represent different circumferential positions. As the R e φ increases, the pressure coefficients at each location do not vary much, except for a slight increase near the upstream border at the high-radius location. The pressure coefficient decreases as the R e φ increases. However, by continuing to increase the R e φ , the pressure coefficient at high radii increases (see the enlargement of Figure 19a). Furthermore, there is a sudden drop in the pressure coefficient at the downstream border (ϕ = π + θ) near the outlet. This is mainly because the fluid in the low-pressure zone is carried upstream and the downstream is filled with high-pressure fluid (similar to Figure 14).
Figure 20 shows the distribution of the pressure coefficients for different Euler numbers, where a, b, c and d represent different circumferential positions. It is clear that an increase in the Euler number does not affect the pressure at low radii. As for the higher radii, the pressure coefficients in and around the low-pressure zone increase significantly as the Euler number increases.
The distribution of pressure coefficients for different θ values is shown in Figure 21, where a, b, c and d represent different circumferential positions. Unlike the E and the R e φ , variations in the θ can affect the pressure coefficient even at low-radius locations, resulting in an increase in the pressure coefficient in the low-pressure region and its immediate vicinity.

3.3.2. Thrust Coefficient

Figure 22 shows the thrust coefficients of the rotor and stator for different cases, where a, b and c show the impact of the R e φ , E and θ respectively. The thrust coefficient characterizes the ratio of the axial force on the disk to the centrifugal force, which is defined as in [25]:
C F = ( p b p ) d s ρ Ω 2 b 4 = a b 2 π ( p b p ) r d r ρ Ω 2 b 4
By definition, C F and C p are related as follows:
C F = ( p b p ) d s ρ Ω 2 b 4 = 1 b 2 ( p b p ) d s ρ Ω 2 b 2 = 1 b 2 C p d s
It is clear that the thrust coefficients of the rotor and the stator are almost equal. E = 0 and θ = 0 mean the blades are intact, while E = 0.2 and θ = 0.1 means a blade is fractured. From Figure 22a, it can be seen that the disk thrust coefficient of the stator increases by 1.7%, 3.3%, 3.8%, and 0.8%, when the R e φ = 2.39 × 10 5 ,   5.12 × 10 5 ,   1.02 × 10 6 ,   and   2.39 × 10 6 , respectively, because of blade fracture. This is mainly due to the fracture of the turbine blades, resulting in a low-pressure area on the disk surface. In addition, the increased ratio first increases and then decreases as the R e φ increase. When the blades are already fractured, the thrust coefficient of the disks decreases significantly as the R e φ increases but as the R e φ continues to increase, the thrust coefficient increases. This is the result of the combined effect of centrifugal force and the low-pressure zone. When the R e φ is small, increasing the R e φ will increase the effect of centrifugal force, making the pressure coefficient gradually converge to zero across the disk (refer to Figure 19), so the thrust coefficient decreases significantly. However, by continuing to increase the R e φ , the pressure coefficient at high radii increases (see the enlargement of Figure 19a). Therefore, continuing to increase the R e φ will increase the thrust coefficient slightly, according to Equation (6). This also means that the thrust coefficient is more sensitive to the R e φ whether blades fracture or not, especially when the R e φ is small. Figure 22b shows that an increase in the Euler number increases the thrust coefficient. Figure 22c shows that an increase in the θ also increases the thrust coefficient, and the effect of the θ on the thrust coefficient is approximately linear.

3.4. Moment Coefficient

Figure 23 shows the moment coefficients of the rotor and the stator on the rotor shaft z for different cases, where a, b and c show the impact of the R e φ , E and θ respectively. The moment coefficient is defined by Equation (7) [25]:
C M = 2 · | M | ρ Ω 2 b 5
From Figure 23, it can be seen that the moment coefficient of the rotor is much greater than that of the stator. In addition, the E (turbine blade fracture) and the θ have almost no effect on the moment coefficient. Figure 23a shows that as the R e φ increases, the moment coefficient of the rotor decreases while that of the stator increases, whether or not blades fracture. This is mainly because as the R e φ increases, the swirl ratio increases, and therefore the moment coefficient of the rotor decreases (according to Han et al. [26]). However, an increase in the swirl ratio increases the relative velocity of the fluid to the stator, and therefore the moment coefficient of the stator increases.
After the turbine blade fracture, a low-pressure zone appears and the pressure distribution in the rotor-stator cavity is no longer symmetrical, so the moment coefficient on the radial direction is not zero. This means that the turbo disk tends to roll over on its side. The moment coefficients for the rotor and the stator on φ = π / 2 + θ / 2 are shown in Figure 24, where a, b and c show the impact of the R e φ , E and θ respectively. It is equal in magnitude to that on the rotating shaft. The moment coefficient in this direction is the largest of all radial directions. The moment coefficients of the stator and the rotor are essentially the same. Figure 24a shows that as the R e φ increases, the moment coefficient decreases. This is because the pressure coefficient gradually converges to zero along the radius as the R e φ increases. From Figure 24b,c, it can be seen that increasing the E and the θ increases the moment coefficient.

4. Conclusions

Turbine blade fracture can lead to engine safety problems, so this paper investigated the effect of turbine blade fracture on the flow in a simple rotor-stator system without a shroud, employing numerical simulation. The results of the numerical simulations were compared with relevant empirical correlations. Specifically, this paper studied the effects of control parameters, such as the rotational Reynolds number, the Euler number, and the range of low-pressure zones on the velocity and pressure fields, as well as the thrust coefficient, the moment coefficient, and heat transfer in a simple rotor-stator cavity. Within the scope of this paper’s research, the following conclusions can be drawn.
  • For the swirl ratio, the effects of the rotational Reynolds number, the Euler number, and the θ are similar. In addition, although the downstream region is more affected than the upstream region, an increase in the Euler number and the θ increases the swirl ratio variation, while an increase in the rotational Reynolds number decreases the swirl ratio variation.
  • Increases in the rotational Reynolds number, the Euler number, and the θ all lead to a more uneven distribution of the flow rate. Furthermore, regardless of the rotational Reynolds number and the Euler number, the flow rate at the upstream border is always smaller than at the downstream border, but an increase in the θ may lead to a more balanced flow rate distribution (there is a critical θ c that makes the flow rate distribution most balanced; θ c 0.2 when C w = 10 , 137 , R e φ = 5.12 × 10 5 , and 0.2 E 0.4 ).
  • Turbine blade fracture causes an increase in the thrust coefficient and is more pronounced at smaller rotational Reynolds numbers. The increase in the thrust coefficient does not exceed 4% when E = 0.2 ,   θ = 0.1 , as discussed in this paper.
  • Changes in the rotational Reynolds number, the Euler number, and the θ have almost no effect on the moment coefficient about the axis of rotation but have a more significant effect on the moment coefficient about the radial direction. The latter will decrease as the rotational Reynolds number increases and increase as the Euler number and the θ increase.

Author Contributions

Funding acquisition, T.Q. and P.L.; writing–original draft, G.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Project of National Science Foundation of China, grant number 61890923. The APC was funded by Project of National Science Foundation of China No. 61890923.

Data Availability Statement

Not applicable.

Acknowledgments

The authors express their sincere gratitude to the support of the National Fund.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ainlet radius, m
Bradius of rotor and stator, m
Saxial spacing between rotor and stator, m
A c effective sealing area, m 2
δ A e , δ A i area of a small orifice where air egresses/ingresses, m 2
C w dimensionless mass flow rate, m / μ b
EEuler number, ( P 2 P 1 ) / 0.5   ρ Ω 2 b 2
Ggap ratio, s/b
m ˙ mass flow rate, kg/s
Nnumber of blades
P 1 , P 2 outlet pressure of rotor-stator cavity, Pa
P * dimensionless pressure difference, ( P 2 p ) / 0.5 ρ Ω 2 r 2
P ( θ ) pressure profile of outlet, Pa
Qvolume flow rate, m3/s
r * dimensionless radius, r/b
R e s rotational Reynolds number based on s, Ω s 2 / ν
R e r radial to rotational Reynolds number, R e r = C w / 2 π G R e φ
R e φ rotational Reynolds number based on b, Ω b 2 / ν
V r , V φ radial and tangential velocity, m/s
Greek
βswirl ratio, V φ / Ω r
θ the range of low-pressure area, rad
Ωrotating velocity, rad/s
νkinematic viscosity, m2/s
ρdensity, kg/m3
λ T turbulent flow parameter, C w / R e φ 0.8
Subscript
r, z, ϕradial, axial, and tangential coordinates, m/m/rad
ccritical
1, 2different location

References

  1. Daily, J.W.; Nece, R.E. Chamber dimension effects on induced flow and frictional resistance of enclosed rotating disks. J. Basic Eng. Trans. ASME 1960, 82, 217–230. [Google Scholar] [CrossRef]
  2. Will, B.-C. Theoretical, Numerical and Experimental Investigation of the Flow in Rotor-Stator Cavities with Application to a Centrifugal Pump; Universitätsbibliothek Duisburg-Essen: Essen, Germany, 2011. [Google Scholar]
  3. Owen, J.M.; Rogers, R.H. Flow and Heat Transfer in Rotating Disc Systems, Vol.1: Rotor-Stator Systems; Research Studies Press LTD.: Taunton, Somerset, UK, 1989; Volume 90, p. 278. [Google Scholar]
  4. Zhang, F.; Adu-Poku, K.A.; Hu, B.; Appiah, D.; Chen, K. Flow theory in the side chambers of the radial pumps: A review. Phys. Fluids 2020, 32, 041301. [Google Scholar] [CrossRef]
  5. Bureau, A.T.S. Power Plant Failures in Turbofan-Powered Aircraft 2008 to 2012; AR-2013-002; Australian Transport Safety Bureau: Canberra, Australian, 2014; pp. 1–34.
  6. Meher-Homji, C.B. Blading vibration and failures in gas turbines: Part A—Blading dynamics and the operating environment. In Proceedings of the ASME 1995 International Gas Turbine and Aeroengine Congress and Exposition, Houston, TX, USA, 5–8 June 1995; pp. 1–11. [Google Scholar]
  7. Hall, J. Reply Refer to A-98-8; National Transportation Safety Board: Washington, DC, USA, 1998; pp. 1–6.
  8. Zhang, S.; Ding, S.; Qiu, T. Aerodynamic Performance Investigation of Turbine in the Event of One Blade Primary Fracture Failure. In Proceedings of the Turbo Expo: Power for Land, Sea, and Air, London, UK, 21–25 September 2020. [Google Scholar]
  9. Federal Aviation Administration. Airworthiness Standards: Aircraft Engines: 33.75 Safety Analysis; Federal Aviation Administration: Washington, DC, USA, 2007.
  10. Bein, M.; Shavit, A.; Solan, A. Nonaxisymmetric flow in the narrow gap between a rotating and a stationary disk. J. Fluids Eng. 1976, 98, 217–223. [Google Scholar] [CrossRef]
  11. Bein, M.; Shavit, A. Nonaxisymmetric Flow in the Narrow Gap Between a Rotating and a Stationary Disk With an Eccentric Source. J. Fluid Eng. 1977, 99, 418–421. [Google Scholar] [CrossRef]
  12. Owen, J.M. Prediction of Ingestion Through Turbine Rim Seals-Part II: Externally Induced and Combined Ingress. J. Turbomach.-Trans. Asme 2011, 133, 031006. [Google Scholar] [CrossRef]
  13. Scobie, J.A.; Sangan, C.M.; Owen, J.M.; Lock, G.D. Review of Ingress in Gas Turbines. J. Eng. Gas Turbines Power-Trans. Asme 2016, 138, 120801. [Google Scholar] [CrossRef]
  14. Owen, J.M. Prediction of Ingestion Through Turbine Rim Seals-Part I: Rotationally Induced Ingress. J. Turbomach.-Trans. Asme 2011, 133, 031005. [Google Scholar] [CrossRef]
  15. Kakade, V.U.; Lock, G.D.; Wilson, M.; Owen, J.M.; Mayhew, J.E. Effect of Radial Location of Nozzles on Heat Transfer in Pre-Swirl Cooling Systems. Proc. Asme Turbo Expo 2009, 3, 1051–1060. [Google Scholar]
  16. Poncet, S.; Chauve, M.P.; Le Gal, P. Turbulent rotating disk flow with inward throughflow. J. Fluid Mech. 2005, 522, 253–262. [Google Scholar] [CrossRef]
  17. Poncet, S.; Chauve, M.P.; Schiestel, R. Batchelor versus Stewartson flow structures in a rotor-stator cavity with throughflow. Phys. Fluids 2005, 17, 075110. [Google Scholar] [CrossRef]
  18. Poncet, S.; Schiestel, R.; Chauve, M.P. Centrifugal flow in a rotor-stator cavity. J. Fluids Eng.-Trans. Asme 2005, 127, 787–794. [Google Scholar] [CrossRef]
  19. Owen, J. An approximate solution for the flow between a rotating and a stationary disk. In Proceedings of the 33rd International Gas Turbine and Aeroengine Congress and Exhibition, Amsterdam, The Netherlands, 5–8 June 1988; pp. 323–332. [Google Scholar]
  20. Yan, Y.; Gord, M.F.; Lock, G.D.; Wilson, M.; Owen, J.M. Fluid dynamics of a pre-swirl rotor-stator system. J. Turbomach.-Trans. Asme 2003, 125, 641–647. [Google Scholar] [CrossRef]
  21. Halila, E.; Lenahan, D.; Thomas, T. Energy Efficient Engine High Pressure Turbine Test Hardware Detailed Design Report; NASA CR-167955; National Aeronautics and Space Adiministration: Washington, DC, USA, 1982; pp. 1–186.
  22. Poncet, S.; Schiestel, R. Numerical modeling of heat transfer and fluid flow in rotor-stator cavities with throughflow. Int. J. Heat Mass. Tran. 2007, 50, 1528–1544. [Google Scholar] [CrossRef] [Green Version]
  23. Da Soghe, R.; Innocenti, L.; Andreini, A.; Poncet, S. Numerical Benchmark of Turbulence Modeling in Gas Turbine Rotor-Stator System. In Proceedings of the ASME Turbo Expo 2010: Turbomachinery: Axial Flow Fan and Compressor Aerodynamics Design Methods, and Cfd Modeling for Turbomachinery, Glasgow, UK, 14–18 June 2010; Volume 7, pp. 771–783. [Google Scholar]
  24. Hu, B.; Brillert, D.; Dohmen, H.J.; Benra, F.-K. Investigation on thrust and moment coefficients of a centrifugal turbomachine. Int. J. Turbomach. Propuls. Power 2018, 3, 9. [Google Scholar] [CrossRef] [Green Version]
  25. Hu, B.; Li, X.S.; Fu, Y.X.; Gu, C.W.; Ren, X.F.; Lu, J.X. Axial Thrust, Disk Frictional Losses, and Heat Transfer in a Gas Turbine Disk Cavity. Energies 2019, 12, 2917. [Google Scholar] [CrossRef] [Green Version]
  26. Hongbiao, H.; Shanqun, G.; Jishen, L.; Yongzhen, Z. Exploring Fluid Resistance of Disk Rotor Based on Boundary Layer Theory. Mech. Sci. Technol. Aerosp. Eng. 2015, 34, 1621–1625. [Google Scholar]
Figure 1. Computational model used in this paper.
Figure 1. Computational model used in this paper.
Aerospace 09 00106 g001
Figure 2. Local mesh.
Figure 2. Local mesh.
Aerospace 09 00106 g002
Figure 3. The β in different circumferential positions when the R e φ = 5.12 × 10 5 at z/s = 0.5.
Figure 3. The β in different circumferential positions when the R e φ = 5.12 × 10 5 at z/s = 0.5.
Aerospace 09 00106 g003
Figure 4. Velocity vector and contour of the β when the R e φ = 5.12 × 10 5 , E = 0.2, θ = 0.1 at z/s = 0.5.
Figure 4. Velocity vector and contour of the β when the R e φ = 5.12 × 10 5 , E = 0.2, θ = 0.1 at z/s = 0.5.
Aerospace 09 00106 g004
Figure 5. β curves with different radius at z/s = 0.5.
Figure 5. β curves with different radius at z/s = 0.5.
Aerospace 09 00106 g005
Figure 6. β curves with different R e φ for C w = 10 , 137 , E = 0.2 , θ = 0.1 at z/s = 0.5.
Figure 6. β curves with different R e φ for C w = 10 , 137 , E = 0.2 , θ = 0.1 at z/s = 0.5.
Aerospace 09 00106 g006
Figure 7. Contour of the β when E = 0.2, θ = 0.1 at z/s = 0.5.
Figure 7. Contour of the β when E = 0.2, θ = 0.1 at z/s = 0.5.
Aerospace 09 00106 g007
Figure 8. β curves for C w = 10137 , R e φ = 5.12 × 10 5 , θ = 0.1 at z/s = 0.5.
Figure 8. β curves for C w = 10137 , R e φ = 5.12 × 10 5 , θ = 0.1 at z/s = 0.5.
Aerospace 09 00106 g008
Figure 9. β curves for C w = 10 , 137 , R e φ = 5.12 × 10 5 , E = 0.2 at z/s = 0.5.
Figure 9. β curves for C w = 10 , 137 , R e φ = 5.12 × 10 5 , E = 0.2 at z/s = 0.5.
Aerospace 09 00106 g009
Figure 10. Dimensionless radial velocity for C w = 10 , 137 , E = 0.2 , θ = 0.1 at z/s = 0.5.
Figure 10. Dimensionless radial velocity for C w = 10 , 137 , E = 0.2 , θ = 0.1 at z/s = 0.5.
Aerospace 09 00106 g010
Figure 11. Dimensionless radial velocity for C w = 10 , 137 , R e φ = 5.12 × 10 5 , θ = 0.1 at z/s = 0.5.
Figure 11. Dimensionless radial velocity for C w = 10 , 137 , R e φ = 5.12 × 10 5 , θ = 0.1 at z/s = 0.5.
Aerospace 09 00106 g011
Figure 12. Dimensionless radial velocity for C w = 10 , 137 , R e φ = 5.12 × 10 5 , E = 0.2 at z/s = 0.5.
Figure 12. Dimensionless radial velocity for C w = 10 , 137 , R e φ = 5.12 × 10 5 , E = 0.2 at z/s = 0.5.
Aerospace 09 00106 g012
Figure 13. Distribution of mass flow rate for C w = 10 , 137 , E = 0.2 , θ = 0.1 at r/b = 1.
Figure 13. Distribution of mass flow rate for C w = 10 , 137 , E = 0.2 , θ = 0.1 at r/b = 1.
Aerospace 09 00106 g013
Figure 14. Pressure coefficient of case D at plane z/s = 0.5.
Figure 14. Pressure coefficient of case D at plane z/s = 0.5.
Aerospace 09 00106 g014
Figure 15. Distribution of mass flow rate for C w = 10 , 137 , R e φ = 5.12 × 10 5 , θ = 0.1 at r/b = 1.
Figure 15. Distribution of mass flow rate for C w = 10 , 137 , R e φ = 5.12 × 10 5 , θ = 0.1 at r/b = 1.
Aerospace 09 00106 g015
Figure 16. Distribution of mass flow rate for C w = 10 , 137 , R e φ = 5.12 × 10 5 , E = 0.2 at r/b = 1.
Figure 16. Distribution of mass flow rate for C w = 10 , 137 , R e φ = 5.12 × 10 5 , E = 0.2 at r/b = 1.
Aerospace 09 00106 g016
Figure 17. Distribution of mass flow rate for C w = 10 , 137 , R e φ = 5.12 × 10 5 , E = 0.4 at r/b = 1.
Figure 17. Distribution of mass flow rate for C w = 10 , 137 , R e φ = 5.12 × 10 5 , E = 0.4 at r/b = 1.
Aerospace 09 00106 g017
Figure 18. Pressure coefficient of different cases at plane z/s = 0.5: (a) case B, (b) case D, (c) case F, (d) case H.
Figure 18. Pressure coefficient of different cases at plane z/s = 0.5: (a) case B, (b) case D, (c) case F, (d) case H.
Aerospace 09 00106 g018
Figure 19. Pressure coefficient for C w = 10 , 137 , E = 0.2 , θ = 0.1 at plane z/s = 0.5: (a) ϕ = 0, (b) ϕ = π, (c) ϕ = π + θ/2, (d) ϕ = π + θ.
Figure 19. Pressure coefficient for C w = 10 , 137 , E = 0.2 , θ = 0.1 at plane z/s = 0.5: (a) ϕ = 0, (b) ϕ = π, (c) ϕ = π + θ/2, (d) ϕ = π + θ.
Aerospace 09 00106 g019
Figure 20. Pressure coefficient for C w = 10 , 137 , R e φ = 5.12 × 10 5 , θ = 0.1 , at plane z/s = 0.5: (a) ϕ = 0, (b) ϕ = π, (c) ϕ = π + θ/2, (d) ϕ = π + θ.
Figure 20. Pressure coefficient for C w = 10 , 137 , R e φ = 5.12 × 10 5 , θ = 0.1 , at plane z/s = 0.5: (a) ϕ = 0, (b) ϕ = π, (c) ϕ = π + θ/2, (d) ϕ = π + θ.
Aerospace 09 00106 g020
Figure 21. Pressure coefficient for C w = 10 , 137 , R e φ = 5.12 × 10 5 , E = 0.2 , at plane z/s = 0.5: (a) ϕ = 0, (b) ϕ = π, (c) ϕ = π + θ/2, (d) ϕ = π + θ.
Figure 21. Pressure coefficient for C w = 10 , 137 , R e φ = 5.12 × 10 5 , E = 0.2 , at plane z/s = 0.5: (a) ϕ = 0, (b) ϕ = π, (c) ϕ = π + θ/2, (d) ϕ = π + θ.
Aerospace 09 00106 g021
Figure 22. Thrust coefficient of rotor and stator: (a) E = 0.2 , θ = 0.1 ; (b) R e φ = 5.12 × 10 5 , θ = 0.1 ; (c) R e φ = 5.12 × 10 5 , E = 0.2 .
Figure 22. Thrust coefficient of rotor and stator: (a) E = 0.2 , θ = 0.1 ; (b) R e φ = 5.12 × 10 5 , θ = 0.1 ; (c) R e φ = 5.12 × 10 5 , E = 0.2 .
Aerospace 09 00106 g022
Figure 23. Moment coefficient of rotor and stator on the rotating shaft (z): (a) E = 0.2 , θ = 0.1 ; (b) R e φ = 5.12 × 10 5 , θ = 0.1 ; (c) R e φ = 5.12 × 10 5 , E = 0.2 .
Figure 23. Moment coefficient of rotor and stator on the rotating shaft (z): (a) E = 0.2 , θ = 0.1 ; (b) R e φ = 5.12 × 10 5 , θ = 0.1 ; (c) R e φ = 5.12 × 10 5 , E = 0.2 .
Aerospace 09 00106 g023
Figure 24. Moment coefficient of rotor and stator on the φ = π / 2 + θ / 2 : (a) E = 0.2 , θ = 0.1 ; (b) R e φ = 5.12 × 10 5 , θ = 0.1 ; (c) R e φ = 5.12 × 10 5 , E = 0.2 .
Figure 24. Moment coefficient of rotor and stator on the φ = π / 2 + θ / 2 : (a) E = 0.2 , θ = 0.1 ; (b) R e φ = 5.12 × 10 5 , θ = 0.1 ; (c) R e φ = 5.12 × 10 5 , E = 0.2 .
Aerospace 09 00106 g024
Table 1. The range of parameters.
Table 1. The range of parameters.
CaseG C w Θ R e φ λ T E
A0.03613,1070.1 2.39 × 10 5 0.650.2
B0.03613,1070.1 5.12 × 10 5 0.350.2
C0.03613,1070.1 1.02 × 10 6 0.200.2
D0.03613,1070.1 2.39 × 10 6 0.100.2
E0.03613,1070.1 5.12 × 10 5 0.350.4
F0.03613,1070.1 5.12 × 10 5 0.350.04
G0.03613,1070.2 5.12 × 10 5 0.350.2
H0.03613,1070.4 5.12 × 10 5 0.350.2
I0.03613,1070 2.39 × 10 5 0.650
J0.03613,1070 5.12 × 10 5 0.350
K0.03613,1070 1.02 × 10 6 0.200
L0.03613,1070 2.39 × 10 6 0.100
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhao, G.; Qiu, T.; Liu, P. Influence of Blade Fracture on the Flow of Rotor-Stator Systems with Centrifugal Superposed Flow. Aerospace 2022, 9, 106. https://doi.org/10.3390/aerospace9020106

AMA Style

Zhao G, Qiu T, Liu P. Influence of Blade Fracture on the Flow of Rotor-Stator Systems with Centrifugal Superposed Flow. Aerospace. 2022; 9(2):106. https://doi.org/10.3390/aerospace9020106

Chicago/Turabian Style

Zhao, Gang, Tian Qiu, and Peng Liu. 2022. "Influence of Blade Fracture on the Flow of Rotor-Stator Systems with Centrifugal Superposed Flow" Aerospace 9, no. 2: 106. https://doi.org/10.3390/aerospace9020106

APA Style

Zhao, G., Qiu, T., & Liu, P. (2022). Influence of Blade Fracture on the Flow of Rotor-Stator Systems with Centrifugal Superposed Flow. Aerospace, 9(2), 106. https://doi.org/10.3390/aerospace9020106

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