Research Papers

An Implicit Harmonic Balance Method in Graphics Processing Units for Oscillating Blades

[+] Author and Article Information
Javier Crespo

Technology and Methods Department,
Industria de TurboPropulsores S.A.,
Alcobendas 28108, Madrid, Spain
e-mail: Javier.Crespo@itp.es

Roque Corral

Advanced Engineering Direction,
Industria de TurboPropulsores S.A.,
Alcobendas 28108, Madrid, Spain
e-mail: Roque.Corral@itp.es

Jesus Pueblas

Technology and Methods Department,
Industria de TurboPropulsores S.A.,
Alcobendas 28108, Madrid, Spain
e-mail: Jesus.Pueblas@itp.es

1Also at the Department of Fluid Dynamics and Aerospace Propulsion of the School of Aeronautics, UPM, Madrid 28040, Spain.

Contributed by the International Gas Turbine Institute (IGTI) of ASME for publication in the JOURNAL OF TURBOMACHINERY. Manuscript received July 20, 2015; final manuscript received October 20, 2015; published online November 17, 2015. Editor: Kenneth C. Hall.

J. Turbomach 138(3), 031001 (Nov 17, 2015) (10 pages) Paper No: TURBO-15-1163; doi: 10.1115/1.4031918 History: Received July 20, 2015; Revised October 20, 2015

An implicit harmonic balance (HB) method for modeling the unsteady nonlinear periodic flow about vibrating airfoils in turbomachinery is presented. An implicit edge-based three-dimensional Reynolds-averaged Navier–Stokes equations (RANS) solver for unstructured grids, which runs both on central processing units (CPUs) and graphics processing units (GPUs), is used. The HB method performs a spectral discretization of the time derivatives and marches in pseudotime, a new system of equations where the unknowns are the variables at different time samples. The application of the method to vibrating airfoils is discussed. It is shown that a time-spectral scheme may achieve the same temporal accuracy at a much lower computational cost than a backward finite-difference method at the expense of using more memory. The performance of the implicit solver has been assessed with several application examples. A speed-up factor of 10 is obtained between the spectral and finite-difference version of the code, whereas an additional speed-up factor of 10 is obtained when the code is ported to GPUs, totalizing a speed factor of 100. The performance of the solver in GPUs has been assessed using the tenth standard aeroelastic configuration and a transonic compressor.

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


Saad, Y. , and Schultz, M. , 1986, “ GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems,” SIAM J. Sci. Stat. Comput., 7(3), pp. 856–869. [CrossRef]
Jameson, A. , 1991, “ Time Dependent Calculations Using Multigrid With Applications to Unsteady Flows Past Airfoils and Wings,” AIAA Paper No. 1991-1596.
Verdon, J. , and Casper, J. , 1984, “ A Linearised Unsteady Aerodynamic Analysis for Transonic Cascades,” J. Fluid Mech., 149, pp. 403–429. [CrossRef]
Sbardella, L. , and Imregun, M. , 2001, “ Linearized Unsteady Viscous Turbomachinery Flows Using Hybrid Grids,” ASME J. Turbomach., 123(3), pp. 568–582. [CrossRef]
Corral, R. , Escribano, A. , Gisbert, F. , Serrano, A. , and Vasco, C. , 2003, “ Validation of a Linear Multigrid Accelerated Unstructured Navier–Stokes Solver for the Computation of Turbine Blades on Hybrid Grids,” AIAA Paper No. 2003-3326.
Hall, K. , 1998, “ Computation of Unsteady Nonlinear Flows in Cascades Using a Harmonic Balance Technique,” Kerrebrock Symposium, A Symposium in Honor of Professor Jack L. Kerrebrock's 70th Birthday, Massachusetts Institute of Technology, Cambridge, MA, Jan. 30.
Hall, K. C. , and Crawley, E. F. , 1989, “ Calculation of Unsteady Flows in Turbomachinery Using the Linearized Euler Equations,” AIAA J., 27(6), pp. 777–787. [CrossRef]
He, L. , and Ning, W. , 1998, “ Computation of Unsteady Flows Around Oscillating Blades Using Linear and Nonlinear Harmonic Euler Methods,” ASME J. Turbomach., 120(3), pp. 508–514. [CrossRef]
McMullen, M. , Jameson, A. , and Alonso, J. J. , 2001, “ Acceleration of Convergence to a Periodic Steady State in Turbomachinery Flows,” AIAA Paper No. 2001-0152.
Frey, C. , Ashcroft, G. , Kersken, H. , and Voigt, C. , 2014, “ A Harmonic Balance Technique for Multistage Turbomachinery Applications,” ASME Paper No. GT2014-25230.
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]
Gopinath, A. , and Jameson, A. , 2005, “ Time Spectral Method for Periodic Unsteady Computations Over Two- and Three-Dimensional Bodies,” AIAA Paper No. 05-1220.
Corrigan, A. , Camelli, F. , Löhner, R. , and Wallin, J. , 2009, “ Running Unstructured Grid CFD Solvers on Modern Graphics Hardware,” AIAA Paper No. 2009-4001.
Brandvik, T. , and Pullan, G. , 2007, “ Acceleration of a Two-Dimensional Euler Solver Using Commodity Graphics Hardware,” Proc. Inst. Mech. Eng., Part C, 221(12), pp. 1745–1748. [CrossRef]
Corral, R. , Crespo, J. , and Gisbert, F. , 2004, “ Parallel Multigrid Unstructured Method for the Solution of the Navier–Stokes Equations,” AIAA Paper No. 2004-0761.
Burgos, M. A. , Corral, R. , and Contreras, J. , 2011, “ Efficient Edge-Based Rotor/Stator Interaction Method,” AIAA J., 49(1), pp. 19–31. [CrossRef]
Canuto, C. , Hussaini, M. , Quarteroni, A. , and Zang, T. , 1988, Spectral Methods in Fluid Dynamics (Springer Series in Computational Physics), Springer-Verlag, New York.
Pueblas, J. , Corral, R. , and Gisbert, F. , 2013, “ Selection of Implicit CFD Techniques for Unstructured Grids and Turbomachinery Applications,” AIAA Paper No. 2013-2575.
Erdos, J. , Alzner, E. , and McNally, W. , 1977, “ Numerical Solution of Periodic Transonic Flow Through a Fan Stage,” AIAA J., 15(11), pp. 1559–1568. [CrossRef]
Koya, M. , and Kotake, S. , 1985, “ Numerical Analysis of Fully Three-Dimensional Periodic Flows Through a Turbine Stage,” ASME J. Eng. Gas Turbines Power, 107(4), pp. 945–952. [CrossRef]
He, L. , 1990, “ An Euler Solution for Unsteady Flows Around Oscillating Blades of Turbomachinery,” ASME J. Turbomach., 112(4), pp. 714–722. [CrossRef]
Gisbert, F. , Corral, R. , and Pastor, G. , 2011, “ Implementation of an Edge-Based Navier-Stokes Solver for Unstructured Grids in Graphics Processing Units,” ASME Paper No. GT2011-46224.
Verdon, J. M. , 1987, “ Linearized Unsteady Aerodynamic Theory,” Manual on Aeroelasticity in Axial-Flow Turbomachines, Vol. 1, AG-298, M. F. Platzer , and F. O. Carta , eds., AGARD, Neuilly sur Seine, France, pp. 2–1, 2–30.
Vega, A. , and Corral, R. , 2016, “ The Low Reduced Frequency Limit of Vibrating Airfoils. Part I: Theoretical Analysis,” ASME J. Turbomach., 138(2), p. 021004.
Zhai, Y. , and Bladh, R. , 2011, “ Aeroelastic Stability Assessment of an Industrial Compressor Blade Including Mistuning Effects,” ASME Paper No. GT2011-45800.
Holzinger, F. , Wartzek, F. , Nestle, M. , Schiffer, M. P. , and Leichtfuβ, S. , 2015, “ Self-Excited Blade Vibration Experimentally Investigated in Transonic Compressors—Acoustic Resonances,” ASME Paper No. GT2015-43618.
Corral, R. , and Crespo, J. , 2011, “ Development of an Edge-Based Harmonic Balance Method for Turbomachinery Flows,” ASME Paper No. GT2011-45170.


Grahic Jump Location
Fig. 1

Amplitude (a) and phase (b) of the distribution along chord of the first harmonic of the pressure coefficient for the tenth standard configuration: M = 0.7, k = 0.5, σ = 0 deg, and α = 2 deg

Grahic Jump Location
Fig. 2

Convergence histories of explicit and implicit methods for the inviscid airfoil case

Grahic Jump Location
Fig. 3

First (a) and second (b) harmonics of the unsteady pressure coefficient distribution for different inlet Mach numbers for the tenth standard configuration: k = 0.5, σ = 0 deg, and α = 2 deg

Grahic Jump Location
Fig. 4

First (a) and second (b) harmonics of the unsteady pressure coefficient distribution for different torsion amplitudes: k = 0.5, σ = 0 deg, and M = 0.8

Grahic Jump Location
Fig. 5

Aerodynamic work per cycle as a function of IBPA (a) and distribution along airfoil chord (b), for different torsion amplitudes: k = 0.5, σ = −180 deg, and M = 0.8

Grahic Jump Location
Fig. 6

Time average ((a) and (b)), first harmonic ((c) and (d)), and second harmonic ((e) and (b)) of the pressure distribution for the HB ((a), (c), and (e)) and the BDF ((b), (c), and (f))

Grahic Jump Location
Fig. 7

Mean value ((a)–(c)), modulus ((d)–(f)), and phase ((g)–(i))) of the first harmonic of the unsteady pressure at h/c = 50% ((a), (d), and (g)), h/c = 75% ((b), (e), and (h)), and h/c = 90% ((c), (f), and (i)) sections

Grahic Jump Location
Fig. 8

Second harmonic of the unsteady pressure distribution along the chord at 90% span

Grahic Jump Location
Fig. 9

Time average of the relative total pressure along the chord over the tip patch (area averaged)

Grahic Jump Location
Fig. 10

Critical aerodynamic damping ratio versus nodal diameter



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