Explosive dynamics of bluff-body-stabilized lean premixed hydrogen flames at blow-off

Two-dimensional direct numerical simulation (DNS) databases of bluff-body-stabilized lean hydrogen flames representative of complicated reactive–diffusive system are analysed using the combined approach of computational singular perturbation (CSP) and tangential stretching rate (TSR) to investigate chemical characteristics in blow-off dynamics. To assess the diagnostic approaches in flame and blow-off dynamics, Damköhler number and TSR variables are applied and compared. Four cases are considered in this study showing different flame dynamics such as the steadily stable mode, local extinction by asymmetric vortex shedding, convective blow-off and lean blow-out. DNS data points in positive explosive eigenvalue conditions were subdivided into four different combinations in TSR and extended TSR space and categorized in four distinct characteristic regions, such as kinetically explosive or dissipative and transport-enhanced or dissipative dynamics. The TSR analysis clearly captures the local extinction point in the complicated vortex shedding and allows an improved understanding of the distinct chemistry-transport interactions occurring in convective blow-off and lean blow-out events.     

Validation of a mixture-averaged thermal diffusion model for premixed lean hydrogen flames

The mixture-averaged thermal diffusion model originally proposed by Chapman and Cowling is validated using multiple flame configurations. Simulations using detailed hydrogen chemistry are done on one-, two-, and three-dimensional flames. The analysis spans flat and stretched, steady and unsteady, and laminar and turbulent flames. Quantitative and qualitative results using the thermal diffusion model compare very well with the more complex multicomponent diffusion model. Comparisons are made using flame speeds, surface areas, species profiles, and chemical source terms. Once validated, this model is applied to three-dimensional laminar and turbulent flames. For these cases, thermal diffusion causes an increase in the propagation speed of the flames as well as increased product chemical source terms in regions of high positive curvature. The results illustrate the necessity for including thermal diffusion, and the accuracy and computational efficiency of the mixture-averaged thermal diffusion model.     

Experimental and numerical study of premixed, lean ethylene flames

Ethylene is a key intermediate in the combustion mechanisms of most practical fuels. It plays also an important role in the formation of aromatic hydrocarbons and soot particules. The latter has motivated many experimental and numerical studies carried out on rich ethylene-air mixtures. Less studies have been devoted to lean mixtures, and the development of strategies based on lean, premixed flames to reduce soot and NO_x production requires additional experimental data in lean conditions. In this work, the chemical structure of lean premixed ethylene-oxygen-nitrogen flames stabilized on a flat-flame burner at atmospheric pressure was determined experimentally. The species mole fraction profiles were also computed by the Premix code (Chemkin II version) and four detailed reaction mechanisms. A very good agreement was observed for the main flame properties: reactants consumption, final products (CO₂, H₂O) and the main intermediates: CO and H₂. Marked differences occurred in the prediction of active intermediate species present in small concentrations. Pathways analyses were performed to identify the origins of these discrepancies. It was shown that the same reactions were involved in the four mechanisms to describe the consumption of ethylene, but with marked differences in their relative importance. C₂H₃ and CH₂HCO are the main radicals formed in this first step and their consumption increases the differences between the mechanisms either by the use of different kinetic data for common reactions or by differences in the nature of the consumption reactions.     

Computational study of lean premixed turbulent flames using RANSPDF and LESPDF methods

A computational study is performed on a series of four piloted, lean, premixed turbulent jet flames. These flames use the Sydney Piloted Premixed Jet Burner (PPJB), and with jet velocities of 50, 100, 150 and 200 m/s are denoted PM150, PM1100, PM1150 and PM1200, respectively. Calculations are performed using the RANSPDF and LESPDF methodologies, with different treatments of molecular diffusion, with detailed chemistry and flamelet-based chemistry modelling, and using different imposed boundary conditions. The sensitivities of the calculations to these different aspects of the modelling are compared and discussed. Comparisons are made to experimental data and to previously-performed calculations. It is found that, given suitable boundary conditions and treatment of molecular diffusion, excellent agreement between the calculations and experimental measurements of the mean and variance fields can be achieved for PM150 and PM1100. The application of a recently developed implementation of molecular diffusion results in a large improvement in the computed variance fields in the LESPDF calculations. The inclusion of differential diffusion in the LESPDF calculations provides insight on the behaviour in the near-field region of the jet, but its effects are found to be confined to this region and to the species CO, OH and H₂. A major discrepancy observed in many previous calculations of these flames is an overprediction of reaction progress in PM1150 and PM1200, and this discrepancy is also observed in the LESPDF calculations; however, a parametric study of the LESPDF mixing model reveals that, with a sufficiently large mixing frequency, calculations of these two flames are capable of yielding improved reaction progress in good qualitative agreement with the mean and RMS scalar measurements up to an x/D of 30. Lastly, the merits of each computational methodology are discussed in light of their computational costs.     

Strain rates,flow patterns and flame surface densities in hydrodynamically unstable,weakly turbulent premixed flames

Recent numerical and experimental studies have unveiled a potentially marked difference between the laminar as well as turbulent propagation of premixed flames exhibiting Darrieus–Landau (DL) (or hydrodynamic) instabilities from flames for which instabilities are inhibited. In this study we utilize two-dimensional numerical simulations of slot burner flames as well as experimental Propane–Air Bunsen flames to analyse differences in turbulent propagation, strain rate and induced flow patterns of hydrodynamically stable and unstable flames. We also investigate the effects of hydrodynamic instability on quantities which are directly related to reaction rate closure models, such as flame surface density and stretch factor. A clear enhancement of turbulent flame speed can be observed for unstable flames, generally mitigated at higher turbulence intensity, which is attributed to a flame area increase induced by the characteristic cusp-like DL-induced corrugation, absent in stable flames, which occurs concurrently and in synergy with turbulent wrinkling. Unstable flames also exhibit, both numerically and experimentally, a different correlation between strain rate and flame curvature and are observed to give rise to a channeling of the induced flow in the fresh mixture. Conditionally averaged flame surface density is also observed to attain smaller values in unstable flames, as a result of the thicker turbulent flame brush, indicating that closure models should incorporate instability-related parameters in addition to turbulence-related parameters.     

The unstable behavior of cellular premixed flames induced by intrinsic instability

The unstable behavior of cellular premixed flames induced by intrinsic instability is studied by two-dimensional unsteady calculations of reactive flows. In the present numerical simulation, the compressible Navier–Stokes equation including a one-step irreversible chemical reaction is employed. We consider two basic types of phenomena to account for the intrinsic instability of premixed flames, i.e., hydrodynamic and diffusive-thermal effects. The hydrodynamic effect is caused by the thermal expansion through the flame front; the diffusive-thermal effect is caused by the preferential diffusion of mass versus heat. A disturbance with several wavelength components is superimposed on a planar flame, and the formation of a cellular flame induced by hydrodynamic and diffusive-thermal effects is numerically simulated. After the cellular-flame formation, the combination and division of cells are observed. The behavior of cellular-flame fronts becomes more unstable when the Lewis number is lower than unity, since the diffusive-thermal effect has a great influence on the unstable behavior. The cell size changes with time, and its average is greater than the critical wavelength and becomes smaller by decreasing the Lewis number. The flame velocity of cellular flames depends strongly on the length of computational domain in the direction tangential to the flame front. As the length of computational domain increases, the flame velocity becomes larger. This is because the long-wavelength components of disturbances play an important role in the shape of cellular flames, i.e., in the flame-surface area.     

Flame front analysis of high-pressure turbulent lean premixed methane–air flames

An experimental study on lean turbulent premixed methane–air flames at high pressure is conducted by using a turbulent Bunsen flame configuration. A single equivalence ratio flame at Φ = 0.6 is explored for pressures ranging from atmospheric pressure to 0.9 MPa. LDA measurements of the cold flow indicate that turbulence intensities and the integral length scale are not sensitive to pressure. Due to the decreased kinematic viscosity with increasing pressure, the turbulent Reynolds numbers increase, and isotropic turbulence scaling relations indicate a large decrease of the smallest turbulence scales. Available experimental results and PREMIX code computations indicate a decrease in laminar flame propagation velocities with increasing pressure, essentially between the atmospheric pressure and 0.5 MPa. The u′/S_L ratio increases therefore accordingly. Instantaneous flame images are obtained by Mie scattering tomography. The images and their analysis show that pressure increase generates small scale flame structures. In an attempt to generalize these results, the variance of the flamelet curvatures, the standard deviation of the flamelet orientation angle, and the flamelet crossing lengths have been plotted against which is proportional to the ratio between the integral and Taylor length scales, and which increases with pressure. These three parameters vary linearly with the ratio between large and small turbulence scales and clearly indicate the strong effect of this parameter on premixed turbulent flame dynamics and structure. An obvious consequence is the increase in flame surface density and hence burning rate with pressure, as confirmed by its direct determination from 2D tomographic images.     

Chemical structure of atmospheric pressure premixed laminar formic acid/hydrogen flames

The work presents an experimental and kinetic modeling study of laminar premixed formic acid [HC(O)OH]/H₂/O₂/Ar flames at different equivalence ratios (φ=0.85, 1.1 and 1.3) stabilized on a flat burner at atmospheric pressure, as well as laminar flame speed of HC(O)OH/O₂/Ar flames (φ=0.5–1.5) at 1 atm. Flame structure as well as laminar flame speed were simulated using three different detailed chemical kinetic mechanisms proposed for formic acid oxidation. The components in the fuel blends show different consumption profiles, namely, hydrogen is consumed slower than formic acid. According to kinetic analysis, the reason of the observed phenomenon is that the studied flames have hydrogen as a fuel but also as an intermediate product formed from HC(O)OH decomposition. Comparison of the measured and simulated flame structure shows that all the mechanisms satisfactorily predict the mole fraction profiles of the reactants, main products, and intermediates. It is noteworthy that the mechanisms proposed by Glarborg et al., Konnov et al. and the updated AramcoMech2.0 adequately predict the spatial variations in the mole fractions of free radicals, such as H, OH O and HO₂. However, some drawbacks of the mechanisms used were identified; in particular, they predict different concentrations of CH₂O. As for laminar flame speed simulations, the Konnov et al. mechanism predicts around two times higher values than in experiment, while the Glarborg et al. and updated AramcoMech2.0 show good agreement with the experimental data.     

Extinction strain rates of premixed ammonia/hydrogen/nitrogen-air counterflow flames

Chemical energy vectors will play a crucial role in the transition of the global energy system, due to their essential advantages in storing energy in form of gaseous, liquid, or solid fuels. Ammonia (NH₃) has been identified as a highly promising candidate, as it is carbon-free, can be stored at moderate pressures, and already has a developed distribution infrastructure. As a fuel NH₃ has poor combustion properties that can be improved by the addition of hydrogen, which can be obtained energy-efficiently by partially cracking ammonia into hydrogen (H₂) and nitrogen (N₂) prior to the combustion process. The resulting NH₃/H₂/N₂ blend leads to significantly improved flame stability and resilience to strain-induced blow-out, despite similar laminar flame properties compared to equivalent methane/air flames. This study reports the first measurements of extinction strain rates, measured using the premixed twin-flame configuration in a laminar opposed jet burner, for two NH₃/H₂/N₂ blends over a range of equivalence ratios. Local strain rates are measured using particle tracking velocimetry (PTV) and are related to the inflow conditions, such that the local strain rate at the extinction point can be approximated. The results are compared with 1D-simulations using three recent kinetic mechanisms for ammonia oxidation. By relating the extinction strain rates to laminar flame properties of the unstretched flame, a comparison of the extinction behaviour of CH₄ and NH₃/H₂/N₂ blends can be made. For lean mixtures, NH₃/H₂/N₂-air flames show a significant higher extinction resistance in comparison to CH₄/air. In addition, a strong non-linear dependence between the resistance to extinction and equivalence ratio for NH₃/H₂/N₂ blends is observed.     

Hydrodynamic and chemical scaling for blow-off dynamics of lean premixed flames stabilized on a meso-scale bluff-body

Direct numerical simulations were conducted to investigate the effect of two parameters, density ratio and laminar flame speed, on the conditions of the onset of local extinction and blow-off of lean premixed flames, stabilized on a meso-scale bluff-body in hydrogen-air and syngas–air mixtures. A total of six simulation cases were considered as isolated comparison of the two parametric effects of the fluid dynamic instability and flame time scale. For all cases under study, the general flame development towards the blow-off limit showed a sequence of five distinct modes, with possible cyclic patterns among the different modes for a range of velocity conditions. The onset of local extinction was observed during the asymmetric vortex shedding and vortex street mode. As the density ratio is decreased, the flow inunder reviewstability is promoted through the increased sinuous mode, and such behavior was properly scaled by the Strouhal number. Although the blow-off velocity is altered by the fluid dynamic effects, the condition for the onset of local extinction and blow-off was mainly dictated by the competition between flow residence time associated with the lateral flame motion and ignition delay of the local mixtures. Time scale analysis supported the validity of the findings across all the cases investigated.     

Phenomenological model of chain-branching premixed flames

In this work, we introduce a global kinetic model that includes fuel, oxygen, products and two radical species involved in the reversible chain-branching, chain-propagation and chain-termination reactions. The model naturally extends the Zeldovich–Liñán and Zeldovich–Barenblatt–Dold models and can be used to describe both premixed and diffusion flames. Here it is applied to the problem of the deflagration wave propagation in the hydrocarbon fuel/air mixture with arbitrary equivalence ratio under the simplifying thermal-diffusive approximation. The conservation equations are solved numerically in order to obtain the velocity and structure of the combustion wave. It is demonstrated that the peak values of the adiabatic flame temperature and deflagration velocity are shifted towards the rich mixture composition if the reverse reactions of product decomposition are taken into account. The dependence of the flame speed and temperature on parameters of the system is analysed. The prospects of further investigation are discussed.     

Pressure disturbances and strained premixed flames

A mathematical model is presented for pressure interactions with premixed flames in a prescribed strained velocity field. A stability analysis is carried out including pressure disturbances and an approximate frequency condition obtained. For positive strain the unsteady analysis indicates that the pulsating instability is suppressed. However, for a converging flow (negative strain), the instability is encouraged. Furthermore, the change of the neutral stability boundary in parameter space is explored, showing that a sharp pressure reduction makes the pulsating instability much more accessible.     

Flame fingers and interactions of hydrodynamic and thermodiffusive instabilities in laminar lean hydrogen flames

Lean hydrogen/air flames are prone to hydrodynamic and thermodiffusive instabilities. In this work, the contribution of each instability mechanism is quantified separately by performing detailed simulations of laminar planar lean hydrogen/air flames with different diffusivity models and equations of state to selectively suppress the hydrodynamic or thermodiffusive instability mechanism.From the analysis of the initial phase of the simulations, the thermodiffusive instability is shown to dominate the flame dynamics. If differential diffusion and, hence, the thermodiffusive instability is suppressed, the flame features a strong reduction of the instability growth rates, whereas if present, a wide range of unstable wave numbers is observed due to the strong destabilizing nature of differential diffusion. When instabilities are fully developed, lean hydrogen/air flames feature the formation of small-scale cellular structures and large-scale flame fingers. While the size of the former is known to be close to the most unstable wave length of a linear stability analysis, this work shows that flame fingers also originate from the thermodiffusive instability and most noteworthy, are not linked to an interaction of the two instability mechanisms. They are stable with respect to external perturbations and feature an enhanced flame propagation as the formation of a central cusp at their tip enables the co-existence of two strongly curved leading edges with high reactivity. The thermodiffusive instability is shown to significantly affect the flames’ consumption speed, while the consumption speed enhancement caused by the hydrodynamic instability is significantly smaller. Further, the surface area increase due to wrinkling is strongly diminished if one of the two instability mechanisms is missing. This is linked to a synergistic interaction between the two mechanisms, as the propagation of flame fingers is enhanced by the presence of the hydrodynamic instability due to a widening of the streamlines ahead of the flame fingers.     

Evaluation of heat release indicators in lean premixed H2/Air cellular tubular flames

Heat release rate in combustion systems must be understood in order to control thermoacoustic instabilities, flame extinction, and heat losses. Traditionally OH chemiluminescence (OH*) is used to trace heat release rate (HRR) in H₂/air flames, but its accuracy as a tracer has not been assessed. Lean premixed H₂/air cellular tubular flames are a good test case to evaluate HRR tracers due to the presence of highly reactive flame cells surrounded by regions of near extinction. Comparing the calculated heat release rate to OH* concentration, one finds that [OH*] profiles correlate with the regions of high reactivity (flame cells) but the correlation fails in the low reactivity regions where the HRR is much higher than the [OH*] value indicates. Alternate HRR tracers including [H] and pixel-by-pixel products of [O₂]x[H], [OH]x[H₂], and [O]x[H₂] are analyzed with detailed numerical simulations. The chosen products derive from the main chain reaction steps that contribute to overall HRR in lean, premixed H₂/air flames. Findings suggest that [H] is an accurate yet simple way of tracking HRR. Planar measurements of HRR are possible if LIF measurements of [H] are improved.     

PDF calculations of piloted premixed jet flames

A series of piloted premixed jet flames with strong finite-rate chemistry effects is studied using the joint velocity-turbulence frequency-composition PDF method. The numerical accuracy of the calculations is demonstrated, and the calculations are compared to experimental data. It is found that all calculations show good agreement with the measurements of mean and rms mixture fraction fields, while the reaction progress is overpredicted to varying degrees depending on the jet velocity. In the calculations of the flame with the lowest jet velocity, the species and temperature show reasonable agreement with the measurements, with the exception of a small region near the centerline where products and temperature are overpredicted and fuel and oxidizer are underpredicted. In the calculations of the flame with the highest jet velocity, however, the overprediction of products and temperature and underprediction of fuel and oxidizer is far more severe. An extensive set of sensitivity studies on inlet boundary conditions, turbulence model constants, mixing models and constants, radiation treatment, and chemical mechanisms is conducted to show that any parameter variation offers little improvement from the base case. To shed light on these discrepancies, diagnostic calculations are performed in which the chemical reactions are artificially slowed. These diagnostic calculations serve to validate the experimental data and to quantify the amount by which the base case calculations overpredict reaction progress. Improved calculations of this flame are achieved only through artificially slowing down the chemical reaction by a factor of about 10. The mixing model behavior in this combustion regime is identified as a likely cause for the observed discrepancy in reaction progress.