Document
NASA/TM—2014-218447
AIAA–2014–3733
Assessing Fan Flutter Stability in Presence of
Inlet Distortion Using One-Way and Two-Way
Coupled Methods
Gregory P. Herrick Glenn Research Center, Cleveland, Ohio
December 2014
NASA STI Program . . . in Pro ¿ le Since its founding, NASA has been dedicated to the • CONFERENCE PUBLICATION. Collected advancement of aeronautics and space science. The papers from scienti ¿ c and technical NASA Scienti ¿ c and Technical Information (STI) conferences, symposia, seminars, or other program plays a key part in helping NASA maintain meetings sponsored or cosponsored by NASA.
this important role.
• SPECIAL PUBLICATION. Scienti ¿ c, The NASA STI Program operates under the auspices technical, or historical information from of the Agency Chief Information Of ¿ cer. It collects, NASA programs, projects, and missions, often organizes, provides for archiving, and disseminates concerned with subjects having substantial NASA’s STI. The NASA STI program provides access public interest.
to the NASA Aeronautics and Space Database and its public interface, the NASA Technical Reports • TECHNICAL TRANSLATION. English- Server, thus providing one of the largest collections language translations of foreign scienti ¿ c and of aeronautical and space science STI in the world. technical material pertinent to NASA’s mission.
Results are published in both non-NASA channels and by NASA in the NASA STI Report Series, which Specialized services also include organizing includes the following report types: and publishing research results, distributing specialized research announcements and feeds, • TECHNICAL PUBLICATION. Reports of providing information desk and personal search completed research or a major signi ¿ cant phase support, and enabling data exchange services.
of research that present the results of NASA programs and include extensive data or theoretical For more information about the NASA STI analysis. Includes compilations of signi ¿ cant program, see the following: scienti ¿ c and technical data and information deemed to be of continuing reference value. • Access the NASA STI program home page at NASA counterpart of peer-reviewed formal http://www.sti.nasa.gov professional papers but has less stringent limitations on manuscript length and extent of • E-mail your question to help@sti.nasa.gov graphic presentations.
• Phone the NASA STI Information Desk at • TECHNICAL MEMORANDUM. Scienti ¿ c 757-864-9658 and technical ¿ ndings that are preliminary or of specialized interest, e.g., quick release • Write to: reports, working papers, and bibliographies that NASA STI Information Desk contain minimal annotation. Does not contain Mail Stop 148 extensive analysis. NASA Langley Research Center Hampton, VA 23681-2199 • CONTRACTOR REPORT. Scienti ¿ c and technical ¿ ndings by NASA-sponsored contractors and grantees.
NASA/TM—2014-218447
AIAA–2014–3733
Assessing Fan Flutter Stability in Presence of
Inlet Distortion Using One-Way and Two-Way
Coupled Methods
Gregory P. Herrick Glenn Research Center, Cleveland, Ohio Prepared for the 50th Joint Propulsion Conference cosponsored by AIAA, ASME, SAE, and ASEE Cleveland, Ohio, July 28–30, 2014 National Aeronautics and Space Administration Glenn Research Center Cleveland, Ohio 44135
December 2014
Level of Review : This material has been technically reviewed by technical management.
Available from NASA STI Information Desk National Technical Information Service Mail Stop 148 5301 Shawnee Road NASA Langley Research Center Alexandria, VA 22312 Hampton, VA 23681-2199 Available electronically at http://www.sti.nasa.gov
Assessing Fan Flutter Stability in Presence of Inlet Distortion Using
One-Way and Two-Way Coupled Methods
Gregory P. Herrick National Aeronautics and Space Administration Glenn Research Center Cleveland, Ohio 44135
Abstract
Concerns regarding noise, propulsive efficiency, and fuel burn are inspiring aircraft designs wherein the propulsive turbomachines are partially (or fully) embedded within the airframe; such designs present serious concerns with regard to aero- dynamic and aeromechanic performance of the compression system in response to inlet distortion. Previously, a preliminary design of a forward-swept high-speed fan exhibited flutter concerns in clean-inlet flows, and the present author then studied this fan further in the presence of off-design distorted in-flows. Continuing this re- search, a three-dimensional, unsteady, Navier-Stokes computational fluid dynamics code is again applied to analyze and corroborate fan performance with clean inlet flow and now with a simplified, sinusoidal distortion of total pressure at the aero- dynamic interface plane. This code, already validated in its application to assess aerodynamic damping of vibrating blades at various flow conditions using a one-way coupled energy-exchange approach, is modified to include a two-way coupled time- marching aeroelastic simulation capability. The two coupling methods are compared in their evaluation of flutter stability in the presence of distorted in-flows.
1 Introduction
Since the advent of turbojet propulsion, “clean” intake flow for the fan (or propul- sor) of the gas turbine engine has been the desired, elusive, goal of airframe, inlet, and propulsion system designers alike in the quest for optimal performance of the propulsion system. As we look to the future, design goals now focus on reduced emissions, reduced fuel consumption, better propulsive efficiency (of the entire air- craft system), and reduced ambient noise. Among the prominent design concepts under development for application to future aircraft is the “Blended Wing Body” (alternatively known as “hybrid wing/body”). Most of these prospective designs feature propulsion systems wherein the propulsive fan, whether coupled directly to a Brayton cycle turbine engine or part of a turbo-electric distributed propulsion system, is partially or fully embedded within the body.
Boundary layer ingestion (BLI) propulsion has the potential for significant re- duction in aircraft fuel burn [1]; previous system studies [2–6] have shown that 5-10% reduction in fuel burn is possible. Recent work [7] has confirmed this potential ben- efit through a system study focused on the propulsion system and its integration into the BWB aircraft. The system study indicated that low-loss inlets and high- N ASA/TM—2014-218447 1 performance, distortion-tolerant turbomachinery are key technologies required to achieve a 3-5% BLI fuel burn benefit for future aircraft relative to a baseline high- performance, pylon-mounted, propulsion system. A robust aeromechanical design with respect to dynamic stresses and flutter stability in the presence of distortion is critical to the successful implementation and application of BLI propulsion systems.
In previous work from Herrick [8], a parallel computational fluid dynamics (CFD) code was modified and customized to begin studying the issues relevant to aerome- chanical response of turbomachinery to distorted inlet flows. In that effort, the aeroelastic analysis employed in the CFD code was “one-way” coupled: The blades’ motion was prescribed and enforced, perturbing the flow field, while the resultant changes in blade surface pressures were not fed back to alter the blades’ motions. In this effort, a “two-way” coupled method of aeroelastic analysis is implemented into the code and applied in the analysis of the fan subject to distorted inlet flows: The motion of the blades perturbs the flow field, and the flow field’s unsteady pressures in turn influence the motion of the blades. Many two-way coupled methods for tur- bomachinery aeroelasticity have been presented in the literature [9–11]; Br´ eard [12] et al. applied their integrated model to the study of forced response due to inlet distortion.
2 Background
The CFD code modified and applied in this study is TURBO, a physics-based sim- ulation tool for multistage turbomachinery. The solver computes the fluid conser- vation laws without ad hoc modeling of any flow phenomena other than models required for turbulence. This code solves the unsteady Reynolds-averaged Navier- Stokes equations and a decoupled k- turbulence model developed by Zhu and Shih [13]. The code is implemented in a portable, scalable form for distributed- memory parallel computers using MPI message passing. The parallel implementa- tion employs domain decomposition and supports general multi-block grids with arbitrary grid-block connectivity. The solution algorithm is a Newton iterative implicit time-accurate scheme with characteristics-based finite-volume spatial dis- cretization. The Newton subiterations are solved using a concurrent block-Jacobi symmetric Gauss-Seidel (BJ-SGS) relaxation scheme. Because all of the fundamen- tal fluid mechanics are computed, the code is capable of capturing the nonlinear characteristics of the flow fields of interest. With the actual modeling of grid move- ment of blade rows in relative motion, this code is capable of computing the unsteady interactions between blade rows. Details of the flow solver are given by Chen and Whitfield [14]. The approach to parallelization for large-scale, complex problems is discussed by Chen and Briley [15]. TURBO has been previously modified for application to study flutter and forced vibration [16–18]. TURBO has been further validated in several previous research efforts regarding aeromechanics [19–26].
The fan of interest is comprised of twenty-two blades. The fan was designed for purposes of noise reduction in a typical podded-engine mounting. Physically, the average radial tip gap as gridded for this study is about 0.28% of blade radial span. The fan was never intended to be part of an embedded propulsion system with distorted inlet flow. Previous aeromechanics analyses [25] of this fan exhibited good corroboration between computational and experimental performance in clean-inlet N ASA/TM—2014-218447 2 flow conditions.
In his initial study of this fan’s aeroelastic behavior in distorted inlet flows, Her- rick [8] modified the code’s total-condition preserved inlet boundary condition to allow for circumferential nonuniformity at the inlet plane. Application of a throttle exit boundary condition, allowing the exit flow to attain equilibrium without arti- ficial constraint, showed that inlet distortions do not fully attenuate upon passage through a single rotor. In the previous work, examination of aeroelastic behavior – using a prescribed motion, one-way coupled approach – in clean inlet conditions gen- erally yielded clear, consistent solutions. The distorted inlet aeroelastic simulations (using the same one-way coupled approach, but with an inlet total pressure distribu- tion approximating that of Inlet A [27]) yielded dramatic variation in aerodynamic damping computations about the annulus, with some annular sectors exhibiting great aeroelastic stability and other annular sectors indicating significant aeroelas- tic instability. These intriguing observations prompt further investigation into the methods used to assess aeroelastic behavior in a distorted flow environment.
This study shall further investigate several issues relevant to the computational assessment of flutter stability for turbomachinery rotors in distorted flow fields: “one-way coupled” versus “two-way coupled” motion of the blades in the unsteady flow field, post-processing/interpretation of computed output (modal force, modal displacement), and initial conditions. The fan is gridded with a single H-block topology for the fan passage, with 214 grid points streamwise, 76 grid points ra- dially, and 58 grid points pitchwise; the blade itself has 81 grid points stream- wise and 70 points in the spanwise direction. No tip grid is present: Instead, a lossless multi-point constraint is applied between the suction and pressure sides in the region between the blade’s radial extent and the casing. Computationally, the fan is preceded by inlet duct of length-ratio 1.387 with respect to blade-tip ax- ial chord and followed by exit duct of length-ratio 2.365 with respect to blade-tip axial chord. Again, the total-condition preserved boundary condition, capable of enforcing circumferentially-varying flow conditions, is applied at the inlet/AIP. The throttle boundary condition is again applied at the exit plane.
To facilitate the goals of this research effort, some fundamental assumptions and simplifications are applied in the study regarding operating speed, inlet distor- tion pattern, blade mode shape, and blade natural frequency. The fan is operated at 10800 rpm, a part-speed condition. From this part-speed condition the fan’s aeroelastic behavior is assessed near op-line and near stall. The mode shape and frequency of interest are artificial; its frequency is dictated to be 220 Hz such that a blade completes one oscillation cycle in the course of traversing precisely 18 blade pitches at 10800 rpm. The CFD mesh of the blade, with deflected geometry, is shown in Fig. 1. The artificial mode shape and frequency, in concert with the part- speed run condition, are chosen with intent to encounter stable and unstable flutter conditions.
The total pressure distortion at the inlet plane is a simple, single-frequency, once-per-revoluton spatial distribution. The simulations are performed rotor-alone, and only one mode is studied, thus the only fundamental frequencies at play are the blades’ vibrational frequency and the rotor’s rotational frequency. The single- frequency, once-per-revolution spatial distortion of total pressure is created by apply- ing a 4% (mean-to-peak) sinusoidal variation upon the (circumferentially-constant, radially-varying) clean flow total pressure profile in a single wave about the in- N ASA/TM—2014-218447 3 let annulus. No circumferential nonuniformity is applied to total temperature or flow angles. Simulations reveal that the sinusoidal distortion yields a 0.12% higher (energy-averaged) integrated total pressure at the inlet plane than the clean, ax- isymmetric pressure distribution. The fan characteristic at part-speed is shown in Fig. 2. The total pressure distributions for clean and distorted in-flow are shown in Fig. 3. CFD simulations are executed with 15 time steps per blade pitch; for clean inlet simulations, 6 Newton subiterations are applied at each time step, while 9 Newton subiterations per time step are used for the sinusoidally-distorted inlet runs.
1.35 Clean Distorted 1.3 Total Pressure Ratio [Pte/Pti] 1.25 31 32 33 34 35 Corrected Mass Flow [kg/s] Figure 2: Speedlines: Clean and distorted Figure 1: Mesh and deflected shape inlets, part-speed.
applied to blade in this research.
(a) Clean Inlet (b) Sinusoidal Inlet Distortion Figure 3: Inlet total pressure distributions for the clean and distorted inlets of this ◦ study. Common colorbar, red is high. Azimuthal orientation: 0 at 12 o’clock, increasing clockwise. Fan rotates counterclockwise.
N ASA/TM—2014-218447 4
3 One-way Coupled Method for Aeroelastic Analysis
3.1 Description and Implementation Flutter stability of the fan is first evaluated using the TURBO aeroelastic analysis code’s one-way coupled method. The energy method [28], or work-per-cycle ap- proach, is used to calculate aerodynamic damping (and hence flutter stability) as described in Ref. 26; Vasanthakumar [29] used a similar approach in flutter analysis of a transonic fan. For a clean inflow prescribed at the inlet boundary, the flow field through the blade row is first calculated with no prescribed blade oscillation.
Starting with this converged steady flow field, blade vibrations are prescribed in a selected mode, frequency ( f ), and nodal diameter pattern (ND) or inter-blade phase angle ( σ ). Note that the number of different possible nodal diameter patterns is equal to the number of blades on the rotor of interest. After the transients in the flow field decay and a periodic flow field is obtained, the work done on the vibrating blade is calculated as follows: ∮ ∫ ∂ X W = − p d A · dt (1) ∂t surf ace where p is the aerodynamic pressure, A is the blade surface area vector, X is the displacement vector on the blade surface, and t denotes time. The aerodynamic damping ratio ( ζ ) can be approximately related to the work-per-cycle ( W ) and the average kinetic energy ( K ) of the blade over one cycle of vibration through the E following expression: W ζ ≈ − (2) 8 K E The schematic for this method is shown in Fig. 4. If aerodynamic damping is negative (i.e., aerodynamic work is positive), flutter may occur. Note that structural damping (material and mechanical) has not been considered.
Calculate new Calculate unsteady aerodynamic flow field on work acting new mesh on blade Calculate new Generate new mesh blade position in accordance and velocity with blade position from prescribed harmonic motion Figure 4: One-way coupling time marching scheme.
N ASA/TM—2014-218447 5 3.2 Application to Distorted In-Flows In previous research, Herrick [8] and Bakhle [26] et al. have utilized the one-way coupled method to assess stability of fans in the presence of distorted in-flows. In each of these endeavors, the temporally- and spatially-averaged damping of all blades – a single value attributed to an idealized “average” blade – was used to gauge the flutter stability of the entire rotor in the presence of the circumferentially-varying flow field.
In the presence of a clean-inlet flow field, this single-value average represents very well the consistent response of each and every blade as the blade traverses the flow field. In Fig. 5(a) the convergence time history of the aerodynamic work on all blades is plotted along with the single-value average; in Fig. 5(b), the annulus is “unrolled”, and the aerodynamic work distribution for each of the twenty-two blades is plotted versus each blade’s respective annular position at the end of each oscillation cycle. In both plots of Fig. 5, the eleven-nodal diameter condition demonstrates the common occurrence whereby the blade may oscillate through a few cycles prior to reaching a converged state. As to be expected for a flow condition with axisymmetric conditions at the inlet plane, Fig. 5(b) shows great consistency about the annulus in a clean flow field. Analyzing blade-to-blade aerodynamic work data in the “unrolled” presentation of Fig. 5(b) proves to be particularly useful when studying a flow condition with inlet distortion.
5e-06 5e-06 0 0 -5e-06 -5e-06 -1e-05 -1e-05 0 ND, Cycle 1 Normalized Aerodynamic Work [-] Normalized Aerodynamic Work [-] 0 ND, Blades 1-22 0 ND, Cycles 2-5 -1.5e-05 -1.5e-05 11 ND, Blades 1-22 11 ND, Cycle 1 0 ND, Average Blade 11 ND, Cycles 2-5 11 ND, Average Blade -2e-05 -2e-05 1 2 3 4 5 0 90 180 270 360 Oscillation Cycle Count Azimuthal Position of Root LE at end of oscillation cycle [deg] (a) Convergence of aerodynamic work for all (b) Convergence of aerodynamic work for all blades versus oscillation cycle count. blades versus azimuthal position.
Figure 5: Time-convergence of aerodynamic work at part-speed, near op-line, clean inlet condition.
In the presence of circumferentially-varying distorted in-flow, the aerodynamic work done on/by the blade – and hence its damping – varies considerably as the blade traverses the flow field, see Fig. 6. Note that the ordinates are of identical scale among all plots in Figures 5 and 6, providing perspective on the constancy of the aerodynamic work about the annulus for the clean inlet condition and the varia- tion of the aerodynamic work about the annulus for the sinusoidally-distorted inlet condition. Observe further the close similarity between the converged “average” blade’s aerodynamic work for both the clean inlet and the sinusoidally-distorted in- let. Recall that the the total pressure distribution of the distorted inlet applied here was specified to be a sinusoidal variation (of magnitude 4%) about the clean inlet’s N ASA/TM—2014-218447 6 radially-varying profile of total pressure. Close examination of the “average” blade’s aerodynamic work quantities reveals the sinusoidally-distorted inlet condition to be marginally/negligibly less stable – but still stable – than the clean inlet condition near op-line.
In instances of inlet distortion wherein the aerodynamic work computations of the one-way coupling method vary from cycle to cycle and from blade to blade about the annulus – and all quantities indicate stability – it is conservative and reasonable to declare the system to be stable; in Fig. 6, the eleven-nodal pattern is clearly stable, and in the work of Bakhle [26] et al. the system was clearly stable under these auspices. However, observe in Fig. 6(b) that a sector spanning approximately ◦ ◦ 150 of the 360 annulus exhibits unstable (positive) aerodynamic work for the zero-nodal diameter pattern. While the “average” blade indicates stability for zero nodal diameters, one cannot conservatively assert that the system is stable with the one-way coupling method consistently indicating instability over a defined sector of the annulus. Similar behavior, not shown here, is observed for for one-, two-, and three-nodal diameter forward traveling waves with inlet distortion near op-line.
Complemenary one-way coupling aeroelastic analyses are performed for the near- stall condition. Again, the “average” blade aerodynamic work computations (and hence aerodynamic dampings) are very similar between clean inlet and sinusoidally- distorted inlet conditions, with the sinusoidally-distorted inlet conditions usually yielding marginally/negligibly less stable dampings than the clean inlet for the part- speed, near-stall condition. For both the clean inlet and sinusodially- distorted inlet, the near-stall “average” blade dampings are less stable than the op-line “average blade” dampings. Notably, the 1 ND FTW (one nodal diameter forward traveing wave) pattern for the part-speed, near-stall condition yields small negative values of damping (i.e., flutter-unstable) for the “average” blade, whereby the sinsusoidally- distorted inlet is marginally/negligibly more unstable. As shown in Figure 7, for both clean inlet and distorted inlet for both one- and two-nodal diameter forward traveling waves, all blades alternate between stable (negative) aerodynamic work and unstable (positive) aerodynamic work as they circumnavigate the annulus. Each black line in Fig. 7(a) represents each individual blade as it traverses the annulus of the clean-inlet flow, while each red line plots the path of individual blades rotating through the distorted-inlet flow, for the one-nodal diameter forward traveling wave.
Figure 7(b) shows similar behavior for the two-nodal diameter forward traveling wave.
Nodal diameter sweeps of “average” blade damping from zero-nodal diameter to three-nodal diameter forward-traveling wave for the part-speed condition – for both near-stall and near op-line, and for both clean inlet and distorted inlet – are plotted in Fig. 8. The eleven-nodal diameter pattern is strongly stable in all studies here; these data points are omitted from Fig. 8 for clarity of the data of nodal diameters zero through three.
4 Two-way Coupled Method for Aeroelastic Analysis
4.1 Description and Implementation The two-way coupled method is implemented via a modal superposition framework similar to Gnesin and Rz¸ adkowski [10], Carstens [9] et al., and Br´ eard [12] et al. In N ASA/TM—2014-218447 7 5e-06 5e-06 0 0 -5e-06 -5e-06 -1e-05 -1e-05 Normalized Aerodynamic Work [-] Normalized Aerodynamic Work [-] 0 ND, Cycle 1 -1.5e-05 -1.5e-05 0 ND, Blades 1-22 0 ND, Cycles 2-5 11 ND, Blades 1-22 11 ND, Cycle 1 0 ND, Average Blade 11 ND, Cycles 2-5 11 ND, Average Blade -2e-05 -2e-05 1 2 3 4 5 0 90 180 270 360 Oscillation Cycle Count Azimuthal Position of Root LE at end of oscillation cycle [deg] (a) Convergence of aerodynamic work for all (b) Convergence of aerodynamic work for all blades versus oscillation cycle count. blades versus azimuthal position.
Figure 6: Time-convergence of aerodynamic work at part-speed, near op-line, sinu- soidal distortion inlet condition.
6e-06 6e-06 4e-06 4e-06 2e-06 2e-06 0 0 -2e-06 -2e-06 Normalized Aerodynamic Work [-] Normalized Aerodynamic Work [-] -4e-06 -4e-06 -6e-06 -6e-06 0 90 180 270 360 0 90 180 270 360 Azimuthal Position of Root LE at end of oscillation cycle [deg] Azimuthal Position of Root LE at end of oscillation cycle [deg] (a) 1 ND FTW (b) 2 ND FTW Figure 7: Aerodynamic work of all blades versus azimuthal position: Near stall, both clean (black) and distorted (red) inlets.
Near-Stall Clean Near-Stall, Distorted Op-line, Clean 10 Op-Line, Distorted Aerodynamic Damping [%] -5 0 1 2 3 Nodal Diameter Figure 8: Aerodynamic damping versus nodal diameter using average-blade aero- dynamic work from one-way coupled method.
N ASA/TM—2014-218447 8 the two-way coupled method, the blade motion and the aerodynamic force acting on the blade each respond to, and act on, one another in sequential fashion in accordance with the equation of motion. Mode shapes are unit-mass normalized and hence modal stiffnesses equate to the squares of the respective frequencies. Modal forces are computed by taking the inner product of the blade surface pressure forces with the blade mode shape.
From the original equation of motion, ¨ ˙ [ M ] { X } + [ C ] { X } + [ K ] { X } = { F } − { f } (3) s T Pre-multiply by the mode shape [Φ] , establish the generalized displacement coor- dinate η , and normalize to unit modal mass: T [Φ] [ M ][Φ] = 1 (4) { X } = [Φ] { η } (5) Applying the conservative simplification of neglecting material and structural damp- ing (i.e., assuming C = 0). yields independent second order ordinary differential equations for each mode i : 2 T ¨ η + ω η = [Φ ] { F − f } (6) i i i i s,i i Equation 6 represents the equation of motion for the blade in terms of its generalized displacement coordinate, η . This equation is time-integrated using a Newmark- β method [30]. This equation applies the modal forces, F , acting on the blades from the existing flow field (with the blades in their existing orientations) to compute the resulting deflections of the blades; these deflections are then applied to the blades to update the grid, and the aerodynamics and modal forces are recomputed and updated; this sequence is shown in Fig. 9.
For each computation and utilization of modal force in Eq. 6, a datum static modal force f – representing the equilibrium loading on the undeformed hot blades s – must be debited such that the blade vibrates about a mean equilibrium position corresponding to the undeformed hot blade. While for axisymmetric clean inlet flow Calculate new Calculate unsteady aerodynamic flow field on forces acting new mesh on blade Calculate new Generate new mesh blade position in accordance and velocity with blade position from aeroelastic equation of motion Figure 9: Two-way coupling time marching scheme.
N ASA/TM—2014-218447 9 the pressure forces acting on the blades are nominally constant as the blade circum- navigates the annulus, these pressure forces acting on the blade can vary significantly about the annulus in the presence of distorted in-flows, hence the computation and utilization of the static modal force demand careful consideration. This variation in unsteady pressures about the annulus is plotted in Fig. 10. Observe the con- siderable variation in blade loading – nearly 25% (mean-to-peak) near op-line – as the blade traverses the annulus in the distorted flow (with an inlet total pressure variation of 4% mean-to-peak). As expected, the blade is more highly loaded near stall than near op-line. While the mean modal force for the distorted flow near op-line is almost equivalent to the nearly-constant modal force for the clean flow at the near op-line condition, there is greater difference between the mean modal force of distorted flow and the nearly-constant modal force of clean flow at the near stall condition. Also, while the variation in modal force within distorted flow near op-line is nearly sinusoidal, the variation in modal force within distorted flow near stall is more disorderly.
In the present implementation into TURBO, the two-way coupled method allows (or requires) user-specification of modal mass, stiffness or frequency of the mode(s), mechanical damping, initial blade displacements and velocities (or prescription of an initial nodal diameter pattern [or interblade phase angle] with displacement am- plitude), and β and γ parameters for the Newmark- β time integration. The user must also specify parameters of the “static modal force” for the simulation: retriev- ing modal force data on rigid blades over a period of one step, one blade pitch, or one revolution; and then sampling that static modal force data as discrete values, blade pitch averages, or full revolution averages. The fluid-structure coupling stencil applied in this study may be classified as conventional serial staggered; improved serial staggered [31] and predictor-corrector stencils are in validation and shall be applied in future work.
4.2 Validation While this two-way coupled method for aeroelastic analysis has been developed and implemented to facilitate greater insight into aeroelastic behavior of turbomachin- ery rotors in the presence of highly unsteady flows, it is instructive to first verify the performance of the algorithm in simpler, more fundamental aeroelastic environ- ments. Given that the one-way coupled analysis procedure deflects the blades in a prescribed sinusoidal motion (irrespective of the aerodynamic modal forces acting on the blades), the two-way coupled analysis procedure should produce identical results when the computed modal forces are set to zero in the solution of Eq. 6. Figure 11 shows the requisite replication of the prescribed motion displacement at the blade tip, suction side, midchord (part-speed, clean inlet, near-stall shown) when using the two-way coupled motion procedure with modal forces neglected; this corroboration verfies that the Newmark- β procedure adequately resolves the displacement-time evolution of the undamped elastic blade in a clean flow, near-stall environment.
Further validation of the method is evident from examination of its performance with its typical modal force coupling enforced. Figure 12 shows the blade’s decreas- ing magnitude of oscillation in the 0 ND pattern near op-line in clean flow. The prescribed motion analysis conveys dampings of 1.64% for zero-nodal diameters (see Fig. 8), while logarithmic-decrement analysis of the tightly-coupled method’s blade N ASA/TM—2014-218447 10 -0.518 0.004 One-way coupling Two-way Coupling, Zero Modal Force -0.519 -0.52 0.0035 -0.521 0.003 -0.522 -0.523 Normalized Modal Force [-] Near Op-Line, Clean, Discrete 0.0025 Near Stall, Clean, Discrete -0.524 Normalized Axial Position of Tip Midchord [-] Near Op-Line, Distorted, Discrete Near Stall, Distorted, Discrete Near Op-Line, Distorted, Mean -0.525 Near Stall, Distorted, Mean 0 200 400 600 800 0.002 Time Step 0 90 180 270 360 Azimuthal Position of Root LE [deg] Figure 11: Blade axial motion (suc- Figure 10: Static modal force: Modal tion side tip, midchord) using both one- force on rigid blades, due to unsteady way prescribed motion procedure and pressure loading of clean and distorted two-way coupled motion procedure with flows. One blade shown.
modal forces zeroed.
motion time histories yields dampings between 1.60% to 1.80%. The positive aero- dynamic damping conservatively indicates flutter stability.
0.0004 0.0002 -0.0002 Normalized Generalized Displacement [-] -0.0004 0 0.01 0.02 0.03 0.04 Time [s] Figure 12: Two-way coupled simulation with decreasing magnitudes of blade os- cillation for all blades, indicative of positive aerodynamic damping. Near op-line, clean-inlet flow, 0 ND initial condition. Static modal force: full revolution, sampled discretely.
4.3 Application to Near-stall and Distorted Flows As discussed previously, the one-way coupling method clearly (and conservatively) shows flutter stability for this fan at part-speed, near op-line, clean inlet conditions.
The flutter-stability question is less clear near stall for both clean and distorted inlet conditions, and there is evidence of periodic instability near op-line with dis- torted inlet conditions. The two-way coupled method shall be applied in further investigation of these fan operating conditions. For each of the following studies, the inlet conditions, exit conditions, fan speed, CFD time step, blade frequency, and blade mode shape are specified identically as for the one-way coupling runs. For N ASA/TM—2014-218447 11 highest fidelity – and with greatest computational expense – the static modal force application in these studies consists of compiling data from one revolution of fan operation with rigid blades, and then sampling this data discretely, matching the annular position of the deflecting blade (the current run) with the annular position of the static modal force database (a previous run with rigid blades but otherwise identical conditions).
4.3.1 Near op-line, sinusoidal inlet distortion For a near op-line condition with a distorted inlet, the zero-nodal diameter pattern exhibits stability by inspection of the “average” blade, but unrolling the annulus reveals that a significant portion of the annulus lends itself to positive (unstable) aerodynamic work computations; see Fig. 6(b). The one-, two-, and three-nodal diameter forward traveling wave patterns (not shown) have similar character in their flutter assessment via the one-way coupling method. Applying the two-way coupled analysis to the zero-nodal diameter pattern near op-line, it is readily evident that the blade is stable under these distorted inlet conditions; see Fig. 13. Comparison of Figures 12 and 13 shows little difference in damping for the zero-nodal diameter pattern, near op-line, between clean and distorted flows; these are nearly-identically stable conditions.
4.3.2 Near stall, clean and distorted inlets Near stall, uncertainties regarding flutter stability arise in both clean and distorted flows with the one-way coupled method, due to the changes in sign of the aerody- namic work from blade to blade and from cycle to cycle; see Fig. 7(a). The one-way coupled method’s “average” blade indicates negative aerodynamic damping (flutter instability) of -0.068% and -0.070% for clean and distorted flows, respectively, near stall (see Fig. 8). In Figures 14(a) and (b), displacements of all blades are plotted for both clean and distorted inlets for 1 ND FTW near stall. Figure 14(a) shows blade-to-blade inconsistency about the rotor with a small global trend of increasing blade displacements. Figure 14(b) shows great blade-to-blade inconsistency about the rotor and a greater global trend in the growth of blade displacements, indicating flutter instability.
The two-nodal diameter forward traveling wave pattern is assessed to be min- imally stable near stall in both clean (+0.49% damping) and distorted (+0.09% damping) flows from the “average” blade computations using the one-way coupled method (see Fig. 8). Like the one-nodal diameter forward traveling wave, this two- nodal diameter pattern shows sign-changes about the annulus in its one-way coupled aerodynamic work computations; see Fig. 7(b). Figures 15(a) and (b) plot the dis- placements of all blades in both clean and distorted flows near stall upon application of the two-way coupled method. Like the 1 ND FTW case in clean flow, Fig. 15(a) shows blade-to- blade nonuniformity of displacements in clean flow; however, for the two-nodal diameter forward traveling wave pattern in clean flow, sporadic increases seem to damp again, as would be expected for the flutter-stable system predicted by the “average” blade using the one-way coupled method. For the 2 ND FTW in distorted flow, Fig. 15(b), the blades again exhibit nonuniformity of displacement about the annulus, but now the displacements seem to grow. While the one-way N ASA/TM—2014-218447 12 coupled method’s “average” blade showed minimal stability in clean and distorted flows for 2 ND FTW near stall, the two-way coupled method shows stability near stall in clean flow but instability near stall in distorted flow.
0.0006 0.0004 0.0002 -0.0002 Normalized Generalized Displacement [-] -0.0004 -0.0006 0 0.01 0.02 0.03 0.04 Time [s] Figure 13: Decreasing magnitudes of oscillation of all blades, indicating flutter sta- bility. Near op-line, distorted inlet, 0 ND initial condition. Static modal force: full revolution sampled discretely.
0.0008 0.0008 0.0006 0.0006 0.0004 0.0004 0.0002 0.0002 0 0 -0.0002 -0.0002 -0.0004 -0.0004 Normalized Generalized Displacement [-] Normalized Generalized Displacement [-] -0.0006 -0.0006 -0.0008 -0.0008 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 Time [s] Time [s] (a) Near stall, clean inlet, 1 ND FTW, all (b) Near stall, distorted inlet, 1 ND FTW, blades. all blades Figure 14: Comparing blade displacements at near-stall, 1 ND FTW, for both clean and distorted inlets. Static modal force: full revolution sampled discretely.
5 Conclusions and Future Work
An initial comparison between one-way and two-way coupled flutter analysis meth- ods implemented in TURBO has been performed. An existing fan, for which TURBO’s one-way coupled flutter analysis method has been experimentally corrob- orated, was used as the computational test article. The fan’s structural dynamics characterstics (mode shape and natural frequency) were altered, and the fan was run at part-speed conditions with intention of simulating stable and unstable flutter environments. The fan was provided with both axisymmetric inlet total pressure as well as a sinusoidal variation of total pressure about the annulus, such that the N ASA/TM—2014-218447 13 0.0008 0.0008 0.0006 0.0006 0.0004 0.0004 0.0002 0.0002 0 0 -0.0002 -0.0002 -0.0004 -0.0004 Normalized Generalized Displacement [-] Normalized Generalized Displacement [-] -0.0006 -0.0006 -0.0008 -0.0008 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 Time [s] Time [s] (a) Near stall, clean inlet, 2 ND FTW, all (b) Near stall, distorted inlet, 2 ND FTW, blades. all blades Figure 15: Comparing blade displacements at near-stall, 2 ND FTW, for both clean and distorted inlets. Static modal force: full revolution sampled discretely.
integrated, mean flow parameters for the distorted flow were equivalent to the clean flow, and thus the circumferential uniformity (or nonuniformity) of the inlet flow would be the only difference between the two flow conditions.
As had been demonstrated previously, TURBO’s one-way coupling method pro- vided consistent work-per-cycle computations for all blades across the annulus in clean flow, while distorted flow caused great variation in aerodynamic work com- putations from blade to blade and from cycle to cycle. The “average” blade can be conservatively applied to assess flutter stability in clean flow, provided that all blades’ aerodynamic work computations indicate stability. Multiple cases were ob- served wherein the “average” blade was stable, but the full annulus of blades would oscillate between stability and instability from blade to blade and from cycle to cycle in an orderly repeating sequence.
The two-way coupling method corroborated the one-way coupling method in assessing flutter stability in clean flow near op-line. The two-way coupling method was then used to investigate the flutter stability of clean flow near stall as well as distorted flow at both near op-line and near stall conditions. In most cases, the two- way coupling corroborated the adequacy of the one-way coupling method’s “average” blade assessment of flutter stability. In general, the distorted flow was slightly less stable (or slightly more unstable) than its clean flow counterpart. However, for 1 ND FTW and 2 ND FTW at near stall conditions, the “average” blade’s (one-way coupling) aerodynamic work oscillated between small values positive and negative in clean flow and among large values positive and negative in distorted flow. The two- way coupled method showed blade-to-blade variation at these conditions, but 1 ND FTW seemed to be minorly unstable in clean flow and more unstable in distorted flow. Most notably, the two-way coupled analysis of 2 ND FTW showed to be stable in clean flow but unstable in distorted flow.
While instabilities seem to have been demonstrated, their magnitude was quite small, and these must be investigated further. Refined gridding, refined time- stepping, and longer simulation times would be conducive to more in-depth study of the fluid mechanics’ role in these potential instabilities. The time-step for these N ASA/TM—2014-218447 14 simulations, though constant between one-way and two-way coupling simulations, was rather coarse. Different fluid-structure coupling stencils exist which could prove critical when operating so closely to neutral stability. All two-way coupled flut- ter simulations discussed herein used a full revolution of “static modal force” data and sampled it discretely because it was thought to offer the highest fidelity (at the highest computational cost) of the fluid-structure interaction; the “static modal force” data is generated by running a rigid-blade simulation with identical param- eters to the subsequent elastic-blade simulation of primary interest. NASA’s work in researching developing BLI propulsion systems continues, and in time there will be experimental aeromechanics data of BLI devices to complement and validate computational method development.
References
1. Smith, L. H.: Wake Ingestion Propulsion Benefit. AIAA Journal of Propulsion and Power , vol. 9, no. 1, Jan–Feb 1993, pp. 74–82.
2. Daggett, D.; Kawai, R.; and Friedman, D.: Blended Wing Body System Studies: Boundary Layer Ingerstion Inlets with Active Flow Control. NASA/CR–2003- 212670, Dec. 2003.
3. Kawai, R.; Friedman, D.; and Serrano, L.: Blended Wing Body (BWB) Bound- ary Layer Ingestion (BLI) Inlet Configuration and System Studies. NASA/CR– 2006-214534, Dec. 2006.
4. Plas, A. P.; Sargeant, M. A.; Madami, V.; Crichton, D.; Greitzer, E. M.; Hynes, T. P.; and Hall, C. A.: Performance of a Boundary Layer Ingesting (BLI) Propulsion System. 45th AIAA Aerospace Sciences Meeting and Exhibit , no.
2007-AIAA-0450, Jan. 2007.
5. Nickol, C. L.: Silent Aircraft Initiative Concept Risk Assessment. NASA/TM– 2008-215112, Feb. 2008.
6. Nickol, C. L.; and McCullers, L. A.: Hybrid Wing Body Configuration System Studies. 47th AIAA Aerospace Sciences Meeting Including The New Horizons Forum and Aerospace Exposition , no. 2009-AIAA-0931, Jan. 2009.
7. Tillman, T. G.; Hardin, L. W.; Moffitt, B. A.; Sharma, O. P.; Lord, W. K.; Berton, J.; and Arend, D.: System-Level Benefits of Boundary Layer Ingesting Propulsion. AIAA 49th Aerospace Sciences Meeting, Jan. 2011. Invited presen- tation.
8. Herrick, G. P.: Effects of Inlet Distortion on Aeromechanical Stability of a Forward-Swept High-Speed Fan. 46th AIAA/ASME/SAE/ASEE Joint Propul- sion Conference & Exhibit , no. 2010-AIAA-6711, July 2010.
9. Carstens, V.; Kemme, R.; and Schmitt, S.: Coupled simulation of flow-structure interaction in turbomachinery. Aerospace Science and Technology , vol. 7, 2003, pp. 298–306.
N ASA/TM—2014-218447 15 10. Gnesin, V.; and Rz¸ adkowski, R.: A Coupled Fluid-Structure Analysis for 3-D Inviscid Flutter of IV Standard Configuration. Journal of Sound and Vibration , vol. 251, no. 2, 2002, pp. 315–327.
11. Marshall, J. G.; and Imregun, M.: An Analysis of the Aeroelastic Behaviour of a Typical Fan Blade with Emphasis on Flutter Mechanisms. International Gas Turbine and Aeroengine Congress & Exposition , no. ASME 96-GT-78.
12. Br´ eard, C.; Vahdati, M.; Sayma, A. I.; and Imregun, M.: An Integrated Time- Domain Aeroelasticity Model for the Prediction of Fan Forced Response due to Inlet Distortion. Journal of Engineering for Gas Turbines and Power , vol. 124, no. 1, 2002, pp. 196–208.
13. Zhu, J.; and Shih, T. H.: CMOTT Turbulence Model for NPARC. NASA/CR– 1997-204143, Aug. 1997.
14. Chen, J.-P.; and Whitfield, D. L.: Navier-Stokes Calcuations for the Unsteady Flowfield of Turbomachinery. 1990-AIAA-0676, Jan. 1990.
15. Chen, J.-P.; and Briley, W. R.: A Parallel Flow Solver for Unsteady Multiple Blade Row Turbomachinery Simulations. Proceedings of ASME TURBO EXPO 2001 , no. GT-2001-0348, June 2001.
16. Srivastava, R.; Bakhle, M. A.; and Keith, Jr., T. G.: Numerical Simulation of Aerodynamic Damping for Flutter Analysis of Turbomachinery Blade Rows.
AIAA Journal of Propulsion and Power , vol. 19, no. 2, Mar. 2003.
17. Srivastava, R.; Bakhle, M. A.; Keith, Jr., T. G.; and Stefko, G. L.: Flutter Analysis of a Transonic Fan. Proceedings of ASME TURBO EXPO 2002 , no.
GT2002-30319, June 2002.
18. Bakhle, M. A.; Liu, J. S.; Panovsky, J.; Keith, Jr., T. G.; and Mehmed, O.: Cal- culation and Correlation of the Unsteady Flowfield in a High Pressure Turbine.
Proceedings of ASME TURBO EXPO 2002 , no. GT2002-30322, June 2002.
19. Panovsky, J.; James, D. K.; and Hassan, K. K.: An Assessment of Fan Flutter Using TURBO. 5th National Turbine Engine High Cycle Fatigue Conference , Universal Technology Corporation, Dayton, Ohio, Apr. 2000.
20. Sanders, A. J.: Stall Flutter Assessment of an Advanced Design Transonic Fan Using TURBO-AE. 6th National Turbine Engine High Cycle Fatigue Confer- ence , Universal Technology Corporation, Dayton, Ohio, Apr. 2001.
21. Panovsky, J.; Liu, J. S.; and Bakhle, M. A.: Comparisons of Experimental and Computational Forced Response in a High Pressure Turbine. 7th National Tur- bine Engine High Cycle Fatigue Conference , Universal Technology Corporation, Dayton, Ohio, Apr. 2002.
22. Kielb, R. E.; Barter, J. W.; Thomas, J. P.; and Hall, K. C.: Blade Excitation by Aerodynamic Instabilities. Proceedings of ASME TURBO EXPO 2003 , no.
GT2003-38634, June 2003.
N ASA/TM—2014-218447 16 23. Sanders, A. J.; Hassan, K. K.; and Rabe, D. C.: Experimental and Numeri- cal Study of Stall Flutter in a Transonic Low-Aspect Ratio Fan Blisk. ASME Journal of Turbomachinery , vol. 126, no. 1, 2004, pp. 166–174.
24. Sanders, A. J.: Nonsynchronous Vibration (NSV) due to a Flow-Induced Aero- dynamic Instability in a Composite Fan Stator. ASME Journal of Turbomachin- ery , vol. 127, no. 2, 2005, pp. 412–421.
25. Bakhle, M. A.; Reddy, T. S. R.; and Stefko, G. L.: Comparisons of Flutter Analyses for an Experimental Fan. NASA/TM–2010-216221, May 2010.
26. Bakhle, M. A.; Reddy, T. S. R.; Herrick, G. P.; Shabbir, A.; and Flo- rea, R.: Aeromechanics Analysis of a Boundary Layer Ingesting Fan. 48th AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit , no. AIAA- 2012-3995, 2012.
27. Berrier, B. L.; Carter, M. B.; and Allan, B. G.: High Reynolds Number Inves- tigation of a Flush-Mounted, S-Duct Inlet With Large Amounts of Boundary Layer Ingestion. NASA/TP–2005-213766, Sept. 2005.
28. Carta, F. O.: Coupled Blade-Disk Shroud Flutter Instabilities in Turbojet En- gine Rotors. ASME Journal of Engineering for Power , July 1967, pp. 419–426.
29. Vasanthakumar, P.: Computation of Aerodynamic Damping for Flutter Analy- sis of a Transonic Fan. Proceedings of ASME TURBO EXPO 2011 , no. GT2011- 45697, June 2011.
30. Newmark, N. M.: A method of computation for structural dynamics. ASCE Journal of Engineering Mechanics , vol. 8, 1959, pp. 67–94.
31. Farhat, C.; and Lesoinne, M.: Two efficient staggered algorithms for the serial and parallel solution of three-dimensional nonlinear transient aeroelastic prob- lems. Computational Methods in Applied Mechanics and Engineering , vol. 182, 2000, pp. 499–515.
N ASA/TM—2014-218447 17