Nonlinear simulations of combustion instabilities with a quasi-1D Navier-Stokes code

As lean premixed combustion systems are more susceptible to combustion instabilities than non-premixed systems, there is an increasing demand for improved numerical design tools that can predict the occurrence of combustion instabilities with high accuracy. The inherent nonlinearities in combustion instabilities can be of crucial importance, and we here propose an approach in which the one-dimensional (1D) Navier-Stokes and scalar transport equations are solved for geometries of variable cross-section. The focus is on attached flames, and for this purpose a new phenomenological model for the unsteady heat release from a flame front is introduced. In the attached flame method (AFM) the heat release occurs over the full length of the flame. The nonlinear code with the use of the AFM approach is validated against analytical results and against an experimental study of thermoacoustic instabilities in oxy-fuel flames by Ditaranto and Hals [Combustion and Flame 146 (2006) 493-512]. The numerical simulations are in accordance with the experimental measurements and the analytical results and both the frequencies and the amplitudes of the resonant acoustic pressure modes are reproduced with good accuracy.     

Phase synchronization and collective interaction of multiple flamelets in a laboratory scale spray combustor

In this paper, we investigate the coupled behvior of the acoustic field in the confinement and the unsteady flame dynamics in a laboratory scale spray combustor. We study this interaction during the intermittency route to thermoacoustic instability when the location of the flame is varied inside the combustor. As the flame location is changed, the synchronization properties of the coupled acoustic pressure and heat release rate signals change from desynchronized aperiodicity (combustion noise) to phase synchronized periodicity (thermoacoustic instability) through intermittent phase synchronization (intermittency). We also characterize the collective interaction between the multiple flamelets anchored at the flame holder and the acoustic field in the system, during different dynamical states observed in the combustor operation. When the signals are desynchronized, we notice that the flamelets exhibit a steady combustion without the exhibition of a prominent feedback with the acoustic field. In a state of intermittent phase synchronization, we observe the existence of a short-term coupling between the heat release rate and the acoustic field. We notice that the onset of collective synchronization in the oscillations of multiple flamelets and the acoustic field leads to the simultaneous emergence of periodicity in the global dynamics of the system. This collective periodicity in both the subsystems causes enhancement of oscillations during epochs of amplitude growth in the intermittency signal. On the contrary, the weakening of the coupling induces suppression of periodic oscillations during epochs of amplitude decay in the intermittency signal. During phase synchronization, we notice a sustained synchronized movement of all flamelets with the periodicity of the acoustic cycle in the system.     

“Slow” active control of combustion instabilities by modification of liquid fuel spray properties

This paper describes an experimental investigation of the feasibility of using “slow” active control approaches, which “instantaneously” change liquid fuel spray properties, to suppress combustion instabilities. The objective of this control approach was to break up the feedback between the combustion process heat release and combustor pressure oscillations that drive the instability by changing the characteristics of the combustion process (e.g., the characteristic combustion time). To demonstrate the feasibility of such control, this study used a proprietary fuel injector (Nanomiser^TM), which can vary its fuel spray properties, to investigate the dependence of acoustics–combustion process coupling, i.e., the driving of combustion instabilities, upon the fuel spray properties. This study showed that by changing the spray characteristics it is possible to significantly damp combustion instabilities. Furthermore, using combustion zone chemiluminescence distributions, which were obtained by Abel’s deconvolution synchronized with measured acoustic data, it has been shown that the instabilities were mostly driven midway between the combustor centerline and wall, a short distance downstream from the flame holder, where the mean axial flow velocity is approximately zero in the vortex near the flame holder. The results of this study strongly suggest that a “slow” active control system that employs controllable fuel injectors could be effectively used to prevent the onset of detrimental combustion instabilities.     

Examination of mode shapes in an unstable model combustor

The coupling between the fluid dynamics, heat addition, and the acoustics of a combustor system determine whether it is prone toward combustion instability. This paper presents results from a benchmark study of the eigenmodes in an unstable experimental combustor. The axisymmetric combustor configuration is representative of a number of practical systems and comprises an injector tube, geometric expansion into a combustion chamber, and a short converging nozzle. Instability limit cycle amplitudes ranged from 5% to nearly 50% of the mean 2.2 MPa pressure. Multiple harmonics were measured for the highly unstable cases. The model combustor was designed to provide a fairly comprehensive set of tested effects: sonic vs subsonic inlets; oxidizer tube lengths that were either quarter-wave, half-wave, or off-resonant acoustic equivalents to the combustion chamber; a significant injector mean flow with Ma∼0.4; and a varied combustion chamber length. The measured mode shape data were analyzed and reduced to provide comparison with results from a linearized one-dimensional Euler model, which included the effects of real boundary conditions, entropy generation, area change, and heat and mass addition, but did not include a model for unsteady heat addition. For low-amplitude instabilities, the measured resonance frequencies agreed with those calculated by the model for the injector tube-combustion chamber system. Resonance frequencies for the high-amplitude oscillation cases corresponded to the first longitudinal frequency of the combustion chamber and its integer multiples. Good quantitative agreement was obtained between computed and measured phase difference profiles, and mode envelopes agreed qualitatively. These results provide a basis for subsequent combustion response studies on the effects of unsteady heat addition.     

Flame-acoustic response measurements in a high-pressure, 42-injector,cryogenic rocket thrust chamber

This work presents measurements of acoustically driven flame dynamics in a 42-element, cryogenic oxygen-hydrogen rocket thrust chamber under supercritical injection conditions. The experiment shows self-excited combustion instabilities for certain operating conditions, and this work describes the nature of the flame dynamics driving the acoustic field, as far as it can be ascertained from state-of-the-art optical measurements. Optical access has been realized in the combustion chamber with both fibre-optical probes and a viewing window. The probes collect point-like measurements of filtered OH* radiation. Their signals were used to calculate the gain and phase of intensity oscillations with respect to acoustic pressure for both stable and unstable operating conditions. Through the window, synchronized high-speed imaging of the flame in filtered OH* and blue radiation wavelengths was collected. The 2D flame response was related to the local acoustic pressure to investigate the distributed intensity and phase relationships. The flame response from OH* measurements is in agreement with the theory of Rayleigh. For stable conditions the oscillations of combustion and pressure were out of phase, whereas for an excited chamber 1T mode the oscillations were closely in phase. The integrated Rayleigh index from blue imaging was not consistent with the OH* results. The reason lies in the depth of field captured by this type of imaging, and must be used in a complementary fashion together with OH* imaging. The flame response values and 2D visualization presented in this work are expected to be of value for the validation of numerical modelling of combustion instabilities.     

Feasibility of “intermittent” active control of combustion instabilities in liquid fueled combustors using a “smart” fuel injector

This paper describes an experimental investigation of the feasibility of an “intermittent” active control approach for suppressing combustion instabilities in liquid fueled combustors. The developed controller employs a “smart” fuel injector that can modify the spray properties in response to changes in combustor operating conditions. This action weakens or breaks up the coupling between the combustion process and combustor acoustic modes oscillations, thus preventing the excitation of large amplitude instabilities. This approach differs significantly from previously proposed active control methods, both in concept and implementation, as it requires only “intermittent” modification of the combustion process by a single control action as opposed to the continuous action required by most other active control methods. The “smart” fuel injector used in this study consisted of a double-staged, air-assisted atomizer in which counter swirling, primary (inner stage) and secondary (outer stage) air streams were supplied to the injector through separate sets of tangentially oriented orifices. Control of the ratio of air mass flow rates supplied to these two stages, by use of a diverter valve, resulted in significant changes in the spray shape and its axial, tangential, and radial velocity components. This variation in spray properties of the “smart” injector was characterized for different values of the inner to outer air flow rate ratio in cold flow tests with a PDPA system. These results were then correlated with the characteristics of the “intermittently” controlled combustor. Measured quantities included the instability amplitudes, axial dependence of the mean and oscillatory heat release amplitudes, and the characteristics of the recirculation zones, which were all shown to depend on the fuel spray properties. The results of this study demonstrate the feasibility of using “smart” fuel injectors with capabilities for varying the combustion process characteristics to reduce the amplitudes of detrimental combustion instabilities in real engines to acceptable levels.     

Fully compressible simulations of the impact of acoustic waves on the dynamics of laminar premixed flames for engine-relevant conditions

Thermo-acoustic instabilities remain problematic in the design of propulsion systems such as gas turbine engines, rocket motors, and ramjets. They arise from the constructive interaction of heat release rate and acoustic pressure oscillations, and can result in increased noise and mechanical fatigue. In the present work, we are concerned with the flame response to the thermodynamic fluctuations that accompany an incident acoustic wave. The objective is to investigate the flame dynamics under engine-relevant conditions using high-fidelity numerical simulations and detailed chemical kinetics. The focus is placed on the combustion of hydrogen and n-heptane, as they are both of practical interest and behave very differently when subjected to acoustic waves. We extract the phase and gain of the unsteady heat release response, which are directly related to the Rayleigh criterion and thus the stability of the system. We highlight the differences between results obtained using the fully compressible Navier-Stokes equations and the low Mach number approximation. The two simulation frameworks agree very well for acoustic wavelengths much larger than the flame thickness. However, they differ significantly at high frequencies. The gain erroneously reaches a plateau under the low Mach number approximation, while it decays to zero using the fully compressible framework. This difference is attributed to the spatial variations in the acoustic pressure, which are not captured by the low Mach number approximation.     

Time-domain characterization of the acoustic damping of a perforated liner with bias flow

Combustion instabilities are caused by the interaction of unsteady heat releases and acoustic waves. To mitigate combustion instabilities, perforated liners, typically subjected to a low Mach number bias flow (a cooling flow through perforated holes), are fitted along the bounding walls of a combustor. They dissipate the acoustic waves by generating vorticity at the rims of perforated apertures. To investigate the absorption of plane waves by a perforated liner with bias flow, a time-domain numerical model of a cylindrical lined duct is developed. The liners' damping mechanism is characterized by using a time-domain "compliance." The development of such time-domain compliance is based on simplified or unsimplified Rayleigh conductivity. Numerical simulations of two different configurations of lined duct systems are performed by combining a 1D acoustic wave model with the compliance model. Comparison is then made between the results from the present models, and those from the experiment and the frequency-domain model of previous investigation [Eldredge and Dowling, J. Fluid Mech. 485, 307-335(2003)]. Good agreement is observed. This confirms that the present model can be used to simulate the propagation and dissipation of acoustic plane waves in a lined duct in real-time.     

Experimental characterization of onset of acoustic instability in a nonpremixed half-dump combustor

This paper reports work on a nonpremixed half-dump combustor, in which methane is injected at the backward-facing step, and mixes and burns with the air flowing past the step in the unsteady recirculation zone. The flow and geometric parameters are widely varied, to gradually change from conditions of low-amplitude noise to excitation of high-amplitude discrete tones. The purpose of the work is to focus on the transition from the former condition to the latter, and to mark the onset of instability. Dimensionless groups such as the Helmholtz and Strouhal numbers are formed based on the observed dominant frequencies, whose variation with the air flow Reynolds number is used to identify the oscillations as those due to the natural acoustic modes or the vortex shedding process. High-speed chemiluminescence imaging reveals shedding of vortical structures in the flame zone. With variation in the conditions, flow-acoustic lock-on and transition from one vortex shedding mode to another is marked by nonlinearity in the corresponding amplitude variations. Such conditions are identified as the onset of instability in terms of the ratio of the flow time scale to the acoustic time scale and mapped against the operating fuel-air equivalence ratio of the combustor.     

Dynamics of spray and swirling flame under acoustic oscillations : A joint experimental and LES investigation

The dynamics of spray swirling flames is investigated by combining experiments on a single sector generic combustor and large eddy simulations of the same configuration. Measurements and calculations correspond to a self-sustained limit cycle operation where combustion coupled by an axial quarter wave acoustic mode induces large amplitude oscillations of pressure in the system. A detailed analysis of the mechanisms controlling the process is carried out first by comparing the measured and calculated spray and flame dynamics. Considering in a second stage that the spray and flame are compact with respect to the acoustic wavelength the analysis can be simplified by defining state variables that are obtained by taking averages over the combustor cross section and representing the behavior of these average quantities as a function of the axial coordinate and time. This reveals a first region in which essentially convective processes prevail. The convective heat release rate then couples further downstream with the pressure field giving rise to positive Rayleigh source terms which feed energy in the axial acoustic mode. In the convective region, the swirl number features oscillations around its mean value with an impact on the flow aerodynamics and flame radial displacement. Fluctuations in the fuel flow rate are initiated at the injector exhaust and likewise convected downstream. The total mass flow rate that exhibits strong convective disturbances is dominated further downstream by the acoustic motion. This information provides new insights on the convective-acoustic coupling that controls the heat release rate disturbances and reveals the time delays governing the combustion oscillation process.     

Impact of turbulence on flame brush development of acoustically excited rod-stabilized flames

Coherent structures, such as those arising from hydrodynamic instabilities or excited by thermoacoustic oscillations, can significantly impact flame structure and, consequently, the nature of heat release. The focus of this work is to study how coherent oscillations of varying amplitudes can impact the growth of the flame brush in a bluff-body stabilized flame and how this impact is influenced by the free stream turbulence intensity of the flow approaching the bluff body. We do this by providing external acoustic excitation at the natural frequency of vortex shedding to simulate a highly-coupled thermoacoustic instability, and we vary the in-flow turbulence intensity using perforated plates upstream of the flame. We use high-speed stereoscopic particle image velocimetry to obtain the three-component velocity field and we use the Mie-scattering images to quantify the behavior of the flame edge. Our results show that in the low-turbulence conditions, presence of high-amplitude acoustic excitation can cause the flame brush to exhibit a step-function growth, indicating that the presence of strong vortical structures close to the flame can suppress flame brush growth. This impact is strongly dependent on the in-flow turbulence intensity and the flame brush development in conditions with higher levels of in-flow turbulence are minimally impacted by increasing amplitudes of acoustic excitation. These findings suggest that the sensitivity of the flow and flame to high-amplitude coherent oscillations is a strong function of the in-flow turbulence intensity.     

黎开管热声振荡的主动控制研究

本文以黎开管内的热声耦合振荡为研究对象,设计基于主动补偿的适应性控制器抑制黎开管内的不稳定燃烧。试验以扬声器为执行机构来改变黎开管的边界条件,从而抑制黎开管内的热声耦合振荡。实时控制效果表明,本文所采用的适应性控制算法能够有效抑制因热声耦合产生的燃烧振荡,为实际动力系统燃烧振荡抑制提供了思路。     

An analytical study of the flame dynamics of a transversely forced asymmetric two-dimensional Bunsen flame

Combustion instabilities in annular combustors are of great interest because of their industrial relevance. Azimuthal acoustic modes, which involve transverse acoustic forcing to flames, have become a key process related to annular combustor instabilities. Transverse mean flow may be a factor that affects azimuthal oscillations. This paper provides an analytical model for a transversely forced two-dimensional Bunsen flame under transverse mean flow. The model is established using a low-amplitude perturbation assumption applied to a G-equation formulation. Forced flame displacement and flame transfer functions (FTFs) are calculated. The results are verified based on numerical solutions of the G-equation. Effects of frequency, transverse mean flow velocity and vertical mean flow velocity on the FTFs are discussed. The symmetric flame without transverse mean flow has a vanishing response to transverse acoustic forcing, while asymmetric flames, which are formed with transverse mean flow, have a bandpass response to transverse forcing. The response at very low and high forcing frequencies is small, with higher transfer function gains only in a certain frequency range. This bandpass response, which is inherently linked to the asymmetry of the flame, is an important factor to account for when considering the flame dynamics related to transverse acoustic effects.     

Impact of symmetry breaking on the Flame Transfer Function of a laminar premixed flame

This work presents a numerical study of the acoustic response of a laminar flame with tunable asymmetry. A V-shaped premixed flame is stabilised in the wake of a cylindrical flame holder that can be rotated. The configuration is symmetric when the flame holder is fixed but increasing its rotation rate breaks the symmetry of the flow. This configuration is submitted to acoustic forcing to measure the effect of rotation of the flame holder on the Flame Transfer Functions. It appears that the asymmetry of the two flame branches changes their respective time delays, resulting in interference in the global unsteady heat release rate fluctuations. Consequently, the Flame Transfer Function exhibits dips and bumps, which are studied via laminar Direct Numerical Simulation. Potential applications for the control of combustion instabilities are discussed.     

Autoignition delay modulation by high-frequency thermoacoustic oscillations in reheat flames

This paper investigates the sensitivity of the autoignition delay in reheat flames to acoustic pulsations associated with high-frequency transverse thermoacoustic oscillations. A reduced order model for the response of purely autoignition-stabilised flames to acoustic disturbances is compared with experimental observations. The experiments identified periodic flame motion associated with high-amplitude transverse limit-cycle oscillations in an atmospheric pressure reheat combustor. This flame motion was assumed to be the result of a superposition of two flame-acoustic coupling mechanisms: autoignition delay modulation by the oscillating acoustic field and displacement and deformation of the flame by the acoustic velocity. The reduced order model coupled to reaction kinetics calculations reveals that a significant portion of the observed flame motion can be attributed to autoignition delay modulation. The ignition position responds instantaneously to the acoustic pressure at the time of ignition, as observed experimentally. The model also provides insight into the importance of the history of acoustic disturbances experienced by the fuel-air mixture prior to ignition. Due to the high-frequency nature of the instability, a fluid particle can experience multiple oscillation cycles before ignition. The ignition delay responds in-phase with the net-acoustic perturbation experienced by a fluid particle between injection and ignition. These findings shed light on the underlying mechanisms of the flame motion observed in experiments and provide useful insight into the importance of autoignition delay modulation as a driving mechanism of high-frequency thermoacoustic instabilities in reheat flames.     

A reduced-order deterministic model describing an intermittency route to combustion instability

Recently, there has been a growing interest in understanding and characterising intermittent burst oscillations that presage the onset of combustion instability. We construct a deterministic model to capture this intermittency route to instability in a bluff-body stabilised combustor by coupling the equations governing vortex shedding and the acoustic wave propagation in a confinement. A feedback mechanism is developed wherein the sound generated due to unsteady combustion affects the vortex shedding. This feedback leads to a variation in the time of impingement of the vortices with the bluff body causing the system to exhibit chaos, intermittency, and limit cycle oscillations. Experimental validation of the model is provided using various precursor measures that quantify the observed intermittent states.     

U-RANS／PDF方法计算钝体后漩涡脱落

本研究发展了U-RANS/PDF混合算法研究湍流和化学反应相互作用对燃烧稳定性的影响,采用有限体积和Monte Carlo相结合的方法在非结构网格中求解相容的U-RANS方程和脉动速度-湍流频率-标量的联合PDF方程.本文对钝体火焰驻定器后冷态流场进行了计算,结果表明此混合算法能够捕捉流场中非稳态的漩涡脱落现象.着重研究了湍流频率模型系数的改变对漩涡脱落频率以及拟序结构在动量输运中的作用的影响.     

Impact of pressure fluctuations on the dynamics of laminar premixed flames

Thermo-acoustic instabilities are problematic in the design of continuous-combustion propulsion systems such as gas turbine engines, rocket motors, jet engine afterburners, and ramjets. Conceptually, the coupling between acoustics and flame dynamics can be divided into two categories: flame area fluctuations and changes in the local flame speed. The latter can be caused by the thermodynamic fluctuations that accompany an acoustic wave. This coupling is the focus of the present work. In this paper, we are concerned with the dynamics of laminar premixed flames involving large hydrocarbon species. Through high-fidelity numerical simulations, we investigate the flame response for a wide range of fuels and acoustic frequencies. The combustion of hydrogen and methane is considered for verification purposes and as baseline cases for comparison with two large hydrocarbon fuels, n-heptane and n-dodecane. We extract the phase and gain of the unsteady heat release response, which are directly related to the Rayleigh criterion and thus the stability of the system. For all fuels, we observe a local peak in the heat release gain. At high frequencies, we find that the fluctuations of the different species mass fractions decrease with the inverse of the acoustic frequency, leading to chemistry being “frozen” in the high-frequency limit. This allows us to predict the flame behavior directly from the steady-state solution.     

Acoustic near-field characteristics of a conical,premixed flame

The occurrence of self-excited pressure oscillations routinely plagues the development of combustion systems. These oscillations are often driven by interactions between the flame and acoustic perturbations. This study was performed to characterize the structure of the acoustic field in the near field of the flame and the manner in which it is influenced by oscillation frequency, combustor geometry, flame length and temperature ratio. The results of these calculations indicate that the acoustic velocity has primarily one- and two-dimensional features near the flame tip and base, respectively. The magnitude of the radial velocity components increases with temperature ratio across the flame, while their axial extent increases with frequency. However, the acoustic pressure has primarily one-dimensional characteristics. They also show that the acoustic field structure exhibits only moderate dependencies upon area expansion and flame temperature ratio for values typical of practical systems. Finally, they show that the local characteristics of the acoustic field, as well as the overall plane-wave reflection coefficient, exhibit a decreasing dependence upon the flame length as the area expansion ratio increases.     

LES studies of the flow in a swirl gas combustor

Environmental and other practical concerns have led to the development of compact gas turbine combustors burning lean mixtures leading to potentially low CO and NO_x emissions. The compact design requires efficient atomization and mixing together with a compact premixed flame. Associated with these requirements are higher temperatures, increased heat transfer, and thermal load, thus increasing the danger of combustion instabilities (causing performance deterioration and excessive mechanical loads), and possible off-design operation. Numerical simulations of reacting flows are well suited to address these issues. To this end, large eddy simulation (LES) is particularly promising. The philosophy behind LES is to explicitly simulate the large scales of the flow and the thermochemistry, affected by boundary conditions whilst modeling only the small scales, including the interaction between the flow and the combustion processes. Here, we examine the flow and the flame in a model gas turbine combustor (General Electric’s lean premixed dry low NO_x LM6000) to evaluate the potential of LES for design studies of engineering applications and to study the effects of the combustor confinement geometry on the flow and on the flame dynamics. Two LES models, a Monotone Integrated LES model with 1 and 2 step Ahrrenius chemistry, and a fractal flame-wrinkling LES model coupled to a conventional one-equation eddy-viscosity subgrid model, are used. Reasonable agreement is found when comparing predictions with experimental data and with other LES computations of the same case. Furthermore, the combustor confinement geometry is found to strongly affect the vortical flow, and hence also the flame and its dynamics.