Research Papers: Gas Turbines: Turbomachinery

Development and Applications of a Coupled Particle Deposition—Dynamic Mesh Morphing Approach for the Numerical Simulation of Gas Turbine Flows

[+] Author and Article Information
Peter R. Forsyth

Southwell Laboratory,
Department of Engineering Science,
University of Oxford,
Oxford OX2 0ES, UK
e-mail: peter.forsyth@eng.ox.ac.uk

David R.H. Gillespie

Southwell Laboratory,
Department of Engineering Science,
University of Oxford,
Oxford OX2 0ES, UK
e-mail: david.gillespie@eng.ox.ac.uk

Matthew McGilvray

Southwell Laboratory,
Department of Engineering Science,
University of Oxford,
Oxford OX2 0ES, UK
e-mail: matthew.mcgilvray@eng.ox.ac.uk

1Corresponding author.

Contributed by the Turbomachinery Committee of ASME for publication in the JOURNAL OF ENGINEERING FOR GAS TURBINES AND POWER. Manuscript received July 9, 2017; final manuscript received July 20, 2017; published online October 3, 2017. Editor: David Wisler.

J. Eng. Gas Turbines Power 140(2), 022603 (Oct 03, 2017) (11 pages) Paper No: GTP-17-1325; doi: 10.1115/1.4037825 History: Received July 09, 2017; Revised July 20, 2017

The presence and accretion of airborne particulates, including ash, sand, dust, and other compounds, in gas turbine engines can adversely affect performance and life of components. Engine experience and experimental work have shown that the thickness of accreted layers of these particulates can become large relative to the engine components on which they form. Numerical simulation to date has largely ignored the effects of resultant changes in the passage geometry due to the build-up of deposited particles. This paper will focus on updating the boundaries of the flow volume geometry by integrating the deposited volume of particulates on the solid surface. The technique is implemented using a novel, coupled deposition-dynamic mesh morphing (DMM) approach to the simulation of particulate-laden flows using Reynolds-averaged Navier–Stokes modeling of the bulk fluid, and Lagrangian-based particulate tracking. On an iterative basis, the particle deposition distributions are used to modify the surface topology by altering the locations of surface nodes, which modifies the mesh. The continuous phase solution and particle tracking are then recalculated. The sensitivity to the modeling time steps employed is explored. An impingement geometry case is used to assess the validity of the technique, and a passage with film cooling holes is interrogated. Differences are seen for all sticking and solid phase motion models employed. At small solid particle sizes, considerable disparity is observed between the particle motion modeling approaches, while the position and level of accretion is altered through the use of a nonisotropic stick and bounce model.

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


Chambers, J. C. A. , 1982, “ The 1982 Encounter of British Airways 77 With the Mt. Galuggung Eruption Cloud,” AIAA Paper No. 85-0097.
Casadevall, T. , 1994, “ The 1989-1990 Eruption of Mt. Redoubt Volcano, Alaska: Impacts on Aircraft Operations,” J. Volcanol. Geotherm. Res., 62(1–4), pp. 301–316. [CrossRef]
Bucknell, A. , Gillespie, D. R. H. , and McGilvray, M. , 2017, “ Heat Transfer in the Core Compressor Under Ice Crystal Icing Conditions,” ASME Paper No. GT2017-63077.
Cardwell, N. , Thole, K. , and Burd, S. , 2010, “ Investigation of Sand Blocking Within Impingement and Film-Cooling Holes,” ASME J. Turbomach., 132(2), p. 021020. [CrossRef]
Prenter, R. , Ameri, A. , and Bons, J. , 2016, “ Deposition of a Cooled Nozzle Guide Vane With Nonuniform Inlet Temperatures,” ASME J. Turbomach., 138(10), p. 101005. [CrossRef]
Lundgreen, R. , Sacco, C. , Prenter, R. , and Bons, J. , 2016, “ Temperature Effects on Nozzle Guide Vane Deposition in a New Turbine Cascade Rig,” ASME Paper No. GT2016-57560.
Whitaker, S. , Prenter, R. , and Bons, J. , 2016, “ The Effect of Freestream Turbulence on Deposition for Nozzle Guide Vanes,” ASME J. Turbomach., 137(12), p. 121001. [CrossRef]
Walsh, W. , Thole, K. , and Joe, C. , 2006, “ Effects of Sand Ingestion on the Blockage of Film-Cooling Holes,” ASME Paper No. GT2006-90067.
Wylie, S. , Bucknell, A. , Forsyth, P. , McGilvray, M. , and Gillespie, D. , 2017, “ Reduction in Flow Parameter Resulting From Volcanic Ash Deposition in Engine-Representative Cooling Passages,” ASME J. Turbomach., 139(3), p. 031008. [CrossRef]
Clum, C. , Casaday, B. , Bokar, E. , and Bons, J. , 2014, “ Particle Deposition in Internal Cooling Cavities of a Nozzle Guide Vane—Part I: Experimental Investigation,” ASME Paper No. GT2014-27150.
Wakeman, T. , and Tabakoff, W. , 1982, “ Measured Particle Rebound Characteristics Useful for Erosion Prediction,” ASME Paper No. 82-GT-170.
Tabakoff, W. , Malak, M. , and Hamed, A. , 1987, “ Laser Measurements of Solid-Particle Rebound Parameters Impacting on 2024 Aluminium and 6Al-4V Titanium Alloys,” AIAA J., 25(5), pp. 721–726. [CrossRef]
Bons, J. , Blunt, R. , and Whitaker, S. , 2015, “ A Comparison of Techniques for Particle Rebound Measurement in Gas Turbine Applications,” ASME Paper No. GT2015-43677.
Reagle, C. , Delimont, J. , Ng, W. , and Ekkad, S. , 2014, “ Study of Microparticle Rebound Characteristics Under High Temperature Conditions,” ASME J. Eng. Gas Turbines Power, 136(1), p. 011501. [CrossRef]
Bons, J. , Prenter, R. , and Whitaker, S. , 2016, “ A Simple Physics-Based Model for Particle Rebound and Deposition in Turbomachinery,” ASME Paper No. GT2016-56697.
Singh, S. , and Tafti, D. , 2013, “ Predicting the Coefficient of Restitution for Particle Wall Impacts in Gas Turbine Components,” ASME Paper No. GT2013-95623.
Whitaker, S. , Peterson, B. , Miller, A. , and Bons, J. , 2016, “ The Effect of Particle Loading, Size, and Temperature on Deposition in a Vane Leading Edge Impingement Cooling Geometry,” ASME Paper No. GT2016-57413.
Venturini, P. , Borello, D. H. K. , and Rispoli, F. , 2012, “ Modelling of Particles Deposition in an Environment Relevant to Solid Fuel Boilers,” Appl. Therm. Eng., 49(1), pp. 131–138. [CrossRef]
Gosman, A. , and Ionnides, E. , 1983, “ Aspects of Computer Simulation of Liquid-Fuelled Combustors,” J. Energy, 7(6), pp. 482–490. [CrossRef]
Suman, A. , Kurz, R. , Aldi, N. , Morini, M. , Brun, K. , Pinelli, M. , and Spina, P. , 2014, “ Quantitative CFD Analyses of Particle Deposition on a Subsonic Axial Compressor Blade—Part I: Particle Zones Impact,” ASME J. Turbomach., 137(2), p. 021009. [CrossRef]
El-Batsh, H. , and Haselbacher, H. , 2000, “ Effect of Turbulence Modelling on Particle Dispersion and Deposition on Compressor and Turbine Blade Surfaces,” ASME Paper No. 2000-GT-0519.
Liu, B. Y. H. , and Agarwal, J. K. , 1974, “ Experimental Observation of Aerosol Deposition in Turbulent Flow,” J. Aerosol Sci., 5(2), pp. 145–155.
Dehbi, A. , 2008, “ Turbulent Particle Dispersion in Arbitrary Wall-Bounded Geometries: A Coupled CFD-Langevin-Equation Based Approach,” Int. J. Multiphase Flow, 34(9), pp. 819–828. [CrossRef]
Forsyth, P. , Gillespie, D. , and McGilvray, M. , 2016, “ Validation and Assessment of the Continuous Random Walk Model for Particle Deposition in Gas Turbine Engines,” ASME Paper No. GT2016-57332.
Prenter, R. , Ameri, A. , and Bons, J. , 2016, “ Computational Simulation of Deposition in a Cooled High-Pressure Turbine Stage With Hot Streaks,” ASME Paper No. GT2016-57815.
Tian, L. , and Ahmadi, G. , 2006, “ Particle Deposition in Turbulent Duct Flows—Comparisons of Different Model Predictions,” J. Aerosol Sci., 38(4), pp. 377–397. [CrossRef]
Ghahramani, E. , Abouali, O. , Emad, H. , and Ahmadi, G. , 2014, “ Numerical Analysis of Stochastic Dispersion of Micro-Particles in Turbulent Flows in a Realistic Model of Human Nasal/Upper Airway,” J. Aerosol Sci., 67(1), pp. 188–206. [CrossRef]
Rybalko, M. , Loth, E. , and Lankford, D. , 2012, “ A Lagrangian Particle Random Walk Model for Hybrid RANS/LES Turbulent Flows,” Powder Technol., 221(1), pp. 105–113. [CrossRef]
Chibbaro, S. , and Minier, J.-P. , 2008, “ Langevin PDF Simulation of Particle Deposition in a Turbulent Pipe Flow,” Aerosol Sci., 39(7), pp. 555–571. [CrossRef]
Mehel, A. , Taniere, A. , Oesterle, A. , and Fontaine, J.-R. , 2010, “ The Influence of an Anisotropic Langevin Dispersion Model on the Prediction of Micro-and Nanoparticle Deposition in Wall-Bounded Turbulent Flows,” J. Aerosol Sci., 41(8), pp. 729–744. [CrossRef]
Cao, Y. , Huang, J. , Xu, Z. , and Yin, J. , 2016, “ Insight Into Rime Ice Accretion on an Aircraft Wing and Corresponding Effects on Aerodynamic Performance,” Aeronaut. J., 120(1229), pp. 1101–1122. [CrossRef]
Shen, Z. , Lin, G. , Yu, J. , Bu, Z. , and Du, C. , 2013, “ Three-Dimensional Numerical Simulation of Ice Accretion at the Engine Inlet,” J. Aircr., 50(2), pp. 635–642. [CrossRef]
Casaday, B. , Clum, C. , and Bons, J. , 2014, “ Particle Deposition in Internal Cooling Cavities of a Nozzle Guide Vane—Part II: Analytical and Computational Modelling,” ASME Paper No. GT2014-27155.
ANSYS, 2011, “ Fluent User Defined Function Manual,” ANSYS Inc., Canonsburg, PA.
ANSYS, 2011, “ Fluent 14 Theory Guide,” ANSYS Inc., Canonsburg, PA.
Rosin, P. , and Rammler, E. , 1933, “ The Laws Governing the Fineness of Powdered Coal,” J. Inst. Fuel, 7(1), pp. 29–36.
Schmuecker, M. , 2012, “ Einfluss von Mineralstaeuben auf Keramische Solarabsorber,” DLR-Sonnenkolloquium, Cologne.
Eppelbaum, L. , 2014, “ Lecture Notes in Earth System Sciences,” Applied Geothermics, Springer, Berlin.


Grahic Jump Location
Fig. 1

Comparison of simulations of Liu and Agarwal [22] pipe flow experiments using the CRW and DRW models. Figure taken from Ref. [24].

Grahic Jump Location
Fig. 2

Diagram demonstrating dynamic mesh morphing approach

Grahic Jump Location
Fig. 3

“Boundary layer smoothing” method demonstrating retention of boundary layer mesh despite significant deformation: (a) undeformed mesh, (b) after one mesh morphing iteration, and (c) after four mesh morphing iterations

Grahic Jump Location
Fig. 4

Computational domain for Clum et al. [10]. Top: schematic. Flow top to bottom from inlet, through impingement holes, onto impingement plate, and exhausted front and back from the lower volume. Bottom: Impingement hole geometry. See original paper for further details.

Grahic Jump Location
Fig. 5

Particle size distribution for impingement simulations, as sampled from Ref. [10]

Grahic Jump Location
Fig. 6

Build up of impingement plate surface due to particle deposition. Top: 15 mg injected. Bottom: 30 mg injected. Contours show δ/Dh, height normalized by impingement hole diameter (0.71 mm).

Grahic Jump Location
Fig. 7

Averaged deposition profiles plotted on Clum et al. [10] data. Contours of injected mass; simulation: 25 mg; experiment: 100 mg.

Grahic Jump Location
Fig. 8

Effect of changing solver precision on solution. Left: single-precision solver. Right: double precision solver. Contours indicate equal mass injection; profiles averaged radially.

Grahic Jump Location
Fig. 9

Effect of boundary conditions (BS: bounce-stick model; AS: all-stick model) on the transient nature of build up of deposition. PRO (projected) indicates one deposition iteration (1 mg actual injection) projected out to 25 mg linearly.

Grahic Jump Location
Fig. 10

Effect of previous deposition on current particle rebound characteristics. Deposition statistics in Table 1: (a) particle impacts on (left) initial, undeformed mesh, and (right) deformed mesh. Color scale particle velocity magnitude (m/s). dp = 4 μm, peak displacement δ/Dh  = 0.32, and (b) normal coefficient of restitution for particle sizes considered (solid lines, 10, 7, 4 μm, outside in), and ensemble-averaged particle velocity at impact (dashed lines, 10, 7, 4 μm, mean flow, left to right).

Grahic Jump Location
Fig. 11

Computational domain of Wylie et al. [9] test piece. Flow right to left, exiting through film cooling holes.

Grahic Jump Location
Fig. 12

Particle size distribution as sampled from Wylie et al. [9]

Grahic Jump Location
Fig. 13

Comparison of numerical and experimental deposition in the film cooling holes (top row; flow in the main channel is left to right) and main channel (bottom row). Both looking into their respective flows. Photographs from preliminary work in Ref. [9]: (a) numerical deposition in film cooling hole, (b) experimental deposition in film cooling hole, (c) numerical deposition in main channel, and (d) experimental deposition in main channel.

Grahic Jump Location
Fig. 14

Mesh in film cooling hole undeformed (left) and substantially deformed (right). Flow moved axially right to left.

Grahic Jump Location
Fig. 15

Comparison of deposited mass (points, left) and consequent RFP (lines, right) for CRW and DRW models

Grahic Jump Location
Fig. 16

Comparison of different ash “compositions” and temperatures. Experimental RFP values from Wylie et al. for 1 g injected ash: Chaiten 14–19%, Eyja (Tw = 1193 K): < 0.5%, Eyja (Tw = 1300 K): 6%.



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