Analysis of an unswept propfan blade with a semiempirical dynamic stall model

The time history response of a propfan wind tunnel model with dynamic stall is studied analytically. The response obtained from the analysis is compared with available experimental data. The governing equations of motion are formulated in terms of blade normal modes which are calculated using the COSMIC-NASTRAN computer code. The response analysis considered the blade plunging and pitching motions. The lift, drag and moment coefficients for angles of attack below the static stall angle are obtained from a quasi-steady theory. For angles above static stall angles, a semiempirical dynamic stall model based on a correction to angle of attack is used to obtain lift, drag and moment coefficients. Using these coefficients, the aerodynamic forces are calculated at a selected number of strips, and integrated to obtain the total generalized forces. The combined momentum-blade element theory is used to calculate the induced velocity. The semiempirical stall model predicted a limit cycle oscillation near the setting angle at which large vibratory stresses were observed in an experiment. The predicted mode and frequency of oscillation also agreed with those measured in the experiment near the setting angle

Numerical simulations of supersonic flow through oscillating cascade sections

A finite difference code was developed for modeling inviscid, unsteady supersonic flow by solution of the compressible Euler equations. The code uses a deforming grid technique to capture the motion of the airfoils and can model oscillating cascades with any arbitrary interblade phase angle. A flat plate cascade is analyzed, and results are compared with results from a small perturbation theory. The results show very good agreement for both the unsteady pressure distributions and the integrated force predictions. The reason for using the numerical Euler code over a small perturbation theory is the ability to model real airfoils that have thickness and camber. Sample predictions are presented for a cascade of loaded airfoils and show appreciable differences in the unsteady surface pressure distributions when compared with the flat plate results

Numerical analysis of supersonic flow through oscillating cascade sections by using a deforming grid

A finite difference code was developed for modeling inviscid, unsteady supersonic flow by solution of the compressible Euler equations. The code uses a deforming grid technique to capture the motion of the airfoils and can model oscillating cascades with any arbitrary interblade phase angle. A flat plate cascade is analyzed, and results are compared with results from a small-perturbation theory. The results show very good agreement for both the unsteady pressure distributions and the integrated force predictions. The reason for using the numerical Euler code over a small-perturbation theory is the ability to model real airfoils that have thickness and camber. Sample predictions are presented for a section of the rotor on a supersonic throughflow compressor designed at NASA Lewis Research Center. Preliminary results indicate that two-dimensional, flat plate analysis predicts conservative flutter boundaries

Aeroelastic Computations of a Compressor Stage Using the Harmonic Balance Method

The aeroelastic characteristics of a compressor stage were analyzed using a computational fluid dynamic (CFD) solver that uses the harmonic balance method to solve the governing equations. The three dimensional solver models the unsteady flow field due to blade vibration using the Reynolds-Averaged Navier-Stokes equations. The formulation enables the study of the effect of blade row interaction through the inclusion of coupling modes between blade rows. It also enables the study of nonlinear effects of high amplitude blade vibration by the inclusion of higher harmonics of the fundamental blade vibration frequency. In the present work, the solver is applied to study in detail the aeroelastic characteristics of a transonic compressor stage. Various parameters were included in the study: number of coupling modes, blade row axial spacing, and operating speeds. Only the first vibration mode is considered with amplitude of oscillation in the linear range. Both aeroelastic stability (flutter) of rotor blade and unsteady loading on the stator are calculated. The study showed that for the stage considered, the rotor aerodynamic damping is not influenced by the presence of the stator even when the axial spacing is reduced by nearly 25 percent. However, the study showed that blade row interaction effects become important for the unsteady loading on the stator when the axial spacing is reduced by the same amount

Stall flutter analysis of propfans

Three semi-empirical aerodynamic stall models are compared with respect to their lift and moment hysteresis loop prediction, limit cycle behavior, easy implementation, and feasibility in developing the parameters required for stall flutter prediction of advanced turbines. For the comparison of aeroelastic response prediction including stall, a typical section model and a plate structural model are considered. The response analysis includes both plunging and pitching motions of the blades. In model A, a correction of the angle of attack is applied when the angle of attack exceeds the static stall angle. In model B, a synthesis procedure is used for angles of attack above static stall angles, and the time history effects are accounted for through the Wagner function

User's Guide for ECAP2D: an Euler Unsteady Aerodynamic and Aeroelastic Analysis Program for Two Dimensional Oscillating Cascades, Version 1.0

This guide describes the input data required for using ECAP2D (Euler Cascade Aeroelastic Program-Two Dimensional). ECAP2D can be used for steady or unsteady aerodynamic and aeroelastic analysis of two dimensional cascades. Euler equations are used to obtain aerodynamic forces. The structural dynamic equations are written for a rigid typical section undergoing pitching (torsion) and plunging (bending) motion. The solution methods include harmonic oscillation method, influence coefficient method, pulse response method, and time integration method. For harmonic oscillation method, example inputs and outputs are provided for pitching motion and plunging motion. For the rest of the methods, input and output for pitching motion only are given

Euler flow predictions for an oscillating cascade using a high resolution wave-split scheme

A compressible flow code that can predict the nonlinear unsteady aerodynamics associated with transonic flows over oscillating cascades is developed and validated. The code solves the two dimensional, unsteady Euler equations using a time-marching, flux-difference splitting scheme. The unsteady pressures and forces can be determined for arbitrary input motions, although only harmonic pitching and plunging motions are addressed. The code solves the flow equations on a H-grid which is allowed to deform with the airfoil motion. Predictions are presented for both flat plate cascades and loaded airfoil cascades. Results are compared to flat plate theory and experimental data. Predictions are also presented for several oscillating cascades with strong normal shocks where the pitching amplitudes, cascade geometry and interblade phase angles are varied to investigate nonlinear behavior

The effects of rotational flow, viscosity, thickness, and shape on transonic flutter dip phenomena

The transonic flutter dip phenomena on thin airfoils, which are employed for propfan blades, is investigated using an integrated Euler/Navier-Stokes code and a two degrees of freedom typical section structural model. As a part of the code validation, the flutter characteristics of the NACA 64A010 airfoil are also investigated. In addition, the effects of artificial dissipation models, rotational flow, initial conditions, mean angle of attack, viscosity, airfoil thickness and shape on flutter are investigated. The results obtained with a Euler code for the NACA 64A010 airfoil are in reasonable agreement with published results obtained by using transonic small disturbance and Euler codes. The two artificial dissipation models, one based on the local pressure gradient scaled by a common factor and the other based on the local pressure gradient scaled by a spectral radius, predicted the same flutter speeds except in the recovery region for the case studied. The effects of rotational flow, initial conditions, mean angle of attack, and viscosity for the Reynold's number studied seem to be negligible or small on the minima of the flutter dip

Symbolic generation of elastic rotor blade equations using a FORTRAN processor and numerical study on dynamic inflow effects on the stability of helicopter rotors

The process of performing an automated stability analysis for an elastic-bladed helicopter rotor is discussed. A symbolic manipulation program, written in FORTRAN, is used to aid in the derivation of the governing equations of motion for the rotor. The blades undergo coupled bending and torsional deformations. Two-dimensional quasi-steady aerodynamics below stall are used. Although reversed flow effects are neglected, unsteady effects, modeled as dynamic inflow are included. Using a Lagrangian approach, the governing equations are derived in generalized coordinates using the symbolic program. The program generates the steady and perturbed equations and writes into subroutines to be called by numerical routines. The symbolic program can operate on both expressions and matrices. For the case of hovering flight, the blade and dynamic inflow equations are converted to equations in a multiblade coordinate system by rearranging the coefficients of the equations. For the case of forward flight, the multiblade equations are obtained through the symbolic program. The final multiblade equations are capable of accommodating any number of elastic blade modes. The computer implementation of this procedure consists of three stages: (1) the symbolic derivation of equations; (2) the coding of the equations into subroutines; and (3) the numerical study after identifying mass, damping, and stiffness coefficients. Damping results are presented in hover and in forward flight with and without dynamic inflow effects for various rotor blade models, including rigid blade lag-flap, elastic flap-lag, flap-lag-torsion, and quasi-static torsion. Results from dynamic inflow effects which are obtained from a lift deficiency function for a quasi-static inflow model in hover are also presented

Analyzing Aeroelasticity in Turbomachines

ASTROP2-LE is a computer program that predicts flutter and forced responses of blades, vanes, and other components of such turbomachines as fans, compressors, and turbines. ASTROP2-LE is based on the ASTROP2 program, developed previously for analysis of stability of turbomachinery components. In developing ASTROP2- LE, ASTROP2 was modified to include a capability for modeling forced responses. The program was also modified to add a capability for analysis of aeroelasticity with mistuning and unsteady aerodynamic solutions from another program, LINFLX2D, that solves the linearized Euler equations of unsteady two-dimensional flow. Using LINFLX2D to calculate unsteady aerodynamic loads, it is possible to analyze effects of transonic flow on flutter and forced response. ASTROP2-LE can be used to analyze subsonic, transonic, and supersonic aerodynamics and structural mistuning for rotors with blades of differing structural properties. It calculates the aerodynamic damping of a blade system operating in airflow so that stability can be assessed. The code also predicts the magnitudes and frequencies of the unsteady aerodynamic forces on the airfoils of a blade row from incoming wakes. This information can be used in high-cycle fatigue analysis to predict the fatigue lives of the blades