Numerical study of methanol flames in laminar forced convective environment using short chemical kinetics mechanism

Due to recent interest in methanol economy, it is seen that a numerical study of combustion of methanol in a comprehensive manner is necessary. Motivated from this interest and based on the studies from literature, a numerical study on prediction of structures of non-premixed methanol-air flames in laminar forced convective environment is reported. Two-dimensional, planar and axisymmetric, computational domains are considered. Corresponding governing equations for conservation of mass, momentum, species and energy have been solved using Ansys FLUENT. The numerical model incorporates multi-component diffusion, variable thermal and physical properties, a short chemical kinetics mechanism with 18 species and 38 elementary reactions, and a non-luminous thermal radiation model. Homogeneous flames in opposed flow and heterogeneous flames in cross-flow and co-flow configurations are studied. For heterogeneous flames, interface conditions at the liquid methanol surface are defined systematically using a user-defined function. Numerical results are validated against the experimental results available in literature. Results in terms of mass burning rates, flow, species and temperature fields have been presented to describe the flame characteristics.     

Turbulence/flame/wall interactions in non-premixed inclined slot-jet flames impinging at a wall using direct numerical simulation

In the present work, three-dimensional turbulent non-premixed oblique slot-jet flames impinging at a wall were investigated using direct numerical simulation (DNS). Two cases are considered with the Damköhler number (Da) of case A being twice that of case B. A 17 species and 73-step mechanism for methane combustion was employed in the simulations. It was found that flame extinction in case B is more prominent compared to case A. Reignition in the lower branch of combustion for case A occurs when the scalar dissipation rate relaxes, while no reignition occurs in the lower branch for case B due to excessive scalar dissipation rate. A method was proposed to identify the flame quenching edges of turbulent non-premixed flames in wall-bounded flows based on the intersections of mixture fraction and OH mass fraction iso-surfaces. The flame/wall interactions were examined in terms of the quenching distance and the wall heat flux along the quenching edges. There is essentially no flame/wall interaction in case B due to the extinction caused by excessive turbulent mixing. In contrast, significant interactions between flames and the wall are observed in case A. The quenching distance is found to be negatively correlated with wall heat flux as previously reported in turbulent premixed flames. The influence of chemical reactions and wall on flow topologies was identified. The FS/U and FC/U topologies are found near flame edges, and the NNN/U topology appears when reignition occurs. The vortex-dominant topologies, FC/U and FS/S, play an increasingly important role as the jet turbulence develops.     

Simulations of laminar non-premixed flames of kerosene with hot combustion products as oxidiser

The effects of hot combustion product dilution in a pressurised kerosene-burning system at gas turbine conditions were investigated with laminar counterflow flame simulations. Hot combustion products from a lean (φ = 0.6) premixed flame were used as an oxidiser with kerosene surrogate as fuel in a non-premixed counterflow flame at 5, 7, 9 and 11 bar. Kerosene-hot product flames, referred to as ‘MILD’, exhibit a flame structure similar to that of kerosene–air flames, referred to as ‘conventional’, at low strain rates. The Heat Release Rate (HRR) of both conventional and MILD flames reflects the pyrolysis of the primary and intermediate fuels on the rich side of the reaction zone. Positive HRR and OH regions in mixture fraction space are of similar width to conventional kerosene flames, suggesting that MILD flames are thin fronts. MILD flames do not exhibit typical extinction behaviour, but gradually transition to a mixing solution at very high rates of strain (above A = 160, 000 s^?1 for all pressures). This is in agreement with literature that suggests heavily preheated and diluted flames have a monotonic S-shaped curve. Despite these differences in comparison with kerosene–air flames, MILD flames follow typical trends as a function of both strain and pressure. Further still, the peak locations of the overlap of OH and CH₂O mass fractions in comparison with the peak HRR indicate that the pixel-by-pixel product of OH- and CH₂O-PLIF signals is a valid experimental marker for non-premixed kerosene MILD and conventional flames.     

A skeletal mechanism for flame inhibition by trimethylphosphate

On the basis of a multi-step kinetic mechanism for flame inhibition by organophosphorus compounds including more than 200 reactions, a skeletal mechanism for flame inhibition by trimethylphosphate was developed. The mechanism consists of 22 irreversible elementary reactions, involving nine phosphorus-containing species. Selection of the crucial steps was performed by analysing P-element fluxes from species to species and by calculating net reaction rates of phosphorus-involving reactions versus the flames zone. The developed mechanism was validated by comparing the modelling results with the measured and simulated (using the starting initial mechanism) speed and the chemical structure of H₂/O₂, CH₄/O₂ and syngas/air flames doped with trimethylphosphate. The mechanism was shown to satisfactorily predict the speed of H₂/O₂/N₂ flames with various dilution ratios, CH₄/air and syngas/air flames doped with trimethylphosphate. The skeletal mechanism was also shown to satisfactorily predict the spatial variation of H and OH radicals and the final phosphorus-containing products of the inhibitor combustion. Further reduction of the skeletal mechanism without modification of the rate constants recommended in the starting mechanism was shown to result in noticeable disagreement of the flame speed and structure.     

Novel insights into mild combustion processes through analyses of hysteresis behavior

In many numerical scientific papers, MILD combustion is defined on the basis of the disappearance of extinction phenomena while varying mixture dilution levels, through the identification of the condition where ignition and extinction collapse in an unique event. Albeit with numerous contributions in elementary reactor configurations (Continuous Stirred Flow Reactors, Opposed-Diffusion flames), operating conditions to achieve this noticeable circumstance do not properly match with experimental evidences, since they lead to extreme mixture dilution levels (higher than 99% for CH₄). In addition, many works suggest the occurrence of extinctions phenomena or instabilities under MILD conditions, symptomatic of extinction/re-ignition phenomena. Simulations in an adiabatic CSTR for CH₄/O₂/N₂ mixtures, varying the N₂ content, were realized with the focus to analyze the system behavior through hysteresis. Results show different hysteresis behaviors from “air” to diluted conditions, demonstrating the occurrence of the condition T_ext=T_ign is not a strict constraint for MILD combustion systems. In addition, the “unstable” branch has characteristic temperatures lower than the ones related to both the ignition and extinction events, enlarging the opportunity to stabilize MILD combustion processes. These results pave the way to the reconsideration of peculiar aspects of MILD combustion processes.     

Rate-ratio asymptotic analysis of non-premixed methane flames

The asymptotic structure of laminar, non-premixed methane flames is analysed using a reduced four-step chemical-kinetic mechanism. Chemical reactions are presumed to take place in two layers: the inner layer and the oxidation layer. In the inner layer the fuel reacts with radicals and the main compounds formed are the intermediate species CO and H₂. These intermediate species are oxidized in the oxidation layer. The structure of the oxidation layer is described by two second-order differential equations: one for CO and the other for H₂. Two limiting cases are considered. At one limit the global step CO+H₂?CO₂+H₂ is presumed to maintain partial equilibrium everywhere in the oxidation layer except in a thin layer adjacent to the inner layer. At the other limit the steady-state approximation is introduced for H₂ everywhere in the oxidation layer except in a thin layer adjacent to the inner layer. This limit, called ‘slow CO oxidation’, has not been analysed previously. The structure of the inner layer is described by two second-order differential equations: one for the fuel and the other for the H radicals. This is a significant improvement over previous models in which either a steady-state approximation is introduced for the H radicals in the inner layer, or the reaction between the fuel and radicals is presumed to be very fast. The chain-breaking elementary reaction CH₃+H+M→CH₄+M is found to have a significant influence on the structure of the inner layer and on the scalar dissipation rates at extinction. The influence of this reaction was either neglected in previous models or was included as a perturbation to the principal elementary reactions taking place to the leading order in the inner layer. Using the results of the asymptotic analysis the scalar dissipation rates at extinction are calculated at a pressure of 1 bar. They are found to agree well with those calculated numerically using a chemical-kinetic mechanism made up of elementary reactions.     

The role of cool-flame chemistry in quasi-steady combustion and extinction of n-heptane droplets

Experiments on the combustion of large n-heptane droplets, performed by the National Aeronautics and Space Administration in the International Space Station, revealed a second stage of continued quasi-steady burning, supported by low-temperature chemistry, that follows radiative extinction of the first stage of burning, which is supported by normal hot-flame chemistry. The second stage of combustion experienced diffusive extinction, after which a large vapour cloud was observed to form around the droplet. In the present work, a 770-step reduced chemical-kinetic mechanism and a new 62-step skeletal chemical-kinetic mechanism, developed as an extension of an earlier 56-step mechanism, are employed to calculate the droplet burning rates, flame structures, and extinction diameters for this cool-flame regime. The calculations are performed for quasi-steady burning with the mixture fraction as the independent variable, which is then related to the physical variables of droplet combustion. The predictions with the new mechanism, which agree well with measured autoignition times, reveal that, in decreasing order of abundance, H₂O, CO, H₂O₂, CH₂O, and C₂H₄ are the principal reaction products during the low-temperature stage and that, during this stage, there is substantial leakage of n-heptane and O₂ through the flame, and very little production of CO₂ with no soot in the mechanism. The fuel leakage has been suggested to be the source of the observed vapour cloud that forms after flame extinction. While the new skeletal chemical-kinetic mechanism facilitates understanding of the chemical kinetics and predicts ignition times well, its predicted droplet diameters at extinction are appreciably larger than observed experimentally, but predictions with the 770-step reduced chemical-kinetic mechanism are in reasonably good agreement with experiment. The computations show how the key ketohydroperoxide compounds control the diffusion-flame structure and its extinction.     

The two-dimensional structure of low strain rate counterflow nonpremixed-methane flames in normal and microgravity

The structure and extinction of low strain rate nonpremixed methane–air flames was studied numerically and experimentally. A time-dependent axisymmetric two-dimensional (2D) model considering buoyancy effects and radiative heat transfer was developed to capture the structure and extinction limits of normal gravity (1-g) and zero gravity (0-g) flames. For comparison with the 2D modelling results, a one-dimensional (1D) flamelet computation using a previously developed numerical code was exercised to provide information on the 0-g flames. A 3-step global reaction mechanism was used in both the 1D and 2D computations to predict the measured extinction limit and flame temperature. Photographic images of flames undergoing the process of extinction were compared with model calculations. The axisymmetric numerical model was validated by comparing flame shapes, temperature profiles, and extinction limits with experiments and with the 1D computational results. The 2D computations yielded insight into the extinction mode and flame structure. A specific maximum heat release rate was introduced to quantify the local flame strength and to elucidate the extinction mechanism. The contribution by each term in the energy equation to the heat release rate was evaluated to investigate the multi-dimensional structure and radiative extinction of the 1-g flames. Two combustion regimes depending on the extinction mode were identified. Lateral heat loss effects and multi-dimensional flame and flow structure were also found. At low strain rates in 1-g flames (‘regime A’), the flame is extinguished from the weak outer edge of the flame, which is attributed to a multi-dimensional flame structure and flow field. At high strain rates, (‘regime B’), the flame extinction initiates near the flame centreline owing to an increased diluent concentration in the reaction zone, similar to the extinction mode of 1D flames. These two extinction modes can be clearly explained by consideration of the specific maximum heat release rate.     

Flame studies of C2 hydrocarbons

The oxidation characteristics of C₂ hydrocarbons were revisited in flames established in the counterflow configuration. Laminar flame speeds of ethane/air, ethylene/air, and acetylene/oxygen/nitrogen mixtures as well as extinction strain rates of non-premixed ethane/air flames were measured using digital particle image velocimetry. The experiments were modeled using three different kinetic models. While the experimental and computed laminar flame speeds agreed closely for all C₂ hydrocarbons under fuel-lean conditions, notable discrepancies were identified under fuel-rich conditions. Using the computed flame structures, insight was provided into the controlling mechanisms that could be responsible for the observed discrepancies. More specifically, the uncertainties associated with the kinetics of the thermal decomposition of the ethyl radical were found to be a potential source of the observed discrepancies for ethane flames. It was shown also by using alternative rate constants for the ethyl radical decomposition, the rate of flame propagation and the extinction propensity are affected notably. Furthermore, the values of the branching ratio of acetylene consumption reactions involving atomic oxygen were found to have a significant effect on the propagation of rich acetylene flames.     

Transported PDF simulation of turbulent CH4/H2 flames under MILD conditions with particle-level sensitivity analysis

Transported probability density function (TPDF) simulation with sensitivity analysis has been conducted for turbulent non-premixed CH₄/H₂ flames of the jet-into-hot-coflow (JHC) burner, which is a typical model to emulate moderate or intense low oxygen dilution combustion (MILD). Specifically, two cases with different levels of oxygen in the coflow stream, namely HM1 and HM3, are simulated to reveal the differences between MILD and hot-temperature combustion. The TPDF simulation well predicts the temperature and species distributions including those of OH, CO and NO for both cases with a 25-species mechanism. The reduced reaction activity in HM1 as reflected in the peak OH concentration is well correlated to the reduced oxygen in the coflow stream. The particle-level local sensitivities with respect to mixing and chemical reaction further show dramatic differences in the flame characteristics. HM1 is less sensitive to mixing and reaction parameters than HM3 due to the suppressed combustion process. Specifically, for HM1 the sensitivities to mixing and chemical reactions have comparable magnitude, indicating that the combustion progress is controlled by both mixing and reaction in MILD combustion. For HM3, there is however a change in the combustion mode: during the flame initialization, the combustion progress is more sensitive to chemical reactions, indicating that finite-rate chemistry is the controlling process during the autoignition process for flame stabilization; at further downstream where the flame has established, the combustion progress is controlled by mixing, which is characteristic of nonpremixed flames. An examination of the particles with the largest sensitivities reveals the difference in the controlling mixtures for flame stabilization, namely, the stoichiometric mixtures are important for HM1, whereas, fuel-lean mixtures are controlling for HM3. The study demonstrates the potential of TPDF simulations with sensitivity analysis to investigate the effects of finite-rate chemistry on the flame characteristics and emissions, and reveal the controlling physio-chemical processes in MILD combustion.     

Thermal and chemical effects of differential diffusion in turbulent non-premixed H2 flames

As a renewable fuel, hydrogen (H₂) may play an increasingly important role in the development and control of piston and gas turbine engines to achieve zero carbon emissions. Predictive modeling of H₂-fueled combustion processes requires a clear understanding of differential diffusion (DD) due to the high diffusivity of H₂. On the assumption that turbulent mixing is a far more dominant process than molecular mixing, DD effects are typically neglected in turbulent combustion simulations to reduce modeling complications. While this assumption is reasonable for hydrocarbon fuels, it is less valid for H₂ combustion, where DD is significant. In this work, two three-dimensional direct numerical simulations of temporally evolving turbulent H₂ jet flames with and without considering DD are performed and compared with laminar flamelet solutions to assess DD effects under turbulent conditions. The emphasis is placed on assessing the suitability of classical mixture fraction Z and Bilger mixture fraction Z_Bilger as conditioning variables for non-premixed turbulent combustion modeling through analyzing DD effects on flame structure, chemical reactions, and tangential diffusion (TD). Furthermore, the persistence of DD effects under turbulent conditions and the suitability of a conventional DD parameter are investigated by comparing the turbulent flames to laminar flamelet solutions. It is found that conditioning the thermochemical state on Z_Bilger helps to capture DD effects and mitigate the relative contribution of TD, which gives Z_Bilger advantages over Z when employing flamelet modeling. Due to close coupling between DD and local chemical reactions, DD can affect the turbulent/laminar flames in the form of thermal effects due to the change in flame temperature, chemical effects due to the change in chemical reactions, and transport effects due to multiple species with varying diffusivities that could result in the difference between Z and Z_Bilger. While the transport effects are suppressed, significant chemical and thermal effects of DD still persist under turbulent conditions, which indicates that the DD parameter is probably unsuitable for comprehensively characterizing and assessing DD effects on the structure of turbulent non-premixed flames.     

Global reduced mechanisms for methane and hydrogen combustion with nitric oxide formation constructed with CSP data

Reduced mechanisms for methane-air and hydrogen-air combustion including NO formation have been constructed with the computational singular perturbation (CSP) method using the fully automated algorithm described by Massias et al. The analysis was performed on solutions of unstrained adiabatic premixed flames with detailed chemical kinetics described by GRI 2.11 for methane and a 71-reaction mechanism for hydrogen including NO _x formation. A 10-step reduced mechanism for methane has been constructed which reproduces accurately laminar burning velocities, flame temperatures and mass fraction distributions of major species for the whole flammability range. Many steady-state species are also predicted satisfactorily. This mechanism is an improvement over the seven-step set of Massias et al, especially for rich flames, because the use of HCNO, HCN and C₂H₂ as major species results in a better calculation of prompt NO. The present 10-step mechanism may thus also be applicable to diffusion flames. A five-step mechanism for lean and hydrogen-rich combustion has also been constructed based on a detailed mechanism including thermal NO. This mechanism is accurate for a wide range of the equivalence ratio and for pressures as high as 40 bar. For both fuels, the CSP algorithm automatically pointed to the same steady-state species as those identified by laborious analysis or intuition in the literature and the global reactions were similar to well established previous methane-reduced mechanisms. This implies that the method is very well suited for the study of complex mechanisms for heavy hydrocarbon combustion.M This article features supplementary data files available from the supplemental page in the online journal.     

LIF measurements and chemical kinetic analysis of methylidyne formation in high-pressure counter-flow partially premixed and non-premixed flames

We report quantitative, spatially resolved, linear laser-induced fluorescence (LIF) measurements of methylidyne concentration ([CH]) in laminar, methane–air, counter-flow partially premixed and non-premixed flames using excitation near 431.5 nm in the A–X (0,0) band. For partially premixed flames, fuel-side equivalence ratios (_B) of 1.45, 1.6 and 2.0 are studied at pressures of 1, 3, 6, 9 and 12 atm. For non-premixed flames, the fuel-side mixture consists of 25% CH₄ and 75% N₂; measurements are obtained at pressures of 1, 2, 3, 4, 5, 6, 9 and 12 atm. The quantitative CH measurements are compared with predictions from an opposed-flow flame code utilizing two GRI chemical kinetic mechanisms (versions 2.11 and 3.0). LIF measurements of [CH] are corrected for variations in the quenching rate coefficient by using major species concentrations and temperatures generated by the code along with suitable quenching cross sections for CH available from the literature. A pathway analysis provides relative contributions from important elementary reactions to the total amount of CH produced at various pressures. Key reactions controlling peak CH concentrations are also identified by using a sensitivity analysis. For the partially premixed flames, measured CH profiles are reproduced reasonably well by GRI 3.0, although some quantitative disagreement exists at all pressures. Two CH radical peaks are observed for _B=1.45 and _B=1.6 at pressures above 3 atm. Peak CH concentrations for the non-premixed flames are significantly underpredicted by GRI 3.0. The latter agrees with previously reported NO concentrations, which are also underpredicted in these same high-pressure counter-flow diffusion flames. PACS 07.35.+k; 42.62.Fi; 82.33.Vx     

Analysis of kinetic mechanism performance in conditional moment closure modelling of turbulent,non-premixed methane flames

This paper presents results obtained from the application of a first-order conditional moment closure approach to the modelling of two methane flames of differing geometries. Predictions are based upon a second-moment turbulence and scalar-flux closure, and supplemented with full and reduced chemical kinetic mechanisms, ranging from a simple 12-step to a complex 1207-step mechanism. Alongside analysis of the full kinetic schemes' performance, is an appraisal of the behaviour of their derivatives obtained using mechanism-reduction techniques. The study was undertaken to analyse the practicality of incorporating kinetic models of varying complexity into calculations of turbulent non-premixed flames, and to make comparison of their performance. Despite extensive studies of the predictive ability of such schemes under laminar flame conditions, systematic evaluations have not been performed for turbulent reacting flows. This paper reflects upon the impact that selection of chemical kinetics has upon subsequent calculations and concludes that, although application of reduced schemes is more than adequate to reproduce experimental data, selection of the parent mechanism is of paramount importance to the prediction of minor species. Although widely used schemes are well documented and validated, their performances vary considerably. Thus, careful consideration must be made to their application and origins during the evaluation of combustion models.     

On transient behavior of non-premixed counter-flow diffusion flames within the REDIM based model reduction concept

In the present work non-stationary behavior of the counter-flow diffusion flame is examined in the context of the recently developed approach of model reduction called REaction–DIffusion Manifolds (REDIM) method. It is a natural extension of the ILDM approach which takes into account both the chemical reaction and the diffusion processes. It has been developed to treat both premixed and non-premixed regimes of combustion. In this work we investigate the ability of the concept to describe transient processes of extinction and re-ignition. A very simple flame configuration and transport model are considered in this current study for the sake of transparency because the main focus is on the transient and non-stationary behavior of flames. H₂/O₂/N₂ combustion system is considered in a non-premixed counter-flow diffusion 1D flame configuration. This study shows how the REDIM concept performs in the transient regimes; it interprets the effect of local extinction and reigniting phenomena using detailed and reduced models. It shows how the unstable/transient behavior of a detailed system can be successfully accounted with the help of the REDIM based reduced model.     

The structure and extinction of nonpremixed methane/nitrous oxide and ethane/nitrous oxide flames

Knowledge of combustion of hydrocarbon fuels with nitrogen-containing oxidizers is a first step in understanding key aspects of combustion of hypergolic and gun propellants. Here an experimental and kinetic-modeling study is carried out to elucidate aspects of nonpremixed combustion of methane (CH₄) and nitrous oxide (N₂O), and ethane (C₂H₆) and N₂O. Experiments are conducted, at a pressure of 1 atm, on flames stabilized between two opposing streams. One stream is a mixture of oxygen (O₂), nitrogen (N₂), and N₂O, and the other a mixture of CH₄ and N₂ or C₂H₆ and N₂. Critical conditions for extinction are measured. Kinetic-modeling studies are performed with the San Diego Mechanism. Experimental data and results of kinetic-modeling show that N₂O inhibits the flame by promoting extinction. Analysis of the flame structure shows that H radicals are produced in the overall chain-branching step 3H₂ + O₂ ? 2H₂O + 2H, in which molecular hydrogen is consumed. Hydrogen is also consumed in the overall step N₂O + H₂ ? N₂ + H₂O where stable products are formed. Inhibition of the flames by N₂O is attributed to competition between these two overall steps.     

Paradigms in turbulent combustion research

The development of the basic conceptual viewpoints, or paradigms, for turbulent combustion in gases over the last 50 years is reviewed. Significant progress has been made. Recent successes in the prediction of pollutant species and extinction/re-ignition phenomena in non-premixed flames are seen as the result of close interaction between experimentalists, theoreticians, and modellers. Premixed turbulent flames seem to be dependent on a much wider range of factors, and predictive capabilities are not so advanced. Implications for large eddy simulation (LES) and partially premixed combustion are outlined.     

Three-dimensional numerical simulations of cellular jet diffusion flames

Recent experimental investigations have demonstrated that the appearance of particular cellular states in circular non-premixed jet flames significantly depends on a number of parameters, including the initial mixture strength, reactant Lewis numbers, and proximity to the extinction limit (Damköhler number). For CO₂-diluted H₂/O₂ jet diffusion flames, these studies have shown that a variety of different cellular patterns or states can form. For given fuel and oxidizer compositions, several preferred states were found to co-exist, and the particular state realized was determined by the initial conditions. To elucidate the dynamics of cellular instabilities, circular non-premixed jet flames are modeled with a combination of three-dimensional numerical simulation and linear stability analysis (LSA). In both formulations, chemistry is described by a single-step, finite-rate reaction, and different reactant Lewis numbers and molecular weights are specified. The three-dimensional numerical simulations show that different cellular flames can be obtained close to extinction and that different states co-exist for the same parameter values. Similar to the experiments, the behavior of the cell structures is sensitive to (numerical) noise. During the transient blow-off process, the flame undergoes transitions to structures with different number of cells, while the flame edge close to the nozzle oscillates in the streamwise direction. For conditions similar to the experiments discussed, the LSA results reveal various cellular instabilities, typically with azimuthal wavenumber m = 1–6. Consistent with previous theoretical work, the propensity for the cellular instabilities is shown to increase with decreasing reactant Lewis number and Damköhler number.     

The influence of water on extinction and ignition of hydrogen and methane flames

The influence of water vapor on critical conditions of extinction and autoignition of premixed and nonpremixed flames is investigated. The fuels tested are hydrogen (H₂) and methane (CH₄). Studies on premixed systems are carried out by injecting a premixed reactant stream made up of fuel, oxygen (O₂), and nitrogen (N₂) from one duct, and an inert-gas stream of N₂ from the other duct. Critical conditions of extinction are measured for various amounts of water vapor added to the premixed reactant stream. The ratio of fuel to oxygen is maintained at a constant value, and the amounts of water vapor and nitrogen are so chosen that the adiabatic temperature remains the same. This ensures that the physical influence of water is the same for all cases. Therefore, changes in values for the critical conditions of extinction are attributed to the chemical influence of water vapor. Studies on nonpremixed systems are carried out by injecting a fuel stream made up of fuel and N₂ from one duct ,and an oxidizer stream made up of O₂ and N₂ from the other duct. Critical conditions of extinction are measured with water vapor added to the oxidizer stream. The concentrations of reactants are so chosen that the adiabatic temperature and the flame position stay the same for all cases. Critical conditions of autoignition are measured by preheating the oxidizer stream of the nonpremixed system. Water vapor is added to the oxidizer stream. Numerical calculations are performed using a detailed chemical-kinetic mechanism and compared with measurements. Experimental and numerical studies show that addition of water makes the premixed and nonpremixed flames easier to extinguish and harder to ignite. The chemical influence of water is attributed to its enhanced chaperon efficiency in three body reactions.