Research Papers

Automatic Design Optimization of Profiled Endwalls Including Real Geometrical Effects to Minimize Turbine Secondary Flows

[+] Author and Article Information
Shahrokh Shahpar

CFD Methods, DSE,
Rolls-Royce plc.,
Derby DE24 8BJ, UK
e-mail: shahrokh.shahpar@rolls-royce.com

Stefano Caloni, Laurens de Prieëlle

CFD Methods, DSE,
Rolls-Royce plc.,
Derby DE24 8BJ, UK

1Corresponding author.

Contributed by the International Gas Turbine Institute (IGTI) of ASME for publication in the JOURNAL OF TURBOMACHINERY. Manuscript received August 29, 2016; final manuscript received November 24, 2016; published online March 7, 2017. Editor: Kenneth Hall.

J. Turbomach 139(7), 071010 (Mar 07, 2017) (11 pages) Paper No: TURBO-16-1214; doi: 10.1115/1.4035510 History: Received August 29, 2016; Revised November 24, 2016

This paper presents a novel optimization methodology based on both adjoint sensitivity analysis and trust-based dynamic response surface modeling to improve the performance of a modern turbine of a large civil aero-engine in the presence of high-fidelity geometry configurations. The system has been applied to the nonaxisymmetric hub and tip endwall optimization of a high-pressure turbine stage making use of multirow 3D simulations, parametric modeling, and rapid meshing of real geometry features such as rim seals and modeling of film cooling flows. It has been shown in previous papers that improvements gained using simplified models of the stage are lost when applying the high-fidelity geometry configuration. New results presented in this paper indicate that controlling the purge flow that exits the disk space through the rim seal at the hub of the main annulus is more significant than the reduction of secondary flows in the main passage. For a given rim sealing mass flow rate and whirl velocity, the nonaxisymmetric endwalls are optimized such that the detrimental impact of the sealing flow on the turbine performance is reduced, and hence, the stage efficiency is significantly increased. The traditional optimization approaches based on evolutionary methods or even sequential modifications for defining the endwalls shape are computationally demanding. Since turbomachinery industry continuously strive to reduce the design cycle time, in particular when high-fidelity 3D computational fluid dynamics (CFD) is used, the main body of this paper outlines the novel methods developed to produce a practical design in a very aggressively short design cycle time.

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


Denton, J. , and Pullan, G. , 2012, “ A Numerical Investigation Into the Sources of Endwall Loss in Axial Flow Turbines,” ASME Paper No. GT2012-69173.
Langston, L. S. , 2001, “ Secondary Flows in Axial Turbines—A Review,” Heat and Mass Transfer in Gas Turbine Systems (Annals of the New York Academy of Sciences, Vol. 943), R. J. Goldstein , ed., New York Academy of Sciences, New York, pp. 11–26.
Rose, M. G. , 1995, “ Gas Turbine Engine Turbine,” U.S. Patent No. 5,466,123.
Harvey, N. W. , Rose, M. G. , Shahpar, S. , Taylor, M. D. , Hartland, J. , and Gregory-Smith, D. G. , 2000, “ Non-Axisymmetric Turbine End Wall Design: Part I. Three-Dimensional Design System,” ASME J. Turbomach., 122(2), pp. 278–285. [CrossRef]
Brennan, G. , Harvey, N. W. , Rose, M. G. , Fomison, N. , and Taylor, M. D. , 2001, “ Improving the Efficiency of the Trent 500 HP Turbine Using Non-Axisymmetric End Walls: Part 1. Turbine Design,” ASME Paper No. 2001-GT-0444.
Harvey, N. , Brennan, G. , Newman, D. A. , and Rose, M. G. , 2002, “ Improving Turbine Efficiency Using Non-Axisymmetric Endwalls: Validation in the Multi-Row Environment and With Low Aspect Ratio Blading,” ASME Paper No. 2002-GT-30337.
Gonzalez, P. , and Lantero, M. , 2006, “ Low Pressure Turbine Design for Rolls-Royce TRENT 900 Turbofan,” ASME Paper No. GT2006-90997.
Ingram, G. , Gregory-Smith, D. , and Harvey, N. , 2005, “ Investigation of a Novel Secondary Flow Feature in a Turbine Cascade With Endwall Profiling,” ASME J. Turbomach., 127(1), pp. 209–214. [CrossRef]
Turgut, O. H. , and Camci, C. , 2013, “ Influence of Leading Edge Fillet and Non-Axisymmetric Contoured Endwall on Turbine NGV Exit Flow Structure and Interactions With the Rim Seal,” ASME Paper No. GT2013-95843.
Schuepbach, P. , Abhari, R. S. , Rose, M. G. , and Gier, J. , 2009, “ Influence of Rim Seal Purge Flow on Performance of an Endwall-Profiled Axial Turbine,” ASME Paper No. GT2009-59653.
Lott, P. T. , Hills, N. J. , Chew, J. W. , Scanlon, T. , and Shahpar, S. , 2009, “ High Pressure Turbine Stage Endwall Profile Optimization for Performance and Rim Seal Effectiveness,” ASME Paper No. 2009-59923.
Chilla, M. , Hodson, H. , and Newman, D. , 2013, “ Unsteady Interaction Between Annulus and Turbine Rim Seal Flows,” ASME J. Turbomach., 135(5), p. 051024. [CrossRef]
Milli, A. , and Shahpar, S. , 2008, “ Full Parametric Design System to Improve the Stage Efficiency of a High-Fidelity HP Turbine Configuration,” ASME Paper No. GT2008-51348.
Aulich, A.-L. , and Goerke, D. , 2013, “ Multidisciplinary Automated Optimization Strategy on a Counter Rotating Fan,” ASME Paper No. GT2013-94259.
Luo, J. , and Liu, F. , 2013, “ Multi-Point Design Optimization of a Transonic Compressor Blade Using an Adjoint Method,” ASME Paper No. GT2013-95155.
Shahpar, S. , and Caloni, S. , 2013, “ Adjoint Optimization of a High Pressure Turbine Stage for a Lean-Burn Combustion System,” 10th European Conference on Turbomachinery Fluid Dynamics and Thermodynamics (ETC 10), Lappeenranta, Finland, Apr. 15–19.
Verstraete, T. , 2012, “ Multidisciplinary Optimization of Turbomachinery Components Using Differential Evolution,” Introduction to Optimization and Multidisciplinary Design in Aeronautics and Turbomachinery (VKI Lecture Series), J. Periaux , and L. S. Verstraete , eds., von Karman Institute, Brussels, Belgium.
Giannakoglou, K. C. , Asouti, V. G. , Kyriacou, S. A. , and Trompoukis, X. S. , 2012, “ Hierarchical, Metamodel-Assisted Evolutionary Algorithms, With Industrial Applications,” Introduction to Optimization and Multidisciplinary Design in Aeronautics and Turbomachinery (VKI Lecture Series), J. Periaux , and L. S. Verstraete , eds., von Karman Institute Brussels, Belgium.
Shahpar, S. , 2012, “ Optimization Strategies Used in Turbomachinery Design From an Industrial Perspective,” Introduction to Optimization and Multidisciplinary Design in Aeronautics and Turbomachinery (VKI Lecture Series), J. Periaux , and L. S. Verstraete , eds., von Karman Institute, Brussels, Belgium.
Shahpar, S. , and Lapworth, L. , 2003, “ PADRAM: Parametric Design and Rapid Meshing System for Turbomachinery Optimization,” ASME Paper No. GT2003-38698.
Milli, A. , and Shahpar, S. , 2012, “ PADRAM: Parametric Design and Rapid Meshing System for Complex Turbomachinery Configurations,” ASME Paper No. GT2012-69030.
Shahpar, S. , 2005, “ SOPHY: An Integrated CFD Based Automatic Design Optimization System,” 17th International Symposium on Air Breathing Engines, Munich, Germany, Sept. 4–9, Paper No. ISABE-2005-1086.
Shahpar, S. , 2007, “ Towards Robust CFD Based Design Optimization of Virtual Engine,” NATO, Athens, Greece, Document No. RTO–AVT-147.
Lapworth, L. , 2004, “ Hydra-CFD: A Framework for Collaborative CFD Development,” International Conference on Scientific and Engineering Computation (IC-SEC), Singapore, June 30–July 2.
Shahpar, S. , 2002, “ SOFT: A New Design and Optimization Tool for Turbomachinery,” Evolutionary Methods for Design, Optimization and Control, K. Giannakoglou , D. T. Tsahalis , J. Periaux , and T. Fogarty , eds., International Center for Numerical Methods in Engineering (CIMNE), Barcelona, Spain.
Jameson, A. , 1988, “ Aerodynamic Design Via Control Theory,” J. Sci. Comput., 3(3), pp. 233–260. [CrossRef]
Jameson, A. , Pierce, N. A. , and Martinelli, L. , 1998, “ Optimum Aerodynamic Design Using the Navier–Stokes Equations,” Theoretical and Computational Fluid Dynamics, Vol. 10, Springer-Verlag, Berlin, pp. 213–237.
Reuther, J. J. , Alonso, J. J. , Jameson, A. , Rimlinger, M. J. , and Saunders , 1999, “ Constrained Multipoint Aerodynamic Shape Optimization Using and Adjoint Formulation and Parallel Computers, Part 1,” J. Aircr., 36(1), pp. 51–60. [CrossRef]
Martins, J. R. R. A. , Alonso, J. J. , and Reuther, J. J. , 2004, “ High-Fidelity Aero-Structural Design Optimization of a Supersonic Business Jet,” J. Aircr., 41(3), pp. 523–530. [CrossRef]
Lohner, R. , Soto, O. , and Yang, C. , 2003, “ An Adjoint-Based Design Methodology for CFD Optimization Problems,” AIAA Paper No. 2003-299.
Mohammadi, B. , and Pironneau, O. , 2001, Applied Shape Optimization for Fluids, Oxford University Press, Oxford, UK.
Giles, M. B. , and Pierce, N. A. , 2000, “ An Introduction to the Adjoint Approach to Design,” Flow Turbul. Combust., 65, pp. 393–415. [CrossRef]
Coral, R. , and Gisbert, F. , 2006, “ Profiled End-Wall Design Using an Adjoint Navier–Stokes Solver,” ASME Paper No. GT2006-90650.
Wang, D. X. , and He, L. , 2008, “ Adjoint Aerodynamic Design Optimization for Blades in Multi-Stage Turbomachines, Part I,” ASME Paper No. GT2008-50208.
Wang, D. X. , and He, L. , 2008, “ Adjoint Aerodynamic Design Optimization for Blades in Multi-Stage Turbomachines, Part II,” ASME Paper No. GT2008-50209.
Wang, D. X. , and He, L. , 2009, “ Concurrent Aerodynamic-Aeromechanic Design Optimization for Turbomachinery Blades Using Adjoint Method,” ASME Paper No. GT2009-59240.
Frey, C. , Kersken, H.-P. , and Nurnberger, D. , 2009, “ The Discrete Adjoint of a Turbomachinery Solver,” ASME Paper No. GT 2009-59062.
Giles, M. B. , Duta, M. C. , Muller, J.-D. , and Pierce, N. A. , 2003, “ Algorithm Developments for Discrete Adjoint Methods,” AIAA J., 41(2), pp. 198–205. [CrossRef]
Marta, A. C. , Shankaran, S. , and Stein, A. , 2010, “ Blade Shape Optimization Using a RANS Discrete Adjoint Solver,” 2nd International Conference on Engineering Optimization, Lisbon, Portugal, Sept. 6–9, Paper No. 1410.
Duta, M. C. , Shahpar, S. , and Giles, M. B. , 2007, “ Turbomachinery Design Optimization Using Automatic Differentiated Adjoint Code,” ASME Paper No. GT2007-28329.
Schittkowski, K. , 1986, “ NLPQL, A Fortran Subroutine for Solving Constrained Non-Linear Programming Problems,” Ann. Oper. Res., 5, pp. 485–500. [CrossRef]
Shahpar, S. , and Caloni, S. , 2012, “ Aerodynamic Optimization of High Pressure Turbines for Lean-Burn Combustion System,” ASME Paper No. GT2012-69228.
Polynkin, A. , Toropov, V. , and Shahpar, S. , 2010, “ Multidisciplinary Optimization of Turbomachinery Based on Metamodel Built by Genetic Programming,” AIAA Paper No. 2010-9397.
Shahpar, S. , 2010, “ High Fidelity Multi-Stage Design Optimization of Multi-Stage Turbine Blades Using a Mid-Range Approximation Method,” AIAA Paper No. 2010-9367.
Popovic, I. , and Hodson, H. P. , 2012, “ The Effects of a Parametric Variation of the Seal Geometry on the Interaction Between Hub Leakage and Mainstream Flows in HP Turbines,” ASME Paper No. GT2012-68025.


Grahic Jump Location
Fig. 1

Schematic flowchart of the SOPHY system

Grahic Jump Location
Fig. 2

Schematic of zooming of trust regions in MAM

Grahic Jump Location
Fig. 3

Full model view and HP + IP turbine rows

Grahic Jump Location
Fig. 4

Computational model of the HP rotor

Grahic Jump Location
Fig. 5

PADRAM mesh, for NGV (top left), zoom on rotor shingling (top right), TE coolant slot (bottom left), and a large nonaxisymmetric endwall mesh perturbation

Grahic Jump Location
Fig. 6

Axial stations set for the rotor endwall profiling

Grahic Jump Location
Fig. 7

Turbine efficiency adjoint sensitivity map—red: raise the platform, green: no sensitivity, and blue: lower the platform

Grahic Jump Location
Fig. 8

Turbine reaction adjoint sensitivity map—red: raise the platform, green: no sensitivity, and blue: lower the platform

Grahic Jump Location
Fig. 9

Contour of change in radius from axisymmetric: (a) adjoint designed and (b) MAM optimized

Grahic Jump Location
Fig. 10

Convergence history of the MAM optimizer

Grahic Jump Location
Fig. 11

(a) Typical convergence history of residuals and (b) oscillation of efficiency

Grahic Jump Location
Fig. 12

Relative total pressure and relative whirl angle at rotor exit plane. Plotted every four iterations of one period.

Grahic Jump Location
Fig. 13

Rotor lift plot at 10% height. Plotted every four iterations of one period.

Grahic Jump Location
Fig. 14

Contribution of various bumps to the overall loss reductions

Grahic Jump Location
Fig. 15

Static pressure above the rim seal and close to the endwall (solid lines), and circumferential variation of the average static pressure over the cavity rim (dashed lines) for one pitch

Grahic Jump Location
Fig. 16

Cavity rim vortex colored by relative total temperature for the (a) datum geometry (red arrow: main annulus gas path) and (b) for the MAM optimum design (blue arrow: purge coolant flow)

Grahic Jump Location
Fig. 17

Contour of relative total temperature in a circumferential plane just behind the FLMP bump: (a) axisymmetric case and (b) MAM optimized

Grahic Jump Location
Fig. 18

Endwall streamlines colored by static pressure ((a) and (b)). Meridional view of “core coolant” and “recirculated” flow ((c) and (d)). Limiting streamlines on the hub endwall, colored by static pressure ((e) and (f)). Left: axisymmetric endwall and right: MAM optimized PEW.



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