Inﬂuence of Blade Fracture on the Flow of Rotor-Stator Systems with Centrifugal Superposed Flow

: 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 ﬂow is studied in this paper. The Euler number E , the rotational Reynolds number Re ϕ , and the low-pressure zone range θ are investigated and, for the ﬁrst 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 Re ϕ is large. As for the distribution of the mass ﬂow rate, there may be a critical θ c at which the other blades are least affected. The θ c would increase as the Re ϕ or the E increase, and the θ c ∼ = 0.2 when C w = 10,137, Re ϕ = 5.12 × 10 5 , and 0.2 ≤ E ≤ 0.4. In addition, the thrust coefﬁcient increases as the E or the θ increases, and the increase in the thrust coefﬁcient does not exceed 4% when the E = 0.2 and the θ = 0.1 in this paper. However, the moment coefﬁcient on the rotating shaft is almost independent of the E and the θ . An increase in the Re ϕ will reduce the effect of turbine blade fracture on the thrust and moment coefﬁcients, when the Re ϕ is small.


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 + 5G 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 rotorstator 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, Re s = Ωs 2 /ν 1; second, the Euler number is sufficiently large, E = (p 2 − p 1 )/0.5ρΩ 2 b 2 Re 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 2 r,1 V 2 r,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 Re ϕ 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.

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.
problem [12,13].This paper only refers to one significant theoretical model for ana this issue: the orifice model.The orifice model was proposed by Owen [14].It is ba two assumptions: (1) that the flow is from a large reservoir through a small nozz (2) that there is a discontinuous flow across an imaginary surface of the actuator di an axial clearance seal, the orifice model is built on an imaginary "orifice ring."Th the tiny regions and , which add up to the clearance area of the seal, and ingress occurs across the separate regions of the orifice ring at the same time.inviscid equations, mass and energy are considered to be continuous inside the se stream tubes for egress and ingress, but there is a pressure discontinuity across the ring.The principal "orifice assumptions" are that ( − ) ⁄ ≪ 1 and , ≪ egress, and vice versa for ingress.When the pressure distribution inside and outs orifice ring is known, the velocity distribution of ingress and egress can be solved orifice model.
For the rotor-stator cavity studied in this paper, the ≫ 1 and the inertia term cannot be ignored, so the model of Bein et al. [10] cannot be used.Additiona cause the radial velocity difference between different radii in the cavity is so sm orifice model cannot be utilized.This paper investigates the effect of turbine blade f on the rotor-stator cavity, employing numerical simulation in terms of swirl ratio flow rate distribution, thrust, and moment coefficients.

Computational Setup
A typical diagram of a rotor-stator system is shown in paper [15], where the of the nozzles is lower than that of the receiver holes and the cold air flow is r outward.Generally, when a turbine blade fractures, structures such as the mort intact.Therefore, this can lead to a low outlet pressure zone in the rotor-stator syst investigate the effect of turbine blade fracture on a rotor-stator system in isolation ple 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 grid used in this paper was meshed by rotating and copying a 2D planar g eliminate the effects of a non-axisymmetric grid on calculation (meshed by ICEM The local grid is shown in Figure 2. Grid setting refers to [22], with 200 radial grid axial grids, 360 tangential grids, a 0.001 mm thickness of the first layer in the wall b ary layer, and a 1.1 growth rate.The total number of grids surpasses ten million, and 1 in almost all the disks in the above cases.
The working medium was air with constant density, viscosity, and heat capacit culated by commercial software, CFX19.1.A nonslip adiabatic wall condition was de for all walls.According to the operating conditions in the pre-swirl system, the wa set to be stationary or rotating, and the rotational speed was set for rotor.The turbu model of SST and a high-order discrete scheme were employed in this article, acco to Poncet et al. [22] and Da [23].To ignore the effect of heat transfer on the flow, t was set as an isothermal fluid.Computational convergence was accepted when the ual of the continuity equation converged to < 10 .

Swirl Ratio
To calculate the axial thrust, it was first necessary to calculate the distribution pressure across the disk surface.The distribution of pressure was, in turn, highly de ent on the flow structure in the cavity.Based on previous studies [1], the flow stru were classified by different distributions of the swirl ratio in the axial direction.Ther the effect of turbine blade fracture on the swirl ratio was analyzed initially, as set o low.
The distribution of the β in the radial direction for different circumferential pos is shown in Figure 3.As can be seen from that figure, β at ϕ = 0 when blade fract 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 RMS < 10 −6 .

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 Re ϕ = 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.
Aerospace 2022, 9, x FOR PEER REVIEW 6 of 25 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 ( * < * ≅ 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 ( * > * ≅ 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, * 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 = 5.12 × 10 , 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.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 ( * < * ≅ 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 ( * > * ≅ 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, * 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 = 5.12 × 10 , 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 Re ϕ increases (from case A/I to D/L), the downstream area of the low-pressure zone is progressively more affected than the upstream area.
erospace 2022, 9, x FOR PEER REVIEW 7 of 2 Figure 5 shows the distribution of the β in the circumferential direction for differen 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 blad fractures.Furthermore, as the increases (from case A/I to D/L), the downstream are of the low-pressure zone is progressively more affected than the upstream area.
-  Another point of interest is that as the increases, the rate of change of the β de creases (see Figures 6 and 7).For example: at r/b = 0.92, the maximum β increases by 180% 102%, 62%, and 50% when = 2.39 × 10 , 5.12 × 10 , 1.02 × 10 , and 2.39 × 10 , re spectively.This is because the influence of centrifugal forces becomes greater as the increases and the circumferential imbalance is gradually suppressed by the centrifuga forces, so that the variations in the β decrease.However, as the continues to increase the centrifugal forces, which are still gradually increasing, slowly plateau.As a conse quence, the rate of change in the β gradually tends to a constant.Another point of interest is that as the Re ϕ increases, the rate of change of the β decreases (see Figures 6 and 7).For example: at r/b = 0.92, the maximum β increases by 180%, 102%, 62%, and 50% when Re ϕ = 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 Re ϕ increases and the circumferential imbalance is gradually suppressed by the centrifugal forces, so that the variations in the β decrease.However, as the Re ϕ 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 nu  The distribution of the β in the circumferential direction for different Euler numbers and different θ numbers is shown in Figures 8 and 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 The distribution of the β in the circumferential direction for different Euler numbers and different θ numbers is shown in Figures 8 and 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.
Aerospace 2022, 9, x FOR PEER REVIEW 9 0.4 (case E in Figure 8).It is worth noting that the change in the maximum and minim values of ≠ 0 or ≠ 0, relative to that of = 0 or = 0, was not very different, is slightly larger downstream, which was probably due to the coriolis force of the coun rotational direction.It is known that when the rotor is stationary, the β will be distrib symmetrically by the middle of the low-pressure zone.However, when the rotor h rotational speed, the fluid near the upstream and downstream borders is subjected reversing tangential coriolis force.As a result, the negative β at downstream borde smaller, while the positive β at upstream borders is somewhat reduced.Another inte ing point in Figure 9 is that at high radii (r/b = 0.92), the peaks of both negative and pos β appearing around the upstream and downstream borders are shifted towards the m dle of the low-pressure zone as the θ increases, compared to a symmetrical case.

Radial Velocity
When the turbine blades are not fractured, the flow in the rotor-stator cavity is isymmetric, and therefore the amount of cold air obtained by each turbine blade is eq However, when the turbine blade fractures and fails, the flow in the rotor-stator cavi no longer symmetrical, resulting in an uneven velocity as well as a mass flow rate di The distribution of the dimensionless radial velocity in the circumferential d for different rotational Reynolds numbers is shown in Figure 10.It is clear that as t increases, the maximum value at the high-radius position (r/b = 0.92) gradually shi the middle of the low-pressure zone towards the downstream border.At the sam the radial velocity distribution becomes more asymmetric and steep, indicating tha velocity is more sensitive to turbine blade fracture when the is high.This is a dent from the tangential velocity distribution.

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 Re ϕ 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 Re ϕ is high.This is also evident from the tangential velocity distribution.
Aerospace 2022, 9, x FOR PEER REVIEW Figures 11 and 12 show the distribution of the dimensionless radial velocity circumferential direction for different E and θ numbers.As can be seen from the f the E and the θ have a similar effect on the radial velocity distribution.The upstream are more affected at the low-radius locations, while the downstream areas are m fected at the high-radius locations.In addition, the high-radius position is more a than the low-radius position.The standard deviation of the radial velocity at high locations is about one order of magnitude larger than that at low radii.For E = 0 standard deviation of the radial velocity at r/b = 0.44, 0.68, and 0.92 is 0.00447, 0.026 0.14153, respectively.One phenomenon is worth noting in Figure 11.An increase i has little effect on the distribution pattern of the radial velocities, which are esse symmetrically distributed.For r/b = 0.92 and E = 0.04, 0.2, and 0.4, the skewness can be used to measure the asymmetry of the distribution) is 5.39, 5.70, and 5.84, tively, while the kurtosis (which can be used to measure the steepness of the distrib is 32.07, 35.63, and 37.34, respectively.Figures 11 and 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.

R PEER REVIEW
12 of 25  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 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.

Mass Flow Rate
Figure 13 shows the distribution of the mass flow rate in the circumferential direction at different 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 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 = 2.39 × 10 , the maximum mass flow rate is 27 times larger than the average value.Meanwhile, as the 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 varies.(This phenomenon is of course more pronounced when the 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 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

Mass Flow Rate
Figure 13 shows the distribution of the mass flow rate in the circumferential direction at different Re ϕ 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 Re ϕ 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 Re ϕ = 2.39 × 10 6 , the maximum mass flow rate is 27 times larger than the average value.Meanwhile, as the Re ϕ 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 Re ϕ varies.(This phenomenon is of course more pronounced when the Re ϕ 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.
therefore more likely to fracture; continuing to increase the   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 therefore more likely to fracture; continuing to increase the will also affect the dow stream blade.Figure 15 shows the distribution of the mass flow rate in the circumferential directi for different Euler numbers.As in the analysis above, the variation in the E does not affe the distribution pattern of the mass flow rate; the mass flow rate remains maximum in t In addition, as the Re ϕ 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 Re ϕ 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 direc for different θ values It is clear from the figure that the maximum values of the mass rate occur in the middle of the low-pressure zone when   = 5.12 × 10 5 , regardles the variation of the θ, and that the maximum values are almost equal.When the θ is s (θ = 0.1), backflow occurs at both the upstream and downstream borders, and the m flow rate at the upstream border is smaller than that at the downstream border; as t increases (θ = 0.2), the backflow zone disappears and the mass flow rate at the upstr 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 Re ϕ = 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, Re ϕ = 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.From the above analysis, it is clear that the   will increase as the   and the E increase.Figure 17 shows the distribution of the mass flow rate in the circumferential direction for different θ when   = 10,137,   = 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   should be slightly greater than 0.2.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 Re ϕ 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, Re ϕ = 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.

𝑤 𝜑
From the above analysis, it is clear that the   will increase as the   and the E increase.Figure 17 shows the distribution of the mass flow rate in the circumferential direction for different θ when   = 10,137,   = 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   should be slightly greater than 0.2.

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.
Based on Equation (1), the pressure along the radius of the disk can be calculated with Equation (2) by Hu et al. [24].
Dimensionless pressure is defined as follows: The pressure coefficient is defined as: 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 Re ϕ 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).

ER REVIEW
18 of 2  Figure 19 shows the distribution of pressure coefficients at different value where a, b, c and d represent different circumferential positions.As the increase the pressure coefficients at each location do not vary much, except for a slight increas Figure 19 shows the distribution of pressure coefficients at different Re ϕ values, where a, b, c and d represent different circumferential positions.As the Re ϕ 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 Re ϕ increases.However, by continuing to increase the Re ϕ , 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   , 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.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 Re ϕ , 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.

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 Re ϕ , 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]:

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 , 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]: By definition, and are related as follows:  It is clear that the thrust coefficients of the rotor and the stator are almost equal.= 0 and = 0 mean the blades are intact, while = 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 = 2.39 × 10 , 5.12 × 10 , 1.02 × 10 , and 2.39 × 10 , 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 increase.When the blades are already fractured, the thrust coefficient of the disks decreases significantly as the increases but as the 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 is small, increasing the 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 , the pressure coefficient at high radii increases (see the enlargement of Figure 19a).Therefore, continuing to increase the will increase the thrust coefficient slightly, according to Equation (6).This also means that the thrust coefficient is more sensitive to the whether blades fracture or not, especially when the 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.By definition, C F and C p are related as follows:

Moment Coefficient
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 Re ϕ = 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 Re ϕ increase.When the blades are already fractured, the thrust coefficient of the disks decreases significantly as the Re ϕ increases but as the Re ϕ 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 Re ϕ is small, increasing the Re ϕ 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 Re ϕ , the pressure coefficient at high radii increases (see the enlargement of Figure 19a).Therefore, continuing to increase the Re ϕ will increase the thrust coefficient slightly, according to Equation (6).This also means that the thrust coefficient is more sensitive to the Re ϕ whether blades fracture or not, especially when the Re ϕ 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.

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 Re ϕ , E and θ respectively.The moment coefficient is defined by Equation ( 7) [25]: coefficient of the rotor decreases while that of the stator increases, whether or not blades fracture.This is mainly because as the 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 Fig- ure 24, where a, b and c show the impact of the , 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 increases, the moment coefficient decreases.This is because the pressure coefficient gradually converges to zero along the radius as the increases.From Figure 24b,c, it can be seen that increasing the E and the θ increases the moment coefficient.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 Re ϕ 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 Re ϕ 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 Re ϕ , 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 Re ϕ increases, the moment coefficient decreases.This is because the pressure coefficient gradually converges to zero along the radius as the Re ϕ increases.From Figure 24b,c, it can be seen that increasing the E and the θ increases the moment coefficient.

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.

1.
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.

2.
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, Re ϕ = 5.12 × 10 5 , and 0.2 ≤ E ≤ 0.4).

3.
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.4.
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.

Figure 1 .
Figure 1.Computational model used in this paper.

Figure 14 .
Figure 14.Pressure coefficient of case D at plane z/s = 0.5.

Figure 14 .
Figure 14.Pressure coefficient of case D at plane z/s = 0.5.

Figure 14 .
Figure 14.Pressure coefficient of case D at plane z/s = 0.5.

Figure 18 .
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 .
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 20
Figure20shows 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 Figure21, where a, b, c and d represent different circumferential positions.Unlike the E and the Re ϕ , 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.

Table 1 .
The range of parameters.