0
Research Papers: Gas Turbines: Turbomachinery

Simulating Bluff-Body Flameholders: On the Use of Proper Orthogonal Decomposition for Wake Dynamics Validation

[+] Author and Article Information
Ryan Blanchard

Virginia Polytechnic Institute and State University,
100S Randolph Hall,
Blacksburg, VA 24061
e-mail: rpberlin@vt.edu

Wing Ng

Virginia Polytechnic Institute and State University,
100S Randolph Hall,
Blacksburg, VA 24061
e-mail: wng@vt.edu

Todd K. Lowe

Virginia Polytechnic Institute and State University,
103-C Randolph Hall,
Blacksburg, VA 24061
e-mail: kelowe@vt.edu

Uri Vandsburger

Virginia Polytechnic Institute and State University,
100S Randolph Hall,
Blacksburg, VA 24061
e-mail: uri@vt.edu

Contributed by the Turbomachinery Committee of ASME for publication in the JOURNAL OF ENGINEERING FOR GAS TURBINES AND POWER. Manuscript received February 24, 2014; final manuscript received March 25, 2014; published online July 2, 2014. Editor: David Wisler.

J. Eng. Gas Turbines Power 136(12), 122603 (Jul 02, 2014) (12 pages) Paper No: GTP-14-1119; doi: 10.1115/1.4027556 History: Received February 24, 2014; Revised March 25, 2014

In this article, we describe a novel use of proper orthogonal decomposition (POD) for validation of the structure of dominant flow features in the wake of a bluff-body flameholder. Large-eddy simulations (LES) using both FLUENT and OpenFOAM were conducted based on experiments with planar particle-image velocimetry (PIV) measurements of the same geometry and conditions. With the vision of extending the LES to reacting flows, a validation process is presented that involves a comparison of experimental and computational results, beginning with single-point mean statistics and then extended to the dynamic modes of the data sets as obtained using POD of the instantaneous flow field results. The results exhibit quantitative agreement between both shapes and mode magnitudes for the first POD modes of the measured and simulated data.

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

References

Macquisten, M. A., and Dowling, A. P., 1993, “Low-Frequency Combustion Oscillations in a Model Afterburner,” Combust. Flame, 94(3), pp. 253–264. [CrossRef]
Giacomazzi, E., Battaglia, V., and Bruno, C., 2004, “The Coupling of Turbulence and Chemistry in a Premixed Bluff-Body Flame as Studied by LES,” Combust. Flame, 138(4), pp. 320–335. [CrossRef]
Shanbhogue, S. J., Husain, S., and Lieuwen, T., 2009, “Lean Blowoff of Bluff Body Stabilized Flames: Scaling and Dynamics,” Prog. Energy Combust. Sci., 35(1), pp. 98–120. [CrossRef]
Cocks, P. A. T., Sankaran, V., and Soteriou, M. C., 2013, “Is LES of Reacting Flows Predictive? Part 1: Impact of Numerics,” 51st AIAA Aerospace Sciences Meeting, Grapevine TX, January 7–10, pp. 1–29 [CrossRef].
Giezendanner, R., Keck, O., Weigand, P., Meier, W., Meier, U., Stricker, W., and Aigner, M., 2003, “Periodic Combustion Instabilities in a Swirl Burner Studied by Phase-Locked Planar Laser-Induced Fluorescence,” Combust. Sci. Tech., 175(4), pp. 721–741. [CrossRef]
Roshko, A., 1952, “On the Development of Turbulent Wakes From Vortex Streets,” Ph.D. thesis, California Institute of Technology, Pasadena, CA.
Yang, J.-T., and Tsai, G.-L., 1992, “The Wake Flow Structure of an Open-Slit Vee Gutter,” Experiment. Therm. Fluid Sci., 5(6), pp. 685–696. [CrossRef]
Karniadakis, E. M., and Triantafyllou, G., 1992, “Three-Dimensional Dynamics and Transition to Turbulence in the Wake of Bluff Objects,” J. Fluid Mech., 238, pp. 1–30. [CrossRef]
Spalart, P. R., and Squires, K. D., 1999, “The Status of Detached-Eddy Simulation for Bluff Bodies,” The Aerodynamics of Heavy Vehicles: Trucks, Buses, and Trains, Vol. 19, R. McCallen, F. Browand, and J. Ross, eds., Springer, Berlin, pp. 29–45 [CrossRef].
Shur, M., 2005, “Three Dimensionality in Reynolds-Averaged Navier–Stokes Solutions Around Two-Dimensional Geometries,” AIAA J., 43(6), pp. 1230–1242 [CrossRef].
Achenbach, E., and Heinecke, E., 1981, “On Vortex Shedding From Smooth and Rough Cylinders in the Range of Reynolds Numbers 6×103 to 5×106,” J. Fluid Mech.109, 239–251. [CrossRef]
West, G. S., and Apelt, C. J., 1982, “The Effects of Tunnel Blockage and Aspect Ratio on the Mean Flow Past a Circular Cylinder With Reynolds Numbers Between 104 and 105,” J. Fluid Mech., 114, pp. 361–377. [CrossRef]
Achenbach, E., 1977, “The Effect of Surface Roughness on the Heat Transfer From a Circular Cylinder to the Cross Flow of Air,” Int. J. Heat Mass Transf., 20(4), pp. 359–369. [CrossRef]
Pearson, K.,1901, “On Lines and Planes of Closest Fit to Systems of Points in Space,” Philosoph. Mag., 2(6), pp. 559–572.
Lumley, J. L., 1967, “The Structure of Inhomogeneous Turbulent Flows,” Atmospheric Turbulence and Radio Wave Propagation, A. M. Yaglom, and V. I. Tatarski, eds., Nauka, Moscow, pp. 166–178.
Gergen, I., and Harmanescu, M., 2012, “Application of Principal Component Analysis in the Pollution Assessment With Heavy Metals of Vegetable Food Chain in the Old Mining Areas,” Chemistry Central J., 6(1), p. 156. [CrossRef]
Smith, T. R., Moehlis, J., and Holmes, P., 2005, “Low-Dimensional Modelling of Turbulence Using the Proper Orthogonal Decomposition: A Tutorial,” Nonlinear Dyn., 41(1–3), pp. 275–307. [CrossRef]
Sirovich, L., and Kirby, M., 1987, “Low-Dimensional Procedure for the Characterization of Human Faces,” J. Opt. Soc. Am. A, 4(3), pp. 519–524. [CrossRef] [PubMed]
Hsieh, T.-H., Chen, J.-J. J., Chen, L.-H., Chiang, P.-T., and Lee, H.-Y., 2011, “Time-Course Gait Analysis of Hemiparkinsonian Rats Following 6-Hydroxydopamine Lesion,” Behav. Brain Res., 222(1), pp. 1–9. [CrossRef] [PubMed]
Berkooz, G., Holmes, P., and Lumley, J. L., 1993, “The Proper Orthogonal Decomposition in the Analysis of Turbulent Flows,” Ann. Rev. Fluid Mech., 25, pp. 539–575 [CrossRef].
Fukunaga, K., 1990, Introduction to Statistical Pattern Recognition, 2nd ed., Academic Press, San Diego, CA, p. 592.
Moin, P., and Moser, R. D., 2006, “Characteristic-Eddy Decomposition of Turbulence in a Channel,” J. Fluid Mech., 200, pp. 471–509. [CrossRef]
Duwig, C., and Iudiciani, P., 2009, “Extended Proper Orthogonal Decomposition for Analysis of Unsteady Flames,” Flow, Turbul. Combust., 84(1), pp. 25–47 [CrossRef].
Böhm, B., Brübach, J., Ertem, C., and Dreizler, A., 2008, “Experiments for Combustion-LES Validation,” Flow, Turbul. Combust., 80(4), pp. 507–529 [CrossRef].
Miranda, A. A., Le Borgne, Y.-A., and Bontempi, G., “New Routes From Minimal Approximation Error to Principal Components,” Neural Process. Lett., 27(3), pp. 1–14 [CrossRef].
Maurel, S., Borée, J., and Lumley, J. L., 2001, “Extended Proper Orthogonal Decomposition: Application to Jet/Vortex Interaction,” Flow, Turbul. Combust., 67(2), pp. 125–136 [CrossRef].
Borée, J., 2003, “Extended Proper Orthogonal Decomposition: A Tool to Analyse Correlated Events in Turbulent Flows,” Exp. Fluids, 35(2), pp. 188–192. [CrossRef]
Iudiciani, P., and Duwig, C., 2011, “Large Eddy Simulation of the Sensitivity of Vortex Breakdown and Flame Stabilisation to Axial Forcing,” Flow, Turbul. Combust., 86(3–4), pp. 639–666. [CrossRef]
Smagornisky, J.1963, “General Circulation Experiments With the Primitive Equations,” Mon. Weather Rev., 91(3), pp. 99–164. [CrossRef]
Meyers, J., and Sagaut, P., 2006, “On the Model Coefficients for the Standard and the Variational Multi-Scale Smagorinsky Model,” J. Fluid Mech., 569, pp. 287–319 [CrossRef].
Germano, M., Piomelli, U., Moin, P., and Cabot, W. H., 1991, “A Dynamic Subgrid-Scale Eddy Viscosity Model,” Phys. Fluids A Fluid Dyn., 3(7), pp. 1760–1765. [CrossRef]
Lilly, D. K., 1992, “A Proposed Modification of the Germano Subgrid-Scale Closure Method,” Phys. Fluids A Fluid Dyn., 4(3), pp. 633–635. [CrossRef]
Ranalli, J. A., and Brien, W. F. O., 2009, “Spatially Resolved Analysis of Flame Dynamics for the Prediction of Thermoacoustic Combustion Instabilities,” Virginia Polytechnic Institute and State University, Blacksburg, VA.
Chatterjee, P., Baumann, W. T., Godfrey, A. G., Saunders, W. R., and West, R. L., 2004, “A Computational Fluid Dynamics Investigation of Thermoacoustic Instabilities in Premixed Laminar and Turbulent Combustion Systems,” Ph.D. dissertation, Virginia Polytechnic Institute and State University, Blacksburg, VA.
Chishty, W. A. “Effects of Thermoacoustic Oscillations on Spray Combustion Dynamics with Implications for Lean Direct Injection Systems Effects of Thermoacoustic Oscillations on Spray Combustion Dynamics With Implications for Lean Direct Injection Systems.” Ph.D. dissertation, Virginia Polytechnic Institute and State University, Blacksburg, VA.
Hendricks, A. G., 2003, “Determination of Flame Dynamics for Unsteady Combustion Systems Using Tunable Diode Laser Absorption Spectroscopy,” Master's thesis, Virginia Polytechnic Institute and State University, Blacksburg, VA.
Mattingly, J. D., 2006, Elements of Propulsion: Gas Turbines and Rockets, American Institute of Aeronautics and Astronautics, Reston, VA, pp. 959 [CrossRef].
Kraichnan, R., 1970, “Diffusion by a Random Velocity Field,” Phys. Fluids, 13(1), pp. 22–31 [CrossRef].
Smirnov, R., Shi, S., and Celik, I., 2001, “Random Flow Generation Technique for Large Eddy Simulations and Particle-Dynamics Modeling,” ASME J. Fluids Eng., 123(2), pp. 359–371. [CrossRef]

Figures

Grahic Jump Location
Fig. 1

Large eddy simulation of nonreacting flow around a bluff body from the current work showing isosurfaces of the Q-criterion. The alternating shedding pattern in the bluff body's wake can be seen.

Grahic Jump Location
Fig. 2

Test section temperature and Mach number limitations from facility constraints

Grahic Jump Location
Fig. 3

Schematic of experiment showing the three sections of the rig and the locations of the fuel injection, vee-gutter flameholder, and windows

Grahic Jump Location
Fig. 4

Transparent view of CAD model of the test section showing the locations of the three windows relative to the vee-gutter flameholder and the PIV laser sheet

Grahic Jump Location
Fig. 5

Isometric view of dimensions of the domain used in the LES models based on the width of the flameholder w = 3.05 cm

Grahic Jump Location
Fig. 6

Plan view dimensions of the domain used in LES models, based on the dimension of the vee-gutter flameholder w = 3.05 cm, which has an included angle of 70 deg

Grahic Jump Location
Fig. 7

Illustration of relative sizes and locations of sampling planes relative to the flameholder: 100 × 100 CFD sample grid points (blue) compared to the effective viewing area for the PIV measurements

Grahic Jump Location
Fig. 8

Mean streamwise velocity plots as measured by Pitot probe and as simulated by FLUENT LES

Grahic Jump Location
Fig. 9

v2 Reynolds stress spectra spectrum from LES model. Note the peak due to vortex shedding at Str = 0.24 and also the decay rate that seems to correspond to the –5/3 power law rule.

Grahic Jump Location
Fig. 10

Contours of the sum of the first three POD modes (left: v’, right: u’) obtained using results from PIV (top), FLUENT (middle), and OpenFOAM (bottom). The profiles shown in Fig. 11 are taken along the lines show on each plot.

Grahic Jump Location
Fig. 11

Comparison between numerical and experimental values of the first three POD modes. Left: v’, right: u’.

Grahic Jump Location
Fig. 12

Contours of u-velocity component of the sum of the first three POD modes of PIV (top), FLUENT (middle), and OpenFOAM (bottom). The profiles shown in Fig. 13 are taken along the lines shown imposed on each figure.

Grahic Jump Location
Fig. 13

Plot of u-velocity component of the sum of the first three POD modes taken along the lines indicated in Fig. 12

Grahic Jump Location
Fig. 14

Plot of v-velocity component of the sum of the first three POD as isolated at y/w = 0 at the midplane, near the wall and near the wall but with EPOD to recover the near-wall modeshapes using the eigenvectors of the midplane analysis

Grahic Jump Location
Fig. 15

Comparison between the baseline domain (left) and the domain with the extended domain (right) with increased distance between the inlet and the vee-gutter to produce a bigger/thicker boundary layer in the vicinity of the vee-gutter

Grahic Jump Location
Fig. 16

Comparison between the sum of the first three POD modes of the baseline simulation (top) and thicker boundary layer configuration (bottom). Contours of the U-component (streamwise) of velocity are shown on the left and V-component contours are shown on the right.

Grahic Jump Location
Fig. 17

Comparison of the v-components (vertical) of velocity of the sum of the first three modes of the baseline simulation and the configuration with the thicker boundary layer

Grahic Jump Location
Fig. 18

Comparison between the sum of the first three POD modes of the baseline simulation (top) and periodic side “wall” configuration (bottom). Contours of the U-component (streamwise) of velocity are shown on the left and V-component contours are shown on the right.

Grahic Jump Location
Fig. 19

Comparison of the first v-components (vertical) of velocity of the sum of the first three modes between baseline simulation and the configuration with periodic sides

Grahic Jump Location
Fig. 20

Comparison between the sum of the first three POD modes of the baseline simulation (top) and turbulent inlet configuration (bottom). Contours of the U-component (streamwise) of velocity are shown on the left and V-component contours are shown on the right.

Grahic Jump Location
Fig. 21

Comparison of the v-components (vertical) of velocity of the sum of the first three modes of the baseline simulation and the configuration with the turbulent inlet

Grahic Jump Location
Fig. 22

Velocity contours of first three modes of baseline (top) and refined grid (bottom) simulations, showing contours of u’ (left) and v’ (right)

Grahic Jump Location
Fig. 23

Comparison of the v-components (vertical) of velocity of the sum of the first three modes of the baseline simulation and the refined grid simulation

Grahic Jump Location
Fig. 24

Pope criterion contours at midplane for baseline (top) and refined mesh (bottom) cases

Grahic Jump Location
Fig. 25

Computation time versus number of snapshots N showing that the computational time grows proportional to N3 as expected (solid line is ~N3)

Grahic Jump Location
Fig. 26

Second-order convergence rate of the first five eigenvalues

Grahic Jump Location
Fig. 27

Convergence of first three POD modes with increasing number of samples showing a first-order convergence rate

Tables

Errata

Discussions

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