Effect of Lewis number on premixed laminar lean-limit flames stabilized on a bluff body

In this paper, we present a study on the effect of Lewis number, Le, on the stabilization and blow-off of laminar lean limit premixed flames stabilized on a cylindrical bluff body. Numerical simulations and experiments are conducted for propane, methane and two blends of hydrogen with methane as fuel gases, containing 20% and 40% of hydrogen by volume, respectively. It is found that the Le?>?1 flame blows-off via convection from the base of the flame (without formation of a neck) when the conditions for flame anchoring are not fulfilled. Le?≤?1 flames exhibit a necking phenomenon just before lean blow-off. This necking of the flame front is a result of the local reduction in mass burning rates causing flame merging and quenching of the thin flame tube formed. The structure of these flames at the necking location is found to be similar to tubular flames. It is found that extinction stretch rates for tubular flames closely match values at the neck location of bluff-body flames of corresponding mixtures, suggesting that excessive flame stretch is directly responsible for blow-off of the studied Le?≤?1 flames. After quenching of the neck, the upstream part forms a steady and stable residual flame in the wake of the bluff body while the downstream part is convected away.     

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.     

2D computations of FREI with cool flames for n-heptane/air mixture

Two-dimensional axisymmetric numerical simulation reproduced flames with repetitive extinction and ignition (FREI) in a micro flow reactor with a controlled temperature profile with a stoichiometric n-heptane/air mixture, which have been observed in the experiment. The ignition of hot flame occurred from consumption reactions of CO that was remained in the previous cycle of FREI. Between extinction and ignition locations of hot flames, several other heat release rate peaks related to cool and blue flames were observed for the first time. After the extinction of the hot flame, cool flame by the low-temperature oxidation of n-heptane appeared first and was stabilized in a low wall temperature region. In the downstream of the stable cool flame, a blue flame by the consumption reactions of cool flame products of CH₂O and H₂O₂ appeared. After that, the hot flame ignition occurred from the remaining CO in the downstream of the blue flame. Then after the next hot flame ignition, the blue flame was swept away by the propagating hot flame. Soon before the hot flame merged with the stable cool flame, the hot flame propagation was intensified by the cool flame. After the hot flame merged with the stable cool flame, the hot flame reacted with the incoming fresh mixture of n-C₇H₁₆ and O₂.     

发射CT法测量层流乙烯扩散火焰中温度与烟黑浓度分布的实验研究

对于非均匀吸收、发射、无散射的轴对称含烟黑火焰对象,常规双色法不再适用。本文基于烟黑辐射特性,利用烟黑单色辐射强度图像信息采用CT算法同时重建含烟黑火焰温度与烟黑浓度分布,对层流乙烯扩散火焰的温度与烟黑容积份额进行测量,得到了较好的结果。     

Experimental comparison of opposed and concurrent flame spread in a forced convective microgravity environment

Flame spread experiments in both concurrent and opposed flow have been carried out in a 5.18-s drop tower with a thin cellulose fuel. Flame spread rate and flame length have been measured over a range of 0–30 cm/s forced flow (in both directions), 3.6–14.7 psia, and oxygen mole fractions 0.24–0.85 in nitrogen. Results are presented for each of the three variables independently to elucidate their individual effects, with special emphasis on pressure/oxygen combinations that result in earth-equivalent oxygen partial pressures (normoxic conditions). Correlations using all three variables combined into a single parameter to predict flame spread rate are presented. The correlations are used to demonstrate that opposed flow flames in typical spacecraft ventilation flows (5–20 cm/s) spread faster than concurrent flow flames under otherwise similar conditions (pressure, oxygen concentration) in nearly all spacecraft atmospheres. This indicates that in the event of an actual fire aboard a spacecraft, the fire is likely to grow most quickly in the opposed mode as the upstream flame spreads faster and the downstream flame is inhibited by the vitiated atmosphere produced by the upstream flame. Additionally, an interesting phenomenon was observed at intermediate values of concurrent forced flow velocity where flow/flame interactions produced a recirculation downstream of the flame, which allowed an opposed flow leading edge to form there.     

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.     

Temperature and CH* measurements and simulations of laminar premixed ethylene jet-wall stagnation flames

New experimental 2D measurements are reported to characterise the flame location, shape and temperature of laminar premixed ethylene jet-wall stagnation flames when the equivalence ratio, exit gas velocity and burner-plate separation distance are varied. Bandpass-filtered optical measurements of the CH* chemiluminescence were used to provide information about the shape and location of the flames. Thin filament pyrometry (TFP) using a 14 µm diameter SiC filament was used to make line measurements of the temperature to reconstruct the full 2D temperature field for the first time in premixed, jet-wall stagnation flames. The comparison of CH* measurements with (intrusive) and without (non-intrusive) the presence of the SiC filament showed that the filament resulted in minimal disturbance of the flame when the filament was placed downstream of the flame front. However, the flame was observed to attach to the filament, resulting in more significant disturbance, when it was placed upstream of the flame front. The flames were simulated using both 1D and 2D models. The 2D simulations were used to provide estimates of the velocity, kinematic viscosity and thermal conductivity required to calculate the gas temperature from the TFP data. The 1D simulations showed excellent agreement with the experimentally observed centreline quantities, but required the strain boundary condition to be fitted in order to match the experimentally observed flame location. The 2D simulations showed excellent agreement without the need for any fitting, and correctly predicted the flame shape, location and temperature as the experimental conditions were varied. A comparison of the set of simulated temperature-residence times along different streamlines showed relatively uniform distributions within each flame. However, the most uniform set of temperature-residence time distributions did not correlate with the flattest flame.     

Framework for simulating stationary spherical flames

Understanding and quantifying the effects of flame stretch rate on the laminar flame speed and flame structure plays an important role from interpreting experimentally-measured laminar burning velocities to characterizing the impact of turbulence on premixed flames. Unfortunately, accounting for these effects often requires an unsteady reacting flow solver and may be computationally expensive. In this work, we propose a mathematical framework to perform simulations of stationary spherical flames. The objective is to maintain the flame at a constant radius (and hence a constant stretch rate) by performing a coordinate change. The governing equations in the new flame-attached frame of reference resemble the original equations for freely-propagating spherical flames. The only difference is the presence of additional source terms whose purpose is to drive the numerical solution to a steady state. These source terms involve one free parameter: the flame stretch rate, which may either be computed in real time or imposed by the user. This parameter controls ultimately the steady state flame radius and the steady state flame speed. That is why, at a given stretch rate, the results of the stationary spherical flame simulations match those of a freely-expanding spherical flame. As an illustration, the dependence of the laminar flame speed on the stretch rate is leveraged to extract Markstein lengths for hydrogen/air mixtures at different equivalence ratios, as well as for hydrocarbon/air mixtures (CH₄ and C₇H₁₆). Numerical predictions are in good agreement with experimental measurements (within experimental uncertainties). Finally, the proposed methodology is implemented in the chemical kinetic software FlameMaster. The use of a dedicated steady-state solver with a non-uniform optimized mesh leads to significant reductions in the computational cost, highlighting that the proposed methodology is ideally suited for other chemical kinetic software such as Chemkin/Premix and Cantera.     

Influence of radiation losses on the stability of premixed flames on a porous-plug burner

Analysis of the planar premixed flames on a porous plug was performed numerically for finite activation energy within the diffusive-thermal model. The paper is focused on the influence of radiation heat loses on the flame standoff distance and its linear stability. We show that the presence of volumetric heat losses limits the range of the mass flow range as well as it can promote the flame instabilities of different kinds, both oscillatory and cellular. The oscillatory instability, which for freely propagating flames can be usually observed for the Lewis number larger than one, in the porous-plug case occurs also for flames with unity and lower than unity Lewis number. For flames with Le < 1 both cellular and oscillatory instabilities can be observed simultaneously in a certain range of the mass flow rate.     

A numerical investigation of stretch effects in counterflow,premixed laminar flames

We use direct numerical simulation of propane/air flames with full chemistry in the geometry of stagnation flow to investigate the effect of different definitions of local flame stretch in the presence of spatially varying velocity gradients. Specifically, we compare simulations with potential- and plug-flow inlet conditions, and show that the widely used definition of upstream stretch leads to unphysical results for flames having the ‘same’ stretch. We then show that a reasonable re-definition of local stretch allows us to produce the ‘same’ flame in the presence of the ‘same’ stretch.     

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.     

How fast can we burn, 2.0

We review the state of the art in measurements and simulations of the behavior of premixed laminar and turbulent flames, subject to differential diffusion, stretch and curvature. The first part of the paper reviews the behavior of premixed laminar flames subject to flow stretch, and how it affects the accuracy of measurements of unstrained laminar flame speeds in stretched and spherically propagating flames. We then examine how flow field stretch and differential diffusion interact with flame propagation, promoting or suppressing the onset of thermodiffusive instabilities. Secondly, we survey the methodology for and results of measurements of turbulent flame speeds in the light of theory, and identify issues of consistency in the definition of mean flame speeds, and their corresponding mean areas. Data for methane at a single operating condition are compared for a range of turbulent conditions, showing that fundamental issues that have yet to be resolved for Bunsen and spherically propagating flames. Finally, we consider how the laminar flame scale response of flames to flow perturbations interacting with differential diffusion leads to very different outcomes to the overall sensitivity of the burning rate to turbulence, according to numerical simulations (DNS). The paper concludes with opportunities for future measurements and model development, including the perennial recommendation for robust archival databases of experimental and DNS results for future testing of models.     

Large scale effects in weakly turbulent premixed flames

In this study we numerically investigate large scale premixed flames in weakly turbulent flow fields. A large scale flame is classified as such based on a reference hydrodynamic lengthscale being larger than a neutral (cutoff) lengthscale for which the hydrodynamic or Darrieus–Landau (DL) instability is balanced by stabilizing diffusive effects. As a result, DL instability can develop for large scale flames and is inhibited otherwise. Direct numerical simulations of both large scale and small scale three-dimensional, weakly turbulent flames are performed at constant Karlovitz and turbulent Reynolds number, using two paradigmatic configurations, namely a statistically planar flame and a slot Bunsen flame. As expected from linear stability analysis, DL instability induces its characteristic cusp-like corrugation only on large scale flames. We therefore observe significant morphological and topological differences as well as DL-enhanced turbulent flame speeds in large scale flames. Furthermore, we investigate issues related to reaction rate modeling in the context of flame surface density closure. Thicker flame brushes are observed for large scale flames resulting in smaller flame surface densities and overall larger wrinkling factors.     

A computational analysis of strained laminar flame propagation in a stratified CH4/H2/air mixture

Propagation of a H₂-added strained laminar CH₄/air flame in a rich-to-lean stratified mixture is numerically studied. The back-support effect, which is known to enhance the consumption speed of a flame propagating into a leaner mixture compared to that into a homogeneous mixture, is evaluated. A new method is devised to characterize unsteady reactant-to-reactant counterflow flames under transiently decreasing equivalence ratio, in order to elucidate the influence of flow strain on the back-support effect. In contrast to the conventional reactant-to-product configurations, the current configuration is more relevant to unsteady stratified flames back-supported by their own combustion products. Moreover, since H₂ distribution downstream of the flame is known to play a crucial role in back-supported CH₄/air flames, the influence of H₂ addition in the upstream mixture is examined. The results suggest that a larger strain rate leads to a larger equivalence ratio gradient at the reaction zone through increased flow divergence, which amplifies the back-support. Meanwhile, since H₂ addition in the upstream mixture does not affect the downstream H₂ content, the relative increase in the consumption speed, i.e. the back-support, is suppressed with larger H₂ addition. Especially, when the upstream H₂ content decreases with the equivalence ratio, the H₂ preferentially diffuses toward the unburned gas, which mitigates H₂ accumulation in the preheat zone and further weakens the back-support.     

Critical temperature and reactant mass flux for radiative extinction of ethylene microgravity spherical diffusion flames at 1 bar

In microgravity combustion, where buoyancy is not present to accelerate the flow field and strain the flame, radiative extinction is of fundamental importance, and has implications for spacecraft fire safety. In this work, the critical point for radiative extinction is identified for normal and inverse ethylene spherical diffusion flames via atmospheric pressure experiments conducted aboard the International Space Station, as well as with a transient numerical model. The fuel is ethylene with nitrogen diluent, and the oxidizer is an oxygen/nitrogen mixture. The burner is a porous stainless-steel sphere. All experiments are conducted at constant reactant flow rate. For normal flames, the ambient oxygen mole fraction was varied from 0.2 to 0.38, burner supply fuel mole fraction from 0.13 to 1, total mass flow rate, ṁ_total, from 0.6 to 12.2 mg/s, and adiabatic flame temperature, T_ad, from 2000 to 2800 K. For inverse flames, the ambient fuel mole fraction was varied from 0.08 to 0.12, burner supply oxygen mole fraction from 0.4 to 0.85, ṁ_total from 2.3 to 11.3 mg/s, and T_ad from 2080 to 2590 K. Despite this broad range of conditions, all flames extinguish at a critical extinction temperature of 1130 K, and a fuel-based mass flux of 0.2 g/m²-s for normal flames, and an oxygen-based mass flux of 0.68 g/m²-s for inverse flames. With this information, a simple equation is developed to estimate the flame size (i.e., location of peak temperature) at extinction for any atmospheric-pressure ethylene spherical diffusion flame given only the reactant mass flow rate. Flame growth, which ultimately leads to radiative extinction if the critical extinction point is reached, is attributed to the natural development of the diffusion-limited system as it approaches steady state and the reduction in the transport properties as the flame temperature drops due to increasing flame radiation with time (radiation-induced growth.)     

Ultra-lean hydrogen-enriched oscillating flames behind a heat conducting bluff-body: Anomalous and normal blow–off

This paper presents a numerical study of ultra-lean hydrogen-methane flames stabilized behind a rectangular, highly conducting metallic bluff body acting as a flame holder. Using high fidelity numerical simulations, we show that lean inverted steady flames exist below normal flammability limits. They have distinct stabilization mechanism from pure methane flames. These flames are blown-off for sufficiently small velocities, a phenomenon called anomalous blow-off. At even leaner conditions oscillating ultra–lean hydrogen-methane flames can be established. These oscillating flames exist within a rather small range of equivalence ratios and inflow velocities, and move to mean locations closer to the flame holder as the reactant flow is increased. We show that the oscillations are associated with the shedding of flame balls from the downstream end of a “residual flame” that remains attached. Unlike their steady counterparts, the oscillating flames exhibit blow-off at both low velocities (anomalous blow-off) and at sufficiently high inflow velocities (normal blow-off). We show that normal blow-off is linked to heat losses to the flame holder.     

Minimum ignition energy and propagation dynamics of laminar premixed cool flames

Laminar premixed cool flames, induced by the coupling of low-temperature chemistry and convective-diffusive transport process, have recently attracted extensive interest in combustion and engine research. In this work, numerical simulations have been conducted using a recently developed open-source reacting flow platform reactingFOAM-SCT, to investigate the minimum ignition energy (MIE) and propagation dynamics of premixed cool flames in a 1D spherical coordinate. Results have shown that when ignition energy is below the MIE of regular hot flames, a class of cool flames could be initiated, which allow much wider flammability limits, both lean and rich, compared to hot flames. Furthermore, the overall cool flame propagation dynamics exhibit intrinsic similarity to those of hot flames, in that, they begin with an ignition kernel propagation regime, followed by two transition regimes, and eventually reach a normal flame propagation regime. However, a spherical expanding cool flame responds completely differently to stretch. Specifically, a regular outwardly propagating hot spherical flame accelerates with increasing stretch rate when the mixture Le < 1 and decelerates when Le > 1. However, it is found that a cool flame always tends to decelerate with increasing stretch rate regardless of mixture composition, exhibiting unique flame aerodynamic characteristic. This research discovers novel features of premixed cool flame initiation and propagation dynamics and sheds light on flame transition, spark-ignition system design, and advanced engine combustion control.     

Response of non-premixed flames to bulk flow perturbations

This paper describes the dynamics of non-premixed flames responding to bulk velocity fluctuations, and compares the dynamics of the flame sheet position and spatially integrated heat release to that of a premixed flame. The space–time dynamics of the non-premixed flame sheet in the fast chemistry limit is described by the stoichiometric mixture fraction surface, extracted from the solution of the

-equation. This procedure has some analogies to premixed flames, where the premixed flame sheet location is extracted from the G = 0 surface of the solution of the G-equation. A key difference between the premixed and non-premixed flame dynamics, however, is the fact that the non-premixed flame sheet dynamics are a function of the disturbance field everywhere, and not just at the reaction sheet, as in the premixed flame problem. A second key difference is that the non-premixed flame does not propagate and so flame wrinkles are convected downstream at the axial flow velocity, while wrinkles in premixed flames convect downstream at a vector sum of the flame speed and axial velocity. With the exception of the flame wrinkle propagation speed, however, we show that that the solutions for the space–time dynamics of the premixed and non-premixed reaction sheets in high velocity axial flows are quite similar. In contrast, there are important differences in their spatially integrated unsteady heat release dynamics. Premixed flame heat release fluctuations are dominated by area fluctuations, while non-premixed flames are dominated by mass burning rate fluctuations. At low Strouhal numbers, the resultant sensitivity of both flames to flow disturbances is the same, but the non-premixed flame response rolls off slower with frequency. Hence, this analysis suggests that non-premixed flames are more sensitive to flow perturbations than premixed flames at O(1) Strouhal numbers.     

A DNS study of the physical mechanisms associated with density ratio influence on turbulent burning velocity in premixed flames

Data obtained in 3D direct numerical simulations of statistically planar, 1D weakly turbulent flames characterised by different density ratios σ are analysed to study the influence of thermal expansion on flame surface area and burning rate. Results show that, on the one hand, the pressure gradient induced within a flame brush owing to heat release in flamelets significantly accelerates the unburned gas that deeply intrudes into the combustion products in the form of an unburned mixture finger, thus causing large-scale oscillations of the burning rate and flame brush thickness. Under the conditions of the present simulations, the contribution of this mechanism to the creation of the flame surface area is substantial and is increased by σ, thus implying an increase in the burning rate by σ. On the other hand, the total flame surface areas simulated at σ = 7.53 and 2.5 are approximately equal. The apparent inconsistency between these results implies the existence of another thermal expansion effect that reduces the influence of σ on the flame surface area and burning rate. Investigation of the issue shows that the flow acceleration by the combustion-induced pressure gradient not only creates the flame surface area by pushing the finger tip into the products, but also mitigates wrinkling of the flame surface (the side surface of the finger) by turbulent eddies. The latter effect is attributed to the high-speed (at σ = 7.53) axial flow of the unburned gas, which is induced by the axial pressure gradient within the flame brush (and the finger). This axial flow acceleration reduces the residence time of a turbulent eddy in an unburned zone of the flame brush (e.g. within the finger). Therefore, the capability of the eddy for wrinkling the flamelet surface (e.g. the side finger surface) is weakened owing to a shorter residence time.