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:
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]:
Furthermore, the error of Equation (1) in a general cascade is
, 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:
It shows the change in blade surface per unit incidence angle (in radians). Hence, VI of the plate cascade follows the equation:
and VI of the isolated plate follows the equation:
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,
, 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:
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:
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:
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
; 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:
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
. Furthermore, the method proposed by Baddoo and Ayton [
19] introduces the opposite correction. Therefore, the correction of Equation (4) changes to:
and the correction of Equation (8) changes to:
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:
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:
namely, the expression of VI at the trailing edge point. Considering the vector representation of blade surface velocity, the
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:
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:
Through the above analysis, the expression of blade surface VI
β of an elliptic cascade with a stagger angle is shown as:
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.
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, , 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.