Absolute and convective instability in gas turbine fuel injectors

Hydrodynamic instabilities in gas turbine fuel injectors help to mix the fuel and air but can sometimes lock into acoustic oscillations and contribute to thermoacoustic instability. This paper describes a linear stability analysis that predicts the frequencies and strengths of hydrodynamic instabilities and identifies the regions of the flow that cause them. It distinguishes between convective instabilities, which grow in time but are convected away by the flow, and absolute instabilities, which grow in time without being convected away. Convectively unstable flows amplify external perturbations, while absolutely unstable flows also oscillate at intrinsic frequencies. As an input, this analysis requires velocity and density fields, either from a steady but unstable solution to the Navier-Stokes equations, or from time-averaged numerical simulations. In the former case, the analysis is a predictive tool. In the latter case, it is a diagnostic tool. This technique is applied to three flows: a swirling wake at Re = 400, a single stream swirling fuel injector at Re - 106, and a lean premixed gas turbine injector with five swirling streams at Re - 106. Its application to the swirling wake demonstrates that this technique can correctly predict the frequency, growth rate and dominant wavemaker region of the flow. It also shows that the zone of absolute instability found from the spatio-temporal analysis is a good approximation to the wavemaker region, which is found by overlapping the direct and adjoint global modes. This approximation is used in the other two flows because it is difficult to calculate their adjoint global modes. Its application to the single stream fuel injector demonstrates that it can identify the regions of the flow that are responsible for generating the hydrodynamic oscillations seen in LES and experimental data. The frequencies predicted by this technique are within a few percent of the measured frequencies. The technique also explains why these oscillations become weaker when a central jet is injected along the centreline. This is because the absolutely unstable region that causes the oscillations becomes convectively unstable. Its application to the lean premixed gas turbine injector reveals that several regions of the flow are hydrodynamically unstable, each with a different frequency and a different strength. For example, it reveals that the central region of confined swirling flow is strongly absolutely unstable and sets up a precessing vortex core, which is likely to aid mixing throughout the injector. It also reveals that the region between the second and third streams is slightly absolutely unstable at a frequency that is likely to coincide with acoustic modes within the combustion chamber. This technique, coupled with knowledge of the acoustic modes in a combustion chamber, is likely to be a useful design tool for the passive control of mixing and combustion instability. Copyright © 2012 by ASME.

Surface-waves, stability, and scattering for a lined duct with flow

We consider a straight cylindrical duct with a steady subsonic axial flow and a reacting boundary (e.g. an acoustic lining). The wave modes are separated into ordinary acoustic duct modes, and surface modes confined to a small neighbourhood of the boundary. Many researchers have used a mass-spring-damper boundary model, for which one surface mode has previously been identified as a convective instability; however, we show the stability analysis used in such cases to be questionable. We investigate instead the stability of the surface modes using the Briggs-Bers criterion for a Flügge thin-shell boundary model. For modest frequencies and wavenumbers the thin-shell has an impedance which is effectively that of a mass-spring-damper, although for the large wavenumbers needed for the stability analysis the thin-shell and mass-spring-damper impedances diverge, owing to the thin shell's bending stiffness. The thin shell model may therefore be viewed as a regularization of the mass-spring-damper model which accounts for nonlocally-reacting effects. We find all modes to be stable for realistic thin-shell parameters, while absolute instabilities are demonstrated for extremely thin boundary thicknesses. The limit of vanishing bending stiffness is found to be a singular limit, yielding absolute instabilities of arbitrarily large temporal growth rate. We propose that the problems with previous stability analyses are due to the neglect of something akin to bending stiffness in the boundary model. Our conclusion is that the surface mode previously identified as a convective instability may well be stable in reality. Finally, inspired by Rienstra's recent analysis, we investigate the scattering of an acoustic mode as it encounters a sudden change from a hard-wall to a thin-shell boundary, using a Wiener-Hopf technique. The thin-shell is considered to be clamped to the hard-wall. The acoustic mode is found to scatter into transmitted and reflected acoustic modes, and surface modes strongly linked to the solid waves in the boundary, although no longitudinal or transverse waves within the boundary are excited. Examples are provided that demonstrate total transmission, total reflection, and a combination of the two. This thin-shell scattering problem is preferable to the mass-spring-damper scattering problem presented by Rienstra, since the thin-shell problem is fully determined and does not need to appeal to a Kutta-like condition or the inclusion of an instability in order to avoid a surface-streamline cusp at the boundary change.

Structural sensitivity of spiral vortex breakdown

Previous numerical simulations have shown that vortex breakdown starts with the formation of a steady axisymmetric bubble and that an unsteady spiralling mode then develops on top of this. We investigate this spiral mode with a linear global stability analysis around the steady bubble and its wake. We obtain the linear direct and adjoint global modes of the linearized Navier-Stokes equations and overlap these to obtain the structural sensitivity of the spiral mode, which identifies the wavemaker region. We also identify regions of absolute instability with a local stability analysis. At moderate swirls, we find that the m=-1 azimuthal mode is the most unstable and that the wavemaker regions of the m=-1 mode lie around the bubble, which is absolutely unstable. The mode is most sensitive to feedback involving the radial and azimuthal components of momentum in the region just upstream of the bubble. To a lesser extent, the mode is also sensitive to feedback involving the axial component of momentum in regions of high shear around the bubble. At an intermediate swirl, in which the bubble and wake have similar absolute growth rates, other researchers have found that the wavemaker of the nonlinear global mode lies in the wake. We agree with their analysis but find that the regions around the bubble are more influential than the wake in determining the growth rate and frequency of the linear global mode. The results from this paper provide the first steps towards passive control strategies for spiral vortex breakdown. © 2013 Cambridge University Press.

Local stability results for the collective behaviors of infinite populations of pulse-coupled oscillators

In this paper, we investigate the behavior of pulse-coupled integrate-and-fire oscillators. Because the stability analysis of finite populations is intricate, we investigate stability results in the approximation of infinite populations. In addition to recovering known stability results of finite populations, we also obtain new stability results for infinite populations. In particular, under a weak coupling assumption, we solve for the continuum model a conjecture still prevailing in the finite dimensional case. © 2011 IEEE.

Analysis of interconnected oscillators by dissipativity theory

This paper employs dissipativity theory for the global analysis of limit cycles in particular dynamical systems of possibly high dimension. Oscillators are regarded as open systems that satisfy a particular dissipation inequality. It is shown that this characterization has implications for the global stability analysis of limit cycle oscillations: i) in isolated oscillators, ii) in interconnections of oscillators, and iii) for the global synchrony analysis in interconnections of identical oscillators. © 2007 IEEE.

What kind of friction model is needed to predict brake squeal?

Predictive models of friction-induced vibration have proved elusive despite decades of research. There are many mechanisms that can cause brake squeal; friction coupled systems can be highly sensitive to small perturbations; and the dynamic properties of friction at the contact zone seem to be poorly understood. This paper describes experimental and theoretical work aimed at identifying the key ingredients of a predictive model. A large-scale experiment was carried out to identify squeal initiations using a pin-on-disc test rig: approximately 30,000 squeal initiations were recorded, covering a very wide range of frequencies. The theoretical model allows for completely general linear systems coupled at a single sliding point by friction: squeal is predicted using a linearised stability analysis. Results will be presented that show that almost all observed squeal events can be predicted within this model framework, but that some subsets require innovative friction modelling: predictions are highly dependent on the particular choice of friction model and its associated parameters. Copyright © 2012 by ASME.

What kind of friction model is needed to predict brake squeal?

Predictive models of friction-induced vibration have proved elusive despite decades of research. There are many mechanisms that can cause brake squeal; friction coupled systems can be highly sensitive to small perturbations; and the dynamic properties of friction at the contact zone seem to be poorly understood. This paper describes experimental and theoretical work aimed at identifying the key ingredients of a predictive model. A large-scale experiment was carried out to identify squeal initiations using a pin-on-disc test rig: approximately 30,000 squeal initiations were recorded, covering a very wide range of frequencies. The theoretical model allows for completely general linear systems coupled at a single sliding point by friction: squeal is predicted using a linearised stability analysis. Results will be presented that show that almost all observed squeal events can be predicted within this model framework, but that some subsets require innovative friction modelling: predictions are highly dependent on the particular choice of friction model and its associated parameters. Copyright © 2012 by ASME.

A differential lyapunov framework for contraction analysis

Lyapunov's second theorem is an essential tool for stability analysis of differential equations. The paper provides an analog theorem for incremental stability analysis by lifting the Lyapunov function to the tangent bundle. The Lyapunov function endows the state-space with a Finsler structure. Incremental stability is inferred from infinitesimal contraction of the Finsler metrics through integration along solutions curves. © 2013 IEEE.

Modal stability theory

This article contains a review of modal stability theory. It covers local stability analysis of parallel flows including temporal stability, spatial stability, phase velocity, group velocity, spatio-temporal stability, the linearized Navier-Stokes equations, the Orr-Sommerfeld equation, the Rayleigh equation, the Briggs-Bers criterion, Poiseuille flow, free shear flows, and secondary modal instability. It also covers the parabolized stability equation (PSE), temporal and spatial biglobal theory, 2D eigenvalue problems, 3D eigenvalue problems, spectral collocation methods, and other numerical solution methods. Computer codes are provided for tutorials described in the article. These tutorials cover the main topics of the article and can be adapted to form the basis of research codes. Copyright © 2014 by ASME.

Linear dimensionality reduction: Survey, insights, and generalizations

© 2015 John P. Cunningham and Zoubin Ghahramani. Linear dimensionality reduction methods are a cornerstone of analyzing high dimensional data, due to their simple geometric interpretations and typically attractive computational properties. These methods capture many data features of interest, such as covariance, dynamical structure, correlation between data sets, input-output relationships, and margin between data classes. Methods have been developed with a variety of names and motivations in many fields, and perhaps as a result the connections between all these methods have not been highlighted. Here we survey methods from this disparate literature as optimization programs over matrix manifolds. We discuss principal component analysis, factor analysis, linear multidimensional scaling, Fisher's linear discriminant analysis, canonical correlations analysis, maximum autocorrelation factors, slow feature analysis, sufficient dimensionality reduction, undercomplete independent component analysis, linear regression, distance metric learning, and more. This optimization framework gives insight to some rarely discussed shortcomings of well-known methods, such as the suboptimality of certain eigenvector solutions. Modern techniques for optimization over matrix manifolds enable a generic linear dimensionality reduction solver, which accepts as input data and an objective to be optimized, and returns, as output, an optimal low-dimensional projection of the data. This simple optimization framework further allows straightforward generalizations and novel variants of classical methods, which we demonstrate here by creating an orthogonal-projection canonical correlations analysis. More broadly, this survey and generic solver suggest that linear dimensionality reduction can move toward becoming a blackbox, objective-agnostic numerical technology.

Development of the 2003 CU-HTK conversational telephone speech transcription system

This paper describes the development of the 2003 CU-HTK large vocabulary speech recognition system for Conversational Telephone Speech (CTS). The system was designed based on a multi-pass, multi-branch structure where the output of all branches is combined using system combination. A number of advanced modelling techniques such as Speaker Adaptive Training, Heteroscedastic Linear Discriminant Analysis, Minimum Phone Error estimation and specially constructed Single Pronunciation dictionaries were employed. The effectiveness of each of these techniques and their potential contribution to the result of system combination was evaluated in the framework of a state-of-the-art LVCSR system with sophisticated adaptation. The final 2003 CU-HTK CTS system constructed from some of these models is described and its performance on the DARPA/NIST 2003 Rich Transcription (RT-03) evaluation test set is discussed.

Automatic transcription of conversational telephone speech

This paper discusses the Cambridge University HTK (CU-HTK) system for the automatic transcription of conversational telephone speech. A detailed discussion of the most important techniques in front-end processing, acoustic modeling and model training, language and pronunciation modeling are presented. These include the use of conversation side based cepstral normalization, vocal tract length normalization, heteroscedastic linear discriminant analysis for feature projection, minimum phone error training and speaker adaptive training, lattice-based model adaptation, confusion network based decoding and confidence score estimation, pronunciation selection, language model interpolation, and class based language models. The transcription system developed for participation in the 2002 NIST Rich Transcription evaluations of English conversational telephone speech data is presented in detail. In this evaluation the CU-HTK system gave an overall word error rate of 23.9%, which was the best performance by a statistically significant margin. Further details on the derivation of faster systems with moderate performance degradation are discussed in the context of the 2002 CU-HTK 10 × RT conversational speech transcription system. © 2005 IEEE.

Calculations and predications of combustion oscillations

The combustion oscillations are the phenomena which we may meet in developing the clean, safe and efficient power and propulsion systems. This paper summarizes authors' systematic work on fuel spray combustion oscillations in the recent years. Combining CFD calculations and stability analysis, a new approach of predicting combustion stabilities was developed. With this approach, detailed flow information and unstable modes can be obtained by CFD and solving perturbation equations, respectively. The results provide the guidance on understanding combustion instability mechanisms and developing the control strategies.

Triggering in a thermoacoustic system with stochastic noise

This paper explores the mechanism of triggering in a simple thermoacoustic system, the Rijke tube. It is demonstrated that additive stochastic perturbations can cause triggering before the linear stability limit of a thermoacoustic system. When triggering from low noise amplitudes, the system is seen to evolve to self-sustained oscillations via an unstable periodic solution of the governing equations. Practical stability is introduced as a measure of the stability of a linearly stable state when finite perturbations are present. The concept of a stochastic stability map is used to demonstrate the change in practical stability limits for a system with a subcritical bifurcation, once stochastic terms are included. The practical stability limits are found to be strongly dependent on the strength of noise.

Stabilizing effect of surrounding gas flow on a plane liquid sheet

The stability of a plane liquid sheet is studied experimentally and theoretically, with an emphasis on the effect of the surrounding gas. Co-blowing with a gas velocity of the same order of magnitude as the liquid velocity is studied, in order to quantify its effect on the stability of the sheet. Experimental results are obtained for a water sheet in air at Reynolds number Rel = 3000 and Weber number W e = 300, based on the half-thickness of the sheet at the inlet, water mean velocity at the inlet, the surface tension between water and air and water density and viscosity. The sheet is excited with different frequencies at the inlet and the growth of the waves in the streamwise direction is measured. The growth rate curves of the disturbances for all air flow velocities under study are found to be within 20 % of the values obtained from a local spatial stability analysis, where water and air viscosities are taken into account, while previous results from literature assuming inviscid air overpredict the most unstable wavelength with a factor 3 and the growth rate with a factor 2. The effect of the air flow on the stability of the sheet is scrutinized numerically and it is concluded that the predicted disturbance growth scales with (i) the absolute velocity difference between water and air (inviscid effect) and (ii) the square root of the shear from air on the water surface (viscous effect).