Research Papers

On the Impact of Geometric Variability on Fan Aerodynamic Performance, Unsteady Blade Row Interaction, and Its Mechanical Characteristics

[+] Author and Article Information
R. Schnell

German Aerospace Center e.V.—DLR,
Institute of Propulsion Technology,
Linder Hoehe,
Cologne 51147, Germany
e-mail: rainer.schnell@dlr.de

T. Lengyel-Kampmann, E. Nicke

German Aerospace Center e.V.—DLR,
Institute of Propulsion Technology,
Linder Hoehe,
Cologne 51147, Germany

1Corresponding author.

Contributed by the International Gas Turbine Institute (IGTI) of ASME for publication in the JOURNAL OF TURBOMACHINERY. Manuscript received October 17, 2013; final manuscript received March 10, 2014; published online April 18, 2014. Assoc. Editor: Seung Jin Song.

J. Turbomach 136(9), 091005 (Apr 18, 2014) (14 pages) Paper No: TURBO-13-1236; doi: 10.1115/1.4027218 History: Received October 17, 2013; Revised March 10, 2014; Accepted March 16, 2014

The focus of the present study is to assess and quantify the uncertainty in predicting the steady and unsteady aerodynamic performance as well as the major mechanical characteristics of a contrarotating turbofan, primarily due to geometric variations stemming from the manufacturing process. The basis of this study is the optically scanned blisk of the first rotor, for which geometric variations from blade to blade are considered. In a first step, selected profile sections of the first rotor were evaluated aerodynamically by applying the 2D coupled Euler/boundary-layer solver mises. Statistical properties of the relevant flow quantities were calculated firstly based on the results of the nine manufactured blades. In a second step, the geometric variations were decomposed into their corresponding eigenforms by means of principal component analysis (PCA). These modes were the basis for carrying out Monte Carlo (MC) simulations in order to analyze in detail the blade's aerodynamic response to the prescribed geometric variations. By means of 3D-computational fluid dynamics (CFD) simulations of the entire fan stage for all the nine scanned rotor 1 blade geometries, the variation of the overall stage performance parameters will be quantified. The impact of the instrumentation will be discussed, here partly doubling the standard deviation of the major performance indicators for the instrumented blades and also triggering a premature laminar/turbulent transition of the boundary layer. In terms of the unsteady blade row interaction, the standard deviation of the resulting blade pressure amplitude shall be discussed based on unsteady simulations, taking advantage of a novel harmonic balance approach. It will be shown that the major uncertainty in terms of the predicted blade pressure amplitude is in the aft part of the front rotor and results from upstream shock/blade interaction. Apart from the aerodynamic performance, an analysis of the mechanical properties in terms of Campbell characteristics and eigenfrequencies was carried out for each of the scanned blades of rotor 1, reflecting the frequency scattering of each eigenmode due to geometric variability.

Copyright © 2014 by ASME
Your Session has timed out. Please sign back in to continue.


Garzon, V. E., 2002, “Probabilistic Aerothermal Design of Compressor Airfoils,” Ph.D. thesis, MIT, Cambridge, MA.
Goodhand, M. N., Miller, R. J., and Lung, H. W., 2012, “The Sensitivity of 2D Compressor Incidence to In-Service Geometric Variation,” ASME Paper No. GT2012-68633. [CrossRef]
Lange, A., Voigt, M., Vogeler, K., Schrapp, H., Johann, E., and Gümmer, V., 2011, “Impact of Manufacturing Variability and Nonaxisymmetry on High-Pressure Compressor Stage Performance,” ASME J. Eng. Gas Turbines Power, 143(3), p. 032504. [CrossRef]
Lamb, C. M., 2005, “Probabilistic Performance-Based Geometric Tolerancing of Compressor Blades,” Master's thesis, MIT, Cambridge, MA.
Giebmanns, A., Backhaus, J., Frey, C., and Schnell, R., 2013, “Compressor Leading Edge Sensitivities and Analysis With an Adjoint Flow Solver,” ASME Paper No. GT2013-94427. [CrossRef]
Aulich, A.-L., Goerke, D., Blocher, M., Nicke, E., and Kocian, F., 2013, “Automated Optimization Strategy on a Counter Rotating Fan,” ASME Paper No. GT2013-94259. [CrossRef]
Schnell, R., Giebmanns, A., Nicke, E., and Dabrock, T. C., 2009, “Aerodynamic Analysis of a Fan for Future Ultra-High-Bypass-Ratio Aero Engines,” 19th International Symposium on Air Breathing Engines (ISABE), Montreal, Canada, September 7–11, Paper No. 2009-1149.
Talbotec, J., and Vernet, M., 2010, “Snecma Counter Rotating Fan Aerodynamic Design: Logic & Tests Results,” 27th Congress of International Council of the Aeronautical Science, Nice, France, September 19–24, Paper ICAS 2010-4.1.2.
Lengyel, T., Voß, C., Schmidt, T., and Nicke, E., 2009, “Design of a Counter Rotating Fan—An Aircraft Engine Technology to Reduce Noise and CO2-Emissions,” 19th International Symposium on Air Breathing Engines (ISABE), Montreal, Canada, September 7–11, Paper No. 2009-1267.
Mileshin, V., Popovyan, A., Korznev, V., Maximov, V., and Khaletskiy, Y., 2009, “C-3A Aero-Acoustic Test Facility of CIAM and Use for VITAL Counter Rotating Fan Tests (CRTF1 Application),” Proceedings of The VITAL Final Workshop, Budapest, Hungary, March 9–10, Paper No. 19.
Garzon, V. E., and Darmofal, D. L., 2003, “Impact of Geometric Variability on Axial Compressor Performance,” ASME J. Turbomach., 125(4), pp. 692–703. [CrossRef]
Holtzhausen, S., Schreiber, S., Schöne, C., Stelzer, R., Heinze, K., and Lange, A., 2009, “Highly Accurate Automated 3D Measuring and Data Conditioning for Turbine and Compressor Blades,” ASME Paper No. GT2009-59902. [CrossRef]
Kumar, A., Nair, P. B., Keane, A. J., and Shahpar, S., 2005, “Probabilistic Performance Analysis of Eroded Compressor Blades,” ASME Paper No. PWR2005-50070. [CrossRef]
Heinze, K., Schrape, S., Voigt, M., and Vogeler, K., 2012, “Probabilistic Endurance Level Analyses of Compressor Blades,” CEAS Aeronaut. J., 2012, 3(1), pp. 55–65. [CrossRef]
Drela, M., and Youngren, H., 1998, A User's Guide to Mises 2.53, MIT, Cambridge, MA.
Goodhand, M. N., and Miller, R. J., 2012, “Compressor Leading Edge Spikes: A New Performance Criterion,” ASME J. Turbomach., 133(2), p. 021006. [CrossRef]
Drela, M., 1998, Mises Implementation of Modified Abu-Ghannam/Shaw Transition Criterion, MIT, Cambridge, MA.
Giebmanns, A., Schnell, R., Steinert, W., Hergt, A., Nicke, E., and Werner-Spatz, C., 2012, “Analysing and Optimising Geometrically Degraded Transonic Fan Blades by Means of 2D and 3D Simulations and Cascade Measurements,” ASME Paper No. GT2012-69064. [CrossRef]
Lengyel-Kampmann, T., Bischoff, A., Meyer, R., and Nicke, E., 2012, “Design of an Economical Counter-Rotating Fan: Comparison of the Calculated and Measured Steady and Unsteady Results,” ASME Turbo Expo, Copenhagen, Denmark, Paper No. GT2012-69587. [CrossRef]
Franke, M., Roeber, T., Kuegeler, E., and Ashcroft, G., 2010, “Turbulence Treatment in Steady and Unsteady Turbomachinery Flows,” Fifth European Conference on Computational Fluid Dynamics (ECCOMAS CFD 2010), Lisbon, Portugal, June 14–17.
Giebmanns, A., Schnell, R., and Nicke, E., 2012, “Numerical Analysis of Typical Degradation Characteristics of a Transonic Fan,” German Aerospace Center (DLR), Cologne, Germany, Internal Report No. DLR-IB-325-03-13.
Schlichting, H., 1982, Boundary Layer Theory (Grenzschichttheorie), 8th ed., G.Braun, ed., Wissenschaftliche Bücherei Verlag, Stuttgart, Germany.
White, F. M., 1991, Viscous Fluid Flow, 2nd ed., McGraw-Hill, New York.
McMullen, M. S., 2003, “The Application of Non-Linear Frequency Domain Methods to the Euler and Navier-Stokes Equations,” Ph.D. thesis, Stanford University, Stanford, CA.
Hall, K. C., Thomas, J. P., and Clark, W. S., 2002, “Computation of Unsteady Nonlinear Flows in Cascades Using a Harmonic Balance Technique,” AIAA J., 40(5), pp. 879–886. [CrossRef]
He, L., 1998, “Efficient Approach for Analysis of Unsteady Viscous Flows in Turbomachines,” AIAA J., 36(11), pp. 2005–2012. [CrossRef]
Frey, C., Ashcroft, G., Kersken, H.-P., and Voigt, C., 2014, “A Harmonic Balance Technique for Turbomachinery Applications,” ASME Paper No. GT2014-25230. [CrossRef]
Sinha, A., Hall, B., Cassenti, B., and Hilber, G., 2008, “Vibratory Parameters of Blades From Coordinate Measurement Machine Data,” ASME J. Turbomach., 130(1), p. 011013. [CrossRef]
Beck, J. A., Brown, J. M., Slater, J. C., and Cross, C. J., 2013, “Probabilistic Mistuning Assessment Using Nominal and Geometry Based Mistuning Methods,” ASME J. Turbomach., 135(5), p. 051004. [CrossRef]
Reddy, T. S. R., Mital, S. K., and Stefko, G. L., 2004, “Probabilistic Aeroelastic Analysis of Turbomachinery Components,” NASA Technical Memorandum No. TM-2004-213063.
Joshi, A. G., and Epureanu, B. I., 2012, “Reduced Order Models for Blade-to-Blade Damping Variability in Mistuned Blisks,” ASME J. Vib. Acoust., 134(5), p. 051015. [CrossRef]
Dhondt, G., 2004, The Finite Element Method for Three-Dimensional Thermomechanical Applications, Wiley, Hoboken, NJ.
Siller, U., and Aulich, M., 2011, “High-Dimensional Constrained Multiobjective Optimization of a Fan Stage,” ASME Paper No. GT2011-45618. [CrossRef]
Cailleau, J. M., Talbotec, J., and Farvacque, B., 2007, “VITAL WP2.4 T2.4.4 CRTF2b Specifications,” EU Project VITAL Deliverable.


Grahic Jump Location
Fig. 1

Contrarotating turbo fan rig CRTF2b installed at CIAM compressor test bed C-3 A [10]

Grahic Jump Location
Fig. 2

Measured and calculated performance characteristic of the CRTF2b fan stage at 100%, 83%, and 54% rotational speed [19]

Grahic Jump Location
Fig. 3

Strain gauge on the rotor of blisk 1 (left) and assembly of rotor 1 and 2 including the booster inlet guide vane (IGV) (right); the two blisks comprised of nine rotor 1 blades and 11 rotor 2 blades

Grahic Jump Location
Fig. 4

Color-coded differences between manufactured and nominal geometry (in mm), including strain gauges' position and corresponding wiring

Grahic Jump Location
Fig. 5

Surface-based CAD model of blade 04 including instrumentation and wiring on the pressure side as it was used for the 3D-CFD and CSM analysis of the real geometry (left) and nodewise distance between clean blade 03 and instrumented blade 04 to quantify sensor dimensions (right)

Grahic Jump Location
Fig. 6

Leading-edge detail for different profile sections at 90% span: nominal (circles), all nine scanned blades of blisk 1 (solid lines), mean manufactured geometry, and ≈1000 geometries created based on PCA modes for MC simulations (solid light gray lines); the manufacturing tolerance band of ±0.1 mm is indicated by dashed circles at several locations along the profile

Grahic Jump Location
Fig. 7

Mode shapes of the first four modes from principal component analysis and defined nondimensional parameters to characterize each mode

Grahic Jump Location
Fig. 8

Outflow angles β2 at different incidence angles β1: results for the nominal geometry (solid line) and the mean manufactured geometry (dashed line) as well as results from Monte Carlo simulations based on PCA modes (plotted are the expectancy values as circles as well as corresponding uncertainties within ±2σ as error bars)

Grahic Jump Location
Fig. 9

Empirical CDFs calculated from the nine samples (bullets) as well as from Monte Carlo simulations (solid lines) in comparison with the CDF from normal distribution (dashed line); the corresponding frequency distributions are shown in Figs. 8(a)8(c)

Grahic Jump Location
Fig. 10

Various flow quantities and corresponding variation within ±2σ at different inflow incidence angles (90% span)

Grahic Jump Location
Fig. 11

Relative standard deviation of different resulting flow parameters from MC simulations, each taking into account a different number of PCA modes to describe the geometric variability

Grahic Jump Location
Fig. 12

Isentropic Mach number distribution at two different incidence angles and its variation due to the geometric variability (based on PCA modes as shown in Fig. 6) as a result from MC simulations; indicated is also the variation of the transition point at high flow incidence within a confidence interval of ±2σ (right)

Grahic Jump Location
Fig. 13

Boundary layer edge velocity in the vicinity of the leading edge and calculated spike diffusion factors according to Ref. [16] for all scanned blade sections

Grahic Jump Location
Fig. 24

Blade-surface static pressure amplitudes (first relative blade passing frequency (BPF) harmonic) resulting from unsteady blade-row interaction

Grahic Jump Location
Fig. 23

Blade-surface static pressure contours at design operating conditions for blade 01 (a), pressure fluctuations (amplitude of the first interaction harmonic) on both blades induced by wake/blade interaction (downstream) (b), and resulting pressure amplitude uncertainty (c)

Grahic Jump Location
Fig. 22

Instantaneous flow field at approximately 70% span highlighting the major interaction phenomena: upstream shock/blade interaction visualized by static pressure contours (left) and downstream wake/blade interaction visualized by contours of eddy viscosity (right)—results from HB computations with five harmonics

Grahic Jump Location
Fig. 21

Comparison between stage and rotor isentropic efficiency standard deviation at three selected operating conditions

Grahic Jump Location
Fig. 20

Boundary layer shape parameter H12 for all nine scanned blades at approximately 65% span (r = 180 mm) for the blade suction side (a) and pressure side (b); nondimensional displacement thickness δ1 as well as skin friction coefficient cf on the pressure side (c)

Grahic Jump Location
Fig. 19

Standard deviation of the rotor 1 total pressure at operating conditions 1–6 (see Fig. 16 for their definition on the 100% speed line); the error bars denote the standard error due to the limited number of samples O(σ/n1/2)

Grahic Jump Location
Fig. 18

Standard deviation of the mass flow rate at operating conditions 1–6 (see Fig. 16 for their definition on the 100% speed line)

Grahic Jump Location
Fig. 17

Standard deviation of the rotor 1 isentropic efficiency at operating conditions 1–6 (see Fig. 16 for their definition on the 100% speed line)

Grahic Jump Location
Fig. 16

Rotor 1 total pressure ratio at 100% rotational speed based on 3D-CFD full-stage simulations for all nine scanned rotor 1 blades and derived standard deviation (±2σ)

Grahic Jump Location
Fig. 15

Rotor 1 isentropic efficiency at 100% rotational speed based on 3D-CFD full-stage simulations for all nine scanned rotor 1 blades and derived standard deviation (±2σ)

Grahic Jump Location
Fig. 14

Computational domain and grid details in an S1-plane on the hub for both rotors close to the leading edge (grid comprising of 6.2 × 106 cells in total, 141 cells in spanwise direction)

Grahic Jump Location
Fig. 25

Measured eigenfrequencies from blade excitation tests in comparison with FEM results as obtained in the present study (black bars, both average over all blades) and results with refined FEM model (permas), including the blade root/disk for the nominal geometry (see Ref. [19] for more details)

Grahic Jump Location
Fig. 26

Campbell diagram for all nine scanned blades of the rotor 1 blisk; vertical dashed lines indicate test relevant operating conditions

Grahic Jump Location
Fig. 27

Standard deviation of the eigenfrequency of the first four eigenmodes based on the results for all scanned blades

Grahic Jump Location
Fig. 28

Standard deviation of the eigenfrequency of the first three eigenmodes—results from blade excitation tests (average over three positions) in comparison with FEM results



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In