Effects of equivalence ratio variation on lean,stratified methane–air laminar counterflow flames

The effects of equivalence ratio variations on flame structure and propagation have been studied computationally. Equivalence ratio stratification is a key technology for advanced low emission combustors. Laminar counterflow simulations of lean methane–air combustion have been presented which show the effect of strain variations on flames stabilized in an equivalence ratio gradient, and the response of flames propagating into a mixture with a time-varying equivalence ratio. ‘Back supported’ lean flames, whose products are closer to stoichiometry than their reactants, display increased propagation velocities and reduced thickness compared with flames where the reactants are richer than the products. The radical concentrations in the vicinity of the flame are modified by the effect of an equivalence ratio gradient on the temperature profile and thermal dissociation. Analysis of steady flames stabilized in an equivalence ratio gradient demonstrates that the radical flux through the flame, and the modified radical concentrations in the reaction zone, contribute to the modified propagation speed and thickness of stratified flames. The modified concentrations of radical species in stratified flames mean that, in general, the reaction rate is not accurately parametrized by progress variable and equivalence ratio alone. A definition of stratified flame propagation based upon the displacement speed of a mixture fraction dependent progress variable was seen to be suitable for stratified combustion. The response times of the reaction, diffusion, and cross-dissipation components which contribute to this displacement speed have been used to explain flame response to stratification and unsteady fluid dynamic strain.     

Diffusive–thermal oscillations of rich premixed hydrogen–air flames in a microflow reactor

In this paper the dynamics of rich hydrogen–air flames in a microflow reactor with controlled temperature of the walls is investigated numerically using the thermal-diffusion model with two-step kinetics in one spatial dimension. It is found that as the parameters of the system are varied the sequence of bifurcation occurs leading to the formation of complex spatio-temporal patterns. These include pulsating, chaotic, mixed-mode and FREI (Flames with Repetitive Extinction and Ignition) oscillations. The critical parameter values for the existence of different dynamical regimes are found in terms of equivalence ratio and flow velocity.     

Measurement of unstable burning velocities of iso-octane–air mixtures at high pressure and the derivation of laminar burning velocities

A new technique is reported for measuring burning velocities at high pressures in the final stages of two inwardly propagating flame kernels in an explosion bomb. The flames were initiated at diametrically opposite spark electrodes, close to the wall, in quiescent mixtures. Measurements of pressure and flame kernel propagation speeds by high-speed photography showed the burning velocities to be elevated above the corresponding laminar burning velocities as a result of the developing flame instabilities. The enhancement increased with increase in pressure and decreased with increase in Markstein number. When the Markstein number was negative, instabilities could be appreciable, as could the enhancement. For the iso-octane–air mixtures investigated, where the mixtures had well-characterised Markstein numbers or critical Peclet numbers at the relevant pressures and temperatures, it was possible to explain the enhancement quantitatively by the spherical explosion flame instability theory of Bechtold and Matalon, provided the critical Peclet number was that observed experimentally, and allowance was made for the changing pressure. With this theoretical procedure, it was possible to derive values of laminar burning velocity from the measured values of burning velocity over a wide range of equivalence ratios, pressures, and temperatures. The values became less reliable at the higher temperatures and pressures as the data on Markstein and critical Peclet numbers became less certain. It was found that with iso-octane as the fuel the laminar burning velocity decreased during isentropic compression.     

Experimental and numerical study of product gas characteristics of ammonia/air premixed laminar flames stabilized in a stagnation flow

Ammonia has widely attracted interest as a potential candidate not only as a hydrogen energy carrier but also as a carbon free fuel for internal combustion engines, such as gas turbines. Because ammonia contains a nitrogen atom in its molecule, nitrogen oxides (NO_x) and other pollutants may be formed when it burns. Therefore, understanding the fundamental product gas characteristics of ammonia/air laminar flames is important for the design of ammonia-fueled combustors to meet stringent emission regulations. In this study, the product gas characteristics of ammonia/air premixed laminar flames for various equivalence ratios were experimentally and numerically investigated up to elevated pressure conditions. In the experiments, a stagnation flame configuration was employed because an ammonia flame can be stabilized by using such a configuration without a pilot flame. The experimental results showed that the maximum NO mole fraction was about 3,500 ppmv, at an equivalence ratio of 0.9 at 0.1 MPa. The NO mole fraction decreased as the equivalence ratio increased. In addition, the maximum value of the NO mole fraction decreased with an increase in mixture pressure. Furthermore, it was experimentally clarified that the simultaneous reduction of NO and unburnt ammonia can be achieved at an equivalence ratio of about 1.06, which is the target equivalence ratio for emission control in rich-lean two-stage ammonia combustors. Comparison of experimental and numerical results showed that even though the reaction mechanisms employed have been optimized for predicting the laminar burning velocity of ammonia/air flames, they failed to satisfactorily predict the measured species in this study. Sensitivity analysis was used to identify elementary reactions that control the species profiles but have negligible effects on the burning velocity. It is considered that these reaction models need to be updated for accurate prediction of product gas characteristics of ammonia/air flames.     

Characterization of the effects of pressure and hydrogen concentration on laminar burning velocities of methane–hydrogen–air mixtures

The aim of the present work was to characterize both the effects of pressure and of hydrogen addition on methane/air premixed laminar flames. The experimental setup consists of a spherical combustion chamber coupled to a classical shadowgraphy system. Flame pictures are recorded by a high speed camera. Global equivalence ratios were varied from 0.7 to 1.2 for the initial pressure range from 0.1 to 0.5 MPa. The mole fraction of hydrogen in the methane + hydrogen mixture was varied from 0 to 0.2. Experimental results were compared to calculations using a detailed chemical kinetic scheme (GRIMECH 3.0). First, the results for atmospheric laminar CH₄/air flames were compared to the literature. Very good agreements were obtained both for laminar burning velocities and for burned gas Markstein length. Then, increasing the hydrogen content in the mixture was found to be responsible for an increase in the laminar burning velocity and for a reduction of the flame dependence on stretch. Transport effects, through the reduction of the fuel Lewis number, play a role in reducing the sensitivity of the fundamental flame velocity to the stretch. Finally, when the pressure was increased, the laminar burning velocity decreased for all mixtures. The pressure domain was limited to 0.5 MPa due to the onset of instabilities at pressures above this value.     

Diffusional-thermal instability in strained diffusion flames with unequal Lewis numbers

Diffusional–thermal instability is analysed for near-extinction counterflow diffusion flames to examine the instability characteristics of strained diffusion flamelets in turbulent flames, with the additional intention of providing a guideline to future experimental investigations. Attention is focused on the linear stability of the instability patterns appearing in the unstrained direction of two-dimensional counterflow diffusion flames, which is treated by the near-equilibrium regime of activation-energy asymptotics with Lewis numbers close to unity. The effects of unequal Lewis numbers for fuel and oxidizer are also taken into account by introducing an effective Lewis number. The resulting formulation describing linear stability of the harmonically decomposed disturbances turns out to be identical to the formulation derived previously for equal fuel and oxidizer Lewis numbers. For effective Lewis numbers less than unity, cellular instability is predicted for the entire range of the equivalence ratio, and the threshold Lewis number maintains a value slightly less than unity. On the other hand, for effective Lewis numbers sufficiently greater than unity, two types of oscillatory instabilities are found. As the effective Lewis number increases from unity, a travelling instability is first encountered for a range of finite wavelengths, and a pulsating instability emerges immediately above the travelling instability. These two types of oscillatory instabilities are predicted only for equivalence ratios sufficiently greater than unity because the threshold Lewis numbers for these instabilities are found to be infinity at unity equivalence ratio. For large values of the equivalence ratio, which is typical of most hydrocarbon flames, oscillatory instabilities are predicted only for flames burning extremely heavy hydrocarbon fuels or for flames heavily diluted by light inert gases.     

Effects of preferential and differential diffusion on the mutual annihilation of two premixed hydrogen–air flames

The unsteady process of upstream head-on quenching of two laminar premixed hydrogen–air flames at different equivalence ratios in one dimension is investigated numerically in the presence of preferential and differential diffusion effects. Important chemical and transport characteristics of the mutual annihilation process are studied during the two primary stages of upstream mutual annihilation, preheat layers' and reaction layers' interactions. Because of the diffusive mobility of the fuel, hydrogen, relative to heat and the oxidizer, preferential and differential diffusion effects result in a shift in the equivalence ratio in the reaction zone to leaner conditions. This shift, in turn, affects the subsequent reaction layers' interactions through qualitative and quantitative changes in the rates of reactants' consumption and radicals' production. Another consequence of this shift is the presence of excess and ‘unburnt’ fuel or oxidizer at the end of the mutual annihilation process. The process of mutual annihilation occurs over time scales that are significantly shorter than characteristic residence times associated with flames.     

Laminar and unstable burning velocities and Markstein lengths of hydrogen–air mixtures at engine-like conditions

Hydrogen offers an attractive alternative to conventional fuels for use in spark ignition engines. It can be burned over a very wide range of equivalence ratios and with considerable exhaust gas recirculation. These help to minimise pumping losses through throttleless operation and oxides of nitrogen (NO_x) production through reduced temperature. Full understanding of hydrogen-fuelled engine operation requires data on the laminar burning rate of hydrogen–air residuals under a wide range of conditions. However, such data are sparse. The present work addresses this need for experimental data. Spherically expanding H₂–air flames were measured at a range of temperatures, pressures, and equivalence ratios and with varying concentrations of residuals of combustion. Unstretched burning velocities, u_l, and Markstein lengths, L_b, were determined from stable flames. At the higher pressures, hydrodynamic and diffusional-thermal instabilities caused the flames to be cellular from inception and prohibited the derivation of values of u_l and L_b. The effect of pressure on the burning rate was demonstrated to have opposing trends when comparing stoichiometric and lean mixtures. The present measurements were compared with those available in the literature, and discrepancies were attributed to neglect, in some works, the effects of stretch and instabilities. From the present measurements, the effects of pressure, temperature, and residual gas concentration on burning velocity are quantified for use in a first step towards a general correlation.     

Modelling of turbulent scalar flux in turbulent premixed flames based on DNS databases

A transport equation for scalar flux in turbulent premixed flames was modelled on the basis of DNS databases. Fully developed turbulent premixed flames were obtained for three different density ratios of flames with a single-step irreversible reaction, while the turbulent intensity was comparable to the laminar burning velocity. These DNS databases showed that the countergradient diffusion was dominant in the flame region. Analyses of the Favre-averaged transport equation for turbulent scalar flux proved that the pressure related terms and the velocity–reaction rate correlation term played important roles on the countergradient diffusion, while the mean velocity gradient term, the mean progress variable gradient term and dissipation terms suppressed it. Based on these analyses, modelling of the combustion-related terms was discussed. The mean pressure gradient term and the fluctuating pressure term were modelled by scaling, and these models were in good agreement with DNS databases. The dissipation terms and the velocity–reaction rate correlation term were also modelled, and these models mimicked DNS well.     

Premixed turbulent flame front structure investigation by Rayleigh scattering in the thin reaction zone regime

Premixed turbulent flames of methane–air and propane–air stabilized on a bunsen type burner were studied using planar Rayleigh scattering and particle image velocimetry. The fuel–air equivalence ratio range was from lean 0.6 to stoichiometric for methane flames, and from 0.7 to stoichiometric for propane flames. The non-dimensional turbulence rms velocity, u′/S_L, covered a range from 3 to 24, corresponding to conditions of corrugated flamelets and thin reaction zones regimes. Flame front thickness increased slightly with increasing non-dimensional turbulence rms velocity in both methane and propane flames, although the flame thickening was more prominent in propane flames. The probability density function of curvature showed a Gaussian-like distribution at all turbulence intensities in both methane and propane flames, at all sections of the flame.The value of the term Dκ, the product of molecular diffusivity evaluated at reaction zone conditions and the flame front curvature, has been shown to be smaller than the magnitude of the laminar burning velocity. This finding questions the validity of extending the level set formulation, developed for corrugated flames region, into the thin reaction zone regime by increasing the local flame propagation by adding the term Dκ to laminar burning velocity.     

Experimental measurements of two-dimensional planar propagating edge flames

The study of edge flames has received increased attention in recent years. This work reports the results of a recent study into two-dimensional, planar, propagating edge flames that are remote from solid surfaces (called here, “free-layer” flames, as opposed to layered flames along floors or ceilings). They represent an ideal case of a flame propagating down a flammable plume, or through a flammable layer in microgravity. The results were generated using a new apparatus in which a thin stream of gaseous fuel is injected into a low-speed laminar wind tunnel thereby forming a flammable layer along the centerline. An airfoil-shaped fuel dispenser downstream of the duct inlet issues ethane from a slot in the trailing edge. The air and ethane mix due to mass diffusion while flowing up towards the duct exit, forming a flammable layer with a steep lateral fuel concentration gradient and smaller axial fuel concentration gradient. We characterized the flow and fuel concentration fields in the duct using hot wire anemometer scans, flow visualization using smoke traces, and non-reacting, numerical modeling using COSMOSFloWorks. In the experiment, a hot wire near the exit ignites the ethane-air layer, with the flame propagating downwards towards the fuel source. Reported here are tests with the air inlet velocity of 25 cm/s and ethane flows of 967-1299 sccm, which gave conditions ranging from lean to rich along the centerline. In these conditions the flame spreads at a constant rate faster than the laminar burning rate for a premixed ethane-air mixture. The flame spread rate increases with increasing transverse fuel gradient (obtained by increasing the fuel flow rate), but appears to reach a maximum. The flow field shows little effect due to the flame approach near the igniter, but shows significant effect, including flow reversal, well ahead of the flame as it approaches the airfoil fuel source.     

Effect of geometrical parameters on thermo-acoustic instability of downward propagating flames in tubes

Experiments and theoretical analysis are presented to clarify the effect of geometrical parameters on thermo-acoustic instability of downward propagating flames in tubes. The experiments reveal that the longer tubes have higher instability compared to shorter tubes and the lower diameter tubes have higher instability compared to higher diameter tubes. The secondary instability leading to turbulent burning is found to be more sensitive to change in geometrical parameters compared to primary instability (oscillating flat flame). The secondary instability is re-stabilized for some intermediate burning velocity conditions even though lower and higher burning velocity conditions show secondary instability. The appearance of such re-stabilization is only observed for some specific lengths of the tube. Present experimental observations pertaining to the effect of geometrical parameters is found to be contradicting the theoretical predictions based on pressure coupling mechanism. To clear the underlying mechanism, analytical growth rate is computed considering velocity coupling mechanism. The computed growth rates correctly predict the effect of geometrical parameters on thermo-acoustic instability of downward propagating flames. This work provides further evidence to believe that the flame -acoustic coupling in downward propagating flames is due to flame area modulation (leading to heat release modulation) through action of acoustic acceleration.     

Insight into fuel isomeric effects on laminar flame propagation of pentanones

Laminar flame propagation was investigated for pentanone isomers/air mixtures (3-pentanone, 2-pentanone and 3-methyl-2-butanone) in a high-pressure constant-volume cylindrical combustion vessel at 393–423 K, 1–10 atm and equivalence ratios of 0.6–1.5, and in a heat flux burner at 393 K, 1 atm and equivalence ratios of 0.6–1.5. Two kinds of methods generally show good agreement, both of which indicate that the laminar burning velocity increases in the order of 3-methyl-2-butanone, 2-pentanone and 3-pentanone. A kinetic model of pentanone isomers was developed and validated against experimental data in this work and in literature. Modeling analysis was performed to provide insight into the flame chemistry of the three pentanone isomers. H-abstraction reactions are concluded to dominate fuel consumption, and further decomposition of fuel radicals eventually produces fuel-specific small radicals. The differences in radical pools are concluded to be responsible for the observed fuel isomeric effects on laminar burning velocity. Among the three pentanone isomers, 3-pentanone tends to produce ethyl and does not prefer to produce methyl and allyl in flames, thus it has the highest reactivity and fastest laminar flame propagation. On the contrary, 3-methyl-2-butanone tends to produce allyl and methyl instead of ethyl, and consequently has the lowest reactivity and slowest laminar flame propagation.     

A flamelet analysis of the burning velocity of premixed turbulent expanding flames

Direct numerical simulation is a very powerful tool to evaluate the validity of new models and theories for turbulent combustion. In this paper, direct numerical simulations of spherically expanding premixed turbulent flames in the corrugated flamelet regime are performed. The flamelet-generated manifold method is used to deal with detailed reaction kinetics. The numerical method is validated for both laminar and turbulent expanding flames. The computational results are analyzed by using an extended flame stretch theory. It is investigated whether this theory is able to describe the influence of flame stretch and curvature on the local burning velocity of the flame. If the full profiles of flame stretch and curvature through the flame front are included in the theory, the local mass burning rate is predicted accurately. The influence of several approximations, which are used in other existing theories, is studied. When flame stretch is assumed to be constant through the flame front or when curvature of the flame front is neglected, the theory fails to predict the local mass burning rate.     

The role of radiative losses in self-extinguishing and self-wrinkling flames

We have developed a general theory of non-adiabatic premixed flames that is valid for flames of arbitrary shape that fully accounts for the hydrodynamic and diffusive-thermal processes, and incorporates the effects of volumetric heat losses. The model is used to describe aspects of experimentally observed phenomena of self-extinguishing (SEFs) and self-wrinkling flames (SWFs), in which radiative heat losses play an important role. SEFs are spherical flames that propagate considerable distances in sub-limit conditions before suddenly extinguishing. Our results capture many aspects of this phenomenon including an explicit determination of flame size and propagation speed at quenching. SWFs are hydrodynamically unstable flames in which cells spontaneously appear on the flame surface once the flame reaches a critical size. Our results yield expressions of the critical flame size at the onset of wrinkling and expected cell size beyond the stability threshold. The various possible burning regimes are mapped out in parameter space.     

近贫燃极限甲烷/空气火焰的层流燃烧速度研究

利用微重力条件下向外传播的球形火焰,对贫燃极限附近甲烷/空气预混火焰的层流燃烧速度进行了测量,得到当量比从0.512(本文微重力实验中测定的可燃极限)到0.601范围内的零拉伸层流燃烧速度,并与前人实验数据和使用3种化学反应动力学模型的计算结果进行了比较.本文实验结果与已有的微重力实验数据非常接近,而其他研究者在常重力...     

Theoretical and experimental studies on mesoscale flame propagation and extinction

Mesoscale flame propagation and extinction of premixed flames in channels are investigated theoretically and experimentally. Emphasis is placed on the effect of wall heat loss and the wall–flame interaction via heat recirculation. At first, an analytical solution of flame speed in mesoscale channels is obtained. The results showed that channel width, flow velocity, and wall thermal properties have dramatic effects on the flame propagation and lead to multiple flame regimes and extinction limits. With the decrease in channel width, there exist two distinct flame regimes, a fast burning regime and a slow burning regime. The existence of the new flame regime and its extended flammability limit render the classical quenching diameter inapplicable. Furthermore, the results showed that at optimum conditions of flow velocity and wall thermal properties, mesoscale flames can propagate faster than the adiabatic flame. Second, numerical simulation with detailed chemistry demonstrated the existence of multiple flame regimes. The results also showed that there is a non-linear dependence of the flame speed on equivalence ratio. Moreover, it is shown that the Nusselt number has a significant impact on this non-linear dependence. Finally, the non-linear dependence of flame speed on equivalence ratio for both flame regimes is measured using a C₃H₈–air mixture. The results are in good agreement with the theory and numerical simulation.     

Cellular instabilities of expanding hydrogen/propane spherical flames at elevated pressures: theory and experiment

An experimental and theoretical investigation of the onset of cellular instabilities on spherically expanding flames in mixtures of hydrogen and propane in air at elevated pressures was conducted. Critical conditions for the onset of instability were measured and mapped out over a range of pressures and mixture compositions. An asymptotic theory of hydrodynamic and diffusional-thermal cell development on flames in mixtures comprised of two scarce fuels burning in air was also formulated. Predicted values of Peclet number, defined as the flame radius at the onset of instability normalized by the flame thickness, were shown to compare favorably with the experimentally measured values.