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).     

Numerical study of sporadic combustion waves in straight channels of different diameters

The dynamics of flames propagating in straight channels filled with a stationary low-Lewis-number premixed gas mixture is studied numerically. A method for determining the propagation velocity of a sporadic combustion wave consisted of separate flame spots is proposed. Dependencies of the sporadic combustion wave propagation velocity, the residual fuel concentration and the number of flame spots on the channel size and the value of radiation heat losses are obtained. Analysis of numerical results show that for the channels of diameter exceeding some value the number of separate cup-like fragments constituting sporadic combustion wave is proportional to the channel cross-sectional area. At smaller diameters, the number of flame spots changes insignificantly and is one or two. It is shown that one of the universal characteristics of the sporadic combustion wave depending only on mixture properties but independent on system geometry is the area necessary to accommodate one reacting spot. Flame velocity which is another fundamental combustion characteristic is found to be almost independent on channel size starting from some critical diameter. This diameter, however, depends on mixture properties or radiative heat loss intensity and corresponds to the sporadic flame containing from several to ten reacting spots. Thus, the main properties of sporadic combustion waves in wide channels can be determined by numerical modeling of the flame propagation in the relatively narrow channels in which the flame consists of 1–10 cup-like fragments.     

An Analytical Study of Hydrodynamic Instability in the Flame. 1. Viscous Gas in the Flame Zone

The problem of the hydrodynamic stability of slow combustion is analytically solved with consideration given to the viscosity of the gas in the flame zone, the temperature dependence of the viscosity, and the dependences of the flame speed on the front curvature according to the Markstein model and on the pressure. The viscous forces in the flame zone alone cannot ensure the stability of the flame at any values of the Reynolds number. These forces act only as amplifiers of the stabilizing factor according to the Markstein model or in the case of a negative dependence of the flame velocity on the pressure. This property of internal friction forces is the more pronounced, the stronger the viscosity increases with the temperature. Thermal expansion is not only a destabilizing factor, leading to an increase in viscosity and other transport coefficients, but also produces a stabilizing effect.     

Flame balls dynamics in divergent channel

A three-dimensional reaction-diffusion model for lean low-Lewis-number premixed flames with radiative heat losses propagating in divergent channel is studied numerically. Effects of inlet gas velocity and heat-loss intensity on flame structure at low Lewis numbers are investigated. It is found that continuous flame front exists at small heat losses and the separate flame balls settled within restricted domain inside the divergent channel at large heat losses. It is shown that the time averaged flame balls coordinate may be considered as important characteristic analogous to coordinate of continuous flame stabilized in divergent channel.     

Fingering instability in forward smolder combustion

We use numerical strategies to examine the linear and nonlinear stability of forward smolder waves in the framework of a simplified thermal–diffusive model, with the hydrodynamic effects completely filtered out. The configuration consists of a horizontal thin solid fuel, over which air blows in the same direction as the smolder front propagation. It is found that, in the absence of convective heat losses, the whole one-dimensional adiabatic solution branch is linearly stable; in contrast, when the convective heat loss effect is taken into account, fingering instability emerges provided the incoming air flow rate is within a narrow range near the one-dimensional extinction limit, a manifestation that is reminiscent of the familiar cellular instability occurring in the context of low-Lewis-number diffusion flames. Accordingly, the fingering instability herein identified in forward smolder combustion is purely thermal–diffusive in nature. Furthermore, a heuristic analysis by drawing an analogy with premixed flame suggests that the occurrence of such fingering instability is the joint consequence of the Lewis number effects and convective heat losses. It is proposed that a Hele–Shaw-type combustion channel may be adopted to experimentally reveal the fingering patterns predicted by current numerical simulations.     

Effect of thermo-hydrodynamic instability on structure and characteristics of filtration combustion wave of solid fuel

Numerical modelling of flame front stability for the inverse wave (with trailing combustion front) of filtration combustion of solid fuel is performed. The problem is treated in terms of dimensionless variables and parameters. It is found that propagation of a plane combustion front becomes unstable under certain conditions. In this case the front spontaneously inclines. The thermo-hydrodynamic mechanism is supposed to be responsible for instability developing. Anisotropic effective mass diffusivity (dispersion) is also taken into account. It turns out that anisotropic diffusivity affects structure and conversion distribution of the inclined combustion front. It is shown that the key parameters determining stability of combustion wave are dimensionless gas flow rate and width of reactor. The range of these parameters corresponding to the stable plane front is determined. It is shown that stability occurs either for small reactor widths (dimensionless values <1), or low gas flow rate (below 0.2). The optimised values of considered dimensionless parameters for maximal productivity are determined.     

Effects of imperfect premixing coupled with hydrodynamic instability on flame propagation

The problem of flame propagation in imperfectly premixed mixtures—mixtures of reactants with variable composition—is considered in this numerical study. We carry out two-dimensional direct numerical simulations of a flame propagating in a globally lean fuel-oxidizer mixture with imposed velocity and composition fluctuations of various intensities. The configuration adopted is that of a flame front interacting with spatially evolving fluctuations, and the characteristic scales of the domain and of the fluctuations imposed are significantly larger than the characteristic thickness of the flame, to account for important flame dynamics such as the hydrodynamic instability. One-step chemistry and Fick’s diffusion law are considered, along with unity Lewis number assumption for all the species. It is observed, in agreement with previous results, that relatively weak fluctuations in composition alone may lead to a large increase in flame length and burning rate. The hydrodynamic instability caused by gas expansion, catalyzed by the composition fluctuations interacting with the flame, is found to be responsible for the flame length enhancement. It is observed as well that the relative importance of this effect diminishes as the velocity fluctuations present become more intense, and that composition fluctuations have a small impact on flame length for these cases. It is additionally found that, with increasing intensity of composition fluctuations, there is eventually a reduction of burning rate per unit length of flame which leads, consequently, to a weak reduction of overall burning rate for the largest velocity fluctuation intensities covered by this study.     

A numerical modeling of the acceleration of a flame by an additional energy input ahead of its front

The acceleration of a flame after an additional energy input ahead of its front was simulated using numerical methods. The combustion of a hydrogen-air mixture in a semiopen channel was considered. The calculations were performed within the framework of a two-dimensional hydrodynamic model of premixed flames, with consideration given to heat transfer, multicomponent diffusion, and chemical kinetics. It was demonstrated that, when the interaction of the flame front with the near-wall boundary layer is taken into account, even a moderate energy input could substantially promote the development of the Landau-Darrieus instability and, possibly, deflagration to detonation transition.     

The study of geometric effects on the explosion front propagation in a horizontal channel with a layer of spherical beads

Experiments were performed in a horizontal channel partially filled with a layer of 12.7 mm ceramic-oxide beads filled with a nitrogen-diluted stoichiometric methane–oxygen mixture, i.e., CH₄ + 2(O₂ + 2/3N₂). Ionization probes and pressure transducers were used to track the explosion front velocity in the 1.22 m long, 76 mm wide and 152 mm high horizontal channel. Schlieren photography and smoked foil techniques are used to gain insight into the explosion front structure. The explosion propagation phenomenon was characterized by the combustion in the bead layer and the unobstructed gap above. It was determined that for a fixed gap height the bead layer thickness had very little effect on the explosion propagation phenomenon. However, for a fixed bead layer height the explosion propagation was strongly influenced by the gap height. The combustion products vented from the bead layer behind the flame propagating in the gap affects the structure of the shock-flame front in the gap and the maximum flame velocity achieved. The coupling between the vented products and the flame velocity in the gap was strongly influenced by the gap height. The gap height also affects the structure of the detonation wave propagating in the gap following DDT that always occurred in the gap. The DDT run-up distance was found to increase with increasing gap height and inversely with initial pressure.     

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.     

Consequence analysis of blast wave from accidental gas explosions

Recently, consequence analyses of accidental gas explosions are often carried out to assess the risk of chemical plants, hazardous-materials sites and new energy systems. In these consequence analyses, it is indispensable to adequately predict the blast-wave (pressure-wave) intensity from gas deflagrations. Some prediction models already exist; however, most of them are based on the theory for explosives and adjusting parameters are needed for evaluating gas deflagrations. In this study, new prediction methods for gas deflagrations were developed. From theoretical analysis of blast-wave generation by a gas deflagration, an evaluation equation of the blast-wave intensity was derived. As the scale of gas deflagration becomes larger, flame front instability (especially hydrodynamic instability) would be more effective and the flame propagating velocity starts to be accelerated. Therefore, the equation was modified considering the effect of flame instability. The evaluations by this modified equation agreed well with the results of large scale experiments. By this analysis, it was found that not only total energy release but also combustion reaction rate has to be introduced into the prediction of gas deflagrations. Using this concept, a modified scale model to predict the blast-wave intensity was developed by improving the previous scale model introducing the term of combustion reaction rate as burning velocity. Furthermore, scale analysis was performed to develop the new scaling law. The universal relationship between scaled distance and overpressure has been realized by this new scaling law for gas deflagrations. In summary, these results provide new methods for accurate prediction of the blast-wave intensity from gas deflagrations.     

DNS analysis of boundary layer flashback in turbulent flow with wall-normal pressure gradient

The presence of swirl in combustion systems produces a marked change in their boundary layer flashback behaviour. Two aspects of swirling flow are investigated in this study: the effect of the swirl-generated wall-normal pressure gradient, and the effect of misalignment between the mean flow direction and the direction of flame propagation. The analysis employs Direct Numerical Simulation (DNS) of fuel-lean premixed hydrogen-air flames in turbulent planar channel flow with friction Reynolds number of 180. The effect of swirl on the flashback process is investigated by imposing a wall-normal pressure gradient profile. Analysis of the DNS data shows how the resulting differences in flow field and flame topology contribute to the differences in the overall flashback speed. Misalignment of the flow and propagation directions leads to asymmetry in the flame shape statistics as streaks of high velocity fluid in the boundary layer cleave into the flame front at an angle, yielding an increase in flame surface density away from the wall. Swirl has a stabilising effect on the turbulent flame front during flashback along the centre-body of a swirling annular flow due to the density stratification across the flame front, and produces a reduction in turbulent consumption speed. However the swirl also sets up a hydrostatic pressure difference that drives the flame forward, and the net effect is that the flashback speed is increased. The dominance of hydrostatic effects motivates development of relatively simple modelling for the effect of swirl on flashback speed. A model accounting for the inviscid momentum balance and for confinement effects is presented which adequately describes the effect of swirl on flashback speed observed in previous experimental studies.     

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.     

聚心火焰在共振腔作用下引发爆轰的数值研究

 基于带化学反应的二维轴对称Euler方程，利用带有Superbee限制函数的波传播算法，对共振腔中的氢气-空气预混气的聚心燃烧进行了数值模拟，讨论了共振腔不同抛物面对起爆的影响。数值结果表明，在开始阶段，燃烧诱导的激波在轴心、火焰和固壁的反射，使火焰失稳，随后共振腔中的抛物壁面上产生一定频率和强度的反射激波，不断穿越火焰，使火焰进一步失稳，加剧了燃烧速度，最终导致爆轰的形成。同时，火焰在与激波的作用过程中，形状扭曲变形，呈封闭端小敞口端大的扁平头部蘑菇云。共振腔抛物面的不同形状引起激波聚焦位置的变化，会影响激波和火焰的相互作用，使起爆提前或推迟，甚至不起爆。     

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.     

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.     

Temporal evolution of flame stretch due to turbulence and the hydrodynamic instability

The temporal evolution of the strain rate on a turbulent premixed flame was measured experimentally using cinema-stereoscopic particle image velocimetry. Turbulence strains a flame due to velocity gradients associated both directly with the turbulence and those caused by the hydrodynamic instability, which are initiated by the turbulence. The development of flame wrinkles caused by both of these mechanisms was observed. Wrinkles generated by the turbulence formed around vortical structures, which passed through the flame and were attenuated. After the turbulent structures had passed, the hydrodynamic instability flow pattern developed and caused additional strain. The hydrodynamic instability also caused the growth of small flame front perturbations into large wrinkles. In the moderately turbulent flame investigated, it was found that the evolution of the strain rate caused by turbulence–flame interactions followed a common pattern involving three temporal regimes. In the first, the turbulence exerted extensive (positive) strain on the flame, creating a wrinkle that had negative curvature (concave towards the reactants). This was followed by a transition period, leading into the third regime in which the flow pattern and strain rate were dominated by the hydrodynamic instability mechanism. It was also found that the magnitudes of the strain rate in the first and third regimes were similar. Hence, the hydrodynamic instability mechanism caused significant strain on a flame and should be included in turbulent combustion models.     

Impact of the gravity field on stability of premixed flames propagating between two closely spaced parallel plates

The stability of a planar flame front propagating between two parallel adiabatic plates inclined at an arbitrary angle is investigated in the frame of narrow-channel approximation. It is demonstrated that buoyancy forces can suppress the hydrodynamic (Darrieus–Landau) and cellular (diffusive-thermal) instabilities for sufficiently large value of the gravity parameter for the case of downward-propagating flames. The stability analysis reveals that in the case of oscillatory diffusive-thermal instability, the flame front cannot be stabilized in the similar way. Finally, the stability results are compared satisfactorily with unsteady numerical simulations.