Theoretical and experimental studies on mesoscale flame propagation and extinction

Mesoscale flame propagation and extinction of premixed flames in channels are investigated theoretically and experimentally. Emphasis is placed on the effect of wall heat loss and the wall–flame interaction via heat recirculation. At first, an analytical solution of flame speed in mesoscale channels is obtained. The results showed that channel width, flow velocity, and wall thermal properties have dramatic effects on the flame propagation and lead to multiple flame regimes and extinction limits. With the decrease in channel width, there exist two distinct flame regimes, a fast burning regime and a slow burning regime. The existence of the new flame regime and its extended flammability limit render the classical quenching diameter inapplicable. Furthermore, the results showed that at optimum conditions of flow velocity and wall thermal properties, mesoscale flames can propagate faster than the adiabatic flame. Second, numerical simulation with detailed chemistry demonstrated the existence of multiple flame regimes. The results also showed that there is a non-linear dependence of the flame speed on equivalence ratio. Moreover, it is shown that the Nusselt number has a significant impact on this non-linear dependence. Finally, the non-linear dependence of flame speed on equivalence ratio for both flame regimes is measured using a C₃H₈–air mixture. The results are in good agreement with the theory and numerical simulation.     

Transition in the propagation mechanism during flame acceleration in porous media

The combustion of stoichiometric hydrogen-air at various initial pressures was investigated in a 7.62 cm square cross-section channel filled with 1.27 cm diameter beads. The flame time-of-arrival and pressure time history along the channel were obtained by ionization probes and piezoelectric pressure transducers. Flame acceleration was found to be very rapid, e.g. at an initial pressure of 45 kPa the flame achieves a velocity of over 600 m/s in roughly 0.3 m. It was determined that at this high speed a well defined planar shock wave precedes a thick reaction zone. It was also shown that there is a transition in the flame propagation mechanism, similar to that observed in an obstacle laden channel [G. Ciccarelli and C. Johansen, The role of shock-flame interactions on flame acceleration in an obstacle laden channel, Proc. 22nd International Colloquium on the Dynamics of Explosions and Reactive Systems, Minsk, 2009]. By varying the initial pressure of the mixture, changes in the axial location of the transition between combustion propagation regimes was also observed. A soot foil technique was used to identify the transition in the propagation mechanism, as well as to provide information concerning the local flow field around the beads and the overall average flow direction.     

Flame propagation in gaseous nitrous oxide

The characteristics of the propagation of a nitrous oxide decomposition flame in a tube with an internal diameter of 70 mm were measured. It was demonstrated that the pattern of flame propagation and the extent of burnout are determined by the convective motion of the flame kernel because of a very slow burning of nitrous oxide. The laminar flame speed estimated from pressure oscillograms and calculated using thermal theory of flame propagation was found to be ～1 cm/s. The critical diameter of flame quenching in channels were measured to decrease from 10 to 4 mm as the pressure was increased from 15 to 20 atm. Because of the possibility of reignition of the fresh mixture behind the flame arrester by the outflowing combustion products, the channel should be significantly longer than 200 mm.     

Hydrodynamic instability of premixed flame propagating in narrow planar channel in the presence of gas flow

The paper analyses the hydrodynamic instability of a flame propagating in the space between two parallel plates in the presence of gas flow. The linear analysis was performed in the framework of a two-dimensional model that describes the averaged gas flow in the space between the plates and the perturbations development of two-dimensional combustion wave. The model includes the parametric dependences of the flame front propagation velocity on its local curvature and on the combustible gas velocity averaged along the height of the channel. It is assumed that the viscous gas flow changes the surface area of the flame front and thereby affects the propagation velocity of the two-dimensional combustion wave. In the absence of the influence of the channel walls on the gas flow, the model transforms into the Darrieus–Landau model of flame hydrodynamic instability. The dependences of the instability growth rate on the wave vector of disturbances, the velocity of the unperturbed gas flow, the viscous friction coefficients and other parameters of the problem are obtained. It is shown that the viscous gas flow in the channel can lead, in some cases, to a significant increase in instability compared with a flame propagating in free space. In particular, the instability increment depends on the direction of the gas flow with respect direction of the flame propagation. In the case when the gas flow moves in the opposite direction to the direction of the flame propagation, the pulsating instability can appear.     

Combustion of n-butanol in a flat minichannel

Combustion of n-butanol was studied experimentally at the two-phase flow in a rectangular channel of 4-mm height. Characteristics of the two-phase flow in channels with close configuration were analyzed, and the conditions of flame spread in the channel were determined. The dependencies of flame propagation velocity on flow velocity of the gas mixture with different oxygen contents were measured.     

Effect of ammonia/oxygen/nitrogen equivalence ratio on spherical turbulent flame propagation of pulverized coal/ammonia co-combustion

Because ammonia is one of the most promising candidates for energy carrier in the future, various applications of ammonia as a fuel are currently considered. One medium for utilizing ammonia is by introducing it to coal-fired boilers. To the best of our knowledge, this paper is the first to report the fundamental mechanism of the flame propagation phenomenon for pulverized coal/ammonia co-combustion. The effects of the equivalence ratio of the ammonia-oxidizer mixture on the flame propagation velocity of pulverized coal/ammonia co-combustion in turbulent fields were clarified by the experiments employing a unique fan-stirred constant volume chamber. The flame propagation velocities of pulverized coal/ammonia co-combustion, pure ammonia combustion, and pure pulverized coal combustion were compared. As expected, the flame propagation velocity of pulverized coal/ammonia was higher than that of the pure pulverized coal combustion for all conditions. However, the comparison of the flame propagation velocities of pulverized coal/ammonia co-combustion and that of the pure ammonia combustion, revealed that whether the flame propagation of the pulverized coal/ammonia was higher than that of the pure ammonia combustion was dependent on the equivalence ratio of the ammonia-oxidizer. This unique feature was explained by a mechanism including three competing effects proposed by the authors. In the ammonia lean condition, the positive effects, which are the strong radiation from the luminous flame and the increment of local equivalence ratio by the addition of volatile matter, are larger than the negative effect, which is the heat absorption by coal particles in preheat zone. In the ammonia rich condition, the effect of an increment of the local equivalence ratio by the addition of volatile matter turns into a negative effect. Consequently, the negative effects overcome the positive effect in the ammonia rich condition resulting in a lower flame propagation velocity of pulverized coal/ammonia co-combustion.     

Formation of dynamic spatial flame structures for gas burning in microchannels with temperature gradients on walls

The thermal-diffusive model was applied to the problem of flame propagation in a microchannel with controlled temperature distribution in the walls; this demonstrated the possibility of formation of oscillating or rotating spatial flame structures, which were described previously in experimental works on microcombustion. Two cases were considered: combustion in a rectangular channel and in the clearance between two disks with radial feeding of premixture. In both cases, the typical across size of the channel was lower than the critical diameter determined with respect to the ambient temperature. The gas flow was assigned and described by the Poiseuille-flow velocity profile. Formation of oscillating flame in a rectangular channel and rotating patterns in a radial channel was observed for a certain range of gas flow rate. At low flow rates beyond this range, repetitive ignition/extinction of flame took place; at high flow rate we observed a steady flame mode. Formation of these special flame structures is related to heat transfer between gas and hot walls of the channel, as well as to velocity maldistribution in the microchannel.     

Characteristic regimes of premixed gas combustion in high-porosity micro-fibrous porous media

Dynamical behaviour of the premixed flame propagating in the inert high-porosity micro-fibrous porous media has been studied numerically. Effects of mixture filtration velocity, equivalence ratio and burner transverse size on the flame structure have been investigated and the regions of existence of different combustion regimes have been determined. It was found that the influence of the hydrodynamic instability on the flame dynamics is significant in the case of the moderate and high filtration velocities and this effect is negligible at the low velocities. At the moderate filtration velocities the effect of hydrodynamic instability manifests in the flame front deformation and in particular in the flame inclination. It was found that the flame can be stabilized within the whole interval of the filtration gas velocity, whereas in the ordinary porous media the standing wave is settled only at fixed value of gas filtration velocity. This finding is in line with recent experimental results on combustion in micro-fibrous porous media (Yang et al., Combust. Sci. Tech. 181 (2009), 1–16). Possible physical interpretation of the flame anchoring effect may be given on the base of present numerical analysis. At the high filtration velocities the hydrodynamic instability manifests itself in periodical appearance of the moving wrinkles on the flame front surface which forms non stationary high temperature trailing spots behind the leading part of the flame front. Such dynamics may be associated with splitting wave structures which were revealed in previous experiments (Yang et al., Combust. Sci. Tech. 181 (2009), 1–16).     

Flame behavior in heated porous sand bed

This paper reports experimental and numerical investigations on the combustion characteristics of a lean methane–air mixture in a heated porous sand bed. The porous bed consisted of sand (SiO₂) particles with a mean particle diameter of 0.56 mm. The horizontally placed quartz tube was heated externally to initiate the combustion reaction in the porous bed combustor. The stabilized flame location curve as a function of averaged mixture velocity was obtained for various equivalence ratios. Contrary to the earlier finding of a C-shape flame stabilization behavior, a new S-shape behavior was observed in the present study. This can be divided into three regimes: high, moderate, and low velocity regimes. In the low velocity regime, flame with very weak luminosity was confirmed and the stabilized flame location moved downstream with the increase of the mixture velocity. For the moderate velocity regime, a stable flat flame was observed and the flame location moved upstream with the increase of the mixture velocity. An oscillatory flame behavior was observed in the high velocity regime. In this oscillatory mode, the flame front oscillated with a characteristic time period of the order on 1 h and increased with the increase of the mixture velocity. In order to further understand these experimental results, one-dimensional computational studies with detailed chemistry and heat transfer mechanisms were carried out. The computational results were in good agreement with experimental observations. The computations showed that solid-to-solid radiation played a significant role in the flame stabilized location. From the examination of the flame structure, it was found that the flame behavior in the low velocity regime was similar to that of the flameless combustion mode.     

Turbulent flame propagation limits of ammonia/methane/air premixed mixture in a constant volume vessel

Ammonia is one of promising energy carriers that can be directly used as carbon-neutral fuel for combustion applications. However, because of the low-burning velocity of ammonia, it is challenging to introduce ammonia to practical combustors those are designed for general hydrocarbon fuels. One of ways to enhance the combustibility of ammonia is by mixing it with other hydrocarbon fuels, such as methane, with a burning velocity is much higher than the burning velocity of ammonia. In this study, we conducted flame propagation experiments of ammonia/methane/air using a fan-stirred constant volume vessel to clarify the effect of methane addition to ammonia on the turbulent flame propagation limit. From experimental results, we constructed the flame propagation maps and clarified the flame propagation limits. The results show that the flame propagation limits were extended with an increase in mixing a fraction of methane to ammonia. Additionally, ammonia/methane/air mixtures with the equivalence ration of 0.9 can propagate at the highest turbulent intensity, even though the peak of the laminar burning velocity is the fuel-rich side because of the diffusional-thermal instability of the flame surface. Furthermore, the Markstein number of the mixture obtained in this research successfully expressed the strength of the diffusional-thermal instability effect on the flame propagation capability. The turbulence Karlovitz number at the flame propagation limit monotonically increases with the decreasing Markstein number.     

基于波传播算法的火焰不稳定性

基于波传播算法构造了多组分反应流的数值格式,利用CH4空气基元反应动力学模型,并采用分离算法,对CH4空气混合物中,入射激波与火焰的相互作用,以及反射激波与火焰的二次作用过程进行了数值模拟.根据计算结果,讨论了激波诱导火焰失稳的发展过程及其特点.结果表明,Helmholtz不稳定、RichtmyerMeshkov不稳定以及反应放热速率对火焰失稳过程有重要影响.计算结果与实验结果进行了比较,对数值方法的有效性进行了验证.     

Characterization of symmetric to non-symmetric flamefront transition in slender microchannels

The purpose of the present work is to analyze propagating two-dimensional flames confined in slender semi open channels, where the combustion process takes place towards the closed end. The study focuses on the calculation of the growth rate of the transition from symmetric to non-symmetric flames propagation identified by Jiménez et al. [1].The combustion cell is initially filled with a stoichiometric mixture of fuel and air at standard conditions. Ignition is induced close to the open end of the channel under planar and gaussian profiles in temperature and species mass fractions which activate a sustained combustion process. The gases inside the chamber, initially stagnant, are accelerated due to the heat generated in the chemical reactions, leading to the development of lateral boundary layers, so that the hot gases exit the channel following a Poiseuille velocity profile. This transverse flow velocity differences are accommodated by means of a symmetric tulip shape formed after a short initial transient.Acoustic waves generated in the ignition process, keep travelling along the channel, bouncing at the walls and interacting with the flame during all the combustion process. Additionally, the flame structure, curved by Darrieus-Landau instability, interacts with the pressure waves triggering small amplitude oscillations (primary oscillation mode), which under certain conditions can transition to higher amplitude oscillations (secondary mode).This transition is observed to be highly dependent not only on the cell geometry, but also in the initial conditions generated by the ignition procedure.The aim of this work is to improve the understanding of this process, complementing the work of Jiménez et al. [1], and to characterize the effect of the channel width in this transition.     

Characteristic analysis of low-velocity gas filtration combustion in an inert packed bed

This paper investigates the low-velocity filtration combustion of lean methane–air mixtures occurring in inert packed beds by using a modified one-temperature model, considering the axial thermal diffusion owing to the convective gas–solid heat transfer. Based on the scaling analysis of various transport terms in different conservation equations, a high-activation energy asymptotic method is applied in the flame zone and results in a set of powerful analytical solutions for combustion macrocharacteristics under the fully developed conditions. These are then combined with the eigenvalue method of the modified one-temperature model in the whole flow region to study the flame behaviour analytically and numerically. Our results have shown that the combustion wave velocity is a key characteristic parameter in the filtration combustion process. Compared with other existing theoretical results, the present analytical solutions demonstrate the intricate relationships among the combustion wave velocity, the flame speed, the peak flame temperature and the effects of the variable thermo-physical properties, and show better prediction performance for the combustion wave velocity, the flame speed and the peak flame temperature. Excellent agreements with experimental results have been observed, especially for very lean filtration combustion with stream-wise propagating combustion fronts.     

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.     

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.     

Influence of water mist on propagation and suppression of laminar premixed flame

The combustion of premixed gas mixtures containing micro droplets of water was studied using one-dimensional approximation. The dependencies of the burning velocity and flammability limits on the initial conditions and on the properties of liquid droplets were analyzed. Effects of droplet size and concentration of added liquid were studied. It was demonstrated that the droplets with smaller diameters are more effective in reducing the flame velocity. For droplets vaporizing in the reaction zone, the burning velocity is independent of droplet size, and it depends only on the concentration of added liquid. With further increase of the droplet diameter the droplets are passing through the reaction zone with completion of vaporization in the combustion products. It was demonstrated that for droplets above a certain size there are two stable stationary modes of flame propagation with transition of hysteresis type. The critical conditions of the transition are due to the appearance of the temperature maximum at the flame front and the temperature gradient with heat losses from the reaction zone to the products, as a result of droplet vaporization passing through the reaction zone. The critical conditions are similar to the critical conditions of the classical flammability limits of flame with the thermal mechanism of flame propagation. The maximum decrease in the burning velocity and decrease in the combustion temperature at the critical turning point corresponds to predictions of the classical theories of flammability limits of Zel'dovich and Spalding. The stability analysis of stationary modes of flame propagation in the presence of water mist showed the lack of oscillatory processes in the frames of the assumed model.     

Geometrical study of spark ignition in two dimensions

Spark ignition, as the first step during the combustion in Otto engines, has a profound impact on the further development of the flame kernel. Over the last ten years growing concern for environment protection, including low emission of pollutants has increased the interest in the numerical simulation of ignition phenomena to guarantee successful flame kernel development even for lean mixtures.

However, the process of spark ignition in a combustible mixture is not yet fully understood. The use of detailed reaction mechanisms, combined with electrodynamical modelling of the spark, is necessary to optimize ignition of lean mixtures.

This work presents simulations of the coupling of flow, chemical reactions and transport with discharge processes in order to investigate the development of a stable flame kernel initiated by an electrical spark. A two-dimensional code to simulate the early stages of flame kernel formation, shortly after the breakdown discharge, has been developed. The model includes Joule heating. The spark plasma channel formed as a consequence of the breakdown is incorporated into the initial conditions. The computations include the initial phase (1–5 µs), which is governed by pressure wave formation, but also the transition to flame propagation. A thorough study of the influence of the electrodes' geometry, i.e. shape and size, and gap width, has been performed for air and a lean H₂–air mixture. Also a detailed methane-air mechanism was chosen as another example including combustion.

Due to the fast expansion of the plasma channel, together with the geometrical complexity of the electrodes, a complicated flow field develops after the emission of a pressure wave by the expanding channel. Special numerical methods, including artificial viscosity, are required to resolve the complicated flow field during these first 1–5 µs. The heat release through chemical reactions and transport processes is almost negligible during this short phase. The second phase, i.e. the development of a propagating flame and the flame kernel expansion, can last up to several milliseconds and is dominated by diffusive processes and chemical reactions. It has been found that the geometry greatly influences the developing flame kernel and the flow field. As soon as chemical reactions begin to contribute significantly to the heat release, the effect becomes smaller.     

The role of duct thickness on the quenching process of premixed flame propagation

The propagation of premixed laminar flame in ducts of circular cross-section considering a thermal-diffusive model is investigated numerically. Heat losses by conduction to the channels walls are taken into account using the thermally thin wall regime. The effects and the relationship between thickness and diameter of the tube with the flame speed propagation are studied and the quenching condition is obtained as a function of the heat-loss parameter. The mathematical model employs the axisymmetric energy and species equations. The calculations are based on a two-step chemistry, with an Arrhenius, energetically neutral, radical production reaction followed by an exothermic radical recombination reaction. For large values of the heat-loss parameter, the wall temperature is close to the free stream temperature and all the heat losses through the wall are convected away. No heat feedback occurs. On the other hand, for small values of the heat-loss parameter, a feedback mechanism occurs by transferring heat from the burned gas to the fresh mixture along the tube wall. For values of the heat-loss parameter of order unity, the heat feedback mechanism is able to sustain the flame propagation and the quenching condition disappears, producing an almost planar flame front as the propagation velocity reduces. For this two-step reaction mechanism, the radical species behaviour at the duct walls seems to have negligible effect on the quenching process.     

A comparative study of H2-air premixed flame in micro combustors with different physical and boundary conditions

A numerical study of H₂-air premixed combustion in the micro channels with a detailed chemical reaction mechanism is performed by solving the two-dimensional fully elliptic governing equations of continuity, momentum, energy and species, coupled with the energy equation in the solid wall. A reference case is defined as the combustion in a cylindrical tube with 0.8 mm inner diameter and 8 mm length with a non-slip wall and a uniform velocity profile at the inlet plane. Different physical and boundary conditions have been applied in order to investigate their respective effects on the flame temperature. The conditions studied in the current paper include the combustor size and geometry, inlet velocity profile, axial heat conduction in the solid wall and slip-wall and temperature jump at the gas–solid interface. It is noted that effects of Knudsen number (slip-wall and temperature jump) on the thermal and fluid field are not very significant in a d = 0.4 mm micro combustor. Furthermore, the qualitative effects of Knudsen number on the flame temperature are analysed. The results of this paper indicate that these various boundary and physical conditions have effects on the flame temperature to different extent and should be carefully monitored when applied for different applications.     

A DNS study on the flame structures and flame stabilization mechanism of a laboratory-scale lean premixed jet flame in crossflow

In the present work, direct numerical simulation (DNS) of a laboratory-scale lean premixed reacting jet flame in crossflow was performed to understand the flame structures and the flame stabilization mechanism. In the DNS, an ethylene-air jet with an equivalence ratio of 0.6 was injected into a hot vitiated crossflow. The jet Reynolds number reaches 6161. The DNS results were compared with those of the experiment with a good agreement. It was found that the windward and leeward branches of the flame show significantly different behaviors. The windward flame branch, appearing lifted and discontinuous, is located in the shear layer regions with high temperature, low vorticity and low scalar dissipation rate. The location of the peak heat release rate shifts to a higher mixture fraction with increasing distance from the jet exit. The leeward branch of the flame anchors in the shear layer near the jet exit. The recirculation zone in the wake of the jet facilitates the stabilization of the leeward flame. The chemical explosive mode analysis (CEMA) and species budget analysis were employed to characterize the local combustion mode. Auto-ignition plays a key role in the stabilization of the windward flame where a large range of extinction is also found due to the high strain rate. In contrast, premixed flame propagation is dominant on the leeward side.