Measurements of burning velocities of dimethyl ether and air premixed flames at elevated pressures

Laminar burning velocities of dimethyl ether (DME) and air premixed flames at elevated pressures up to 10 atm were measured by using a newly developed pressure-release type spherical bomb. The measurement system was validated using laminar burning velocities of methane–air flames. A comparison with the previous experimental data shows an excellent agreement and demonstrates the accuracy and reliability of the present experimental system. The measured flame speeds of DME–air flames were compared with the previous experimental data and the predictions using the full and reduced mechanisms. At atmospheric pressure, the measured laminar burning velocities of DME–air flames are in reasonable agreement with the previous data from spherical bomb method, but are much lower than both predictions and the experimental data of the PIV based counterflow flame measurements. The laminar burning velocities of DME–air flames at 2, 6, and 10 atm were also measured. It was found that flame speed decreases considerably with the increase of pressure. Moreover, the measured flame speeds are also lower than the predictions at high pressures. In addition, experiments showed that at high pressures the rich DME–air flames are strongly affected by the hydrodynamic and thermal-diffusive instabilities. Markstein lengths and the overall reaction order at different equivalence ratios were extracted from the flame speed data at elevated pressures. Sensitivity analysis showed that reactions involving methyl and formyl radicals play an important role in DME–air flame propagation and suggested that systematic modification of the reactions rates associated with methyl and formyl formations are necessary to reduce the discrepancies between predictions and measurements.     

Influences of heat flux on extinction characteristics of steady/unsteady premixed stagnation flames

This paper investigates the extinction characteristics of premixed stagnation flames (PSFs) with controlled heat losses and flow disturbances. The low-frequency air flow pulsations that imitate the operational transients in practical combustors were specially introduced. The tunable diode-laser absorption spectroscopy (TDLAS) measurement was applied to obtain the temperature profile and wall heat flux. It is found that, for steady flame with a fixed equivalence ratio, the extinction stretch rate dramatically increases as the wall heat flux decreases. The extinction criterion is summarized as a global Karlovitz number of 0.57 by establishing a relationship between the global and local stretch rates. Numerical simulations reveal that the local extinction Karlovitz number of steady PSFs is approximately 1.0 regardless of the conditions such as heat flux and equivalence ratio. Further experiments present that the air pulsations with a repetition of ～5 Hz significantly deteriorate the flame stability. Particularly, for unsteady perturbed flames, the extinction stretch rate exhibits a nonlinear trend, yielding two regimes with discrepant sensitivities to wall heat flux. The unsteady simulation then highlights a local stretch rate overshoot in the presence of pulsation. It is caused by the time delay between the inlet velocity and flame front movement that eventually leads to poor flame stability. Moreover, in the high heat-flux regime, a smaller local stretch rate overshoot results in the weak dependence of extinction limits on heat fluxes.     

Flame speed and self-similar propagation of expanding turbulent premixed flames

In this Letter we present turbulent flame speeds and their scaling from experimental measurements on constant-pressure, unity Lewis number expanding turbulent flames, propagating in nearly homogeneous isotropic turbulence in a dual-chamber, fan-stirred vessel. It is found that the normalized turbulent flame speed as a function of the average radius scales as a turbulent Reynolds number to the one-half power, where the average radius is the length scale and the thermal diffusivity is the transport property, thus showing self-similar propagation. Utilizing this dependence it is found that the turbulent flame speeds from the present expanding flames and those from the Bunsen geometry in the literature can be unified by a turbulent Reynolds number based on flame length scales using recent theoretical results obtained by spectral closure of the transformed G equation.     

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.     

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.     

Minimum ignition energy and propagation dynamics of laminar premixed cool flames

Laminar premixed cool flames, induced by the coupling of low-temperature chemistry and convective-diffusive transport process, have recently attracted extensive interest in combustion and engine research. In this work, numerical simulations have been conducted using a recently developed open-source reacting flow platform reactingFOAM-SCT, to investigate the minimum ignition energy (MIE) and propagation dynamics of premixed cool flames in a 1D spherical coordinate. Results have shown that when ignition energy is below the MIE of regular hot flames, a class of cool flames could be initiated, which allow much wider flammability limits, both lean and rich, compared to hot flames. Furthermore, the overall cool flame propagation dynamics exhibit intrinsic similarity to those of hot flames, in that, they begin with an ignition kernel propagation regime, followed by two transition regimes, and eventually reach a normal flame propagation regime. However, a spherical expanding cool flame responds completely differently to stretch. Specifically, a regular outwardly propagating hot spherical flame accelerates with increasing stretch rate when the mixture Le < 1 and decelerates when Le > 1. However, it is found that a cool flame always tends to decelerate with increasing stretch rate regardless of mixture composition, exhibiting unique flame aerodynamic characteristic. This research discovers novel features of premixed cool flame initiation and propagation dynamics and sheds light on flame transition, spark-ignition system design, and advanced engine combustion control.     

Effect of H-atom concentration on soot formation in premixed ethylene/air flames

Soot formation is a major challenge in the development of clean and efficient combustion systems based on hydrocarbon fuels. Fundamental understanding of the reaction mechanism leading to soot formation can be obtained by investigating the role of key reactive species such as atomic hydrogen taking part in soot formation pathways. In this study, two-dimensional laser induced incandescence (LII) measurements using λ?=?1064?nm laser have been used to measure soot volume fraction (f_V) in a series of rich ethylene (C₂H₄)/air flames, stabilized over a McKenna burner fitted with a flame stabilizing metal disc. Moreover, a comparison of UV (λ?=?283?nm), visible (λ?=?532?nm) and IR (λ?=?1064?nm) laser excited LII measurements of soot is discussed. Recently developed, femtosecond two-photon laser-induced fluorescence (fs-TPLIF) technique has been applied for obtaining spatially resolved H-atom concentration ([H]) profiles under the same flame conditions. The structure of the flames has also been determined using hydroxyl radical (OH) planar laser induced fluorescence (PLIF) imaging. The results indicate an inverse dependence of f_V on [H] for a range of C₂H₄/air rich flames up to an equivalence ratio, Φ?=?3.0. Although an absolute relationship between [H] and f_V cannot be easily derived owing to the multiple steps involving H and other intermediate species in soot formation pathways, the present study demonstrates the feasibility to couple [H] and f_V obtained using advanced optical techniques for soot formation studies.     

Cavity ring-down measurements of seeded NO in premixed atmospheric-pressure H2/air and CH4/air flames

Quantitative aspects of using cavity ring-down absorption spectroscopy near 226 nm for measurements of NO mole fractions in premixed atmospheric-pressure flames are discussed. Measurements in methane–air flames showed strong broadband absorption near 226 nm by hot CO₂ molecules, precluding using the cavity ring-down method in these flames at atmospheric pressure. In hydrogen–air flames, the broadband absorption at this wavelength was substantially lower. Absorption cross sections derived from non-seeded cavity ring-down spectra suggest that absorption by water is the major contribution to the background in these flames. The detectability limit for NO by cavity ring-down measurements in hydrogen–air flames using the current setup is estimated to be 10 ppm. Effects of the cold boundary layer on the measured NO mole fraction were accounted for by measuring the radial distributions of temperature and NO mole fraction using coherent anti-Stokes Raman scattering and laser-induced fluorescence (LIF), respectively. Measurements performed in seeded stoichiometric and lean hydrogen–air flames showed no reburning at temperatures above 1750 K, demonstrating the adequacy of using these flames for calibration of LIF measurements. At lower temperatures, the mole fraction of NO in the hot gases was up to 30% lower than that expected from the degree of seeding in the cold gases. PACS 42.62.Fi; 42.68.Ca; 82.33.Vx     

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.     

Unsteady propagation of premixed methane/propane flames in a mesoscale disk burner of variable-gaps

Unsteady flame propagation of air-premixed methane and propane flames was investigated in a new mesoscale disk burner, of which disk-gap could be precisely varied. To begin with, the quenching disk-gaps on the flammability limits were measured. In most cases, with the slight increase of the disk-gap, cellular flame structures could be generated. The initiation of such cellular structures could be explained by the thermally induced hydrodynamic instability, and it could be enhanced if the Lewis number was sufficiently small. When the disk-gap was sufficiently larger than a critical value (approximately 1.5 times the quenching distance), the cellular flame structure changed into a smooth one in the azimuthal direction. With a further increase of the disk-gap, the flame propagation velocities approached to constant values. These values were comparable to the laminar burning velocities except for the propane-rich conditions, in which much larger propagation velocities were observed. The flame stretch effects (coupled with Le-effects) within a narrow space were suspected as the reason. The structural transition of the premixed flame could be investigated successfully through various disk-gaps, from the smallest quenching-scale to the ordinary large scales via whole mesoscales including Hele–Shaw scales.     

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.     

Direct numerical simulations of rich premixed turbulent n-dodecane/air flames at diesel engine conditions

Rich premixed turbulent n-dodecane/air flames at diesel engine conditions are analyzed using direct numerical simulations. The conditions correspond to a parametric variation of the Engine Combustion Network Spray A (pressure 60 atm; oxidizer oxygen level and temperature 21% and 900 K, respectively; fuel temperature 363 K). Three simulations with equivalence ratios of 3, 5, and 7 are performed with a Karlovitz number (Ka, based on flame time) of order 100 to match the estimated Ka of the rich premixed combustion region in Spray A. At these conditions, the reference laminar flames exhibit a complex structure which involves both low-temperature chemistry (LTC) and high-temperature chemistry over a wide range of length scales. In the presence of turbulence, the flame structure is strongly affected in physical space and the reaction zone exhibits a very complex structure in which broken, distributed, and thin regions co-exist, especially for the leanest case. However, the contribution of the LTC pathway is only weakly affected by turbulence. In progress variable space, the mean flame structure, including the chemical source terms, is found to match remarkably well that of the corresponding unity Lewis number laminar flame, particularly for the $?=$ 3 and 5 cases. This behavior is attributed to the strong turbulent mixing occurring throughout the flames/reaction zones, which suppresses differential diffusion effects. Nevertheless, large conditional fluctuations around the mean chemical source terms are identified. These are found to correlate very well with radical species mass fractions such as OH. In addition, a similar functional dependence is obtained from counterflow laminar flames. As such, it appears from these results that laminar flame models have a potential to be used to represent the thermochemical state of rich premixed turbulent flames under diesel engine conditions.     

Burning velocity correlation of methane/air turbulent premixed flames at high pressure and high temperature

Turbulent burning velocities for methane/air mixtures at pressures ranging from atmospheric pressure up to 1.0 MPa and mixture temperatures of 300 and 573 K were measured, which covers the typical operating conditions of premixed-type gas-turbine combustors. A bunsen-type flame stabilized in a high-pressure chamber was used, and OH-PLIF visualization was performed with the pressure and mixture temperature being kept constant. In addition to a burner with an outlet diameter of 20 mm for the high-pressure experiments, a large-scale burner with an outlet diameter of 60 mm was used at atmospheric pressure to extend the turbulence Reynolds number based on the Taylor microscale, R_λ, as a common parameter to compare the pressure and temperature effects. It was confirmed that R_λ over 100 could be attained and that u′/S_L could be extended even at atmospheric pressure. Based on the contours of the mean progress variable c = 0.1 determined using OH-PLIF images, turbulent burning velocity was measured. S_T/S_L was also found to be greatly affected by pressure for preheated mixtures at 573 K. The bending tendency of the S_T/S_L curves with u′/S_L was seen regardless of pressure and mixture temperature and the R_λ region where the bending occurs corresponded well to the region where the smallest scale of flame wrinkling measured as a fractal inner-cutoff approaches the characteristic flame instability scale and becomes almost constant. A power law of S_T/S_L with (P/P₀)(u′/S_L) was clearly seen when S_T was determined using c = 0.1 contours, and the exponent was close to 0.4, indicating agreement with the previous results using the mean flame cone method and the significant pressure effects on turbulent burning velocity.     

Improvement of two-photon laser induced fluorescence measurements of H- and O-atoms in premixed methane/air flames

3 photodissociation is shown to be an important photolytic source of H atoms in the reaction zone of the methane flames. At 226 nm, an efficient energy transfer between O(³P) and N₂ is established from the observation of O-resonant emissions from the second positive system of N₂. The subsequent rate of O(³P) depletion appears to be essentially “controlled” by the O(³P) concentration and is probably only minimally temperature dependent. Finally quenching rate coefficients for O(³P) by water and nitrogen are deduced from quenching rates measurements performed in CH₄/O₂/N₂ and CH₄/O₂/Ar flames. Received: 7 November 1996/Revised version: 28 January 1997     

Measurements of three-dimensional mean flame surface area ratio in turbulent premixed Bunsen flames

A data processing scheme with particular emphasis on proper flame contour smoothing is developed and applied to measure the three-dimensional mean flame surface area ratio in turbulent premixed flames. The scheme is based on the two-sheet imaging technique such that the mean flame surface area ratio is an average within a window covering a finite section of the turbulent flame brush. This is in contrast to the crossed-plane tomograph technique which applies only to a line. Two sets of Bunsen flames have been investigated in this work with the turbulent Reynolds number up to 4000 and the Damköhler number ranging from less than unity to close to 10. The results show that three-dimensional effects are substantial. The measured three-dimensional mean flame surface area ratio correlates well with a formula similar to the Zimont model for turbulent burning velocity but with different model constants. Also, the mean flame surface area ratio displays a weak dependency on turbulence intensity but a strong positive dependency on the turbulence integral length scale.     

Structure and propagation of two-dimensional,partially premixed,laminar flames in diesel engine conditions

We investigate the influence of inflow velocity (V_in) and scalar dissipation rate (χ) on the flame structure and stabilisation mechanism of steady, laminar partially premixed n-dodecane edge flames stabilised on a convective mixing layer. Numerical simulations were performed for three different χ profiles and several V_in (V_in = 0.2 to 2.5m/s). The ambient thermochemical conditions were the same as the Engine Combustion Network’s (ECN) Spray A flame, which in turn represents conditions in a typical heavy duty diesel engine. The results of a combustion mode analysis of the simulations indicate that the flame structure and stabilisation mechanism depend on V_in and χ. For low V_in the flame is attached. Increasing V_in causes the high-temperature chemistry (HTC) flame to lift-off, while the low-temperature chemistry (LTC) flame is still attached. A unique speed S_R associated with this transition is defined as the velocity at which the lifted height has the maximum sensitivity to changes in V_in. This transition velocity is negatively correlated with χ. Near $V_{in}=S_R$ a tetrabrachial flame structure is observed consisting of a triple flame, stabilised by flame propagation into the products of an upstream LTC branch. Further increasing the inlet velocity changes the flame structure to a pentabrachial one, where an additional HTC ignition branch is observed upstream of the triple flame and ignition begins to contribute to the flame stabilisation. At large V_in, the LTC is eventually lifted, and the speed at which this transition occurs is insensitive to χ. Further increasing V_in increases the contribution of ignition to flame stabilisation until the flame is completely ignition stabilised. Flow divergence caused by the LTC branch reduces the χ at the HTC branches making the HTC more resilient to χ. The results are discussed in the context of identification of possible stabilisation modes in turbulent flames.     

Flame front analysis of high-pressure turbulent lean premixed methane–air flames

An experimental study on lean turbulent premixed methane–air flames at high pressure is conducted by using a turbulent Bunsen flame configuration. A single equivalence ratio flame at Φ = 0.6 is explored for pressures ranging from atmospheric pressure to 0.9 MPa. LDA measurements of the cold flow indicate that turbulence intensities and the integral length scale are not sensitive to pressure. Due to the decreased kinematic viscosity with increasing pressure, the turbulent Reynolds numbers increase, and isotropic turbulence scaling relations indicate a large decrease of the smallest turbulence scales. Available experimental results and PREMIX code computations indicate a decrease in laminar flame propagation velocities with increasing pressure, essentially between the atmospheric pressure and 0.5 MPa. The u′/S_L ratio increases therefore accordingly. Instantaneous flame images are obtained by Mie scattering tomography. The images and their analysis show that pressure increase generates small scale flame structures. In an attempt to generalize these results, the variance of the flamelet curvatures, the standard deviation of the flamelet orientation angle, and the flamelet crossing lengths have been plotted against which is proportional to the ratio between the integral and Taylor length scales, and which increases with pressure. These three parameters vary linearly with the ratio between large and small turbulence scales and clearly indicate the strong effect of this parameter on premixed turbulent flame dynamics and structure. An obvious consequence is the increase in flame surface density and hence burning rate with pressure, as confirmed by its direct determination from 2D tomographic images.     

Issues arising in the construction of QSSA mechanisms: the case of reduced n-heptane/air models for premixed flames

A model reduction methodology, based on the quasi steady-state approximation (QSSA), is employed for the construction of reduced mechanisms in the case of an n-heptane/air premixed flame. Several issues related to the construction of these reduced mechanisms are discussed; such as the influence of the size of the starting skeletal mechanism, the stiffness reduction, and the truncation/simplification of (i) the expressions of the global rates and (ii) the steady-state relations. The starting point for the reduction is two skeletal mechanisms that involve 177/768 and 66/326 species/reactions, respectively [J. Prager, H.N. Najm, M. Valorani, and D.A. Goussis, Skeletal mechanism generation with CSP and validation for premixed n-heptane flames, Proc. Combust. Inst. 32 (2009), pp. 509–517] and which were derived from the detailed mechanism of Curran et al. [H.J. Curran, P. Gaffuri, W.J. Pitz, and C.K. Westbrook, A comprehensive modeling study of iso-octane oxidation, Combust. Flame 129 (2002), pp. 253–280], which involves 561/2538 species/reactions. From these two skeletal mechanisms, a number of reduced mechanisms of various sizes are produced and analysed. The validity of the reduced mechanism with the minimum size is demonstrated by considering its accuracy regarding the mass fractions of major and minor species, the temperature, and the flame speed, over a wide range of equivalence ratios and pressures.     

Turbulent flame propagation limits of ammonia/methane/air premixed mixture in a constant volume vessel

Ammonia is one of promising energy carriers that can be directly used as carbon-neutral fuel for combustion applications. However, because of the low-burning velocity of ammonia, it is challenging to introduce ammonia to practical combustors those are designed for general hydrocarbon fuels. One of ways to enhance the combustibility of ammonia is by mixing it with other hydrocarbon fuels, such as methane, with a burning velocity is much higher than the burning velocity of ammonia. In this study, we conducted flame propagation experiments of ammonia/methane/air using a fan-stirred constant volume vessel to clarify the effect of methane addition to ammonia on the turbulent flame propagation limit. From experimental results, we constructed the flame propagation maps and clarified the flame propagation limits. The results show that the flame propagation limits were extended with an increase in mixing a fraction of methane to ammonia. Additionally, ammonia/methane/air mixtures with the equivalence ration of 0.9 can propagate at the highest turbulent intensity, even though the peak of the laminar burning velocity is the fuel-rich side because of the diffusional-thermal instability of the flame surface. Furthermore, the Markstein number of the mixture obtained in this research successfully expressed the strength of the diffusional-thermal instability effect on the flame propagation capability. The turbulence Karlovitz number at the flame propagation limit monotonically increases with the decreasing Markstein number.