Spatiotemporal measurements of flame stretch and propagation rates for lean and rich CH4/air premixed flames interacting with a turbulent-wake

A 1.5 m long turbulent-wake combustion vessel with a 0.15 m × 0.15 m cross-sectional area is proposed for spatiotemporal measurements of curvature, strain, dilatation and burning rates along a freely downward-propagating premixed flame interacting with a parallel row of staggered vortex pairs having both compression (negative) and extension (positive) strains simultaneously. The wanted wake is generated by rapidly withdrawing an electrically-controlled, horizontally-oriented sliding plate of 5 mm thickness for flame–wake interactions. Both rich and lean CH₄/air flames at the equivalence ratios  = 1.4 and  = 0.7 with nearly the same laminar burning velocity are studied, where flame–wake interactions and their time-dependent velocity fields are obtained by high-speed, high-resolution DPIV and laser-tomography. Correlations among curvature, strain, stretch, and dilatation rates along wrinkled flame fronts at different times are measured and thus their influences on front propagation rates can be analyzed. It is found that strain-related effects have significant influence on front propagation rates of rich CH₄/air (diffusionally stable) flames even when the curvature weights more in the total stretch than the strain rate does. The local propagation rates along the wrinkled flame front are more intense at negative strain rates corresponding to positive peak dilatation rates but the global propagation rate averaged along the rich flame front remains constant during all period of flame–wake interaction. For lean CH₄/air (diffusionally unstable) flames, the curvature becomes a dominant parameter influencing the structure and propagation of the wrinkled flame front, where both local and global propagation rates increase significantly with time, showing unsteady flame propagation. These experimental results suggest that the theory of laminar flame stretch can be applicable to a more complex flame–wake interaction involving unsteadiness and multitudinous interactions between vortices.     

Mechanism of combustion dynamics in a backward-facing step stabilized premixed flame

Combustion dynamics leading to thermoacoustic instability in a rearward-facing step stabilized premixed flame is experimentally examined with the objective of investigating the fluid dynamic mechanism that drives heat release rate fluctuations, and how it couples with the acoustic field. The field is probed visually, using linear photodiode arrays that capture the spatiotemporal distribution of CH* and OH*; an equivalence ratio monitor; and a number of pressure sensors. Results show resonance between the acoustic quarter wave mode of the combustion tunnel and a fluid dynamic mode of the wake. Under unstable conditions, the flame is convoluted around a large vortex that extends several step heights downstream. During a typical cycle, while the velocity is decreasing, the vortex grows, and the flame extends downstream around its outer edge. As the velocity reaches its minimum, becoming mostly negative, the vortex reaches its maximum size, and the flame collides with the upper wall; its leading edge folds, trapping reactants pockets, and its trailing edge propagates far upstream of the step. In the next phase, while the velocity is increasing, the heat release grows rapidly as trapped reactant’ pockets are consumed by flames converging towards their centers, and the upstream flame is dislodged back downstream. The heat release rate reaches its maximum halfway into the velocity rise period, leading the maximum velocity by about 90°. In this quarter-wave mode, the pressure leads the velocity by 90° as well, that is, it is in phase with the heat release rate. Numerical modeling results support this mechanism. Equivalence ratio contribution to the instability mechanism is shown to be minor, i.e., heat release dynamics are governed by the cyclical formation of the wake vortex and its interaction with the flame.     

Subgrid combustion modeling of 3-D premixed flames in the thin-reaction-zone regime

Large eddy simulation (LES) is used to investigate three-dimensional (3D) lean premixed turbulent methane–air flames in the thin-reaction-zone regime. In this regime, the Kolmogorov scale is smaller than the preheat zone thickness, but larger than the reaction zone thickness. Past numerical studies of similar flames were primarily direct numerical simulation either in two-dimensions or using the artificially thickened flame approach in 3D. For an LES the effect of small (unresolved) scales on the scalar field must be, modeled accurately to capture the correct flame structure. A subgrid combustion model based on the linear-eddy-mixing (LEM) model is used within an LES framework (called LEM–LES hereafter) to capture the 3D flame-structure of the highly stretched premixed flames. A finite-rate, one-step methane–air chemistry with a non-unity Lewis number formulation is used in this study. The simulated flame structure resembles flames experimentally studied in the thin-reaction-zone regime. Even though the preheat zone is broadened by the penetration of small eddies, the chemical reaction zone remains thin and localized. This feature is captured properly in the current LEM–LES approach. The flame structure and other statistics such as the flame area evolution, curvature, and strain-rate statistics computed using the LEM–LES are also in good agreement with the past DNS studies.     

Temporal evolution of flame stretch due to turbulence and the hydrodynamic instability

The temporal evolution of the strain rate on a turbulent premixed flame was measured experimentally using cinema-stereoscopic particle image velocimetry. Turbulence strains a flame due to velocity gradients associated both directly with the turbulence and those caused by the hydrodynamic instability, which are initiated by the turbulence. The development of flame wrinkles caused by both of these mechanisms was observed. Wrinkles generated by the turbulence formed around vortical structures, which passed through the flame and were attenuated. After the turbulent structures had passed, the hydrodynamic instability flow pattern developed and caused additional strain. The hydrodynamic instability also caused the growth of small flame front perturbations into large wrinkles. In the moderately turbulent flame investigated, it was found that the evolution of the strain rate caused by turbulence–flame interactions followed a common pattern involving three temporal regimes. In the first, the turbulence exerted extensive (positive) strain on the flame, creating a wrinkle that had negative curvature (concave towards the reactants). This was followed by a transition period, leading into the third regime in which the flow pattern and strain rate were dominated by the hydrodynamic instability mechanism. It was also found that the magnitudes of the strain rate in the first and third regimes were similar. Hence, the hydrodynamic instability mechanism caused significant strain on a flame and should be included in turbulent combustion models.     

An experimental study of partially premixed flame sound

The occurrence of oscillating combustion and combustion instability has led to resurgence of interest in the causes, mechanisms, suppression, and control of combustion noise. Noise generated by enclosed flames is of greater practical interest but is more complicated than that by open flames, which itself is not clearly understood. Studies have shown that different modes of combustion, premixed and non-premixed, differ in their sound generation characteristics. However, there is lack of understanding of the region bridging these two combustion modes. This study investigates sound generation by partially premixed flames. Starting from a non-premixed flame, air was gradually added to achieve partial premixing while maintaining the fuel flow rate constant. Methane, ethylene, and ethane partially premixed flames were studied with hydrogen added for flame stabilization. The sound pressure generated by methane partially premixed flames scales with M⁵ compared to M³ for turbulent non-premixed methane flames. Also, the sound pressure generated by partially premixed flames of ethane and ethylene scales as M^4.5. With progressive partial premixing, spectra level increases at all frequencies with a greater increase in the high-frequency region compared to the low-frequency region; flames develop a peak and later a constant level plateau in the low frequency region. The partially premixed flames of methane, ethylene, and ethane generate a similar SPL as a function of equivalence ratio when the fuel volume flow rate is matched. However, when fuel mass flow rate is matched, the ethane and ethylene flames produce a similar SPL, which is lower than that produced by the methane flame.     

Experimental and LES investigation of premixed methane/air flame propagating in a tube with a thin obstacle

In this paper, an experimental and numerical investigation of premixed methane/air flame dynamics in a closed combustion vessel with a thin obstacle is described. In the experiment, high-speed video photography and a pressure transducer are used to study the flame shape changes and pressure dynamics. In the numerical simulation, four sub-grid scale viscosity models and three sub-grid scale combustion models are evaluated for their individual prediction compared with the experimental data. High-speed photographs show that the flame propagation process can be divided into five stages: spherical flame, finger-shaped flame, jet flame, mushroom-shaped flame and bidirectional propagation flame. Compared with the other sub-grid scale viscosity models and sub-grid scale combustion models, the dynamic Smagorinsky–Lilly model and the power-law flame wrinkling model are better able to predict the flame behaviour, respectively. Thus, coupling the dynamic Smagorinsky–Lilly model and the power-law flame wrinkling model, the numerical results demonstrate that flame shape change is a purely hydrodynamic phenomenon, and the mushroom-shaped flame and bidirectional propagation flame are the result of flame–vortex interaction. In addition, the transition from “corrugated flamelets” to “thin reaction zones” is observed in the simulation.     

Measurements of flame orientation and scalar dissipation in turbulent partially premixed methane flames

Three turbulent flames were studied using a new experimental facility developed at Sandia National Laboratories. Line imaging of Raman and Rayleigh scattering and CO laser-induced fluorescence (LIF) yielded information on all major species, temperature, mixture fraction, and a 1D surrogate measure of scalar dissipation. Simultaneously, crossed planar OH LIF imaging provided information on the instantaneous flame orientation, allowing estimation of the full 3D (flame-normal) scalar dissipation rate. The three flames studied were methane–air piloted jet flames (Sandia flames C, D, and E), which cover a range in Reynolds number from 13,400 to 33,600. The statistics of the instantaneous flame orientation are examined in the different flames, with the purpose of studying the prevailing kinematics of isoscalar contours. The 1D and 3D results for scalar dissipation rate are examined in detail, both in the form of conditional averages and in the form of probability density functions. The effect of overall strain and Reynolds number on flame suppression and eventual extinction is also investigated, by examining the doubly conditional statistics of temperature in the form of S-shaped curves. This latter analysis reveals that double conditioning of temperature on both mixture fraction and scalar dissipation does not collapse the data from these flames onto the same curve at low scalar dissipation rates, as might be expected from simple flamelet concepts.     

Experimental and numerical study of a conical turbulent partially premixed flame

The structure and dynamics of a turbulent partially premixed methane/air flame in a conical burner were investigated using laser diagnostics and large-eddy simulations (LES). The flame structure inside the cone was characterized in detail using LES based on a two-scalar flamelet model, with the mixture fraction for the mixing field and level-set G-function for the partially premixed flame front propagation. In addition, planar laser induced florescence (PLIF) of CH and chemiluminescence imaging with high speed video were performed through a glass cone. CH and CH₂O PLIF were also used to examine the flame structures above the cone. It is shown that in the entire flame the CH layer remains very thin, whereas the CH₂O layer is rather thick. The flame is stabilized inside the cone a short distance above the nozzle. The stabilization of the flame can be simulated by the triple-flame model but not the flamelet-quenching model. The results show that flame stabilization in the cone is a result of premixed flame front propagation and flow reversal near the wall of the cone which is deemed to be dependent on the cone angle. Flamelet based LES is shown to capture the measured CH structures whereas the predicted CH₂O structure is somewhat thinner than the experiments.     

Experimental and numerical study of premixed, lean ethylene flames

Ethylene is a key intermediate in the combustion mechanisms of most practical fuels. It plays also an important role in the formation of aromatic hydrocarbons and soot particules. The latter has motivated many experimental and numerical studies carried out on rich ethylene-air mixtures. Less studies have been devoted to lean mixtures, and the development of strategies based on lean, premixed flames to reduce soot and NO_x production requires additional experimental data in lean conditions. In this work, the chemical structure of lean premixed ethylene-oxygen-nitrogen flames stabilized on a flat-flame burner at atmospheric pressure was determined experimentally. The species mole fraction profiles were also computed by the Premix code (Chemkin II version) and four detailed reaction mechanisms. A very good agreement was observed for the main flame properties: reactants consumption, final products (CO₂, H₂O) and the main intermediates: CO and H₂. Marked differences occurred in the prediction of active intermediate species present in small concentrations. Pathways analyses were performed to identify the origins of these discrepancies. It was shown that the same reactions were involved in the four mechanisms to describe the consumption of ethylene, but with marked differences in their relative importance. C₂H₃ and CH₂HCO are the main radicals formed in this first step and their consumption increases the differences between the mechanisms either by the use of different kinetic data for common reactions or by differences in the nature of the consumption reactions.     

An experimental study of the rich premixed ethylbenzene flame at low pressure

A slightly sooting premixed ethylbenzene flame with an equivalence ratio of 1.90 was investigated at low pressure (4.0 kPa) using molecular-beam mass spectrometry (MBMS) and tunable synchrotron vacuum ultraviolet (VUV) photoionization. Basing on the ionization threshold measurements of photoionization efficiency (PIE) spectra, combustion intermediates up to C₁₉H₁₂ were identified, including a number of radicals and isomeric species. Mole fraction profiles of observed flame species were evaluated from the measurements of burner scan at the photon energies near ionization thresholds. Besides, the flame temperature profile was measured by a Pt/Pt-13%Rh thermocouple. From the intermediate identification and mole fraction measurements, the degradation of ethylbenzene, as well as the formation of some interested polycyclic aromatic hydrocarbons (PAHs), was discussed in detail. It is suggested that the formation of most typical PAHs observed in this work can be related to the H-abstraction/C₂H₂-addition (HACA) mechanism. Furthermore, the high concentration levels of intermediates in this flame is ascribed to the weak C-C bonds in the sidechain of ethylbenzene, which provides a potential explanation of the high sooting tendencies of ethylbenzene and other monocyclic aromatic fuels with complex sidechain structure. This study is anticipated to be constructive for combustion investigations of aromatic fuels, and the discussion is hoped to be helpful for further modeling studies concerning PAHs formation in combustion process.     

An experimental and kinetic modeling study of a premixed nitromethane flame at low pressure

A premixed nitromethane/oxygen/argon flame at low pressure (4.67 kPa) has been investigated using tunable vacuum ultraviolet (VUV) photoionization and molecular-beam mass spectrometry. About 30 flame species including hydrocarbons, oxygenated and nitrogenous intermediates have been identified by measurements of photoionization efficiency spectra. Mole fraction profiles of the flame species have been determined by scanning burner position at some selected photon energies. The results indicate that N₂ and NO are the major nitrogenous products in the nitromethane flame. Compared with previous studies on nitromethane combustion, a number of unreported intermediates, including C₃H₄, C₄H₆, C₄H₈, C₂H₂O, C₂H₄O, CH₃CN, H₂CNHO, C₃H₃N and C₃H₇N, are observed in this work. Based on our experimental results and previous modeling studies, a detailed oxidation mechanism including 69 species and 314 reactions has been developed to simulate the flame structure. Despite some small discrepancies, the predictions by the modeling study are in reasonable agreement with the experimental results.     

An experimental and kinetic modeling study of premixed nitroethane flames at low pressure

An experimental and kinetic modeling study is reported on three premixed nitroethane/oxygen/argon flames at low pressure (4.655 kPa) with the equivalence ratios (Φ) of 1.0, 1.5 and 2.0. Over 30 flame species were identified with tunable synchrotron vacuum ultraviolet photoionization mass spectrometry, with their mole fractions quantified as the function of the height above burner. The flame temperature profiles were measured with a Pt–6%Rh/Pt–30%Rh thermocouple. A detailed kinetic mechanism with 115 species and 730 reactions was proposed and validated against experimental results. The computed predictions have shown satisfactory agreement with the experimental results. Basing on the rate-of-production analysis, the reaction pathways that feature the combustion of nitroethane were revealed, including the primary decomposition of C–N bond fission, the oxidation of C2 and C1 hydrocarbons and the formation of nitrogenous species. The presence of NO₂ and NO has been proved to be important for these processes.     

An experimental and kinetic modeling study on nitric oxide formation in premixed C3 alcohols flames

This study provides new quantitative NO concentrations measurements in n-propanol + air and i-propanol + air flames together with a new combustion kinetic model. The heat flux method was employed to stabilize propyl alcohols flames and the initial gas conditions were set to 323 K, 1 atm, and Φ=0.7–1.4. Saturated laser-induced fluorescence was employed to measure NO concentration in the post-combustion region. The presented and literature models, namely the POLIMI and Bohon et al. (2018) kinetic mechanisms, were assessed against new experimental data. Experimental results showed a higher NO formation in the thermal zone for n-propanol flames, whereas i-propanol flames indicate a higher amount of NO formed at fuel-rich conditions. Overall among the tested models, the present mechanism exhibited the best agreement in emulating NO experimental profiles; conversely, numerical simulations from the POLIMI model showed significant inconsistencies at fuel-rich conditions and the Bohon et al. (2018) model was unable to reproduce the measured data, notably underpredicting experimental values at all investigated conditions. However, the present model manifested some uncertainties in reproducing NO formation in the prompt region; therefore, in connection with this important aspect, the new experimental data obtained in this work will provide a valid support to further develop more reliable kinetic models.     

Structural response of different Lewis number premixed flames interacting with a toroidal vortex

Simultaneous measurements of temperature, CH* and OH* chemiluminescent species are carried out to explore the impact of stretch rate and curvature on the structure of premixed flames. The configuration of an initially flat premixed flame interacting with a toroidal vortex is selected for the present study and reasons for this choice are discussed. Lewis number effects are assessed by comparing methane and propane flames. It is emphasized that the flame structure experiences very strong variations. In particular, the flame is shrunk (broadened) in the initial (final) period of the interaction with the vortex where strain rate (curvature) contribution of the stretch rate is predominant. By further analysing independently the thickness of the preheat and reaction zones, it is shown that for propane flames, not only the former but also the latter is significantly altered in zones where the flame curvature is negative. Changes in the reaction zone properties are further emphasized using CH* and OH* radicals. It is demonstrated that higher thermal diffusivity plays a significant role around curved regions, in which the enhanced diffusion of heat leads to a strong increase of CH* compared to OH* intensity. As an overall conclusion, this study suggests that it would be interesting to reassess the internal flame structure at lower and moderate Karlovitz numbers since changes might appear for a moderate vortex intensity with typical size much larger than the flame thickness.     

An investigation of flame zones and burning velocities of laminar unconfined methane-oxygen premixed flames

Numerical and experimental investigations of unconfined methane-oxygen laminar premixed flames are presented. In a lab-scale burner, premixed flame experiments have been conducted using pure methane and pure oxygen mixtures having different equivalence ratios. Digital photographs of the flames have been captured and the radial temperature profiles at different axial locations have been measured using a thermocouple. Numerical simulations have been carried out with a C2 chemical mechanism having 25 species and 121 reactions and with an optically thin radiation sub-model. The numerical results are validated against the experimental and numerical results for methane-air premixed flames reported in literature. Further, the numerical results are validated against the results from the present methane-oxygen flame experiments. Visible regions in digital flame photographs have been compared with OH isopleths predicted by the numerical model. Parametric studies have been carried out for a range of equivalence ratios, varying from 0.24 to 1.55. The contours of OH, temperature and mass fractions of product species such as CO, CO₂ and H₂O, are presented and discussed for various cases. By using the net methane consumption rate, an estimate of the laminar flame speed has been obtained as a function of equivalence ratio.     

A numerical study on the effect of CO addition on extinction limits and NO x formation in lean counterflow CH4/air premixed flames

The effect of CO addition on extinction and NO _x formation in lean premixed counterflow CH₄/air flames was investigated by numerical simulation. Detailed chemistry and complex thermal and transport properties were employed. A method that gradually switched off the initial reactions of NO formation from different routes was used to analyse the variation of NO formation mechanism. The results indicate that the addition of certain amount of CO increases the strain extinction limits and reduces the radiation extinction limits. As a result, the lean flammability limit of CH₄/air premixed flame is extended to leaner side by the addition of CO. The formation of NO in a flame is increased with the addition of CO at a constant equivalence ratio. For an ultra-lean flame, the increase in the formation of NO is mainly because of the increase in the contribution from the NNH intermediate route, while for a near stoichiometric flame, this increase is mainly attributable to the rise in the contribution from the thermal route. With the fraction of added CO being gradually increased, the formation of NO₂ in a flame first decreases and then increases at a given equivalence ratio. The addition of CO reduces the formation N₂O in an ultra-lean flame, while affects little on the formation of N₂O in a near stoichiometric flame.     

Strain rates,flow patterns and flame surface densities in hydrodynamically unstable,weakly turbulent premixed flames

Recent numerical and experimental studies have unveiled a potentially marked difference between the laminar as well as turbulent propagation of premixed flames exhibiting Darrieus–Landau (DL) (or hydrodynamic) instabilities from flames for which instabilities are inhibited. In this study we utilize two-dimensional numerical simulations of slot burner flames as well as experimental Propane–Air Bunsen flames to analyse differences in turbulent propagation, strain rate and induced flow patterns of hydrodynamically stable and unstable flames. We also investigate the effects of hydrodynamic instability on quantities which are directly related to reaction rate closure models, such as flame surface density and stretch factor. A clear enhancement of turbulent flame speed can be observed for unstable flames, generally mitigated at higher turbulence intensity, which is attributed to a flame area increase induced by the characteristic cusp-like DL-induced corrugation, absent in stable flames, which occurs concurrently and in synergy with turbulent wrinkling. Unstable flames also exhibit, both numerically and experimentally, a different correlation between strain rate and flame curvature and are observed to give rise to a channeling of the induced flow in the fresh mixture. Conditionally averaged flame surface density is also observed to attain smaller values in unstable flames, as a result of the thicker turbulent flame brush, indicating that closure models should incorporate instability-related parameters in addition to turbulence-related parameters.     

Particle formation in opposed-flow diffusion flames of ethylene: An experimental and numerical study

An experimental and numerical study on particles inception and growth is performed in opposed-flow diffusion flames of ethylene and air characterized by different sooting tendencies. Spectrally resolved UV-visible laser induced fluorescence, laser induced incandescence and laser light scattering measurements are used to characterize different classes of combustion-generated compounds. A detailed kinetic model accounting for both gas-phase and particle formation is used. Comparison between experimental results and numerical predictions gives a qualitative view of the mechanism of particle formation in opposed-flow flames.Particle inception is the result of both chemical growth and coagulation of aromatic compounds. In the region close to the flame front where the temperature is relatively high and radicals are abundant, the particle inception is due to a chemical growth mechanism by which aromatic molecules add aromatic radicals leading to the formation of biphenyl-like structures. The growth process continues as high-molecular mass aromatics are moved away from the flame zone towards the stagnation plane by the addition of acetylene and other aromatics forming particles of increasing sizes. Graphitization of these particles and thermal annealing lead to the formation of soot particles. At relatively lower temperatures, found across the stagnation plane, particles growth still occurs and it is mainly due to a process of physical coagulation of PAHs.The experimental and numerical results obtained in this work demonstrate and explain the sensibility of inception and growth of particles to radical concentration and temperature in opposed-flow flame configurations.