Structure of partially premixed n-heptane–air counterflow flames

To avoid the complexities associated with the droplet/vapor transport and nonuniform evaporation processes, a fundamental investigation of liquid fuel combustion in idealized configurations is very useful. An experimental–computational investigation of prevaporized n-heptane nonpremixed and partially premixed flames established in a counterflow burner is described. There is a general agreement between various facets of our nonpremixed flame measurements and the literature data. The partially premixed flames are characterized by a double flame structure. This becomes more distinct as the strain rate decreases and partial premixing increases, which also increases the separation distance between the two reaction zones. The peak partially premixed flame temperature increases with increasing premixing of the fuel stream. The peak CO₂ and H₂O concentrations are relatively insensitive to partial premixing. The CO and H₂ peak concentrations on the premixed side increase as the fuel-side equivalence ratio decreases. These species are transported to the nonpremixed reaction zone where they oxidize. The C₂ species have peaks in the premixed reaction zone. The concentrations of olefins are ten times larger than those of the corresponding paraffins. The oxidizer is present in partially premixed flames throughout the combustion system and there are no regions characterized by simultaneous high temperature and high fuel concentration. As a result, pyrolysis reactions leading to soot formation are greatly diminished.     

An apparatus-independent extinction strain rate in counterflow flames

Resistance to extinction by stretch is a key property of any flame, and recent work has shown that this property controls the overall structure of several important types of turbulent flames. Multiple definitions of the critical strain rate at extinction (ESR) have been presented in the literature. However, even if the same definition is used, different experiments report different extinction strain rates for flames burning the same fuel-air mixture at very similar temperatures using similarly constructed opposed-flow instruments. Here we show that at extinction, all these flames are essentially identical, so one would expect that each would be assigned the same value of a parameter representing its intrinsic resistance-to-stretch-induced-extinction, regardless of the specifics of the experimental apparatus. A similar situation arises in laminar flame speed measurements since different apparatuses could result in different strain rate distributions. In that instance, the community has agreed to report the unstretched laminar flame speed, and methods have been developed to translate the experimental (stretched) flame speed into a universal unstretched laminar flame speed. We propose an analogous method for translating experimental measurements for stretch-induced extinction into an unambiguous and apparatus-independent quantity (ESR_∞) by extrapolating to infinite opposing burner separation distance. The uniqueness of the flame at extinction is shown numerically and supported experimentally for twin premixed, single premixed, and diffusion flames at Lewis numbers greater than and less than one. A method for deriving ESR_∞ from finite-boundary experimental studies is proposed and demonstrated for methane and propane experimental diffusion and premixed single flame data. The two values agree within the range of ESR differences typically observed between experimental measurements and simulation results for the traditional ESR definition.     

Extinction and flame bifurcations of stretched dimethyl ether premixed flames

Extinction limits and flame bifurcation of lean premixed dimethyl ether–air flames are numerically investigated using the counterflow flame with a reduced chemistry. Emphasis is paid to the combined effect of radiation and flame stretch on the extinction and flammability limits. A method based on the reaction front is presented to predict the Markstein length. The predicted positive Markstein length agrees well with the experimental data. The results show that flow stretch significantly reduces the flame speed and narrows the flammability limit of the stretched dimethyl ether–air flame. It is found that the combined effect of radiation and flow stretch results in a new flame bifurcation and multiple flame regimes. At an equivalence ratio slightly higher than the flammability limit of the planar flame, the distant flame regime appears at low stretch rates. With an increase in the equivalence ratio, in addition to the distant flame, a weak flame isola emerges at moderate stretch rates. With a further increase in the equivalence ratio, the distant flame and the weak flame branches merge together, resulting in the splitting of the weak flame branch into two weak flame branches, one at low stretch and the other at high stretch. Flame stability analysis demonstrates that the high stretch weak flame is also stable. Furthermore, a K-shaped flammability limit diagram showing various flame regimes and their extinction limits is obtained.     

Effects of radiation heat loss on laminar premixed ammonia/air flames

Ammonia (NH₃) direct combustion is attracting attention for energy utilization without CO₂ emissions, but fundamental knowledge related to ammonia combustion is still insufficient. This study was designed to examine effects of radiation heat loss on laminar ammonia/air premixed flames because of their very low flame speeds. After numerical simulations for 1-D planar flames with and without radiation heat loss modeled by the optically thin model were conducted, effects of radiation heat loss on flame speeds, flame structure and emissions were investigated. Simulations were also conducted for methane/air mixtures as a reference. Effects of radiation heat loss on flame speeds were strong only near the flammability limits for methane, but were strong over widely diverse equivalence ratios for ammonia. The lower radiative flame temperature suppressed the thermal decomposition of unburned ammonia to hydrogen (H₂) at rich conditions. The equivalence ratio for a low emission window of ammonia and nitric oxide (NO) in the radiative condition shifted to a lower value than that in the adiabatic condition.     

Experimental and numerical determination of heat release in counterflow premixed laminar flames

This paper presents an experimental and numerical study of heat release in atmospheric laminar counterflow premixed flames. The measurements are based on simultaneous planar laser-induced fluorescence (PLIF) of OH and HCHO. These measurements are compared to numerical results obtained using detailed chemistry and multicomponent transport properties. A low Mach number formulation along the stagnation streamline is employed to describe the reactive flow. The conservation equations are completed with CHEMKIN and EGLIB packages. They are solved using finite differences, Newton iterations, and an adaptive gridding technique. The comparison is done along the burner axis for both, maximum heat release location and heat release profile width. It is shown that the product of OH and HCHO concentrations yields a result closely related to the heat release. These comparisons lead to the conclusion that the experimental method used seems to be a good tool for the determination of heat release in flames.     

A numerical investigation of stretch effects in counterflow,premixed laminar flames

We use direct numerical simulation of propane/air flames with full chemistry in the geometry of stagnation flow to investigate the effect of different definitions of local flame stretch in the presence of spatially varying velocity gradients. Specifically, we compare simulations with potential- and plug-flow inlet conditions, and show that the widely used definition of upstream stretch leads to unphysical results for flames having the ‘same’ stretch. We then show that a reasonable re-definition of local stretch allows us to produce the ‘same’ flame in the presence of the ‘same’ stretch.     

Effect of flame thickness variation on the mass burning rate of premixed counterflow flames with an oscillating strain rate

The mass-based stretch rate is used to study the response of premixed axisymmetric counterflow flames subject to an oscillating strain rate. Integral analysis is used to estimate the mass burning rate of the oscillating counterflow flames. From this study it can be concluded that the flame responds in a nonlinear manner. With an increase of the applied strain frequencies, it is found that unsteady stretch effects arising due to flame thickness variations become significant and the mass-based stretch rate is able to capture these nonlinear effects. The inclusion of these unsteady stretch effects in the mass-based stretch helps the integral analysis to predict the mass-burning rate of oscillating flames more accurately.     

Chemical structure of atmospheric pressure premixed laminar formic acid/hydrogen flames

The work presents an experimental and kinetic modeling study of laminar premixed formic acid [HC(O)OH]/H₂/O₂/Ar flames at different equivalence ratios (φ=0.85, 1.1 and 1.3) stabilized on a flat burner at atmospheric pressure, as well as laminar flame speed of HC(O)OH/O₂/Ar flames (φ=0.5–1.5) at 1 atm. Flame structure as well as laminar flame speed were simulated using three different detailed chemical kinetic mechanisms proposed for formic acid oxidation. The components in the fuel blends show different consumption profiles, namely, hydrogen is consumed slower than formic acid. According to kinetic analysis, the reason of the observed phenomenon is that the studied flames have hydrogen as a fuel but also as an intermediate product formed from HC(O)OH decomposition. Comparison of the measured and simulated flame structure shows that all the mechanisms satisfactorily predict the mole fraction profiles of the reactants, main products, and intermediates. It is noteworthy that the mechanisms proposed by Glarborg et al., Konnov et al. and the updated AramcoMech2.0 adequately predict the spatial variations in the mole fractions of free radicals, such as H, OH O and HO₂. However, some drawbacks of the mechanisms used were identified; in particular, they predict different concentrations of CH₂O. As for laminar flame speed simulations, the Konnov et al. mechanism predicts around two times higher values than in experiment, while the Glarborg et al. and updated AramcoMech2.0 show good agreement with the experimental data.     

The effect of hydrogen addition on flammability limit and NOx emission in ultra-lean counterflow CH4/air premixed flames

The effect of hydrogen addition to ultra lean counterflow CH₄/air premixed flames on the extinction limits and the characteristics of NO_x emission was investigated by numerical simulation. Detailed chemistry and complex thermal and transport properties were employed. The results show that the addition of hydrogen can significantly enlarge the flammable region and extend the flammability limit to lower equivalence ratios. If the equivalence ratio is kept constant, the addition of hydrogen increases the emission of NO in a flame due to the enhancement in the rate of the NNH or N₂O intermediate NO formation routes. The addition of hydrogen causes a monotonic decrease in the formation of NO₂ and N₂O, except flames near the extinction limits, where the emission of NO₂ and N₂O first increases, and then decreases with the increase in the fraction of hydrogen. Overall, hydrogen enrichment technology allows stable combustion under ultra lean conditions, resulting in significant CO₂ and NO emission reduction.     

Characteristic patterns of thermodiffusively unstable premixed lean hydrogen flames

Large-scale two-dimensional numerical simulations of thermodiffusively unstable, lean, premixed hydrogen flames have been performed using detailed finite rate chemistry to analyze flame intrinsic scales. The simulations feature a long integration time and large domain sizes to rule out effects of confinement on the dynamics of the flame front. For sufficiently large domain sizes, the total consumption speed of the flame is found to become independent of the domain size. An assessment of the characteristic scales of the flame front corrugation reveals the existence of a smallest and a largest flame intrinsic length scale. The smallest length manifests itself by local cusps, which lead to the formation of characteristic cells along the flame front. Their size is remarkably close to the most unstable wavelength predicted by a linear stability analysis of the flame front evolution in the linear regime. Independently of the domain size, a specific largest flame intrinsic structure, here referred to as flame finger, emerges from the interaction of multiple small-scale cusps. Thermodiffusively unstable flames are found to periodically form and destroy these flame fingers, but the formation of a global cusp that is known to emerge for purely hydrodynamically unstable flames is suppressed. The finite size of the largest scale fingers is explained by an instability in their movement. As they proceed towards the unburnt mixture, they tend to tilt and move laterally, thereby eventually being incorporated again by the rest of the flame. This behavior arises from the interaction of the flame fingers and the diverging velocity field ahead of them. Finally, the effect of equivalence ratio and unburnt gas temperature is investigated showing that flame fingers are found to develop only in case of a thermodiffusively unstable flame.     

Explosive dynamics of bluff-body-stabilized lean premixed hydrogen flames at blow-off

Two-dimensional direct numerical simulation (DNS) databases of bluff-body-stabilized lean hydrogen flames representative of complicated reactive–diffusive system are analysed using the combined approach of computational singular perturbation (CSP) and tangential stretching rate (TSR) to investigate chemical characteristics in blow-off dynamics. To assess the diagnostic approaches in flame and blow-off dynamics, Damköhler number and TSR variables are applied and compared. Four cases are considered in this study showing different flame dynamics such as the steadily stable mode, local extinction by asymmetric vortex shedding, convective blow-off and lean blow-out. DNS data points in positive explosive eigenvalue conditions were subdivided into four different combinations in TSR and extended TSR space and categorized in four distinct characteristic regions, such as kinetically explosive or dissipative and transport-enhanced or dissipative dynamics. The TSR analysis clearly captures the local extinction point in the complicated vortex shedding and allows an improved understanding of the distinct chemistry-transport interactions occurring in convective blow-off and lean blow-out events.     

Visible chemiluminescence of ammonia premixed flames and its application for flame diagnostics

We report a spatially resolved spectroscopic study of the visible chemiluminescence emission from different premixed ammonia-air-oxygen flames stabilized on a laminar flat flame burner, with equivalence ratio ranging from 0.7 to 1.35 and an O₂/N₂ ratio of 0.4. In the reaction zone of the observed flames, the visible emission was recognized to be the chemiluminescence of excited NH₂* radicals, while in the post-flame zone, two types of chemiluminescence were observed: NO₂* chemiluminescence dominated in the fuel-lean flames and NH₂* chemiluminescence dominated in the fuel-rich flames. The high-resolution spectra of the NO₂* and NH₂* chemiluminescence in the visible region (400-700 nm) were recorded. The intensity of both spectra increased gradually with wavelength. However, the NO₂*-chemiluminescence spectrum appeared to be continuous and unstructured, while the NH₂*-chemiluminescence spectrum consisted of groups of distinct emission lines. Based on the spectral feature, the ratios of the integrated chemiluminescence intensities over the 598-603 nm wavelength range to the intensities over the 586-592 nm range and 447-453 nm range were used to sense equivalence ratio. In addition, slightly different colors of the fuel-lean and fuel-rich flames were observed, due to the fact that NO₂* chemiluminescence had a relatively stronger signal in the blue region than NH₂* chemiluminescence. The difference was used to infer flame equivalence ratio using the flame images recorded by a RGB digital camera, where the ratio of the signal from the red channel to the signal from the blue channel was calculated.