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 10⁵–10⁶ 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.

A detailed understanding of liquid propellant combustion is necessary for the development of improved and more reliable propulsion systems. This article describes experimental investigations aimed at providing such a fundamental basis for design and engineering of combustion components. It reports recent applications of imaging techniques to cryogenic combustion at high pressure. The flame structure is investigated in the transcritical range where the pressure exceeds the critical pressure of oxygen (p > p _c (O₂ = 5.04MPa)) but the temperature of the injected liquid oxygen is below its critical value . Data obtained from imaging of OH* radicals emission, CH* radicals emission in the case of LOx/GCH₄ flames and backlighting provide a detailed view of the flame structure for a set of injection conditions. The data may be used to guide numerical modelling of transcritical flames and the theoretical and numerical analysis of the stabilization process. Calculations of the flame edge are used to illustrate this aspect. Results obtained may also be employed to devise engineering modelling tools and methodologies for component development aimed at improved efficiency and augmented reliability.

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.

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 spiralling mode). The most energetic POD mode pair is identified as the precessing vortex core. By analysing the fast Fourier transform (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 Kelvin–Helmholtz 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.