Research Papers

Uncertainty Quantification: A Stochastic Method for Heat Transfer Prediction Using LES

[+] Author and Article Information
M. Carnevale

Dipartimento di Energetica,
Università degli Studi di Firenze,
via di S. Marta, 3,
Firenze 50139, Italia;
Whittle Laboratory,
1 JJ Thompson Evenue,
University of Cambridge,
Cambridge CB3 0DY, UK

F. Montomoli

Whittle Laboratory,
1 JJ Thompson Evenue,
University of Cambridge,
Cambridge CB3 0DY, UK;
Division of Mechanical, Medical, and Aerospace Engineering,
University of Surrey,
Guildford GU2 7XH, UK

A. D’Ammaro

Whittle Laboratory,
1 JJ Thompson Evenue,
University of Cambridge,
Cambridge CB3 0DY, UK

F. Martelli

Dipartimento di Energetica,
Università degli Studi di Firenze,
via di S. Marta, 3,
Firenze 50139, Italia

Contributed by the International Gas Turbine Institute (IGTI) of ASME for publication in the JOURNAL OF TURBOMACHINERY. Manuscript received August 2, 2012; final manuscript received August 31, 2012; published online June 28, 2013. Editor: David Wisler.

J. Turbomach 135(5), 051021 (Jun 28, 2013) (8 pages) Paper No: TURBO-12-1162; doi: 10.1115/1.4007836 History: Received August 02, 2012; Revised August 31, 2012

In computational fluid dynamics (CFD), it is possible to identify namely two uncertainties: epistemic, related to the turbulence model, and aleatoric, representing the random-unknown conditions such as the boundary values and or geometrical variations. In the field of epistemic uncertainty, large eddy simulation (LES and DES) is the state of the art in terms of turbulence closures to predict the heat transfer in internal channels. The problem is still unresolved for the stochastic variations and how to include these effects in the LES studies. In this paper, for the first time in literature, a stochastic approach is proposed to include these variations in LES. By using a classical uncertainty quantification approach, the probabilistic collocation method is coupled to numerical large eddy simulation (NLES) in a duct with pin fins. The Reynolds number has been chosen as a stochastic variable with a normal distribution. The Reynolds number is representative of the uncertainties associated with the operating conditions, i.e., velocity and density, and geometrical variations such as the pin fin diameter. This work shows that assuming a Gaussian distribution for the Reynolds number of ±25%, it is possible to define the probability to achieve a specified heat loading under stochastic conditions, which can affect the component life by more than 100%. The same method, applied to a steady RANS, generates a different level of uncertainty. New methods have been proposed based on the different level of aleatoric uncertainties which provides information on the epistemic uncertainty. This proves, for the first time, that the uncertainties related to the unknown conditions, aleatoric, and those related to the physical model, epistemic, are strongly interconnected. This result, which is idealized for this specific issue, can be extrapolated, and has direct consequences in uncertainty quantification science and not only in the gas turbine world.

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


Kassinos, S. C., Alnger, A. C., Iaccarino, G., and Moin, P., 2007, Complex Effects in Large Eddy Simulation, Springer, New York.
Rozati, A., and Tafti, D. K., 2008, “Large-Eddy Simulations of Leading Edge Film Cooling: Analysis of Flow Structures, Effectiveness, and Heat Transfer Coefficient,” Int. J. Heat Fluid Flow, 29(1), pp. 1–17. [CrossRef]
Tafti, D. K., Zhang, X., Huang, W., and Wang, G., 2000, “Large-Eddy Simulations of Flow and Heat Transfer in Complex Three-Dimensional Multilouvered Fins,” ACRC Tech. Report Air Conditioning and Refrigeration Center, Report No. TR-165.
Iaccarino, G., 2008, Quantification of Uncertainty in Flow Simulations Using Probabilistic Methods, (VKI Lectures Series), von Karman Institute, Rhode-St-Genèse, Belgium.
Walters, R. W., and Huyse, L., 2002, “Uncertainty Analysis for Fluid Mechanics With Application,” NASA/CR Report No. 2002-211449.
Montomoli, F., Massini, M., and Salvadori, S., 2010, “Geometrical Uncertainty in Turbomachinery: Tip Gap and Fillet Radius,” Comput. Fluids, 46(1), pp. 362–368. [CrossRef]
Montomoli, F., Massini, M., Salvadori, S., and Martelli, F., 2012, “Geometrical Uncertainty and Film Cooling: Fillet Radii,” ASME J. Turbomach.,134(1), p. 011019. [CrossRef]
D’Ammaro, A., and Montomoli, F., 2011, “Geometrical Uncertainty and Film Cooling,” Comput. Fluids, 71, pp. 320–326. [CrossRef]
Tucker, P. G., 2011, “Computation of Unsteady Turbomachinery Flows: Part 2—LES and Hybrids,” Prog. Aerosp. Sci., 47(7), pp. 546–569. [CrossRef]
Tatang, M. A., 1995, “Direct Incorporation of Uncertainty in Chemical and Environmental Engineering System,” Ph.D. thesis, MIT, Cambridge, MA.
Han, J., Dutta, S., and Ekkad, S. V., 2000, Gas Turbine Heat Transfer and Cooling Technology, Taylor & Francis, New York.
Ames, F. E., and Dvorak, L. A., 2006, “Turbulent Transport in Pin Fin Arrays: Experimental Data and Predictions,” ASME J. Turbomach., 128(1), pp. 71–81. [CrossRef]
Ames, F. E., Nordquist, C. A., and Klennert, L. A., 2007, “Endwall Heat Transfer Measurements in a Staggered Pin Fin Array With an Adiabatic Pin,” Proceedings of the GT 2007 ASME Turbo Expo, Montreal, Canada, May 14–17, ASME Paper No. GT2007-27432. [CrossRef]
Delibra, G., Borello, D., Hanjalic, K., and Rispoli, F., 2010, “An LES Insight Into Convective Mechanism of Heat Transfer in a Wall-Bounded Pin Matrix,” Proceedings of the International Heat Transfer Conference (IHTC14), Washington, DC, Aug. 8–13, ASME Paper No. IHTC14-23205. [CrossRef]
Delibra, G., Hanjalic, K., Borello, D., and RispoliF., 2009, “Vortex Structures and Heat Transfer in a Wall-Bounded Pin Matrix: LES With a RANS Wall-Treatment,” Int. J. Heat Fluid Flow, 31(5), pp. 740–753. [CrossRef]
Delibra, G., Borello, D., Hanjalic, K., and Rispoli, F., 2009, “URANS of Flow and Endwall Heat Transfer in a Pinned Passage Relevant to Gas-Turbine Blade Cooling,” Int. J. Heat Fluid Flow, 30(3), pp. 549–560. [CrossRef]
Simoneau, R. J., and Van Fossen, G. J., Jr., 1984, “Effect of Location in an Array on Heat Transfer to a Short Cylinder in Crossflow,” ASME J. Heat Transfer, 106(1), pp. 183–190. [CrossRef]
Adami, P., Martelli, F., and Michelassi, V., 2000, “Three Dimensional Investigations for Axial Turbines by an Implicit Unstructured Multi-Block Flow Solver,” Proceedings of the ASME Turbo Expo 2000, Munich, Germany, ASME Paper No. GT2000-0636.
Montomoli, F., Adami, P., and Martelli, F., 2009, “A Finite-Volume Method for the Conjugate Heat Transfer in Film Cooling Devices,” J. Power Energy, 223(A2), pp. 191–200. [CrossRef]
Salvadori, S., Adami, P., and Martelli, F., 2011, “On the Implementation of a Phase Lag Approach for Multi-Row Simulations,” Proceedings of the 10th International Symposium on Experimental Computational Aerothermodynamics of Internal Flows, Brussels, Belgium, July 4–7, Paper No. ISAIF10-063.
Salvadori, S., Bernardini, C., Adami, P., and Martelli, F., 2009, “Turbulence and Transition Modelling in Transonic Turbine Stage,” Proceedings of the 19th ISABE Conference, Montreal, Canada, September 7–11.
Bernardini, C., Carnevale, M., Salvadori, S., and Martelli, F., 2011, “On the Assessment of an Unstructured Finite-Volume DES/LES Solver for Turbomachinery Application,” WSEAS Trans. Fluid Mech., 6(3), pp. 160–173.
Roe, P. L., 1986, “Characteristic-Based Scheme for the Euler Equations,” Ann. Rev. Fluid Mech., 18, pp. 337–365. [CrossRef]
Anderson, W. K., and Bonhaus, D. L., 1994, “An Implicit Upwind Algorithm for Computing Turbulent Flows on Unstructured Grids,” Comput. Fluids, 23(1), pp. 1–21. [CrossRef]
Wilcox, D. C., 1993, Turbulence Modelling for CFD, DCW Industries Inc., La Canada, CA.
Menter, F. R., 1992, “Improved Two-Equation k-ω Turbulence Models for Aerodynamic Flow,” NASA Technical Memorandum 1992, Technical Report No. 103975.
Dahlquist, G., and Bjork, A., 1974, Numerical Methods, Prentice-Hall, Englewood Cliffs, NJ.
Calvo, M., Franco, J. M., and Randez, L. A., 2004, “A New Minimum Storage Runge-Kutta Scheme for Computational Acoustics,” J. Comput. Phys., 201(1), pp. 1–12. [CrossRef]
Thornber, B., Mosedale, A., Drikakis, D., Youngs, D., and Williams, R. J. R., 2008, “An Improved Reconstruction Method for Compressible Flows With Low Mach Number Features,” J. Comput. Phys., 227, pp. 4873–4894. [CrossRef]
Mary, I., and Sagaut, P., 2002, “Large Eddy Simulation of Flow Around an Airfoil Near Stall,” AIAA J., 40(6), pp. 1139–1145. [CrossRef]
Garnier, E., Mossi, M., Sagaut, P., and Davile, M., 1999, “On the Use of Shock-Capturing Schemes for Large-Eddy Simulation,” J. Comput. Phys., 153(2), pp. 272–311. [CrossRef]
Belardini, E., Adami, P., and Martelli, F., 2001, Development of an Unsteady Parallel Approach for 3D Stator-Rotor Interaction, Science and Supercomputing at CINECA, Casalecchio di Reno, Bologna, Italy,, pp. 338–355.
Ames, F. E., Dvorak, L. A., and Morrow, M. J., 2005, “Turbulent Argumentation of Internal Convection Over Pins in Staggered Pin-Fin Arrays,” ASME J. Turbomach., 127, pp. 183–190. [CrossRef]
Xiu, D., 2010, Numerical Methods for Stochastic Computations: A Spectral Methods Approach, Princeton University Press, Princeton, NJ.
Churchill, S. W., and Usagi, R., 1972, “A General Expression for the Correlation of Rates of Transfer and Other Phenomena,” AIChE J., 18(6), pp. 1121–1128. [CrossRef]


Grahic Jump Location
Fig. 2

Mesh details: left—overall view; right—pin mesh

Grahic Jump Location
Fig. 3

Pressure coefficient pin 1, pin 2, pin 3, comparison CFD with experimental data

Grahic Jump Location
Fig. 4

Comparison of mean velocity U/Vmax in the midplane along the line (a–a) for rows 1–3

Grahic Jump Location
Fig. 5

Spectrum of velocity signal

Grahic Jump Location
Fig. 6

Nusselt average distribution

Grahic Jump Location
Fig. 7

Statistical distribution of Nusselt number and σ: NLES calculation

Grahic Jump Location
Fig. 8

Statistical distribution of Nusselt number and σ: RANS calculation

Grahic Jump Location
Fig. 10

φ-parameter ration distribution in streamwise direction for RANS closure

Grahic Jump Location
Fig. 9

φ-parameter ration distribution in streamwise direction for NLES closure




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