Next Article in Journal
Flow Characterization of the UTSA Hypersonic Ludwieg Tube
Next Article in Special Issue
Experimental Firing Test Campaign and Nozzle Heat Transfer Reconstruction in a 200 N Hybrid Rocket Engine with Different Paraffin-Based Fuel Grain Lengths
Previous Article in Journal
A Tool for the Design of Turbomachinery Disks for an Aero-Engine Preliminary Design Framework
Previous Article in Special Issue
A Simplified Thermal Analysis Model for Regeneratively Cooled Rocket Engine Thrust Chambers and Its Calibration with Experimental Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Velocity Increment on Incidence Angle near the Leading Edge of the Compressor Cascade

1
School of Energy and Power Engineering, Beihang University, Beijing 100083, China
2
National Key Laboratory of Science & Technology on Aero-Engine Aero-Thermodynamics, Beihang University, Beijing 100083, China
3
Research Institute of Aero-Engine, Beihang University, Beijing 100083, China
*
Author to whom correspondence should be addressed.
Aerospace 2023, 10(5), 461; https://doi.org/10.3390/aerospace10050461
Submission received: 14 February 2023 / Revised: 6 May 2023 / Accepted: 15 May 2023 / Published: 16 May 2023
(This article belongs to the Special Issue Fluid-Dynamics and Heat Transfer in Aerospace Propulsion Systems)

Abstract

:
The geometry of a compressor leading edge has an important effect on the aerodynamic performance at an off-designed incidence angle. The current geometric design methods of the leading edge are usually developed based on the flow characteristics at the designed incidence angle. However, few research focuses on the quantitative rules of the leading edge flow characteristics at the off-designed incidence angle in a compressor cascade. This situation restricts the further optimization and development of the leading edge geometry design method. In this paper, starting from the research of a potential cascade theory, the singularity point, where the surface velocity approaches infinity in the leading edge region, is eliminated by applying the characteristic that the ratio of the velocity increasement on the incidence angle in the plate cascade and the isolated plate flow is finite. Secondly, the equivalent pitch lengths based on 1/cos(β) and VI caused by a diffuser deceleration in the cascade passage were employed to correct the effect of the stagger angle. Finally, by introducing the isolated flow around the thick airfoil and considering the influence of the camber line geometry, a model of the variation of the surface velocity near the leading edge under the off-designed incidence angle, named the velocity increment on incidence angle, is derived from any compressor cascade. Hence, the relation between the off-designed incidence angle and the designed incidence angle of the surface velocity in a cascade blade is established, and it depends only on the geometrical parameters. Through a verification using numerical calculations and experimental measurement, the explicit formula for the velocity increment on incidence angle proposed in this paper has high precision near the leading edge.

1. Introduction

The geometry of a leading edge has an important effect on the aerodynamic performance of the axial flow compressor blade at an off-designed incidence angle, such as that showed by Carter in his experimental study whereby the smaller the blade leading edge size, the smaller the blade loss and the larger the operational incidence range [1]. Unlike the circular leading edge, which is widely used, Walraevens and Cumpsty confirmed that the boundary layer behind the elliptic leading edge was healthier, and they indicated that the difference was due to the laminar separation bubbles caused by the velocity spike around the leading edge [2]. With the deepening of the research, Shi et al. indicated that the entire integral entropy generation in the boundary layer increased sharply when the laminar separation bubble moved upstream to the leading edge [3], and Lambert et al. observed that the size of the laminar separation bubble increased with a larger incidence angle [4]. Given the importance of a leading edge spike, Goodhand et al. quantified the spike intensity [5], and they presented a continuous-curvature leading edge to eliminate a spike at the designed incidence angle. Following this research, Xiang et al. [6] discussed that the continuity of the first derivative of curvature can restrain the separation of the boundary layer and delay the transition process. Liu et al. [7] and Song et al. [8] also confirmed that a well-designed continuous-curvature leading edge can eliminate spikes at the designed incidence angle and improve the blade performance by using the CST method and B-spline method, respectively.
The current problem of the blade leading edge research is how to generate the best geometry with the largest operational incidence range. The current geometry generation methods are usually applied for the single designed flow condition, then, the off-designed aerodynamic performance is obtained by a numerical method, and finally, the optimal geometry with the targeted composite performance could be selected after several iterations [9]. Many optimization methods have been proposed, such as, a genetic algorithm [10], particle swarm optimization [11], a surrogate model with a response surface [12], deep learning [13] or reinforcement learning [14], et al. However, in those methods, the off-designed aerodynamic performance appears passively, and more resources are consumed in the iterations. A multipoint inverse design method posed by Michael is another alternative road, and it considers the ideal velocity distribution on the blade surface under different working conditions [15,16]. However, their generation of airfoil still relies on iterations.
Unlike the increasingly complex optimization methods, the theoretical innovation could simplify the blade generation process. However, for the blade under off-designed incidence angles, there are no explicit quantitative rules of the flow characteristics for a universal compressor cascade. If the off-designed performance could be actively embodied in the design procedure by utilizing the aerodynamic relation between the off-designed and the designed flow condition, the iterations in the process of leading edge geometry generation can be reduced or even avoided. However, in current studies, the surface velocity under the off-designed incidence angle is usually obtained from CFD results or implicit potential flow solution. The rise of the CFD method has brought about a large number of flow field details, but these flow field details cannot be used to explore the general rules under different flow conditions. Before CFD was applied in engineering, the analytical conformal transformation method based on potential flow theory was widely used as a flow analysis tool in blade profile design [17,18]. However, the potential theories for cascade flow are all implicit equations and have remained stagnant for decades.
To understand the basic mathematical rules of the cascade flow field in more detail, the potential flow theory with explicit equations for cascade should be proposed. Recently, along with the development in applied mathematics, the integral problem with boundary conditions was applied to the cascade flows. In this way, Baddoo and Ayton obtained an explicit solution for the cascade flows with a small stagger angle, small turning angle, thin thickness, and sharp trailing edge [19]. The explicit solution can directly reflect the mathematical rules of various geometric parameters on the flow field in the compressor cascade, but the accuracy is poor near the blade leading edge because of the infinity tendency at the leading edge point. Moreover, the explicit solution should be extended to a real cascade blade with a larger stagger angle and larger thickness.
This paper will study the steady, incompressible, non-viscous and non-rotational potential cascade flow theory, supplemented with numerical and experimental verification, in order to find the aerodynamic relation between the off-designed and the designed flow condition. Finally, the explicit mathematical rules of the surface velocity increasement on the incidence angle for universal compressor cascades would be presented.

2. Potential Analysis of Simplified Cascade

For the two-dimensional cascade of general incompressible axial compressor studied in this paper, different values can be selected for pitch length, t, inlet metal angle, β1, outlet metal angle, β2, camber line geometry and thickness distribution when the chord length, c, is specified. The cascade geometry and parameter definitions are shown in Figure 1, in which α1 is the inlet flow angle, α2 is the outlet flow angle, β is the stagger angle, and the incidence angle is defined as i = α1β. When the velocity on the blade surface is represented by a vector, the positive direction flows clockwise. To consider the high degree of freedom in cascade geometry, the study of the general rules of universal cascade should first simplify the two-dimensional cascades properly, and then add the influence of different geometry parameters into the simplified model.
In this paper, the study of plate cascades will be firstly introduced and then extended to elliptic straight cascades with stagger angles. Finally, the influence of actual thickness distribution and camber line geometry will be considered to form a model of the variation of the surface velocity near the leading edge under the off-designed incidence angle, which is derived from any compressor cascade. Furthermore, the flow fields calculated with the conformal transformation method would be treated as the true value, as they have mathematical accuracy in potential flow [20].

2.1. Potential Analysis of Plate Cascade

From Baddoo and Ayton’s study [19], when the trailing edge always conforms to the Kutta condition, the explicit expression of the influence of the incidence angle on the axial velocity of the blade surface is shown in the following equation:
u i ± x = i U e π t s i n h π 1 x / t s i n h π 1 + x / t
where, the value of x is −1 to 1, hence, the chord length is always 2. Since the influence of the stagger angle and thickness was ignored in the derivation procedure, Equation (1) is essentially the relationship of the surface velocity of the cascade blade on the incidence angle, the pitch length, and the axial position. When the pitch length tends to infinity, Equation (1) is equivalent to the incidence term of the flow around an isolated plate [21]:
u i ± x = i U 1 x 1 + x
Furthermore, the error of Equation (1) in a general cascade is  O ( ε 2 ) , where  ε  is approximately equal to the ratio of thickness height to chord length.
In this paper, consider the definition of the velocity increment on incidence angle, VI, as the derivative of the velocity to incidence angle, as described in the equation:
V I = 1 U d u d i
It shows the change in blade surface per unit incidence angle (in radians). Hence, VI of the plate cascade follows the equation:
V I p l a t e , c a s c a d e = e π t s i n h π 1 x / t s i n h π 1 + x / t
and VI of the isolated plate follows the equation:
V I p l a t e , i s o l a t e d = 1 x 1 + x
The VI distribution of plate cascade with axial position x under different pitch length is shown in Figure 2. For VI at any pitch length, the value is maximum at the leading edge point and tends to infinity, and then decreases rapidly at the axial position toward the trailing edge. When the pitch length tends to infinity, the value of Equation (4) is completely equivalent to the value of Equation (5).

2.2. Potential Analysis of Elliptic Straight Cascade

Evidently, VI at the leading edge point of any thick blade will not take an infinite value, which shows the application limitations of the analytical explicit expressions shown in Equations (4) and (5). However, the ratio of VI in the plate cascade to VI on an isolated plate,  V I p l a t e , c a s c a d e / V I p l a t e , i s o l a t e d , maintains to be finite, as shown in Figure 3. Hence, the defect whereby the current explicit analytical solution cannot be applied to the leading edge region has been solved.
For the flow around an ellipse with a major axis of 2 on the x axis, the surface velocity under the Kutta condition follows the equation:
U θ = 2 U sin θ 2 cos θ 2 i ( A R + 1 ) A R 2 sin 2 θ + cos 2 θ
in which  θ  is the angular parameter (the axial position x = cos( θ )) and AR is the axial ratio of the ellipse. Hence, VI of the isolated ellipse could be expressed as:
V I e l l i p s e , i s o l a t e d = 2 sin θ 2 sin θ 2 i A R + 1 A R 2 sin 2 θ + cos 2 θ
and VI of the elliptic straight cascades could be established from the flow around the isolated ellipse by utilizing the relationship between plate cascade and flow around the isolated plate, which is expressed as:
V I e l l i p s e , c a s c a d e = V I p l a t e , c a s c a d e V I p l a t e , i s o l a t e d · V I e l l i p s e , i s o l a t e d
The comparison between the value obtained from Equation (8) for an elliptical cascade with t = 2 and the true value is shown in Figure 4. In Figure 4a, the inferred value is almost the same to the true value; it shows the availability of Equation (8). Figure 4b shows the specific difference between the two values with different axial ratios; it indicates that the leading edge region has the largest difference, and the closer the blade surface position is to the trailing edge, the smaller is the difference. At the same time, with the same pitch length, the smaller the axial ratio of the ellipse, the larger the difference. However, the axial ratio does not affect the difference at the leading edge point. The errors between the inferred value and true value are not larger than  O ( ε 2 ) ; therefore, the assumption made by Equation (8) is reasonable.

2.3. Potential Analysis of Elliptic Straight Cascade with Stagger Angle

For the elliptic cascade with t = 2, AR = 10, the effect of the stagger angle is shown in Figure 5, while the inlet flow angles change, following the stagger angles to keep the front stagnation point on the leading edge point. Figure 5a indicates that the larger the stagger angle, the larger the value of VI near the leading edge region, and the larger the stagger angle is, the larger the difference between VI values of the pressure surface and suction surface in the blade body. Evidently, this is the influence of the change of inlet flow angle under the large stagger angle, which leads to the relative expansion and deceleration of the flow passage. The distribution of VI along the blade surface forms a closed ring, and the suction surface branch is always located in the lower half of the closed ring, that is, VI is always positive near the leading edge of the suction surface; the back half is negative, and the pressure surface is always positive. While the surface velocity is scalar, the distribution of VI is plotted in Figure 5b. In the figure, the upper part is the suction surface side, and the lower part is the pressure surface side. That is, with the increase in the incidence angle, the absolute velocity of the pressure surface always decreases, while the velocity of the suction surface increases in the leading edge region and decreases in the latter part of the blade.

2.3.1. Equivalent Pitch Length on Stagger Angle

The method proposed by Baddoo and Ayton to correct the stagger angle is to consider only the change of chord line projection [19], which has a large error in the actual cascade. However, using the equivalent pitch length in:
t e q = t cos β
we can correct the influence of the stagger angle near the leading edge, as shown in Figure 6. In that figure, the corrected value of VI on the leading edge point matches the true value well, and the correction error is also smaller than  O ( ε 2 ) . Furthermore, the method proposed by Baddoo and Ayton [19] introduces the opposite correction. Therefore, the correction of Equation (4) changes to:
V I p l a t e , c a s c a d e = e π t e q s i n h π 1 x / t e q s i n h π 1 + x / t e q
and the correction of Equation (8) changes to:
V I s t a g g e r , c o r r = e π t e q s i n h π 1 x / t e q s i n h π 1 + x / t e q · 1 + x 1 x · 2 sin θ 2 sin θ 2 i A R + 1 A R 2 sin 2 θ + cos 2 θ

2.3.2. Diffuser Deceleration on Stagger Angle

The outlet velocity of cascade is only related to the inlet and outlet flow angle, and the three are constrained by the flow conservation equation:
U o u t = U i n c o s ( α 1 ) c o s ( α 2 )
under the incompressible assumption. Due to the Kutta condition, the rear stagnation point of the flow is always at the trailing edge point, so the outlet flow angle almost does not change with the inlet flow angle within a certain range. Therefore, when the inlet flow angle changes, the diffuser degree of the flow passage changes. In this case, the derivative of the outlet velocity changing with the inlet angle of attack satisfies:
d U o u t U i n d i = V I T E = s i n ( α 1 ) c o s ( α 2 )
namely, the expression of VI at the trailing edge point. Considering the vector representation of blade surface velocity, the  V I T E  of pressure surface and suction surface should be equal in magnitude and opposite in sign.
In the body region from the leading edge to the trailing edge, the average local flow angle is between the inlet and outlet flow angles. As the turning of airflow needs the driver of pressure difference and the velocity increasement represents the velocity difference driven by the pressure difference, the ratio of the flow turning angle to the inlet flow angle in the body region could be described with:
T x = 1 1 x V I β , c o r r ( s ) d s 1 1 V I β , c o r r ( s ) d s
In this definition, at the leading edge, T(−1) = 1, and at the trailing edge, T(1) = 0. Therefore, VI caused by a diffuser deceleration, caused by the changing of the inlet flow angle, satisfies:
V I d i f f u s i o n = [ 1 T x ] s i n ( α 1 ) c o s ( α 2 )
Through the above analysis, the expression of blade surface VIβ of an elliptic cascade with a stagger angle is shown as:
V I β = V I s t a g g e r , c o r r + V I d i f f u s i o n
When t = 2, AR = 10, β = 30°, the calculated value of VI in the elliptic cascade has a good matching with the true value, as shown in Figure 7a. Figure 7b shows the distribution of the estimated error of Equation (16) under different stagger angles. It can be seen that the error of Equation (16) gradually increases with the increase in the stagger angle, but the error is always at least one order of magnitude smaller than the value of VI when the stagger angle is not too large. Therefore, Equation (16) can be used as the estimation of VI of the cascades with a stagger angle.
It can be seen that diffuser deceleration has a greater impact on the blade body, while the impact on the leading edge is not significant. Therefore, Equation (16) proposed in this paper has been able to achieve a high-precision estimation of the leading edge region.

3. Extension in Actual Compressor Cascade

3.1. Coordinate System on the Blade Surface

In the above studies, the abscissa of VI is the axial position, while the leading and trailing edge points of the ellipse are located on the (−1, 0) and (1, 0), respectively. When the coordinates are out of the range [−1, 1], Equation (16) would be unsolvable. In the actual blade with a significant curvature in the camber line, the normalization position value should be in the section, [0, 1], strictly, to adapt the x-axis value of [−1, 1] in Equation (16). However, the x-axis value or the normalized chordwise value of the points near the leading edge could be smaller than 0, and the normalized arc length coordinate value then has a discontinuity of derivative near the leading edge. Therefore, the camber projection normalization method is used to define the normalized position of any blade position. As shown in Figure 8, the arc length on the blade surface between any point and the leading edge point is s; the arc length on the camber line between the projection point and the leading edge point is sp; the arc length on the blade surface between the leading edge point and the trailing edge point is S0; and the arc length on the camber line between leading edge point and the trailing edge point is scamber. Hence, the normalized camber projection coordinate, sc, is described as sc = sp/scamber, and the normalized surface arc length coordinate is written as s/S0. As the normalized surface arc length coordinate has a higher proportion near the leading edge, it is used in the experiment to show the details in the field-of-view area.
The normalization process of the camber projection of the blade surface points is actually the process of solving blade thickness distribution and camber geometry. At this time, two points on the suction surface and pressure surface with the same scamber are symmetric along the coordinate axis of the blade thickness distribution. On the thickness distribution curve, the symmetry of the suction surface and pressure surface can be maintained from the leading edge point to the trailing edge point. Therefore, the camber projection normalization method has the best symmetry on the whole blade.

3.2. Effect of Thickness Distribution

The thickness distribution curve of a real blade is not strictly elliptic, and the curvature radius at the front and trailing edges is usually larger than the curvature radius at the end of the ellipse’s major axis when the maximum thickness is the same. Based on the curvature radius formula of the ellipse, r = a/AR2, in which r is the curvature radius at the end of the major axis and a is the length of the major axis, a definition of the equivalent axial ratio based on chord length and the curvature radius at leading edge point is shown as:
A R r = c / 2 r L E
Moreover, another definition of the equivalent axial ratio based on chord length and the maximum thickness is shown as:
A R t = c / H m a x
In the potential flow, VI near the leading edge of a straight cascade with any thickness distribution should be similar to that of an elliptic cascade with the axial ratio equal to ARr, and VI in the blade body region should be similar to that of an elliptic cascade with the axial ratio equal to ARt. The black line in Figure 9a shows a straight cascade with a stagger angle of 30° and a maximum thickness of 5% chord length, but with a leading edge curvature radius of 0.5% chord length. In this case, its ARr is 10 and ARt is 20. Furthermore, the standard ellipse with AR equal to 20 is shown as the blue line and the ellipse with AR equal to 10 is shown as the red line. The comparison of VI among the three geometries is shown in Figure 9b, which indicates that the elliptical cascade with ARr is the most similar to the freedom thickness straight cascade in the leading edge region, while the elliptical cascade with ARt is the most similar to the freedom thickness straight cascade in the body region. When β changes in Figure 9c, the three leading edges keep the relative trends. Corresponding to the research content of this paper, in the actual cascade, ARr should be used as the equivalent axial ratio to calculate VI near the leading edge.

3.3. Effect of Turning Angle of the Camber Line

The curvature of the actual camber line of the blade will affect the diffuser strength in the blade body region. If the tangent angle of the camber line at any position is defined as βx, then according to Equations (13)–(16), VI value caused by diffuser at any position of the cascade is:
V I d i f f u s i o n = 1 x V I β , c o r r ( s ) d s 1 1 V I β , c o r r ( s ) d s · s i n ( α 1 ) c o s ( α 1 + ( β x β 1 ) α 2 α 1 β 2 β 1 )
However, VIdiffusion is very small near the leading edge, meaning that the effect of the channel diffuser deceleration is not significant near the leading edge. Therefore, the influence of camber geometry on VI near the leading edge is mainly reflected in the change of the cascade stagger angle, and VIdiffusion could be ignored when the research is focused on the leading edge region.

4. Numerical and Experimental Verification

The above analysis of VI starts from the potential flow solution of the plate cascade, then studies the VI distribution rule of elliptical cascade, and finally builds a VI distribution expression of an arbitrary compressor cascade by adding the influence of the thickness distribution and the camber line geometry of real blade. Due to the influence of the boundary layer, the error of the potential flow analysis results is inevitable in practice. However, the boundary layer near the blade leading edge is usually thin and has a weak influence on the mainstream flow field. Therefore, the analysis results in this paper are of high practical value. In this section, the expression of VI obtained in this paper will be verified by numerical calculations and experimental measurement.

4.1. Numerical Verification

The cascade bases on a typical compressor stator blade mid-span profile was used to analyze the difference between the above theoretical equation and CFD calculation near the blade leading edge, whose inlet and outlet metal angles are 51° and 6°, respectively, and the solidity is 1.44, as shown in Figure 10a. MISES was employed as the CFD tool to solve the flow field of cascade near the designed inlet flow angle (45°). In order to calculate VI, two incidence angles of ±0.5° were used, and the difference of blade surface velocity between them was used to calculate the CFD solution of VI. As the central difference scheme was adopted, the value of VI reached the second-order accuracy and met the test requirements. Figure 10b shows the difference in the values of VI between the CFD solution and the theory estimation from Equations (11) and (16). It can be seen that  V I β  from Equation (16) can effectively predict VI in a large range of the blade surface, and there is a certain error only in the middle area of the blade suction surface. The error in the middle area of the suction surface is mainly due to the fact that the main flow velocity is higher than the incoming flow velocity, so the denominator of the relative velocity should be the local main flow velocity rather than the incoming flow velocity.  V I s t a g g e r , c o r r  from Equation (11) can only accurately predict the increment of angular attack velocity near the leading edge, as shown in Figure 10c. For the study of the leading edge region,  V I s t a g g e r , c o r r  can also be used to effectively analyze the distribution of VI and benefit from its very simple calculation process.

4.2. Experimental Verification

In order to further verify the accuracy of the theoretical analysis results and numerical calculation results, high spatial resolution PIV measurements have been carried out in the leading edge area to verify the changing process of blade surface velocity at different incidence angles. The experiment involved in this paper were carried out in the low-speed planar cascade facility located in Beihang University with a high-quality inlet environment (less than 1.6% turbulence) and a 420 × 120 mm2 outlet cross-sectional area. Basic parameters of the cascade facility are shown in Table 1, in which Reynolds number is calculated based on chord length and inlet velocity.
For the PIV measurement setup, a 140 mJ/Pulse Nd:YAG laser was employed, and the light sheet thickness was approximately 1 mm. ImagerProX-4M camera with a Nikon Micro 105 mm lens and a resolution of 2048 × 2048 pixels (at 4 Hz recording rate) were used for image capturing, and the magnification factor of the lens was approximately 1.0. DEHS and a Laskin nozzle-based particle generator were used to generate the seeding particles. The time interval of two laser beams was 2.5 μs, and the averaged displacement of the particles is 4.5 pixels, while the maximum displacement of the particles is 9.8 pixels. Figure 11a shows the schematic of the blade and the camera range of the PIV measurement. The physical size of the camera field-of-view shown in Figure 11 is 14.5 × 14.5 mm2; it measured 0~8% in chord length near the leading edge. For the captured particle images shown in Figure 11b, the particle diameter is approximately 1.5–3.0 pixels, the particle distribution density is approximately 0.021–0.048 particles/pixels, and about 600 snapshots were obtained for each incidence angle condition in the experiment. While processing the pre-processed transformed images, three passes of 32 × 32 pixels and then three passes of 12 × 12 pixels interrogation windows are employed to achieve the WIDIM algorithm with a 50% overlap. As the displacement uncertainty for this algorithm is 0.05 pixels, the measured velocity uncertainty is about 1.1%. The details of the post-processing procedure was based on the method by Liu and Xu [22].
In the experiment, velocity fields in the area of interest under six different incidence angles were obtained, and the velocity at the main flow boundary was obtained by identifying the boundary layer, as shown in Figure 12. Since the boundary layer near the leading edge at the designed incidence angle is thin, the velocity at the main flow boundary is essentially equivalent to the isentropic velocity of the blade surface. As can be seen from Figure 12, local separation bubbles appear in the leading edge area at 9.5° incidence angle, and it indicates that the boundary layer could be ignored only in the case that the incidence angle is smaller than 8.5°.
By considering the flow angle as an independent variable, the measured velocity is plotted in Figure 13a. In the figure, the higher line in the left region represents the normalized position being farther back, and all of the lines are close to the straight line. Figure 13b shows the linearity of lines in Figure 13a in different flow positions and different incidence angle ranges (black line for the small incidence angle range, while the red line for the large incidence angle range). It can be seen that in the range of the small incidence angle, the surface velocity distribution ultimately presents a linear relationship with the change of the incidence angle. Hence, when there is no serious deterioration of the boundary layer near the leading edge, the slope of the lines in Figure 13a has a physical significance, called VI, as mentioned above.
For different intervals of incidence angles, such as [0, 4°], [0, 6°], [0, 7.5°], the value of VI measured by PIV in the leading edge region is calculated and plotted as the black dot in Figure 14. In this figure, the red line is the value calculated by CFD, and the blue line is the value calculated with Equation (16). In the experiment, VI obtained by CFD is very close to the experimental results, while VI calculated from Equation (16) is slightly different from the experimental results. However, considering that Equation (16) is an explicit analytical formula, its application is simple, so Equation (16) can be used as an effective and rapid estimation of VI near the leading edge of the blade.
Furthermore, both the theoretical analysis and experimental research were carried out in incompressible flow. For a subsonic compressible environment, compressibility correction can be used to further extend the calculation process in this paper [23,24,25].

5. Conclusions

In this paper, the influence of the off-designed incidence angle on the flow in the leading edge region is studied in the compressor cascade environment, and the concept of the velocity increment on incidence angle, VI, is raised to present the surface velocity variation with the off-designed incidence angle. With the analysis of potential theory, the explicit equation for VI is derived, and it is verified by numerical calculations and experimental measurement. The conclusions drawn are:
By employing the ratio of the value of VI between the plate cascade and the isolated plate, the infinity tendency of the current explicit solution at the leading edge point is eliminated. For the stagger angle of the cascade, the larger the stagger angle, the larger the value of VI near the leading edge region and the larger the difference between VI values of the pressure surface and suction surface in the blade body. Therefore, the equivalent pitch lengths based on 1/cos(β) and VI caused by diffuser deceleration in the cascade passage were employed to correct the effect of the stagger angle in the explicit equation in a simplified elliptic cascade.
For the real blade profile, VI near the leading edge is similar to that of an elliptic cascade with an axial ratio equal to the equivalent axial ratio,  c / 2 r L E , and the VI value caused by camber turning can be ignored near the leading edge region. Finally, the explicit equation of VI derived in this paper depends only on the geometrical parameters of the cascade and blade.
Through numerical calculations, the equation of VI has a good coincidence accuracy near the leading edge, trailing edge, and also in the blade body region of the pressure surface. By measuring the velocity field in the leading edge area with PIV, the linear rule of the blade surface velocity on the incidence angle in the leading edge region is verified, and the solving accuracies of the equation of VI are analyzed and verified.
According to the above research, the equation of VI is suitable for a wide range of incidence angles near the leading edge of a compressor blade cascade, which greatly improve the accuracy of an explicit analytical solution of blades under all working conditions.
For feasible further research, the equation of VI could be used in a design method with lesser iterations for blade geometry by considering flow characteristics under off-designed conditions, and it could also provide a basis for research on pollution deposit positions in aeroengines by providing the calculation method for a front-stagnation position in all flow conditions. It may also be helpful in the calculation of unsteady blade force of axial compressors.

Author Contributions

Conceptualization, X.X. and X.Y.; methodology, X.X. and X.Y.; software, X.X.; validation, B.L. and X.Y.; formal analysis, X.X. and X.Y.; investigation, X.X., B.L., X.Y. and G.A.; resources, X.Y. and G.A.; data curation, X.X. and X.Y.; writing—original draft preparation, X.X.; writing—review and editing, G.A.; visualization, X.X.; supervision, B.L.; project administration, B.L. and X.Y.; funding acquisition, B.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. 52276025), the National Science and Technology Major Project (2017-II-0001-0013, J2019-II-0004-0024, J2019-II-0003-0023), the Advanced Jet Propulsion Innovation Center/AEAC (funding number HKCX2022-01-008), and the Fundamental Research Funds for the Central Universities.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Abbreviations

SymbolDescription
α1Inlet flow angle
α2Outlet flow angle
ARAxial ratio of the ellipse
A R r Equivalent   axial   ratio   based   on   r L E
A R t Equivalent axial ratio based on the maximum thickness
β1Inlet metal angle
β2Outlet metal angle
βStagger angle
cChord length
H m a x Maximum thickness of the blade
iIncidence angle
r L E Curvature radius at leading edge point
sArc length on the blade from leading edge point
scamberArc length on the camber line from leading edge point
scNormalized camber projection length
S0Arc length on the blade from leading edge to trailing edge
T x Ratio of flow turning angle on incidence angle
tPitch length
t e q Equivalent pitch length
U Velocity (scalar)
V I Velocity increment on incidence angle
V I p l a t e , c a s c a d e VI value of the plate cascade
V I p l a t e , i s o l a t e d VI value of the isolated plate
V I e l l i p s e , i s o l a t e d VI value of the isolated ellipse
V I e l l i p s e , c a s c a d e VI value of the elliptic straight cascade
V I s t a g g e r , c o r r VI value of the elliptic staggered cascade
V I T E VI value at trailing edge point
V I d i f f u s i o n VI value caused by diffuser deceleration
V I β Comprehensive VI value

References

  1. Carter, A.D.S. Blade Profiles for Axial-Flow Fans, Pumps, Compressors, Etc. Proc. Inst. Mech. Eng. 1961, 175, 775–806. [Google Scholar] [CrossRef]
  2. Walraevens, R.E.; Cumpsty, N.A. Leading Edge Separation Bubbles on Turbomachine Blades. J. Turbomach. 1995, 117, 115–125. [Google Scholar] [CrossRef]
  3. Shi, L.; Ma, H. Effect of Incidence Angle on Entropy Generation in the Boundary Layers on the Blade Suction Surface in a Compressor Cascade. Entropy 2019, 21, 1049. [Google Scholar] [CrossRef]
  4. Lambert, A.; Yarusevych, S. Effect of angle of attack on vortex dynamics in laminar separation bubbles. Phys. Fluids 2019, 31, 064105. [Google Scholar] [CrossRef]
  5. Goodhand, M.N.; Miller, R.J. Compressor Leading Edge Spikes: A New Performance Criterion. J. Turbomach. 2011, 133, 021006. [Google Scholar] [CrossRef]
  6. Xiang, S.; Avital, E.; Rezaienia, M.A.; Paul, G.; Korakianitis, T. Computational Methods for Investigation of Surface Curvature Effects on Airfoil Boundary Layer Behavior. J. Algorithms Comput. Technol. 2017, 11, 68–82. [Google Scholar] [CrossRef]
  7. Liu, B.J.; Yuan, C.X.; Yu, X.J. Effects of leading-edge geometry on aerodynamic performance in controlled diffusion airfoil. J. Propuls. Technol. 2013, 34, 890–897. [Google Scholar]
  8. Song, Y.; Gu, C. Effects of Curvature Continuity of Compressor Blade Profiles on Their Performances. In Proceedings of the ASME Turbo Expo 2014: Turbine Technical Conference and Exposition, Düsseldorf, Germany, 16–20 June 2014. [Google Scholar]
  9. Li, Z.H.; Zheng, X.Q. Review of design optimization methods for turbomachinery aerodynamics. Prog. Aerosp. Sci. 2017, 93, 1–23. [Google Scholar] [CrossRef]
  10. Trigg, M.A.; Tubby, G.R.; Sheard, A.G. Automatic Genetic Optimization Approach to Two-Dimensional Blade Profile Design for Steam Turbines. J. Turbomach. 1999, 121, 11–17. [Google Scholar] [CrossRef]
  11. Li, Y.; Wei, K.; Yang, W.; Wang, Q. Improving wind turbine blade based on multi-objective particle swarm optimization. Renew. Energy 2020, 161, 525–542. [Google Scholar] [CrossRef]
  12. Tosin, S.; Friedrichs, J.; Sperling, J.; Kožulović, D. Aerodynamic Optimization of Turboprop Turbine Blades Using a Response Surface Methodology Based Algorithm. In Proceedings of the ASME 2014 12th International Conference on Nanochannels, Microchannels, and Minichannels, Chicago, IL, USA, 3–7 August 2014. [Google Scholar]
  13. Zeinalzadeh, A.; Pakatchian, M.R. Evaluation of novel-objective functions in the design optimization of a transonic rotor by using deep learning. Eng. Appl. Comput. Fluid Mech. 2021, 15, 561–583. [Google Scholar] [CrossRef]
  14. Qin, S.; Wang, S.Y.; Wang, L.Y.; Wang, C.; Sun, G.; Zhong, Y.J. Multi-Objective Optimization of Cascade Blade Profile Based on Reinforcement Learning. Appl. Sci. 2021, 11, 106. [Google Scholar] [CrossRef]
  15. Selig, M.S. Multipoint Inverse Design of an Infinite Cascade of Airfoils. AIAA J. 1994, 32, 774–782. [Google Scholar] [CrossRef]
  16. Selig, M.S.; Maughmer, M.D. Multipoint Inverse Airfoil Design Method based on Conformal Mapping. AIAA J. 1992, 30, 1162–1170. [Google Scholar] [CrossRef]
  17. Theodorsen, T.; Garrick, I.E. General Potential Theory of Arbitrary Wing Sections; No. 452; National Advisory Committee for Aeronautics: Hampton, VA, USA, 1933.
  18. Hall, W.S.; Thwaites, B. On the Calculation of Cascade Flows; C.P. No. 806; Her Majesty’s Stationery Office: London, UK, 1965. [Google Scholar]
  19. Baddoo, P.J.; Ayton, L.J. Potential Flow through a Cascade of Aerofoils: Direct and Inverse Problems. Proc. R. Soc. A Math. Phys. Eng. Sci. 2018, 474, 20180065. [Google Scholar] [CrossRef] [PubMed]
  20. Pollard, D.; Wordsworth, J. A Comparison of Two Methods for Predicting the Potential Flow around Arbitrary Airfoils in Cascade; C.P. No. 618; Her Majesty’s Stationery Office: London, UK, 1963. [Google Scholar]
  21. Thwaites, B. Incompressible Aerodynamics: An Account of the Theory and Observation of the Steady Flow of Incompressible Fluid Past Aerofoils, Wings, and Other Bodies; Clarendon Press: Oxford, UK, 1960. [Google Scholar]
  22. Liu, B.; Xu, X.; Yu, X.; An, G. Method for utilizing PIV to investigate high curvature and acceleration boundary layer flows around the compressor blade leading edge. Chin. J. Aeronaut. 2022, 35, 72–88. [Google Scholar] [CrossRef]
  23. Hess, R.V.; Gardner, C.S. Study by the Prandtl-Glauert Method of Compressibility Effects and Critical Mach Number for Ellipsoids of Various Aspect Ratios and Thickness Ratios; National Aeronautics and Space Administration: Washington, DC, USA, 1949.
  24. Karman, T.V. Compressibility Effects in Aerodynamics. J. Spacecr. Rocket. 2003, 40, 992–1011. [Google Scholar] [CrossRef]
  25. Laitone, E.V. New Compressibility Correction for Two-Dimensional Subsonic Flow. J. Aeronaut. Sci. 2015, 18, 350. [Google Scholar] [CrossRef]
Figure 1. Geometry parameters of the two-dimensional planer cascade.
Figure 1. Geometry parameters of the two-dimensional planer cascade.
Aerospace 10 00461 g001
Figure 2. Value of VI for plate cascade with different pitch length and axial position.
Figure 2. Value of VI for plate cascade with different pitch length and axial position.
Aerospace 10 00461 g002
Figure 3. Ratio of VI in the plate cascade to VI on an isolated plate.
Figure 3. Ratio of VI in the plate cascade to VI on an isolated plate.
Aerospace 10 00461 g003
Figure 4. True value and the value from Equation (8) of VI in an elliptic straight cascade. (a) t = 2, AR = 10, (b) t = 2.
Figure 4. True value and the value from Equation (8) of VI in an elliptic straight cascade. (a) t = 2, AR = 10, (b) t = 2.
Aerospace 10 00461 g004
Figure 5. VI distribution of an elliptic straight cascade with different stagger angles under t = 2, AR = 10. (a) VI calculated from vector velocity, (b) VI calculated from scalar velocity.
Figure 5. VI distribution of an elliptic straight cascade with different stagger angles under t = 2, AR = 10. (a) VI calculated from vector velocity, (b) VI calculated from scalar velocity.
Aerospace 10 00461 g005
Figure 6. Stagger effect on VI at leading edge, t = 2, AR = 10 [19].
Figure 6. Stagger effect on VI at leading edge, t = 2, AR = 10 [19].
Aerospace 10 00461 g006
Figure 7. Estimated VI of the elliptic straight cascades with stagger angle based on Equation (16). (a) t = 2, AR = 10, β = 30°, (b) t = 2, AR = 10.
Figure 7. Estimated VI of the elliptic straight cascades with stagger angle based on Equation (16). (a) t = 2, AR = 10, β = 30°, (b) t = 2, AR = 10.
Aerospace 10 00461 g007
Figure 8. Definitions in camber projection normalization method.
Figure 8. Definitions in camber projection normalization method.
Aerospace 10 00461 g008
Figure 9. Effect of ARt and ARr for a freedom thickness straight cascade. (a) Geometry, (b) value of VI when β = 30°, (c) value of VI with various β at leading edge.
Figure 9. Effect of ARt and ARr for a freedom thickness straight cascade. (a) Geometry, (b) value of VI when β = 30°, (c) value of VI with various β at leading edge.
Aerospace 10 00461 g009
Figure 10. Numerical verification with a typical compressor stator blade cascade. (a) Geometry, (b) value of VI, (c) differences between calculated values and CFD result.
Figure 10. Numerical verification with a typical compressor stator blade cascade. (a) Geometry, (b) value of VI, (c) differences between calculated values and CFD result.
Aerospace 10 00461 g010
Figure 11. Schematic of the camera field-of-view on the cascade. (a) Experimental setup, (b) sample image.
Figure 11. Schematic of the camera field-of-view on the cascade. (a) Experimental setup, (b) sample image.
Aerospace 10 00461 g011
Figure 12. Measured velocity at the main flow boundary.
Figure 12. Measured velocity at the main flow boundary.
Aerospace 10 00461 g012
Figure 13. Relationship between the surface velocity and the change of incidence angle. (a) Velocity distribution on the incidence angle at a fixed point, (b) linearity of lines.
Figure 13. Relationship between the surface velocity and the change of incidence angle. (a) Velocity distribution on the incidence angle at a fixed point, (b) linearity of lines.
Aerospace 10 00461 g013
Figure 14. Comparison of VI Values from CFD, Equation (16), and measured by PIV.
Figure 14. Comparison of VI Values from CFD, Equation (16), and measured by PIV.
Aerospace 10 00461 g014
Table 1. Basic parameters of the cascade facility.
Table 1. Basic parameters of the cascade facility.
ParametersValue
chord, c (mm)180
solidity1.6
inlet metal angle,  β 1  (°)47
outlet metal angle,  β 2  (°)6
designed inlet flow angle (°)43
inlet Mach number0.13
Reynolds number5.45 × 105
pitch, t (mm)112.5
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xu, X.; Liu, B.; Yu, X.; An, G. Velocity Increment on Incidence Angle near the Leading Edge of the Compressor Cascade. Aerospace 2023, 10, 461. https://doi.org/10.3390/aerospace10050461

AMA Style

Xu X, Liu B, Yu X, An G. Velocity Increment on Incidence Angle near the Leading Edge of the Compressor Cascade. Aerospace. 2023; 10(5):461. https://doi.org/10.3390/aerospace10050461

Chicago/Turabian Style

Xu, Xiaobin, Baojie Liu, Xianjun Yu, and Guangfeng An. 2023. "Velocity Increment on Incidence Angle near the Leading Edge of the Compressor Cascade" Aerospace 10, no. 5: 461. https://doi.org/10.3390/aerospace10050461

APA Style

Xu, X., Liu, B., Yu, X., & An, G. (2023). Velocity Increment on Incidence Angle near the Leading Edge of the Compressor Cascade. Aerospace, 10(5), 461. https://doi.org/10.3390/aerospace10050461

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop