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.

Numerical simulation of the upward propagation of a flame in a vertical tube filled with a very lean mixture.

Upwardpropagation of a premixed flame in averticaltubefilled with a very leanmixture is simulated numerically using a single irreversible Arrhenius reaction model with infinitely high activation energy. In the absence of heat losses and preferential diffusion effects, a curved flame with stationary shape and velocity close to those of an open bubble ascending in the same tube is found for values of the fuel mass fraction above a certain minimum that increases with the radius of the tube, while the numerical computations cease to converge to a stationary solution below this minimum mass fraction. The vortical flow of the gas behind the flame and in its transport region is described for tubes of different radii. It is argued that this flow may become unstable when the fuel mass fraction is decreased, and that this instability, together with the flame stretch due to the strong curvature of the flame tip in narrow tubes, may be responsible for the minimum fuel mass fraction. Radiation losses and a Lewis number of the fuel slightly above unity decrease the final combustion temperature at the flame tip and increase the minimum fuel mass fraction, while a Lewis number slightly below unity has the opposite effect.

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.

On the formation of nitrogen oxides in the combustion processes of hydrogen in air

A research programme is being carried out at the Institute Nacional de Tecnica Aeroespacial of Spain, on several aspects of the formation of nitrogen oxides in continuous flow combustion systems, considering hydrogen and hydrocarbons as fuels. The research programme is fundamentally oriented on the basic aspects of the problem, although it also includes the study of the influence on the formation process of several operational and design variables of the combusters, such as type of fuels, fuel/air ratio, degree of mixing in premixed type flames, existence of droplets as compared with homogeneous combustion.This problem of nitrogen oxides formation is receiving lately great attention, specially in connection with automobile reciprocating engines and aircraft gas turbines. This is due to the fact of the increasing frequency and intensity of photochemical hazes or smog, typical of urban areas submitted to strong solar radiation, which are originated by the action on organic compounds of the oxidants resulting from the photochemical decomposition of nitrogen dioxide N02. In the combustion process almost all nitrogen oxides are in form of NO. This nitric oxide reacts with the oxygen of the air and forms N02, this reaction only taking place in or near the exhaust of tne motors, since the N0-02 reaction becomes frozen for the concentration existing in the atmosphere.

Numerical analysis of deflagration initation by a hot jet

Numerical simulations of axisymmetric reactive jets with one-step Arrhenius kinetics are used to investigate the problem of deflagration initiation in a premixed fuel–air mixture by the sudden discharge of a hot jet of its adiabatic reaction products. For the moderately large values of the jet Reynolds number considered in the computations, chemical reaction is seen to occur initially in the thin mixing layer that separates the hot products from the cold reactants. This mixing layer is wrapped around by the starting vortex, thereby enhancing mixing at the jet head, which is followed by an annular mixing layer that trails behind, connecting the leading vortex with the orifice rim. A successful deflagration is seen to develop for values of the orifice radius larger than a critical value a c in the order of the flame thickness of the planar deflagration δL. Introduction of appropriate scales provides the dimensionless formulation of the problem, with flame initiation characterised in terms of a critical Damköhler number Δc=(a d/δL)2, whose parametric dependence is investigated. The numerical computations reveal that, while the jet Reynolds number exerts a limited influence on the criticality conditions, the effect of the reactant diffusivity on ignition is much more pronounced, with the value of Δc increasing significantly with increasing Lewis numbers. The reactant diffusivity affects also the way ignition takes place, so that for reactants with the flame develops as a result of ignition in the annular mixing layer surrounding the developing jet stem, whereas for highly diffusive reactants with Lewis numbers sufficiently smaller than unity combustion is initiated in the mixed core formed around the starting vortex. The analysis provides increased understanding of deflagration initiation processes, including the effects of differential diffusion, and points to the need for further investigations corporating detailed chemistry models for specific fuel–air mixtures.

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.

Source terms for calculations of vaporizing and burning fuel sprays with non-unity Lewis numbers in gases with temperature-dependent thermal conductivities

Liquid-fueled burners are used in a number of propulsion devices ranging from internal combustion engines to gas turbines. The structure of spray flames is quite complex and involves a wide range of time and spatial scales in both premixed and non-premixed modes (Williams 1965; Luo et al. 2011). A number of spray-combustion regimes can be observed experimentally in canonical scenarios of practical relevance such as counterflow diffusion flames (Li 1997), as sketched in figure 1, and for which different microscalemodelling strategies are needed. In this study, source terms for the conservation equations are calculated for heating, vaporizing and burning sprays in the single-droplet combustion regime. The present analysis provides extended formulation for source terms, which include non-unity Lewis numbers and variable thermal conductivities.

Flamelet structures in spray ignition

In typical liquid-fueled burners the fuel is injected as a high-velocity liquid jet that breaks up to form the spray. The initial heating and vaporization of the liquid fuel rely on the relatively large temperatures of the sourrounding gas, which may include hot combustion products and preheated air. The heat exchange between the liquid and the gas phases is enhanced by droplet dispersion arising from the turbulent motion. Chemical reaction takes place once molecular mixing between the fuel vapor and the oxidizer has occurred in mixing layers separating the spray flow from the hot air stream. Since in most applications the injection velocities are much larger than the premixed-flame propagation velocity, combustion stabilization relies on autoignition of the fuel-oxygen mixture, with the combustion stand-off distance being controlled by the interaction of turbulent transport, droplet heating and vaporization, and gas-phase chemical reactions. In this study, conditions are identified under which analyses of laminar flamelets canshed light on aspects of turbulent spray ignition. This study extends earlier fundamental work by Liñan & Crespo (1976) on ignition in gaseous mixing layers to ignition of sprays. Studies of laminar mixing layers have been found to be instrumental in developing un-derstanding of turbulent combustion (Peters 2000), including the ignition of turbulent gaseous diffusion flames (Mastorakos 2009). For the spray problem at hand, the configuration selected, shown in Figure 1, involves a coflow mixing layer formed between a stream of hot air moving at velocity UA and a monodisperse spray moving at velocity USUA. The boundary-layer approximation will be used below to describe the resulting sl ender flow, which exhibits different igniting behaviors depending on the characteristics of t he fuel. In this approximation, consideration of the case U A = U S enables laminar ignition distances to be related to ignition times of unstrained spray flamelets, thereby pro viding quantitative information of direct applicability in regions of low scala r dissipation-rate in turbulent reactive flows (see the discussion in pp. 181–186 of Peters (2000)) . This report is organized as follows. Effects of droplet dispersion dynamics on ignition of sprays in turbulent mixing layers are discussed in Section 2. The formulation f or ignition in laminar mixing layers is outlined in Sections 3 and 4. The results are presented in Section 5. In Section 6, the mixture-fraction field and associated scalar dissipat ion rates for spray ignition are discussed. Finally, some brief conclusions are drawn in Section 7.