Large eddy simulation of premixed hydrogen/methane/air flame propagation in a closed duct

In this paper, large eddy simulation (LES) is performed to investigate the propagation characteristics of premixed hydrogen/methane/air flames in a closed duct. In LES, three stoichiometric hydrogen/methane/air mixtures with hydrogen fractions (volume fractions) of 0, 50% and 100% are used. The numerical results have been verified by comparison with experimental data. All stages of flame propagation that occurred in the experiment are reproduced qualitatively in LES. For fuel/air mixtures with hydrogen fractions of 0 and 50%, only four stages of “tulip” flame formation are observed, but when the hydrogen fraction is 100%, the distorted “tulip” flame appears after flame front inversion. In the acceleration stage, the LES and experimental flame speed and pressure dynamic coincide with each other, except for a hydrogen fraction of 0. After “tulip” flame formation, all LES and experimental flame propagation speeds and pressure dynamics exhibit the same trends for hydrogen fractions of 0 and 100%. However, when the hydrogen fraction is 50%, a slight periodic oscillation appears only in the experiment. In general, the different structures displayed in the flame front during flame propagation can be attributed to the interaction between the flame front, the vortex and the reverse flow formed in the unburned and burned zones.     

Dynamics of premixed hydrogen/air flame in a closed combustion vessel

The dynamics of a premixed hydrogen/air flame propagating in a closed vessel is investigated using high-speed schlieren cinematography, pressure measurement and numerical simulation. A dynamically thickened flame approach with a 19-step detailed chemistry is employed in the numerical simulation to model the premixed combustion. The schlieren photographs show that a remarkable distorted tulip flame is initiated after a classical tulip flame has been fully produced. A second distorted tulip flame is generated with a cascade of indentations created in succession before the vanishing of the first one. The flame dynamics observed in the experiments is well reproduced in the numerical simulation. The burnt region near the flame front is entirely dominated by a reverse flow during the formation of the distorted tulip flame. The distorted tulip flame can be formed in the absence of vortex motion. The pressure wave leads to periodic flame deceleration and plays an essential role in the distorted tulip formation. The numerical results corroborate the mechanism that the distorted tulip flame formation is a manifestation of Taylor instability.     

An experimental study of premixed hydrogen/air flame propagation in a partially open duct

High-speed schlieren cinematography and pressure records are used to investigate the dynamics of premixed hydrogen/air flame propagation and pressure build up in a partially open duct with an opening located in the upper wall near the right end of the duct. This work provides basic understanding of flame behaviors and the effects of opening ratio on the combustion dynamics. The flame behaves differently under different opening conditions. The opening ratio has an important influence on the flame propagation and pressure dynamics. When the opening ratio α ≤ 0.075 a significant distorted tulip flame can be formed after the full formation of a classical tulip flame. The propagation speed of flame leading tip increases with the opening ratio. The coupling of flame front with the pressure wave is strong at low opening ratio. Both the pressure growth rate and oscillation amplitude inside the duct increases as the opening ratio decreases. The formation times of tulip and distorted tulip flames and the corresponding distances of flame front increase with the increase of the opening ratio.     

Experimental study on the behaviors and shape changes of premixed hydrogen-air flames propagating in horizontal duct

The behaviors and shape changes of premixed hydrogen-air flames at various equivalence ratios propagating in half-open and closed horizontal ducts are experimentally investigated using high-speed schlieren imaging and pressure sensors. The study shows that the premixed hydrogen-air flame undergoes more complex shape changes and exhibits more distinct characteristics than that of other gaseous fuels. One of the outstanding findings is that obvious distortion happens to tulip flame after its full formation when equivalence ratio ranges from 0.84 to 4.22 in the closed duct. The salient tulip flame distortions are specially scrutinized and distinguished from the classical tulip collapse and disappearance. The dynamics of distorting tulip flame is different from that of classical tulip flame. The normal tulip flame can be reproduced after the first distortion followed by another distortion. The initiation of flame shape changes coincides with the deceleration both of pressure rise and flame front speed for flames with tulip distortions. And the formation and dynamics of tulip/distorting tulip flames depend on the mixture composition.     

On the transition from a highly turbulent curved flame into a tulip flame

Experimental and numerical investigations of premixed flame propagation behaviour associated with vortex interactions due to planar pressure waves crossing a curved flame front have been carried out. The resulting “tulip flame” formation in such a closed tube has been studied by Schlieren visualization. The “tulip flame” phenomenon was observed only in closed tubes, while cellular flame fronts appeared in half-open tubes. A physical model has been developed and implemented in a discrete vortex method combined with a flame tracking algorithm. The numerical method has been applied to model and understand the processes that cause the flame to change from a curved to a tulip shape. The results of the simulation are in good agreement with the experimental observations. We find that the rotational flows causing the tulip formation in our experimental case originate from the baroclinic effect — an interaction of non-parallel density and pressure gradients. Pressure waves were generated ahead of the accelerating and highly turbulent flame front. In closed tubes the pressure waves were reflected and crossed the curved flame front. As a result we saw the “tulip flame”. Within 0.5 ms, the flame front velocity reversed from about 50 m/s to about −20 m/s.     

Hydrogen-air deflagrations in open atmosphere: Large eddy simulation analysis of experimental data

The largest known experiment on hydrogen-air deflagration in the open atmosphere has been analysed by means of the large eddy simulation (LES). The combustion model is based on the progress variable equation to simulate a premixed flame front propagation and the gradient method to decouple the physical combustion rate from numerical peculiarities. The hydrodynamic instability has been partially resolved by LES and unresolved effects have been modelled by Yakhot's turbulent premixed combustion model. The main contributor to high flame propagation velocity is the additional turbulence generated by the flame front itself. It has been modelled based on the maximum flame wrinkling factor predicted by Karlovitz et al. theory and the transitional distance reported by Gostintsev with colleagues. Simulations are in a good agreement with experimental data on flame propagation dynamics, flame shape, and outgoing pressure wave peaks and structure. The model is built from the first principles and no adjustable parameters were applied to get agreement with the experiment.     

Propagation of premixed hydrogen-air flame initiated by a planar ignition in a closed tube

Numerical simulations were used to study the dynamics of premixed flames propagating after planar ignition in a closed tube filled with stoichiometric hydrogen-air mixture. The two-dimensional fully compressible reactive Navier–Stokes equations coupled to a calibrated chemical-diffusive model were solved using a high-order numerical method and adaptive mesh refinement. The results show that the flame evolves from an initially planar flame to a double-cusped tulip flame, subsequently to a multi-cusped tulip flame, and finally to a series of distorted tulip flames (DTFs). The DTF forms one after another until the end of combustion. The initial flame lips of the double-cusped tulip flame are produced due to the stretching effect of nonuniform flow caused by the wall friction. The multi-cusped tulip flame forms as secondary cusps are created on the leading flame tips near the sidewalls. The formation of DTFs here is thought to be closely connected to pressure waves generated in the flame propagation process. The first DTF is caused by the combined effects of the vortex motions and the Rayleigh–Taylor (RT) instability driven by pressure waves, while the subsequent DTFs form due to reverse flows and RT instability. Nevertheless, both the vortex motions and reverse flows are essentially induced by the interactions between pressure waves and flow fields. Furthermore, the numerical results were compared to that in the case with a semicircular ignition. It was found that although there are significant differences in the early flame acceleration and tulip formation stages between the two differently shaped ignitions, the dynamics of DTFs are substantially consistent.     

Experimental and numerical study on premixed hydrogen/air flame propagation in a horizontal rectangular closed duct

Hydrogen is a promising energy in the future, and it is desirable to characterize the combustion behavior of its blends with air. The premixed hydrogen/air flame microstructure and propagation in a horizontal rectangular closed duct were recorded using high-speed video and Schlieren device. Numerical simulation was also performed on Fluent CFD code to compare with the experimental result. A tulip flame is formed during the flame propagating, and then the tulip flame formation mechanism was proposed based on the analysis. The induced reverse flow and vortex motion were observed both in experiment and simulation. The interactions among the flame, reverse flow and vortices in the burned gas change the flame shape and ultimately it develops into a tulip flame. During the formation of the tulip flame, the tulip cusp slows down and stops moving after its slightly forward moving, and then, it starts to move backward and keeps on a longer time, after that, it moves forward again. The structure of the tulip flame is becoming less stable with its length decreasing in flame propagation direction. The flame thickness increases gradually which is due to turbulence combustion.     

Effects of ignition location on premixed hydrogen/air flame propagation in a closed combustion tube

The dynamics of premixed hydrogen/air flame ignited at different locations in a finite-size closed tube is experimentally studied. The flame behaves differently in the experiments with different ignition positions. The ignition location exhibits an important impact on the flame behavior. When the flame is ignited at one of the tube ends, the heat losses to the end wall reduce the effective thermal expansion and moderate the flame propagation and acceleration. When the ignition source is at a short distance off one of the ends, the tulip flame dynamics closely agrees with that in the theory. And both the tulip and distorted tulip flames are more pronounced than those in the case with the ignition source placed at one of the ends. Besides, the flame–pressure wave coupling is quite strong and a second distorted tulip flame is generated. When the ignition source is in the tube center, the flame propagates in a much gentler way and the tulip flame can not be formed. The flame oscillations are weaker since the flame–pressure wave interaction is weaker.     

A filtered tabulated chemistry model for LES of premixed combustion

A new modeling strategy called F-TACLES (Filtered Tabulated Chemistry for Large Eddy Simulation) is developed to introduce tabulated chemistry methods in Large Eddy Simulation (LES) of turbulent premixed combustion. The objective is to recover the correct laminar flame propagation speed of the filtered flame front when subgrid scale turbulence vanishes as LES should tend toward Direct Numerical Simulation (DNS). The filtered flame structure is mapped using 1-D filtered laminar premixed flames. Closure of the filtered progress variable and the energy balance equations are carefully addressed in a fully compressible formulation. The methodology is first applied to 1-D filtered laminar flames, showing the ability of the model to recover the laminar flame speed and the correct chemical structure when the flame wrinkling is completely resolved. The model is then extended to turbulent combustion regimes by including subgrid scale wrinkling effects in the flame front propagation. Finally, preliminary tests of LES in a 3-D turbulent premixed flame are performed.     

Experimental study and three-dimensional simulation of premixed hydrogen/air flame propagation in a closed duct

An experimental and numerical study of premixed hydrogen/air flame propagation in a closed duct is presented. High-speed schlieren photography is used in the experiment to record the changes in flame shape and location. The pressure transient during the combustion is measured using a pressure transducer. A dynamic thickened flame model is applied to model the premixed combustion in the numerical simulation. The four stages of the flame dynamics observed in the experiment are well reproduced in the numerical simulation. The oscillations of the flame speed and pressure growth, induced by the pressure wave, indicate that the pressure wave plays an important role in the combustion dynamics. The predicted pressure dynamics in the numerical simulation is also in good agreement with that in the experiment. The close correspondence between the numerical simulation and experiment demonstrate that the TF approach is quite reliable for the study of premixed hydrogen/air flame propagation in the closed duct. It is shown that the flame wrinkling is important for the flame dynamics at the later stages.     

Impacts of plate slits on flame acceleration of premixed methane/air in a closed tube

The paper aims at revealing the interaction of various numbers of premixed methane/air jet flames in a closed duct. In the experiment, a high-speed video camera and pressure transducers are used to study the flame structure and pressure dynamics. In the numerical simulations, large eddy simulation (LES) with Power-Law combustion model is employed to investigate the interaction between the moving flame and vortices induced by the thin plate. The results demonstrate that the flame propagation for all plate configurations can be divided into four typical stages, i.e. hemispherical flame, finger-shaped flame, jet flame and bidirectional propagation flame. For three plate configurations, the jet flames merge together under the effect of the vortices, and the more slits with the same blockage ratio (BR) do not mean the stronger deflagration. It is observed that the peaks of flame tip speed and pressure growth rate decrease with the increase of the number of slits. The sub-grid scale combustion model, Power-Law model, coupled with sub-grid scale viscosity model, dynamic Smagorinsky-Lilly eddy viscosity model can well reproduce the flame propagation. By analyzing the numeric flow structure, the flame propagation mechanism of premixed methane/air flame propagation in a tube with various slits can be explained in the view of pure hydrodynamics.     

Dynamics of premixed hydrogen-air flame propagation in the duct with pellets bed

Hydrogen, as the promising clean alternative energy in the future, is in the spotlight now all over the world. However, its flammable and explosive hazards should be highly considered during its practical application. In this study, the experiments are performed to study premixed hydrogen-air flame propagation in the duct with pellets bed, especially for fuel-rich condition. High-speed schlieren photography is employed to capture flame front development during the experiments. As well as the pressure transducer, is used to track the pressure buildup in the flame propagation process. Different diameters of pellets and different concentrations of gas mixture are considered in this experimental study. The typical evolutions about the tulip flame are similar in all cases, although the tulip flame formation time caused by the laminar flame speed are different. The flame propagation velocity is pretty enhanced in fuel-lean mixture under the effect of large diameter pellets bed, but it is significantly suppressed in fuel-rich conditions. While for the small diameter pellets (d = 3 mm), the suppression effect on flame propagation and pressure is obtained over a wider range of equivalence ratios, especially a better suppression effect is generated near the stoichiometric condition.     

Effects of the external turbulence on centrally-ignited spherical unstable Ch4/H2/air flames in the constant-volume combustion bomb

In the present study, we conducted experiments to investigate the effects of external turbulence on the development of spherical H₂/CH₄/air unstable flames developments at two different equivalence ratios associated with different turbulent intensities using a spherical constant-volume turbulent combustion bomb and high speed schlieren photography technology. Flame front morphology and acceleration process were recorded and different effects of weak external turbulent flow field and intrinsic flame instability on the unstable flame propagation were compared. Results showed the external turbulence has a great influence on the unstable flame propagation under rich fuel conditions. For fuel-lean premixed flames, however, the effects of external turbulence on the morphology of the cellular structure on the flame front was not that obvious. Critical radius decreased firstly and then kept almost unchanged with the augment of the turbulence intensity. This indicated the dominating inhibiting effect of flame stretch on the turbulent premixed flame at the initial stage of the flame front development. Beyond the critical radius, the acceleration exponent was found increasing with the enhancement of initial turbulence intensity for fuel-lean premixed flames. For fuel-rich conditions, however, the initial turbulence intensity had little effect on acceleration exponent. In order to evaluate the important impact of the intrinsic flame instability and external turbulent flow field for spherical propagating premixed flames, intrinsic flame instability scale and average diameter of vortex tube were calculated. Intrinsic flame instability scale decreased greatly and then stayed unchanged with the propagation of the flame front. The comparison between intrinsic flame instability scale and average diameter of vortex tube demonstrated that the external turbulent flow filed will be more important for the evolution of wrinkle structure in the final stage of the flame propagation, when the turbulence intensity was more than 0.404 m/s.     

Experimental study of hydrogen/air premixed flame propagation in a closed channel with inhibitions for safety consideration

An experimental study of hydrogen/air premixed flame propagation in a closed rectangular channel with the inhibitions (N₂ or CO₂) was conducted to investigate the inhibiting effect of N₂ and CO₂ on the flame properties during its propagation. Both Schlieren system and the pressure sensor were used to capture the evolution of flame shape and pressure changes in the channel. It was found that both N₂ and CO₂ have considerable inhibiting effect on hydrogen/air premixed flames. Compared with N₂, CO₂ has more prominent inhibition, which has been interpreted from thermal and kinetic standpoints. In all the flames, the classic tulip shape was observed. With different inhibitor concentration, the flame demonstrated three types of deformation after the classic tulip inversion. A simple theoretical analysis has also been conducted to indicate that the pressure wave generated upon the first flame-wall contact can affect the flame deformation depending on its meeting moment with the flame front. Most importantly, the meeting moment is always behind the start of tulip inversion, which suggests the non-dominant role of pressure wave on this featured phenomenon.     

用光学纹影摄影术观察分析定容燃烧室内氢燃料的燃烧过程

论述了采用纹影摄影术和高速摄影法观察分析氢气和空气预混合燃料在定容燃烧室内的火花点火燃烧过程,定性地分析了预混合氢气燃料的火焰形态和变化过程,以及燃烧室内的初始压力和空燃比对火焰传播速度及其燃烧压力的影响,通过采用纹影摄影术方法,初步揭示了预混合氢气燃料在定容燃烧室内燃烧时火焰初期紊流产生的机理,以及由开始的层流状火焰发展到湍流状火焰的过程,研究结果表明,预混合氢气燃料燃烧的火焰传播速度及燃烧压力明显地受燃烧室内的初始压力和空燃比的影响。     

Effect of obstacle arrangement on premixed hydrogen flame: Eddy-dissipation concept model based numerical simulation

In this paper, computational fluid dynamics (CFD) numerical simulation is used to analyze and discuss the horizontal propagation process of premixed hydrogen flame with obstacles. A total of three different obstacle channel arrangements at the blocking ratio of 0.5, which will affect the explosion flame and pressure development. The results show that the premixed flame is affected by flow instabilities and vortices when propagating through the obstacle channel, thereby distorting the flame. The vortices outside the flame boundary are more conducive to the acceleration of the flame. The continuous acceleration and synergistic promotion of the flame is more prominent due to the existence of the channel in the central axis of flame propagation, and the maximum velocity even achieved 307.91  m/s. The degree of the wrinkle of flame increases with the number of obstacle channels. The flame propagation process is always accompanied by pressure variations, and the dynamic pressure builds up at the flame front and intensifies periodically. But the downstream pressure gradually increases as the number of obstacle channels increases. CFD simulation of the explosion process clearly reveals the changing trends and interactions of explosion characteristic factors.     

Effects of non-equidiffusion on unsteady propagation of hydrogen-enriched methane/air premixed flames

A Large Eddy Simulation (LES) model was developed to simulate the unsteady propagation of hydrogen-enriched methane/air premixed flames around toroidal vortices. Although the LES model does not take into account the non-equidiffusive effects associated with the hydrogen presence (preferential diffusion and non-unity Lewis number), it gives good predictions of experimental data previously obtained for lean mixtures with hydrogen mole fraction in the fuel (hydrogen plus methane) varying from 0 to 0.5. In particular, for each fuel composition, size and velocity of the toroidal vortex generated ahead of the propagating flame front are well reproduced along with the evolution of the flame shape and structure resulting from the interaction with the vortex. The negligible role played by the non-equidiffusive effects has been attributed to the fact that, at the conditions investigated, the characteristic time of hydrogen diffusion is one order of magnitude higher than the characteristic time of flame roll-up around the vortex.     

Numerical investigation on the effects of reaction orders on the flame propagation dynamic behaviors for premixed gas in a closed tube

In this paper, the premixed flame propagation in a closed tube is surveyed using Computational Fluid Dynamics. The propagation characteristics of premixed flame are obtained coupling a single-step reaction mechanism with a laminar flame model. Three single-step reaction mechanisms are established with different reaction orders for hydrocarbon fuels. This study is to establish a wider range of reaction mechanisms and represent actual experimental conditions better. The numerical simulation results demonstrate that reaction orders can affect the tulip flame development. As the flame spreads, the tulip flame fronts become wrinkled. When the reaction order is 2, there are more wrinkles in the flame front and the degree of wrinkles is more obvious. Reaction orders also affect the flame tip velocity and the flame skirt velocity. The main reason is that laminar flame speeds are significantly different. When the reaction orders are 1.5 and 2, laminar flame speeds are mainly affected by temperature, which respectively increase by about 25% and 75%. When the reaction order is 1, the pressure is crucial for the variation of laminar flame speed. The laminar flame speed decreases by about 33%.     

Experimental investigation on the self-acceleration of 10%H2/90%CO/air turbulent expanding premixed flame

The self-acceleration characteristics of a syngas/air mixture turbulent premixed flame were experimentally evaluated using a 10% H₂/90% CO/air mixture turbulent premixed flame by varying the turbulence intensity and equivalence ratio at atmospheric pressure and temperature. The propagation characteristics of the turbulent premixed flame including the variation in the flame propagation speed and turbulent burning velocity of the syngas/air mixture turbulent premixed flame were evaluated. In addition, the effect of the self-acceleration characteristics of the turbulent premixed flame was also evaluated. The results show that turbulence gradually changes the radius of the premixed flame from linear growth to nonlinear growth. With the increase of turbulence intensity, the formation of a cellular structure of the flame front accelerated, increasing the flame propagation speed and burning speed. In the transition stage, the acceleration exponent and fractal excess of the turbulent premixed flame decreased with increasing equivalence ratio and increased with increasing turbulence intensity at an equivalence ratio of 0.6. The acceleration exponent was always greater than 1.5.