Experimental and numerical investigation of premixed flame propagation with distorted tulip shape in a closed duct

High-speed schlieren photography, pressure records and large eddy simulation (LES) model are used to study the shape changes, dynamics of premixed flame propagation and pressure build up in a closed duct. The study provides further understanding of the interaction between flame front, pressure wave and combustion-generated flow, especially when the flame acquires a “distorted tulip” shape. The Ulster multi-phenomena LES premixed combustion model is applied to gain an insight into the phenomenon of “distorted tulip” flame and explain the experimental observations. The model accounts for the effects of flow turbulence, turbulence generated by flame front itself, selective diffusion, and transient pressure and temperature on the turbulent burning velocity. The schlieren images show that the flame exhibits a salient “distorted tulip” shape with two secondary cusps superimposed onto the two original tulip lips. This curious flame shape appears after a well-pronounced classical tulip flame is formed. The dynamics of “distorted tulip” flame observed in the experiment is well reproduced by LES. The numerical simulations show that large-scale vortices are generated in the burnt gas after the formation of a classical tulip flame. The vortices remain in the proximity of the flame front and modify the flow field around the flame front. As a result, the flame front in the original cusp and near the sidewalls propagates faster than that close to the centre of the original tulip lips. The discrepancy in the flame propagation rate finally leads to the formation of the “distorted tulip” flame. The LES model validated previously against large-scale hydrogen/air deflagrations is successfully applied in this study to reproduce the dynamics of flame propagation and pressure build up in the small-scale duct. It is confirmed that grid resolution has an influence to a certain extent on the simulated combustion dynamics after the flame inversion.     

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.     

Experimental study on the propagation characteristics of hydrogen/methane/air premixed flames in a narrow channel

This work is focused on the explosion characteristics of premixed gas containing different volume fractions of hydrogen in a narrow channel (1000 mm × 50 mm × 10 mm) under the circumstance of stoichiometric ratio. The ignition positions were set in the closed end and the middle of the pipeline respectively. The results showed that when the gas was ignited at the pipeline closed end, the propagating flame was tulip structure for different premixed gas. When the hydrogen volume fraction was less than 40%, the flame propagation speed increased significantly with the rise of hydrogen volume fraction, and the overpressure peak also appeared obviously in advance. However, when the volume fraction of hydrogen was more than 40%, the increase of flame propagation speed and the overpressure peak occurrence time varied slightly. Furthermore, when the ignition position was placed in the middle of the pipeline, the flame propagation speed propagating to the opening end was much faster than that propagating to the closing end, and there was no tulip shape when the flame propagates to the opening end. The flame propagating to the closed end appeared tulip shape under the influence of airflow, and high-frequency flame oscillation occurred during the propagation. This work shows that the hydrogen volume fraction and ignition position significantly affected the flame structure, flame front speed, and explosion overpressure.     

Experimental study on explosion characteristics of syngas with different ignition positions and hydrogen fraction

The influence of different ignition positions and hydrogen volume fractions on the explosion characteristics of syngas is studied in a rectangular half-open tube. Three ignition positions were set at the axis of the tube, which are 0 mm, 600 mm and 1100 mm away from the closed end, respectively. A range of hydrogen volume fraction (φ) from 10% to 90% were concerned. Experimental results show that different ignition positions and hydrogen volume fraction have important influence on flame propagation structure. When ignited at 600 mm from the closed end on the tube axis, distorted tulip flame forms when flame propagates to the closed end. The formations of the tulip flame and the distorted tulip flame are accompanied by a change in the direction of the flame front propagation. The flame propagation structure and pressure are largely affected by the ignition position and the hydrogen volume fraction. At the same ignition position, flame propagation speed increases with the growing of hydrogen volume fraction. And the pressure oscillates more severe as the ignition location is closer to the open end. And pressure oscillations bring two different forms. The first form is that the pressure has a periodic oscillation. The amplitude of the pressure oscillation gradually increases. It takes several cycles from the start of the oscillation to the peak. For the second form, the pressure reaches the peak of the oscillation in the first cycle of the start to the oscillation.     

Experimental study of premixed syngas/air flame deflagration in a closed duct

The propagation behaviour of a deflagration premixed syngas/air flame over a wide range of equivalence ratios is investigated experimentally in a closed rectangular duct using a high-speed camera and pressure transducer. The syngas hydrogen volume fraction, φ, ranges from 0.1 to 0.9. The flame propagation parameters such as flame structure, propagation time, velocity and overpressure are obtained from the experiment. The effects of the equivalence ratio and hydrogen fraction on flame propagation behaviour are examined. The results indicate that the hydrogen fraction in a syngas mixture greatly influences the flame propagation behaviour. When φ, the hydrogen fraction, is ≥0.5, the prominently distorted tulip flame can be formed in all equivalence ratios, and the minimum propagation time can be obtained at an equivalence ratio of 2.0. When φ < 0.5, the tulip flame distortion only occurs in a hydrogen fraction of φ = 0.3 with an equivalence ratio of 1.5 and above. The minimum flame propagation time can be acquired at an equivalence ratio of 1.5. The distortion occurs when the maximum flame propagation velocity is larger than 31.27 m s^?1. The observable oscillation and stepped rise in the overpressure trajectory indicate that the pressure wave plays an important role in the syngas/air deflagration. The initial tulip distortion time and the plane flame formation time share the same tendency in all equivalence ratios, and the time interval between them is nearly constant, 4.03 ms. This parameter is important for exploring the quantitative theory or models of distorted tulip flames.     

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.     

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.     

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.     

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.     

Numerical study of hydrogen/methane buoyant fires using FireFoam

Hydrogen/methane buoyant fires with various hydrogen volume fractions ranging from 0% to 20% were numerically studied in this paper. The modified eddy dissipation concept combustion model for multi-fuels in the large eddy simulation (LES) framework was employed for combustion, and especially the infinitely fast rate based on “global” concept was improved. Combined with the weighted sum of gray gas model for emission/absorption coefficient, the finite volume discrete ordinates model was used to compute the radiative heat transfer. The predicted centerline temperature, velocity, and flame height are in good consistence with the measured data. Furthermore, the detailed analysis was conducted on the dependency of the parameters such as centerline temperature and velocity, flame height, and soot volume fraction on hydrogen volume concentration.     

Effect of N2 and CO2 on explosion behavior of syngas/air mixtures in a closed duct

The explosion behavior of syngas/air mixtures under the effect of N₂ and CO₂ addition is experimentally investigated in three cases of N₂ and CO₂ volume fractions (0, 20% and 40%). Tests are performed for syngas/air mixtures with varying equivalent ratios (from 0.8 to 2.5) and hydrogen fractions (from 25% to 75%). The effects of N₂ and CO₂ addition on flame structure evolution, flame speed and overpressure histories are analyzed. The results showed that the tulip shaped flames appear in all cases regardless of whether N₂ and CO₂ are added. After flame inversion, the appearance of tulip shaped flame distortion can be observed in syngas/air without N₂ and CO₂ addition and meanwhile the oscillations are seen in the flame front position and speed trajectories. The flame distortion becomes less pronounced with N₂ and CO₂ addition, and the oscillation amplitude of the flame front position and speed reduce accordingly. Both addition of N₂ or CO₂ decrease the flame speed and the maximum overpressure. Therefore, it increases the time required for flame arriving to the discharge vent. Whereas CO₂ has evidently better inhibition performance for syngas/air explosion.     

Laminar burning velocities and flame instabilities of diluted H2/CO/air mixtures under different hydrogen fractions

An experimental study was conducted using outwardly propagating flame to evaluate the laminar burning velocity and flame intrinsic instability of diluted H₂/CO/air mixtures. The laminar burning velocity of H₂/CO/air mixtures diluted with CO₂ and N₂ was measured at lean equivalence ratios with different dilution fractions and hydrogen fractions at 0.1 MPa; two fitting formulas are proposed to express the laminar burning velocity in our experimental scope. The flame instability was evaluated for diluted H₂/CO/air mixtures under different hydrogen fractions at 0.3 MPa and room temperature. As the H₂ fraction in H₂/CO mixtures was more than 50%, the flame became more unstable with the decrease in equivalence ratio; however, the flame became more stable with the decrease in equivalence ratio when the hydrogen fraction was low. The flame instability of 70%H₂/30%CO premixed flames hardly changed with increasing dilution fraction. However, the flames became more stable with increasing dilution fraction for 30%H₂/70%CO premixed flames. The variation in cellular instability was analyzed, and the effects of hydrogen fraction, equivalence ratio, and dilution fraction on diffusive-thermal and hydrodynamic instabilities were discussed.     

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.     

Experimental study on the Influence of Hydrogen Fraction on Self-acceleration of H2/CO/air laminar premixed flame

To investigate self-acceleration propagation characteristics of a laminar premixed flame, an experimental study of H₂/CO/air mixtures with various hydrogen fractions and equivalence ratios was conducted. The acceleration exponent and fractal excess were defined to quantitatively investigated flame self-acceleration in the transition and saturation stages. Also, the influence of flame inherent instabilities on the acceleration exponent in the transition stage were investigated. The results indicate that with an increase in the hydrogen fraction, the first and second critical radius decreased, the proportion of the transition (saturation) stage in the whole flame propagation process decreased (increased), and the acceleration exponent and fractal excess of the transition and saturation stages increased. Because of the limits of flame radius and different degrees of pulsation in the saturation stage, the acceleration exponent and fractal excess at the saturation stage measured do not show obvious regularity; the values are less than 1.5 and 0.33, respectively. When the hydrogen fraction in syngas is changed, the acceleration exponent in the transition stage showed a nonlinear decreasing trend with an increase in the effective Le number. The hydrodynamic instability usually increased with a decrease in flame thickness, and the acceleration exponent in the transition stage increased.     

Flame characteristics of methane/air with hydrogen addition in the micro confined combustion space

This paper investigated methane/air flame characteristics with hydrogen addition in micro confined combustion space experimentally and computationally. The focus is on the effect of hydrogen addition on the methane/air flame stabilization, the onset of flame with repetitive extinction and ignition (FREI), and the global flame quenching in decreasing continuously combustion space. Furthermore, the effects of hydrogen addition on the flame temperature and the local equivalence ratio distribution were analyzed systematically using numerical simulations. In addition, the effects of hydrogen addition on the concentrations of OH and H radicals, and the critical scalar dissipation rate of local flame extinction were discussed. With a higher hydrogen ratio, the mixing is faster, and the flame is smaller. When the micro confined space is narrower, the heat loss to the combustor walls has a higher impact on the flames. The flames with higher hydrogen ratios have therefore lower peak flame temperatures and lower concentrations of H and OH radicals. The results show that hydrogen addition can effectively widen the stable combustion range of methane/air flames in the micro confined space by about 20% when the hydrogen addition ratio reaches 50%. The frequency and the maximum propagation velocity of FREI flames can be increased as well. The quenching distance of methane/hydrogen/air flames decreases nearly linearly with the increase of hydrogen ratio. This is attributed to the higher critical scalar dissipation rate of local flame extinction in flames with a higher hydrogen ratio.     

Experimental study of explosion dynamics of syngas flames in the narrow channel

This article introduced the experimental study of the propagation of a syngas premixed flame in a narrow channel. The structural evolution, flame front position and velocity characteristics of lean and rich premixed flames were investigated at different hydrogen volume fractions as the flame was ignited at the open end of the pipe and propagated to the closed end. The comparative study of the syngas fuel characteristics, flame oscillation frequency and overpressure oscillation frequency prove that the syngas explosion flame oscillation in the narrow passage has a strong coupling relationship with overpressure and fuel heat release rate. The results was shown that the flame structure was strongly influenced by the hydrogen volume fraction of the syngas and the fuel concentration. The distorted tulip flame only appears in lean mixture. At 30% of hydrogen volume fraction, the flame exhibits intense and unstable propagation, manifested as the reciprocating and alternating movement of the flame front. As the volume fraction of hydrogen increased, the velocity of flame propagation and the frequency of oscillation increased. When the hydrogen volume fraction γ ≥ 0.4 at the equivalence ratio of Φ = 0.8, the pressure oscillation amplitude gradually increases and reaching the peak after 200–320 ms. Significantly, when γ = 0.3, the pressure peak increases abnormally. This work can provide support for the safe use of syngas in industry by experimental study of various explosion parameters in the narrow channel.     

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.     

Laminar burning velocity and Markstein length of nitrogen diluted natural gas/hydrogen/air mixtures at normal,reduced and elevated pressures

Flame propagation of premixed nitrogen diluted natural gas/hydrogen/air mixtures was studied in a constant volume combustion bomb under various initial pressures. Laminar burning velocities and Markstein lengths were obtained for the diluted stoichiometric fuel/air mixtures with different hydrogen fractions and diluent ratios under various initial pressures. The results showed that both unstretched flame speed and unstretched burning velocity are reduced with the increase in initial pressure (except when the hydrogen fraction is 80%) as well as diluent ratio. The velocity reduction rate due to diluent addition is determined mainly by hydrogen fraction and diluent ratio, and the effect of initial pressure is negligible. Flame stability was studied by analyzing Markstein length. It was found that the increase of initial pressure and hydrogen fraction decreases flame stability and the flame tends to be more stable with the addition of diluent gas. Generally speaking, Markstein length of a fuel with low hydrogen fraction is more sensitive to the change of initial pressure than that of a one with high hydrogen fraction.     

Effects of radiation and compression on propagating spherical flames of methane/air mixtures near the lean flammability limit

Large discrepancies between the laminar flame speeds and Markstein lengths measured in experiments and those predicted by simulations for ultra-lean methane/air mixtures bring a great concern for kinetic mechanism validation. In order to quantitatively explain these discrepancies, a computational study is performed for propagating spherical flames of lean methane/air mixtures in different spherical chambers using different radiation models. The emphasis is focused on the effects of radiation and compression. It is found that the spherical flame propagation speed is greatly reduced by the coupling between thermal effect (change of flame temperature or unburned gas temperature) and flow effect (inward flow of burned gas) induced by radiation and/or compression. As a result, for methane/air mixtures near the lean flammability limit, the radiation and compression cause large amounts of under-prediction of the laminar flame speeds and Markstein lengths extracted from propagating spherical flames. Since radiation and compression both exist in the experiments on ultra-lean methane/air mixtures reported in the literature, the measured laminar flame speeds and Markstein lengths are much lower than results from simulation and thus cannot be used for kinetic mechanism validation.     

Flame characteristics of premixed H2-air mixtures explosion venting in a spherical container through a duct

The explosion venting duct can effectively reduce the hazard degree of a gas explosion and conduct the venting energy to the safe area. To investigate the flame quantitative propagation law of explosion venting with a duct, the effects of hydrogen fraction and explosion venting duct length on jet flame propagation characteristics of premixed H₂-air mixtures were analyzed through experiment and simulation. The experiment results under initial conditions of room temperature and 1 atm show that when hydrogen fraction was high enough, part of the unburned hydrogen was mixed with air again to reach an ignitable concentration, resulting in the secondary combustion was easier produced and the duration of the secondary flame increased. With the increase of venting duct length, the flame front distance and propagation velocity increased. Meanwhile, the spatial distribution of pressure field and temperature field, and the propagation process and mechanism of the flame venting with a duct were analyzed using FLUENT software. The variation of the pressure wave and the pressure reflection oscillation law in the explosion venting duct was captured. Therefore, in the industrial explosion venting design with a duct, the hazard caused by the coupling of venting pressure and venting flame under different fractions should be considered comprehensively.