Time-resolved measurements of hydroxyl in stable and extinguishing partially premixed turbulent opposed-jet flames

Quantitative hydroxyl time-series measurements from a set of stable and extinguishing turbulent opposed-flow partially premixed CH₄/air flames are used to investigate the effect of Reynolds number and fuel-side equivalence ratio on the structure of turbulent partially premixed flames. The hydroxyl (OH) integral time scale, computed from the autocorrelation function, is used to characterize OH fluctuations and is found to reach a minimum at the axial location of peak OH. Analyses of the duration of and period between bursts in the OH time series are used to examine the dynamics of flame-front motion. In general, with increasing Reynolds number (Re), the distribution in OH burst times shifts towards smaller time scales. A hydroxyl intermittency parameter is also defined from the bursts to quantify the presence or absence of OH. For flames with the same fuel-side equivalence ratio, the hydroxyl intermittency at peak OH remains almost constant when going from stable to extinguishing flames. However, histograms portray an increase in burst separation times for flames displaying occasional extinction events. Hydroxyl time series for a partially premixed flame at a fuel-side equivalence ratio of 2.0 and Re = 6650 are synthesized by using mixture-fraction simulations based on calculated state relationships for OH versus mixture fraction (f). The laminar-flamelet model is employed to explore relations between OH and f so as to predict trends in mixture-fraction time scales.“Time-Series Measurements in Turbulent Opposed-Jet Flames" is submitted for consideration as a full length article to Flow Turbulence and Combustion.     

An experimental investigation of the interaction of swirl flow with partially premixed disk stabilized propane flames

The present work describes the experimental investigation of reacting wakes established through fuel injection and staged premixing with air in an axisymmetric double cavity arrangement, formed along three concentric disks, and stabilized in the downstream vortex region of the afterbody. The burner assembly is operated with a co-flow of swirling air, aerodynamically introduced upstream of the burner exit plane, to allow for the study of the interaction between the resulting partially premixed recirculating afterbody flames with the surrounding swirl. At low swirl the primary afterbody disk stabilizes the partially premixed annular jet in the downstream reacting wake formation region. As swirl increases, a system of two successive vortices emerges along the axis of the developing wake; the primary afterbody vortex is cooperating with an adjacent, swirl induced, central recirculation zone and this combination further promotes turbulent mixing in the hot wake.Complementary measurements of the counterpart isothermal turbulent velocity fields provided important information on the near wake aerodynamics under the interaction of the variable swirl and the double cavity produced annular jet stabilized by the afterbody. Under reacting conditions, measurements of turbulent velocities, temperatures and statistics together with an evaluation of the exhaust emissions were performed using LDV, thin digitally-compensated thermocouples and gas analyzers. A selected number of lean and ultra-lean flames were investigated by regulating the injected fuel and the air supply ratio, while the influence of the variation of the imposed swirl on wake development, flame characteristics and emission performance was documented for constant fuel injections. The differences and similarities between the present partially premixed stabilizer and other types of axisymmetric configurations are also highlighted and discussed.     

Simultaneous OH and Schlieren visualization of premixed flames at the lean blow-out limit

The structure of an air-propane premixed flame was studied experimentally at the lean flammability limit, using Schlieren photography synchronized with OH-imaging done with the Planar Laser Induced Fluorescence (PLIF) technique. The flame was studied in a wide range of fuel equivalence ratios. Various steps in the process of the flame destabilization were investigated, including partial lift-off, stable lift-off, and final blow-out conditions. The flame structure was visualized for each stage showing the transition from a flame held at the nozzle to a flame held by the flow structures. In order to study the latter conditions in more detail the flame was acoustically excited at the preferred mode frequency generating large, stable, coherent structures in the core region. The modified flame structure was visualized to understand the interaction between the flame and vortical flow dynamics.It is shown that for the flow conditions when the flame cannot be stabilized at the nozzle, a new anchoring point is reached at the location of the initial vortex roll-up in the jet shear layer. At this point the flow reversal and transition to turbulence produce stagnation points with relatively low local velocities and velocity gradients where the flame can be stabilized. When the flame jet is being forced at the jet most unstable frequency, large coherent structures are formed and the flame is stabilized intermittently on these vortices.     

Imaging of premixed flames in microgravity

A laser schlieren system which uses video recording and digital images analysis has been developed and applied successfully to microgravity combustion experiments performed in a drop-tower. The optical system and the experiment are installed within a small package which is subjected to free-fall. The images are recorded on video tape and are digitized and analyzed by a computer-controlled image processor. The experimental results include laminar and turbulent premixed conical flames in microgravity, normal positive gravity (upward), and reverse gravity (downward). The procedures to extract frequency information from the digitized images are described. Many gross features of the effects of gravity on premixed conical flames are found. Flames that ignite easily in normal gravity fail to ignite in microgravity. Buoyancy driven instabilities associated with an interface formed between the hot products and the cold surrounding air is the mechanism through which gravity influences premixed laminar and turbulent flames. In normal gravity, this causes the flame to flicker. In reverse gravity, -g, and microgravity, g, the interface is stable and flame flickering ceases. The flickering frequencies of +g flames vary with changing upstream boundary conditions. The absence of flame flickering in g suggest that g flames would be less sensitive to these changes.This work is supported by NASA Microgravity Sciences and Applications Divisions under contract No. C-32000-R through the U.S. Department of Energy Contract No. DE-AC03-76F00098. Technical support is provided by NASA Lewis Research Center. Project Scientist is Dr. Karen J. Weiland. The authors would like to acknowledge Dr. Liming Zhou for his contribution to early testing of the schlieren system, and to Mr. Gray Hubbard for writing the image analysis software     

Structures and stabilization of low calorific value gas turbulent partially premixed flames in a conical burner

Experiments are carried out on partially premixed turbulent flames stabilized in a conical burner. The investigated gaseous fuels are methane, methane diluted with nitrogen, and mixtures of CH₄, CO, CO₂, H₂ and N₂, simulating typical products from gasification of biomass, and co-firing of gasification gas with methane. The fuel and air are partially premixed in concentric tubes. Flame stabilization behavior is investigated and significantly different stabilization characteristics are observed in flames with and without the cone. Planar laser induced fluorescence (LIF) imaging of a fuel-tracer species, acetone, and OH radicals is carried out to characterize the flame structures. Large eddy simulations of the conical flames are carried out to gain further understanding of the flame/flow interaction in the cone. The data show that the flames with the cone are more stable than those without the cone. Without the cone (i.e. jet burner) the critical jet velocities for blowoff and liftoff of biomass derived gases are higher than that for methane/nitrogen mixture with the same heating values, indicating the enhanced flame stabilization by hydrogen in the mixture. With the cone the stability of flames is not sensitive to the compositions of the fuels, owing to the different flame stabilization mechanism in the conical flames than that in the jet flames. From the PLIF images it is shown that in the conical burner, the flame is stabilized by the cone at nearly the same position for different fuels. From large eddy simulations, the flames are shown to be controlled by the recirculation flows inside cone, which depends on the cone angle, but less sensitive to the fuel compositions and flow speed. The flames tend to be hold in the recirculation zones even at very high flow speed. Flame blowoff occurs when significant local extinction in the main body of the flame appears at high turbulence intensities.     

Tulip flames: changes in shape of premixed flames propagating in closed tubes

 The experimental results that are the subject of this communication provide high-speed schlieren images of the closed-tube flame shape that has come to be known as the tulip flame. The schlieren images, along with in-chamber pressure records, help demonstrate the effects of chamber length, equivalence ratio, and igniter geometry on formation of the tulip flame. The pressure/time records show distinct features which correlate with flame shape changes during the transition to tulip. The measurements indicate that the basic tulip flame formation is a robust phenomenon that depends on little except the overall geometry of the combustion vessel. Received: 22 April 1997/Accepted: 7 July 1997     

Dynamics of PIV seeding particles in turbulent premixed flames

Particle image velocimetry (PIV) estimates the fluid velocity field measuring the displacement of small dispersed particles between two successive instants separated by a small time interval. The accuracy of the measurements depends on the ability of the particles to accommodate their velocity to the fluid fluctuations. When the fluid is subjected to extreme accelerations, the small but finite inertia prevents the particles from following the fluid, originating a substantial relative velocity. This effect is shown to be crucial for applications of PIV to turbulent premixed combustion, particularly in the product region at locations just behind the instantaneous flame front. The issuing inaccuracy may easily spoil the estimate of certain statistical observables which are of crucial importance in the theory of turbulent premixed combustion. By exploiting the direct numerical simulation of a model air/methane flame, a suitable criterion for proper particle seeding is validated and compared with the corresponding experiments with a combined PIV/OH-LIF (laser-induced fluorescence) system. The proposed parameter, the flamelet Stokes number, depends on particle properties and thermochemical conditions of the flame and substantially restricts the particle dimensions required for a reliable estimate of the relevant flow statistics.     

Correlation coefficients of the fluctuating density in turbulent premixed flames

 Two-point density measurements by laser induced Rayleigh scattering are used in this study to fully characterise the scalar field in a Bunsen type turbulent premixed flame. The two points are separated within the flame brush in the axial or radial directions. Correlation coefficients are obtained by comparing the evolution of one-point density fluctuations in time or the two-point density fluctuations in both space and time. Time and length scales of the scalar field, and the mean convection velocity of the turbulent scalar structures are deduced from these correlation coefficients. Time scales are calculated from the auto-correlation coefficients, length scales are determined from the space correlation coefficients and the mean convection velocity of the scalar structures in the axial direction is deduced from the space–time correlation coefficients. The relevance of these results for analysing and modelling the structure of turbulent premixed flames is discussed. Received: 30 April 1996 / Accepted: 2 September 1997     

Experimental investigation of dynamics of premixed acetylene–air flames in a micro-combustor

This paper describes a detailed experimental study performed to investigate the flame propagation behaviour of premixed flames in micro-channels. A novel, modular, stackable micro-combustor was developed for this purpose. For a chosen planar channel geometry, the flow condition and the mixture equivalence ratio of premixed acetylene–air were varied to investigate various modes of operation. Three different modes of operation were observed; they were (i) stable periodic operation – consisting of ignition, flame propagation, flame extinction, and re-ignition, (ii) a-periodic operation, and (iii) anchored flame condition. The present work also aims to provide quantitative information on the dynamics of premixed acetylene–air flames propagating inside micro-channels. A novel measurement approach based on OH^* chemiluminescence measurements employing a single photomultiplier unit was developed for this purpose. The data recorded were post processed using an in-house developed MATLAB code to evaluate the mean flame propagation speed measured between three different spatial locations along the length of the micro-channel. The results from the flame propagation speed measurements performed during ‘periodic’ mode of operation indicated that the flame travelled at higher propagation speed in the mid-length region of the channel compared to that at the initial entry point, suggesting flame acceleration. This flame acceleration could be attributed to a situation where the flame experienced different local equivalence ratio conditions at different upstream locations. The results suggest that after completion of a cycle of operation consisting of ignition, flame propagation and flame extinction, the fresh mixture that filled the channel was diluted with the exhaust gas from the previous cycle. This pocket of diluted mixture convected downstream with time, thus enabling the spatial variation in local equivalence ratio along the micro-channel.     

Experimental analysis of flashback in lean premixed swirling flames: upstream flame propagation

Simultaneous 10-kHz OH-PLIF and 20-kHz two-component PIV were made in conjunction with wide-field 20-kHz flame luminescence imaging of an unconfined, swirling, lean premixed, bluff-body stabilized flame during flashback. Flashback was induced by increasing the stoichiometry or swirl number or reducing the Reynolds number. A detailed stability regime was prepared and compared to predictions. Analysis of the time-correlated flame history inside the exit nozzle during flashback and non-flashback flame events led to a new hypothesis for the flashback mechanism.     

Numerical study on the instability of premixed plane flames in the three-dimensional field

The instability of premixed plane flames in the three-dimensional (3-D) field is investigated by means of the numerical simulation. We show numerically that infinitesimal disturbances superimposed on the flames grow exponentially with time, as predicted in the linear analysis, and obtain the growth rates of disturbances depending on the absolute values of the wave-number vectors. The growth rates of the 3-D flames are consistent with those of the two-dimensional (2-D) flames. The hydrodynamic effect has a destabilizing influence on the instability of flames, and the diffusive-thermal effect has a destabilizing/stabilizing influence for Le < 1/Le > 1. Moreover, we produce the hexagonal cellular structure of the flame front not only for Le < 1 but also for Le > 1, where the spacing between cells in flames for Le < 1 is small compared to that for Le > 1. The spacing of the 3-D flames is 2/√3 times as long as the cell size of the 2-D flames.     

A novel method for chemistry tabulation of strained premixed/stratified flames based on principal component analysis

The principal component analysis (PCA) is used to analyze the high dimensional chemistry data of laminar premixed/stratified flames under strain effects. The first few principal components (PCs), which have larger contribution ratios, are chosen as the tabulated scalars to build the look-up chemistry table. Prior tests show that the strained premixed flame structure can be well reconstructed. To highlight the physical meanings of the tabulated scalars in the stratified flames, a modified PCA method is developed, where the mixture fraction is used to replace one of the PCs with the highest correlation coefficient with a mixture fraction. The other two tabulated scalars are then modified with the Schmidt orthogonalization. The modified tabulated scalars not only have clear physical meanings, but also contain passive scalars. The PCA method has good commonality, and can be extended for building the thermo-chemistry table including strain rate effects when different fuels are used.     

Experimental analysis of flashback in lean premixed swirling flames: conditions close to flashback

Swirling lean premixed flames are of practical relevance due to their potential for low nitric oxide (NOx) emissions. Unfortunately, these flames have various drawbacks. One critical attribute is the possibility for flashback of the reacting flow into the nozzle. Advanced numerical simulations should be able in the future to predict the transition from stable flames to flashback. For a better understanding of the process itself and for validation of numerical simulation a well-documented generic benchmark experiment is needed. This study presents a burner configuration that has already been studied extensively in the past. By minor geometrical adaptations, and via the possibility to vary the swirl intensity in a controlled way, the transition from stable flames to flashback is now accessible to detailed characterisation using advanced laser diagnostics. In a first part of this study the different states of the flame were classified. In the second part, both a stable and a precessing flame very close to flash back were compared in terms of flow and scalar field. The variation of the swirl intensity on the flame is discussed. Because the flame is strongly influenced by its inflow conditions additional velocity measurements inside the nozzle were carried out. This is of special importance for subsequent numerical simulations to match the experimental conditions. The quantitative investigation of the flame during flashback is subjected to consecutive experiments where planar laser diagnostics at high repetition rates will be exploited.     

Detailed measurements of equivalence ratio modulations in premixed flames using laser Rayleigh scattering and absorption spectroscopy

Equivalence ratio non-uniformities may give rise to some of the instabilities observed in modern lean premixed combustion systems. The present work intends to investigate the influence of equivalence ratio perturbations on the dynamics of premixed flames. A burner equipped with a secondary injection system is used to generate equivalence ratio perturbations which are convected by the flow and impinge on a conical flame. Two laser-diagnostics, based on Rayleigh scattering and hydrocarbon infrared absorption, respectively, are employed to give insight into the spatial and temporal evolution of the mixture composition field. Rayleigh scattering images also reveal the flame front dynamics providing an indication on the response of a weakly turbulent flame subject to mixture composition inhomogeneities. Laser light absorption provides a time resolved signal which is used to estimate the equivalence ratio perturbation level. A theoretical model based on the G-equation is used to interpret the experimental data and compare the relative effects of velocity and equivalence ratio perturbations.     

Measurement of propagation speeds in adiabatic cellular premixed flames of CH4 + O2 + CO2

Experimental measurements of the propagation speed of adiabatic flames of methane + oxygen + carbon dioxide are presented. The oxygen content O₂/(O₂ + CO₂) in the artificial air was 31.55% and 35%. Non-stretched flames were stabilized on a perforated plate burner at atmospheric pressure. A heat flux method was used to determine propagation speeds under conditions when the net heat loss of the flame is zero. Under specific experimental conditions the flames become cellular; this leads to significant modification of the flame propagation speed. The onset of cellularity was observed throughout the stoichiometric range of the mixtures studied. Measurements in cellular flames are presented and compared with those for laminar flat flames. Cellularity disappeared when the flames became only slightly sub-adiabatic. Visual and photographic observations of the flames were performed to quantify their cellular structure. Increasing the oxygen content in the artificial air and increasing the temperature of the burner plate led to increase of the number of cells observed.     

Determination of limit cycles for a modified van der Pol oscillator

The determination of amplitude and period of limit cycles is a crucial question in non-linear mechanics. In this paper, a van der Pol oscillator containing a periodic potential is considered. Amplitude and period of limit cycles are calculated by He’s variational method and Krylov–Bogoliubov–Mitropolsky (KBM) method.     

A-priori testing of an eddy viscosity model for the density-weighted subgrid scale stress tensor in turbulent premixed flames

In this study, we report on the direct measurement of the density-weighted subgrid scale (SGS) stress tensor in turbulent premixed flames. In large-eddy simulations (LES), this unresolved tensor is typically modelled using eddy viscosity approaches. Additionally to the direct measurement, we provide a pure experimentally based a-priori test of the commonly used eddy viscosity model suggested by Smagorinsky. For two turbulent premixed V-shaped methane–air flames, a statistical analysis is presented where the correlation between the directly measured SGS stress tensor and the eddy viscosity model following Smagorinsky is tested. The measurement strategy is based on the application of a dual-plane stereo-PIV technique which enables the measurement of the 3D flow field in two parallel planes. This allows the determination of velocities as well as velocity gradients in all three directions. Here, a vector resolution of 118 μm was achieved. For a priori testing, the data are subjected to a spatial filtering procedure that reproduces the application of the filter function in LES. The calculation of velocity gradients is performed after the application of this spatial averaging. Additionally to the velocity field, the flame front position is deduced from the clearly observable step in the tracer particle number density between burnt and unburnt regions of the flame. This facilitates the direct single-shot-based evaluation of all components of the density-weighted SGS stress tensor. Additionally, the model expressions related to these terms can be determined, which is done in this first study for the static Smagorinsky model. With that, the instantaneous local comparison between directly measured stress terms and modelled terms is possible, based on the instantaneous local evaluation procedure. The measurement procedure is described, and first results are presented and discussed. They show a rather poor performance of the static form of the Smagorinsky model (with fixed Smagorinsky constant). Our future aims are to use the directly measured SGS data for the a-priori comparison with more advanced models.     

Effects of operating pressure on flame oscillation and emission characteristics in a partially premixed swirl combustor

The influence of varying combustor pressure on flame oscillation and emission characteristics in the partially premixed turbulent flame were investigated. In order to investigate combustion characteristics in the partially premixed turbulent flame, the combustor pressure was controlled in the range of −30 to 30 kPa for each equivalence ratio (Φ = 0.8-1.2). The r.m.s. of the pressure fluctuations increased with decreasing combustor pressure for the lean condition. The combustor pressure had a sizeable influence on combustion oscillation, whose dominant frequency varied with the combustor pressure. Combustion instabilities could be controlled by increasing the turbulent intensity of the unburned mixture under the lean condition. An unstable flame was caused by incomplete combustion; hence, EICO greatly increased. Furthermore, EINO_x simply reduced with decreasing combustor pressure at a rate of 0.035 g/10 kPa. The possibility of combustion control on the combusting mode and exhaust gas emission was demonstrated.     

Assessment of the performances of sub-grid scalar flux models for premixed flames with different global Lewis numbers: A Direct Numerical Simulation analysis

The statistical behaviours of sub-grid flux of reaction progress variable has been assessed for premixed turbulent flames with global Lewis number Le (=thermal diffusivity/mass diffusivity) ranging from 0.34 to 1.2 using a Direct Numerical Simulation (DNS) database of freely propagating statistically planar flames. It is known that the sub-grid scalar flux shows counter-gradient transport when the velocity jump across the flame due to heat release overcomes the effects of turbulent velocity fluctuation. The results show that the sub-grid scalar flux components exhibit counter-gradient transport for all cases considered here. The extent of counter-gradient transport increases with increasing filter width Δ and decreasing value of Le. This is due to the fact that flames with Le ≪ 1 (e.g. Le = 0.34) exhibit thermo-diffusive instabilities, which in turn increases the extent of counter-gradient transport. The effects of heat release and flame normal acceleration weaken with increasing Le. Several established algebraic models have been assessed in comparison to the sub-grid scalar flux components extracted from explicitly filtered DNS data in terms of their correlation coefficients at the vector level and their mean variation conditional on the Favre-filtered progress variable. The gradient transport closure does neither capture the quantitative nor the qualitative behaviour of the different sub-grid scalar flux components for all filter widths in all cases considered here. Models which account for local flame normal acceleration perform better, especially when the flame remains completely unresolved. In particular those models that account for the alignment of local resolved velocity and scalar gradients by using a tensor diffusivity, perform relatively better than the other alternative models irrespective of Le.