Structure of propagating edge diffusion flames in hydrocarbon fuel jets

Direct numerical simulations with a C₃-chemistry model have been performed to investigate the transient behavior and internal structure of flames propagating in an axisymmetric fuel jet of methane, ethane, ethylene, acetylene, or propane in normal earth gravity (1g) and zero gravity (0g). The fuel issued from a 3-mm-i.d. tube into quasi-quiescent air for a fixed mixing time of 0.3 s before it was ignited along the centerline where the fuel–air mixture was at stoichiometry. The edge of the flame formed a vigorously burning peak reactivity spot, i.e., reaction kernel, and propagated through a flammable mixture layer, leaving behind a trailing diffusion flame. The reaction kernel broadened laterally across the flammable mixture layer and possessed characteristics of premixed flames in the direction of propagation and unique flame structure in the transverse direction. The reaction kernel grew wings on both fuel and air sides to form a triple-flame-like structure, particularly for ethylene and acetylene, whereas for alkanes, the fuel-rich wing tended to merge with the main diffusion flame zone, particularly methane. The topology of edge diffusion flames depend on the properties of fuels, particularly the rich flammability limit, and the mechanistic oxidation pathways. The transit velocity of edge diffusion flames, determined from a time series of calculated temperature field, equaled to the measured laminar flame speed of the stoichiometric fuel–air mixtures, available in the literature, independent of the gravity level.     

Influences of stoichiometry on steadily propagating triple flames in counterflows

Most studies of triple flames in counterflowing streams of fuel and oxidizer have been focused on the symmetric problem in which the stoichiometric mixture fraction is 1/2. There then exist lean and rich premixed flames of roughly equal strengths, with a diffusion flame trailing behind from the stoichiometric point at which they meet. In the majority of realistic situations, however, the stoichiometric mixture fraction departs appreciably from unity, typically being quite small. With the objective of clarifying the influences of stoichiometry, attention is focused on one of the simplest possible models, addressed here mainly by numerical integration. When the stoichiometric mixture fraction departs appreciably from 1/2, one of the premixed wings is found to be dominant to such an extent that the diffusion flame and the other premixed flame are very weak by comparison. These curved, partially premixed flames are expected to be relevant in realistic configurations. In addition, a simple kinematic balance is shown to predict the shape of the front and the propagation velocity reasonably well in the limit of low stretch and low curvature.     

Effect of ignition energy on the uncertainty in the determination of laminar flame speed using outwardly propagating spherical flames

The initial propagation processes of expanding spherical flames of CH₄/N₂/O₂/He mixtures at different ignition energies were investigated experimentally and numerically to reduce the effect of ignition energy on the accurate determination of laminar flame speeds. The experiments were conducted in a constant-volume combustion bomb at initial pressures of 0.07???0.7?MPa, initial temperatures of 298???398?K, and equivalence ratios of 0.9???1.3 with various Lewis numbers. The A-SURF program was employed to simulate the corresponding flame propagation processes. The results show that elevating the ignition energy increases the initial flame propagation speed and expands the range of flame trajectory which is affected by ignition energy, but the increase rates of the speed and range decrease with the ignition energy. Based on the trend of the minimum flame propagation speed during the initial period with the ignition energy, the minimum reliable ignition energy (MRIE) is derived by considering the initial flame propagation speed and energy conservation. It is observed that MRIE first decreases and then increases with the increasing equivalence ratio and monotonously decreases with increasing initial pressure and temperature. As the Lewis number rises, MRIE increases. The results also suggest that during the data processing of the spherical flame experiment, the accuracy of determination of laminar flame speeds can be enhanced when taking the flame radius influenced by MRIE as the lower limit of the flame radius range. Then the flame radius influenced by MRIE was defined as RFR. It can also be found that there exist nonlinear relationships between RFR and the equivalence ratio and Lewis number, and the RFR decreases with increasing initial pressure and temperature.     

Influence of advective heat flux on extinction scalar dissipation rate and velocity in negative edge flames

A counterflow flame geometry, which has previously been experimentally shown to produce stable negative edge flames, was studied using numerical simulations. In this geometry, the flame edge is formed off the counterflow centreline owing to a local increase in scalar dissipation rate. Hot products from the stable nonpremixed flame on the centreline flow through the edge at velocities of ～ 1–5 m/s. The size of the counterflow burner and the gas flowrates are varied in the simulations to alter the flame strength and velocity at the flame edge. The advection of products through the edge is shown to extend the flame extinction to higher scalar dissipation rates than required for centreline extinction. For high velocities, the scalar dissipation rate required for flame extinction can be related to the centreline extinction value by considering only the effect of energy addition to the flame edge via advection. However, for lower edge flame velocities, the effects of increased thermal and species diffusion through the edge must also be included. Since the advection at the edge is a product of both the local velocity and temperature gradient, a single correlation between the scalar dissipation rate and the negative edge flame velocity does not exist.     

On the determination of laminar flame speed from low-pressure and super-adiabatic propagating spherical flames

The outwardly propagating spherical flame (OPF) method is popularly used to measure the laminar flame speed (LFS). Recently, great efforts have been devoted to improving the accuracy of the LFS measurement from OPF. In the OPF method, several assumptions are made. For examples, the burned gas is assumed to be static and in chemical equilibrium. However, these assumptions may not be satisfied under certain conditions. Here we consider low-pressure and super-adiabatic propagating spherical flames, for which chemical non-equilibrium exists and the burned gas may not be static. The objective is to assess the chemical non-equilibrium effects on the accuracy of LFS measurement from the OPF method. Numerical simulations considering detailed chemistry and transport are conducted. Stoichiometric methane/air flames at sub-atmospheric pressures and methane/oxygen flames at different equivalence ratios are considered. At low pressures, broad heat release zone is observed and the burned gas cannot quickly reach the adiabatic flame temperature, indicating the existence of chemical non-equilibrium of burned gas. Positive flow in the burned gas is identified and it is shown to become stronger at lower initial pressure. Consequently, the LFS measurement from OPF at low pressures is not accurate if the burned gas is assumed to be static and at chemical equilibrium. For super-adiabatic spherical flames, the burned gas speed is found to be negative due to the local temperature overshoot at the flame front. Such negative speed of burned gas can also reduce the accuracy of LFS measurement. It is recommended that the direct method measuring both flame propagation speed and flow speed of unburned gas should be used to determine the LFS at low pressures or for mixtures with super-adiabatic flame temperature.     

Measurement and simulation of partially-premixed cellular tubular flames

Experimental and numerical simulation results are reported of partially-premixed cellular tubular flames. Parametric measurements across stretch rate and equivalence ratio are taken by chemiluminescent imaging and are presented for the first time. Select hybrid cases with both cellular and non-cellular flame structures are examined with laser-induced spontaneous Raman scattering. Results are spatially resolved in two dimensions and radial interpolations of reaction and extinction zones are compared to numerical simulations using multicomponent transport and detailed chemical kinetics. Experimental cell structures and extinction zones are well predicted by numerical simulation, with discrepancies of temperature and H₂O and temperature primarily observed in locations with moderate and high mole fractions of CO₂. A novel cellular structure, denoted as a “split-cell” flame, is reported for the first time with both chemiluminescent imaging and Raman scattering. Results indicate that partially-premixed flames are valuable as experimental and numerical benchmarks to advance fundamental combustion research.     

Experimental measurements of geometric properties of turbulent stratified flames

Combustion under stratified conditions is common in many systems. However, relatively little is known about the structure and dynamics of turbulent stratified flames. Two-dimensional imaging diagnostics are applied to premixed and stratified V-flames at a mean equivalence ratio of 0.77, and low turbulent intensity, within the corrugated flame range. The present results show that stratification affects the mean turbulent flame speed, structure and geometric properties. Stratification increases the flame surface density above the premixed flame levels in all cases, with a maximum reached at intermediate levels of stratification. The flame surface density (FSD) of stratified flames is higher than that of premixed flames at the same mean equivalence ratio. Under the present conditions, the FSD peaks at a stratification ratio around 3.0. The FSD curves for stratified flames are further skewed towards the product side. The distribution of flame curvature in stratified flames is broader and more symmetric relative to premixed flames, indicating an additional mechanism of curvature generation, which is not necessarily due to cusping. These experiments indicate that flame stratification affects the intrinsic behaviour of turbulent flames and suggest that models may need to be revised in the light of the current evidence.     

An assessment of large eddy simulations of premixed flames propagating past repeated obstacles

This paper presents an assessment of Large Eddy Simulations (LES) in calculating the structure of turbulent premixed flames propagating past solid obstacles. One objective of the present study is to evaluate the LES simulations and identify the drawbacks in accounting the chemical reaction rate. Another objective is to analyse the flame structure and to calculate flame speed, generated overpressure at different time intervals following ignition of a stoichiometric propane/air mixture. The combustion chamber has built-in repeated solid obstructions to enhance the turbulence level and hence increase the flame propagating speed. Various numerical tests have also been carried out to determine the regimes of combustion at different stages of the flame propagation. These have been identified from the calculated results for the flow and flame characteristic parameters. It is found that the flame lies within the ‘thin reaction zone’ regime which supports the use of the laminar flamelet approach for modelling turbulent premixed flames. A submodel to calculate the model coefficient in the algebraic flame surface density model is implemented and examined. It is found that the LES predictions are slightly improved owing to the calculation of model coefficient by using submodel. Results are presented and discussed in this paper are for the flame structure, position, speed, generated pressure and the regimes of combustion during all stages of flame propagation from ignition to venting. The calculated results are validated against available experimental data.     

Structure and propagation of two-dimensional,partially premixed,laminar flames in diesel engine conditions

We investigate the influence of inflow velocity (V_in) and scalar dissipation rate (χ) on the flame structure and stabilisation mechanism of steady, laminar partially premixed n-dodecane edge flames stabilised on a convective mixing layer. Numerical simulations were performed for three different χ profiles and several V_in (V_in = 0.2 to 2.5m/s). The ambient thermochemical conditions were the same as the Engine Combustion Network’s (ECN) Spray A flame, which in turn represents conditions in a typical heavy duty diesel engine. The results of a combustion mode analysis of the simulations indicate that the flame structure and stabilisation mechanism depend on V_in and χ. For low V_in the flame is attached. Increasing V_in causes the high-temperature chemistry (HTC) flame to lift-off, while the low-temperature chemistry (LTC) flame is still attached. A unique speed S_R associated with this transition is defined as the velocity at which the lifted height has the maximum sensitivity to changes in V_in. This transition velocity is negatively correlated with χ. Near $V_{in}=S_R$ a tetrabrachial flame structure is observed consisting of a triple flame, stabilised by flame propagation into the products of an upstream LTC branch. Further increasing the inlet velocity changes the flame structure to a pentabrachial one, where an additional HTC ignition branch is observed upstream of the triple flame and ignition begins to contribute to the flame stabilisation. At large V_in, the LTC is eventually lifted, and the speed at which this transition occurs is insensitive to χ. Further increasing V_in increases the contribution of ignition to flame stabilisation until the flame is completely ignition stabilised. Flow divergence caused by the LTC branch reduces the χ at the HTC branches making the HTC more resilient to χ. The results are discussed in the context of identification of possible stabilisation modes in turbulent flames.     

An apparatus-independent extinction strain rate in counterflow flames

Resistance to extinction by stretch is a key property of any flame, and recent work has shown that this property controls the overall structure of several important types of turbulent flames. Multiple definitions of the critical strain rate at extinction (ESR) have been presented in the literature. However, even if the same definition is used, different experiments report different extinction strain rates for flames burning the same fuel-air mixture at very similar temperatures using similarly constructed opposed-flow instruments. Here we show that at extinction, all these flames are essentially identical, so one would expect that each would be assigned the same value of a parameter representing its intrinsic resistance-to-stretch-induced-extinction, regardless of the specifics of the experimental apparatus. A similar situation arises in laminar flame speed measurements since different apparatuses could result in different strain rate distributions. In that instance, the community has agreed to report the unstretched laminar flame speed, and methods have been developed to translate the experimental (stretched) flame speed into a universal unstretched laminar flame speed. We propose an analogous method for translating experimental measurements for stretch-induced extinction into an unambiguous and apparatus-independent quantity (ESR_∞) by extrapolating to infinite opposing burner separation distance. The uniqueness of the flame at extinction is shown numerically and supported experimentally for twin premixed, single premixed, and diffusion flames at Lewis numbers greater than and less than one. A method for deriving ESR_∞ from finite-boundary experimental studies is proposed and demonstrated for methane and propane experimental diffusion and premixed single flame data. The two values agree within the range of ESR differences typically observed between experimental measurements and simulation results for the traditional ESR definition.     

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.     

The effects of dimethyl ether and ethanol on benzene and soot formation in ethylene nonpremixed flames

Soot volume fractions, C1-C12 hydrocarbon concentrations, and gas temperature were measured in ethylene/air nonpremixed flames with up to 10% dimethyl ether (CH₃OCH₃) or ethanol (CH₃CH₂OH) added to the fuel. The measurement techniques were laser-induced incandescence, photoionization mass spectroscopy, and thermocouples. Oxygenated hydrocarbons have been proposed as soot-reducing fuel additives, and nonpremixed flames are good laboratory-scale models of the fuel-rich reaction zones where soot forms in many full-scale combustion devices. However, addition of both dimethyl ether and ethanol increased the maximum soot volume fractions in the ethylene flames studied here, even though ethylene is a much sootier fuel than either oxygenate. Furthermore, dimethyl ether produced a larger increase in soot even though neat dimethyl ether flames produce less soot than neat ethanol flames. The detailed species measurements suggest that the oxygenates increase soot concentrations because they decompose to methyl radical, which promotes the formation of propargyl radical (C₃H₃) through C1 + C2 addition reactions and consequently the formation of benzene through propargyl self-reaction. Dimethyl ether has a stronger effect than ethanol because it decomposes more completely to methyl radical. Ethylene does not decompose to methyl, so its flames are particularly sensitive to this mechanism; the alkane-based fuels used in most practical fuels do decompose to methyl radical, so the mechanism will be much less important for practical devices.     

Stabilization, propagation and instability of tribrachial triple flames

A tribrachial (or triple) flame is one kind of edge flame that can be encountered in nonpremixed mixing layers, consisting of a lean and a rich premixed flame wing together with a trailing diffusion flame all extending from a single point. The flame could play an important role on the characteristics of various flame behaviors including lifted flames in jets, flame propagation in two-dimensional mixing layers, and autoignition fronts. The structure of tribrachial flame suggests that the edge is located along the stoichiometric contour in a mixing layer due to the coexistence of all three different types of flames. Since the edge has a premixed nature, it has unique propagation characteristics. In this review, the propagation speed of tribrachial flames will be discussed for flames propagating in mixing layers, including the effects of concentration gradient, velocity gradient, and burnt gas expansion. Based on the tribrachial edge structure observed experimentally in laminar lifted flames in jets, the flame stabilization characteristics including liftoff height, reattachment, and blowout behaviors and their buoyancy-induced instability will be explained. Various effects on liftoff heights in both free and coflow jets including jet velocity, the Schmidt number of fuel, nozzle diameter, partial premixing of air to fuel, and inert dilution to fuel are discussed. Implications of edge flames in the modeling of turbulent nonpremixed flames and the stabilization of turbulent lifted flames in jets are covered.     

Effect of the reversibility of the chemical reaction on triple flames

We examine a new aspect of triple flames, namely the effect of the reversibility of the chemical reaction on flame propagation. The study is carried out in the configuration of the two-dimensional strained mixing layer formed between two opposing streams of fuel and oxidiser. The chemical reaction is modelled as a single reversible reaction following an Arrhenius law in the forward and backward directions. The problem is formulated within the constant-density (thermo-diffusive) approximation, the main non-dimensional parameters relevant to this study being a reversibility parameter R (equal to zero in the irreversible case), a non-dimensional measure of the strain rate ? and a stoichiometric parameter S. We provide analytical (asymptotic) expressions for the local burning speed of the triple flame in terms of ?, S, and R. In particular we describe how the propagation speed of the front is decreased by an increase in R and how the location of its leading edge is affected by reversibility. For example, it is found that the leading edge shifts towards the fuel stream for S > 1 and towards the oxidiser if S < 1, as R is increased. A detailed numerical study is conducted covering all propagation regimes ranging from weakly stretched positively propagating (ignition) fronts to thick negatively propagating (extinction) fronts. In the weakly stretched cases we show that the numerics are in good agreement with the asymptotic findings. Furthermore, the results allow the determination of the domains of the distinct propagation regimes, mainly in the terms of R and ?. In line with our physical intuition, it is found that reversibility reduces the domain of ignition fronts and promotes extinction. The results provide a systematic investigation which can be considered as a first step when considering a more realistic chemistry, or poorly explored aspects (such as the existence of a temperature gradient in the unburnt mixture), when analyzing triple flames.     

Influence of edge propagation on downward flame spread over three-dimensional PMMA samples

This systematic experimental study measures the effect of flame propagation along vertical edges on the overall downward spread of flames using Polymethyl Methacrylate (PMMA). Samples with a wide range of regular cross-sections – from triangular through octagonal – as well as irregular ones, have been used to test a large variation of internal angles. A MATLAB-based tool was used to calculate instantaneous spread rate for central and edge flames. The edge flame is shown to significantly enhance the spread rate, as much as five times, in respect to samples with no edges. This amplification is shown to depend primarily on the internal angle at the edge (the smaller the angle, the faster the flame) and fuel thickness, and not on other factors such as aspect ratio or cross-sectional area. Using a phenomenological argument, the edge propagation rate is correlated to the spread rate over an equivalent cylindrical fuel (the limiting shape with infinite edges) with an effective radius obtained from the geometry of the edges and the diffusion length scale of the solid phase. A formula for flame spread over cylindrical fuel from the literature is used to link the new results to existing models. Both thick and thin limits are shown to encompass the wide range of three-dimensional spread rate data within the effective radius (the independent variable), which can be determined from the known parameters. Based on these results, different types of cross-sectional areas can be sorted in the order of their inherent fire safety characteristics.     

Detailed measurements of local scalar front structures in stagnation-type turbulent premixed flames

Local scalar front structures of OH mole fraction, reaction progress variable, and its three-dimensional gradient have been measured in stagnation-type turbulent premixed flames. The reaction progress variable front is observed to change with increasing turbulence from parallel iso-scalar contours but reduced progress variable gradients, called the lamella-like front, to disrupted non-parallel iso-contours that deviate substantially from those of wrinkled laminar flamelets, called the non-flamelet front. This transition is attributed to the different scales of interaction between the flame internal structure and a spectrum of turbulence extending from the integral scale to the Kolmogorov scale. The lamella-like front pattern occurs when the length scales of interaction are smaller than the laminar flame thickness but the time scales are greater than the flame residence time. The non-flamelet front pattern occurs when the length scales of interaction are greater than the laminar flame thickness but the time scales are smaller than the flame residence time. This difference corresponds to the change of combustion regime from complex-strain flame front to turbulent flame front on a revised regime diagram. A correlation is also proposed for the turbulent flame brush thickness as a function of turbulent Reynolds number and heat release parameter. The heat release parameter is considered to arise from the non-passive effects of flame-surface wrinkling.     

Polydisperse effects in jet spray flames

A laminar jet polydisperse spray diffusion flame is analysed mathematically for the first time using an extension of classical similarity solutions for gaseous jet flames. The analysis enables a comparison to be drawn between conditions for flame stability or flame blow-out for purely gaseous flames and for spray flames. It is found that, in contrast to the Schmidt number criteria relevant to gas flames, droplet size and initial spray polydispersity play a critical role in determining potential flame scenarios. Some qualitative agreement for lift-off height is found when comparing predictions of the theory and sparse independent experimental evidence from the literature.