Abstract

The damage due to particulate matter ingestion by propulsion gas turbine engines can be significant, impacting the operability and performance of plant components. Here, we focus on the axial compressor whose blades become damaged when operated in dusty/sandy environments, resulting in significant performance degradation. In this work, CFD studies are performed to model the effects of airfoil damage on the first-stage rotor blading of a GE T700–401C compressor. We use thermoplastic additive manufacturing to construct representative physical models of three damage morphologies—ballistically bent/curved leading edges, cragged erosion of leading edges, and eroded leading/tailing edges at outer span locations. The resultant damaged plastic geometries, and a baseline undamaged configuration are then optically scanned and incorporated into sublayer resolved full stage, unsteady RANS analyses. Boundary conditions are imposed that conform to damaged compressor operation protocols, and this iterative process for accommodating corrected mass flow and off-design powering is presented. The results for the three damaged and one undamaged configuration are studied in terms of compressible wave field and secondary/tip flows, spanwise performance parameter distributions and efficiency. A method to estimate the effect of rotor damage on engine SFC is presented. The code, modeling, and meshing strategies pursued here are consistent with a validation study carried out for NASA Rotor 37 — these results are briefly included, and provide confidence in the predictions of the T700 geometry studied. The results provide quantitative comparisons of, and insight into, the physical mechanisms associated with damaged compressor performance degradation.

1 Introduction

In certain environments, solid particles such as sand, dust, ice, and volcanic ash can be ingested by a gas turbine engine in significant quantities. The higher inertia of some of these particles can cause them to deviate from the gas path streamlines and impact the fore-stages of axial fans and compressors [1]. The damage to the blades is caused by direct impingement of large particles as well as the impact and deposition of finer particles due to secondary flows through the blade passage [2]. In the case of aeroengines, particle ingestion into the core occurs during all phases of flight operation, however, it is most critical during takeoff and landing where particle/dust densities are typically higher. Damage modes associated with particle impact with compressor blades and endwalls include pitting, material loss, and curling of blade leading and trailing edges, increased effective tip clearances, and erosive blade surface roughness. The aerodynamic performance of the compressor blades is affected by these geometry changes. This, in turn, affects the efficiency and operating range of the compressor.

A number of research groups have studied these various damage modes in axial compressors. Balan and Tabakoff [3] studied the effect of sand erosion on axial compressor cascade aerodynamics and performance deterioration. They observed particularly severe erosion at the rotor leading edge and on the pressure surface. Tabakoff et al. [4] simulated and investigated the effects of single stage erosion on overall compressor performance. Sallee [5] reported significant in-service damage to various components, including the high pressure compressor for the JT9D engine. Also reported was an assessment of the effect of compressor operation time on rotor blade damage. Ghenaiet et al. [6] assessed the damage to an axial fan due to sand ingestion through a qualitative study of erosion patterns and quantitative measurements of degradation in aerodynamic characteristics. Batcho et al. [7] performed experiments on the Pratt and Whitney TF33 turbofan and J57 turbojet engines to assess particle ingestion damage on the performance of the axial compressors, and determined the increased susceptibility to surge at low power settings.

Several groups have studied compressor damage using computational fluid dynamics (CFD). Li and Sayma [8] assessed the effect of tip curl damage on the stall margin of a transonic axial compressor. They observed that the presence of one damaged blade prevented the stall cells from growing and improved stability. Fedechkin et al. [9] modeled different “nicks” and bending of the first stage rotor blade leading edge in a four-stage axial compressor. They determined limits for each of these damage modes for observable performance degradation at multiple operating conditions. Suzuki and Yamamoto [10] performed two-phase CFD analysis of a single stage axial compressor. The particle trajectories and impaction were used to model erosion depth explicitly used in a subsequent CFD analysis. They reported localized damage patterns and performance deterioration consistent with experimental results.

The research above focuses on ballistically impaired geometry. Other researchers have focused their studies on less severe damage modes including increased surface roughness and fouling/deposition. Suder et al. [11] performed laser anemometry experiments and some attendant CFD analysis on roughened and thickened compressor blading and observed significant performance degradation (>5% reduction in pressure ratio and adiabatic efficiency). Most of these losses were attributed to leading edge thickening/roughening and transition effects. Morini et al. [12] performed CFD simulations of NASA Stage 37 to assess the effect of compressor fouling by imposing different combinations of added thickness and surface roughness to the baseline geometry. They observed that the geometric modifications led to decreased mass flowrate, efficiency, and stagnation pressure ratio. Gbadebo et al. [13] performed experimental and CFD studies on an artificially roughened low-speed compressor stator. Total pressure traverses and surface flow visualization showed a significant roughness induced increase in the size of the hub endwall corner separation, this leading to attendant increased loss, blockage, and deviation. Syverud et al. [14,15] studied the impact of salt deposits on a jet engine compressor and again observed significant deterioration on stage performance. They also quantified the distribution of the salt fouling on the rotor and stator blades within the machine. Walton et al. [16] also quantified the significant leading edge surface roughness of ex-service compressor blades, their research focusing on the metrology of the roughness rather than its performance impact.

In this work, we study several ballistic compressor blade damage morphologies observed in the field. Specifically, we aim to study local aerodynamic and attendant performance implications of representative geometric models of three of these modes—ballistically bent/curved leading edges, spanwise cragged erosion of leading edges, and eroded leading/tailing edges at outer spanwise locations. We employ thermoplastic additive manufacturing (AM) to build a number of baseline undamaged stage 1 rotor blades, and then apply heat and tooling treatments to obtain the damaged geometries which are then optically scanned and incorporated into sublayer resolved Reynolds-averaged Navier–Stokes (RANS) analysis.

The paper is organized as follows: First, we summarize the technical approaches associated with geometry generation, CFD meshing, and CFD modeling. We present results for the well-studied NASA Rotor 37 configuration to ascribe credibility to the modeling we perform for the T700 stage. In axial compressors, damage to the first stage rotor significantly impacts the instantaneous and average flow field in the rotor and downstream passages. Accordingly, we perform a full stage unsteady analysis in order to avoid the significant washing out of local aerodynamic changes encountered in a steady analysis. We present an approach for appropriately accommodating operating point variation for a damaged compressor. Next, we present comparisons of the flow field and performance predictions between the undamaged and three damaged blades. A short summary of how results such as these can be used to predict overall engine performance is included.

2 Technical Approach

The geometry of the first stage compressor rotor was obtained from GE and NAVAIR for the T700–401C helicopter engine [17]. Figure 1 shows several views of the stage 1 rotor blade geometry.

Fig. 1
Two views of the T700 first stage rotor blade
Fig. 1
Two views of the T700 first stage rotor blade
Close modal

2.1 Damaged Blades.

Since the damaged blades from the T700 experimental tests are forthcoming, an AM based approach was used to represent damage for the CFD studies. The blades were fabricated with a Creality3D Ender-3 V2 3D printer using polylactic acid plastic which is easily amenable to local mechanical and thermal deformation without compromising the shape and aerodynamic design of the blade. The damaged plastic blades were then 3D optically scanned to produce point clouds. The point clouds were then processed in solidworks [18] to produce 3D CAD models to be used for simulation.

Three damage modes were chosen for this study, designated here as, Cragged, Eroded, and Curled. These damage modes have been observed by many workers including Walton et al. [16], Aust and Pons [19], and Zuniga and Osvaldo [20]. Cragged and Eroded damages arise from the surface abrasion/wear of the blades over long exposure times. Curled tip damage arises due to impact from singular foreign large scale object ingestion. These and other damage modes are explored in Aust and Pons [19] and the definitions, descriptions, and photography in that reference were used to adapt geometries studied herein.

Figure 2 shows a compressor blade with representative cragged leading edge damage, and the corresponding 3D-printed blade used in this study. The blade leading edge has notches of depth 1 mm each at 20%, 35%, 46%, 65%, and 84% spans. The notches were made and then mechanically smoothed.

Fig. 2
Original [19] and 3D-printed blades with Cragged damage
Fig. 2
Original [19] and 3D-printed blades with Cragged damage
Close modal

Figure 3 shows a blade with leading and trailing edge chord material loss, and the corresponding 3D-printed blade. At the leading edge, there is blade material loss from 75% span to 25% chord length at 100% span. At the trailing edge, there is blade material loss from 70% span to 30% chord length at 100% span.

Fig. 3
Original [19] and 3D-printed blades with the Eroded damage
Fig. 3
Original [19] and 3D-printed blades with the Eroded damage
Close modal

Figure 4 shows a blade with a curled tip, and the corresponding 3D-printed blade. The blade tip is bent 45 deg toward the suction side along a line connecting the leading edge at 75% span to 10% percent chord at the tip. Further details of the 3D-printed damaged blades are available in the first author’s forthcoming Ph.D. thesis [21].

Fig. 4
Original [19] and 3D-printed blades with the Curled damage
Fig. 4
Original [19] and 3D-printed blades with the Curled damage
Close modal

2.2 Computation Domain and Mesh.

A computational domain with four rotor blades and six stator blades is considered for each case, and is shown in Fig. 5. This ratio corresponds very closely to the integer blade-vane-count in the T700 first stage. We use 4:6 rather than 2:3 to enable analysis of damage scenarios where all blades do not exhibit the same damage. In this work, we use three undamaged and 1 damaged blade for each four-blade annulus section.

Fig. 5
A view of the full stage computational domain
Fig. 5
A view of the full stage computational domain
Close modal

The mesher in star-ccm+ (v2021.1) [22] was used to create a sublayer resolved unstructured polyhedral mesh with prism layers at the blade surfaces and endwalls. A structured grid built in Pointwise [23] was used to extend the outlet domain one axial chord downstream of the stator trailing edge plane. This is done to reduce overall mesh count and to facilitate accurate spanwise averages at the stage exit. The tip clearance gap was explicitly resolved using 31 prism layers. Views of the mesh used are shown in Figs. 6 and 7.

Fig. 6
A 3D view of the mesh used for the simulation
Fig. 6
A 3D view of the mesh used for the simulation
Close modal
Fig. 7
Top view of the mesh at 50% span
Fig. 7
Top view of the mesh at 50% span
Close modal

2.3 Simulation Parameters and Boundary Conditions.

The RANS simulations were performed in star-ccm+ (v2021.1) [22]. All cases were run time accurate, using the available coupled solver, with second-order space and time discretization. Perfect gas air was the working fluid. The Menter-kω SST turbulence model [24] was used. An exactly conservative sliding interface was defined between the rotor and stator meshes. The model was run in the absolute frame-of-reference, with rotor mesh motion. A physical time-step of 1.37 × 10−6 s was used, corresponding to 50 time-steps per rotor passing period. Each simulation was run for 4000 rotor passing periods to achieve near statistical stationarity for the four rotor blades (1 damaged) : six stator vane configurations studied.

Rotationally periodic boundaries were used. At the inlet, radial profiles of the stagnation pressure and stagnation temperature were specified based on information provided by GE. Inlet turbulence intensity and length scale were specified to accommodate expected test conditions. At the outlet, hub static pressure was specified, and a simplified radial equilibrium pressure distribution imposed. The blade and endwall surfaces were defined as adiabatic non-slip walls.

A grid independence study was performed on a single undamaged rotor blade passage to determine an appropriate mesh for the CFD analyses. A 3.50 × 106 cell mesh and a finer 7.03 × 106 cell mesh were run. Predicted stagnation pressure ratio, stagnation temperature ratio, and adiabatic efficiency varied by less than 0.1% between the two meshes. Accordingly, we have used the meshing parameters of the 3.50 × 106 cell mesh for all damaged and undamaged rotors blade and stator passages in this paper. This led to full stage cell counts of approximately 3.4 × 107 cells for the 4:6 rotor:stator stage simulations presented below.

2.4 Target Operating Conditions.

Damaged blading implies off-design compressor operation. In order to accommodate real-world off-design operational conditions (i.e., pilot control), several considerations were made regarding the target conditions for the simulations with the damaged rotor blades. The propulsion plant output power needs to be maintained for adequate aircraft performance. So, in the T700 turboshaft application, the mass flowrate to, and the shaft speed of, the power turbine is maintained. Hence, for the stage with damaged blades, the target mass flowrate needs to match that with the undamaged rotor operating on-design. Therefore, compressor shaft speed must be increased. In addition, the stagnation enthalpy rises across the entire six-stage compressor (Δh0,1−6) should also match that of the undamaged compressor.

We note that the percentage of overall compressor stagnation enthalpy rise supplied by the damaged 1st stage rotor is lower than the undamaged rotor. Accordingly, since the mass flowrate and stagnation enthalpy change across the entire compressor is assumed unchanged, the increment in rotation speed means that the subsequent undamaged stages of the compressor perform more work. Hence, downstream of the damaged first stage, application of the Euler turbomachinery equation yields for each rotor
(1)
where subscripts i, and f denote the initial/on-design and final/off-design conditions respectively.

Additionally, as per meanline velocity triangles for these rotors, increasing the compressor shaft speed Ω leads to a change in the tangential component of the absolute velocity (vθ) at constant mass flowrate (vz remains the same) as shown in Fig. 8.

Fig. 8
Velocity triangles at design point (left), off-design point (at increased Ω) (right)
Fig. 8
Velocity triangles at design point (left), off-design point (at increased Ω) (right)
Close modal
A linear relationship between flow turning ratio (Δvθ,fvθ,i) and compressor shaft speed ratio (Ωfi) for the rotors is determined through a study of a single undamaged stage 1 rotor blade within the operational range of compressor shaft speed Ω. We illustrate this in Fig. 9 which plots the flow turning ratio versus compressor shaft speed ratio for five shaft speeds in the range of interest. Linear regression yields
(2)
Considering the undamaged compressor, the stagnation enthalpy change across the five axial and one centrifugal stages is given by Δh0,1−6,und and the stagnation enthalpy change for the damaged compressor can be expressed as
(3)
Hence, the target stagnation enthalpy change for the damaged first-stage rotor can be computed using Eqs. (1), (2), and (3) as
(4)
where in Eq. (4) the design intent stagnation enthalpy rise for each stage is known/supplied from the engine manufacturer (i.e., Δh0,1−6,und, Δh0,2−6,und) and subscripts i and f denote the initial/undamaged and final/damaged conditions respectively.
Fig. 9
Δvθ,fΔvθ,i versus ΩfΩi for five selected compressor shaft speeds in range of interest and linear fit
Fig. 9
Δvθ,fΔvθ,i versus ΩfΩi for five selected compressor shaft speeds in range of interest and linear fit
Close modal

In summary, for the damaged stage simulations, the rotation speed and the outlet hub static pressure are adjusted, during iteration, such that the converged mass flowrate matches that of the on-design mass flowrate, and the stagnation enthalpy rise for the damaged stage satisfies Eq. (4).

3 Results

3.1 Method Validation.

Experimental data from T700 field tests are forthcoming, therefore the CFD approach was validated against the well-studied NASA Rotor 37 transonic axial compressor. Rotor 37 was designed and tested by Reid and Moore [25], and retested by Suder [26]. The same meshing and simulation approaches presented above were used to analyze the rotor (here, single blade row, steady-state, relative frame-of-reference simulation) at a given experimental operating point near peak efficiency, i.e., at 98% of the experimental choked mass flowrate. Here, results are compared to the experimental measurements of Suder [26] and the CFD results of Ameri [27] and Bruna and Turner [28]. Table 1 shows a comparison of the overall performance parameters. Stagnation pressure ratio (SPR), stagnation temperature ratio (STR), and adiabatic efficiency predictions are in very close agreement with the experimental data and other CFD results. Figure 10 shows predicted contours of instantaneous relative Mach number at 70% span, illustrating the fairly complicated shock system, attendant boundary layer interaction, and wake field in this flow.

Fig. 10
Rotor 37 instantaneous relative Mach number contours at 70% span
Fig. 10
Rotor 37 instantaneous relative Mach number contours at 70% span
Close modal
Table 1

Rotor 37 overall performance parameters

SuderAmeriBruna and TurnerPresent
m˙/m˙c0.9820.9800.9850.980
SPR2.0912.0582.0872.063
STR1.2671.2671.2701.270
ηtt0.8780.8600.8660.852
SuderAmeriBruna and TurnerPresent
m˙/m˙c0.9820.9800.9850.980
SPR2.0912.0582.0872.063
STR1.2671.2671.2701.270
ηtt0.8780.8600.8660.852

Figures 11 and 12 show the spanwise distribution of mass weighted STR and enthalpy weighted SPR at the outlet. Again, good agreement is observed between the present results and the experiment and previous CFD analyses. Discrepancies between the CFD models and experiments are likely associated with the outer wall thermal boundary condition [28], hub leakage [29], and turbulence modeling [30]. Further discussion of the Rotor 37 analysis is beyond the scope of this paper. These results justify that the meshing and modeling methods used are suitable for the parametric study of the T700 stage carried out below.

Fig. 11
Rotor 37 spanwise distribution of stagnation temperature ratio at the outlet
Fig. 11
Rotor 37 spanwise distribution of stagnation temperature ratio at the outlet
Close modal
Fig. 12
Rotor 37 spanwise distribution of stagnation pressure ratio at the outlet
Fig. 12
Rotor 37 spanwise distribution of stagnation pressure ratio at the outlet
Close modal

3.2 Undamaged Stage.

Results for the baseline undamaged T700 first stage are presented first. Like Rotor 37, the first stage rotor is transonic and a complicated compressible wave field arises, as illustrated by the instantaneous relative Mach number contours at 90% span shown in Fig. 13. This and all subsequent instantaneous visualizations are obtained at the same physical time-step (i.e., relative rotor-stator positions). The bow shock propagates upstream on the suction side and intersects the adjacent suction surface aft of midchord as a lambda shock which decelerates the boundary layer there. Expansion waves/acceleration are observed on the pressure side. Downstream of the shock-boundary-layer interaction a weak normal shock appears and at this near casing location separates the pressure side boundary layer. The shock system in the blade passage clearly gives rise to significant boundary layer and wake thickening.

Fig. 13
Instantaneous relative Mach number contours at 90% span for the Undamaged case
Fig. 13
Instantaneous relative Mach number contours at 90% span for the Undamaged case
Close modal

Figure 14 shows a front view of instantaneous relative Mach number contours at (x/c)hub = 0.27 downstream of the hub leading edge. The spanwise variation of the shock location is evident. A low momentum region at the casing is induced by the impinging shock there. The undamaged rotor has a tip clearance of 0.4% of span, and the tip clearance flow is visualized in Fig. 15, using relative frame streamlines integrated through the time-averaged flow field. Double leakage above the adjacent blade is observed as the tip clearance streamlines migrate across the blade passage. Figure 16 shows a top view of the same streamlines, with a 0.98 span surface contoured by time-averaged relative velocity. In both views, we observe the interaction of the tip vortex with the passage shock, which is seen to redirect and decelerate the vortex. The redirected vortex structure develops as it migrates across the blade passage and impinges at the pressure side of the adjacent blade aft of midchord. These observations are consistent with the Rotor 37 studies of Yamada et al. [31].

Fig. 14
Instantaneous relative Mach number contours at (x/c)hub = 0.27: Undamaged
Fig. 14
Instantaneous relative Mach number contours at (x/c)hub = 0.27: Undamaged
Close modal
Fig. 15
Tip clearance gap streamlines of Undamaged case with contours of relative velocity at several axial planes in the passage
Fig. 15
Tip clearance gap streamlines of Undamaged case with contours of relative velocity at several axial planes in the passage
Close modal
Fig. 16
Streamlines in the tip clearance gap region with contours of relative velocity at 98% span
Fig. 16
Streamlines in the tip clearance gap region with contours of relative velocity at 98% span
Close modal

3.3 Damaged Configurations.

We present qualitative and quantitative comparisons of the damaged configurations with the baseline Undamaged case. The Cragged configuration shown in Fig. 2, is considered first. Figure 17 shows instantaneous relative Mach number contours at 90% span. In this figure, and all subsequent blade-to-blade representations, the damaged blade is labeled. Comparing the two contour plots, modest differences are observed, and this is somewhat expected since damage to the rotor blade is confined to the leading edge along the span. Slight differences in shock and boundary layer fields are observed, particularly in the passages on either side of the damaged blade. Downstream of the bow shock on the suction side, we observe an increased acceleration of the flow compared to the undamaged case. Although relative Mach numbers are supersonic, the axial Mach number does not exceed one, so upstream influence of the Cragged passage flow arises, leading to off-design physics in all four of the blade passages. The wake downstream of the cragged blade is thinner and dissipates more quickly at this location.

Fig. 17
Instantaneous relative Mach number contours at 90% span for the Cragged damage case
Fig. 17
Instantaneous relative Mach number contours at 90% span for the Cragged damage case
Close modal

Figure 18 shows instantaneous relative Mach number contours at 90% span for the Eroded damage case. Here, we can clearly observe the impact of chord loss at the leading and trailing edges, and the resultant modifications to the flow field. The leading edge and its attendant bow shock are further downstream, and a richer shock-expansion field arises. We see a stronger near-normal shock in the suction side passage adjacent to the damaged blade. There is significant difference in the flow physics in all the blade passages. Also, the shock emanating from the blade’s pressure surface, is significantly stronger than in the Undamaged case, and impinges on the wake of the adjacent blade, further downstream than in the Undamaged case. The wake behind the Eroded blade is small and quickly dissipates at this outer span location.

Fig. 18
Instantaneous relative Mach number contours at 90% span for the Eroded damage case
Fig. 18
Instantaneous relative Mach number contours at 90% span for the Eroded damage case
Close modal

Figure 19 shows instantaneous relative Mach number contours at 90% span for the Curled damage case. The bow shock is much stronger since the leading edge curl acts as a much larger leading edge. A reverse flow region arises immediately downstream of the curl feature. The flow field within the damaged blade pressure side passage is more complex, with multiple shock expansion interaction features arising. On the suction side, much more expansion is observed. The wakes of all four blade passages are significantly impacted.

Fig. 19
Instantaneous relative Mach number contours at 90% span for the Curled damage case
Fig. 19
Instantaneous relative Mach number contours at 90% span for the Curled damage case
Close modal

In Figs. 20 and 21, instantaneous relative stagnation pressure contours are plotted at (x/c)hub = 0.66 downstream of the hub leading edge for the damaged and undamaged cases. In the undamaged blade passages, a high loss region is observed at the casing-suction surface corner associated with the shock-tip clearance vortex physics discussed above (Fig. 15). A number of additional observations apply. First, a low p0,rel region arises on the suction surface near mid-span for all passages/cases. This high loss region is due to boundary layer interaction with the adjacent blade’s bow shock. The size and shape of this loss region are significantly impacted by all three of the damaged blade geometries. The curled tip and eroded blades in particular give rise to significantly larger regions of high loss adjacent to the blade surfaces in passages on either side of the damage.

Fig. 20
Instantaneous relative stagnation pressure contours at (x/c)hub = 0.66; Undamaged (top), Eroded (bottom)
Fig. 20
Instantaneous relative stagnation pressure contours at (x/c)hub = 0.66; Undamaged (top), Eroded (bottom)
Close modal
Fig. 21
Instantaneous relative stagnation pressure contours at (x/c)hub = 0.66; Cragged (top), Curled (bottom)
Fig. 21
Instantaneous relative stagnation pressure contours at (x/c)hub = 0.66; Cragged (top), Curled (bottom)
Close modal

Figures 22 and 23 show the tip clearance flow for the Eroded and Curled cases visualized with the same viewing angle and seeding locations as Fig. 15. Comparing the three streamline plots, they are generally similar, however some differences can be observed. In the Eroded case, the tip vortex remains more closely confined to the blade, not extending all the way to the pressure side of the adjacent blade. Unlike the Undamaged case, double leakage over the adjacent blade is not observed. In conjunction with Fig. 20, we see evidence of the tip vortex developing along the shroud creating a broad region of low momentum fluid. In the Curled case, we observe a significantly larger tip vortex, with migration across the blade passage. In conjunction with Fig. 21, we see the attendant regions of low momentum fluid near the shroud, closer to the pressure side of the adjacent blade. In this case, we observe minimal double leakage over the adjacent blade. The streamlines for the Cragged case are not included here, since the near-tip geometry is similar to the undamaged blade, and the streamlines are thereby quite similar.

Fig. 22
Tip clearance gap streamlines of Eroded case with contours of relative velocity at several axial planes in the passage
Fig. 22
Tip clearance gap streamlines of Eroded case with contours of relative velocity at several axial planes in the passage
Close modal
Fig. 23
Tip clearance gap streamlines of Curled case with contours of relative velocity at several axial planes in the passage
Fig. 23
Tip clearance gap streamlines of Curled case with contours of relative velocity at several axial planes in the passage
Close modal

Stage performance parameters are summarized in Table 2. These are defined based on computational domain outlet and inlet boundary values. These parameters have been normalized with the results for the Undamaged case. Specifically, enthalpy averaged SPR, mass weighted STR, and adiabatic efficiency (based on mass weighted enthalpies) are tabulated. As discussed in Sec. 2.4, specifying appropriate operating conditions for the damaged runs involves specifying a higher shaft speed and lower rotor 1 stagnation enthalpy rise, Δh0,1, and these are listed as well. All three damaged rotors exhibit significant reductions in the stagnation pressure and temperatures ratios, and adiabatic efficiency. Consistent with the observations above, the cragged rotor exhibits modest performance deterioration, followed by more significant reductions for the curled blade, and yet more severe impact for the eroded damage.

Table 2

Comparison of overall performance parameters

CaseUnd.Crag.Curl.Erod.
RPM/RPMund1.0001.0011.0031.004
Δh01h01,und1.0000.9750.9320.912
SPR/SPRund1.0000.9810.9490.932
STR/STRund1.0000.9950.9870.983
ηtt/ηtt,und1.0000.9910.9730.966
CaseUnd.Crag.Curl.Erod.
RPM/RPMund1.0001.0011.0031.004
Δh01h01,und1.0000.9750.9320.912
SPR/SPRund1.0000.9810.9490.932
STR/STRund1.0000.9950.9870.983
ηtt/ηtt,und1.0000.9910.9730.966

Figure 24 shows the spanwise distribution of normalized adiabatic efficiency at an axial plane 0.34 stator axial chord downstream of the stator trailing edge plane. All profile data have been normalized using the mid-span values of the Undamaged case. For all cases, the normalized adiabatic efficiency decreases across the span, from approximately 1.05 near the hub to 0.9 at 80% span, then dropping precipitously to near 0.6 at the casing, this decrease arising due to increasing Mrel and shock strength with radius, tip clearance flow, and relative casing rotation.

Fig. 24
Spanwise distribution of normalized adiabatic efficiency at the stator outlet
Fig. 24
Spanwise distribution of normalized adiabatic efficiency at the stator outlet
Close modal

Figures 25 and 26 show the spanwise distributions of mass-weighted STR and enthalpy-weighted SPR, also at an axial plane 0.34 stator axial chord downstream of the stator trailing edge plane. Figure 27 shows the spanwise profile of absolute flow angle at the rotor-stator interface. Flow angles are also referenced to the mid-span Undamaged case. A number of observations are forthcoming. As shown in Fig. 25, all four cases exhibit decreased stage work and rotor turning with span across the inner 80% span, a general characteristic of the stage design. All three of the damaged blades give rise to significant reductions in STR and attendant rotor under-turning across most of the span. This redistributive unloading is observed for all stages despite the damage itself being restricted to outer spans for the curled and eroded blades. As shown in Fig. 26, the stagnation pressure rise decreases correspondingly with the increasing severity of the damage morphology.

Fig. 25
Spanwise distribution of normalized stagnation temperature ratio at the stator outlet
Fig. 25
Spanwise distribution of normalized stagnation temperature ratio at the stator outlet
Close modal
Fig. 26
Spanwise distribution of normalized stagnation pressure ratio at the stator outlet
Fig. 26
Spanwise distribution of normalized stagnation pressure ratio at the stator outlet
Close modal
Fig. 27
Spanwise distribution of normalized absolute flow angle at the rotor-stator interface
Fig. 27
Spanwise distribution of normalized absolute flow angle at the rotor-stator interface
Close modal

The damaged stage performance profiles show generally consistent spanwise trends with the undamaged stage, since the qualitative physics of the damaged passage flow field remain largely the same (boundary/operating conditions, wave fields, secondary flows), and only one blade of every four is damaged in this study. However, there is less work done and increased losses across the entire span for all three damage modes. Damage severity can again be rank ordered from Cragged, which shows only modest deterioration of the four parameters considered, across the span, followed by Curled and Eroded, which exhibit highly compromised performance deterioration.

3.4 Impact on Specific Fuel Consumption.

In Sec. 2.4, a method was presented for setting appropriate operating and boundary conditions for a damaged compressor CFD analysis, i.e., as a pilot would operate the engine under these circumstances. Maintaining full power with significantly reduced stage 1 compressor efficiency of course detrimentally impacts overall engine specific fuel consumption (SFC). Estimating this powerplant SFC reduction is of importance to aircraft operators in that it impacts mission duration, fuel costs, and off-platform refurbishment planning.

Accordingly, a standard gas turbine cycle deck was written to estimate the impact of the stage 1 rotor damage modes studied here on expected turboshaft SFC. The goal is to have an easily adaptable/upgradable tool to assess compressor and overall engine performance using clean and damaged blade measurements (T700/NAVAIR), and present CFD. Referring to Fig. 28, a side-view cutaway of the T700 compressor is shown. For this machine we “split” the compressor into stage 1, stages 2–5, and the centrifugal stage, per the numbering shown in the figure. Using this numbering, performance stacking relations for incremental stagnation pressures and temperatures within the compressor are given in Eqs. (5)(10), leading to overall compressor stagnation pressure and temperature ratios Eqs. (11)(12).
(5)
(6)
(7)
(8)
(9)
(10)
(11)
(12)
By invoking reasonable/known assumptions for the efficiencies of the inlet, turbine, burner, and exhaust engine sections, as well as assumptions for fuel type, bypass ratio, and freestream Mach number (both zero here), and turbine/exhaust specific heat ratio, one can derive an expression for SFC following [32] (for example). The inputs to the model are the predicted or measured stage 1 stagnation pressure ratio, CPRS1, and adiabatic efficiency, ηS1, and the measured, computed (or design intent) values for the other compressor stages, CPRS2−5, CPRS3, ηS2−5, ηS3.
Fig. 28
Side cut-away view of T700 compressor section [33]
Fig. 28
Side cut-away view of T700 compressor section [33]
Close modal

In Fig. 29, predicted values of SFC, normalized by the Undamaged model value, are plotted for the four configurations. As the nature of the damage becomes more severe, the normalized SFC increases as expected, to a maximum of just over a 1.4% increase for the Eroded damaged stage. This is attended by an overall reduction in CPR (also normalized) of 6.8% as plotted in Fig. 30.

Fig. 29
Normalized SFC versus normalized stage 1 adiabatic efficiency
Fig. 29
Normalized SFC versus normalized stage 1 adiabatic efficiency
Close modal
Fig. 30
Normalized SFC versus normalized compressor pressure ratio
Fig. 30
Normalized SFC versus normalized compressor pressure ratio
Close modal

4 Conclusion

CFD studies were performed to model the effects of ballistic airfoil damage on axial compressor blading. Physically representative geometric models of three damage modes were constructed and analyzed using a RANS approach which was validated here against Rotor 37. A scheme to implement boundary and operating conditions consistent with damaged compressor operation was presented. Full stage unsteady calculations for the three damaged and one undamaged configuration were studied in terms of compressible wave field and secondary/tip flows, spanwise performance parameter distributions, and efficiency. A simple scheme to estimate the effect of rotor damage on engine SFC was presented. The results enabled the rank ordering of the various damage modes in terms of expected impact on compressor and engine performance. The present study invoked three representative damage morphology geometries, and a particular damaged blade count (every fourth blade). The validated best-estimate CFD methods pursued here are likely to provide good engineering estimates for similar damage morphologies.

Acknowledgment

This material is based upon work supported by the Office of Naval Research under grant N00014-19-1-2232, with Dr. Steven Martens as Technical Monitor.

Conflict of Interest

There are no conflicts of interest.

Data Availability Statement

The datasets generated and supporting the findings of this article are obtainable from the corresponding author upon reasonable request.

Nomenclature

r =

radius

m˙ =

mass flowrate

h0 =

specific stagnation enthalpy

p0,k =

stagnation pressure at inlet of rotor stage k

vz =

absolute velocity in axial direction

vθ =

absolute velocity in tangential direction

T0,k =

stagnation temperature at stage k inlet

m˙c =

mass flow rate at choke

crag =

cragged rotor 1 stage

curl =

curled rotor 1 stage

erod =

eroded rotor 1 stage

und =

undamaged rotor 1 stage

Δh0,k =

stagnation enthalpy rise for stage k

Δh0,1,tgt =

target Δh0,k for stage 1

Δvθ,k =

flow turning for rotor stage k

CPR =

compressor pressure ratio

SFC =

specific fuel consumption

SPR =

stagnation pressure ratio

STR =

stagnation temperature ratio

Greek Symbols

α =

absolute flow angle

ηtt =

adiabatic efficiency

Ω =

compressor shaft speed

Superscripts and Subscripts

a =

station 2a at outlet of stage 1

b =

station 2b at outlet of stage 5

f =

final/off-design condition

i =

initial/on-design condition

mag =

magnitude

mid =

mid-span value

rel =

relative frame-of-reference

SC =

centrifugal compressor stage

S1 =

axial compressor stage 1

S2 − 5 =

axial compressor stages 2–5

und =

undamaged rotor 1 stage

2 =

station 2 at compressor inlet

3 =

station 3 at compressor outlet

References

1.
El-Sayed
,
A.
,
2017
,
Aircraft Propulsion and Gas Turbine Engines
, 2nd ed.,
CRC Press
,
Boca Raton, FL
, pp.
927
1057
.
2.
Neilson
,
J.
, and
Gilchrist
,
A.
,
1968
, “
Erosion by a Stream of Solid Particles
,”
Wear
,
11
(
2
), pp.
111
122
.
3.
Balan
,
C.
, and
Tabakoff
,
W.
,
1984
, “
Axial Flow Compressor Performance Deterioration
,”
20th Joint Propulsion Conference
,
Cincinnati, OH
,
June 11–13
, p.
1280
.
4.
Tabakoff
,
W.
,
Lakshminarasimha
,
A. N.
, and
Pasin
,
M.
,
1990
, “
Simulation of Compressor Performance Deterioration Due to Erosion
,”
ASME J. Turbomach.
,
112
(
1
), pp.
78
83
.
5.
Sallee
,
G.
,
1978
, “
Performance Deterioration Based on Existing (Historical) Data: JT9D Jet Engine Diagnostics Program
,” Technical Report NASA-CR-135448,
NASA
.
6.
Ghenaiet
,
A.
,
Tan
,
S. C.
, and
Elder
,
R. L.
,
2004
, “
Experimental Investigation of Axial Fan Erosion and Performance Degradation
,”
Proc. Inst. Mech. Eng., Part A: J. Power Energy
,
218
(
6
), pp.
437
450
.
7.
Batcho
,
P. F.
,
Moller
,
J. C.
,
Padova
,
C.
, and
Dunn
,
M. G.
,
1987
, “
Interpretation of Gas Turbine Response Due to Dust Ingestion
,”
ASME J. Eng. Gas. Turbines Power.
,
109
(
3
), pp.
344
352
.
8.
Li
,
Y.
, and
Sayma
,
A.
,
2012
, “
Effects of Blade Damage on the Performance of a Transonic Axial Compressor Rotor
,”
Turbo Expo: Power for Land, Sea, and Air
,
Copenhagen, Denmark
,
June 11–15
, Vol.
8
, pp.
2427
2437
.
9.
Fedechkin
,
K. S.
,
Kuzmin
,
M. V.
, and
Marchukov
,
E. U.
,
2020
, “
The Design Analysis of the Damaged Compressor Blade Impact on the Aerodynamic Parameters
,”
Turbo Expo: Power for Land, Sea, and Air
,
Virtual, Online
,
Sept. 21–25
, Vol.
2A
.
10.
Suzuki
,
M.
, and
Yamamoto
,
M.
,
2011
, “
Numerical Simulation of Sand Erosion Phenomena in a Single-Stage Axial Compressor
,”
J. Fluid Sci. Technol.
,
6
(
1
), pp.
98
113
.
11.
Suder
,
K. L.
,
Chima
,
R. V.
,
Strazisar
,
A. J.
, and
Roberts
,
W. B.
,
1995
, “
The Effect of Adding Roughness and Thickness to a Transonic Axial Compressor Rotor
,”
ASME J. Turbomach.
,
117
(
4
), pp.
491
505
.
12.
Morini
,
M.
,
Pinelli
,
M.
,
Spina
,
P. R.
, and
Venturini
,
M.
,
2009
, “
CFD Simulation of Fouling on Axial Compressor Stages
,”
Turbo Expo: Power for Land, Sea, and Air
,
Orlando, FL
,
June 8–12
, Vol.
5
, pp.
331
342
.
13.
Gbadebo
,
S. A.
,
Hynes
,
T. P.
, and
Cumpsty
,
N. A.
,
2004
, “
Influence of Surface Roughness on Three-Dimensional Separation in Axial Compressors
,”
ASME J. Turbomach.
,
126
(
4
), pp.
455
463
.
14.
Syverud
,
E.
, and
Bakken
,
L. E.
,
2006
, “
The Impact of Surface Roughness on Axial Compressor Performance Deterioration
,”
Turbo Expo: Power for Land, Sea, and Air
,
Barcelona, Spain
,
May 8–11
, Vol.
5
, pp.
491
501
.
15.
Syverud
,
E.
,
Brekke
,
O.
, and
Bakken
,
L. E.
,
2007
, “
Axial Compressor Deterioration Caused by Saltwater Ingestion
,”
ASME J. Turbomach.
,
129
(
1
), pp.
119
126
.
16.
Walton
,
K.
,
Blunt
,
L.
,
Fleming
,
L.
,
Goodhand
,
M.
, and
Lung
,
H.
,
2014
, “
Areal Parametric Characterisation of Ex-service Compressor Blade Leading Edges
,”
Wear
,
321
, pp.
79
86
.
17.
GE Aviation
. “
T700-401C/-701C Turboshaft Engines Data Sheet
.” www.ge.com/aviation
18.
Dassault Systèmes
. “
SolidWorks 2020 Essentials
.” https://www.solidworks.com/sw/support
19.
Aust
,
J.
, and
Pons
,
D.
,
2019
, “
Taxonomy of Gas Turbine Blade Defects
,”
Aerospace
,
6
(
5
), p.
58
.
20.
Zuniga
,
V.
, and
Osvaldo
,
M.
,
2007
, “
Analysis of Gas Turbine Compressor Fouling and Washing on Line
,” Ph.D. thesis,
Cranfield University
,
Cranfield, UK
. http://hdl.handle.net/1826/2448
21.
Chirayath
,
E.
,
2023
, “
CFD Modeling of Particle Ingestion Damage and Its Impact on Multistage Axial Compressor Performance
,” Ph.D. thesis,
Pennsylvania State University
,
State College, PA
.
22.
Siemens
, “
STAR-CCM+ v2021.1
,” User Manual.
23.
Pointwise, Inc.
, “
Pointwise V18.2 R2, User Manual
.”
24.
Menter
,
F. R.
,
1994
, “
Two-Equation Eddy-Viscosity Turbulence Models for Engineering Applications
,”
AIAA J.
,
32
(
8
), pp.
1598
1605
.
25.
Reid
,
L.
, and
Moore
,
R. D.
,
1978
, “
Performance of Single-Stage Axial-Flow Transonic Compressor With Rotor and Stator Aspect Ratios of 1.19 and 1.26, Respectively, and With Design Pressure Ration of 1.82
,” Technical Report NASA-TP-1338,
NASA
.
26.
Suder
,
K. L.
,
1996
, “
Experimental Investigation of the Flow Field in a Transonic, Axial Flow Compressor With Respect to the Development of Blockage and Loss
,” Technical Report NASA-TM-107310,
NASA
.
27.
Ameri
,
Ali
,
2009
, “
NASA Rotor 37 CFD Code Validation Using Glenn-HT Code
,” Technical Report NASA/CR-2010-216235,
NASA
.
28.
Bruna
,
D.
, and
Turner
,
M. G.
,
2013
, “
Isothermal Boundary Condition at Casing Applied to the Rotor 37 Transonic Axial Flow Compressor
,”
ASME J. Turbomach.
,
135
(
3
), p.
034501
.
29.
Shabbir
,
A.
,
Celestina
,
M. L.
,
Adamczyk
,
J. J.
, and
Strazisar
,
A. J.
,
1997
, “
The Effect of Hub Leakage Flow on Two High Speed Axial Flow Compressor Rotors
,”
Turbo Expo: Power for Land, Sea, and Air
,
Orlando, FL
,
June 2–5
, Vol.
1
, p.
V02AT32A071
.
30.
Denton
,
J. D.
,
1997
, “
Lessons From Rotor 37
,”
J. Thermal Sci.
,
6
(
1
), pp.
1
13
.
31.
Yamada
,
K.
,
Furukawa
,
M.
,
Nakano
,
T.
,
Inoue
,
M.
, and
Funazaki
,
K.
,
2004
, “
Unsteady Three-Dimensional Flow Phenomena Due to Breakdown of Tip Leakage Vortex in a Transonic Axial Compressor Rotor
,”
Turbo Expo: Power for Land, Sea, and Air
,
Vienna, Austria
,
June 14–17
, Vol.
5
, pp.
515
526
.
32.
Hill
,
P. G.
, and
Peterson
,
C. R.
,
1992
,
Mechanics and Thermodynamics of Propulsion
, 2nd ed.,
Addison-Wesley
,
Reading, MA
.
33.
Cutawaycreations
,
2018
, “
GE Black Hawk T700 Engine
,” Retrieved from https://youtu.be/-Ox1wMXBo-A