Amplitude-Dependent Damping and Driving Rates of High-Frequency Thermoacoustic Oscillations in a Lab-Scale Lean-Premixed Gas Turbine Combustor

2021 ◽  
Author(s):  
Thomas Hofmeister ◽  
Thomas Sattelmayer
Keyword(s):  
Flow Field ◽  
Gas Turbine ◽  
Limit Cycle ◽  
High Frequency ◽  
Mean Flow ◽  
Gas Turbine Combustor ◽  
Cfd Simulations ◽  
The Mean ◽  
Amplitude Dependency ◽  
Thermoacoustic Oscillations

Abstract This paper presents numerical investigations of the amplitude-dependent stability behavior of thermoacoustic oscillations at screech level frequencies in a lean-premixed, swirl-stabilized, lab-scale gas turbine combustor. A hybrid Computational Fluid Dynamics / Computational AeroAcoustics (CFD / CAA) approach is applied to individually compute thermoacoustic damping and driving rates for various acoustic amplitude levels at the combustors' first transversal (T1) eigenfrequency. Forced CFD simulations with the Unsteady Reynolds-Averaged Navier-Stokes (URANS) equations mimic the real combustor's rotating T1 eigenmode. An increase of the forcing amplitude over time allows observation of the amplitude-dependent flow field and flame evolution. In accordance with measured OH*-chemiluminescence images, a pulsation amplitude-dependent flame contraction is reproduced in the CFD simulations. At several amplitude levels, period-averaged flow fields are then denoted as reference states, which serve as inputs for the CAA part. There, eigenfrequency simulations with linearized flow equations are performed with the Finite Element Method (FEM). The outcomes are damping and driving rates as a response to the amplitude-dependency of the mean flow field. It is found that driving due to flame-acoustics interactions governs a weak amplitude-dependency, which agrees with experimentally based studies at the authors' institute. This disqualifies the perception of heat release saturation as the root-cause for limit-cycle oscillations in this high-frequency thermoacoustic system. Instead, significantly increased dissipation due to the interaction of acoustically induced vorticity perturbations with the mean flow is identified, which may explain the formation of a limit-cycle.

Amplitude-Dependent Damping and Driving Rates of High-Frequency Thermoacoustic Oscillations in a Lab-Scale Lean-Premixed Gas Turbine Combustor

2021 ◽  
Author(s):  
Thomas Hofmeister ◽  
Thomas Sattelmayer
Keyword(s):  
Flow Field ◽  
Gas Turbine ◽  
High Frequency ◽  
Mean Flow ◽  
Gas Turbine Combustor ◽  
Cfd Simulations ◽  
Root Cause ◽  
The Mean ◽  
Amplitude Dependency ◽  
Thermoacoustic Oscillations

Abstract This paper presents the numerical investigations of amplitude-dependent stability behavior of thermoacoustic oscillations at screech level frequencies in a lean-premixed, atmospheric, swirl-stabilized, lab-scale gas turbine combustor. A hybrid Computational Fluid Dynamics / Computational AeroAcoustics (CFD / CAA) approach is applied to individually compute thermoacoustic damping and driving rates for various acoustic amplitude levels at the combustors’ first transversal (T1) eigenfrequency. Harmonically forced CFD simulations with the Unsteady Reynolds-Averaged Navier-Stokes (URANS) equations mimic the real combustor’s rotating T1 eigenmode. A slow and monotonous increase of the forcing amplitude over time allows observation of the amplitude-dependent flow field and flame evolution. In accordance with measured OH*-chemiluminescence images, a pulsation amplitude-dependent flame contraction is reproduced in the CFD simulations, where acoustically induced backflow at the combustion chamber inlet is identified as the root-cause of this phenomenon. At several amplitude levels, period-averaged flow fields are then denoted as reference states, which serve as inputs for the CAA part. There, eigenfrequency simulations with linearized flow equations are performed with the Finite Element Method (FEM). The outcomes are damping and driving rates as a response to the amplitude dependency of the mean flow field, which combined give the net thermoacoustic growth rate. It is found that driving due to flame-acoustics interactions only governs a weak amplitude dependency, which agrees with prior, experimentally based studies at the authors’ institute. This disqualifies the perception of heat release saturation as the root-cause for limit-cycle oscillations — at least in this high-frequency thermoacoustic system. Instead, significantly increased dissipation due to the interaction of acoustically induced vorticity perturbations with the mean flow is identified, which may explain the formation of a limit-cycle.

Pulsation-Amplitude-Dependent Flame Dynamics of High-Frequency Thermoacoustic Oscillations in Lean-Premixed Gas Turbine Combustors

2017 ◽  
Author(s):  
Frederik M. Berger ◽  
Tobias Hummel ◽  
Bruno Schuermans ◽  
Thomas Sattelmayer
Keyword(s):  
Heat Release ◽  
Gas Turbine ◽  
Growth Rates ◽  
High Frequency ◽  
Individual Growth ◽  
Mean Flow ◽  
Flame Dynamics ◽  
Pulsation Amplitude ◽  
The Mean ◽  
Thermoacoustic Oscillations

This paper presents the experimental investigation of pulsation-amplitude-dependent flame dynamics associated with transverse thermoacoustic oscillations at screech level frequencies in a generic gas turbine combustor. Specifically, the flame behavior at different levels of pulsation amplitudes is assessed and interpreted. Spatial dynamics of the flame are measured by imaging the OH* chemiluminescence signal synchronously to the dynamic pressure at the combustor’s face plate. First, linear thermoacoustic stability states, modal dynamics, as well as flame-acoustic phase relations are evaluated. It is found that the unstable acoustic modes converge into a predominantly rotating character in the direction of the mean flow swirl. Furthermore, the flame modulation is observed to be in phase with the acoustic pressure at all levels of the oscillation amplitude. Second, distributed flame dynamics are investigated by means of visualizing the mean and oscillating heat release distribution at different pulsation amplitudes. The observed flame dynamics are then compared against numerical evaluations of the respective amplitude-dependent thermoacoustic growth rates, which are computed using analytical models in the fashion of a non-compact flame-describing function. While results show a nonlinear contribution for the individual growth rates, the superposition of flame deformation and displacements balances out to a constant flame driving. This latter observation contradicts the state-of-the-art perception of root-causes for limit-cycle oscillations in thermoacoustic gas turbine systems, for which the heat release saturates with increasing amplitudes. Consequently, the systematic observations and analysis of amplitude-dependent flame modulation shows alternative paths to the explanation of mechanisms that might cause thermoacoustic saturation in high frequency systems.

Pulsation-Amplitude-Dependent Flame Dynamics of High-Frequency Thermoacoustic Oscillations in Lean-Premixed Gas Turbine Combustors

2017 ◽  
Vol 140 (4) ◽  
Author(s):  
Frederik M. Berger ◽  
Tobias Hummel ◽  
Bruno Schuermans ◽  
Thomas Sattelmayer
Keyword(s):  
Heat Release ◽  
Gas Turbine ◽  
Growth Rates ◽  
High Frequency ◽  
Individual Growth ◽  
Mean Flow ◽  
Flame Dynamics ◽  
Pulsation Amplitude ◽  
The Mean ◽  
Thermoacoustic Oscillations

This paper presents the experimental investigation of pulsation-amplitude-dependent flame dynamics associated with transverse thermoacoustic oscillations at screech level frequencies in a generic gas turbine combustor. Specifically, the flame behavior at different levels of pulsation amplitudes is assessed and interpreted. Spatial dynamics of the flame are measured by imaging the OH⋆ chemiluminescence (CL) signal synchronously to the dynamic pressure at the combustor's face plate. First, linear thermoacoustic stability states, modal dynamics, and flame-acoustic phase relations are evaluated. It is found that the unstable acoustic modes converge into a predominantly rotating character in the direction of the mean flow swirl. Furthermore, the flame modulation is observed to be in phase with the acoustic pressure at all levels of the oscillation amplitude. Second, distributed flame dynamics are investigated by means of visualizing the mean and oscillating heat release distribution at different pulsation amplitudes. The observed flame dynamics are then compared against numerical evaluations of the respective amplitude-dependent thermoacoustic growth rates, which are computed using analytical models in the fashion of a noncompact flame-describing function. While results show a nonlinear contribution for the individual growth rates, the superposition of flame deformation and displacement balances out to a constant flame driving. This latter observation contradicts the state-of-the-art perception of root-causes for limit-cycle oscillations in thermoacoustic gas turbine systems, for which the heat release saturates with increasing amplitudes. Consequently, the systematic observations and analysis of amplitude-dependent flame modulation shows alternative paths to the explanation of mechanisms that might cause thermoacoustic saturation in high frequency systems.

A Time-Domain Network Model for Nonlinear Thermoacoustic Oscillations

2008 ◽  
Author(s):  
Simon R. Stow ◽  
Ann P. Dowling
Keyword(s):  
Heat Release ◽  
Frequency Domain ◽  
Limit Cycle ◽  
Network Model ◽  
Time Domain ◽  
Acoustic Waves ◽  
Mean Flow ◽  
The Mean ◽  
The Time Domain ◽  
Thermoacoustic Oscillations

Lean premixed prevaporised (LPP) combustion can reduce NOx emissions from gas turbines, but often leads to combustion instability. Acoustic waves produce fluctuations in heat release, for instance by perturbing the fuel-air ratio. These heat fluctuations will in turn generate more acoustic waves and in some situations linear oscillations grow into large amplitude self-sustained oscillations. The resulting limit cycles can cause structural damage. Thermoacoustic oscillations will have a low amplitude initially. Thus linear models can describe the initial growth and hence give stability predictions. An unstable linear mode will grow in amplitude until nonlinear effects become sufficiently important to achieve a limit cycle. While the frequency of the linear mode can often provide a good approximation to that of the resulting limit cycle, linear theories give no prediction of its resulting amplitude. In previous work, we developed a low-order frequency-domain method to model thermoacoustic limit cycles in LPP combustors. This was based on a ‘describing function’ approach and is only applicable when there is a dominant mode and the main nonlinearity is in the combustion response to flow perturbations. In this paper that method is extended into the time domain. The main advantage of the time-domain approach is that limit-cycle stability, the influence of harmonics, and the interaction between different modes can be simulated. In LPP combustion, fluctuations in the inlet fuel-air ratio have been shown to be the dominant cause of unsteady combustion: these occur because velocity perturbations in the premix ducts cause a time-varying fuel-air ratio, which then convects downstream. If the velocity perturbation becomes comparable to the mean flow, there will be an amplitude-dependent effect on the equivalence ratio fluctuations entering the combustor and hence on the rate of heat release. Since the Mach number is low, the velocity perturbation can be comparable to the mean flow, with even reverse flow occurring, while the disturbances are still acoustically linear in that the pressure perturbation is still much smaller than the mean. Hence while the combustion response to flow velocity and equivalence ratio fluctuations must be modelled nonlinearly, the flow perturbations generated as a result of the unsteady combustion can be treated as linear. In developing a time-domain network model for nonlinear thermoacoustic oscillations an initial frequency-domain calculation is performed. The linear network model, LOTAN, is used to categorise the combustor geometry by finding the transfer function for the response of flow perturbations (at the fuel injectors, say) to heat-release oscillations. This transfer function is then converted into the time domain through an inverse Fourier transform to obtain the Green’s function, which thus relates unsteady flow to heat release at previous times. By combining this with a nonlinear flame model (relating heat release to unsteady flow at previous times) a complete time-domain solution can be found by stepping forward in time. If an unstable mode is present, its amplitude will initially grow exponentially (in accordance with linear theory) until saturation effects in the flame model become significant, and eventually a stable limit cycle will be attained. The time-domain approach enables determination of the limit-cycle. In addition, the influence of harmonics and the interaction and exchange of energy between different modes can be simulated. These effects are investigated for longitudinal and circumferential instabilities in an example combustor system and results are compared to frequency-domain limit-cycle predictions.

A Time-Domain Network Model for Nonlinear Thermoacoustic Oscillations

2009 ◽  
Vol 131 (3) ◽  
Author(s):  
Simon R. Stow ◽  
Ann P. Dowling
Keyword(s):  
Heat Release ◽  
Frequency Domain ◽  
Limit Cycle ◽  
Network Model ◽  
Time Domain ◽  
Acoustic Waves ◽  
Mean Flow ◽  
The Mean ◽  
The Time Domain ◽  
Thermoacoustic Oscillations

Lean premixed prevaporized (LPP) combustion can reduce NOx emissions from gas turbines but often leads to combustion instability. Acoustic waves produce fluctuations in heat release, for instance, by perturbing the fuel-air ratio. These heat fluctuations will in turn generate more acoustic waves and in some situations linear oscillations grow into large-amplitude self-sustained oscillations. The resulting limit cycles can cause structural damage. Thermoacoustic oscillations will have a low amplitude initially. Thus linear models can describe the initial growth and hence give stability predictions. An unstable linear mode will grow in amplitude until nonlinear effects become sufficiently important to achieve a limit cycle. While the frequency of the linear mode can often provide a good approximation to that of the resulting limit cycle, linear theories give no prediction of its resulting amplitude. In previous work, we developed a low-order frequency-domain method to model thermoacoustic limit cycles in LPP combustors. This was based on a “describing-function” approach and is only applicable when there is a dominant mode and the main nonlinearity is in the combustion response to flow perturbations. In this paper that method is extended into the time domain. The main advantage of the time-domain approach is that limit-cycle stability, the influence of harmonics, and the interaction between different modes can be simulated. In LPP combustion, fluctuations in the inlet fuel-air ratio have been shown to be the dominant cause of unsteady combustion: These occur because velocity perturbations in the premix ducts cause a time-varying fuel-air ratio, which then convects downstream. If the velocity perturbation becomes comparable to the mean flow, there will be an amplitude-dependent effect on the equivalence ratio fluctuations entering the combustor and hence on the rate of heat release. Since the Mach number is low, the velocity perturbation can be comparable to the mean flow, with even reverse flow occurring, while the disturbances are still acoustically linear in that the pressure perturbation is still much smaller than the mean. Hence while the combustion response to flow velocity and equivalence ratio fluctuations must be modeled nonlinearly, the flow perturbations generated as a result of the unsteady combustion can be treated as linear. In developing a time-domain network model for nonlinear thermoacoustic oscillations an initial frequency-domain calculation is performed. The linear network model, LOTAN, is used to categorize the combustor geometry by finding the transfer function for the response of flow perturbations (at the fuel injectors, say) to heat-release oscillations. This transfer function is then converted into the time domain through an inverse Fourier transform to obtain Green’s function, which thus relates unsteady flow to heat release at previous times. By combining this with a nonlinear flame model (relating heat release to unsteady flow at previous times) a complete time-domain solution can be found by stepping forward in time. If an unstable mode is present, its amplitude will initially grow exponentially (in accordance with linear theory) until saturation effects in the flame model become significant, and eventually a stable limit cycle will be attained. The time-domain approach enables determination of the limit cycle. In addition, the influence of harmonics and the interaction and exchange of energy between different modes can be simulated. These effects are investigated for longitudinal and circumferential instabilities in an example combustor system and the results are compared with frequency-domain limit-cycle predictions.

scholarly journals Instability wave models for the near-field fluctuations of turbulent jets

2011 ◽  
Vol 689 ◽  
pp. 97-128 ◽  
Author(s):  
K. Gudmundsson ◽  
Tim Colonius
Keyword(s):  
Flow Field ◽  
Large Scale ◽  
Microphone Array ◽  
Near Field ◽  
Turbulent Jets ◽  
Instability Wave ◽  
Mean Flow ◽  
Velocity Fluctuations ◽  
The Mean ◽  
Scale Structures

AbstractPrevious work has shown that aspects of the evolution of large-scale structures, particularly in forced and transitional mixing layers and jets, can be described by linear and nonlinear stability theories. However, questions persist as to the choice of the basic (steady) flow field to perturb, and the extent to which disturbances in natural (unforced), initially turbulent jets may be modelled with the theory. For unforced jets, identification is made difficult by the lack of a phase reference that would permit a portion of the signal associated with the instability wave to be isolated from other, uncorrelated fluctuations. In this paper, we investigate the extent to which pressure and velocity fluctuations in subsonic, turbulent round jets can be described aslinearperturbations to the mean flow field. The disturbances are expanded about the experimentally measured jet mean flow field, and evolved using linear parabolized stability equations (PSE) that account, in an approximate way, for the weakly non-parallel jet mean flow field. We utilize data from an extensive microphone array that measures pressure fluctuations just outside the jet shear layer to show that, up to an unknown initial disturbance spectrum, the phase, wavelength, and amplitude envelope of convecting wavepackets agree well with PSE solutions at frequencies and azimuthal wavenumbers that can be accurately measured with the array. We next apply the proper orthogonal decomposition to near-field velocity fluctuations measured with particle image velocimetry, and show that the structure of the most energetic modes is also similar to eigenfunctions from the linear theory. Importantly, the amplitudes of the modes inferred from the velocity fluctuations are in reasonable agreement with those identified from the microphone array. The results therefore suggest that, to predict, with reasonable accuracy, the evolution of the largest-scale structures that comprise the most energetic portion of the turbulent spectrum of natural jets, nonlinear effects need only be indirectly accounted for by considering perturbations to the mean turbulent flow field, while neglecting any non-zero frequency disturbance interactions.

An Investigation of Flow Fields Over Multi-Element Aerofoils

2001 ◽  
Vol 124 (1) ◽  
pp. 154-165 ◽  
Author(s):  
S. R. Maddah ◽  
H. H. Bruun
Keyword(s):  
Flow Field ◽  
Computational Study ◽  
Separated Flow ◽  
Mean Flow ◽  
Model Configuration ◽  
Attached Flow ◽  
The Mean ◽  
Flap Deflection ◽  
Comparative Results ◽  
Lift And Drag

This paper presents results obtained from a combined experimental and computational study of the flow field over a multi-element aerofoil with and without an advanced slat. Detailed measurements of the mean flow and turbulent quantities over a multi-element aerofoil model in a wind tunnel have been carried out using stationary and flying hot-wire (FHW) probes. The model configuration which spans the test section 600mm×600mm, is made of three parts: 1) an advanced (heel-less) slat, 2) a NACA 4412 main aerofoil and 3) a NACA 4415 flap. The chord lengths of the elements were 38, 250 and 83 mm, respectively. The results were obtained at a chord Reynolds number of 3×105 and a free Mach number of less than 0.1. The variations in the flow field are explained with reference to three distinct flow field regimes: attached flow, intermittent separated flow, and separated flow. Initial comparative results are presented for the single main aerofoil and the main aerofoil with a nondeflected flap at angles of attacks of 5, 10, and 15 deg. This is followed by the results for the three-element aerofoil with emphasis on the slat performance at angles of attack α=10, 15, 20, and 25 deg. Results are discussed both for a nondeflected flap δf=0deg and a deflected flap δf=25deg. The measurements presented are combined with other related aerofoil measurements to explain the main interaction of the slat/main aerofoil and main aerofoil/flap both for nondeflected and deflected flap conditions. These results are linked to numerically calculated variations in lift and drag coefficients with angle of attack and flap deflection angle.

Measurement of the Mean Flow Field in a Smooth Rotating Channel With Coriolis and Buoyancy Effects

2017 ◽  
Author(s):  
Ruquan You ◽  
Haiwang Li ◽  
Zhi Tao ◽  
Kuan Wei
Keyword(s):  
Flow Field ◽  
Coriolis Force ◽  
Indium Tin Oxide ◽  
Buoyancy Force ◽  
Flow Direction ◽  
Density Ratio ◽  
Mean Flow ◽  
The Mean ◽  
Static Conditions ◽  
Rotating Channel

The mean flow field in a smooth rotating channel was measured by particle image velocimetry under the effect of buoyancy force. In the experiments, the Reynolds number, based on the channel hydraulic diameter (D) and the bulk mean velocity (Um), is 10000, and the rotation numbers are 0, 0.13, 0.26, 0.39, 0.52, respectively. The four channel walls are heated with Indium Tin Oxide (ITO) heater glass, making the density ratio (d.r.) about 0.1 and the maximum value of buoyancy number up to 0.27. The mean flow field was simulated on a 3D reconstruction at the position of 3.5<X/D<6.5, where X is along the mean flow direction. The effect of Coriolis force and buoyancy force on the mean flow was taken into consideration in the current work. The results show that the Coriolis force pushes the mean flow to the trailing side, making the asymmetry of the mean flow with that in the static conditions. On the leading surface, due to the effect of buoyancy force, the mean flow field changes considerably. Comparing with the case without buoyancy force, separated flow was captured by PIV on the leading side in the case with buoyancy force. More details of the flow field will be presented in this work.

scholarly journals Dynamic Response of Fuel Nozzles for Liquid-Fueled Gas Turbine Combustors

1996 ◽  
Author(s):  
Mounir Ibrahim ◽  
Terry Sanders ◽  
Douglas Darling ◽  
Michelle Zaller
Keyword(s):  
Gas Turbine ◽  
Cone Angle ◽  
Outer Region ◽  
Test Fluid ◽  
Mean Flow ◽  
Spray Cone Angle ◽  
Spray Cone ◽  
Gas Turbine Combustors ◽  
The Mean ◽  
Piezoelectric Drive

To imitate resonances that might occur in the fuel delivery system of gas turbine combustors, the incoming liquid streams of two pressure swirl nozzles were perturbed using a piezoelectric driver. Frequencies of perturbations examined were from 3 to 20 kHz, and water was used as the test fluid. A video camera and a Phase Doppler Particle Analyzer (PDPA) were used to study the effect of perturbations on the mean flow quantities of the sprays. Various lighting arrangements were used for the video photography: back lighting, front lighting, a strobe synchronized with the input to the piezoelectric, and a laser sheet oriented along the midplane of the sprays. The study showed that the piezoelectric drive had an effect an the spray system at discrete frequencies. At these particular frequencies, by increasing the input voltage, it was found that the piezoelectric drive affected the atomization in the following ways: (1) the mean flow rate decreased, (2) the spray cone angle decreased, (3) the break up length decreased, (4) the peak of the spatial distribution of the mean droplet size decreased, and (5) the mean droplet sizes and velocities increased near the spray center line and decreased in the outer region of the spray. A hysteresis effect of the drive frequency on the spray cone angle was observed. The results indicated that more fundamental research is needed to gain an in-depth understanding of the physical processes induced in the spray by the piezoelectric drive.

Sign in / Sign up

Export Citation Format

Share Document