Response of two acoustically excited turbulent premixed flames to an imposed phase lag

This paper describes an experimental investigation into the interactions that occur between two acoustically forced lean turbulent premixed flames for an induced phase lag. Phase-averaged FSD from cinematographic OH-PLIF measurements and global heat release measurements were obtained for a range phase lags (ψ S) and amplitudes (A) as a function of flame separation distance, S. The effect of bringing two flames closer together causes jet merging, which alters the vortex flame interactions that drive the thermo-acoustic response. To simulate circumferential modes a phase lag was introduced, which affected the flame dynamics in the region of flame-flame interaction, with these changes dependent on S. For moderate separation distances, the flame structure becomes increasingly asymmetric inducing a very small transverse oscillation. However, for moderate phase lags (φ s ≤ 20) the magnitude of these changes and their subsequent influence on the thermo-acoustic response was found to be slight in comparison with changes in S. Copyright © 2012 by the American Institute of Aeronautics and Astronautics, Inc.

The behaviour of laminar stratified methane/air flames in counterflow

Flames propagating through a mixture with a gradient of equivalence ratio have been previously demonstrated to travel faster or slower than their equivalent premixed flames. The present study aims to numerically investigate the response of strained laminar methane-air flames to such gradients. The flames are simulated in a counterflow configuration where a premixed reactant stream at equivalence ratio φR opposes a hot equilibrium stream at equivalence ratio φP. Premixed and stratified flames are compared with respect to the equivalence ratio φ* and the corresponding gradient ∇φ* at the point of peak heat release rate, for three strain rates, a=50, 300 and 500s-1 and a range of φ*. The effect of different stratification levels is also investigated by varying the ratio of φP to φR, Θ. Results indicate that, as long as flames stabilize within the diffusion layer and Θ>1, increased heat release rate Q is seen throughout the progress variable space in comparison to the premixed state. In contrast, an attenuation of heat release rate is seen for Θ<1. The enhancement (or attenuation) of heat release varies monotonically with Θ. The effect of stratification on flame behavior becomes more pronounced as the strain rate increases. The present study reveals the mechanisms for the propagation of quasi-steady stratified flames under lean and rich conditions: stratified flames are primarily dominated by the diffusion of heat under lean conditions, and diffusion of H2 under rich conditions. Thanks to species and thermal support, stratified flames continue to burn beyond the premixed lean and rich flammability limits. Further investigation on the unsteady response of flames to the fluctuating equivalence ratio implies that the steady results represent the unsteady response well, as long as φ* and ∇φ* are similar in both steady and unsteady cases. © 2013 The Combustion Institute.

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.