Effect of differential diffusion on turbulent lean premixed hydrogen enriched flames through structure analysis

This study presents the flame structure influenced by the differential diffusion effects and evaluates the structural modifications induced by the turbulence, thus to understand the coupling effects of the diffusively unstable flame fronts and the turbulence distortion. Lean premixed CH₄/H₂/air flames were conducted using a piloted Bunsen burner. Three hydrogen fractions of 0, 30% and 60% were adopted and the laminar flame speed was kept constant. The turbulence was generated with a single-layer perforated plate, which was combined with different bulk velocities to obtain varied turbulence intensities. Quasi-laminar flames without the plate were also performed. Explicit flame morphology was obtained using the OH-PLIF. The curvature, flame surface density and turbulent burning velocity were measured. Results show that the preferential transport of hydrogen produces negatively curved cusps flanked with positively curved bulges, which are featured by skewed curvature pdfs and consistent with the typical structure caused by the Darrieus-Landau instability. Prevalent bulge-cusp like wrinkles remain with relatively weak turbulence. However, stronger turbulence can break the bulges to be finer, and induce random positively curved cusps, therefore to destroy the bulge-cusp structures. Evident positive curvatures are generated in this process modifying the skewed curvature pdfs to be more symmetric, while the negative curvatures are not affected seriously. From low to high turbulence intensities, the hydrogen addition always strengthens the flame wrinkling. The augmentation of flame surface density and turbulent burning velocity with hydrogen is even more obvious at higher turbulence intensity. It is suggested that the differential diffusion can persist and even be strengthened with strong turbulence.     

Effect of hydrogen enrichment on flame broadening of turbulent premixed flames in thin reaction regime

An experimental study to identify the effect of hydrogen enrichment and differential diffusion on the flame broadening is conducted. Turbulent lean premixed flames in the Broadened Preheat–Thin Reaction (BP-TR) regime are obtained. The flames are stabilized on a Bunsen burner and CH₄/H₂/air mixtures are adopted with three hydrogen fractions of 0, 30% and 60%. The preheat zone and heat release zone are captured with the multi-species Planar Laser-Induced Fluorescence (PLIF) of OH and CH₂O radicals. Flame thicknesses of the preheat and heat release layers are measured. Results show broadened preheat zone and thin heat release layers for the flames, as predicted by the BP-TR regime. The preheat zone thickness can be increased to about 3–6 times compared to the laminar preheat thickness. An apparently decreased preheat zone thickness with hydrogen addition is observed. The differential diffusion is anticipated to locally thicken the heat release zone along the flame front. The mean heat release thickness is nearly not affected by the turbulence or hydrogen addition.     

Lewis number effect on the propagation of premixed flames in narrow adiabatic channels: Symmetric and non-symmetric flames and their linear stability analysis

The propagation of premixed flames with different Lewis numbers in a planar channel subject to a Poiseuille flow is considered within the diffusive-thermal model for steady and time-dependent cases. It was found that, depending on the Lewis number and the flow rate, symmetric and non-symmetric flames are possible. The existence of multiple steady solutions in cases of the low Lewis number is demonstrated. The time-dependent simulations carried out for high Lewis number flames also showed the symmetric and non-symmetric oscillatory solutions.Linear stability analysis of two-dimensional steady-states was performed using a practical method developed in the paper and applied to calculate the main eigenvalue. It was shown that for symmetric flames with a low Lewis number the increase in the flow rate leads to a loss of stability with subsequent formation of non-symmetric solutions. For flames with a high Lewis number the Poiseuille flow produces a stabilization effect. The results of the stability analysis were successfully compared with the results of direct numerical simulations.     

Heat release rate variations in high hydrogen content premixed syngas flames at elevated pressures: Effect of equivalence ratio

Three-dimensional direct numerical simulations with detailed chemistry were performed to investigate the effect of equivalence ratio on spatial variations of the heat release rate and flame markers of hydrogen/carbon monoxide syngas expanding spherical premixed flames under turbulent conditions at elevated pressures. The flame structures and the heat release rate were analysed and compared between fuel-lean, stoichiometric and fuel-rich centrally ignited spherical flames. The equivalence ratio changes the balance among thermo-diffusive effects, Darrieus–Landau instability and turbulence, leading to different flame dynamics and the heat release rate distribution, despite exhibiting similar cellular and wrinkling flames. The Darrieus–Landau instability is relatively insensitive to the equivalence ratio while the thermo-diffusive process is strongly affected by the equivalence ratio. As the thermo-diffusive effect increases as the equivalence ratio decreases, the fuel-lean flame is more unstable than the fuel-rich flame with the stoichiometric flame in between, under the joint effects of the thermo-diffusive instability and the Darrieus–Landau instability. The local heat release rate and curvature display a positive correlation for the lean flame, no correlation for the stoichiometric flame, and negative correlation for the rich flame. Furthermore, for the fuel-lean flame, the low and high heat release rate values are found in the negative and positive curvature zones, respectively, while for the fuel-rich flame, the opposite trends are found. It is found that heat release rate markers based on species concentrations vary strongly with changing equivalence ratio. The results suggest that the HCO, HO₂ concentrations and product of OH and CH₂O concentrations show good correlation with the local heat release rate for H₂/CO premixed syngas-air stoichiometric flame under turbulent conditions at elevated pressures.     

Experimental study of Markstein number effects on laminar flamelet velocity in turbulent premixed flames

Effects of turbulent flame stretch on mean local laminar burning velocity of flamelets, , were investigated experimentally in an explosion vessel at normal temperature and pressure. In this context, the wrinkling, At/Al, and the burning velocity, ut, of turbulent flames were measured simultaneously. With the flamelet assumption the mean local laminar burning velocity of flamelets, , was calculated for different turbulence intensities. The results were compared to the influence of stretch on spherically expanding laminar flames. For spherically expanding laminar flames the stretched laminar burning velocity, un, varied linearly with the Karlovitz stretch factor, yielding Markstein numbers that depend on the mixture composition. Six different mixtures with positive and negative Markstein numbers were investigated. The measurements of the mean local laminar burning velocity of turbulent flamelets were used to derive an efficiency parameter, I, which reflects the impact of the Markstein number and turbulent flame stretch—expressed by the turbulent Karlovitz stretch factor—on the local laminar burning velocity of flamelets. The results showed that the efficiency is reduced with increasing turbulence intensity and the reduction can be correlated to unsteady effects.     

A numerical study of the stability of one-dimensional laminar premixed flames in inert porous media

This work presents a numerical study of the stabilization diagram of methane/air premixed flames in a finite porous media foam with a uniform ambient temperature. A set of steady computations are considered, using a 1D numerical model that takes into account solid and gas energy equations as well as chemistry and radiation models. The present results show that both stable and unstable solutions, for upper and lower flames, exist either at the surface or submerged in the porous matrix. The influence of the 1D computational domain, boundary conditions, and gas/solid interface treatment on the stability of the calculated flames is also discussed. A linearized version of the discrete-ordinates radiation model is included in the linear stability analysis to discuss the influence of radiation on the stability of the flames. The full stabilization diagram and the linear stability analysis provide information on the stability of the flames, pointing to the existence of unstable upstream surface flames as well as unstable submerged flames on the downstream part of the porous media.     

Effect of high hydrogen enrichment on the outer-shear-layer flame of confined lean premixed CH4/H2/air swirl flames

In this study, we investigated the H₂-induced transition of confined swirl flames from the “V” to “M” shape. H₂-enriched lean premixed CH₄/H₂/air flames with H₂ fractions up to 80% were conducted. The flame structure was obtained with Planar Laser-Induced Fluorescence (PLIF) of the OH radical. Flow fields were measured with Particle Image Velocimetry (PIV). It was observed that the flame tip in the outer shear layer gradually propagated upstream and finally anchored to the injector with the hydrogen fractions increase, yielding the transition from the “V” to “M” flame. We examined the flame structures and the flame flow dynamics during the transition. The shape transition was directly related to the evolution of the corner flame along the outer shear layer. With H₂ addition, the outer recirculation zone first appeared downstream where the corner flame started to propagate upstream; then, the recirculation zone expanded upward to form a stable “M” flame gradually. The flow straining was observed to influence the stabilization of the outer shear layer flame significantly. This study can be useful for the understanding of recirculation-stabilized swirling flames with strong confinement. The flame structure and the flow characteristics of flames with a high H₂ content are also valuable for model validation.     

Effects of H2 and H preferential diffusion and unity Lewis number on superadiabatic flame temperatures in rich premixed methane flames

The structures of freely propagating rich CH₄/air and CH₄/O₂ flames were studied numerically using a relatively detailed reaction mechanism. Species diffusion was modeled using five different methods/assumptions to investigate the effects of species diffusion, in particular H₂ and H, on superadiabatic flame temperature. With the preferential diffusion of H₂ and H accounted for, significant amount of H₂ and H produced in the flame front diffuse from the reaction zone to the preheat zone. The preferential diffusion of H₂ from the reaction zone to the preheat zone has negligible effects on the phenomenon of superadiabatic flame temperature in both CH₄/air and CH₄/O₂ flames. It is therefore demonstrated that the superadiabatic flame temperature phenomenon in rich hydrocarbon flames is not due to the preferential diffusion of H₂ from the reaction zone to the preheat zone as recently suggested by Zamashchikov et al. [V.V. Zamashchikov, I.G. Namyatov, V.A. Bunev, V.S. Babkin, Combust. Explosion Shock Waves 40 (2004) 32]. The suppression of the preferential diffusion of H radicals from the reaction zone to the preheat zone drastically reduces the degree of superadiabaticity in rich CH₄/O₂ flames. The preferential diffusion of H radicals plays an important role in the occurrence of superadiabatic flame temperature. The assumption of unity Lewis number for all species leads to the suppression of H radical diffusion from the reaction zone to the preheat zone and significant diffusion of CO₂ from the postflame zone to the reaction zone. Consequently, the degree of superadiabaticity of flame temperature is also significantly reduced. Through reaction flux analyses and numerical experiments, the chemical nature of the superadiabatic flame temperature phenomenon in rich CH₄/air and CH₄/O₂ flames was identified to be the relative scarcity of H radical, which leads to overshoot of H₂O and CH₂CO in CH₄/air flames and overshoot of H₂O in CH₄/O₂ flames.     

Effects of Damköhler number on flame extinction and reignition in turbulent non-premixed flames using DNS

Results from a parametric study of flame extinction and reignition with varying Damköhler number using direct numerical simulation are presented. Three planar, non-premixed ethylene jet flames were simulated at a constant Reynolds number of 5120. The fuel and oxidizer stream compositions were varied to adjust the steady laminar extinction scalar dissipation rate, while maintaining constant flow and geometric conditions. Peak flame extinction varies from approximately 40% to nearly global blowout as the Damköhler number decreases. The degree of extinction significantly affects the development of the jets and the degree of mixing of fuel, oxidizer, and combustion products prior to reignition. The global characteristics of the flames are presented along with an analysis of the modes of reignition. It is found that the initially non-premixed flame undergoing nearly global extinction reignites through premixed flame propagation in a highly stratified mixture. A progress variable is defined and a budget of key terms in its transport equation is presented.     

Pressure influence on the flame front curvature of turbulent premixed flames: comparison between experiment and theory

In gas turbines, lean premixed combustion is executed in strongly turbulent flow fields and under high-pressure to allow large thermal loads within small-size combustors. Previous research on turbulent premixed flames has revealed the vital importance of flame-vortex interactions, but most of these investigations have been performed only at atmospheric pressure disregarding the large pressure dependency of the flame front dynamics. We report about spatially high-resolved laser-induced predissociation fluorescence imaging of OH (OH-LIPF) in premixed, high-pressure bluff-body stabilized methane/air flames. For each of the two measurement series with different equivalence ratio (φ = 0.7 and φ = 1.0), the planar flame topology at different pressures (0.1 to 1.1 MPa) but constant exit velocity was detected and stored for analysis. As the pressure was increased, the flame front contour of both equivalence ratios became strongly wrinkled with formation of highly curved flame front elements. For quantification of this phenomenon, the probability density function of flame curvature was evaluated with definition of the mean curvature radius as representative folding scale. To discuss different mechanisms of flame front disturbances according to their relevance, the flame curvature is compared with characteristic turbulence scales of the flow field and with the expected folding scale derived with Sivashinsky‘s formulation of linear flame instability theory. Significant changes become obvious especially if the pressure is increased up to 0.5 MPa. The mean curvature radius decreases distinctly and can be linked to the decreasing size of the Taylor length. Additionally, the formation of highly convoluted flame front elements is enforced by the increasing flame instability behavior. As the results show, the flame stoichiometry has a strong impact on the flame front topology at increasing pressures due to the differences of their flame dynamics.     

Hysteresis loops of methane catalytic partial oxidation for hydrogen production under the effects of varied Reynolds number and Damköhler number

Hysteresis loops of catalytic partial oxidation of methane (CPOM) for hydrogen production under the effects of varied Reynolds number and Damköhler number are investigated numerically in this study. The physical phenomena are predicted using the indirect mechanism, which consists of the total oxidation (or combustion), steam reforming and CO₂ reforming of methane in a catalyst bed. Numerical results reveal that, when the Damköhler number is relatively low, a hysteresis loop of CPOM from varying Reynolds number develops. Increasing the Damköhler number leads to the loop shifting toward the regime of high Reynolds number. However, once the Damköhler number is large to a certain extent, the chemical reactions are always exhibited for the Reynolds number less than 2000. A closed loop is thus not observed. Alternatively, for a given Reynolds number, an ignited Damköhler number and an extinguished Damköhler number can be obtained. Accordingly, three different regions in the plot of Damköhler number versus Reynolds number are identified. Physically, when the role played by Damköhler number on CPOM is much more important than by the Reynolds number (Region I), the thermal effect governs the chemical reactions. In contrast, if the Reynolds number plays a key role in determining the CPOM (Region III), the chemically frozen flow prevails over the catalyst bed. When the residence times of the total oxidation and convection in the catalyst bed are in an equivalent state (Region II), CPOM is characterized by a dual-solution, rendering the hysteresis loops. From the distributions of ignited and extinguished Damköhler numbers, the catalytic reactor and operation of partial oxidation of methane and other fuels can be designed accordingly.     

A numerical study on the effects of pressure and gravity in laminar ethylene diffusion flames

The effects of pressure and gravity on sooting characteristics and flame structure were studied numerically in coflow ethylene–air laminar diffusion flames between 0.5 and 5 atm. Computations were performed by solving the unmodified and fully-coupled equations governing reactive, compressible, gaseous mixtures which include complex chemistry, detailed radiation heat transfer, and soot formation/oxidation. Soot formation/oxidation was modeled using an acetylene-based, semi-empirical model which has been verified with previously published experimental data to correctly capture many of the observed trends at normal-gravity. Calculations for each pressure considered were performed for both normal- and zero-gravity conditions to help separate the effects of pressure and buoyancy on soot formation. Based on the numerical predictions, pressure and gravity were observed to significantly influence the flames through their effects on buoyancy and reaction rates. The zero-gravity flames have higher soot concentrations, lower temperatures and broader soot-containing zones than normal-gravity flames at the same pressure. The zero-gravity flames were also found to be longer and wider. Differences were observed between the two levels of gravity when pressure was increased. The zero-gravity flames displayed a stronger dependence of the maximum soot yield on pressure from 0.5 to 2 atm and a weaker dependence from 2 to 5 atm as compared to the normal-gravity flames. In addition, flame diameter decreased with increasing pressure under normal-gravity while it increased with pressure in the zero-gravity cases. Changing the prescribed wall boundary condition from fixed-temperature to adiabatic significantly altered the numerical predictions at 5 atm. When the walls were assumed to be adiabatic, peak soot volume fractions and temperatures increased in both the zero- and normal-gravity flames, emphasizing the importance of heat conduction to the burner rim on flame structure.     

Effect of elliptic burner geometry and air equivalence ratio on the nitric oxide emissions from turbulent hydrogen flames

An experimental study on turbulent hydrogen flames from circular and elliptic burners with varying degrees of premixedness (diffusion, fuel-rich, stoichiometric, and fuel-lean) is presented. Flame stability, visible flame height, flame radiation, global nitric oxide (NO) concentration, and inflame temperature and NO concentration profiles were measured. We found that the elliptic burner flames had lower liftoff velocity, were shorter, and radiated less heat to the surrounding as compared to circular burner flames. Global NO concentration decreased with an increase in air equivalence ratio for both circular and elliptic burner flames. Peak in-flame NO concentration along the flame centerline increased with a decrease in air equivalence ratio. Elliptic burner flames produced higher peak in-flame temperatures. Overall, the elliptic burner flames produced less peak NO as compared to circular burner flames at all air equivalence ratios except zero (diffusion flames) in accordance with the global emission measurements.     

The interaction between internal heat loss and external heat loss on the extinction of stretched spray flames with nonunity Lewis number

The influences of flow stretch, preferential diffusion, internal heat transfer and external heat loss on the extinction of dilute spray flames propagating in a stagnation-point flow are analyzed using activation energy asymptotics. A completely prevaporized mode and a partially prevaporized mode of flame propagation are identified. The internal heat transfer, associated with the liquid fuel loading and the initial droplet size of the spray, provides heat loss for rich sprays but heat gain for lean sprays. The flow stretch respectively weakens and intensifies the burning intensity of the lean methanol-spray flame (Le>1) and rich methanol-spray flame (Le<1). Results show that the Le>1 flame can be extinguished with or without external heat loss. Flame extinction characterized by a C-shaped curve is dominated by the external heat loss or the flow stretch. For the Le<1 flame without external heat loss, no extinction occurs under the influence of flow stretch. However, the Le<1 flame with completely prevaporized fuel sprays enduring a small amount of flow stretch can be extinguished by the external heat loss and this behavior is characterized by a C-shaped curve. Note that the W-shaped extinction curve is mainly governed by the internal heat loss. Flame extinction characterized by a W-shaped curve occurs when the Le<1 spray flame with external heat loss endures a positive stretch and experiences a partially prevaporized spray with sufficiently large liquid fuel loading and droplet size.     

Effects of strain rate,turbulence, reactant stoichiometry and heat losses on the interaction of turbulent premixed flames with stoichiometric counterflowing combustion products

Intense strain, turbulence, heat transfer, and mixing with combustion products can affect premixed flames in practical combustion devices. These effects are systematically studied in turbulent premixed CH₄/N₂/O₂ flames using a reactant versus product counterflow system and independently varying bulk strain rate, turbulent Reynolds number, equivalence ratio of the reactant mixture, and temperature of the stoichiometric counterflowing combustion products. The flow field and the turbulent flames are investigated using particle image velocimetry (PIV) measurements and laser-induced fluorescence (LIF) imaging of OH. The OH-LIF images are used to identify the interface between the counterflowing streams, referred to here as the gas mixing layer interface (GMLI). The flame response for different flow conditions is compared in terms of the probability of localized extinction along the GMLI, the turbulent flame brush thickness, and flame position relative to the GMLI, by using an OH-LIF-based progress variable. The probability of localized extinction at the GMLI increases as the separation between the turbulent flame brush and the GMLI decreases. Flame fronts in the vicinity of the GMLI are more likely to extinguish as a result of heat losses, dilution of the reaction zone by the product stream, and large local strain rates. A higher probability of localized extinction at the GMLI is induced by either a larger bulk strain rate or a slower flame speed. As the turbulent Reynolds number increases, the corresponding increase in turbulent flame brush thickness enhances the interactions of the flame fronts with the GMLI. Heat losses are substantially less significant for cases in which the turbulent flame brush is sufficiently separated from the GMLI. For flames in close proximity to the GMLI, the effects of the product stream on the flame front differ for lean and rich reactant mixtures. These disparities are attributed in part to differences in the ignitibility of the reactant mixtures by the hot product stream.     

Effects of H2 or CO2 addition, equivalence ratio, and turbulent straining on turbulent burning velocities for lean premixed methane combustion

Using hydrogen or carbon dioxide as an additive, we investigate the bending effect of turbulent burning velocities (ST/SL) over a wide range of turbulent intensities (u^′/SL) up to 40 for lean premixed methane combustion at various equivalence ratios (?), where SL is the laminar burning velocity. Experiments are carried out in a cruciform burner, in which a sizable downward-propagating premixed CH₄/diluent/air flame interacts with intense isotropic turbulence in the central region without influences of ignition and unwanted turbulence from walls. Simultaneous measurements using the pressure transducer and pairs of ion-probe sensors at various positions of the burner show that effects of gas velocities and pressure rise due to turbulent combustion on ST of lean CH₄/H₂/air flames can be neglected, confirming the accuracy of the ST data. Results with increasing hydrogen additions (δ=10, 20, and 30% in volume) show that the bending of ST/SL vs u^′/SL plots is diminished when compared to data with δ=0, revealing that high reactivity and diffusivity of hydrogen additives help the reaction zone remaining thin even at high u^′/SL. In contrast, the bending effect is strongly promoted when CO₂ is added due to radiation heat losses. This leads to lower values of ST/SL at fixed u^′/SL and ?, where the slope n can change signs from positive to negative at sufficiently large u^′/SL, suggesting that the reaction zone is no longer thin. All ST data with various δ can be well approximated by a general correlation (ST−SL)/u^′=0.17Da^0.43, covering both corrugated flamelet and distributed regimes with very small data scatter, where Da is the turbulent Damköhler number. These results are useful in better understanding how turbulence and diluents can influence the canonical structures of turbulent premixed flames and thus turbulent burning rates.     

A 5 kWt catalytic burner for PEM fuel cells: Effect of fuel type,fuel content and fuel loads on the capacity of the catalytic burner

For proton exchange membrane fuel cell systems (PEMFC) integrated with fuel processors, the calorific value of reformate gases produced during the start-up phase must be recovered. An appropriate exhaust after treatment system has crucial importance for PEMFC systems. Catalytic combustion is a promising alternative regarding its total oxidation capability of low calorific value gases at low temperatures, thereby reducing environmentally hazardous emissions. The aim of the study is to develop an after treatment system using a catalytic burner with a nominal capacity of 5 kW_t, which is also adaptive to partial loads of PEM fuel cell capacity. Fuel type, fuel composition and fuel loads are important parameters determining the operating window of the catalytic burner. Precious metal based catalysts, as proved to be the most active catalysts for the oxidation of hydrocarbons, can withstand temperatures of about 1073 K without exhibiting a rapid deactivation. This is the main barrier dictating the operating window and thereby determining the capacity of the burner. In this work, 1.5% natural gas (NG) alone was found to be the upper limit to control the catalyst bed temperature below 1073 K. In the case of catalytic combustion of hydrogen–NG mixture, 7% of hydrogen with NG up to 0.6% could be totally oxidized below 1073 K. Within the experimented ranges of fuel loads, between 2.5 kW_t and 5.5 kW_t, the temperature of the catalyst bed was seen to increase with increasing the fuel load at constant fuel percentages. It has been observed that fuel type was another parameter affecting the exhaust gas temperature.