1. Introduction
The demand for fuel-efficient, lightweight, and high-performance commercial aircraft engines has driven the manufacturers towards designing lighter blades with higher pressure ratios per stage making them prone to flutter. This aeroelastic instability is mostly seen in the first stages of the compressor and the last stages of the turbine. While the rotating stall of a compressor has been studied extensively from the aerodynamic perspective, the study of flutter/aeroelasticity provides another view for understanding the complex flow phenomenon of compressor stalls.
Flutter is a self-excited and self-sustained aeroelastic instability phenomenon where the unsteady pressure field of the flow caused due to the vibration of the structure adds energy to the structure, which is beyond the damping capacity of the system, causing an exponential increase in the oscillation amplitude. For example, in an ideal stable system, the aerodynamic forces acting on a structural element exposed to fluid flow are balanced out by the inertial and elastic forces, dissipating the energy imparted by the fluid to the structure. In turbomachinery, flutter occurs when energy imparted by the flow on the blade over a cycle due to the blade’s vibration at the natural frequency is sufficient to overcome the damping capability of the blade.
Flutter due to its inherent nonlinear nature and dependency on both the aerodynamic and structural characteristics of the blade made computational fluid dynamics (CFD) difficult to implement and computationally expensive. The advancements in integrated CFD methods and affordable computational power have given headway to significant progress in developing methods for flutter prediction in turbomachinery. Casoni and Benini [
1] in their study divided the existing methods for predicting flutter in turbomachinery into two broad categories, namely the classical method and the coupled method. In the classical method, the fluid and structural domain is computed separately, while in the coupled method the fluid and structural domains are considered as a single domain and there is a transfer of data between them, i.e., any changes in the fluid affect the structure and vice versa. Several studies [
2,
3,
4,
5,
6,
7,
8,
9,
10,
11,
12,
13] have been conducted using the uncoupled and coupled methods for accurately predicting flutter in turbomachinery. For example, the study by Zheng and Yang [
2] used a coupled fluid–structure method based on the time-domain solution of the fluid–structure where the aerodynamic and structural equations were marched simultaneously to predict flutter of a transonic fan. The method used was not only able to predict the stability of the blade but also able to compute the amplitude and frequency of the vibration. Another example of fully coupled flutter prediction is the study by Chen et al. [
3], where fully coupled numerical methodology was developed to calculate the flow-structure interaction and to predict transonic NACA 64A010 airfoil flutter. The study was not only able to capture the limit cycle oscillation (LCO) but also obtained the flutter boundary with the transonic dip. This paper focuses on the use of the unsteady uncoupled energy method to predict the flutter of a transonic compressor in regions of near-peak efficiency and stall. The energy method is a classical method developed by Carta [
4] based on calculating the algebraic sum of the work done by the aerodynamic forces on the blade oscillating with a pre-determined amplitude. Even though the fully coupled method is known for its accuracy and ability to capture coupled flutter characteristics, it can be computationally expensive.
The CFD energy method compared to the coupled methods of flutter prediction is efficient, accurate and computationally inexpensive as the fluid and structural domains are treated separately [
5,
6]. Chahine et al. [
5] in their study compared the coupled and decoupled fan flutter prediction methods for varying values of mass ratio and blade stiffness of a transonic three-dimensional rotor. The study concluded that the results of both methods were comparable for all the cases of isolated variation in mass ratio and stiffness from baseline and the decoupled/energy method was able to predict the stability for a maximum frequency shift value of 1.2%. The work done by Li and Sheng [
7] checked the extent, to which the energy method could be used to predict stall flutter in a transonic blade, and included a structural model with two degrees of freedom combined with an unsteady aerodynamic model. The stall flutter boundaries for two rotors were predicted with the energy method and were found to agree well with the experimental values. The work by Clark and Hall [
8] used a time-linearized Navier–Stokes computational method for predicting the onset of flutter in a viscous flow through a two-dimensional cascade to predict stall flutter by decomposing the flow into a nonlinear spatially varying mean flow with the addition of a small perturbation, and the comparison of the CFD results to experimental values of stall flutter were found to be in agreement.
Based on the above-mentioned literature survey, the main aim of this study is to investigate the application of the energy method for the prediction of stall flutter onset in a transonic compressor. The NASA Rotor 67 transonic compressor was chosen for this analysis due to the wide-range availability of its experimental and CFD data (in literature). The energy method is an uncoupled method of flutter study, which, in fact, is convenient due to its lower consumption of computational resources and time efficiency. In addition, the energy method has already been previously used for flutter prediction; however, the stall-related flutter onset prediction by using the energy method has not been fully addressed. Therefore, in this study, stall flutter onset prediction was carried out in a stepwise manner. In the first step, a modal analysis was carried out to evaluate the natural frequencies and mode shapes. The modal analysis results (of this study) were compared with the previous studies [
9] on the NASA Rotor 67. In the next step, CFD simulations were carried out for the evaluation of the compressor performance, which was close to the experimental results of Strazisar et al. [
14,
15]. In the final step, transient CFD simulations were performed to determine the blade aerodynamic damping for two cases: (i) near-peak efficiency, (ii) and stall condition. From the obtained values of aerodynamic damping, the system stability was evaluated. At the near-peak compressor efficiency, no flutter-related instability was detected, which was also proven in the previous studies [
6,
9]; however, for the compressor stall condition, system instability was detected by a negative aerodynamic damping value at a backward-traveling wave (negative nodal diameter). This stepwise compressor stall flutter analysis successfully demonstrated the applicability of the computational inexpensive energy method. Further validations of the stall flutter onset prediction along with the flow phenomenon causing such system instability are planned in the future.
This study is part of the series of current [
16,
17,
18,
19,
20] and future studies under the European Latvian Council of Science 4R-UAV project and the National Research Program Project: Sports, which is funded by the Ministry of Science and Education of Latvia with the aim to develop innovative and environmentally friendly aerospace and sports equipment at Riga Technical University, Latvia.
3. Results
The results of the FEM and CFD analyses are presented in this section. At first, the modal analysis was performed, and then a steady-state and transient analysis was performed to observe the onset of the flutter in the compressor blade at peak efficiency and stall condition.
3.1. Modal Analysis
The modal analysis was conducted to obtain the mode shapes of natural frequencies of the NASA Rotor 67 compressor blade. For this, the blade was prestressed by taking into consideration the centrifugal loading at compressor design speed, while the hub was set as a fixture (with 0 DOF). In accordance with the actual manufacturing practices, fully solid titanium blade material was chosen for the analysis. The mechanical properties of the chosen material are given in
Table 4, while the obtained first three natural frequencies and the associated mode shapes for the rotor blade are shown in
Figure 6,
Figure 7 and
Figure 8.
It can be observed in
Figure 6 for the first bending mode, the leading edge at the blade tip possesses the maximum displacement with an even displacement distribution from hub to tip. Similar to the first mode, the second bending mode also exhibited maximum displacement at the leading edge; however, slightly higher displacement was observed across the width of the blade compared to the first mode, as shown in
Figure 7. The third mode in
Figure 8 is the torsional mode with maximum displacement at both the leading and trailing edge. In this study, the obtained natural frequency values were compared (with a maximum deviation of 0.5%) with the frequency values obtained in the study by Ren et al. [
9], where the modal analysis was conducted on the NASA Rotor 67 as one of the steps to analyze couple-mode flutter in turbomachinery.
The natural frequency deviations are negligible. The natural frequency values obtained by the authors and the comparitive study are illustrated in
Table 5. The mode shapes from the modal analysis of the blade along with its natural frequencies, maximum displacement, and node coordinates (in the form of cartesian co-ordinates) were extracted from Ansys Mechanical using custom APDL code in CSV format to be imposed as mesh motion for the blade in Ansys CFX.
3.2. CFD Analysis; Compressor Performance
The steady-state CFD analysis was carried out and validated by comparing the compressor performance (obtained in this study) with the experimental data of Strazisar et al. [
14,
15] for the NASA Rotor 67. The CFD-based steady-state simulations were conducted for a fixed inlet pressure of 1 atm with variable outlet pressure using the k-ω SST turbulence model.
Figure 9 and
Figure 10 illustrate the compressor performance map where the variations of isentropic efficiency and compressor pressure ratios can be observed against the mass flow rate (normalized for choke conditions). The simulated results satisfy the experimental performance trends with minor deviations (less than 3%). These steady-state numerical validations set the base for initializing the transient cases, which help in achieving a faster convergence of the TBR simulations for predicting the onset of flutter, as discussed in
Section 2.3 (Computational Scheme).
3.3. Flutter Prediction
The TBR modelling with Fourier transformation were used to predict flutter onset in the Rotor 67 for the first three modes of natural vibrations and the nodal diameter in the range of −1 to +4 for compressor near-peak efficiency and stall condition. After the TBR simulation achieved transient periodic solution, the energy method was employed to calculate the aerodynamic damping in order to determine the stability of the system (blade) at each nodal diameter and natural blade vibration frequencies. The work on the blade in each vibration cycle was calculated using Equation (6), where
is the fluid pressure,
is the velocity to the predetermined displacements,
is the surface area of the blade,
is the surface normal unit vector, T is the period of one vibration cycle (
, and
is the time at the beginning of the vibration cycle [
23].
The nondimensionalized aerodynamic damping was obtained by dividing the work on the blade by the blade kinetic energy.
where
is the maximum displacement amplitude and
is the natural frequency of the mode (considered). In the energy method, positive aerodynamic damping designates a damped system for a particular nodal diameter; on the other hand, a negative value of aerodynamic damping indicates an undamped vibration, which is referred to as an unstable system with the possible onset of flutter.
For this study, flutter onset was analyzed for the two conditions, i.e., near-peak and stall conditions, as explained below:
Aerodynamic damping (near-peak efficiency):
The aerodynamic damping values for the near-peak efficiency region were calculated for the first three natural frequencies of the Rotor 67 (determined by modal analysis in
Section 3.1).
Figure 11,
Figure 12 and
Figure 13 present the aerodynamic damping of the first three natural frequencies for a range of nodal diameters (NDs). For the first mode (
Figure 11), the nodal diameter 1 is the least stable (lowest aerodynamic damping) and nodal diameter 4 the most stable (maximum aerodynamic damping). For the second mode (
Figure 12), a similar trend was observed i.e., highest damping at ND 4 (most stable) while lowest damping at ND 1 (least stable). In the case of the third mode, ND 0 is most stable while ND 1 is least stable, as shown in
Figure 13. To summarize, all three modes had a consistent trend of positive aerodynamic damping manifesting the overall stability of the rotor in the near-peak region, which means that there was no flutter predicted at the near-peak efficiency. This trend has also been observed in previous studies [
6,
9].
Aerodynamic Damping (Stall Condition):
The main aim of this study was to predict the onset of flutter in the compressor stall condition. Even though the energy method has been previously used in predicting the compressor blade flutter, its implementation for the stall-related flutter prediction is rather overlooked and rare, as literature mainly covers the flutter prediction in stable compressor-operating regions. Therefore, in this study, the aerodynamic damping values for the stall region were specifically evaluated.
Figure 14,
Figure 15 and
Figure 16 illustrate the aerodynamic damping of the first three natural frequencies for a range of NDs. The existence of a negative value of aerodynamic damping can be observed in the first mode at ND −1 (IBPA 16.36 degrees), which, in fact, depicts the system’s instability. The other two modes (second and third) show a stable tendency within all considered NDs.
4. Discussion and Conclusions
This paper explored a relatively “untouched” topic of flutter onset prediction in transonic compressor blade stall conditions using the energy method. As mentioned previously, the energy method is well explored and offers a convenient way of blade flutter analysis that is computationally inexpensive and can be conducted without advanced high-performance computing (HPC) systems. For careful analysis, each step of this study was either verified or compared with the known data. For example, the modal analysis, the compressor performance, and the flutter analysis (near-peak efficiency region) were thoroughly conducted and verified. Furthermore, in this study, the ability of the energy method to predict the stall flutter onset was also successfully explored. The energy method used in this study predicted the stall flutter onset at −1 nodal diameter. The authors strongly believe that the energy method is sufficient to predict the stall flutter onset with the least computational cost; however, further investigations at this nodal diameter (−1) are inevitable, and a thorough fluid flow phenomenon related to this should be explored further. While the energy method successfully predicted the stall flutter onset, further investigations for the expanded validity of the energy method are mandatory. In addition, the causes of the stall flutter onset and its related flow mechanism should also be investigated in the future. Such aeroelastic (aerodynamic and structural) analyses will also be helpful in understanding the classical question of the relationship between stall inception and flutter, which may include the computationally expensive coupled analysis.
The following conclusions can be drawn from the study:
The energy method successfully predicted the stall flutter onset in the transonic compressor.
At compressor near-peak efficiency, the energy method predicted the system stability; however, in stall condition, the energy method detected system instability at least for one nodal diameter.
Significantly, the prediction of stall flutter onset by the energy method was carried out using the FEM and transient CFD analyses without including any additional or uncommon methods. This computationally inexpensive and time-efficient method provides a wide range of applications for the compressor flutter prediction for future studies.
In the future, further investigations into the expanded validity of the energy method are planned. The relationship between the cause of the stall flutter inception and the related stall flow phenomenon will also be investigated in future studies. The knowledge attained from this study for aerodynamic simulations will be used for designing and simulating enhanced aerodynamic surfaces of sports equipment as a part of the VPP Sports Project.