Nonlinear phenomena in thermoacoustics: A comparison between single-mode and multi-mode methods

Nonlinear analysis of thermoacoustic instability is essential for prediction of frequencies and amplitudes of limit cycles. In frequency domain analyses, a quasi-linear transfer function between acoustic velocity and heat release rate perturbations, called the flame describing function (FDF), is obtained from a flame model or experiments. The FDF is a function of the frequency and amplitude of velocity perturbations but only contains the heat release response at the forcing frequency. While the gain and phase of the FDF provide insight into the nonlinear dynamics of the system, the accuracy of its predictions remains to be verified for different types of nonlinearity. In time domain analyses, the governing equations of the fully coupled problem are solved to find the time evolution of the system. One method is to discretize the governing equations using a suitable basis, such as the natural acoustic modes of the system. The number of modes used in the discretization alters the accuracy of the solution. In our previous work we have shown that predictions using the FDF are almost exactly the same as those obtained from the time-domain using only one mode for the discretization. We call this the single-mode method. In this paper we compare results from the single-mode and multi-mode methods, applied to a thermoacoustic system of a premixed flame in a tube. For some cases, the results differ greatly in both amplitude as well as frequency content. This study shows that the contribution from higher and subharmonics to the nonlinear dynamics can be significant and must be considered for an accurate and comprehensive analysis of thermoacoustic systems. Hence multi-mode simulations are necessary, and the single-mode method or the FDF may be insufficient to capture some of the complex nonlinear behaviour in fhermoacoustics.

Experimental investigation of dynamic response of acoustically forced turbulent premixed CH4/CO2/air flames

Increasing demand for energy and continuing increase in environmental as well as financial cost of use of fossil fuels drive the need for utilization of fuels from sustainable sources for power generation. Development of fuel-flexible combustion systems is vital in enabling the use of sustainable fuels. It is also important that these sustainable combustion systems meet the strict governmental emission legislations. Biogas is considered as one of the viable sustainable fuels that can be used to power modern gas turbines: However, the change in chemical, thermal and transport properties as well as change in Wobbe index due to the variation of the fuel constituents can have a significant effect on the performance of the combustor. It is known that the fuel properties have strong influence on the dynamic flame response; however there is a lack of detailed information regarding the effect of fuel compositions on the sensitivity of the flames subjected to flow perturbations. In this study, we describe an experimental effort investigating the response of premixed biogas-air turbulent flames with varying proportions of CH4 and CO2 to velocity perturbations. The flame was stabilized using a centrally placed conical bluff body. Acoustic perturbations were imposed to the flow using loud speakers. The flame dynamics and the local heat release rate of these acoustically excited biogas flames were studied using simultaneous measurements of OH and H2CO planar laser induced fluorescence. OH* chemiluminescence along with acoustic pressure measurements were also recorded to estimate the total flame heat release modulation and the velocity fluctuations. The measurements were carried out by keeping the theoretical laminar flame speed constant while varying the bulk velocity and the fuel composition. The results indicate that the flame sensitivity to perturbations increased with increased dilution of CH4 by CO2 at low amplitude forcing, while at high amplitude forcing conditions the magnitude of the flame response was independent of dilution.

Complex chemistry DNS of n-heptane spray autoignition at high pressure and intermediate temperature conditions

Direct Numerical Simulations (DNS) of turbulent n-heptane sprays autoigniting at high pressure (P=24bar) and intermediate air temperature (Tair=1000K) have been performed to investigate the physical mechanisms present under conditions where low-temperature chemistry is expected to be important. The initial turbulence in the carrier gas, the global equivalence ratio in the spray region, and the initial droplet size distribution of the spray were varied. Results show that spray ignition exhibits a spotty nature, with several kernels developing independently in those regions where the mixture fraction is close to its most reactive value ξMR (as determined from homogeneous reactor calculations) and the scalar dissipation rate is low. Turbulence reduces the ignition delay time as it promotes mixing between air and the fuel vapor, eventually resulting in lower values of scalar dissipation. High values of the global equivalence ratio are responsible for a larger number of ignition kernels, due to the higher probability of finding regions where ξ=ξMR. Spray polydispersity results in the occurrence of ignition over a wider range of mixture fraction values. This is a consequence of the inhomogeneities in the mixing field that characterize these sprays, where poorly mixed rich spots are seen to alternate with leaner ones which are well-mixed. The DNS simulations presented in this work have also been used to assess the applicability of the Conditional Moment Closure (CMC) method to the simulation of spray combustion. CMC is found to be a valid method for capturing spray autoignition, although care should be taken in the modelling of the unclosed terms appearing in the CMC equations. © 2013 The Combustion Institute.

Flow field measurements of pulverized coal combustion using optical diagnostic techniques

This paper demonstrates the application of laser Doppler velocimetry (LDV) and particle image velocimetry (PIV) techniques to a particle-laden reacting flow of pulverized coal. A laboratory-scale open-type annular burner is utilized to generate velocity profiles of coal particles and micrometric alumina particles. Pair-wise two-component LDV measurements and high-speed stereo PIV measurements provide three-dimensional velocity components of the flow field. A detailed comparison of velocities for alumina and coal particle seeding revealed differences attributed to the wide size distribution of coal particles. In addition, the non-spherical shape and high flame luminosity associated with coal particle combustion introduces noise to the Mie scatter images. The comparison of mean and RMS velocities measured by LDV and PIV techniques showed that PIV measurements are affected by the wide size distribution of coal particles, whereas LDV measurements become biased toward the velocity of small particles, as signals from large particles are rejected. This small-particle bias is also reflected in the spectral characteristics for both techniques, which are in good agreement within the range of frequencies accessible. PIV measurements showed an expected lack of response of large coal particles to the turbulence fluctuations. The overall good agreement between LDV and PIV measurements demonstrates the applicability of the high-speed PIV technique to a particle-laden, high luminosity coal flame while highlighting some of its limitations. © 2013 Springer-Verlag Berlin Heidelberg.

Time resolved measurement of cold start hc concentration using the fast FID

Understanding mixture formation phenomena during the first few cycles of an engine cold start is extremely important for achieving the minimum engine-out emission levels at the time when the catalytic converter is not yet operational. Of special importance is the structure of the charge (film, droplets and vapour) which enters the cylinder during this time interval as well as its concentration profile. However, direct experimental studies of the fuel behaviour in the inlet port have so far been less than fully successful due to the brevity of the process and lack of a suitable experimental technique. We present measurements of the hydrocarbon (HC) concentration in the manifold and port of a production SI engine using the Fast Response Flame Ionisation Detector (FRFID). It has been widely reported in the past few years how the FRFID can be used to study the exhaust and in-cylinder HC concentrations with a time resolution of a few degrees of crank angle, and the device has contributed significantly to the understanding of unburned HC emissions. Using the FRFID in the inlet manifold is difficult because of the presence of liquid droplets, and the low and fluctuating pressure levels, which leads to significant changes in the response time of the instrument. However, using recently developed procedures to correct for the errors caused by these effects, the concentration at the sampling point can be reconstructed to align the FRFID signal with actual events in the engine. © 1996 Society of Automotive Engineers, Inc.

Effect of operating conditions and fuel type on crevice hc emissions: Model results and comparison with experiments

A one-dimensional model for crevice HC post-flame oxidation is used to calculate and understand the effect of operating parameters and fuel type (propane and isooctane) on the extent of crevice hydrocarbon and the product distribution in the post flame environment. The calculations show that the main parameters controlling oxidation are: bulk burned gas temperatures, wall temperatures, turbulent diffusivity, and fuel oxidation rates. Calculated extents of oxidation agree well with experimental values, and the sensitivities to operating conditions (wall temperatures, equivalence ratio, fuel type) are reasonably well captured. Whereas the bulk gas temperatures largely determine the extent of oxidation, the hydrocarbon product distribution is not very much affected by the burned gas temperatures, but mostly by diffusion rates. Uncertainties in both turbulent diffusion rates as well as in mechanisms are an important factor limiting the predictive capabilities of the model. However, it seems well suited to sensitivity calculations about a baseline. Copyright © 1999 Society of Automotive Engineers, Inc.

Development of a time and space resolved sampling probe diagnostic for engine exhaust hydrocarbons

In order to understand how unburned hydrocarbons emerge from SI engines and, in particular, how non-fuel hydrocarbons are formed and oxidized, a new gas sampling technique has been developed. A sampling unit, based on a combination of techniques used in the Fast Flame Ionization Detector (FFID) and wall-mounted sampling valves, was designed and built to capture a sample of exhaust gas during a specific period of the exhaust process and from a specific location within the exhaust port. The sampling unit consists of a transfer tube with one end in the exhaust port and the other connected to a three-way valve that leads, on one side, to a FFID and, on the other, to a vacuum chamber with a high-speed solenoid valve. Exhaust gas, drawn by the pressure drop into the vacuum chamber, impinges on the face of the solenoid valve and flows radially outward. Once per cycle during a specified crank angle interval, the solenoid valve opens and traps exhaust gas in a storage unit, from which gas chromatography (GC) measurements are made. The port end of the transfer tube can be moved to different locations longitudinally or radially, thus allowing spatial resolution and capturing any concentration differences between port walls and the center of the flow stream. Further, the solenoid valve's opening and closing times can be adjusted to allow sampling over a window as small as 0.6 ms during any portion of the cycle, allowing resolution of a crank angle interval as small as 15°CA. Cycle averaged total HC concentration measured by the FFID and that measured by the sampling unit are in good agreement, while the sampling unit goes one step further than the FFID by providing species concentrations. Comparison with previous measurements using wall-mounted sampling valves suggests that this sampling unit is fully capable of providing species concentration information as a function of air/fuel ratio, load, and engine speed at specific crank angles. © Copyright 1996 Society of Automotive Engineers, Inc.

Initial phase modelling in numerical explosion applied to process safety

The utilisation of computational fluid dynamics (CFD) in process safety has increased significantly in recent years. The modelling of accidental explosion via CFD has in many cases replaced the classical Multi Energy and Brake Strehlow methods. The benefits obtained with CFD modelling can be diminished if proper modelling of the initial phase of explosion is neglected. In the early stages of an explosion, the flame propagates in a quasi-laminar regime. Proper modelling of the initial laminar phase is a key aspect in order to predict the peak pressure and the time to peak pressure. The present work suggests a modelling approach for the initial laminar phase in explosion scenarios. Findings are compared with experimental data for two classical explosion test cases which resemble the common features in chemical process areas (confinement and congestion). A detailed analysis of the threshold for the transition from laminar to turbulent regime is also carried out. The modelling is implemented in a fully 3D Navier-Stokes compressible formulation. Combustion is treated using a laminar flamelet approach based on the Bray, Moss and Libby (BML) formulation. A novel modified porosity approach developed for the unstructured solver is also considered. Results agree satisfactorily with experiments and the modelling is found to be robust. © 2013 The Institution of Chemical Engineers.

Analysis of high pressure premixed flames using Equivalent Reactor Networks for predicting NOx emissions

This paper describes a computational study of lean premixed high pressure methane-air flames, using Computational Fluid Dynamics (CFD) together with a reactor network approach. A detailed chemical reaction mechanism is employed to predict pollutant concentrations, placing emphasis on nitrogen oxide emissions. The reacting flow field is divided into separate zones in which homogeneity of the physical and chemical conditions prevails. The defined zones are interconnected forming an Equivalent Reactor Network (ERN). Three flames are examined for which experimental data is available. Flame A is characterised by an equivalence ratio of 0.43 while Flames B and C are richer with equivalence ratios of 0.5 and 0.56 respectively. Computations are performed for a range of operating conditions, quantifying the effect in the emitted NOx levels. Model predictions are compared against the available experimental data. Sensitivity analysis is performed to investigate the effect of the network size, in order to define the optimum number of reactors for accurate predictions of the species mass fractions. © 2012 Elsevier Ltd. All rights reserved.

Low-order modelling of the response of ducted flames in annular geometries

The adoption of lean premixed prevaporised combustion systems can reduce NOx emissions from gas turbines, but unfortunately also increases their susceptibility to thermoacoustic instabilities. Initially, acoustic waves can produce heat release fluctuations by a variety of mechanisms, often by perturbing the equivalence ratio. If correctly phased, heat release fluctuations can subsequently generate more acoustic waves, which at high amplitude can result in significant structural damage to the combustor. The prediction of this phenomenon is of great industrial interest. In previous work, we have coupled a physics based, kinematic model of the flame with a network model to provide the planar acoustic response necessary to close the feedback loop and predict the onset and amplitude of thermoacoustic instabilities in a lab-scale, axisymmetric single burner combustor. The advantage of a time domain approach is that the modal interaction, the influence of harmonics, and flame saturation can be investigated. This paper extends this approach to more realistic, annular geometries, where both planar and circumferential modes must be considered. In lean premixed prevaporised combustors, fluctuations in equivalence ratio have been shown to be a dominant cause of unsteady combustion. These can occur, for example, due to velocity perturbations in the premix ducts, which can lead to equivalence ratio fluctuations at the fuel injectors, which are subsequently convected downstream to the flame surfaces. Here, they can perturb the heat release by locally altering the flame speed, enthalpy of combustion, and, indirectly, the flame surface area. In many gas turbine designs, particularly aeroengines, the geometries are composed of a ring of premix ducts linking a plenum and an annular combustor. The most unstable modes are often circumferential modes. The network model is used to characterise the flow response of the geometry to heat fluctuations at an appropriate location, such as the fuel injectors. The heat release at each flame holder is determined in the time domain using the kinematic flame model derived, as a function of the flow perturbations in the premix duct. This approach is demonstrated for an annular ring of burners on a in a simple geometry. The approach is then extended to an industrial type gas turbine combustor, and used to predict the limit cycle amplitudes. Copyright © 2012 by ASME.

Nonlinear phenomena in thermoacoustic systems with premixed flames

Nonlinear analysis of thermoacoustic instability is essential for prediction of frequencies, amplitudes and stability of limit cycles. Limit cycles in thermoacoustic systems are reached when the energy input from driving processes and energy losses from damping processes balance each other over a cycle of the oscillation. In this paper an integral relation for the rate of change of energy of a thermoacoustic system is derived. This relation is analogous to the well-known Rayleigh criterion in thermoacoustics, but can be used to calculate the amplitudes of limit cycles, as well as their stability. The relation is applied to a thermoacoustic system of a ducted slot-stabilized 2-D premixed flame. The flame is modelled using a nonlinear kinematic model based on the G-equation, while the acoustics of planar waves in the tube are governed by linearised momentum and energy equations. Using open-loop forced simulations, the flame describing function (FDF) is calculated. The gain and phase information from the FDF is used with the integral relation to construct a cyclic integral rate of change of energy (CIRCE) diagram that indicates the amplitude and stability of limit cycles. This diagram is also used to identify the types of bifurcation the system exhibits and to find the minimum amplitude of excitation needed to reach a stable limit cycle from another linearly stable state, for single- mode thermoacoustic systems. Furthermore, this diagram shows precisely how the choice of velocity model and the amplitudedependence of the gain and the phase of the FDF influence the nonlinear dynamics of the system. Time domain simulations of the coupled thermoacoustic system are performed with a Galerkin discretization for acoustic pressure and velocity. Limit cycle calculations using a single mode, as well as twenty modes, are compared against predictions from the CIRCE diagram. For the single mode system, the time domain calculations agree well with the frequency domain predictions. The heat release rate is highly nonlinear but, because there is only a single acoustic mode, this does not affect the limit cycle amplitude. For the twenty-mode system, however, the higher harmonics of the heat release rate and acoustic velocity interact resulting in a larger limit cycle amplitude. Multimode simulations show that in some situations the contribution from higher harmonics to the nonlinear dynamics can be significant and must be considered for an accurate and comprehensive analysis of thermoacoustic systems. Copyright © 2012 by ASME.

Obtaining bifurcation diagrams with a thermoacoustic network model

Linear techniques can predict whether the non-oscillating (steady) state of a thermoacoustic system is stable or unstable. With a sufficiently large impulse, however, a thermoacoustic system can reach a stable oscillating state even when the steady state is also stable. A nonlinear analysis is required to predict the existence of this oscillating state. Continuation methods are often used for this but they are computationally expensive. In this paper, an acoustic network code called LOTAN is used to obtain the steady and the oscillating solutions for a horizontal Rijke tube. The heat release is modelled as a nonlinear function of the mass flow rate. Several test cases from the literature are analysed in order to investigate the effect of various nonlinear terms in the flame model. The results agree well with the literature, showing that LOTAN can be used to map the steady and oscillating solutions as a function of the control parameters. Furthermore, the nature of the bifurcation between steady and oscillating states can be predicted directly from the nonlinear terms inside the flame model. Copyright © 2012 by ASME.

On Mixing and Shaking : Structure and Dynamics of Turbulent Stratified Flames

The drive for low emission combustion systems encourages applications using premixed flames. Yet in many applications, considerations of flame stability or mixing times lead to systems with neither premixed nor diffusion flames, which are often called technically premixed or stratified flames. In this talk we discuss the current state of understanding of the effect of mixing and extent of stratification on the structure, microstructure and dynamics of selected turbulent stratified flames. Over the past few years, a significant database of scalar and velocity data has been built to analyze the effects of unmixedness on local and global flame structure. Microscale studies of the flame structures show in detail how the effect of local stratification affects (or not!) the flame structure, flame surface density and scalar dissipation rates, and production of selected species. The experiments place exacting demands on current spectroscopic diagnostics, and reveal the progress and limits to our understanding of turbulent flames in general. The dynamics of stratified flames with respect to instabilities is also shown to be very rich, as the particular shape of the flames and the stabilization points are is significantly affected by the fuel distribution, modifying the rate and location of heat release, and thus the coupling with the surrounding acoustics and determining the onset of self-excitations.