Propagation of nonstationary curved and stretched premixed flames with multistep reaction mechanisms

The propagation speed of a thin premixed flame disturbed by an unsteady fluid flow of a larger scale is considered. The flame may also have a general shape but the reaction zone is assumed to be thin compared to the flame thickness. Unlike in preceding publications, the presented asymptotic analysis is performed for a general multistep reaction mechanism and, at the same time, the flame front is curved by the fluid flow. The resulting equations define the propagation speed of disturbed flames in terms of the properties of undisturbed planar flames and the flame stretch. Special attention is paid to the near-equidiffusion limit. In this case, the flame propagation speed is shown to depend on the effective Zeldovich number Z(f) , and the flame stretch. Unlike the conventional Zeldovich number, the effective Zeldovich number is not necessarily linked directly to the activation energies of the reactions. Several examples of determining the effective Zeldovich number for reduced combustion mechanisms are given while, for realistic reactions, the effective Zeldovich number is determined from experiments. Another feature of the present approach is represented by the relatively simple asymptotic technique based on the adaptive generalized curvilinear system of coordinates attached to the flame (i.e., intrinsic disturbed flame equations [IDFE]).

NOx and CO emissions and soot presence in partially premixed acoustically excited flames

The pulsating combustion process has attracted interest in current research because its application in energy generation can offer several advantages, such as fuel economy, reduced pollutants formation, increased rate of convective heat transfer and reduced investment, when compared with other new techniques of combustion. An experimental study has been conducted with the objective of investigating the effects of combustion driven acoustic oscillations in the emission rates of combustion gases, especially carbon monoxide and nitrogen oxides, and soot presence in partial premixed flames in confined partially premixed liquefied petroleum gas flames. The results basically showed that a more uniform fuel/air mixture due to the presence of an acoustic field increases the NOx emissions in operations close to stoichiometric equivalence ratios and the frequency is the most important parameter. Carbon monoxide and soot reduced significantly.

A unified model of flames as gasdynamic discontinuities

Viewed on a hydrodynamic scale, flames in experiments are often thin so that they may be described as gasdynamic discontinuities separating the dense cold fresh mixture from the light hot burned products. The original model of a flame as a gasdynamic discontinuity was due to Darrieus and to Landau. In addition to the fluid dynamical equations, the model consists of a flame speed relation describing the evolution of the discontinuity surface, and jump conditions across the surface which relate the fluid variables on the two sides of the surface. The Darrieus-Landau model predicts, in contrast to observations, that a uniformly propagating planar flame is absolutely unstable and that the strength of the instability grows with increasing perturbation wavenumber so that there is no high-wavenumber cutoff of the instability. The model was modified by Markstein to exhibit a high-wavenumber cutoff if a phenomenological constant in the model has an appropriate sign. Both models are postulated, rather than derived from first principles, and both ignore the flame structure, which depends on chemical kinetics and transport processes within the flame. At present, there are two models which have been derived, rather than postulated, and which are valid in two non-overlapping regions of parameter space. Sivashinsky derived a generalization of the Darrieus-Landau model which is valid for Lewis numbers (ratio of thermal diffusivity to mass diffusivity of the deficient reaction component) bounded away from unity. Matalon & Matkowsky derived a model valid for Lewis numbers close to unity. Each model has its own advantages and disadvantages. Under appropriate conditions the Matalon-Matkowsky model exhibits a high-wavenumber cutoff of the Darrieus-Landau instability. However, since the Lewis numbers considered lie too close to unity, the Matalon-Matkowsky model does not capture the pulsating instability. The Sivashinsky model does capture the pulsating instability, but does not exhibit its high-wavenumber cutoff. In this paper, we derive a model consisting of a new flame speed relation and new jump conditions, which is valid for arbitrary Lewis numbers. It captures the pulsating instability and exhibits the high-wavenumber cutoff of all instabilities. The flame speed relation includes the effect of short wavelengths, not previously considered, which leads to stabilizing transverse surface diffusion terms.

Stability of planar flames as gasdynamic discontinuities

The stability of a steadily propagating planar premixed flame has been the subject of numerous studies since Darrieus and Landau showed that in their model flames are unstable to perturbations of any wavelength. Moreover, the instability was shown to persist even for very small wavelengths, i.e. there was no high-wavenumber cutoff of the instability. In addition to the Darrieus-Landau instability, which results from thermal expansion, analysis of the diffusional thermal model indicates that premixed flames may exhibit cellular and pulsating instabilities as a consequence of preferential diffusion. However, no previous theory captured all the instabilities including a high-wavenumber cutoff for each. In Class, Matkowsky & Klimenko (2003) a unified theory is proposed which, in appropriate limits and under appropriate assumptions, recovers all the relevant previous theories. It also includes additional new terms, not present in previous theories. In the present paper we consider the stability of a uniformly propagating planar flame as a solution of the unified model. The results are then compared to those based on the models of Darrieus-Landau, Sivashinsky and Matalon-Matkowsky. In particular, it is shown that the unified model is the only model to capture the Darrieus-Landau, cellular and pulsating instabilities including a high-wavenumber cutoff for each.

Flammability conditions for ultra-lean hydrogen premixed combustion based on flame-ball analyses

It has been reasoned that the structures of strongly cellular flames in very lean mixtures approach an array of flame balls, each burning as if it were isolated, thereby indicating a connection between the critical conditions required for existence of steady flame balls and those necessary for occurrence of self-sustained premixed combustion. This is the starting assumption of the present study, in which structures of near-limit steady sphericosym-metrical flame balls are investigated with the objective of providing analytic expressions for critical combustion conditions in ultra-lean hydrogen-oxygen mixtures diluted with N2 and water vapor. If attention were restricted to planar premixed flames, then the lean-limit mole fraction of H2 would be found to be roughly ten percent, more than twice the observed flammability limits, thereby emphasizing the relevance of the flame-ball phenomena. Numerical integrations using detailed models for chemistry and radiation show that a onestep chemical-kinetic reduced mechanism based on steady-state assumptions for all chemical intermediates, together with a simple, optically thin approximation for water-vapor radiation, can be used to compute near-limit fuel-lean flame balls with excellent accuracy. The previously developed one-step reaction rate includes a crossover temperature that determines in the first approximation a chemical-kinetic lean limit below which combustión cannot occur, with critical conditions achieved when the diffusion-controlled radiation-free peak temperature, computed with account taken of hydrogen Soret diffusion, is equal to the crossover temperature. First-order corrections are found by activation-energy asymptotics in a solution that involves a near-field radiation-free zone surrounding a spherical flame sheet, together with a far-field radiation-conduction balance for the temperature profile. Different scalings are found depending on whether or not the surrounding atmosphere contains wáter vapor, leading to different analytic expressions for the critical conditions for flame-ball existence, which give results in very good agreement with those obtained by detailed numerical computations.

Local burning velocity in a Bunsen jet flame

A PIV-based system has been set-up for the simultaneous measurement of the local burning velocity of premixed flames and the flame stretch due to the flame front curvature and the incoming flow strain rate. For moderately short jet flames, these measurements allow an indirect determination of the Markstein length, according to Clavin and Joulin (C–J) theory. For tall flames, the flame curvature becomes relatively large in a region around the tip where the C–J theory breaks down. However, our experiments confirm the appearance of a new linear relation between burning velocity and curvature at the flame tip. This relation defines a new proportionality factor which is probably associated to the evolution from rounded tips to slender tips when the jet velocity is increased.

Emissions from the premixed combustion of gasified polyethylene

An investigation was conducted on pollutants emitted from steady-state, steady-flow gasification and combustion of polyethylene (PE) in a two-stage furnace. The polymer, in pulverized form, was first pyrolyzed at 1000 degrees C, and subsequently, its gaseous pyrolyzates were burned, upon mixing with air at high temperatures (900-1100 degrees C). The motivation for this indirect type of burning PE was to attain nominally premixed combustion of the pyrolyzate gases with air, thereby achieving lower pollutant emissions than those emanating from the direct burning of the solid PE polymer. This work assessed the effluents of the two-stage furnace and examined the effects of the combustion temperature, as well as the polymer feed rate and the associated fuel/air equivalence ratio (0.3 < phi < 1.4). It was found that, whereas the yield of pyrolysis gas decreased with an increasing polymer feed rate, its composition was nearly independent of the feed rate. CO2 emissions peaked at an equivalence ratio near unity, while the CO emissions increased with an increasing equivalence ratio. The total light volatile hydrocarbon and semivolatile polycyclic aromatic hydrocarbon (PAH) emissions of combustion increased with an increasing equivalence ratio. The generated particulates were mostly submicrometer in size. Overall, PAH and soot emissions from this indirect burning of PE were an order of magnitude lower than corresponding emissions from the direct burning of the solid polymer, obtained previously in this laboratory using identical sampling and analytical techniques. Because pyrolysis of this polymer requires a nominal heat input that amounts to only a diminutive fraction of the heat released during its combustion, implementation of this technique is deemed advantageous.

Temperature measurements by oh lif and chemiluminescence kinetic modeling for ethanol flames

OH LIF-thermometry was applied to premixed ethanol flames at atmospheric pressure in a burner for three flame conditions. Flame temperatures were simulated from energy equation with PREMIX code of CHEMKIN software package for comparison. A kinetic modeling based on a model validated through chemiluminescence measurements and on a set of reactions for nitrogen chemistry was evaluated. Marinov's mechanism was also tested. Sensitivity analysis was performed for fuel-rich flame condition with Φ = 1.34. Simulated temperatures from both reaction mechanisms evaluated were higher than experimental values. However, the proposed kinetic modeling resulted in temperature profiles qualitatively very close to the experimental.

Gas concentration and temperature excited Delft turbulent jet in acoustically flames

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

CH and C2 radicals characterization in natural gas turbulent diffusion flames

This paper reports the construction of an axisymmetric nonpremixed piloted jet burner, with well-defined initial and boundary conditions, known as the Delft burner, to assess turbulence-chemistry interaction in non-premixed turbulent flames. Detailed experimental information is described, involving hot-wire anemometry, thin-wire thermocouples and chemiluminescence visualization measurements. Radial profile of the axial mean velocity indicates excellent agreement between flow patterns developed within Delft installation and the one described herein. Chemiluminescence emissions from CH and C2 free-radicals were acquired with a CCD camera. Tomography reconstruction analysis was utilised to compare radical emissions and temperature spatial distributions. There was a strong dependence between temperature and CH/C 2 emissions. This is an indication that these radicals can be used in flame front studies.

Experimental aspects of partially premixed pulsating combustion

The pulsating combustion process has won interest in current research due to indications that its application in energy generation can offer several advantages, such as: fuel economy, reduced pollutants formation, increased rate of convective heat transfer and reduced investment, when compared with conventional techniques. An experimental study has been conducted with the objective of investigating the effects of combustion driven acoustic oscillations in the emission rates of combustion gases, especially carbon monoxide and nitrogen oxides. The experiments were conducted in a water-jacketed 1-m long by 25-cm internal diameter stainless steel vertical tube. The combustor operated with liquefied petroleum gas (LPG) in both oscillatory and non oscillatory conditions, under the same input conditions. Part of the reactant mixture was excited acoustically, before the burner exit, by a speaker positioned strategically. The burner was aligned with the chamber longitudinal axis and positioned at its bottom. The experiments were conducted for 0.16 g/s of LPG burning in stoichiometric equivalence ratio. The main conclusions were: a) the pulsating combustion process produces more uniform fuel/air profile than the non pulsating process, b) close to stoichiometric equivalence ratio the pulsating combustion process generates higher rates of NO x; c) the frequency has a strong influence in NO x emission, but the pressure amplitude has a weak influence; d) the presence of the acoustic field may change drastically the combustion gas emissions in diffusion flames, but in pre-mixed flames the influence is not as strong.

Dynamics of thermal ignition of spray flames in mixing layers

Conditions are identified under which analyses of laminar mixing layers can shed light on aspects of turbulent spray combustion. With this in mind, laminar spray-combustion models are formulated for both non-premixed and partially premixed systems. The laminar mixing layer separating a hot-air stream from a monodisperse spray carried by either an inert gas or air is investigated numerically and analytically in an effort to increase understanding of the ignition process leading to stabilization of high-speed spray combustion. The problem is formulated in an Eulerian framework, with the conservation equations written in the boundary-layer approximation and with a one-step Arrhenius model adopted for the chemistry description. The numerical integrations unveil two different types of ignition behaviour depending on the fuel availability in the reaction kernel, which in turn depends on the rates of droplet vaporization and fuel-vapour diffusion. When sufficient fuel is available near the hot boundary, as occurs when the thermochemical properties of heptane are employed for the fuel in the integrations, combustion is established through a precipitous temperature increase at a well-defined thermal-runaway location, a phenomenon that is amenable to a theoretical analysis based on activation-energy asymptotics, presented here, following earlier ideas developed in describing unsteady gaseous ignition in mixing layers. By way of contrast, when the amount of fuel vapour reaching the hot boundary is small, as is observed in the computations employing the thermochemical properties of methanol, the incipient chemical reaction gives rise to a slowly developing lean deflagration that consumes the available fuel as it propagates across the mixing layer towards the spray. The flame structure that develops downstream from the ignition point depends on the fuel considered and also on the spray carrier gas, with fuel sprays carried by air displaying either a lean deflagration bounding a region of distributed reaction or a distinct double-flame structure with a rich premixed flame on the spray side and a diffusion flame on the air side. Results are calculated for the distributions of mixture fraction and scalar dissipation rate across the mixing layer that reveal complexities that serve to identify differences between spray-flamelet and gaseous-flamelet problems.

Modeling of Turbulent Natural Gas and Biogas Flames of the Delft Jet-in-Hot-Coflow Burner: Effects of Coflow Temperature, Fuel Temperature and Fuel Composition on the Flame Lift-Off Height

The structure of a turbulent non-premixed flame of a biogas fuel in a hot and diluted coflow mimicking moderate and intense low dilution (MILD) combustion is studied numerically. Biogas fuel is obtained by dilution of Dutch natural gas (DNG) with CO2. The results of biogas combustion are compared with those of DNG combustion in the Delft Jet-in-Hot-Coflow (DJHC) burner. New experimental measurements of lift-off height and of velocity and temperature statistics have been made to provide a database for evaluating the capability of numerical methods in predicting the flame structure. Compared to the lift-off height of the DNG flame, addition of 30 % carbon dioxide to the fuel increases the lift-off height by less than 15 %. Numerical simulations are conducted by solving the RANS equations using Reynolds stress model (RSM) as turbulence model in combination with EDC (Eddy Dissipation Concept) and transported probability density function (PDF) as turbulence-chemistry interaction models. The DRM19 reduced mechanism is used as chemical kinetics with the EDC model. A tabulated chemistry model based on the Flamelet Generated Manifold (FGM) is adopted in the PDF method. The table describes a non-adiabatic three stream mixing problem between fuel, coflow and ambient air based on igniting counterflow diffusion flamelets. The results show that the EDC/DRM19 and PDF/FGM models predict the experimentally observed decreasing trend of lift-off height with increase of the coflow temperature. Although more detailed chemistry is used with EDC, the temperature fluctuations at the coflow inlet (approximately 100K) cannot be included resulting in a significant overprediction of the flame temperature. Only the PDF modeling results with temperature fluctuations predict the correct mean temperature profiles of the biogas case and compare well with the experimental temperature distributions.