Matthew P. Juniper

Triggering in the horizontal Rijke tube: non-normality, transient growth and bypass transition

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 selfsustained 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 5 –10 6 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.

STRUCTURE AND DYNAMICS OF CRYOGENIC FLAMES AT SUPERCRITICAL PRESSURE

ABSTRACT 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 (O2 = 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/GCH4 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.

The effect of confinement on the stability of two-dimensional shear flows

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 stability of ducted compound flows and consequences for the geometry of coaxial injectors

A two-dimensional wake-like compound flow, formed by a low-speed stream embedded within a high-speed flow, is examined in this article. The range of absolutely unstable flow in parameter space greatly increases when such a flow is confined within a duct. Parameters studied here are: the density ratio, which is from 0.1 to 1000; the velocity ratio, which varies from co-flow to counter-flow; and the ratio of the duct width to the width of the central jet. Absolutely unstable flows permit perturbations to propagate upstream, and can lead to self-sustained global oscillations similar to the vortex shedding process which takes place in the wake of a bluff body. This theoretical situation models the wake-like behaviour of a two-fluid coaxial injector with a recessed central tube. The aerodynamic destabilizing mechanism is of primary importance whereas the stabilizing mechanisms, which are not considered here, are of secondary importance. The conclusions from this analysis of a ducted compound flow can explain why one observes self-sustained oscillations in recessed coaxial injectors

Adjoint algorithms for the Navier-Stokes equations in the low Mach number limit

This paper describes a derivation of the adjoint low Mach number equations and their implementation and validation within a global mode solver. The advantage of using the low Mach number equations and their adjoints is that they are appropriate for flows with variable density, such as flames, but do not require resolution of acoustic waves. Two versions of the adjoint are implemented and assessed: a discrete-adjoint and a continuous-adjoint. The most unstable global mode calculated with the discrete-adjoint has exactly the same eigenvalue as the corresponding direct global mode but contains numerical artifacts near the inlet. The most unstable global mode calculated with the continuous-adjoint has no numerical artifacts but a slightly different eigenvalue. The eigenvalues converge, however, as the timestep reduces. Apart from the numerical artifacts, the mode shapes are very similar, which supports the expectation that they are otherwise equivalent. The continuous-adjoint requires less resolution and usually converges more quickly than the discrete-adjoint but is more challenging to implement. Finally, the direct and adjoint global modes are combined in order to calculate the wavemaker region of a low density jet.

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.

Azimuthal instabilities in annular combustors: Standing and spinning modes

This theoretical study investigates spinning and standing modes in azimuthally symmetric annular combustion chambers. Both modes are observed in experiments and simulations, and an existing model predicts that spinning modes are the only stable state of the system. We extend this model to take into account the effect that the acoustic azimuthal velocity, u, has on the flames, and propose a phenomenological model based on experiments performed on transversely forced flames. This model contains a parameter, δ, that quantifies the influence that the transversal excitation has on the fluctuating heat release. For small values of δ, spinning modes are the only stable state of the system. In an intermediate range of δ, both spinning and standing modes are stable states. For large values of δ, standing modes are the only stable state. This study shows that a flame’s response to azimuthal velocity fluctuations plays an important role in determining the type of thermoacoustic oscillations found in annular combustors.

Sensitivity analysis of a time-delayed thermo-acoustic system via an adjoint-based approach

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 time-delay 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.

The full impulse response of two-dimensional jet/wake flows and implications for confinement

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.

The effect of confinement on the stability of viscous planar jets and wakes

This theoretical study examines confined viscous planar jet/wake flows with continuous velocity profiles. These flows are characterized by the shear, confinement, Reynolds number and shear-layer thickness. The primary aim of this paper is to determine the effect of confinement on viscous jets and wakes and to compare these results with corresponding inviscid results. The secondary aim is to consider the effect of viscosity and shear-layer thickness. A spatio-temporal analysis is performed in order to determine absolute/convective instability criteria. This analysis is carried out numerically by solving the Orr-Sommerfeld equation using a Chebyshev collocation method. Results are produced over a large range of parameter space, including both co-flow and counter-flow domains and confinements corresponding to 0.1 < h 2 /h 1 < 10, where the subscripts 1 and 2 refer to the inner and outer streams, respectively. The Reynolds number, which is defined using the channel width, takes values between 10 and 1000. Different velocity profiles are used so that the shear layers occupy between 1/2 and 1/24 of the channel width. Results indicate that confinement has a destabilizing effect on both inviscid and viscous flows. Viscosity is found always to be stabilizing, although its effect can safely be neglected above Re = 1000. Thick shear layers are found to have a stabilizing effect on the flow, but infinitely thin shear layers are not the most unstable; having shear layers of a small, but finite, thickness gives rise to the strongest instability.