International audienceThe decomposition of experimental data into dynamic modes using a data-based algorithm is applied to Schlieren snapshots of a helium jet and to time-resolved PIV-measurements of an unforced and harmonically forced jet. The algorithm relies on the reconstruction of a low-dimensional inter-snapshot map from the available flow field data. The spectral decomposition of this map results in an eigenvalue and eigenvector representation (referred to as dynamic modes) of the underlying fluid behavior contained in the processed flow fields. This dynamic mode decomposition allows the breakdown of a fluid process into dynamically revelant and coherent structures and thus aids in the characterization and quantification of physical mechanisms in fluid flow. © 2010 Springer-Verlag

With a sufficiently large impulse, a thermoacoustic system can reach self-sustained oscillations even when it is linearly stable, a process known as triggering. In this paper, a procedure is developed to find the lowest initial energy that can trigger self-sustained oscillations, as well as the corresponding initial state. This is known as the ‘most dangerous’ initial state. The procedure is based on adjoint looping of the nonlinear governing equations, combined with an optimization routine. It is developed for a simple model of a thermoacoustic system, the horizontal Rijke tube, and can be extended to more sophisticated thermoacoustic models. It is observed that the most dangerous initial state grows transiently towards an unstable periodic solution before growing to a stable periodic solution. The initial energy required to trigger these self-sustained oscillations is much lower than the energy of the oscillations themselves and slightly lower than the lowest energy on the unstable periodic solution. It is shown that this transient growth arises due to non-normality of the governing equations. This is analogous to the sequence of events observed in bypass transition to turbulence in fluid mechanical systems and has the same underlying cause. The most dangerous initial state is calculated as a function of the heat-release parameter. It is found that self-sustained oscillations can be reached over approximately half the linearly stable domain. Transient growth in real thermoacoustic systems is 105–106 times greater than that in this simple model. One practical conclusion is that, even in the linearly stable regime, it may take very little initial energy for a real thermoacoustic system to trigger to high-amplitude self-sustained oscillations through the mechanism described in this paper.

Nine decades of rocket engine and gas turbine development have shown that thermoacoustic oscillations are difficult to predict, but can usually be eliminated with relatively small ad hoc design changes. These changes can, however, be ruinously expensive to devise. This review explains why linear and nonlinear thermoacoustic behaviour is so sensitive to parameters such as operating point, fuel composition, and injector geometry. It shows how non-periodic behaviour arises in experiments and simulations and discusses how fluctuations in thermoacoustic systems with turbulent reacting flow, which are usually filtered or averaged out as noise, can reveal useful information. Finally, it proposes tools to exploit this sensitivity in the future: adjoint-based sensitivity analysis to optimize passive control designs, and complex systems theory to warn of impending thermoacoustic oscillations and to identify the most sensitive elements of a thermoacoustic system.

We apply adjoint-based sensitivity analysis to a time-delayed thermo-acoustic system: a Rijke tube containing a hot wire. We calculate how the growth rate and frequency of small oscillations about a base state are affected either by a generic passive control element in the system (the structural sensitivity analysis) or by a generic change to its base state (the base-state sensitivity analysis). We illustrate the structural sensitivity by calculating the effect of a second hot wire with a small heat release parameter. In a single calculation, this shows how the second hot wire changes the growth rate and frequency of the small oscillations, as a function of its position in the tube. We then examine the components of the structural sensitivity in order to determine the passive control mechanism that has the strongest influence on the growth rate. We find that a force applied to the acoustic momentum equation in the opposite direction to the instantaneous velocity is the most stabilizing feedback mechanism. We also find that its effect is maximized when it is placed at the downstream end of the tube. This feedback mechanism could be supplied, for example, by an adiabatic mesh. We illustrate the base-state sensitivity by calculating the effects of small variations in the damping factor, the heat-release timedelay coefficient, the heat-release parameter, and the hot wire location. The successful application of sensitivity analysis to thermo-acoustics opens up new possibilities for the passive control of thermo-acoustic oscillations by providing gradient information that can be combined with constrained optimization algorithms in order to reduce linear growth rates.

It has been shown recently that the instability of a two-dimensional wake increases when it is confined in the transverse direction by two flat plates. Confinement causes the transition from convective to absolute instability to occur at lower values of shear. This paper examines this effect comprehensively and concludes that it is due to the constructive interaction of modes with zero group velocity in the wake (or jet) and in the surrounding flow. Maximum instability occurs when the wavenumber of the fundamental mode in the wake (or jet) matches that of the fundamental mode in the surrounding flow. Other regions of high instability occur when the harmonics of one mode interact with the fundamental of the other. This effect is examined at density ratios from 0.001 to 1000. At each density ratio, the confinement which causes maximum absolute instability can be predicted. This study also shows that it is vital to examine the wavenumber of absolutely unstable modes in order to avoid over-predicting the absolute instability. In some situations this wavenumber is vanishingly small and the mode must be discounted on physical grounds.

The large-scale coherent motions in a realistic swirl fuel injector geometry are analysed by direct numerical simulations (DNS), proper orthogonal decomposition (POD), and linear global modes. The aim is to identify the origin of instability in this turbulent flow in a complex internal geometry.The flow field in the nonlinear simulation is highly turbulent, but with a distinguishable coherent structure: the precessing vortex core (a spiraling mode).The most energetic POD mode pair is identified as the precessing vortex core. By analysing the FFT of the time coefficients of the POD modes, we conclude that the first four POD modes contain the coherent fluctuations. The remaining POD modes (incoherent fluctuations) are used to form a turbulent viscosity field, using the Newtonian eddy model.The turbulence sets in from convective shear layer instabilities even before the nonlinear flow reaches the other end of the domain, indicating that equilibrium solutions of the Navier-Stokes are never observed. Linear global modes are computed around the mean flow from DNS, applying the turbulent viscosity extracted from POD modes. A slightly stable discrete m = 1 eigenmode is found, well separated from the continuous spectrum, in very good agreement with the POD mode shape and frequency. The structural sensitivity of the precessing vortex core is located upstream of the central recirculation zone, identifying it as a spiral vortex breakdown instability in the nozzle. Furthermore, the structural sensitivity indicates that the dominant instability mechanism is the KelvinHelmholtz instability at the inflection point forming near vortex breakdown. Adjoint modes are strong in the shear layer along the whole extent of the nozzle, showing that the optimal initial condition for the global mode is localized in the shear layer.We analyse the qualitative influence of turbulent dissipation in the stability problem (eddy viscosity) on the eigenmodes by comparing them to eigenmodes computed without eddy viscosity. The results show that the eddy viscosity improves the complex frequency and shape of global modes around the fuel injector mean flow, while a qualitative wavemaker position can be obtained with or without turbulent dissipation, in agreement with previous studies.This study shows how sensitivity analysis can identify which parts of the flow in a complex geometry need to be altered in order to change its hydrodynamic stability characteristics.Key Words: † Email address for correspondence: outi-leena.tammisola@nottingham.ac.uk , and the exit velocity from the swirler (U E ). The coordinate system is also defined: the origin is at the centerline at the swirler exit location. The relative dimensions of the geometry are the same as in the numerical simulation, except that the numerical domain is longer in the downstream direction.

Thermoacoustic systems can oscillate self-excitedly, and often non-periodically, owing to coupling between unsteady heat release and acoustic waves. We study a slot-stabilized two-dimensional premixed flame in a duct via numerical simulations of a G-equation flame coupled with duct acoustics. We examine the bifurcations and routes to chaos for three control parameters: (i) the flame position in the duct, (ii) the length of the duct and (iii) the mean flow velocity. We observe period-1, period-2, quasi-periodic and chaotic oscillations. For certain parameter ranges, more than one stable state exists, so mode switching is possible. At intermediate times, the system is attracted to and repelled from unstable states, which are also identified. Two routes to chaos are established for this system: the period-doubling route and the Ruelle-Takens-Newhouse route. These are corroborated by analyses of the power spectra of the acoustic velocity. Instantaneous flame images reveal that the wrinkles on the flame surface and pinch-off of flame pockets are regular for periodic oscillations, while they are irregular and have multiple time and length scales for quasi-periodic and aperiodic oscillations. This study complements recent experiments by providing a reduced-order model of a system with approximately 5000 degrees of freedom that captures much of the elaborate nonlinear behaviour of ducted premixed flames observed in the laboratory.

In this theoretical study, a linear spatio-temporal analysis is performed on unconfined and confined inviscid jet/wake flows in order to determine whether they are absolutely or convectively unstable. The impulse response is evaluated in the entire outer fluid, rather than just at the point of impulse, over a wide range of density ratios. This confirms that the dominant saddle point can validly migrate into the plane of diverging eigenfunctions. This reveals that, at certain density ratios and shear numbers, the response can grow upstream in some directions with a cross-stream component, even though it decays directly upstream of, and at, the point of impulse. This type of flow is convectively unstable when unconfined, but becomes absolutely unstable when confined. Other effects of confinement are described in a previous paper. Together, these articles have important implications for the design of fuel injectors, which often employ confined shear flows at high Reynolds number and large density ratios to generate strong mixing in combustion chambers.