FPCAS3D User's guide: A three dimensional full potential aeroelastic program, version 1

The FPCAS3D computer code has been developed for aeroelastic stability analysis of bladed disks such as those in fans, compressors, turbines, propellers, or propfans. The aerodynamic analysis used in this code is based on the unsteady three-dimensional full potential equation which is solved for a blade row. The structural analysis is based on a finite-element model for each blade. Detailed explanations of the aerodynamic analysis, the numerical algorithms, and the aeroelastic analysis are not given in this report. This guide can be used to assist in the preparation of the input data required by the FPCAS3D code. A complete description of the input data is provided in this report. In addition, six examples, including inputs and outputs, are provided

Flutter Version of Propulsion Aeroelasticity Code Completed

NASA s Advanced Subsonic Technology (AST) Program seeks to develop new technologies to increase the fuel efficiency of commercial aircraft engines, improve the safety of engine operation, and reduce emissions and engine noise. For new designs of ducted fans, compressors, and turbines to achieve these goals, a basic aeroelastic requirement is that there should be no flutter or high resonant blade stresses in the operating regime. For verifying the aeroelastic soundness of the design, an accurate prediction/analysis code is required. Such a three-dimensional viscous propulsion aeroelastic code, named TURBO-AE, is being developed at the NASA Lewis Research Center. The development and verification of the flutter version of the TURBO-AE code (version 4) has been completed. Validation of the code is partially complete

Propulsion Aeroelastic Analysis Developed for Flutter and Forced Response

The NASA Glenn Research Center at Lewis Field develops new technologies to increase the fuel efficiency of aircraft engines, improve the safety of engine operation, reduce emissions, and reduce engine noise. With the development of new designs for fans, compressors, and turbines to achieve these goals, the basic aeroelastic requirements are that there should be no flutter (self-excited vibrations) or high resonant blade stresses (due to forced response) in the operating regime. Therefore, an accurate prediction and analysis capability is required to verify the aeroelastic soundness of the designs. Such a three-dimensional viscous propulsion aeroelastic analysis capability has been developed at Glenn with support from the Advanced Subsonic Technology (AST) program. This newly developed aeroelastic analysis capability is based on TURBO, a threedimensional unsteady aerodynamic Reynolds-averaged Navier-Stokes turbomachinery code developed previously under a grant from Glenn. TURBO can model the viscous flow effects that play an important role in certain aeroelastic problems such as flutter with flow separation, flutter at high loading conditions near the stall line (stall flutter), flutter in the presence of shock and boundary-layer interaction, and forced response due to wakes and shock impingement. In aeroelastic analysis, the structural dynamics representation of the blades is based on normal modes. A finite-element analysis code is used to calculate these in-vacuum vibration modes and the associated natural frequencies. In an aeroelastic analysis using the TURBO code, flutter and forced response are modeled as being uncoupled. To calculate if a blade row will flutter, one prescribes the motion of the blade to be a harmonic vibration in a specified in-vacuum normal mode. An aeroelastic analysis preprocessor is used to generate the displacement field required for the analysis. The work done by aerodynamic forces on the vibrating blade during a cycle of vibration is calculated. If this work is positive, the blade is dynamically unstable, since it will extract energy from the flow, leading to an increase in the blade s oscillation amplitude. The forced-response excitations on a blade row are calculated by modeling the flow through two adjacent blade rows using the TURBO code. The blades are assumed to be rigid. As an option, a single blade row can be modeled with the upstream blade row influence represented by a time-varying disturbance (gust) at the inlet boundary. The unsteady forces on a blade row from such analyses are used in a structural analysis along with the blade structural dynamics characteristics and aerodynamic damping associated with blade vibration to calculate the resulting dynamic stresses on the blade

Flutter Stability of the Efficient Low Noise Fan Calculated

The TURBO-AE aeroelastic code has been used to verify the flutter stability of the Efficient Low Noise Fan (ELNF), which is also referred to as the trailing-edge blowing fan. The ELNF is a unique technology demonstrator being designed and fabricated at the NASA Glenn Research Center for testing in Glenn's 9-by-15-Foot Low-Speed Wind Tunnel. In the ELNF, air can be blown out of slots near the trailing edges of the fan blades to fill in the wakes downstream of the rotating blades. This filling of the wakes leads to a reduction of the rotor-stator interaction (tone) noise that results from the interaction of wakes with the downstream stators. The ELNF will demonstrate a 1.6-EPNdB1 reduction in tone noise through wake filling, without increasing the broadband noise. Furthermore, the reduced blade row interaction will decrease the possibility of forced response and enable closer spacing of blade rows, thus reducing engine length and weight. During the design of the ELNF, the rotor blades were checked for flutter stability using the detailed aeroelastic analysis capability of the three-dimensional Navier-Stokes TURBOAE code. The aeroelastic calculations were preceded by steady calculations in which the blades were not allowed to vibrate. For each rotational speed, as the back-pressure was increased, the mass flow rate decreased, and the operating point moved along the constant speed characteristic (speed-line) from choke to stall as shown on the fan map. The TURBO-AE aeroelastic analyses were performed separately for the first two vibration modes (bending and torsion) and covered the complete range of interblade phase angles or nodal diameters at which flutter can occur. The results indicated that the ELNF blades would not encounter flutter at takeoff conditions. The calculations were then repeated for a part-speed condition (70-percent rotational speed), and the results again showed no flutter in the operating region. On the fan map (shown), the predicted flutter point at part speed condition was located beyond the stall line, which means that the ELNF will not encounter flutter since it will never operate beyond the stall line. All the calculations done so far have been for the nonblowing case, and selected calculations will be repeated with air blowing from the trailing edge of the fan

Flutter Stability Verified for the Trailing Edge Blowing Fan

The TURBO-AE aeroelastic code has been used to verify the flutter stability of the trailing edge blowing (TEB) fan, which is a unique technology demonstrator being designed and fabricated at the NASA Glenn Research Center for testing in Glenn s 9- by 15-Foot Low-Speed Wind Tunnel. Air can be blown out of slots near the trailing edges of the TEB fan blades to fill in the wakes downstream of the rotating blades, which reduces the rotor-stator interaction (tone) noise caused by the interaction of wakes with the downstream stators. The TEB fan will demonstrate a 1.6-EPNdB reduction in tone noise through wake filling. Furthermore, the reduced blade-row interaction will decrease the possibility of forced-response vibrations and enable closer spacing of blade rows, thus reducing engine length and weight. The detailed aeroelastic analysis capability of the three-dimensional Navier-Stokes TURBO-AE code was used to check the TEB fan rotor blades for flutter stability. Flutter calculations were first performed with no TEB flow; then select calculations were repeated with TEB flow turned on

Harmonic Balance Computations of Fan Aeroelastic Stability

A harmonic balance (HB) aeroelastic analysis, which has been recently developed, was used to determine the aeroelastic stability (flutter) characteristics of an experimental fan. To assess the numerical accuracy of this HB aeroelastic analysis, a time-domain aeroelastic analysis was also used to determine the aeroelastic stability characteristics of the same fan. Both of these three-dimensional analysis codes model the unsteady flowfield due to blade vibrations using the Reynolds-averaged Navier-Stokes (RANS) equations. In the HB analysis, the unsteady flow equations are converted to a HB form and solved using a pseudo-time marching method. In the time-domain analysis, the unsteady flow equations are solved using an implicit time-marching approach. Steady and unsteady computations for two vibration modes were carried out at two rotational speeds: 100 percent (design) and 70 percent (part-speed). The steady and unsteady results obtained from the two analysis methods compare well, thus verifying the recently developed HB aeroelastic analysis. Based on the results, the experimental fan was found to have no aeroelastic instability (flutter) at the conditions examined in this study

Fan Flutter Computations Using the Harmonic Balance Method

An experimental forward-swept fan encountered flutter at part-speed conditions during wind tunnel testing. A new propulsion aeroelasticity code, based on a computational fluid dynamics (CFD) approach, was used to model the aeroelastic behavior of this fan. This threedimensional code models the unsteady flowfield due to blade vibrations using a harmonic balance method to solve the Navier-Stokes equations. This paper describes the flutter calculations and compares the results to experimental measurements and previous results from a time-accurate propulsion aeroelasticity code

Unsteady Flowfield in a High-Pressure Turbine Modeled by TURBO

Forced response, or resonant vibrations, in turbomachinery components can cause blades to crack or fail because of the large vibratory blade stresses and subsequent high-cycle fatigue. Forced-response vibrations occur when turbomachinery blades are subjected to periodic excitation at a frequency close to their natural frequency. Rotor blades in a turbine are constantly subjected to periodic excitations when they pass through the spatially nonuniform flowfield created by upstream vanes. Accurate numerical prediction of the unsteady aerodynamics phenomena that cause forced-response vibrations can lead to an improved understanding of the problem and offer potential approaches to reduce or eliminate specific forced-response problems. The objective of the current work was to validate an unsteady aerodynamics code (named TURBO) for the modeling of the unsteady blade row interactions that can cause forced response vibrations. The three-dimensional, unsteady, multi-blade-row, Reynolds-averaged Navier-Stokes turbomachinery code named TURBO was used to model a high-pressure turbine stage for which benchmark data were recently acquired under a NASA contract by researchers at the Ohio State University. The test article was an initial design for a high-pressure turbine stage that experienced forced-response vibrations which were eliminated by increasing the axial gap. The data, acquired in a short duration or shock tunnel test facility, included unsteady blade surface pressures and vibratory strains

Cascade flutter analysis with transient response aerodynamics

Two methods for calculating linear frequency domain aerodynamic coefficients from a time marching Full Potential cascade solver are developed and verified. In the first method, the Influence Coefficient, solutions to elemental problems are superposed to obtain the solutions for a cascade in which all blades are vibrating with a constant interblade phase angle. The elemental problem consists of a single blade in the cascade oscillating while the other blades remain stationary. In the second method, the Pulse Response, the response to the transient motion of a blade is used to calculate influence coefficients. This is done by calculating the Fourier Transforms of the blade motion and the response. Both methods are validated by comparison with the Harmonic Oscillation method and give accurate results. The aerodynamic coefficients obtained from these methods are used for frequency domain flutter calculations involving a typical section blade structural model. An eigenvalue problem is solved for each interblade phase angle mode and the eigenvalues are used to determine aeroelastic stability. Flutter calculations are performed for two examples over a range of subsonic Mach numbers

Linearized Aeroelastic Solver Applied to the Flutter Prediction of Real Configurations

A fast-running unsteady aerodynamics code, LINFLUX, was previously developed for predicting turbomachinery flutter. This linearized code, based on a frequency domain method, models the effects of steady blade loading through a nonlinear steady flow field. The LINFLUX code, which is 6 to 7 times faster than the corresponding nonlinear time domain code, is suitable for use in the initial design phase. Earlier, this code was verified through application to a research fan, and it was shown that the predictions of work per cycle and flutter compared well with those from a nonlinear time-marching aeroelastic code, TURBO-AE. Now, the LINFLUX code has been applied to real configurations: fans developed under the Energy Efficient Engine (E-cubed) Program and the Quiet Aircraft Technology (QAT) project. The LINFLUX code starts with a steady nonlinear aerodynamic flow field and solves the unsteady linearized Euler equations to calculate the unsteady aerodynamic forces on the turbomachinery blades. First, a steady aerodynamic solution is computed for given operating conditions using the nonlinear unsteady aerodynamic code TURBO-AE. A blade vibration analysis is done to determine the frequencies and mode shapes of the vibrating blades, and an interface code is used to convert the steady aerodynamic solution to a form required by LINFLUX. A preprocessor is used to interpolate the mode shapes from the structural dynamics mesh onto the computational fluid dynamics mesh. Then, LINFLUX is used to calculate the unsteady aerodynamic pressure distribution for a given vibration mode, frequency, and interblade phase angle. Finally, a post-processor uses the unsteady pressures to calculate the generalized aerodynamic forces, eigenvalues, an esponse amplitudes. The eigenvalues determine the flutter frequency and damping. Results of flutter calculations from the LINFLUX code are presented for (1) the E-cubed fan developed under the E-cubed program and (2) the Quiet High Speed Fan (QHSF) developed under the Quiet Aircraft Technology project. The results are compared with those obtained from the TURBO-AE code. A graph of the work done per vibration cycle for the first vibration mode of the E-cubed fan is shown. It can be seen that the LINFLUX results show a very good comparison with TURBO-AE results over the entire range of interblade phase angle. The work done per vibration cycle for the first vibration mode of the QHSF fan is shown. Once again, the LINFLUX results compare very well with the results from the TURBOAE code