Modeling the initial flame acceleration in an obstructed channel using large eddy simulation

The propagation and acceleration of a flame surface past obstructions in a closed square channel was investigated using large eddy simulation. The dynamic Smagorinsky–Lilly subgrid model and the Boger flame surface density combustion model were used. The geometry is essentially two-dimensional with fence-type obstacles distributed on the top and bottom surfaces, equally spaced along the channel length at the channel height. Flame propagation, however, is three dimensional as ignition occurs at a point at the center of the channel cross-section. The effect of obstacle blockage ratio on the development of the flame structure was investigated by varying the obstacle height. Three-dimensional cases were simulated from the initiation of a combustion kernel through spark ignition to the acceleration of the flame front at speeds up to 80 m/s. The transition from laminar flame propagation to turbulent flame propagation within the “thin reaction zone” regime was observed in the simulations. By analyzing the development of the three dimensional flame surface and unburned gas flow field, the formation of several flame structures observed experimentally are explained. Global quantities such as the total flame area and centerline flame velocity were ascertained and compared to the experimental data. High amplitude oscillations in the centerline flame velocity were found to occur from a combination of the unburned gas flow field and fluctuations in the volumetric burning rate.     

Numerical simulation of flame propagation and localized preflame autoignition in enclosures

A novel computational approach based on the coupled 3D Flame-Tracking–Particle (FTP) method is used for numerical simulation of confined explosions caused by preflame autoignition. The Flame-Tracking (FT) technique implies continuous tracing of the mean flame surface and application of the laminar/turbulent flame velocity concepts. The Particle method is based on the joint velocity–scalar probability density function approach for simulating reactive mixture autoignition in the preflame zone. The coupled algorithm is supplemented with the database of tabulated laminar flame velocities as well as with reaction rates of hydrocarbon fuel oxidation in wide ranges of initial temperature, pressure, and equivalence ratio. The main advantage of the FTP method is that it covers both possible modes of premixed combustion, namely, frontal and volumetric. As examples, combustion of premixed hydrogen–air, propane–air, and n-heptane–air mixtures in enclosures of different geometry is considered. At certain conditions, volumetric hot spots ahead of the propagating flame are identified. These hot spots transform to localized exothermic centers giving birth to spontaneous ignition waves traversing the preflame zone at very high apparent velocities, i.e., nearly homogeneous preflame explosion occurs. The abrupt pressure rise results in the formation of shock waves producing high overpressure peaks after reflections from enclosure walls.     

Effect of ignition position on the run-up distance to DDT for hydrogen-air explosions

The method described in this paper enabled reliable and accurate positioning of an overdriven detonation by calculation of shock wave velocities (detonation and retonation) for hydrogen explosions in a closed 18 m long horizontal DN150 pipe. This enabled an empirical correlation between the ignition position and the run-up distance to DDT to be determined. It was shown that the initial ability of the flame to expand unobstructed and the piston-like effect of burnt gas expanding against the closed end of the tube contributed to initial flame acceleration and hence were able to affect the run-up distance to overdriven detonation. Flame speeds and rates of initial pressure rise were also used to explain how these two competing effects were able to produce a minimum in the run-up distance to DDT. The shortest run-up distance to DDT, relative to the ignition position, for this pipe and gas configuration was found when the ignition position was placed 5.6 pipe diameters (or 0.9 m) from the closed pipe end. The shortest run-up distance to DDT relative to the end of the pipe was recorded when the ignition source was placed 4.4 pipe diameters or 0.7 m from the pipe end.     

Numerical simulation of premixed methane–air deflagration in a semi-confined obstructed chamber

In this paper, large eddy simulation coupled with a turbulent flame speed cloure (TFC) subgrid combustion model has been utilized to simulate premixed methane–air deflagration in a semi-confined chamber with three obstacles mounted inside.The computational results are in good agreement with published experimental data, including flame structures, pressure time history and flame speed. The attention is focused on the flame flow field interaction, pressure dynamics, as well as the mechanism of obstacle-induced deflagration. It is found that there is a positive feedback mechanism established between the flame propagation and the flow field. The pressure time history can be divided into four stages and the pseudo-combustion concept is proposed to explain the pressure oscillation phenomenon. The obstacle-induction mechanism includes direct effect and indirect effect, but do not always occur at the same time.     

Numerical study on the spontaneous ignition of pressurized hydrogen during its sudden release into the tube with varying lengths and diameters

Fuel cell vehicles (FCV) and other hydrogen systems with pressurized hydrogen has a safety hazard of spontaneous ignition during its sudden release into the tube. Tube parameter is a key factor affecting the spontaneous ignition of pressurized hydrogen. In this paper, a numerical study on the spontaneous ignition of pressurized hydrogen during its sudden release into the tube with varying lengths and diameters is conducted. The models of Large Eddy Simulation (LES), Eddy Dissipation Concept (EDC), Renormalization Group (RNG), 10-step like opening process of burst disk and 18-step detailed hydrogen combustion mechanism are employed. 6 cases are simulated based on the previous experiments. Numerical results show that the possibility of spontaneous ignition of pressurized hydrogen increases inside the longer and thinner tubes, which agrees with the experimental results. The increasing of tube length has little influence on the shock wave formation and propagation inside the tube. However, there exists critical tube lengths for the generation of Mach disk and the normal shock wave: the maximum and minimum distances for the generation of the Mach disk in 10 mm diameter tube are 7.8 and 6.7 mm, respectively. As for the normal shock wave, these critical values are 22.1 and 19.4 mm, respectively. In addition, the formation times and initial positions of Mach disk and normal shock wave are delayed inside the thicker tube. Due to the shock-affected time increases with the increasing of tube length, the temperature could rise to the critical ignition temperature and triggers the spontaneous ignition due to the sufficient tube length even though the less hydrogen/air mixture and the contact surface with lower temperature is produced inside the thicker tube. Finally, a simple time scale analysis is conducted.     

An experimental study of flame acceleration and deflagration to detonation transition in representative process piping

The paper summarizes the results of experimental tests and accompanying analyses to investigate the factors that govern flame acceleration and potential transition to detonation in a relatively long unobstructed piping system. The overall aim of the work was to obtain sufficient experimental data so as to be able to develop and evaluate methodologies for classifying and predicting potential detonation flame acceleration and deflagration to detonation transition (DDT) hazard in industrial process pipes and mixtures. The present results show that the flame acceleration process in an unobstructed pipe exhibit three distinct phases: an initial establishment phase; a second rapid acceleration phase and a final transition to detonation phase. Test results with ethylene indicate that the acceleration process is not sensitive to initial pressure (all other parameters remaining constant) but can be sensitivity to initial pipe wall temperature or possibly mixture humidity. The presence of bends increases the local rate of turbulent combustion, an effect attributed to the additional turbulence generated downstream of the bend. For straight pipes, detonation was only observed to develop for hydrogen–air and ethylene–air mixtures. Detonation was not observed with methane, propane or acetone as fuel in the present piping apparatus.     

Solution adaptive CFD simulation of premixed flame propagation over various solid obstructions

This paper presents the results of a number of calculations carried out in order to simulate combustion past obstacles of different shape and blockage ratio. The obstacle shapes considered are circles, squares, triangles and flat plates. Two-dimensional simulations are carried out with the McNEWT code. The code solves the reacting flow field with a laminar flamelet model on an unstructured mesh. Adaptive mesh refinement is applied so that the flame front is accompanied by mesh refinement throughout the calculation domain. A transition from laminar to turbulent combustion induced by passage past the obstacle is seen in the simulations. Evidence for the transition is found in the change in flame shape, flame speed and pressure. The simulations are compared with experimental data and there is good agreement between experiment and simulation.     

Modeling of detonation processes in chemically active bubble systems at normal and elevated initial pressures

The detonation processes in chemically active bubble liquids under elevated initial pressures are investigated theoretically. It is shown, that supersonic regimes of wave propagation can exist, if the initial pressure is relatively high and the volume fraction of the bubbles is relatively small. Characteristic values of the bubble detonation wave pressure at sub- and supersonic regimes differ by an order of magnitude.The principal possibility of detonation wave structure transformation in the case of propagation in the mixture with high initial pressure and longitudinal gradient of bubble volume fraction is predicted. The leading shock may transform into a smooth wave of compression.The Chapman-Jouguet conditions for self-sustaining supersonic bubble detonation wave is obtained.A model of shock induced single bubble dynamics and ignition taking into account the real properties of the liquid, inter-phase transition processes, mechanical mixing of phases, ignition delay and continuous shift of chemical equilibrium have been described. Calculations for the oxygen containing bubble in liquid cyclohexane have been performed.     

Flame propagation through aluminum particle cloud in a combustion tube

The propagation of a flame is investigated experimentally and theoretically for a large, horizontal combustion tube containing a mixture of air and aluminum powder with pre-existing turbulence. One end of the tube is closed and the other is connected to a large dump-tank. Twenty dispersion systems are used on the tube to produce a uniform suspension of aluminum dust in the tube with a mean diameter of 6 μm. The characteristics of a flame front from the ignitors at the closed end are measured using photodiodes and the development of pressure is monitored by transducers. Experimental results revealed the entire process of an accelerating flame and the development of shock waves. A set of conservation equations for two-phase turbulent combustion flow is derived, using the two-fluid model, k–ε model, Hinze–Tchen model and EBU-Arrhenius model for turbulent combustion. The SIMPLE scheme usually applied to the homogeneous turbulent combustion is extended to fit this two-phase, reactive behavior. The results of calculations show the positive feedback coupling among combustion, expansion and turbulence during flame propagation. Computed and measured results are generally in good agreement.     

A sneak peek into the phenomenology of fuel mist explosions: The key role of vapor fractions

The modern world depends greatly on hydrocarbons, which are ubiquitous, indispensable fuels used in nearly every existing industry. Although important, their use may trigger dangerous incidents, whether in their production, handling, storage, or transporting phase, especially when aerosolized. In light of proposing a standard procedure to assess the flammability and explosivity of fuel mists, a new test method was established based on the EN 14034 standards series. For the previous purposes, a gravity-fed mist generation system was designed and employed in a modified 20 L explosion vessel. This test method allowed the determination of the ignition sensitivity of several fuels. In addition, their explosion severity was represented by the explosion overpressure P_ex, and the rate of pressure rise dP/dt_ex, two thermo-kinetic parameters determined with a specifically developed control system and custom software. Nonetheless, a noticeable difference in the ignition sensitivity and the explosion severity was perceived when changing suppliers or petroleum cuts of some fuels. Moreover, sensitivity studies showed that both the droplet size distribution and the temperature of the droplets play a significant role in fuel mist explosion. These parameters can be directly related to the vapor fraction surrounding a droplet during its ignition. Consequently, this study focuses on the influence of varying the composition of three well-known and abundantly used fuels. Different petroleum cuts were introduced in different fractions into isooctane, Jet A1 aviation fuel, and diesel fuel mixtures, which were then aerosolized into a uniformly distributed turbulent mist cloud and ignited using spark ignitors of 100 J. Subsequently, complementary tests were executed in a vertical flame propagation tube coupled with a high-speed video camera allowing the visualization of the flame and the determination of the spatial flame velocity, and a tentative estimation of the laminar burning velocity. The latter was also estimated from the pressure-time evolution in the 20 L sphere using existing correlations. Indeed, the determination of the laminar burning velocity can be useful in modeling such accidents. Finally, highlighting the essential role of the mist and vapor fraction during their ignition has led to a better understanding of their explosion mechanisms.     

爆炸形成过程中火焰加速的试验研究

为预防和控制工业爆炸事故,并为脉冲爆轰发动机的研究提供理论指导,分析火焰加速导致的燃烧转爆轰过程的影响因素。采用爆轰管探讨障碍物的阻塞比、混合物的组成、初始压力和点火能等4个因素对爆炸性气体火焰速度和爆轰压力的影响规律。试验结果表明:障碍物的存在能大大提高火焰速度和爆轰压力;爆轰压力随管内障碍物阻塞比的增大先变大后减小,并在阻塞比为0.498,燃料种类为天然气,化学当量比为1时达到最大;爆轰压力还随混合气体初始压力的增大和点火能的提高而增大。选择适宜的条件可大大提高火焰加速速率,促进燃烧向爆轰过程转变。     

Structure and flame speed of dilute and dense layered coal-dust explosions

Multidimensional time-dependent simulations were performed to study the interaction of a stoichiometric methane–air detonation with layers of coal dust. The simulations solved equations representing a Eulerian kinetic-theory-based granular multiphase model applicable to dense and dilute particle volume fractions. These equations were solved using a high-order Godunov-based method for compressible fluid dynamics. Two dust layer concentrations were considered: loose with an initial volume fraction of 1%, and dense with an initial volume fraction of 47%. Each layer was simulated with two types of dust: reactive coal and inert ash. Burning of the coal particles results in a coupled complex consisting of an accelerating shock leading a coal-dust flame. The overall structure of the shock–flame complex resembles that of a premixed fast flame with length scales on the order of several meters. The large length scales are direct results of time needed to lift, mix, heat, and autoignite the particle. The flame speeds are large and much larger than the gas-phase velocity. Large spikes of flame speed are characteristic of the 47% case. These spikes and high flame speed are caused by pockets of coal dust autoigniting ahead of the flame. The flame is choked in the 1% case due to the gas-phase products exceeding the sonic velocity with respect to the flame. The 47% case is choked due to attenuation of pressure waves as they propagate through particles. Inert layers of dust substantially reduce the overpressure, impulse, and speed produced by propagating blast wave. The results also show that loose layers of dust are far more dangerous than dense layers. The shock and flame are more strongly coupled for loose layers, propagate at higher velocity, and produce large overpressures and impulses.     

Experimental study on detonation flame penetrating through flame arrester

In-line detonation flame arresters are important safety apparatus to prevent group tank fires caused by the spreading of fire through vapor connection lines. In this study, a DN50 experimental apparatus aimed at the detonation flame penetration characteristics and failure mechanisms in a flame arrester was set up, and a series of experiments were carried out with 6.6% C2H4 and air mixture. Pressure, and velocity of flame penetrating through flame arrester housing and filters were analyzed. Experimental results showed that the attenuation of pressure and velocity was proportional to the thickness of the filters. Two failure modes of the fire-extinguishing process in the flame arrester were captured directly with a high-speed camera. In Mode I, the detonation flame could go straight through the flame arrester filters when the filters were too thin. In Mode II, when the filters were not sufficiently thick, the remained shock wave pressure of detonation flame was still several times of the initial pressure and could rise sharply at the downstream contraction section, resulting in that the flammable gas at the downstream transition section could be compressed and reignited even the flame had been extinguished by filters. These conclusions are helpful to reveal the nature of failure modes of fire-extinguishing process and design flame arresters with high fire-resisting performance by structure improved.     

Flamelet surface density modelling of turbulent deflagrating flames in vented explosions

The performance of two reaction rate models based on the laminar flamelet concept have been examined by calculating the behaviour of turbulent flame deflagration inside a semi-confined explosion tube. The models formulate the mean rate of reaction as a function of a transport equation for the flamelet surface density. The difference in the models is in modelling the source/sink terms of the flamelet surface density transport equation. The models are validated using laser diagnostics of flame deflagration in methane–air flammable mixture. The predictions are compared with experimental results for propagation, pressure history and flame speed. Sensitivity to cross-flow effects are investigated through comparison between two- and three-dimensional calculations. The numerically simulated results show that experimental trends are well reproduced by both models.     

Effect of concentration and ignition position on vented methane–air explosions

Experiments were conducted in a 1 m³ vessel with a top vent to investigate the effect of methane concentration and ignition position on pressure buildup and flame behavior. Three pressure peaks (p₁, p₂, and P_ext) and two types of pressure oscillations (Helmholtz and acoustic oscillations) were observed. The rupture of vent cover results in p₁ that is insensitive to methane concentration and ignition position. Owing to the interaction between acoustic wave and the flame, p₂ forms in the central and top ignition explosions when the methane–air mixture is near–stoichiometric. When the methane–air mixture is centrally ignited, p₂ first increases and then decreases with an increase in the methane concentration. The external explosion-induced P_ext is observed only in the bottom ignition explosions with an amplitude of several kilopascals. Under the current experimental conditions, flame–acoustic interaction leads to the most serious explosions in central ignition tests. Methane concentration and ignition position have little effect on the frequency of Helmholtz and acoustic oscillations; however, the Helmholtz oscillation lasts longer and first decreases and then increases as the methane concentration increases for top ignition cases. The ignition position significantly affects the Taylor instability of the flame front resulting from the Helmholtz oscillation.     

Detonation transition criteria from the interaction of supersonic shock-flame complexes with different shaped obstacles

The present study is an experimental investigation of the last stages of the deflagration-to-detonation transition. A fast flame following a lead shock was generated by passing a detonation wave through a perforated plate. The shock flame complex then interacts with an obstacle of different shape. We study the influence of the obstacle shape on the transition mechanism to a detonation. The obstacles studied are a single round or square obstacle, a flat plate, a C-shaped and an H-shaped obstacle. The experiments were performed in a thin transparent channel permitting high speed schlieren visualization. Stoichiometric propane-oxygen was investigated at sub-atmospheric conditions. For each obstacle configuration, the initial pressure was changed to modify the flame burning velocity and the Mach number of the leading shock. The burning velocity prior to the interaction was measured experimentally from the displacement velocity of the flame in the videos. This required estimating the speed of the gas ahead of the flame. A linear correction to the speed immediately behind the lead shock was applied using the shock change equations and the measured pressure gradient behind the lead shock in order to account for the non-steadiness of the lead shock and viscous losses to the walls. Three main findings were that the obstacle shape had a minimal influence on the critical flame strength required for transition, although obstacles with a forward facing cavity were able to suppress the transition by isolating the re-initiation event inside the cavity. The main transition mechanism for all geometries was the enhancement of the flame burning velocity through the flame interaction with the shock reflected on the obstacle leading to Richtmyer-Meshkov instability. Finally, it was found that the flame burning velocity of the initial flame required for transition was closely approximated by the Chapman-Jouguet burning velocity. Consistent with the visual observations, this supports the view that transition is favored when the flame is in phase with the acoustic waves, and strong internal pressure waves can be amplified.     

障碍物对油气-空气混合气体泄压爆炸火焰传播特性影响

为了研究障碍物对油气泄压爆炸火焰传播特性的影响规律,进行了不同数量障碍物工况下的对比实验,并利用纹影仪和高速摄影仪记录了火焰传播过程,针对障碍物对火焰形态、火焰锋面位置及火焰传播速度的影响规律进行了研究,结果表明:圆柱体障碍物会导致油气泄压爆炸火焰形态产生褶皱和弯曲变形,诱导层流火焰向湍流火焰转变,加速火焰的传播,对油气泄压爆炸火焰的初始传播形态有显著影响;随着障碍物数量的增多,火焰锋面传播距离点火端的最大距离增大,但到达最远距离的时间减少;障碍物能够增强火焰的传播速度,尤其对障碍物下游火焰影响最为显著,随着障碍物数量的增多,火焰传播的最大速度也随之增大,但达到最大火焰传播速度的时间却随之减少;障碍物的存在增大了油气泄压爆炸过程外部爆炸压力,并且随着障碍物数量的增多,外部爆炸压力峰值增长幅度增大。