Numerical studies of flames in wide tubes: stability limits of curved stationary flames

Flame dynamics in wide tubes with ideally adiabatical and slip walls is studied by means of direct numerical simulations of the complete set of hydrodynamical equations including thermal conduction, fuel diffusion, viscosity, and chemical kinetics. Stability limits of curved stationary flames in wide tubes and the hydrodynamic instability of these flames (the secondary Darrieus-Landau instability) are investigated. The stability limits found in the present numerical simulations are in a very good agreement with the previous theoretical predictions. It is obtained that close to the stability limits the secondary Darrieus-Landau instability results in an extra cusp at the flame front. It is shown that the curved flames subject to the secondary Darrieus-Landau instability propagate with velocity considerably larger than the velocity of the stationary flames.     

Nonlinear equation for curved nonstationary flames and flame stability

A time-dependent nonlinear equation for a nonstationary curved flame front of an arbitrary expansion coefficient is derived under the assumptions of a small but finite flame thickness and weak nonlinearity. On the basis of the derived equation, stability of two-dimensional curved stationary flames propagating in tubes with ideally adiabatic and slip walls is studied. The stability analysis shows that curved stationary flames become unstable for sufficiently wide tubes. The obtained stability limits are in a good agreement with the results of numerical simulations of flame dynamics and with semiqualitative stability analysis of curved stationary flames. Possible outcomes of the obtained instability at the nonlinear stage are discussed. The instability may result in extra wrinkles at a flame front close to the stability limits and in self-turbulization of the flame far from the limits. The self-turbulization can also be interpreted as a fractal structure. The fractal dimension of a flame front and velocity of a self-turbulized flame are evaluated.     

Effect of vorticity production on the structure and velocity of curved flames

A nonlinear equation describing curved stationary flames with arbitrary gas expansion theta = rho(fuel)/rho(burnt) is obtained in a closed form without an assumption of weak nonlinearity. The equation respects all conservation laws and takes into account vorticity production in the flame. In the scope of the asymptotic expansion for theta -->1, the new equation solves the problem of stationary flame propagation with accuracy of the sixth order in theta - 1. Its analytical solutions give the flame velocity in tubes of arbitrary width, which agrees with available results of direct numerical modeling.     

Flame fingers and interactions of hydrodynamic and thermodiffusive instabilities in laminar lean hydrogen flames

Lean hydrogen/air flames are prone to hydrodynamic and thermodiffusive instabilities. In this work, the contribution of each instability mechanism is quantified separately by performing detailed simulations of laminar planar lean hydrogen/air flames with different diffusivity models and equations of state to selectively suppress the hydrodynamic or thermodiffusive instability mechanism.From the analysis of the initial phase of the simulations, the thermodiffusive instability is shown to dominate the flame dynamics. If differential diffusion and, hence, the thermodiffusive instability is suppressed, the flame features a strong reduction of the instability growth rates, whereas if present, a wide range of unstable wave numbers is observed due to the strong destabilizing nature of differential diffusion. When instabilities are fully developed, lean hydrogen/air flames feature the formation of small-scale cellular structures and large-scale flame fingers. While the size of the former is known to be close to the most unstable wave length of a linear stability analysis, this work shows that flame fingers also originate from the thermodiffusive instability and most noteworthy, are not linked to an interaction of the two instability mechanisms. They are stable with respect to external perturbations and feature an enhanced flame propagation as the formation of a central cusp at their tip enables the co-existence of two strongly curved leading edges with high reactivity. The thermodiffusive instability is shown to significantly affect the flames’ consumption speed, while the consumption speed enhancement caused by the hydrodynamic instability is significantly smaller. Further, the surface area increase due to wrinkling is strongly diminished if one of the two instability mechanisms is missing. This is linked to a synergistic interaction between the two mechanisms, as the propagation of flame fingers is enhanced by the presence of the hydrodynamic instability due to a widening of the streamlines ahead of the flame fingers.     

Turbulent propagation of premixed flames in the presence of Darrieus–Landau instability

We investigate the role played by hydrodynamic instability in the wrinkled flamelet regime of turbulent combustion, where the intensity of turbulence is small compared to the laminar flame speed and the scale large compared to the flame thickness. To this end the Michelson–Sivashinsky (MS) equation for flame front propagation in one and two spatial dimensions is studied in the presence of uncorrelated and correlated noise representing a turbulent flow field. The combined effect of turbulence intensity, integral scale, and an instability parameter related to the Markstein length are examined and turbulent propagation speed monitored for both stable planar flames and corrugated flames for which the planar conformation is unstable. For planar flames a particularly simple scaling law emerges, involving quadratic dependence on intensity and a linear dependence on the degree of instability. For corrugated flames we find the dependence on intensity to be substantially weaker than quadratic, revealing that corrugated flames are more resilient to turbulence than planar flames. The existence of a threshold turbulence intensity is also observed, below which the corrugated flame in the presence of turbulence behaves like a laminar flame. We also analyze the conformation of the flame surface in the presence of turbulence, revealing primary, large-scale wrinkles of a size comparable to the main corrugation. When the integral scale is much smaller than the characteristic corrugation length we observe, in addition to primary wrinkles, secondary small-scale wrinkles contaminating the surface. The flame then acquires a multi-scale, self-similar conformation, with a fractal dimension, for one-dimensional flames, plateauing at 1.23 for large intensities. The existence of an intermediate integral scale is also found at which the turbulent speed is maximized. When two-dimensional flames are subject to turbulence, the primary wrinkling patterns give rise to polyhedral-cellular structures which bear a very close resemblance to those observed in experiments on hydrodynamically unstable expanding spherical flames.     

Numerical simulation of flames as gas-dynamic discontinuities

The dynamics of thin premixed flames is computationally studied within the context of a hydrodynamic theory. A level-set method is used to track down the flame, which is treated as a free-boundary interface. The flow field is described by the incompressible Navier–Stokes equations, with different densities for the burnt and unburnt gases, supplemented by singular source terms that properly account for thermal expansion effects. The numerical scheme has been tested on several benchmark problems and was shown to be stable and accurate. In particular, the propagation of a planar flame front and the dynamics of hydrodynamically unstable flames were successfully simulated. This includes recovering the planar front in narrow domains, the Darrieus–Landau linear growth rate for long waves of small amplitude, and the nonlinear development of cusp-like structures predicted by the Michelson–Sivashinsky equation for a small density change. The stationary flame of a Bunsen burner with uniform and parabolic outlet flows were also simulated, showing in particular a careful mapping of the flow field. Finally, the evolution of a hydrodynamically unstable flame was studied for finite amplitude disturbances and realistic values of thermal expansion. These results, which constitute one of the main objectives of this study, elucidate the effect of thermal expansion on flame dynamics.     

The unstable behavior of cellular premixed flames induced by intrinsic instability

The unstable behavior of cellular premixed flames induced by intrinsic instability is studied by two-dimensional unsteady calculations of reactive flows. In the present numerical simulation, the compressible Navier–Stokes equation including a one-step irreversible chemical reaction is employed. We consider two basic types of phenomena to account for the intrinsic instability of premixed flames, i.e., hydrodynamic and diffusive-thermal effects. The hydrodynamic effect is caused by the thermal expansion through the flame front; the diffusive-thermal effect is caused by the preferential diffusion of mass versus heat. A disturbance with several wavelength components is superimposed on a planar flame, and the formation of a cellular flame induced by hydrodynamic and diffusive-thermal effects is numerically simulated. After the cellular-flame formation, the combination and division of cells are observed. The behavior of cellular-flame fronts becomes more unstable when the Lewis number is lower than unity, since the diffusive-thermal effect has a great influence on the unstable behavior. The cell size changes with time, and its average is greater than the critical wavelength and becomes smaller by decreasing the Lewis number. The flame velocity of cellular flames depends strongly on the length of computational domain in the direction tangential to the flame front. As the length of computational domain increases, the flame velocity becomes larger. This is because the long-wavelength components of disturbances play an important role in the shape of cellular flames, i.e., in the flame-surface area.     

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.     

Effect of geometrical parameters on thermo-acoustic instability of downward propagating flames in tubes

Experiments and theoretical analysis are presented to clarify the effect of geometrical parameters on thermo-acoustic instability of downward propagating flames in tubes. The experiments reveal that the longer tubes have higher instability compared to shorter tubes and the lower diameter tubes have higher instability compared to higher diameter tubes. The secondary instability leading to turbulent burning is found to be more sensitive to change in geometrical parameters compared to primary instability (oscillating flat flame). The secondary instability is re-stabilized for some intermediate burning velocity conditions even though lower and higher burning velocity conditions show secondary instability. The appearance of such re-stabilization is only observed for some specific lengths of the tube. Present experimental observations pertaining to the effect of geometrical parameters is found to be contradicting the theoretical predictions based on pressure coupling mechanism. To clear the underlying mechanism, analytical growth rate is computed considering velocity coupling mechanism. The computed growth rates correctly predict the effect of geometrical parameters on thermo-acoustic instability of downward propagating flames. This work provides further evidence to believe that the flame -acoustic coupling in downward propagating flames is due to flame area modulation (leading to heat release modulation) through action of acoustic acceleration.     

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.     

Large scale effects in weakly turbulent premixed flames

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

Strain rates,flow patterns and flame surface densities in hydrodynamically unstable,weakly turbulent premixed flames

Recent numerical and experimental studies have unveiled a potentially marked difference between the laminar as well as turbulent propagation of premixed flames exhibiting Darrieus–Landau (DL) (or hydrodynamic) instabilities from flames for which instabilities are inhibited. In this study we utilize two-dimensional numerical simulations of slot burner flames as well as experimental Propane–Air Bunsen flames to analyse differences in turbulent propagation, strain rate and induced flow patterns of hydrodynamically stable and unstable flames. We also investigate the effects of hydrodynamic instability on quantities which are directly related to reaction rate closure models, such as flame surface density and stretch factor. A clear enhancement of turbulent flame speed can be observed for unstable flames, generally mitigated at higher turbulence intensity, which is attributed to a flame area increase induced by the characteristic cusp-like DL-induced corrugation, absent in stable flames, which occurs concurrently and in synergy with turbulent wrinkling. Unstable flames also exhibit, both numerically and experimentally, a different correlation between strain rate and flame curvature and are observed to give rise to a channeling of the induced flow in the fresh mixture. Conditionally averaged flame surface density is also observed to attain smaller values in unstable flames, as a result of the thicker turbulent flame brush, indicating that closure models should incorporate instability-related parameters in addition to turbulence-related parameters.     

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

Reaction propagation modes in millimeter-scale tubes for ethylene/oxygen mixtures

Effects of tube diameter and equivalence ratio on reaction front propagations of ethylene/oxygen mixtures in capillary tubes were experimentally analyzed using high speed cinematography. The inner diameters of the tubes investigated were 0.5, 1, 2 and 3 mm. The flame was ignited at the center of the 1.5 m long smooth tube under ambient pressure and temperature before propagated towards the exits in the opposite directions. A total of five reaction propagation scenarios, including deflagration-to-detonation transition followed by steady detonation wave transmission (DDT/C–J detonation), oscillating flame, steady deflagration, galloping detonation and quenching flame, were identified. DDT/C–J detonation mode was observed for all tubes for equivalence ratios in the vicinity of stoichiometry. The velocity for the steady detonation wave propagation was approximately Chapman–Jouguet velocity for 1, 2, and 3 mm I.D. tubes; however, a velocity deficit of 5% was found for the case in 0.5 mm I.D. tube. For leaner mixtures, an oscillating flame mode was found for tubes with diameters of 1 to 3 mm, and the reaction front travelled in a steady deflagrative flame mode with velocities around 2–3 m/s when the mixture equivalence ratio becomes even leaner. Galloping detonation wave propagation was the dominant mode for the fuel lean regime in the 0.5 mm I.D. tube. For rich mixtures beyond the detonation limits, a fast flame followed by flame quenching was observed.     

Effect of velocity gradient on propagation speed of tribrachial flames in laminar coflow jets

The effect of velocity gradient on the propagation speed of tribrachial flame edge has been investigated experimentally in laminar coflow jets for propane fuel. It was observed that the propagation speed of tribrachial flame showed appreciable deviations at various jet velocities in high mixture fraction gradient regime. From the similarity solutions, it was demonstrated that the velocity gradient varied significantly during the flame propagation. To examine the effect of velocity gradient, detail structures of tribrachial flames were investigated from OH LIF images and Abel transformed images of flame luminosity. It was revealed that the tribrachial point was located on the slanted surface of the premixed wing, and this slanted angle was correlated with the velocity gradient along the stoichiometric contour. The temperature field was visualized qualitatively by the Rayleigh scattering image. The propagation speed of tribrachial flame was corrected by considering the direction of flame propagation with the slanted angle and effective heat conduction to upstream. The corrected propagation speed of tribrachial flame was correlated well. Thus, the mixture fraction gradient together with the velocity gradient affected the propagation speed.     

预混气体燃烧火焰闪烁现象分析

在低速射流的预混火焰和扩散火焰中都存在火焰闪烁现象。对扩散火焰,其机理已比较明确,是由于浮力诱导引起的一种水力学不稳定性。而对预混火焰闪烁现象则存在水力学不稳定性和热驱动不稳定性两种观点。本文根据水力学不不稳定性观点,把预混火焰的闪烁现象看成是包围火焰锋面的已燃混气层中内、外区间在垂直方向上的相对脉动,应用Ｋｅｌｖｉｎ－Ｈｅｌｍｈｏｌｔｚ不稳定性机理进行了分析,获得了火焰闪烁频率与重力和压力的关系式,并与已有的结果作了对比。     

Velocity of weakly turbulent flames of finite thickness

The velocity increase of a weakly turbulent flame of finite thickness is investigated using analytical theory developed in previous papers. The obtained velocity increase depends on the flow parameters: on the turbulent intensity, on the turbulent spectrum and on the characteristic length scale. It also depends on the thermal and chemical properties of the burning matter: thermal expansion, the Markstein number and the temperature dependence of transport coefficients. It is shown that the influence of the finite flame thickness is especially strong close to the resonance point, when the wavelength of the turbulent harmonic is equal to the cut off wavelength of the Darrieus–Landau instability. The velocity increase is almost independent of the Prandtl number. On the contrary, the Markstein number is one of the most important parameters controlling the velocity increase. The relative role of the external turbulence and the Darrieus–Landau instability for the velocity increase is studied for different parameters of the flow and the burning matter. The velocity increase for turbulent flames in methane and propane fuel mixtures is calculated for different values of the equivalence ratio. The present theoretical results are compared with previous experiments on turbulent flames. In order to perform the comparison, the theoretical results of the present paper are extrapolated to the case of a strongly corrugated flame front using the ideas of self-similar flame dynamics. The obtained theoretical results are in a reasonable agreement with the experimental data, taking into account the uncertainties of both the theory and the experiments. It is shown that in many experiments on turbulent flames the Darrieus–Landau instability is more important for the flame velocity than the external turbulence.     

Premixed turbulent flame front structure investigation by Rayleigh scattering in the thin reaction zone regime

Premixed turbulent flames of methane–air and propane–air stabilized on a bunsen type burner were studied using planar Rayleigh scattering and particle image velocimetry. The fuel–air equivalence ratio range was from lean 0.6 to stoichiometric for methane flames, and from 0.7 to stoichiometric for propane flames. The non-dimensional turbulence rms velocity, u′/S_L, covered a range from 3 to 24, corresponding to conditions of corrugated flamelets and thin reaction zones regimes. Flame front thickness increased slightly with increasing non-dimensional turbulence rms velocity in both methane and propane flames, although the flame thickening was more prominent in propane flames. The probability density function of curvature showed a Gaussian-like distribution at all turbulence intensities in both methane and propane flames, at all sections of the flame.The value of the term Dκ, the product of molecular diffusivity evaluated at reaction zone conditions and the flame front curvature, has been shown to be smaller than the magnitude of the laminar burning velocity. This finding questions the validity of extending the level set formulation, developed for corrugated flames region, into the thin reaction zone regime by increasing the local flame propagation by adding the term Dκ to laminar burning velocity.     

Exact equation for curved stationary flames with arbitrary gas expansion

An exact equation describing freely propagating stationary flames with arbitrary values of the gas expansion coefficient is obtained. This equation respects all conservation laws at the flame front, and provides a consistent nonperturbative account of the effect of vorticity produced by the curved flame on the front structure. It is verified that the new equation is in agreement with the approximate equations derived previously in the case of weak gas expansion.