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.     

Asymptotic analysis to the effect of temperature gradient on the propagation of triple flames

We study asymptotically in this paper the influence of the temperature gradient across the mixing layer on the propagation triple flames formed inside a porous wall channel. The study begins by formulating the problem mathematically using the thermo-diffusive model and then presents a thorough asymptotic analysis of the problem in the limit of large activation energy and thin flames. Analytical formulae for the local burning speed, the flame shape and the propagation speed in terms of the temperature gradient parameter have been derived. It was shown that varying the feed temperatures can significantly enhance the burning of the reactants up to a critical threshold, beyond which no solutions can be obtained. In addition, the study showed that increasing the temperature at the boundaries will modify the usual triple structure of the flame by inverting the upper premixed branch and extending it to the boundary, which may have great implications on the safety of the adopted combustion chambers.     

Lift-off characteristics of triple flame with concentration gradient

The concentration gradient and uniform mean velocity of a triple flame in a mixing layer were studied using a multi-slot burner, which can stabilize the lift-off flame especially at a very small concentration gradient. Flame stabilization conditions were examined, and the lift-off heights of the triple flame were measured for methane and propane flames. A hot-wire anemometer was used to measure the velocity distributions. Mass spectroscopy (for methane) and Rayleigh scattering (for propane) were used to measure the concentration gradients. OH radical distribution was measured by laser-induced fluorescence (LIF), and in-stream velocity variation was measured with particle-image velocimetry (PIV). Maximum in-stream temperatures were measured using the coherent anti-Stokes Raman scattering (CARS) technique. Lift-off heights of triple flames have minimum values during the increase of concentration gradient, and the propagation velocity of triple flames reaches its maximum at a critical concentration gradient. This is caused by three factors: velocity distribution upstream, flammable limit of premixed gas, and reaction of diffusion flame. The critical concentration gradient, which maximizes the propagation velocity is suggested as a new criterion of transition from a premixed flame to a triple flame.     

Dynamics of triple-flames in ignition of turbulent dual fuel mixture: A direct numerical simulation study

Pilot-ignited dual fuel combustion involves a complex transition between the pilot fuel autoignition and the premixed-like phase of combustion, which is challenging for experimental measurement and numerical modelling, and not sufficiently explored. To further understand the fundamentals of the dual fuel ignition processes, the transient ignition and subsequent flame development in a turbulent dimethyl ether (DME)/methane-air mixing layer under diesel engine-relevant conditions are studied by direct numerical simulations (DNS). Results indicate that combustion is initiated by a two-stage autoignition that involves both low-temperature and high-temperature chemistry. The first stage autoignition is initiated at the stoichiometric mixture, and then the ignition front propagates against the mixture fraction gradient into rich mixtures and eventually forms a diffusively-supported cool flame. The second stage ignition kernels are spatially distributed around the most reactive mixture fraction with a low scalar dissipation rate. Multiple triple flames are established and propagate along the stoichiometric mixture, which is proven to play an essential role in the flame developing process. The edge flames gradually get close to each other with their branches eventually connected. It is the leading lean premixed branch that initiates the steady propagating methane-air flame. The time required for the initiation of steady flame is substantially shorter than the autoignition delay time of the methane-air mixture under the same thermochemical condition. Temporal evolution of the displacement speed at the flame front is also investigated to clarify the propagation characteristics of the combustion waves. Cool flame and propagation of triple flames are also identified in this study, which are novel features of the pilot-ignited dual fuel combustion.     

重力对扩散射流火焰动态特性的影响

本文探讨重力对扩散射流火焰动态特性的影响规律。结果表明,火焰闪烁现象是一种浮力诱导不稳定性,在浮力消失或反向重力场中,不存在这种不稳定性现象,闪烁频率与燃料射流速度无直接关系,但涡的大小随燃料射流速度的增大而增大。存在触发火焰闪烁的临界高度,闪烁频率与重力成平方根关系式。反向重力情况下,也存在浮力稳定型平面火焰,它反映了浮力与火焰的耦合作用。     

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.     

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.     

Triple flame: Inherent asymmetries and pentasectional character

A two-dimensional triple-flame numerical model of a laminar combustion process reflects flame asymmetric structural features that other analytical models do not generate. It reveals the pentasectional character of the triple flame, composed of the central pure diffusion-flame branch and the fuel-rich and fuel-lean branches, each of which is divided into two sections: a near-stoichiometric section and a previously unreported near-flammability-limits section with combined diffusion and premixed character. Results include propagation velocity, fuel and oxidiser mass fractions, temperature and reaction rates. Realistic stoichiometric ratios and reaction orders match experimental planar flame characteristics. Constant density, a one-step reaction, and a mixture fraction gradient at the inlet as the simulation parameter are imposed. The upstream equivalence ratio or the upstream reactant mass fractions are linear or hyperbolic functions of the transverse coordinate. The use here of experimental kinetics data differs from previous analytical works and results in flame asymmetry and different flammability limits. Upstream mixture composition gradient affects propagation velocity, flame curvature, diffusion flame reaction rate, and flammability limits. Flammability limits extend beyond those of a planar flame due to transverse heat and mass diffusion causing the pentasectional character.     

The role of unequal diffusivities in ignition and extinction fronts in strained mixing layers

We have studied flame propagation in a strained mixing layer formed between a fuel stream and an oxidizer stream, which can have different initial temperatures. Allowing the Lewis numbers to deviate from unity, the problem is first formulated within the framework of a thermo-diffusive model and a single irreversible reaction. A compact formulation is then derived in the limit of large activation energy, and solved analytically for high values of the Damköhler number. Simple expressions describing the flame shape and its propagation velocity are obtained. In particular, it is found that the Lewis numbers affect the propagation of the triple flame in a way similar to that obtained in the studies of stretched premixed flames. For example, the flame curvature determined by the transverse enthalpy gradients in the frozen mixing layer leads to flame-front velocities which grow with decreasing values of the Lewis numbers.

The analytical results are complemented by a numerical study which focuses on preferential-diffusion effects on triple flames. The results cover, for different values of the fuel Lewis number, a wide range of values of the Damköhler number leading to propagation speeds which vary from positive values down to large negative values     

A comparison of computational and experimental lift-off heights of coflow laminar diffusion flames

As a sensitive marker of changes in flame structure, the number densities of excited-state CH (denoted CH^*), and excited-state OH (denoted OH^*) are imaged in coflow laminar diffusion flames. Measurements are made both in normal gravity and on the NASA KC-135 reduced-gravity aircraft. The spatial distribution of these radicals provides information about flame structure and lift-off heights that can be directly compared with computational predictions. Measurements and computations are compared over a range of buoyancy and fuel dilution levels. Results indicate that the lift-off heights and flame shapes predicted by the computations are in excellent agreement with measurement for both normal gravity (1g) and reduced gravity flames at low dilution levels. As the fuel mixture is increasingly diluted, however, the 1g lift-off heights become underpredicted. This trend continues until the computations predict stable flames at highly dilute fuel mixtures beyond the 1g experimental blow-off limit. To better understand this behavior, an analysis was performed, which indicates that the lift-off height is sensitive to the laminar flame speed of the corresponding premixed mixture at the flame edge. By varying the rates of two key “flame speed” controlling reactions, we were able to modify the predicted lift-off heights so as to be in closer agreement with the experiments. The results indicate that reaction sets that work well in low dilution systems may need to be modified to accommodate high dilution flames.     

Behaviour of a confined fire located in an unventilated zone

The behaviour of a fire in an enclosure is studied for a configuration where the fuel source is located in the upper hot unventilated zone trapped by a soffit. The experimental study, undertaken in a laboratory-scale compartment with a fuel source above the level of a soffit, included the determination of the parameters (ventilation factor, rate of fuel supply) controlling the combustion or leading to extinction. Measurements (PIV, thermocouples, gas sampling and analysis) were performed to propose a hypothesis on the structure of the flame (flame stabilisation mechanisms, premixed or diffusion types). Video photography is used to determine the area covered by the flames. This information is used as a criterion to identify the combustion regimes. The results show that the gaseous fuel is diluted in the combustion products (CP) in the upper layer and that a recirculatory motion is formed, driven by buoyancy forces, which enhances the mixing of fuel and CP. These then travel horizontally towards the vent along the interface between the lower fresh air and upper zones, and are premixed with the convected air in the enclosure, before entering the reaction zone and being burnt. The flame stabilises at the interface between the upper hot and lower ventilated layers in the compartment. The observed “ghosting flame” is stabilised by a triple flame if the flame speed of the premixed flame is higher than the natural convection velocity induced in the compartment. The flame stability is quantified by a criterion based on the area of the horizontal flame. It has been observed that the combustion is controlled by the available mass fuel flux at the reaction zone if the ventilation is sufficient. This information is essential for the modelling of the phenomena involved in fires with such an underventilated fuel source.     

Heat recirculation effects on flame propagation and flame structure in a mesoscale tube

Heat recirculation effects on flame propagation and flame structure are theoretically and experimentally examined in a mesoscale tube as the simplest model of heat-recirculating burners. Solutions for steady propagation are obtained using a one-dimensional two-temperature approximation. The results show that the low heat diffusivities of common solid materials permit significant heat recirculation through the wall only for a slowly-propagating condition, otherwise the flame behaves almost like a freely-propagating nonadiabatic flame. This limited heat recirculation sharply pinches and stretches two well-known branches of the freely-propagating nonadiabatic flame, resulting in the appearance of two slow-propagation branches. On the upper slow-propagation branch flames can reach superadiabatic temperatures and on the lower one, which is stretched from the classical unstable lower branch, flames can be stable. As the tube inner diameter decreases, another burning regime appears where flames are barely sustained by the heat recirculation. Further reduction of the tube inner diameter makes no flame exist. It is also revealed that a flame in a mesoscale tube has two length scales, i.e. the conventional flame thickness and a convective preheat zone thickness, and that the latter should be much larger than the former for significant heat recirculation. It is theoretically predicted that a heat-recirculating, even superadiabatic, flame with positive propagation velocity against the gas flow can exist in a mesoscale tube. It is also found that a flame transition from one branch to another in a given tube is well described by only one dimensionless parameter. Finally, these theoretical results show good qualitative agreements with experiments, especially for the transition behaviours.     

Stabilization mechanism of lifted flame edge in the near field of coflow jets for diluted methane

The stabilization mechanism of lifted flames in the near field of coflow jets has been investigated experimentally and numerically for methane fuel diluted with nitrogen. The lifted flames were observed only in the near field of coflow jets until blowout occurred in the normal gravity condition. To elucidate the stabilization mechanism for the stationary lifted flames of methane having the Schmidt number smaller than unity, the behavior of the flame in the buoyancy-free condition, and unsteady propagation characteristics after ignition were investigated numerically at various conditions of jet velocity. It has been found that buoyancy plays an important role for flame stabilization of lifted flames under normal gravity, such that the flame becomes attached to the nozzle in microgravity. The stabilization mechanism is found to be due to the variation of the propagation speed of the lifted flame edge with axial distance from the nozzle in the near field of the coflow as compared to the local flow velocity variation at the edge.     

The role of radiative losses in self-extinguishing and self-wrinkling flames

We have developed a general theory of non-adiabatic premixed flames that is valid for flames of arbitrary shape that fully accounts for the hydrodynamic and diffusive-thermal processes, and incorporates the effects of volumetric heat losses. The model is used to describe aspects of experimentally observed phenomena of self-extinguishing (SEFs) and self-wrinkling flames (SWFs), in which radiative heat losses play an important role. SEFs are spherical flames that propagate considerable distances in sub-limit conditions before suddenly extinguishing. Our results capture many aspects of this phenomenon including an explicit determination of flame size and propagation speed at quenching. SWFs are hydrodynamically unstable flames in which cells spontaneously appear on the flame surface once the flame reaches a critical size. Our results yield expressions of the critical flame size at the onset of wrinkling and expected cell size beyond the stability threshold. The various possible burning regimes are mapped out in parameter space.     

Role of hydroxyl production and heat release in the two-zone fuel-rich adiabatic dimethyl ether/air flames at atmospheric pressure

Fuel-rich laminar adiabatic flames of premixed dimethyl ether/air mixtures at a high initial temperature and atmospheric pressure have been studied by numerical simulation and sensitivity analysis. These flames, having two heat release zones, are of great interest as an unusual and little-studied subject. We have investigated the chemical processes occurring in the two zones and analysed the mechanism of heat release in the flame. It has been found that the key reactions that have a significant influence on the flame speed are those involving dimethyl ether and the products of its incomplete oxidation. Calculation of the heat release rate confirms the presence of two heat release zones in the flame. A comparison of the reactions making a major contribution to the heat release with those significantly affecting the flame speed indicates that the main factor determining the flame speed is the formation of hydroxyls, rather than heat release. Analysis of the flame speed sensitivity shows that in the case of a two-zone structure of the flame, its speed is mainly determined by the reactions taking place in the low-temperature zone. That is, the cool zone with a higher temperature gradient is the leading one.     

Vortex dynamics in different combustion regions of a cavity-based scramjet

A hybrid RANS/LES study of a cavity-based scramjet was performed and reasonable agreements were found between simulation results and experimental measurements. In the current case, the flame was stabilized by the subsonic cavity shear layer and propagated downstream into the supersonic flow. The vortex dynamic in the flow, mixing, and combustion regions was comparatively investigated. The averaged vorticity in the combustion regions was lower by 55% compared to the mixing region, primarily due to dilatation as a result of the heat release. Furthermore, the combustion zone was decomposed into four regions based on premixed/diffusion flame and subsonic/supersonic combustion. Then the vorticity and its transport in the four regions were compared. The averaged vorticity in the premixed combustion regions was only slightly larger than that in the diffusion combustion regions. However, the averaged heat release rate was nearly 3 times larger in the premixed regions, leading to higher contributions of dilatation and baroclinic torque in the premixed regions, with an overall weak positive impact on the vorticity generation. In the subsonic combustion regions, the vorticity was three times larger than that in the supersonic combustion regions, despite similar heat release rates on average. It could be explained by the relatively large magnitude of dilatation and baroclinic torque in the supersonic flow. Vortex stretching and dilatation were comparable in the supersonic flame but the former became two times larger than the latter in the subsonic flame. Moreover, the baroclinic torque had larger contributions than diffusion in the supersonic flame whereas the opposite trend was found in the subsonic flame. The results highlight that the subsonic combustion regions in the cavity shear layer and near the walls significantly contribute to the vortex dynamics and mixing process, in addition to flame stabilization.     

Analytical predictions of shapes of laminar diffusion flames in microgravity and earth gravity

Flame shape is an important observed characteristic of flames that can be used to scale flame properties such as heat release rates and radiation. Flame shape is affected by fuel type, oxygen levels in the oxidiser, inverse burning and gravity. The objective of this study is to understand the effect of high oxygen concentrations, inverse burning, and gravity on the predictions of flame shapes. Flame shapes are obtained from recent analytical models and compared with experimental data for a number of inverse and normal ethane flame configurations with varying oxygen concentrations in the oxidiser and under earth gravity and microgravity conditions. The Roper flame shape model was extended to predict the complete flame shapes of laminar gas jet normal and inverse diffusion flames on round burners. The Spalding model was extended to inverse diffusion flames. The results show that the extended Roper model results in reasonable predictions for all microgravity and earth gravity flames except for enhanced oxygen normal diffusion flames under earth gravity conditions. The results also show trends towards cooler flames in microgravity that are in line with past experimental observations. Some key characteristics of the predicted flame shapes and parameters needed to describe the flame shape using the extended Roper model are discussed.     

Effect of non-zero relative velocity on the flame speed of two-phase laminar flames

A numerical study of one-dimensional n-heptane/air spray flames is presented. The objective is to evaluate the flame propagation speed in the case where droplets evaporate inside the reaction zone with possibly non-zero relative velocity. A Direct Numerical Simulation approach for the gaseous phase is coupled to a discrete particle Lagrangian formalism for the dispersed phase. A global two-step n-heptane/air chemical mechanism is used. The effects of initial droplet diameter, overall equivalence ratio, liquid loading and relative velocity between gaseous and liquid phases on the laminar spray flame speed and structure are studied. For lean premixed cases, it is found that the laminar flame speed decreases with increasing initial droplet diameter and relative velocity. On the contrary, rich premixed cases show a range of diameters for which the flame speed is enhanced compared to the corresponding purely gaseous flame. Finally, spray flames controlled by evaporation always have lower flame speeds. To highlight the controlling parameters of spray flame speed, approximate analytical expressions are proposed, which give the correct trends of the spray flame propagation speed behavior for both lean and rich mixtures.