Validation of a Physics Based Low Order Thermo Acoustic Model of a Liquid Fueled‐‐‐ Gas Turbine‐‐‐ Combustor‐‐‐ and its Application for Predicting‐‐‐ Combustion Driven Oscillations
A dissertation submitted to the Graduate School of the University of Cincinnati
in partial fulfilment of the requirements for the degree of Doctor of Philosophy (Ph.D.)
In the Department of Aerospace Engineering and Engineering Mechanics of the College of Engineering and Applied Sciences
By
Michael Knadler
B.S., Aerospace Engineering, University of Cincinnati, 2011
November 2017
Committee Chair: Dr. Jongguen Lee, Ph.D.
Abstract
This research validates a physics based model for the thermo-acoustic behavior of a liquid-fueled gas turbine combustor as a tool for diagnosing the cause of combustion oscillations. A single nozzle, acoustically tunable gas turbine combustion rig fueled with Jet-A was built capable of operating in the unsteady combustion regime. A parametric study was performed with the experimental rig to determine the operating conditions resulting in thermoacoustic instabilities.
The flame transfer function has been determined for varying fuel injection and flame stabilization arrangements to better understand the feedback loop concerning the heat release and acoustics. The acoustic impedance of the boundaries of the combustion system was experimentally determined. The results were implemented in a COMSOL Multiphysics model as complex impedance boundary conditions at the inlet and exit and a source term to model the flame and heat release. The validity of that model was determined based on an eigenvalue study comparing both the frequency and growth rate of the eigenvalues with the experimentally measured frequencies and pressures of the stable and unstable operating conditions. The model demonstrated that it can accurately predict the instability of the examined operating conditions.
The model also closely predicted the frequency of instability and demonstrated the usefulness of including the experimentally determined acoustic boundary conditions over idealized sound hard boundaries.
i
© Copyright by Michael Knadler 2017 All Rights Reserved
ii
Acknowledgments:
I would first like to thank my advisor, Dr. Jongguen Lee, for his year of guidance and patience in leading me through my research.
I would also like to thank my committee members, Dr. Jay Kim, Dr. Kwanwoo Kim, and Dr. Mark Turner for their insight and advice throughout my dissertation research and studies.
Also in need of recognition are my fellow graduate students who help me both in completing my research, Arda Cakmakci and Thomas Caley, and keeping me sane in and out of the lab, Jun Hee Han.
Of course nothing would get done in the lab without Curt Fox who has helped me with most every piece of technical equipment I have needed.
And finally, Kiev, for always reminding me to never work too hard.
iii
TABLE OF CONTENTS
1. Introduction ...... 1 1.1. The Thermoacoustic Problem ...... 1 1.2. Objectives ...... 1 2. Theory of Thermoacoustic Instability ...... 3 2.1. The Rayleigh Criterion for Thermoacoustic Instability ...... 3 2.2. Mathematical Model for Thermoacoustic Instability ...... 8 2.3. One Dimensional Wave Theory for Acoustic Perturbations ...... 11 3. Previous Approaches to Thermoacoustic Modeling ...... 16 3.1. Computational Fluid Dynamics (CFD) ...... 16 3.2. Thermoacoustic Network Model...... 17 4. Finite Element Modeling for Thermoacoustic Instabilities ...... 26 4.1. Acoustic Wave Coefficients ...... 26 4.1.1. Determining Acoustic Wave Coefficients ...... 26 4.1.2. Acoustic Wave Coefficient Error Analysis ...... 28 4.2. Acoustic Impedance Measurements and Boundary Conditions ...... 29 4.3. Flame Modelling ...... 39 5. Experimental Details ...... 54 5.1. Single Nozzle Acoustically Tunable Gas Turbine Combustion Rig Setup ...... 54 5.1.1. Inline Heater ...... 55 5.1.2. Air Siren ...... 55 5.1.3. Inlet Plenum ...... 57 5.1.4. Fuel Nozzle and Swirler ...... 57 5.1.5. Combustion Chamber and Pressure Vessel ...... 61 5.1.6. Transition Tube ...... 62 5.1.7. Pressure Screws ...... 62 5.2. Dynamic Pressure Sensor Setup ...... 63 5.3. Flame Transfer Function Setup ...... 67 5.4. High Speed Camera Setup ...... 72 5.5. COMSOL Model Development ...... 72 6. Results ...... 78 6.1. Stability Map ...... 78 6.2. Flame Transfer Function Measurement ...... 83
iv
6.3. High Speed Flame Imaging ...... 86 6.3.1. High Speed Color Imaging ...... 87 6.3.2. OH* ICCD Imaging ...... 89 6.4. Time Delay Determination ...... 94 6.5. Acoustic Boundary Condition Measurement ...... 98 6.6. Combustor Temperature Profile ...... 101 6.7. Eigenfrequency Study ...... 102 7. Conclusion ...... 110 8. Bibliography ...... 112 A. APPENDIX A: Flush to Recess Mounting Calibration and MATLAB Code ...... 116 B. APPENDIX B: Flame Transfer Function MATLAB Code ...... 125 C. APPENDIX C: Three-line Pyrometry Theory and Calibration ...... 131 D. APPENDIX D: Impedance Calculation MATLAB Code...... 134
v
LIST OF FIGURES
Figure 2.1: Thermodynamic interpretation of the Rayleigh criterion. Heat addition in phase with pressure (red) and out of phase with pressure (green) ...... 6 Figure 2.2: Block diagram representation of feedback loop between acoustic fluctuations and heat release ...... 8 Figure 2.3: Acoustic wave propagation in a duct ...... 13 Figure 3.1: (a) Schematic of a model gas turbine combustor with acoustic waves A and B upstream and downstream of the flame (b) Block diagram of thermoacoustic network model ... 18 Figure 3.2: (left) Upstream and downstream acoustic variables with relation to thermoacoustic network element (right) Reimann invariants upstream and downstream with respective directions of flow ...... 19 Figure 3.3: Converging-diverging nozzle ...... 20 Figure 3.4: Combustor with premix duct used to develop transfer matrices presented in Equations 3.2.22 and 3.2.23 ...... 23 Figure 4.1: Flow through an orifice showing vena contracta effect ...... 32 Figure 4.2: Orifice showing vena contracta producing edge vortices ...... 32 Figure 4.3: Effect of bias flow on acoustic impedance as determined by Jing and Sun [29] ...... 34 Figure 4.4: Su's [38] impedance model and its dependence on frequency ...... 34 Figure 4.5: Testud's [39] model for predicting orifice whistling. Negative indicates whistling potential...... 36 Figure 4.6: RANS simulation confirming orifice whistling potential model ...... 36 Figure 4.7: Converging-diverging nozzle with shock used to model choked orifice impedance . 37 Figure 4.8: Magnitude and phase of reflection coefficient comparison between numerical (circle) and analytical (solid) results of Stow [43] ...... 38 Figure 4.9: Impedance comparison between finite difference solver, Euler solver, and analytical results with a compaact nozzle assumption [46] ...... 39 Figure 4.10: Schematic of combustor used to derive T22 flame model ...... 41 Figure 4.11: Comparison of measured (red crosses) and modeled (blue line) values for transfer matrix element T22 for various temperatures [48] ...... 45 Figure 4.12: (left) Absolute value and (right) phase of flame transfer matrix elements comparison for measured (black square), Rankine-Hugoniot approximation (red line), and OH* chemiluminescence measurement (blue diamonds) results ...... 47 Figure 4.13: (left) Effect of β on first resonant frequency of combustor comparing exact solution (solid line) and one term Galerkan approximation (dashed line), (right) effect of time delay on growth rate for varying β [13] ...... 50 Figure 4.14: Effect of using temperature jump on first resonant frequency for single temperature jump (solid line), uniform temperature rise (dashed line), and five smaller jump approximations (-.-) ...... 51 Figure 5.1: Overview of single nozzle, acoustically tunable, gas turbine combustion rig ...... 54 Figure 5.2: Cross sectional view of combustion rig ...... 55 Figure 5.3: Inlet section connected to air modulation siren and siren bypass line ...... 56 Figure 5.4: Air modulation siren ...... 56 Figure 5.5: 3D printed and cross sectional CAD view of swirler used in testing ...... 58 Figure 5.6: Swirler mounted on effusion plate viewed from upstream/inlet side (left) and downstream/combustor side (right). Splash plate not included...... 60
vi
Figure 5.7: Downstream/Combustor side view of swirler mounted with splash plate attached and flush with flared exit of swirler ...... 60 Figure 5.8: Pressure drop vs. main air flow rate for open and closed effusion hole cases ...... 61 Figure 5.9: Piezoelectric pressure transducer seated in wave guide mount for recessed mounting ...... 63 Figure 5.10: Inlet orifice impedance test setup ...... 64 Figure 5.11: Schematic of exit impedance setup ...... 66 Figure 5.12: Exit orifice impedance testing setup ...... 66 Figure 5.13: Microphone mounting on transition tube for exit orifice impedance testing ...... 67 Figure 5.14: Optical arrangement of the three line pyrometer...... 69 Figure 5.15: Time traces of downstream pressure,, three line pyrometer emissions, and OH* chemiluminescence emissions showing in-phase relationship ...... 70 Figure 5.16: (left) Schematic of OH* chemiluminescence setup focusing emissions from the side window of the combustor onto the PMT sensor and (right) actual experimental hardware in place for testing ...... 71 Figure 5.17: COMSOL geometry with elements representing (from left to right) inlet, nozzle, combustor, flame zone, and transition tube ...... 73 Figure 5.18: COMSOL geometry for the 0.0508m-long transition tube case with fine mesh visible ...... 74 Figure 6.1: Acoustic pressure fluctuations in the combustor for the 0.254m-long transition tube setup ...... 79 Figure 6.2: Acoustic pressure fluctuations in the combustor for the 0.508m-long transition tube setup ...... 79 Figure 6.3: Combustor pressure fluctuation spectrum for Case 1 ...... 81 Figure 6.4: Combustor pressure fluctuation spectrum for Case 2 ...... 82 Figure 6.5: Combustor pressure fluctuation spectrum for Case 3 ...... 82 Figure 6.6: Combustor pressure fluctuation spectrum for Case 4 ...... 83 Figure 6.7: Exemplary plots used in determining the gain of the FTF for (top) 1500SLPM case fordced at 400Hz and (bottom) 2000SLPM case forced at 400Hz ...... 84 Figure 6.8: Gain of the FTF for 2000SLPM and 1500SLPM cases ...... 85 Figure 6.9: Phase of the FTF for 2000SLPM and 1500SLPM cases ...... 86 Figure 6.10: Averaged high speed color camera images for (top left) Case 1 (top right) Case 2 and (bottom) Case 4 showing pressure fluctuations in the inlet (P1') and combustor (Pcomb'), frequency of fluctuations ...... 88 Figure 6.11: Averaged ICCD images for (top left) Case 1 (top right) Case 2 and (bottom) Case 4 with locations of peak and centroid intensity ...... 90 Figure 6.12: Average ICCD images of 2000SLPM 20in case under forcing from 200Hz, 400Hz, 600Hz, and 800Hz ...... 91 Figure 6.13: Average ICCD images of 1500SLPM 20in case under forcing from 200Hz, 400Hz, 600Hz, and 800Hz ...... 91 Figure 6.14: Phase averaged ICCD images for the 2000SLPM 20in case under natural instability ...... 93 Figure 6.15: Intensity contours for the phase averaged ICCD images for the 2000SLPM 20in case under natural instability ...... 93 Figure 6.16: Movement of centroid location downstream of nozzle for the 2000SLPM 20in case under 200Hz forcing by the siren...... 94
vii
Figure 6.17: Video FFT analyzed OH* emission peaks for (top left) 200Hz, (top right) 400Hz, (bottom left) 600Hz, and (bottom right) 800Hz ...... 97 Figure 6.18: Real component to impedance for upstream and downstream orifices for both 1500SLPM and 2000SLPM cases ...... 99 Figure 6.19: Imaginary component to impedance for upstream and downstream orifices for both 1500SLPM and 2000SLPM cases ...... 99 Figure 6.20: Reflection coefficient for upstream and downstream orifices for both 1500SLPM and 2000SLPM cases ...... 100 Figure 6.21: Percent error in acoustic pressure amplitude predictions for upstream and downstream orifices for both 1500SLPM and 2000SLPM cases ...... 100 Figure 6.22: Temperature profiles used in COMSOL model ...... 102 Figure 6.23: Mode shape for Case 1 comparing predicted COMSOL values to experimentally determined acoustic pressure amplitudes in inlet and combustor with error bars indicating one standard deviation ...... 108 Figure 6.24: Mode shape for Case 2 comparing predicted COMSOL values to experimentally determined acoustic pressure amplitudes in inlet and combustor with error bars indicating one standard deviation ...... 109
viii
LIST OF TABLES
Table 5.1: Swirler Geometry ...... 58 Table 5.2: Measured Effective Area and Projected Mass Flow Rates for Ten SN0.6 Swirlers [54] ...... 59 Table 5.3: Inlet transducer mount locations ...... 65 Table 5.4: Exit orifice impedance test transducer locations ...... 66 Table 5.6: Mesh density impact on eigenfrequency of rig model without flow or source present74 Table 5.7: COMSOL grid impact study results with experimental operating conditions in place for 2000SLPM 0.508m-long transition tube case ...... 75 Table 5.8: Exemplary input parameters for the COMSOL model ...... 76 Table 6.1: Equivalence ratios for stability map operating conditions ...... 78 Table 6.2: Time delay calculations based on known convection distances and calculated velocities ...... 96 Table 6.3: Results of final COMSOL model compared to experimentally determined acoustic pressure fluctuations in the combustor and their respective frequencies ...... 102 Table 6.4: Effect of impedance boundary conditions on converged eigenfrequencies in the COMSOL model ...... 103 Table 6.5: Effect of temperature gradient on converged eigenfrequencies in COMSOL model 104 Table 6.6: Model prediction with additional time delays ...... 105 Table 6.7: Eigenfrequency predictions using visible spectrum emissions in calculating flame zone and time delay ...... 107
ix
1. Introduction
1.1. The Thermoacoustic Problem
In a drive to produce more efficient and lower NOx emission gas turbine engines techniques have been employed in engine design that increase the system’s risk for self-sustained combustion instabilities. These included the use of lean premixed prevaporized (LPP) combustion that along with reducing local hot spots in the combustion chamber, leading to lower
NOx emissions, reduce both flame anchoring and dampening [1]. Several consequences of LPP combustion lead to this increased susceptibility for thermoacoustic oscillations including a lower dilution air supply, leaner equivalence ratio, and a pressure antinode located near the flame [2].
These oscillations, under the right conditions, can cause undesirable pressure fluctuations in the combustor of such amplitude as to significantly increase the wear and tear on engine components leading to unsafe operating conditions and failure. In addition, thermoacoustic instabilities are not just a local phenomenon. Acoustics of the entire combustion system, from inlet boundary conditions through outlet boundary conditions and everything in between through which the acoustic waves propagate, can affect the combustor’s stability [3]. Although local flow characteristics can play an important role, a global view of the system must be taken to properly diagnose thermoacoustic instabilities. This often inhibits the further development of stable low- emission combustors and so a greater understanding of the mechanisms which initiate and sustain the instabilities and various operating conditions is needed.
1.2. Objectives
The objective of this study is to develop and validate a theoretical finite element model capable of predicting the formation of combustion driven pressure oscillations through a detailed understanding of the physical response of the flame to periodic disturbances, the acoustic
1 conditions at the boundaries, and the propagation of acoustic waves in the air supply and combustor. This knowledge is currently limited and the validation of these combustion dynamics models is not found.
2
2. Theory of Thermoacoustic Instability
2.1. The Rayleigh Criterion for Thermoacoustic Instability
Unstable combustion is the self-sustained combustion oscillations at the resonant frequency of the combustion chamber which are the result of closed loop coupling between unsteady heat release and pressure fluctuations. These flame dynamics have been observed and studies by interested researchers for well over the past two centuries. The first recorded observation was of what became known as the “singing” flame by Higgins in 1777 [4]. Higgins and other researchers discovered that anchoring a flame on a smaller diameter fuel line inside a larger diameter tube would produce an audible level of sound by exciting the fundamental or harmonic modes of the larger tube, similar to the excitation of an organ pipe. Years later a “dancing” flame was observed by LeConte where a flame would pulse in sync to beats of music (if only the flame could learn to act then it would a true triple threat!). Around the same time, Rijke discovered that sound could be generated in a vertical tube, open on both ends, by heating a metal gauze placed exactly at one quarter the distance from the bottom to top of the tube. This setup, known as the
Rijke tube, comes about as a result of acoustic velocity fluctuations in the bottom and top halves of the tube being opposite each other and the heat source being placed at the L/4 location causing the velocity fluctuations to lead (in the bottom half) or lag (in the top half) by 90 o
These phenomena were first described by Lord Rayleigh in 1878 [5]. In his own words the thermoacoustic phenomenon is described as follows: “If heat be periodically communicated to, and abstracted from, a mass of air vibrating in a cylinder bounded by a piston, the effect produced will depend upon the phase of the vibration at which the transfer of heat takes place. If heat be given to the air at the moment of greatest condensation or to be taken from it at the moment of greatest rarefaction, the vibration is encouraged. On the other hand, if heat be given
3 at the moment of greatest rarefaction, or extracted at the moment of greatest condensation, the vibration is discouraged”. In practical term applied to gas turbines, the acoustic waves cause fluctuations in heat release by altering the fuel and/or air mass flow rate. The fluctuations in resulting heat release then amplify the acoustic waves which results in self-sustained instabilities.
When these heat release fluctuations and pressure fluctuations are in phase the system will be naturally unstable and when they are out of phase natural dampening will occur. More specifically, any upstream flow perturbation can perturb the heat release further downstream at the flame. When the heat release fluctuates, then volume expansion of the flame fluctuates as well. As this cycle continues the volume of the expansion fluctuations of the flame will release acoustic waves throughout the combustion chambers at a frequency dependent upon the initial perturbation frequency. The acoustic pressure wave will propagate all along the combustor, eventually reflecting back upstream and causing acoustic velocity fluctuations in the flame.
This theory can be formulated into an inequality integrating the relationship between pressure and heat release fluctuations and the natural dissipation of the system.