# A Novel Approach to Evaluate the Benefits of Casing Treatment in Axial Compressors

Authors: Guillaume Legras; Isabelle Trébinjac; Nicolas Gourdain; Xavier Ottavy; Lionel Castillon
Journal: International Journal of Rotating Machinery
 (2012)

---

## Abstract

Passive control devices based on casing treatments have already shown their capability to improve the flow stability in axial compressors. However, their optimization remains complex due to a partial understanding of the related physical mechanisms. In order to quantitatively assess the interaction between slots and the blade tip flow, the present paper develops a novel analysis methodology based on a control-volume approach located in the rotor tip region. This methodology may be used for analyzing the casing treatment based on both axi- and non-axisymmetric slots design. The second issue of the paper focuses on the application of the current approach to better understand the effects of axi- and non-axisymmetric grooves in three different axial compressors which differ by the flow regime (subsonic/transonic) and the smooth casing shape (cylindrical/concave). Numerical simulations are performed, and results of the current approach with and without casing treatments are compared.


---

## Body

## 1. Introduction

Today, it is mandatory for compressor designers to improve performance in terms of efficiency and operating range characterized by the stall margin at low mass flow rate. One of the main difficulties encountered in this process is that compressor stall is not always controlled through normal aerodynamic design. Thus stall prevention techniques must be used, and one promising technology known to bring substantial stability for tip critical compressor rotor is casing treatment (CT) [1]. This passive control device consists of slots within the rotor casing and presents various types of geometries: circumferential grooves [2–9], non-axisymmetric slot-type CT [10–14] and self-recirculating flow channels [15–17], and honeycomb [18]. Hathaway [19] provides an extensive overview of the research studies over the last 50 years that attempt to uncover the physics behind the improvement in stall margin.However, this understanding is still not complete since experimental measurements in the near casing region are elaborate. Moreover, postprocessing numerical results using 2D or 3D views hardly allow a quantification and can cause in some cases a misunderstanding of the flow mechanisms. Therefore, CT design requires analysis methods that need to be as easy as possible to handle for aeroengine designers. In this context, Shabbir and Adamczyk [3] proposed an approach based on a budget analysis of the steady axial momentum equation close to the rotor casing. This methodology provides further insight into the flow mechanisms relevant to compressor stability systems. In fact, their paper shows the approach is innovative to quantitatively ascertain the influence of a design on the near casing flow and provides guidance on groove design along the axial direction. However, the Shabbir and Adamczyk methodology is restrained to the knowledge of changes in the balance of the steady axial momentum equation while CT interaction with the main flow is strongly complex (tridimensional and unsteady), thus requiring information coming from each of the Navier-Stokes equations.The current paper aims at further understanding and quantitatively diagnosing complex flow mechanisms such as those induced by CT. First, the paper described an extended Shabbir and Adamczyk model (ESA model) that calculates the budget analysis of the Navier-Stokes set of (un)steady equations. Then, the model is used to investigate the flow mechanisms induced by CT in three different axial compressors. The idea is here to have a large test range in order to study the correlation between CT and the flow regime (subsonic and transonic), CT and the original smooth wall (SW) geometry (cylindrical or conical), and finally the main difference between axi- and non-axisymmetric CT. Results are based on numerical 3D calculations. In order to apprehend this novel approach, the control-volume analysis is hereafter limited to the axial momentum equation.


## 2. Flow Analysis Methodology

### 2.1. Equations

A control volume fixed in time and located in the near casing flow is retained. This will provide a quantitative understanding of the relevant fluid mechanisms associated with SW and CT configurations. Thus, the balance between the various terms which appear in the Navier-Stokes equations is analyzed based on its finite volume formulation. For simplicity reasons, the equations are considered in the relative rotor frame and in cylindrical coordinates. Using the divergence theorem, the Navier-Stokes equations in integral form can be written as follows:(1)V∫∂t∂Q∂t=-∮∂Ω[BdAr+CdAθ+DdAz]+VT,
where Ar, Aθ, and Az are the projection areas of the control volume, V the constant volume of the control domain, Q the conservative variables, B, C, and D vectors resulting from the development of the advective and diffusive fluxes, and T the forces per unit volume (usually named source terms). These vectors are recalled in (2). The stresses in (2) include both viscous and Reynolds stresses. In the same manner, the heat flux takes into account the heat flux and the enthalpy turbulent diffusion flux. This last term and the Reynolds stresses are both approximated by an eddy viscosity model:
(2)Q=[ρρWrρWθρWzρEr],B=[ρWr(ρWr2+Ps)-τrrρWθWr-τθrρWzWr-τzr(ρEr+Ps-τrr)Wr-(Wθτrθ+Wzτrz+qr)],C=[ρWθρWrWθ-τrθ(ρWθ2+Ps)-τθθρWzWθ-τzθ(ρEr+Ps-τθθ)Wθ-(Wrτθr+Wzτθz+qθ)],D=[ρWzρWrWz-τrzρWθWz-τθz(ρWzz2+Ps)-τzz(ρEr+Ps-τzz)Wz-(Wrτzr+Wθτzθ+qz)],T=[0ρWθ2-Ps-τθθr+2ρωWθ+ρω2rρWθWr-τrθr-2ρωWr0ρω2rWθ].


### 2.2. Numerical Resolution

Since the objective is to understand the balance of the various terms using numerical CFD results, the semidiscretised in space for uncoupled time/space integration formulation of (1) is considered. This procedure corresponds to the finite volume method. For an individual basic hexahedral cell, this formulation is written as follows:(3)∂Q∂t=-1V[∑l=16F⋅NΣl-VT]=-1VR(Q),
where l designates the lth face bounding the cell, F the numerical approximation of the exact flux (including the below tensors B, C, and D), NΣl the external (nonunit) normal to the face Σl=[Arl,Aθl,Azl]T, and R the numerical modelling residual of variable Q.In order to access the information of all individual terms, a computation postprocessing tool has been developed. Since the volume does not depend on time, the source termsT can be directly determined using (2). Concerning the time derivative terms, the choice is done to determine those terms at instant n by computing the opposite of the numerical modelling residual at the same instant:
(4)V∂Qn∂t=-R(Qn).
Therefore, only the convective fluxes have to be computed (here with a 2nd-order centered Jameson scheme [21].Since most CFD codes solve the Navier-Stokes equations in the cartesian reference frame, a specific treatment was done to recover terms of (2) in cylindrical coordinates. To do so, the fluxes constitutive of terms in (2) are build up through projections of the cartesian fluxes into the cylindrical reference frame. This strategy, instead of applying a cylindrical spatial scheme, permits to ensure the same numerical modelling residual than that in the cartesian equations. Equation (5) presents the transformation matrix P for rotation by an angle of θ over the longitudinal direction (i.e., z). Equation (6) recalls the expressions for projections from cartesian to cylindrical coordinates of the face areas A, the relative velocity vector W, the stress tensor τ, and the heat flux vector q:
(5)P=[cos⁡θ-sin⁡θ0sin⁡θcos⁡θ0001](6)Arθz=P-T⋅Axyz,Wrθz=P-T⋅Wxyz,τrθz=P-T⋅τxyz⋅P,qrθz=P-T⋅qxyz.For sake of clarity, the nomenclature used by Shabbir and Adamczyk [3] is chosen. The operator Δ()=∑l=16() is introduced and characterizes the balance of flux on an individual basic cell. For example, the term Δ(ρWzWrAr) appearing in the axial momentum equation corresponds to the transport of the axial momentum across the radial faces of a basic grid cell.The objective is to understand the balance of the various terms by using data on a mesh grid. Thus, it is obvious that the control volume retained is based on this mesh. In consequence, the current approach is extended to a control volume composed of multiple grid cells. The operator∑r∑θ∑z()=∑r,θ,z() is introduced and realizes the cumulative sum on each individual grid cell. For example, the unsteady axial momentum equation can be written as follows:(7)∑r,θ,zΔ(ρWzWrAr)+∑r,θ,zΔ(ρWzWθAθ)+∑r,θ,zΔ(ρWz2Az)+∑r,θ,zΔ(PsAz)+∑r,θ,z(Fλ2)+∑r,θ,z(Fλ4)-∑r,θ,zΔ(τrzAr)-∑r,θ,zΔ(τθzAθ)-∑r,θ,zΔ(τzzAz)=R(ρWz)=V∂ρWz∂t,
where Fλ2 and Fλ4 correspond, respectively, to the 2nd and 4th numerical scalar artificial dissipation fluxes. These have been added to the equation due to the use of the 2nd-order centered Jameson scheme [21] in the numerical simulations. Therefore, the novel method allows the quantification of artificial dissipation fluxes on the balance of the axial momentum equation. Moreover, it is possible to localize regions and flow structures where these fluxes are active. Finally, notice that each term of (7) is homogeneous to a force per unit volume ([kg·m−2·s−1]).


### 2.3. Application of the ESA Model

#### 2.3.1. On Steady Flows

In the case of steady flow problems, time derivative terms are null leading to a balance of the advective and viscous fluxes and forces per volume unit. In practice, the precision of this equilibrium depends on the value of the numerical residualR. For numerical data, it tends to reach zero at convergence according to the machines’ precision. Equation (3) can be simplified as follows:(8)R(Q)≈0⟺∑i=16F⋅NΣi-VTΩ≈0.


#### 2.3.2. On Unsteady Flows

Concerning unsteady flow problems, time derivative terms need to be determined since they can strongly impact the equilibrium of the equations. Those are taken equal to-R(Q), which is consistent only if the result analyzed is a time-consistent solution (usually a periodic solution).The current approach can be applied both on unsteady and time-averaged solution allowing access of specific information of the numerical modelling residual. In fact, assuming that any flow variableQ can be decomposed into its mean term Q¯ plus its deterministic component Q′ (stochastic component Q′′ is neglected), that is, Q=Q¯+Q′, and by replacing each variable in (3), the residual operator can be written as(9)R(Q)¯=R(Q¯)+R(Q′)¯=0,
where R(Q)¯ denotes the time-averaged balance of equation and tends to zero since the unsteady problem is periodic in time. R(Q¯) designates the residual applied to the time-averaged flow field. R(Q′)¯ indicates the time-averaged effect of the unsteadiness (also called lumped deterministic source terms by Ratzlaff et al. [23]).


## 2.1. Equations

A control volume fixed in time and located in the near casing flow is retained. This will provide a quantitative understanding of the relevant fluid mechanisms associated with SW and CT configurations. Thus, the balance between the various terms which appear in the Navier-Stokes equations is analyzed based on its finite volume formulation. For simplicity reasons, the equations are considered in the relative rotor frame and in cylindrical coordinates. Using the divergence theorem, the Navier-Stokes equations in integral form can be written as follows:(1)V∫∂t∂Q∂t=-∮∂Ω[BdAr+CdAθ+DdAz]+VT,
where Ar, Aθ, and Az are the projection areas of the control volume, V the constant volume of the control domain, Q the conservative variables, B, C, and D vectors resulting from the development of the advective and diffusive fluxes, and T the forces per unit volume (usually named source terms). These vectors are recalled in (2). The stresses in (2) include both viscous and Reynolds stresses. In the same manner, the heat flux takes into account the heat flux and the enthalpy turbulent diffusion flux. This last term and the Reynolds stresses are both approximated by an eddy viscosity model:
(2)Q=[ρρWrρWθρWzρEr],B=[ρWr(ρWr2+Ps)-τrrρWθWr-τθrρWzWr-τzr(ρEr+Ps-τrr)Wr-(Wθτrθ+Wzτrz+qr)],C=[ρWθρWrWθ-τrθ(ρWθ2+Ps)-τθθρWzWθ-τzθ(ρEr+Ps-τθθ)Wθ-(Wrτθr+Wzτθz+qθ)],D=[ρWzρWrWz-τrzρWθWz-τθz(ρWzz2+Ps)-τzz(ρEr+Ps-τzz)Wz-(Wrτzr+Wθτzθ+qz)],T=[0ρWθ2-Ps-τθθr+2ρωWθ+ρω2rρWθWr-τrθr-2ρωWr0ρω2rWθ].


## 2.2. Numerical Resolution

Since the objective is to understand the balance of the various terms using numerical CFD results, the semidiscretised in space for uncoupled time/space integration formulation of (1) is considered. This procedure corresponds to the finite volume method. For an individual basic hexahedral cell, this formulation is written as follows:(3)∂Q∂t=-1V[∑l=16F⋅NΣl-VT]=-1VR(Q),
where l designates the lth face bounding the cell, F the numerical approximation of the exact flux (including the below tensors B, C, and D), NΣl the external (nonunit) normal to the face Σl=[Arl,Aθl,Azl]T, and R the numerical modelling residual of variable Q.In order to access the information of all individual terms, a computation postprocessing tool has been developed. Since the volume does not depend on time, the source termsT can be directly determined using (2). Concerning the time derivative terms, the choice is done to determine those terms at instant n by computing the opposite of the numerical modelling residual at the same instant:
(4)V∂Qn∂t=-R(Qn).
Therefore, only the convective fluxes have to be computed (here with a 2nd-order centered Jameson scheme [21].Since most CFD codes solve the Navier-Stokes equations in the cartesian reference frame, a specific treatment was done to recover terms of (2) in cylindrical coordinates. To do so, the fluxes constitutive of terms in (2) are build up through projections of the cartesian fluxes into the cylindrical reference frame. This strategy, instead of applying a cylindrical spatial scheme, permits to ensure the same numerical modelling residual than that in the cartesian equations. Equation (5) presents the transformation matrix P for rotation by an angle of θ over the longitudinal direction (i.e., z). Equation (6) recalls the expressions for projections from cartesian to cylindrical coordinates of the face areas A, the relative velocity vector W, the stress tensor τ, and the heat flux vector q:
(5)P=[cos⁡θ-sin⁡θ0sin⁡θcos⁡θ0001](6)Arθz=P-T⋅Axyz,Wrθz=P-T⋅Wxyz,τrθz=P-T⋅τxyz⋅P,qrθz=P-T⋅qxyz.For sake of clarity, the nomenclature used by Shabbir and Adamczyk [3] is chosen. The operator Δ()=∑l=16() is introduced and characterizes the balance of flux on an individual basic cell. For example, the term Δ(ρWzWrAr) appearing in the axial momentum equation corresponds to the transport of the axial momentum across the radial faces of a basic grid cell.The objective is to understand the balance of the various terms by using data on a mesh grid. Thus, it is obvious that the control volume retained is based on this mesh. In consequence, the current approach is extended to a control volume composed of multiple grid cells. The operator∑r∑θ∑z()=∑r,θ,z() is introduced and realizes the cumulative sum on each individual grid cell. For example, the unsteady axial momentum equation can be written as follows:(7)∑r,θ,zΔ(ρWzWrAr)+∑r,θ,zΔ(ρWzWθAθ)+∑r,θ,zΔ(ρWz2Az)+∑r,θ,zΔ(PsAz)+∑r,θ,z(Fλ2)+∑r,θ,z(Fλ4)-∑r,θ,zΔ(τrzAr)-∑r,θ,zΔ(τθzAθ)-∑r,θ,zΔ(τzzAz)=R(ρWz)=V∂ρWz∂t,
where Fλ2 and Fλ4 correspond, respectively, to the 2nd and 4th numerical scalar artificial dissipation fluxes. These have been added to the equation due to the use of the 2nd-order centered Jameson scheme [21] in the numerical simulations. Therefore, the novel method allows the quantification of artificial dissipation fluxes on the balance of the axial momentum equation. Moreover, it is possible to localize regions and flow structures where these fluxes are active. Finally, notice that each term of (7) is homogeneous to a force per unit volume ([kg·m−2·s−1]).


## 2.3. Application of the ESA Model

### 2.3.1. On Steady Flows

In the case of steady flow problems, time derivative terms are null leading to a balance of the advective and viscous fluxes and forces per volume unit. In practice, the precision of this equilibrium depends on the value of the numerical residualR. For numerical data, it tends to reach zero at convergence according to the machines’ precision. Equation (3) can be simplified as follows:(8)R(Q)≈0⟺∑i=16F⋅NΣi-VTΩ≈0.


### 2.3.2. On Unsteady Flows

Concerning unsteady flow problems, time derivative terms need to be determined since they can strongly impact the equilibrium of the equations. Those are taken equal to-R(Q), which is consistent only if the result analyzed is a time-consistent solution (usually a periodic solution).The current approach can be applied both on unsteady and time-averaged solution allowing access of specific information of the numerical modelling residual. In fact, assuming that any flow variableQ can be decomposed into its mean term Q¯ plus its deterministic component Q′ (stochastic component Q′′ is neglected), that is, Q=Q¯+Q′, and by replacing each variable in (3), the residual operator can be written as(9)R(Q)¯=R(Q¯)+R(Q′)¯=0,
where R(Q)¯ denotes the time-averaged balance of equation and tends to zero since the unsteady problem is periodic in time. R(Q¯) designates the residual applied to the time-averaged flow field. R(Q′)¯ indicates the time-averaged effect of the unsteadiness (also called lumped deterministic source terms by Ratzlaff et al. [23]).


## 2.3.1. On Steady Flows

In the case of steady flow problems, time derivative terms are null leading to a balance of the advective and viscous fluxes and forces per volume unit. In practice, the precision of this equilibrium depends on the value of the numerical residualR. For numerical data, it tends to reach zero at convergence according to the machines’ precision. Equation (3) can be simplified as follows:(8)R(Q)≈0⟺∑i=16F⋅NΣi-VTΩ≈0.


## 2.3.2. On Unsteady Flows

Concerning unsteady flow problems, time derivative terms need to be determined since they can strongly impact the equilibrium of the equations. Those are taken equal to-R(Q), which is consistent only if the result analyzed is a time-consistent solution (usually a periodic solution).The current approach can be applied both on unsteady and time-averaged solution allowing access of specific information of the numerical modelling residual. In fact, assuming that any flow variableQ can be decomposed into its mean term Q¯ plus its deterministic component Q′ (stochastic component Q′′ is neglected), that is, Q=Q¯+Q′, and by replacing each variable in (3), the residual operator can be written as(9)R(Q)¯=R(Q¯)+R(Q′)¯=0,
where R(Q)¯ denotes the time-averaged balance of equation and tends to zero since the unsteady problem is periodic in time. R(Q¯) designates the residual applied to the time-averaged flow field. R(Q′)¯ indicates the time-averaged effect of the unsteadiness (also called lumped deterministic source terms by Ratzlaff et al. [23]).


## 3. Description of the Investigated Axial Compressors

In order to understand the CT effects on axial compressors features, three axial compressors are numerically investigated with and without CT: NASA Rotor 37 [24, 25], CREATE [26], and CBUAA [22, 27]. They differ by the geometry of the smooth casing (cylindrical versus noncylindrical), the flow regime (subsonic versus transonic), and the CT geometry incorporated (circumferential grooves or slot-type CT). These differences are resumed in Table 1. The following paragraphs further describe the design features of these three compressors.Table 1
Differences between the axial compressors investigated.
Investigated compressorsNASA Rotor 37CREATE S1-R2CBUAAGeometry of the casingCylindrical—XXNoncylindricalX——Geometry of CTCircumferentialXX—Slot type——XType of regimeSubsonic—X—TransonicX—X


### 3.1. NASA Rotor 37

The well-known transonic experimental compressor NASA Rotor 37 (Figure1) was designed and tested at the NASA Lewis Research Center [24, 25]. The main characteristics are recalled in Table 2. More details about the compressor design, extensive results, and comparisons of numerical and experimental data may be found in the following references: Suder and Celestina [28], Chima [29], Denton [30], and Dunham [20]. The NASA Rotor 37 has been chosen since it is known for the occurrence of the blade tip located rotating stall phenomena at operating close to surge [31]. Note that the original configuration experimentally tested does not include any CT (i.e., experimental data are only available for the SW configuration). The implementation of casing grooves is only investigated through numerical simulations.Table 2
Design features of the NASA Rotor 37 [20].
Number of blades36Tip radius at leading edge252 mmAspect ratio1.19Hub-tip radius ratio0.70Tip solidity1.288Tip clearance height0.356 mmRotation speed17188.7 rpmTip speed457 m/sTotal pressure ratio2.106Adiabatic efficiency0.877Design mass flow20.188 kg·s−1Choked mass flow20.93±0.14 kg·s−1Figure 1
Cross-section of the NASA Rotor 37 [28].The CT geometry is taken from Legras et al. [8] and consists of 6 circumferential slots (Figure 5(a)). Each groove is 3 mm wide (10.9% of rotor tip axial chord) with a gap of 1.5 mm between each cavity. The height-to-width ratio is 0.7 : 1. The first groove starts at around x/C=5% of the tip chord upstream of the rotor tip leading edge. Concerning the design, the objective was not to optimize the CT geometry but to reach a comprehensive description of the change in flow properties generated by the grooves. Therefore, it was found interesting to investigate the impact of a slot that surmounts the blade tip leading edge, even if its usefullness is not obvious.


### 3.2. Multistage Compressor CREATE

CREATE (Compresseur de Recherche pour l’Etude des Effets Aérodynamiques et Technologiques, Figure2) is a research 3.5-stage axial compressor designed by SNECMA and tested at the Ecole Centrale de Lyon in LMFA research laboratory. This subsonic compressor is dedicated to aerothermal and aerodynamic studies. Its geometry and rotation speed are representative of high-pressure compressor median-rear blocks of modern turbojet engine. Its design features are recalled in Tables 3 and 4. More details about the compressor design can be found in Touyeras and Villain [26]. Thanks to the large amount of measurements, this experimental compressor has been the subject of different research works focused on rotor/stator interactions [32, 33] and CFD methods calibration and improvement [34–36].Table 3
Design features of CREATE.
Cylindrical outer casing diameter0.52 mRotation speed11543 rpmTip rotor speed313 m·s−1Table 4
Number of blades of CREATE rows.
RowIGVR1S1R2S2R3S3Number of blades per row (for2π)3264968011280128Number of blades for2π/162465758Figure 2
Cross-section of the CREATE high-pressure compressor and measurement planes.It was numerically seen at nominal speed that the stability of the untreated casing compressor is limited by the tip clearance effects on stages 2 and 3 [36]. As a result, it was chosen to implement on both stages a CT in the form of five axisymmetric grooves of equal width and equally distributed in the axial direction. Details of the design study can be found in Perrot et al. [4]. Note that, for confidentiality reasons, slots dimensions are not mentioned, and the following results are normalized. However, for the purpose of the current study, the CT flow analysis is only focused on the grooves implemented in the R2 rotor row. Furthermore, notice that the numerical results are based on Legras et al. [37]. Performances are obtained on the overall configuration R1-S3, whereas the CTs implemented in R2 are investigated using modelisation of the S1-R2 configuration.


### 3.3. CBUAA Compressor Test Case

The Compressor of the Beijing University of Aeronautics & Astronautics (CBUAA) is a single-stage transonic axial compressor (Figure3). A brief description of the test rotor can be seen in Table 5. In the experiment, the stator is placed far downstream of test rotor; thus the isolated rotor environment can be established. Experimental and numerical studies of the test case can be found in Lin et al. [22] and Ning and Xu [27].Table 5
Design features of the CBUAA [22].
Blade number17Tip diameter355.8 mmDesign speed22000 rpmHub-to-tip ratio0.565Aspect ratio0.956Mass flow at design point13.4 kg·s−1Total pressure ratio at design speed1.6 (stage)Adiabatic efficiency at design speed0.88 (stage)Relative Mach number of tip rotor speed1.404Figure 3
Cross-section view of the test section [22].The slot-type CT is illustrated in Figure4 [22, 27]. The total slot number is 153 (9 times the blade number). The slots are skewed to have a 45 degrees angle so that their openings face the pressure side of the blades. The axial length of the slots is 79% of the blade tip axial chord (with 39% exposed to the blade tip). The ratio of the open area of the slots to the casing annulus area is 0.8. The radial depth of the slots is 9 mm. Moreover, the CT has a recess chamber and its radial depth is 9 mm.Figure 4
Sketch of the CBUAA slot-type casing treatment [22].Chimera meshes used for modelling circumferential CT.
(a)
NASA Rotor 37(b)
CREATE


## 3.1. NASA Rotor 37

The well-known transonic experimental compressor NASA Rotor 37 (Figure1) was designed and tested at the NASA Lewis Research Center [24, 25]. The main characteristics are recalled in Table 2. More details about the compressor design, extensive results, and comparisons of numerical and experimental data may be found in the following references: Suder and Celestina [28], Chima [29], Denton [30], and Dunham [20]. The NASA Rotor 37 has been chosen since it is known for the occurrence of the blade tip located rotating stall phenomena at operating close to surge [31]. Note that the original configuration experimentally tested does not include any CT (i.e., experimental data are only available for the SW configuration). The implementation of casing grooves is only investigated through numerical simulations.Table 2
Design features of the NASA Rotor 37 [20].
Number of blades36Tip radius at leading edge252 mmAspect ratio1.19Hub-tip radius ratio0.70Tip solidity1.288Tip clearance height0.356 mmRotation speed17188.7 rpmTip speed457 m/sTotal pressure ratio2.106Adiabatic efficiency0.877Design mass flow20.188 kg·s−1Choked mass flow20.93±0.14 kg·s−1Figure 1
Cross-section of the NASA Rotor 37 [28].The CT geometry is taken from Legras et al. [8] and consists of 6 circumferential slots (Figure 5(a)). Each groove is 3 mm wide (10.9% of rotor tip axial chord) with a gap of 1.5 mm between each cavity. The height-to-width ratio is 0.7 : 1. The first groove starts at around x/C=5% of the tip chord upstream of the rotor tip leading edge. Concerning the design, the objective was not to optimize the CT geometry but to reach a comprehensive description of the change in flow properties generated by the grooves. Therefore, it was found interesting to investigate the impact of a slot that surmounts the blade tip leading edge, even if its usefullness is not obvious.


## 3.2. Multistage Compressor CREATE

CREATE (Compresseur de Recherche pour l’Etude des Effets Aérodynamiques et Technologiques, Figure2) is a research 3.5-stage axial compressor designed by SNECMA and tested at the Ecole Centrale de Lyon in LMFA research laboratory. This subsonic compressor is dedicated to aerothermal and aerodynamic studies. Its geometry and rotation speed are representative of high-pressure compressor median-rear blocks of modern turbojet engine. Its design features are recalled in Tables 3 and 4. More details about the compressor design can be found in Touyeras and Villain [26]. Thanks to the large amount of measurements, this experimental compressor has been the subject of different research works focused on rotor/stator interactions [32, 33] and CFD methods calibration and improvement [34–36].Table 3
Design features of CREATE.
Cylindrical outer casing diameter0.52 mRotation speed11543 rpmTip rotor speed313 m·s−1Table 4
Number of blades of CREATE rows.
RowIGVR1S1R2S2R3S3Number of blades per row (for2π)3264968011280128Number of blades for2π/162465758Figure 2
Cross-section of the CREATE high-pressure compressor and measurement planes.It was numerically seen at nominal speed that the stability of the untreated casing compressor is limited by the tip clearance effects on stages 2 and 3 [36]. As a result, it was chosen to implement on both stages a CT in the form of five axisymmetric grooves of equal width and equally distributed in the axial direction. Details of the design study can be found in Perrot et al. [4]. Note that, for confidentiality reasons, slots dimensions are not mentioned, and the following results are normalized. However, for the purpose of the current study, the CT flow analysis is only focused on the grooves implemented in the R2 rotor row. Furthermore, notice that the numerical results are based on Legras et al. [37]. Performances are obtained on the overall configuration R1-S3, whereas the CTs implemented in R2 are investigated using modelisation of the S1-R2 configuration.


## 3.3. CBUAA Compressor Test Case

The Compressor of the Beijing University of Aeronautics & Astronautics (CBUAA) is a single-stage transonic axial compressor (Figure3). A brief description of the test rotor can be seen in Table 5. In the experiment, the stator is placed far downstream of test rotor; thus the isolated rotor environment can be established. Experimental and numerical studies of the test case can be found in Lin et al. [22] and Ning and Xu [27].Table 5
Design features of the CBUAA [22].
Blade number17Tip diameter355.8 mmDesign speed22000 rpmHub-to-tip ratio0.565Aspect ratio0.956Mass flow at design point13.4 kg·s−1Total pressure ratio at design speed1.6 (stage)Adiabatic efficiency at design speed0.88 (stage)Relative Mach number of tip rotor speed1.404Figure 3
Cross-section view of the test section [22].The slot-type CT is illustrated in Figure4 [22, 27]. The total slot number is 153 (9 times the blade number). The slots are skewed to have a 45 degrees angle so that their openings face the pressure side of the blades. The axial length of the slots is 79% of the blade tip axial chord (with 39% exposed to the blade tip). The ratio of the open area of the slots to the casing annulus area is 0.8. The radial depth of the slots is 9 mm. Moreover, the CT has a recess chamber and its radial depth is 9 mm.Figure 4
Sketch of the CBUAA slot-type casing treatment [22].Chimera meshes used for modelling circumferential CT.
(a)
NASA Rotor 37(b)
CREATE


## 4. Compressors Modelisation

### 4.1. Description of the Numerical Methods

All numerical simulations are carried out using theelsA software developed by ONERA and CERFACS [38]. The code solves the Favre-Reynolds-averaged Navier-Stokes equations on multiblock structured meshes using a cell-centered finite-volume approach. It also allows the use of the Chimera method dedicated to modelize complex geometries typically generated by technological effects [39]. Thus, the Chimera method is here used to model the CT.Computations were run with a 2nd-order centered Jameson scheme [21] for the estimation of convective fluxes. The time marching is performed by an efficient implicit time integration scheme based on the backward Euler scheme and a scalar lower-upper (lu) symmetric successive overrelaxation (SSOR) method. Concerning the numerical approach, it differs depending on the investigated compressor case (see Table 6).Table 6
Numerical approaches performed according to the compressor investigated.
Numerical approachesInvestigated compressorsNASA Rotor 37 axi. CTCREATE S1-R2 axi. CTCBUAA non axi. CTSteadyX——Unsteady (phase-lagged)—XXFor the NASA Rotor 37 configuration, circumferential grooves are axisymmetric. They can be considered as a continuing circumferential modification of the casing and then be computed similarly to SW configuration in a conventional steady-state approach.Concerning the CREATE S1-R2 configuration, the use of an unsteady approach with phase-lagged boundary conditions at the rows’ interfaces and on the periodic boundaries enables to consider only one single-blade passage for each row [40, 41].For CBUAA case with slot-type CT, simulations used an unsteady overset grids’ method coupled with a phase-lagged condition at the interface between rotor and CT slot Chimera blocks [42]. This procedure reduces the computational domain to one single-blade passage and one CT slot.The turbulence is modelized by the two-equation modelk-ω proposed by Wilcox [43] for both the NASA Rotor 37 and the CBUAA, and the two-equation model of k-ε low-Reynolds of Launder-Sharma for the CREATE case. The flow is assumed to be fully turbulent since the mean Reynolds number based on the blade chord is approximately 5.106.Numerical validation and other modelling details are given, respectively, for(i)
the NASA Rotor 37 in Legras et al. [8, 9],(ii)
the CREATE case in Legras et al. [37],(iii)
the CBUAA in Legras et al. [44], Castillon and Legras [42].


### 4.2. Meshing Strategies

For all compressors, the flow domain is discretized with a low-Reynolds multiblock structured approach. Each SW configuration counts approximately a total of 1.2–1.5 million grid points. The meshes are clustered towards the solid boundaries in order to reach the resolution requirement ofy+≈1 (the size of the first layer is approximatively 1 μm). They are characterized by 89 grid points in the spanwise direction. The tip leakage region is discretized using a “O-H’’ grid topology with 25 points in the radial direction.Concerning the CT configurations, they are based on the SW meshes at which “H’’ Chimera blocks modelling the casing grooves have been added (Figures5 and 6).Figure 6
Chimera mesh used for modelling the slot-type CT of the CBUAA case.


## 4.1. Description of the Numerical Methods

All numerical simulations are carried out using theelsA software developed by ONERA and CERFACS [38]. The code solves the Favre-Reynolds-averaged Navier-Stokes equations on multiblock structured meshes using a cell-centered finite-volume approach. It also allows the use of the Chimera method dedicated to modelize complex geometries typically generated by technological effects [39]. Thus, the Chimera method is here used to model the CT.Computations were run with a 2nd-order centered Jameson scheme [21] for the estimation of convective fluxes. The time marching is performed by an efficient implicit time integration scheme based on the backward Euler scheme and a scalar lower-upper (lu) symmetric successive overrelaxation (SSOR) method. Concerning the numerical approach, it differs depending on the investigated compressor case (see Table 6).Table 6
Numerical approaches performed according to the compressor investigated.
Numerical approachesInvestigated compressorsNASA Rotor 37 axi. CTCREATE S1-R2 axi. CTCBUAA non axi. CTSteadyX——Unsteady (phase-lagged)—XXFor the NASA Rotor 37 configuration, circumferential grooves are axisymmetric. They can be considered as a continuing circumferential modification of the casing and then be computed similarly to SW configuration in a conventional steady-state approach.Concerning the CREATE S1-R2 configuration, the use of an unsteady approach with phase-lagged boundary conditions at the rows’ interfaces and on the periodic boundaries enables to consider only one single-blade passage for each row [40, 41].For CBUAA case with slot-type CT, simulations used an unsteady overset grids’ method coupled with a phase-lagged condition at the interface between rotor and CT slot Chimera blocks [42]. This procedure reduces the computational domain to one single-blade passage and one CT slot.The turbulence is modelized by the two-equation modelk-ω proposed by Wilcox [43] for both the NASA Rotor 37 and the CBUAA, and the two-equation model of k-ε low-Reynolds of Launder-Sharma for the CREATE case. The flow is assumed to be fully turbulent since the mean Reynolds number based on the blade chord is approximately 5.106.Numerical validation and other modelling details are given, respectively, for(i)
the NASA Rotor 37 in Legras et al. [8, 9],(ii)
the CREATE case in Legras et al. [37],(iii)
the CBUAA in Legras et al. [44], Castillon and Legras [42].


## 4.2. Meshing Strategies

For all compressors, the flow domain is discretized with a low-Reynolds multiblock structured approach. Each SW configuration counts approximately a total of 1.2–1.5 million grid points. The meshes are clustered towards the solid boundaries in order to reach the resolution requirement ofy+≈1 (the size of the first layer is approximatively 1 μm). They are characterized by 89 grid points in the spanwise direction. The tip leakage region is discretized using a “O-H’’ grid topology with 25 points in the radial direction.Concerning the CT configurations, they are based on the SW meshes at which “H’’ Chimera blocks modelling the casing grooves have been added (Figures5 and 6).Figure 6
Chimera mesh used for modelling the slot-type CT of the CBUAA case.


## 5. Overall Performances

Figures7, 8, and 9 present the measured and calculated performances at design speed, respectively, for the NASA Rotor 37, CREATE R1-S3 configuration, and CBUAA.Predicted NASA Rotor 37 characteristics at design speed.
(a)
Pressure rise(b)
Isentropic efficiencyCREATE characteristics relative to nominal operating point (experimental versus numerical simulations).
(a)
Pressure rise(b)
Isentropic efficiencyCBUAA characteristics at design speed (experimental versus numerical simulations).
(a)
Pressure rise(b)
Isentropic efficiency


### 5.1. Results of the NASA Rotor 37

In Figure7, the mass flow rates are normalized by the experimental choked mass flow (20.93 kg·s−1). The pressure ratio-mass flow curves show a fairly good agreement between numerical and experimental results for the SW configuration. However, the efficiency curve resulting from the simulation is clearly shifted which leads to underestimate the efficiency whatever the mass flow.The numerical results of the CT configuration show only slight differences with the SW case in terms of performance (Figure7). As observed by Legras et al. [8, 9], those small differences are all the more surprising that important radial flow exchanges occur between the 2nd and the 4th slots and the main flow. It is interesting to note that nearly all numerical studies that investigate CT for the NASA Rotor 37 present a relative insensitivity of the slots on the performance (circumferential grooves [45–47]; semicircular slots [31]). The application of the current methodology applied on the near casing flow region can be useful to ascertain the flow mechanisms, to quantify the fluxes and finally to understand why the implementation of CT is nonefficient in the NASA Rotor 37.


### 5.2. Results of the CREATE R1-S3 Configuration

Overall performances presented in Figure8 are normalized by the SW values at nominal operating conditions. Experimental uncertainties are 0.46%, 0.17%, and 0.32%, respectively, for the mass flow rate, the pressure ratio, and the isentropic efficiency. Concerning the steady numerical results, the limit of stability is estimated considering the last converged calculation.Results in Figure8 show that the shapes of the pressure rise and isentropic efficiency curves are correctly represented by the simulations. The main information of Figure 8 is that the numerical simulation is able to predict the trend observed in the measurement that CT improves performances compared to the SW case. The relative shifts in performance between both cases at design operating condition speed are approximately Δηis,max=+1%, ΔRpimax=+2%. For the steady simulations of the whole compressor being validated, let us analyze the unsteady calculations of the S1-R2 configuration.Moreover, Legras et al. [37] present local flow comparisons between CFD results and LDA measurements done in the rotor tip clearance and inside the CT of CREATE compressor. They proved that the CFD model correctly predicts the main characteristics of the flow in this region and particularly the rotor tip vortex, the flow inside the CT, and the flow interactions between the CT and the throughflow (i.e., radial flow exchanges).


### 5.3. Results of the CBUAA

Rotor characteristics of both SW and CT configurations obtained at 98% speed by the experiment [22] and predicted by unsteady simulations are compared in Figure 9. The numerical curves shapes are in fairly good agreement with the measurements. However, numerical characteristics are shifted due to a larger mass flow rate at chocked conditions than the measured value (as predicted by Lin et al. [22] and Ning and Xu [27]). Simulations overestimate total pressure as well as efficiency magnitudes. The stable operating range of the SW configuration is correctly predicted. Concerning CT case, simulations clearly underestimate the operating range due to the onset of stall point at higher mass flow rate than the experiment. Finally, the calculated efficiency curves of the CT configuration show lower magnitude than the SW case one in the SW operating range, whereas the magnitude of the total pressure CT curve is similar. These differences have also been observed by other authors who numerically investigated the present test case [22, 27].


### 5.4. Discussion of the Effect of CT Design on Overall Performances

The comparison of the three compressors performances indicates that the implementation of axisymmetric CT induces a low stall margin improvement and no penalty in compressor efficiency. On the contrary, non-axisymmetric CT greatly improves the stall margin but reduces the compressor efficiency. This observation suggests that the more the CT is efficient to enhance stall margin, the more the compressor is subject to a decrease of its efficiency. This observation is in agreement with results observed in the literature by Hathaway [19].Finally, based on the comparisons between experimental and numerical results, it was considered that the numerical models are able to reveal the overall flow mechanisms occurring near the outer casing. This is essential for the objective of the present work aiming at understanding the local influence of CT slots.


## 5.1. Results of the NASA Rotor 37

In Figure7, the mass flow rates are normalized by the experimental choked mass flow (20.93 kg·s−1). The pressure ratio-mass flow curves show a fairly good agreement between numerical and experimental results for the SW configuration. However, the efficiency curve resulting from the simulation is clearly shifted which leads to underestimate the efficiency whatever the mass flow.The numerical results of the CT configuration show only slight differences with the SW case in terms of performance (Figure7). As observed by Legras et al. [8, 9], those small differences are all the more surprising that important radial flow exchanges occur between the 2nd and the 4th slots and the main flow. It is interesting to note that nearly all numerical studies that investigate CT for the NASA Rotor 37 present a relative insensitivity of the slots on the performance (circumferential grooves [45–47]; semicircular slots [31]). The application of the current methodology applied on the near casing flow region can be useful to ascertain the flow mechanisms, to quantify the fluxes and finally to understand why the implementation of CT is nonefficient in the NASA Rotor 37.


## 5.2. Results of the CREATE R1-S3 Configuration

Overall performances presented in Figure8 are normalized by the SW values at nominal operating conditions. Experimental uncertainties are 0.46%, 0.17%, and 0.32%, respectively, for the mass flow rate, the pressure ratio, and the isentropic efficiency. Concerning the steady numerical results, the limit of stability is estimated considering the last converged calculation.Results in Figure8 show that the shapes of the pressure rise and isentropic efficiency curves are correctly represented by the simulations. The main information of Figure 8 is that the numerical simulation is able to predict the trend observed in the measurement that CT improves performances compared to the SW case. The relative shifts in performance between both cases at design operating condition speed are approximately Δηis,max=+1%, ΔRpimax=+2%. For the steady simulations of the whole compressor being validated, let us analyze the unsteady calculations of the S1-R2 configuration.Moreover, Legras et al. [37] present local flow comparisons between CFD results and LDA measurements done in the rotor tip clearance and inside the CT of CREATE compressor. They proved that the CFD model correctly predicts the main characteristics of the flow in this region and particularly the rotor tip vortex, the flow inside the CT, and the flow interactions between the CT and the throughflow (i.e., radial flow exchanges).


## 5.3. Results of the CBUAA

Rotor characteristics of both SW and CT configurations obtained at 98% speed by the experiment [22] and predicted by unsteady simulations are compared in Figure 9. The numerical curves shapes are in fairly good agreement with the measurements. However, numerical characteristics are shifted due to a larger mass flow rate at chocked conditions than the measured value (as predicted by Lin et al. [22] and Ning and Xu [27]). Simulations overestimate total pressure as well as efficiency magnitudes. The stable operating range of the SW configuration is correctly predicted. Concerning CT case, simulations clearly underestimate the operating range due to the onset of stall point at higher mass flow rate than the experiment. Finally, the calculated efficiency curves of the CT configuration show lower magnitude than the SW case one in the SW operating range, whereas the magnitude of the total pressure CT curve is similar. These differences have also been observed by other authors who numerically investigated the present test case [22, 27].


## 5.4. Discussion of the Effect of CT Design on Overall Performances

The comparison of the three compressors performances indicates that the implementation of axisymmetric CT induces a low stall margin improvement and no penalty in compressor efficiency. On the contrary, non-axisymmetric CT greatly improves the stall margin but reduces the compressor efficiency. This observation suggests that the more the CT is efficient to enhance stall margin, the more the compressor is subject to a decrease of its efficiency. This observation is in agreement with results observed in the literature by Hathaway [19].Finally, based on the comparisons between experimental and numerical results, it was considered that the numerical models are able to reveal the overall flow mechanisms occurring near the outer casing. This is essential for the objective of the present work aiming at understanding the local influence of CT slots.


## 6. CT Flow Mechanisms Analysis according to the ESA Model

This section is dedicated to the analysis of the flow interaction between rotor and CT at the last SW stable operating point. To do so, the ESA model is used to get further insight into the understanding of the phenomenon. In practical terms, a control volume analysis is here performed on the blade tip flow and comparisons between SW and CT configurations are done.


### 6.1. Equation Analyzed

As it has been observed in previous studies, the CT interaction with the blade tip flow is mainly tridimensional and unsteady. In order to handle the results interpretation, the scope of the CT flow mechanisms’ study is restrained to the analysis of the balance of the axial momentum (7). Thus, results in terms of axial forces acting on the blade tip flow are expected. This equation is chosen since it describes the global flow across the rotor channel, specially its pressure rise.Here, we notice that the CFD computations presented before has been performed with steady as well as unsteady approaches. Thus, for the NASA Rotor 37, steady axial momentum balance have been computed. On the contrary, for the CREATE and CBUAA cases, instantaneous axial momentum balances have been obtained for instants equally spread over a respective known temporal period:(i)
for CREATE, 25 instants equally spread over one upstream stator blade S1 passing temporal period (TS1);(ii)
for CBUAA, 17 instants spread over one CT slot passing temporal period (TCT).For these two last compressors, the contribution of the time derivative term∂ρWz/∂t can be computed, which cannot be done for NASA Rotor 37. In order to permit comparisons between the three compressors, CREATE and CBUAA instantaneous axial momentum balances are time-averaged on their respective temporal periods. In other words, NASA steady axial momentum balance is here compared to CREATE and CBUAA time-averaged axial momentum balances.


### 6.2. Definition of the Control Volume

For all configurations, the control volume surrounds the rotor blade tip flow field as shown in Figure10. It is circumferentially delimited by the rotor blade pitch, axially extended just upstream and downstream of the rotor blade tip, and radially bounded by the 5th/25 to the 20th/25 grid layers modelling the rotor tip clearance. Therefore, due to the Chimera technique in the grooved configurations, the control volume is not attached to the shroud. In fact, the two last radial layers normal to the casing are used for interpolation between grooves and the blade passage blocks. These are not taken into account in the control volume since they can mislead the Navier-Stokes equations equilibrium.Views of the control volume located the NASA Rotor 37 tip clearance.
(a)
3D view(b)
Meridian view


### 6.3. Global Flow Mechanisms Acting on the Near Casing Flow

Results of the axial momentum balance for the SW and CT configurations are presented, respectively, for NASA Rotor 37 in Figure11, CREATE in Figure 12, and CBUAA in Figure 13. This kind of representation, as histograms, provides a macroscopic view of the different axial forces acting on the control volumes (i.e., the near casing flow). As a consequence, one has to take care of the results interpretation since histograms do not describe the local influence of CT. This will be further done in the next section titled “Longitudinal evolutions of the axial forces and CT contribution’’.Figure 11
NASA Rotor 37 with and without CT: balance of the axial momentum equation.Figure 12
CREATE-R2 with and without CT: balance of the axial momentum equation.Figure 13
CBUAA with and without CT: balance of the axial momentum equation.Before discussing the physical analysis, it is instructive to comment on the magnitude of the numerical terms of the balance. Results show that the numerical modelling residualsR(ρWz) tend to zero leading to a good precision of the balance of the equation. For the CREATE and CBUAA cases, it signifies that the time-averaged temporal derivative term (i.e., =-R(ρWz)¯) obtained from unsteady simulations converges to zero. This observation lends support that the unsteady calculations correctly reach a time-periodic state. Moreover, it can be observed that scalar artificial viscosity fluxes (Fλ2 and Fλ4) are insignificant in spite of the presence of a strong shock wave located in the blade tip region of the NASA Rotor 37 and the CBUAA.


#### 6.3.1. SW Results

General Overview
For all SW configurations, results highlight the presence of four main axial forces being applied on the control volume:(i)
the axial pressure force∑Δ(PsAz) that corresponds to the axial flow pressure rise across the rotor.(ii)
The axial force due to the transport of the axial momentum across the longitudinal faces∑Δ(ρWz2Az). This term is similar to a dynamic pressure force that acts on the same surfaces than the pressure force. Its value relates to the flow axial acceleration or deceleration. In fact, by simply differentiating the terms between the axial outlet and inlet of the control volume (i.e., [(ρWz2Az)outlet-(ρWz2Az)inlet]) and assuming that Wz is mainly positive in the volume (which is not locally true due to the presence of the tip leakage, casing boundary layer,…), a negative (resp. positive) force corresponds to a deceleration (resp. acceleration) of the flow. The same reasonment can be done with Wz<0 that leads to the inverse conclusion. Results in Figures 11, 12, and 13 indicate that the incoming flow decelerates along the axial direction. This observation is coherent with the general aerodynamic of compressor since the casing boundary layer and the tip leakage flow greatly influence the incoming flow. Furthermore, the term ∑Δ(ρWz2Az) gives information on the axial flow blockage occurring in the control volume.(iii)
The axial shear force acting on the radial faces of the control volumes-∑Δ(τrzAr) induced by casing and blade tip boundary layers as well as the tip leakage flow.(iv)
The axial force due to the transport of the axial momentum across the radial faces∑Δ(ρWzWrAr). This term expresses the influence of the radial flow on the balance of axial forces. This force acts on the same faces than the axial shear force -∑Δ(τrzAr). Thus, it is believed that a change of one of these forces can strongly influence the other.All other terms are zero due to periodicity (forces acting onAθ faces) or to the radial thickness of the control volume. It is worth noticing that the forces applied on the near casing flow correspond to those observed in similar works done by Shabbir and Adamczyk [3].These axial forces can be distinguished depending on whether they are acting in the same direction of the flow (z>0) or in opposite direction of the flow (z<0). This information is given by the sign of the terms. A sketch describing the nomenclature of the axial forces acting on the control volume is presented in Figure 22. The results show for all compressors some coherent trends:(i)
the axial force due to∑Δ(PsAz) always acts in the opposite direction of the flow advance since the term is positive;(ii)
the axial forces due to∑Δ(ρWz2Az) and -∑Δ(τrzAr) always act in the same direction as the flow advance since the terms are negative.This observation is consistent with the general aerodynamic of axial compressor. The casing boundary layer (represented by-∑Δ(τrzAr)) and the upstream flow advance (described by ∑Δ(ρWz2Az)) balance the adverse pressure gradient (∑Δ(PsAz)) created by the axial compressor. If the boundary layer separates from the casing and/or if the flow advance is too low (like when the mass flow rate is too low), the pressure rise cannot be held, thus leading to an unstable operating flow condition.The different compressors differ by the magnitude of these axial forces and the direction of the∑Δ(ρWzWrAr).Influence of the Casing Geometry
Figures11, 12, and 13 clearly highlight the impact of the annulus shape (cylindrical versus conical). In fact, both CREATE and CBUAA compressors, which have a cylindrical casing, show similar results characterized by(i)
same relative balance between∑Δ(PsAz) and the opposite axial forces ∑Δ(ρWz2Az) and -∑Δ(τrzAr),(ii)
insignificant influence of∑Δ(ρWzWrAr).On the contrary, the NASA Rotor 37 that has a noncylindrical casing shows the occurrence of another axial force that is due to ∑Δ(ρWzWrAr). This force comes to contribute with the adverse pressure force ∑Δ(PsAz) to counteract the flow advance (i.e., z<0). This result is clearly obvious since the conical casing forces the flow to have a radial velocity component Wr<0. The near casing flow deflection is large enough to affect the balance of axial momentum, and it manifests itself through the axial force ∑Δ(ρWzWrAr).Influence of the Flow Regime
The comparison between CREATE and CBUAA (Figures12 and 13), which, respectively, works in subsonic and transonic flow regime, indicates no major difference in the axial momentum balance. Furthermore, the comparison between NASA Rotor 37 and CBUAA (Figures 11 and 13) proves that the shape of the rotor casing has a greater influence than the flow regime.


#### 6.3.2. CT Results

General Overview
The results of the CT configurations in Figures11, 12, and 13 show the same type of forces acting on the blade tip flow than SW cases. However, the influence of CT is here revealed by changes of the axial forces contribution.For all compressors, the magnitude of pressure rise across the rotor∑Δ(PsAz) follows the predicted overall performances (see Figures 7, 8, and 9).Due to the radial flow exchanges between the grooves and the throughflow (see [9, 37, 44] for images of the flow exchanges), CT is expected to contribute to balance the axial momentum equation with additional terms. As previously assumed in the SW analysis, the radial flow exchanges might influence the axial momentum equation through the force ∑Δ(ρWzWrAr). In fact, results are in agreement with this hypothesis. They clearly show that the main impact of the CT, whether it is axisymmetric or not, is revealed by the changes of ∑Δ(ρWzWrAr) and the viscous force -∑Δ(τrzAr) due to the casing boundary layer. As previously announced, these two forces might influence each other since they both act on the same faces of the control volume. In fact, results reveal that the axial force ∑Δ(ρWzWrAr) changes proportionnaly to -∑Δ(τrzAr). These observations can be interpreted as follows: the radial fluid exchanges between the CT slots and the rotor channel tend to replace the casing boundary layer in its role to balance the adverse pressure gradient. This suggests that the Wr magnitude is a key point to the CT flow mechanisms. If the reasonment is pushed further, one explanation of the CT efficiency to improve the stall margin might come from its capability to perform an axial force due to the radial flow exchanges (through the axial force ∑Δ(ρWzWrAr)) more robust to the adverse pressure gradient than the casing boundary layer (through the axial force -∑Δ(τrzAr)). Moreover, since the shear axial force can induce a loss of efficiency, the use of radial flow exchanges as an alternative to casing boundary layer effect is very interesting in order to increase the global efficiency of the compressor rotor. As a consequence, the main difference between axisymmetric and non-axisymmetric CT can be evaluated depending on the Wr evolution.The results of the NASA Rotor 37 (Figure11) and CREATE (Figure 12) indicate that circumferential CT contribute with additional axial force ∑Δ(ρWzWrAr)<0 acting on the flow advance. As expected, it is the effect of the radial flow exchanges and this observation is in agreement with results of Shabbir and Adamczyk [3]. This adding term ∑Δ(ρWzWrAr)<0 is accompanied with a proportional decrease of the viscous axial force -∑Δ(τrzAr). It can be relied to the slight gain in efficiency observed in CREATE performances (Figure 8). As a matter of fact, the main role of axisymmetric CT is to turn or drive the force ∑Δ(ρWzWrAr) more along the main flow direction, that is, z>0. This is consistent with the general conclusions that previous studies have found that the reverse tip leakage flow through the tip gap got pushed into the grooves by the re-injecting flow (in main flow direction) from the grooves.On the contrary, the CBUUA slot-type CT contributes by additional force∑Δ(ρWzWrAr)>0 acting in opposite direction to the flow as well as ∑Δ(PsAz) (Figure 13). Moreover, the viscous axial force -∑Δ(τrzAr) has increased proportionally, thus explaining the slight decrease in efficiency observed in CBUAA performances (Figure 9). This first result reveals that slot-type CT behaves in a different manner than circumferential CT. However, as mentioned below histograms present a global and not a local influence of the CT. For this particular CT geometry, a local analysis has to be done to further understand the flow mechanisms as it will be seen in the next paragraphs.For both axi- and non-axisymmetric grooves, the balance change is made possible through the modification of the radial velocityWr induced by flow exchanges between grooves and the throughflow. It further supports that Wr is quite important for the flow stabilization in the tip region. This observation has to be confirmed by the local analysis of the CT flow mechanisms.


#### 6.3.3. Unsteady Flow Mechanisms

Figures14 and 15 present, respectively, for CREATE and CBUAA the temporal evolution of the unsteady axial momentum balance. This analysis allows us to gain a quantitative and qualitative knowledge of the contribution of the time derivative term ∂ρWz/∂t.Figure 14
CREATE-R2 with and without CT: temporal evolution of the main terms of the unsteady axial momentum equation.Figure 15
CBUAA with and without CT: temporal evolution of the main terms of the unsteady axial momentum equation.Results of the CREATE-R2
For the CREATE case, unsteady results are presented over the upstream stator S1 passing temporal period. Plots in Figure14 show a sinusoidal time-dependent evolution of the time-derivative term influenced by the stator S1 blade passage frequency TS1. Its time-averaged value tends to zero, in agreement with results shown in Figure 12. Concerning its temporal fluctuations, amplitudes are particularly pronounced for the CT case with higher values of a factor 2 compared to the SW. Moreover, results show that ∂ρWz/∂t amplitudes are mainly driven by those of ∑Δ(PsAz) in the SW case and those of ∑Δ(ρWzWrAr) in the CT case. In fact, CT deadens the magnitude of the ∑Δ(PsAz) fluctuations by a factor 2 which is believed to greatly enhance the flow stability (because the pressure gradient opposes the flow). The increase in ∑Δ(ρWzWrAr) fluctuation, which drives the time-dependent term, proves that groove mechanisms strongly respond to upstream flow unsteadiness.Results of the CBUAA
Concerning the CBUAA case, unsteady results are presented over the CT slot passing temporal periodTCT. Plots in Figure 15 indicate that all terms present a “quasi-steady’’ behaviour, supporting that unsteadiness of the CT interaction with the throughflow is insignificant. This observation is especially confirmed by the zero magnitude over TCT of the time-derivative term ∂ρWz/∂t. This “quasi-steady’’ behaviour has already been observed by Lin et al. [22], and Ning and Xu [27]. Indeed, the unsteady time scale of the discrete treatment slots is therefore considerably small compared to that of the blade passage [27]; thus the unsteadiness of blade passage is weak. One reason to this phenomena might come from the ratio between the slot area and the area between the individual slots (which is very small here). Despite the high temporal resolution of the CT slot passing temporal period (TCT = 170 iterations), this particular test case fails to meet expectations in terms of analysis of unsteady flow interaction between the CT slots and the rotor channel. However, the effort is now focused on the analysis of local changes of the axial momentum balance in order to get more qualitative and quantitative information on the contribution of the CT slots.


### 6.4. Longitudinal Evolutions of the Axial Forces and CT Contribution

Since axial forces due to-∑Δ(τrzAr), ∑Δ(ρWzWrAr), and ∑Δ(ρWz2Az) highlight the change between SW and CT cases, it is interesting to investigate their respective axial evolution to further understand the flow interaction between CT and the throughflow. Figures 16 and 18 present the cumulative sum along the axial direction of these terms, respectively, for the NASA Rotor 37 and CBUAA. Figure 17 recalled from Legras et al. [37] the cumulative sum for the CREATE compressor of both steady and unsteady axial momentum balances. In the current study, only the unsteady curves are analyzed. Further comments on the difference between the steady and unsteady results are discussed in Legras et al. [37]. These diagrams provide information on the region where flow is highly constrained and on individual grooves contribution. The shaded bands denote the axial location of the casing grooves.Figure 16
NASA Rotor 37 with and without CT: cumulative sum in the axial direction of the principal terms involved in the axial momentum equation.CREATE-R2 with and without CT: cumulative sum in the axial direction of the principal terms involved in the axial momentum equation [37].
(a)(b)(c)Figure 18
CBUAA with and without CT: cumulative sum in the axial direction of the principal terms involved in the axial momentum equation.


#### 6.4.1. SW Results

Results of all SW configurations show that the efforts evolve in amplitude in a nonlinear manner along the axial direction.The curve of viscous term-Δ(τrzAr) begins to grow near the leading edge. To relate these trends to the flow physics, Figure 19(a) presents at the mid-height of the CBUAA rotor tip gap, the flow field of -Δ(τrzAr). It can be observed that the growth of -Δ(τrzAr) can be relied firstly to the tip leakage vortex (referred as region A in Figure 19(a)) and secondly to the interaction of the detached shock with the shroud boundary layers (in a less pronounced manner).CBUAA: Flow field denoted at midgap by negative values of-Δ(τzrAr).
(a)
Smooth wall(b)
Grooved wallConcerning the curve of the termΔ(ρWz2Az), it begins to grow in negative value before the blade leading edge for the NASA Rotor 37 and CBUAA cases. Contour plots of Δ(ρWz2Az) in the CBUAA rotor tip gap presented. Figure 20(a) reveals that this growing trend is linked to the flow acceleration before the blade detached shock. Next, the curve slightly increases in positive values. This observation relates to the incoming flow deceleration due to its interaction with the tip leakage flow. Further downstream in the blade channel, the curve increases in negative values due to the presence of the tip leakage flow.CBUAA: Flow field denoted at midgap by negative values ofΔ(ρWz2Az).
(a)
Smooth wall(b)
Grooved wallThe curve of termΔ(ρWzWrAr) axially evolves in “quasi-’’ opposite magnitude of Δ(ρWz2Az) curve. In fact, Wr acts here as a “coefficient of amplification.’’ This observation is confirmed by Figure 21(a), where contour plots are very similar to those of Δ(ρWz2Az) in Figure 20(a).CBUAA: Flow field denoted at midgap by negative values ofΔ(ρWzWrAr).
(a)
Smooth wall(b)
Grooved wallFigure 22
Nomenclature and sketch of the axial forces acting on a control volume located in the rotor tip clearance.To resume, the SW results highlight strong axial force magnitude due to the tip leakage flow and, when a detached blade shock exists, due to the interaction between leakage flow and blade passage shock. However, the detached blade shock has a lower influence in term of force magnitude than the tip leakage flow. It explains why the axial momentum balances of CREATE (subsonic regime) and CBUAA (transonic regime) are very similar. Concerning the difference ofΔ(ρWzWrAr) and Δ(ρWz2Az) between NASA Rotor 37 and CBUAA, it can clearly be explained by the casing shape. In fact, the rapid growth of Δ(ρWzWrAr) and Δ(ρWz2Az) curves in NASA Rotor 37 corresponds to the location of the greatest slope of the casing form.


#### 6.4.2. CT Results

Results of the CT configurations clearly highlight the local contribution of the slots.For both axi- and non-axisymmetric CT, the curve shapes of the viscous term-Δ(τrzAr) are very similar to those of the SW cases, except that CT delays further downstream the beginning of the growth rate. Since this term can be relied to the tip leakage flow, this observation can be explained by an attenuation of the spread of the leakage flow perpendicular to the blade. This is confirmed when comparing Figures 19(a) and 19(b). Moreover, referring to the CBUAA histograms in Figure 13, CT configuration shows higher magnitude of viscous term -Δ(τrzAr). In fact, contour plots in Figure 19(b) indicate the presence of a second region of high viscous axial force denoted B and located along the rear part of the CT slots. As seen by Legras et al. [44], this part of the slot is responsible for flow bleeding, whereas at the front part the recirculated flow is reinjected into the throughflow. This observation suggests that the flow bleeding mechanism induces strong shear layer that can increase viscous losses, thus explaining the efficiency losses observed with the use of the slot-type CT.The CT effects are even more highlighted by the curve ofΔ(ρWzWrAr). In fact, one can notice that axisymmetric CTs in NASA Rotor 37 and CREATE cases create an axial force Δ(ρWzWrAr)<0 in opposite direction to the flow advance. This explains the change observed in the previous histograms (Figures 11 and 12). This creation of axial force Δ(ρWzWrAr)<0 can be relied to the bleeding effect of the grooves since Wr>0 has been observed in previous studies [9, 37]. Moreover, this kind of analysis permits to identify the groove effectiveness depending on its capability to create this previous axial force, thus a bleeding effect. For the NASA Rotor 37, results show that the 2nd, 3rd, and 4th grooves are useful at the last stable point numerically seen for SW configuration. For the CREATE, Legras et al. [37] conclude that the 3rd groove is efficient at the last stable point numerically seen for SW configuration.Concerning the CBUAA slot-type CT, the bleeding effect along the rear part of the CT slots creates an axial force∑Δ(ρWzWrAr)<0 directed along the flow advance. This observation is confirmed by Figure 21(a) and is refered as zone D. As a matter of fact, the bleeding mechanisms of slot-type CT have a relatively similar effect than observed with circumferential CT. This is consistent with the general conclusions that previous studies have found that the reverse tip leakage flow through the tip gap got pushed into the grooves by the reinjecting flow (in main flow direction) from the grooves. On the contrary, the flow blowing mechanism along the front part of the CT induces axial force of ∑Δ(ρWzWrAr)>0 acting in opposite direction to the flow advance (refered as zone C in Figure 21(a)). In fact, the reinjected flow tends to radially deviate the incoming flow.Concerning the term ofΔ(ρWz2Az), curves mainly evolve in opposite magnitude to the term ∑Δ(ρWzWrAr). Moreover, results of CBUAA with slot-type CT indicate that the blowing mechanism at the CT front part creates an axial force Δ(ρWz2Az)<0 in the direction to the flow advance (see Figure 20(b), zone C). This observation suggests that the near casing flow is further accelerated. This is consistent with previous studies that the reinjected flow energizes the near casing flow. On the contrary, the flow bleeding mechanism creates an axial force Δ(ρWz2Az)>0 in the opposite direction to the flow advance. This is consistent since the near casing flow bled lost its axial velocity in order to be radially deflected into the CT slots.


### 6.5. Recommendation on the Useful Interpretation of the CT Flow Mechanisms Using the ESA Model

The analysis of the slot-type CT of the CBUAA compressor shows the limitation of the use of the histograms as a tool to investigate CT flow mechanisms. In fact, the histograms present a global overview of the forces being applied on a rotor tip control volume. They tend to “smooth” the local effect of CT flow mecha nisms (mainly both bleeding and blowing mechanisms). For an analysis of axisymmetric CT, histograms can be useful, whereas for slot-type CT the histograms prove to be not sufficient. Thus, one has to take care of the interpretation of the histograms. The main interest of the model comes from the analysis of the longitudinal evolutions of the axial forces since it permits to reveal the local CT flow mechanisms as well as the CT individual grooves’ contribution. So, it is recommended to perform both global and local analysis using the ESA model in order to further understand the flow mechanisms induced by CT.


## 6.1. Equation Analyzed

As it has been observed in previous studies, the CT interaction with the blade tip flow is mainly tridimensional and unsteady. In order to handle the results interpretation, the scope of the CT flow mechanisms’ study is restrained to the analysis of the balance of the axial momentum (7). Thus, results in terms of axial forces acting on the blade tip flow are expected. This equation is chosen since it describes the global flow across the rotor channel, specially its pressure rise.Here, we notice that the CFD computations presented before has been performed with steady as well as unsteady approaches. Thus, for the NASA Rotor 37, steady axial momentum balance have been computed. On the contrary, for the CREATE and CBUAA cases, instantaneous axial momentum balances have been obtained for instants equally spread over a respective known temporal period:(i)
for CREATE, 25 instants equally spread over one upstream stator blade S1 passing temporal period (TS1);(ii)
for CBUAA, 17 instants spread over one CT slot passing temporal period (TCT).For these two last compressors, the contribution of the time derivative term∂ρWz/∂t can be computed, which cannot be done for NASA Rotor 37. In order to permit comparisons between the three compressors, CREATE and CBUAA instantaneous axial momentum balances are time-averaged on their respective temporal periods. In other words, NASA steady axial momentum balance is here compared to CREATE and CBUAA time-averaged axial momentum balances.


## 6.2. Definition of the Control Volume

For all configurations, the control volume surrounds the rotor blade tip flow field as shown in Figure10. It is circumferentially delimited by the rotor blade pitch, axially extended just upstream and downstream of the rotor blade tip, and radially bounded by the 5th/25 to the 20th/25 grid layers modelling the rotor tip clearance. Therefore, due to the Chimera technique in the grooved configurations, the control volume is not attached to the shroud. In fact, the two last radial layers normal to the casing are used for interpolation between grooves and the blade passage blocks. These are not taken into account in the control volume since they can mislead the Navier-Stokes equations equilibrium.Views of the control volume located the NASA Rotor 37 tip clearance.
(a)
3D view(b)
Meridian view


## 6.3. Global Flow Mechanisms Acting on the Near Casing Flow

Results of the axial momentum balance for the SW and CT configurations are presented, respectively, for NASA Rotor 37 in Figure11, CREATE in Figure 12, and CBUAA in Figure 13. This kind of representation, as histograms, provides a macroscopic view of the different axial forces acting on the control volumes (i.e., the near casing flow). As a consequence, one has to take care of the results interpretation since histograms do not describe the local influence of CT. This will be further done in the next section titled “Longitudinal evolutions of the axial forces and CT contribution’’.Figure 11
NASA Rotor 37 with and without CT: balance of the axial momentum equation.Figure 12
CREATE-R2 with and without CT: balance of the axial momentum equation.Figure 13
CBUAA with and without CT: balance of the axial momentum equation.Before discussing the physical analysis, it is instructive to comment on the magnitude of the numerical terms of the balance. Results show that the numerical modelling residualsR(ρWz) tend to zero leading to a good precision of the balance of the equation. For the CREATE and CBUAA cases, it signifies that the time-averaged temporal derivative term (i.e., =-R(ρWz)¯) obtained from unsteady simulations converges to zero. This observation lends support that the unsteady calculations correctly reach a time-periodic state. Moreover, it can be observed that scalar artificial viscosity fluxes (Fλ2 and Fλ4) are insignificant in spite of the presence of a strong shock wave located in the blade tip region of the NASA Rotor 37 and the CBUAA.


### 6.3.1. SW Results

General Overview
For all SW configurations, results highlight the presence of four main axial forces being applied on the control volume:(i)
the axial pressure force∑Δ(PsAz) that corresponds to the axial flow pressure rise across the rotor.(ii)
The axial force due to the transport of the axial momentum across the longitudinal faces∑Δ(ρWz2Az). This term is similar to a dynamic pressure force that acts on the same surfaces than the pressure force. Its value relates to the flow axial acceleration or deceleration. In fact, by simply differentiating the terms between the axial outlet and inlet of the control volume (i.e., [(ρWz2Az)outlet-(ρWz2Az)inlet]) and assuming that Wz is mainly positive in the volume (which is not locally true due to the presence of the tip leakage, casing boundary layer,…), a negative (resp. positive) force corresponds to a deceleration (resp. acceleration) of the flow. The same reasonment can be done with Wz<0 that leads to the inverse conclusion. Results in Figures 11, 12, and 13 indicate that the incoming flow decelerates along the axial direction. This observation is coherent with the general aerodynamic of compressor since the casing boundary layer and the tip leakage flow greatly influence the incoming flow. Furthermore, the term ∑Δ(ρWz2Az) gives information on the axial flow blockage occurring in the control volume.(iii)
The axial shear force acting on the radial faces of the control volumes-∑Δ(τrzAr) induced by casing and blade tip boundary layers as well as the tip leakage flow.(iv)
The axial force due to the transport of the axial momentum across the radial faces∑Δ(ρWzWrAr). This term expresses the influence of the radial flow on the balance of axial forces. This force acts on the same faces than the axial shear force -∑Δ(τrzAr). Thus, it is believed that a change of one of these forces can strongly influence the other.All other terms are zero due to periodicity (forces acting onAθ faces) or to the radial thickness of the control volume. It is worth noticing that the forces applied on the near casing flow correspond to those observed in similar works done by Shabbir and Adamczyk [3].These axial forces can be distinguished depending on whether they are acting in the same direction of the flow (z>0) or in opposite direction of the flow (z<0). This information is given by the sign of the terms. A sketch describing the nomenclature of the axial forces acting on the control volume is presented in Figure 22. The results show for all compressors some coherent trends:(i)
the axial force due to∑Δ(PsAz) always acts in the opposite direction of the flow advance since the term is positive;(ii)
the axial forces due to∑Δ(ρWz2Az) and -∑Δ(τrzAr) always act in the same direction as the flow advance since the terms are negative.This observation is consistent with the general aerodynamic of axial compressor. The casing boundary layer (represented by-∑Δ(τrzAr)) and the upstream flow advance (described by ∑Δ(ρWz2Az)) balance the adverse pressure gradient (∑Δ(PsAz)) created by the axial compressor. If the boundary layer separates from the casing and/or if the flow advance is too low (like when the mass flow rate is too low), the pressure rise cannot be held, thus leading to an unstable operating flow condition.The different compressors differ by the magnitude of these axial forces and the direction of the∑Δ(ρWzWrAr).Influence of the Casing Geometry
Figures11, 12, and 13 clearly highlight the impact of the annulus shape (cylindrical versus conical). In fact, both CREATE and CBUAA compressors, which have a cylindrical casing, show similar results characterized by(i)
same relative balance between∑Δ(PsAz) and the opposite axial forces ∑Δ(ρWz2Az) and -∑Δ(τrzAr),(ii)
insignificant influence of∑Δ(ρWzWrAr).On the contrary, the NASA Rotor 37 that has a noncylindrical casing shows the occurrence of another axial force that is due to ∑Δ(ρWzWrAr). This force comes to contribute with the adverse pressure force ∑Δ(PsAz) to counteract the flow advance (i.e., z<0). This result is clearly obvious since the conical casing forces the flow to have a radial velocity component Wr<0. The near casing flow deflection is large enough to affect the balance of axial momentum, and it manifests itself through the axial force ∑Δ(ρWzWrAr).Influence of the Flow Regime
The comparison between CREATE and CBUAA (Figures12 and 13), which, respectively, works in subsonic and transonic flow regime, indicates no major difference in the axial momentum balance. Furthermore, the comparison between NASA Rotor 37 and CBUAA (Figures 11 and 13) proves that the shape of the rotor casing has a greater influence than the flow regime.


### 6.3.2. CT Results

General Overview
The results of the CT configurations in Figures11, 12, and 13 show the same type of forces acting on the blade tip flow than SW cases. However, the influence of CT is here revealed by changes of the axial forces contribution.For all compressors, the magnitude of pressure rise across the rotor∑Δ(PsAz) follows the predicted overall performances (see Figures 7, 8, and 9).Due to the radial flow exchanges between the grooves and the throughflow (see [9, 37, 44] for images of the flow exchanges), CT is expected to contribute to balance the axial momentum equation with additional terms. As previously assumed in the SW analysis, the radial flow exchanges might influence the axial momentum equation through the force ∑Δ(ρWzWrAr). In fact, results are in agreement with this hypothesis. They clearly show that the main impact of the CT, whether it is axisymmetric or not, is revealed by the changes of ∑Δ(ρWzWrAr) and the viscous force -∑Δ(τrzAr) due to the casing boundary layer. As previously announced, these two forces might influence each other since they both act on the same faces of the control volume. In fact, results reveal that the axial force ∑Δ(ρWzWrAr) changes proportionnaly to -∑Δ(τrzAr). These observations can be interpreted as follows: the radial fluid exchanges between the CT slots and the rotor channel tend to replace the casing boundary layer in its role to balance the adverse pressure gradient. This suggests that the Wr magnitude is a key point to the CT flow mechanisms. If the reasonment is pushed further, one explanation of the CT efficiency to improve the stall margin might come from its capability to perform an axial force due to the radial flow exchanges (through the axial force ∑Δ(ρWzWrAr)) more robust to the adverse pressure gradient than the casing boundary layer (through the axial force -∑Δ(τrzAr)). Moreover, since the shear axial force can induce a loss of efficiency, the use of radial flow exchanges as an alternative to casing boundary layer effect is very interesting in order to increase the global efficiency of the compressor rotor. As a consequence, the main difference between axisymmetric and non-axisymmetric CT can be evaluated depending on the Wr evolution.The results of the NASA Rotor 37 (Figure11) and CREATE (Figure 12) indicate that circumferential CT contribute with additional axial force ∑Δ(ρWzWrAr)<0 acting on the flow advance. As expected, it is the effect of the radial flow exchanges and this observation is in agreement with results of Shabbir and Adamczyk [3]. This adding term ∑Δ(ρWzWrAr)<0 is accompanied with a proportional decrease of the viscous axial force -∑Δ(τrzAr). It can be relied to the slight gain in efficiency observed in CREATE performances (Figure 8). As a matter of fact, the main role of axisymmetric CT is to turn or drive the force ∑Δ(ρWzWrAr) more along the main flow direction, that is, z>0. This is consistent with the general conclusions that previous studies have found that the reverse tip leakage flow through the tip gap got pushed into the grooves by the re-injecting flow (in main flow direction) from the grooves.On the contrary, the CBUUA slot-type CT contributes by additional force∑Δ(ρWzWrAr)>0 acting in opposite direction to the flow as well as ∑Δ(PsAz) (Figure 13). Moreover, the viscous axial force -∑Δ(τrzAr) has increased proportionally, thus explaining the slight decrease in efficiency observed in CBUAA performances (Figure 9). This first result reveals that slot-type CT behaves in a different manner than circumferential CT. However, as mentioned below histograms present a global and not a local influence of the CT. For this particular CT geometry, a local analysis has to be done to further understand the flow mechanisms as it will be seen in the next paragraphs.For both axi- and non-axisymmetric grooves, the balance change is made possible through the modification of the radial velocityWr induced by flow exchanges between grooves and the throughflow. It further supports that Wr is quite important for the flow stabilization in the tip region. This observation has to be confirmed by the local analysis of the CT flow mechanisms.


### 6.3.3. Unsteady Flow Mechanisms

Figures14 and 15 present, respectively, for CREATE and CBUAA the temporal evolution of the unsteady axial momentum balance. This analysis allows us to gain a quantitative and qualitative knowledge of the contribution of the time derivative term ∂ρWz/∂t.Figure 14
CREATE-R2 with and without CT: temporal evolution of the main terms of the unsteady axial momentum equation.Figure 15
CBUAA with and without CT: temporal evolution of the main terms of the unsteady axial momentum equation.Results of the CREATE-R2
For the CREATE case, unsteady results are presented over the upstream stator S1 passing temporal period. Plots in Figure14 show a sinusoidal time-dependent evolution of the time-derivative term influenced by the stator S1 blade passage frequency TS1. Its time-averaged value tends to zero, in agreement with results shown in Figure 12. Concerning its temporal fluctuations, amplitudes are particularly pronounced for the CT case with higher values of a factor 2 compared to the SW. Moreover, results show that ∂ρWz/∂t amplitudes are mainly driven by those of ∑Δ(PsAz) in the SW case and those of ∑Δ(ρWzWrAr) in the CT case. In fact, CT deadens the magnitude of the ∑Δ(PsAz) fluctuations by a factor 2 which is believed to greatly enhance the flow stability (because the pressure gradient opposes the flow). The increase in ∑Δ(ρWzWrAr) fluctuation, which drives the time-dependent term, proves that groove mechanisms strongly respond to upstream flow unsteadiness.Results of the CBUAA
Concerning the CBUAA case, unsteady results are presented over the CT slot passing temporal periodTCT. Plots in Figure 15 indicate that all terms present a “quasi-steady’’ behaviour, supporting that unsteadiness of the CT interaction with the throughflow is insignificant. This observation is especially confirmed by the zero magnitude over TCT of the time-derivative term ∂ρWz/∂t. This “quasi-steady’’ behaviour has already been observed by Lin et al. [22], and Ning and Xu [27]. Indeed, the unsteady time scale of the discrete treatment slots is therefore considerably small compared to that of the blade passage [27]; thus the unsteadiness of blade passage is weak. One reason to this phenomena might come from the ratio between the slot area and the area between the individual slots (which is very small here). Despite the high temporal resolution of the CT slot passing temporal period (TCT = 170 iterations), this particular test case fails to meet expectations in terms of analysis of unsteady flow interaction between the CT slots and the rotor channel. However, the effort is now focused on the analysis of local changes of the axial momentum balance in order to get more qualitative and quantitative information on the contribution of the CT slots.


## 6.3.1. SW Results

General Overview
For all SW configurations, results highlight the presence of four main axial forces being applied on the control volume:(i)
the axial pressure force∑Δ(PsAz) that corresponds to the axial flow pressure rise across the rotor.(ii)
The axial force due to the transport of the axial momentum across the longitudinal faces∑Δ(ρWz2Az). This term is similar to a dynamic pressure force that acts on the same surfaces than the pressure force. Its value relates to the flow axial acceleration or deceleration. In fact, by simply differentiating the terms between the axial outlet and inlet of the control volume (i.e., [(ρWz2Az)outlet-(ρWz2Az)inlet]) and assuming that Wz is mainly positive in the volume (which is not locally true due to the presence of the tip leakage, casing boundary layer,…), a negative (resp. positive) force corresponds to a deceleration (resp. acceleration) of the flow. The same reasonment can be done with Wz<0 that leads to the inverse conclusion. Results in Figures 11, 12, and 13 indicate that the incoming flow decelerates along the axial direction. This observation is coherent with the general aerodynamic of compressor since the casing boundary layer and the tip leakage flow greatly influence the incoming flow. Furthermore, the term ∑Δ(ρWz2Az) gives information on the axial flow blockage occurring in the control volume.(iii)
The axial shear force acting on the radial faces of the control volumes-∑Δ(τrzAr) induced by casing and blade tip boundary layers as well as the tip leakage flow.(iv)
The axial force due to the transport of the axial momentum across the radial faces∑Δ(ρWzWrAr). This term expresses the influence of the radial flow on the balance of axial forces. This force acts on the same faces than the axial shear force -∑Δ(τrzAr). Thus, it is believed that a change of one of these forces can strongly influence the other.All other terms are zero due to periodicity (forces acting onAθ faces) or to the radial thickness of the control volume. It is worth noticing that the forces applied on the near casing flow correspond to those observed in similar works done by Shabbir and Adamczyk [3].These axial forces can be distinguished depending on whether they are acting in the same direction of the flow (z>0) or in opposite direction of the flow (z<0). This information is given by the sign of the terms. A sketch describing the nomenclature of the axial forces acting on the control volume is presented in Figure 22. The results show for all compressors some coherent trends:(i)
the axial force due to∑Δ(PsAz) always acts in the opposite direction of the flow advance since the term is positive;(ii)
the axial forces due to∑Δ(ρWz2Az) and -∑Δ(τrzAr) always act in the same direction as the flow advance since the terms are negative.This observation is consistent with the general aerodynamic of axial compressor. The casing boundary layer (represented by-∑Δ(τrzAr)) and the upstream flow advance (described by ∑Δ(ρWz2Az)) balance the adverse pressure gradient (∑Δ(PsAz)) created by the axial compressor. If the boundary layer separates from the casing and/or if the flow advance is too low (like when the mass flow rate is too low), the pressure rise cannot be held, thus leading to an unstable operating flow condition.The different compressors differ by the magnitude of these axial forces and the direction of the∑Δ(ρWzWrAr).Influence of the Casing Geometry
Figures11, 12, and 13 clearly highlight the impact of the annulus shape (cylindrical versus conical). In fact, both CREATE and CBUAA compressors, which have a cylindrical casing, show similar results characterized by(i)
same relative balance between∑Δ(PsAz) and the opposite axial forces ∑Δ(ρWz2Az) and -∑Δ(τrzAr),(ii)
insignificant influence of∑Δ(ρWzWrAr).On the contrary, the NASA Rotor 37 that has a noncylindrical casing shows the occurrence of another axial force that is due to ∑Δ(ρWzWrAr). This force comes to contribute with the adverse pressure force ∑Δ(PsAz) to counteract the flow advance (i.e., z<0). This result is clearly obvious since the conical casing forces the flow to have a radial velocity component Wr<0. The near casing flow deflection is large enough to affect the balance of axial momentum, and it manifests itself through the axial force ∑Δ(ρWzWrAr).Influence of the Flow Regime
The comparison between CREATE and CBUAA (Figures12 and 13), which, respectively, works in subsonic and transonic flow regime, indicates no major difference in the axial momentum balance. Furthermore, the comparison between NASA Rotor 37 and CBUAA (Figures 11 and 13) proves that the shape of the rotor casing has a greater influence than the flow regime.


## 6.3.2. CT Results

General Overview
The results of the CT configurations in Figures11, 12, and 13 show the same type of forces acting on the blade tip flow than SW cases. However, the influence of CT is here revealed by changes of the axial forces contribution.For all compressors, the magnitude of pressure rise across the rotor∑Δ(PsAz) follows the predicted overall performances (see Figures 7, 8, and 9).Due to the radial flow exchanges between the grooves and the throughflow (see [9, 37, 44] for images of the flow exchanges), CT is expected to contribute to balance the axial momentum equation with additional terms. As previously assumed in the SW analysis, the radial flow exchanges might influence the axial momentum equation through the force ∑Δ(ρWzWrAr). In fact, results are in agreement with this hypothesis. They clearly show that the main impact of the CT, whether it is axisymmetric or not, is revealed by the changes of ∑Δ(ρWzWrAr) and the viscous force -∑Δ(τrzAr) due to the casing boundary layer. As previously announced, these two forces might influence each other since they both act on the same faces of the control volume. In fact, results reveal that the axial force ∑Δ(ρWzWrAr) changes proportionnaly to -∑Δ(τrzAr). These observations can be interpreted as follows: the radial fluid exchanges between the CT slots and the rotor channel tend to replace the casing boundary layer in its role to balance the adverse pressure gradient. This suggests that the Wr magnitude is a key point to the CT flow mechanisms. If the reasonment is pushed further, one explanation of the CT efficiency to improve the stall margin might come from its capability to perform an axial force due to the radial flow exchanges (through the axial force ∑Δ(ρWzWrAr)) more robust to the adverse pressure gradient than the casing boundary layer (through the axial force -∑Δ(τrzAr)). Moreover, since the shear axial force can induce a loss of efficiency, the use of radial flow exchanges as an alternative to casing boundary layer effect is very interesting in order to increase the global efficiency of the compressor rotor. As a consequence, the main difference between axisymmetric and non-axisymmetric CT can be evaluated depending on the Wr evolution.The results of the NASA Rotor 37 (Figure11) and CREATE (Figure 12) indicate that circumferential CT contribute with additional axial force ∑Δ(ρWzWrAr)<0 acting on the flow advance. As expected, it is the effect of the radial flow exchanges and this observation is in agreement with results of Shabbir and Adamczyk [3]. This adding term ∑Δ(ρWzWrAr)<0 is accompanied with a proportional decrease of the viscous axial force -∑Δ(τrzAr). It can be relied to the slight gain in efficiency observed in CREATE performances (Figure 8). As a matter of fact, the main role of axisymmetric CT is to turn or drive the force ∑Δ(ρWzWrAr) more along the main flow direction, that is, z>0. This is consistent with the general conclusions that previous studies have found that the reverse tip leakage flow through the tip gap got pushed into the grooves by the re-injecting flow (in main flow direction) from the grooves.On the contrary, the CBUUA slot-type CT contributes by additional force∑Δ(ρWzWrAr)>0 acting in opposite direction to the flow as well as ∑Δ(PsAz) (Figure 13). Moreover, the viscous axial force -∑Δ(τrzAr) has increased proportionally, thus explaining the slight decrease in efficiency observed in CBUAA performances (Figure 9). This first result reveals that slot-type CT behaves in a different manner than circumferential CT. However, as mentioned below histograms present a global and not a local influence of the CT. For this particular CT geometry, a local analysis has to be done to further understand the flow mechanisms as it will be seen in the next paragraphs.For both axi- and non-axisymmetric grooves, the balance change is made possible through the modification of the radial velocityWr induced by flow exchanges between grooves and the throughflow. It further supports that Wr is quite important for the flow stabilization in the tip region. This observation has to be confirmed by the local analysis of the CT flow mechanisms.


## 6.3.3. Unsteady Flow Mechanisms

Figures14 and 15 present, respectively, for CREATE and CBUAA the temporal evolution of the unsteady axial momentum balance. This analysis allows us to gain a quantitative and qualitative knowledge of the contribution of the time derivative term ∂ρWz/∂t.Figure 14
CREATE-R2 with and without CT: temporal evolution of the main terms of the unsteady axial momentum equation.Figure 15
CBUAA with and without CT: temporal evolution of the main terms of the unsteady axial momentum equation.Results of the CREATE-R2
For the CREATE case, unsteady results are presented over the upstream stator S1 passing temporal period. Plots in Figure14 show a sinusoidal time-dependent evolution of the time-derivative term influenced by the stator S1 blade passage frequency TS1. Its time-averaged value tends to zero, in agreement with results shown in Figure 12. Concerning its temporal fluctuations, amplitudes are particularly pronounced for the CT case with higher values of a factor 2 compared to the SW. Moreover, results show that ∂ρWz/∂t amplitudes are mainly driven by those of ∑Δ(PsAz) in the SW case and those of ∑Δ(ρWzWrAr) in the CT case. In fact, CT deadens the magnitude of the ∑Δ(PsAz) fluctuations by a factor 2 which is believed to greatly enhance the flow stability (because the pressure gradient opposes the flow). The increase in ∑Δ(ρWzWrAr) fluctuation, which drives the time-dependent term, proves that groove mechanisms strongly respond to upstream flow unsteadiness.Results of the CBUAA
Concerning the CBUAA case, unsteady results are presented over the CT slot passing temporal periodTCT. Plots in Figure 15 indicate that all terms present a “quasi-steady’’ behaviour, supporting that unsteadiness of the CT interaction with the throughflow is insignificant. This observation is especially confirmed by the zero magnitude over TCT of the time-derivative term ∂ρWz/∂t. This “quasi-steady’’ behaviour has already been observed by Lin et al. [22], and Ning and Xu [27]. Indeed, the unsteady time scale of the discrete treatment slots is therefore considerably small compared to that of the blade passage [27]; thus the unsteadiness of blade passage is weak. One reason to this phenomena might come from the ratio between the slot area and the area between the individual slots (which is very small here). Despite the high temporal resolution of the CT slot passing temporal period (TCT = 170 iterations), this particular test case fails to meet expectations in terms of analysis of unsteady flow interaction between the CT slots and the rotor channel. However, the effort is now focused on the analysis of local changes of the axial momentum balance in order to get more qualitative and quantitative information on the contribution of the CT slots.


## 6.4. Longitudinal Evolutions of the Axial Forces and CT Contribution

Since axial forces due to-∑Δ(τrzAr), ∑Δ(ρWzWrAr), and ∑Δ(ρWz2Az) highlight the change between SW and CT cases, it is interesting to investigate their respective axial evolution to further understand the flow interaction between CT and the throughflow. Figures 16 and 18 present the cumulative sum along the axial direction of these terms, respectively, for the NASA Rotor 37 and CBUAA. Figure 17 recalled from Legras et al. [37] the cumulative sum for the CREATE compressor of both steady and unsteady axial momentum balances. In the current study, only the unsteady curves are analyzed. Further comments on the difference between the steady and unsteady results are discussed in Legras et al. [37]. These diagrams provide information on the region where flow is highly constrained and on individual grooves contribution. The shaded bands denote the axial location of the casing grooves.Figure 16
NASA Rotor 37 with and without CT: cumulative sum in the axial direction of the principal terms involved in the axial momentum equation.CREATE-R2 with and without CT: cumulative sum in the axial direction of the principal terms involved in the axial momentum equation [37].
(a)(b)(c)Figure 18
CBUAA with and without CT: cumulative sum in the axial direction of the principal terms involved in the axial momentum equation.


### 6.4.1. SW Results

Results of all SW configurations show that the efforts evolve in amplitude in a nonlinear manner along the axial direction.The curve of viscous term-Δ(τrzAr) begins to grow near the leading edge. To relate these trends to the flow physics, Figure 19(a) presents at the mid-height of the CBUAA rotor tip gap, the flow field of -Δ(τrzAr). It can be observed that the growth of -Δ(τrzAr) can be relied firstly to the tip leakage vortex (referred as region A in Figure 19(a)) and secondly to the interaction of the detached shock with the shroud boundary layers (in a less pronounced manner).CBUAA: Flow field denoted at midgap by negative values of-Δ(τzrAr).
(a)
Smooth wall(b)
Grooved wallConcerning the curve of the termΔ(ρWz2Az), it begins to grow in negative value before the blade leading edge for the NASA Rotor 37 and CBUAA cases. Contour plots of Δ(ρWz2Az) in the CBUAA rotor tip gap presented. Figure 20(a) reveals that this growing trend is linked to the flow acceleration before the blade detached shock. Next, the curve slightly increases in positive values. This observation relates to the incoming flow deceleration due to its interaction with the tip leakage flow. Further downstream in the blade channel, the curve increases in negative values due to the presence of the tip leakage flow.CBUAA: Flow field denoted at midgap by negative values ofΔ(ρWz2Az).
(a)
Smooth wall(b)
Grooved wallThe curve of termΔ(ρWzWrAr) axially evolves in “quasi-’’ opposite magnitude of Δ(ρWz2Az) curve. In fact, Wr acts here as a “coefficient of amplification.’’ This observation is confirmed by Figure 21(a), where contour plots are very similar to those of Δ(ρWz2Az) in Figure 20(a).CBUAA: Flow field denoted at midgap by negative values ofΔ(ρWzWrAr).
(a)
Smooth wall(b)
Grooved wallFigure 22
Nomenclature and sketch of the axial forces acting on a control volume located in the rotor tip clearance.To resume, the SW results highlight strong axial force magnitude due to the tip leakage flow and, when a detached blade shock exists, due to the interaction between leakage flow and blade passage shock. However, the detached blade shock has a lower influence in term of force magnitude than the tip leakage flow. It explains why the axial momentum balances of CREATE (subsonic regime) and CBUAA (transonic regime) are very similar. Concerning the difference ofΔ(ρWzWrAr) and Δ(ρWz2Az) between NASA Rotor 37 and CBUAA, it can clearly be explained by the casing shape. In fact, the rapid growth of Δ(ρWzWrAr) and Δ(ρWz2Az) curves in NASA Rotor 37 corresponds to the location of the greatest slope of the casing form.


### 6.4.2. CT Results

Results of the CT configurations clearly highlight the local contribution of the slots.For both axi- and non-axisymmetric CT, the curve shapes of the viscous term-Δ(τrzAr) are very similar to those of the SW cases, except that CT delays further downstream the beginning of the growth rate. Since this term can be relied to the tip leakage flow, this observation can be explained by an attenuation of the spread of the leakage flow perpendicular to the blade. This is confirmed when comparing Figures 19(a) and 19(b). Moreover, referring to the CBUAA histograms in Figure 13, CT configuration shows higher magnitude of viscous term -Δ(τrzAr). In fact, contour plots in Figure 19(b) indicate the presence of a second region of high viscous axial force denoted B and located along the rear part of the CT slots. As seen by Legras et al. [44], this part of the slot is responsible for flow bleeding, whereas at the front part the recirculated flow is reinjected into the throughflow. This observation suggests that the flow bleeding mechanism induces strong shear layer that can increase viscous losses, thus explaining the efficiency losses observed with the use of the slot-type CT.The CT effects are even more highlighted by the curve ofΔ(ρWzWrAr). In fact, one can notice that axisymmetric CTs in NASA Rotor 37 and CREATE cases create an axial force Δ(ρWzWrAr)<0 in opposite direction to the flow advance. This explains the change observed in the previous histograms (Figures 11 and 12). This creation of axial force Δ(ρWzWrAr)<0 can be relied to the bleeding effect of the grooves since Wr>0 has been observed in previous studies [9, 37]. Moreover, this kind of analysis permits to identify the groove effectiveness depending on its capability to create this previous axial force, thus a bleeding effect. For the NASA Rotor 37, results show that the 2nd, 3rd, and 4th grooves are useful at the last stable point numerically seen for SW configuration. For the CREATE, Legras et al. [37] conclude that the 3rd groove is efficient at the last stable point numerically seen for SW configuration.Concerning the CBUAA slot-type CT, the bleeding effect along the rear part of the CT slots creates an axial force∑Δ(ρWzWrAr)<0 directed along the flow advance. This observation is confirmed by Figure 21(a) and is refered as zone D. As a matter of fact, the bleeding mechanisms of slot-type CT have a relatively similar effect than observed with circumferential CT. This is consistent with the general conclusions that previous studies have found that the reverse tip leakage flow through the tip gap got pushed into the grooves by the reinjecting flow (in main flow direction) from the grooves. On the contrary, the flow blowing mechanism along the front part of the CT induces axial force of ∑Δ(ρWzWrAr)>0 acting in opposite direction to the flow advance (refered as zone C in Figure 21(a)). In fact, the reinjected flow tends to radially deviate the incoming flow.Concerning the term ofΔ(ρWz2Az), curves mainly evolve in opposite magnitude to the term ∑Δ(ρWzWrAr). Moreover, results of CBUAA with slot-type CT indicate that the blowing mechanism at the CT front part creates an axial force Δ(ρWz2Az)<0 in the direction to the flow advance (see Figure 20(b), zone C). This observation suggests that the near casing flow is further accelerated. This is consistent with previous studies that the reinjected flow energizes the near casing flow. On the contrary, the flow bleeding mechanism creates an axial force Δ(ρWz2Az)>0 in the opposite direction to the flow advance. This is consistent since the near casing flow bled lost its axial velocity in order to be radially deflected into the CT slots.


## 6.4.1. SW Results

Results of all SW configurations show that the efforts evolve in amplitude in a nonlinear manner along the axial direction.The curve of viscous term-Δ(τrzAr) begins to grow near the leading edge. To relate these trends to the flow physics, Figure 19(a) presents at the mid-height of the CBUAA rotor tip gap, the flow field of -Δ(τrzAr). It can be observed that the growth of -Δ(τrzAr) can be relied firstly to the tip leakage vortex (referred as region A in Figure 19(a)) and secondly to the interaction of the detached shock with the shroud boundary layers (in a less pronounced manner).CBUAA: Flow field denoted at midgap by negative values of-Δ(τzrAr).
(a)
Smooth wall(b)
Grooved wallConcerning the curve of the termΔ(ρWz2Az), it begins to grow in negative value before the blade leading edge for the NASA Rotor 37 and CBUAA cases. Contour plots of Δ(ρWz2Az) in the CBUAA rotor tip gap presented. Figure 20(a) reveals that this growing trend is linked to the flow acceleration before the blade detached shock. Next, the curve slightly increases in positive values. This observation relates to the incoming flow deceleration due to its interaction with the tip leakage flow. Further downstream in the blade channel, the curve increases in negative values due to the presence of the tip leakage flow.CBUAA: Flow field denoted at midgap by negative values ofΔ(ρWz2Az).
(a)
Smooth wall(b)
Grooved wallThe curve of termΔ(ρWzWrAr) axially evolves in “quasi-’’ opposite magnitude of Δ(ρWz2Az) curve. In fact, Wr acts here as a “coefficient of amplification.’’ This observation is confirmed by Figure 21(a), where contour plots are very similar to those of Δ(ρWz2Az) in Figure 20(a).CBUAA: Flow field denoted at midgap by negative values ofΔ(ρWzWrAr).
(a)
Smooth wall(b)
Grooved wallFigure 22
Nomenclature and sketch of the axial forces acting on a control volume located in the rotor tip clearance.To resume, the SW results highlight strong axial force magnitude due to the tip leakage flow and, when a detached blade shock exists, due to the interaction between leakage flow and blade passage shock. However, the detached blade shock has a lower influence in term of force magnitude than the tip leakage flow. It explains why the axial momentum balances of CREATE (subsonic regime) and CBUAA (transonic regime) are very similar. Concerning the difference ofΔ(ρWzWrAr) and Δ(ρWz2Az) between NASA Rotor 37 and CBUAA, it can clearly be explained by the casing shape. In fact, the rapid growth of Δ(ρWzWrAr) and Δ(ρWz2Az) curves in NASA Rotor 37 corresponds to the location of the greatest slope of the casing form.


## 6.4.2. CT Results

Results of the CT configurations clearly highlight the local contribution of the slots.For both axi- and non-axisymmetric CT, the curve shapes of the viscous term-Δ(τrzAr) are very similar to those of the SW cases, except that CT delays further downstream the beginning of the growth rate. Since this term can be relied to the tip leakage flow, this observation can be explained by an attenuation of the spread of the leakage flow perpendicular to the blade. This is confirmed when comparing Figures 19(a) and 19(b). Moreover, referring to the CBUAA histograms in Figure 13, CT configuration shows higher magnitude of viscous term -Δ(τrzAr). In fact, contour plots in Figure 19(b) indicate the presence of a second region of high viscous axial force denoted B and located along the rear part of the CT slots. As seen by Legras et al. [44], this part of the slot is responsible for flow bleeding, whereas at the front part the recirculated flow is reinjected into the throughflow. This observation suggests that the flow bleeding mechanism induces strong shear layer that can increase viscous losses, thus explaining the efficiency losses observed with the use of the slot-type CT.The CT effects are even more highlighted by the curve ofΔ(ρWzWrAr). In fact, one can notice that axisymmetric CTs in NASA Rotor 37 and CREATE cases create an axial force Δ(ρWzWrAr)<0 in opposite direction to the flow advance. This explains the change observed in the previous histograms (Figures 11 and 12). This creation of axial force Δ(ρWzWrAr)<0 can be relied to the bleeding effect of the grooves since Wr>0 has been observed in previous studies [9, 37]. Moreover, this kind of analysis permits to identify the groove effectiveness depending on its capability to create this previous axial force, thus a bleeding effect. For the NASA Rotor 37, results show that the 2nd, 3rd, and 4th grooves are useful at the last stable point numerically seen for SW configuration. For the CREATE, Legras et al. [37] conclude that the 3rd groove is efficient at the last stable point numerically seen for SW configuration.Concerning the CBUAA slot-type CT, the bleeding effect along the rear part of the CT slots creates an axial force∑Δ(ρWzWrAr)<0 directed along the flow advance. This observation is confirmed by Figure 21(a) and is refered as zone D. As a matter of fact, the bleeding mechanisms of slot-type CT have a relatively similar effect than observed with circumferential CT. This is consistent with the general conclusions that previous studies have found that the reverse tip leakage flow through the tip gap got pushed into the grooves by the reinjecting flow (in main flow direction) from the grooves. On the contrary, the flow blowing mechanism along the front part of the CT induces axial force of ∑Δ(ρWzWrAr)>0 acting in opposite direction to the flow advance (refered as zone C in Figure 21(a)). In fact, the reinjected flow tends to radially deviate the incoming flow.Concerning the term ofΔ(ρWz2Az), curves mainly evolve in opposite magnitude to the term ∑Δ(ρWzWrAr). Moreover, results of CBUAA with slot-type CT indicate that the blowing mechanism at the CT front part creates an axial force Δ(ρWz2Az)<0 in the direction to the flow advance (see Figure 20(b), zone C). This observation suggests that the near casing flow is further accelerated. This is consistent with previous studies that the reinjected flow energizes the near casing flow. On the contrary, the flow bleeding mechanism creates an axial force Δ(ρWz2Az)>0 in the opposite direction to the flow advance. This is consistent since the near casing flow bled lost its axial velocity in order to be radially deflected into the CT slots.


## 6.5. Recommendation on the Useful Interpretation of the CT Flow Mechanisms Using the ESA Model

The analysis of the slot-type CT of the CBUAA compressor shows the limitation of the use of the histograms as a tool to investigate CT flow mechanisms. In fact, the histograms present a global overview of the forces being applied on a rotor tip control volume. They tend to “smooth” the local effect of CT flow mecha nisms (mainly both bleeding and blowing mechanisms). For an analysis of axisymmetric CT, histograms can be useful, whereas for slot-type CT the histograms prove to be not sufficient. Thus, one has to take care of the interpretation of the histograms. The main interest of the model comes from the analysis of the longitudinal evolutions of the axial forces since it permits to reveal the local CT flow mechanisms as well as the CT individual grooves’ contribution. So, it is recommended to perform both global and local analysis using the ESA model in order to further understand the flow mechanisms induced by CT.


## 7. Conclusion

This paper has presented a generalized method of the Shabbir and Adamczyk approach for uncovering the flow mechanisms induced by casing treatment (CT) using CFD numerical simulations. Compared to Shabbir’s approach, this method permits to analyze any CT geometry, thanks to a budget analysis of each of the Navier-Stokes (un)steady momentum equations. In the current paper, only the equilibrium of the axial momentum equations is analyzed since it describes the global flow across the rotor channel, specially its pressure rise. The novel method has been used to investigate the flow interaction between axi and non-axisymmetric CT and the throughflow for three different axial compressors. They differ by the flow regime (subsonic/transonic), the untreated casing shape (cylindrical/conical), and the CT implemented (axi-/non axisymmetric CT).The comparison of the three compressors indicates that axisymmetric CT is characterized by its flow bleeding effect that creates an axial forceΔ(ρWzWrAr)<0. It is shown that the benefit of this force on the flow stability is more dependant on the casing shape than the flow regime. Concerning slot-type CT, it is characterized by both a bleeding effect at the rear part, which creates an axial force Δ(ρWzWrAr)<0, and a blowing effect at the front part that energizes the near casing flow by creating axial force Δ(ρWz2Az)<0. However, its strong bleeding effect is responsible of additional viscous losses explaining the efficiency losses observed with the use of slot-type CT. These conclusions are rather difficult to be stated from the only histograms. Thus, it is strongly recommended to focus on the axial distribution of the axial forces to further understand the interaction between the CT and the flow. Furthermore, since both axi- and non-axisymmetric CT create an axial force Δ(ρWzWrAr)<0, it lends support that the radial velocity component Wr is of a main interest for the stability of the tip flow.The current analysis methodology proved its potential to help in the designing of CT. However, there are still many questions on the understanding of the CT mechanisms that the model can answer. The main perspectives concern(i)
the analysis of tangential momentum: since the tangential force is the main force of compressors, it could be interesting to investigate the influence of CT,(ii)
the analysis of the balance at part speed,(iii)
the resolution of the momentum equation in the streamline reference frame rather than cylindrical.Finally, it can also be easily used to uncover the flow mechanisms in order to design a large panel of flow passive control devices: film cooling, and boundary layer aspiration in the turbomachinery field as well as for aircraft or helicopter applications.



---

*Source: 975407-2012-06-14.xml*