Effects of gas compressibility on the dynamics of premixed flames in long narrow adiabatic channels

We examine the dynamics of premixed flames in long, narrow, adiabatic channels focusing, in particular, on the effects of gas compressibility on the propagation. Recognising the importance of the boundary conditions, we examine and compare three cases: flame propagation in channels open at both ends, where the pressure must adjust to the ambient pressure at both ends and the expanding gas is allowed to leave the channel freely, and flame propagation in channels that remain closed at one of the two ends, where the burned/unburned gas remains trapped between the flame and one of the two walls. Earlier studies have shown that a flame accelerates when travelling down a narrow channel as a result of the combined effects of wall friction and thermal expansion. In the present work we show that compressibility effects enhance the transition to fast accelerating flames in channels open at both ends and in channels closed at the ignition end. In both situations, the accelerating flames could reach values that, depending on the effective Mach number, are as large as fifty times the laminar flame speed. In contrast, when the channel is closed at the far end, the acceleration is limited and the propagation speed is damped as the flame approaches the far boundary. Moreover, we show that, in channels closed at their ignition end, the flame in sufficiently long channels evolves into a steadily propagating compression-driven flame. The propagation speed of these flames depends exponentially on the constant-volume equilibrium temperature, which is higher than the (constant pressure) adiabatic flame temperature, and is therefore larger than for ordinary isobaric flames. Fast propagating compression waves cannot emerge in channels that remain open at their ignition end because of the reduced pressure forced by the open boundary.     

半开口管道中的氢/空气火焰加速和压力发展过程

本文研究了氢／空气预混火焰在半开口管道中的火焰加速现象和压力发展过程．结果表明,重复布置的障碍物对火焰速度和压力提升产生显著的影响．火焰传播状态随着氢气当量比的变化而发生改变,在氢气当量比约为０．３４时,火焰速度出现第一次跃变;随着氢气当量比进一步提高,火焰速度发生第二次跃变,即由爆燃转为爆轰．发生爆轰时氢气当量比的范围随着阻塞比的不同而发生变化．     

Flame acceleration in long narrow open channels

We study the propagation of premixed flames in long but finite channels, when the mixture is ignited at one end and both ends remain open and exposed to atmospheric pressure. Thermal expansion produces a continuous flow of burned gas directed away from the flame and towards the end of the channel where ignition took place. Owing to viscous drag, the flow is retarded at the walls and accelerated in the center, producing a pressure gradient that pushes the unburned gas ahead of the flame towards the other end of the channel. As a result the flame accelerates when it travels from end to end of the channel. The total travel time depends on the length of the channel and is proportional to γ^?1ln(1 + γ), where γ is the heat release parameter.     

Flame acceleration in a narrow channel with flow compressibility and diverging or converging walls

We study the effects of non-parallel (diverging or converging) channel walls on flame propagation and acceleration in planar and cylindrical narrow channels, closed at the ignition end and open at the other, accounting for thermal expansion in both the zero Mach number and weakly compressible flow limits. For parallel channel walls, previous work has shown that thermal expansion induces an axial flow in the channel, which can significantly increase the propagation speed and acceleration of the flame. In this study, we consider examples of diverging/converging linear walls, although our asymptotic analysis is also valid for curved walls. The slope of the channel walls is chosen so that the magnitude of the thermal-expansion induced flow through the channel obtained for parallel walls is modified at leading-order, thereby influencing the leading-order flame propagation. For zero Mach number flows, the effect of the diverging/converging channel walls is moderate. However, for weakly compressible flows, the non-parallel walls directly affect the rate at which pressure diffuses through the channel, significantly inhibiting flame acceleration for diverging walls, whereas the flame acceleration process is enhanced for converging walls. We consider several values of the compressibility factor and channel wall slopes. We also show that the effect of a cylindrical channel geometry can act to significantly enhance flame acceleration relative to planar channels. The study reveals several new physical insights on how non-parallel channel walls can influence the ability of flames to accelerate by modifying the flow and pressure distribution induced by thermal expansion.     

Premixed flame response to equivalence ratio perturbations

This paper studies the heat-release oscillation response of premixed flames to oscillations in reactant stream fuel/air ratio. Prior analyses have studied this problem in the linear regime and have shown that heat release dynamics are controlled by the superposition of three processes: flame speed, heat of reaction, and flame surface area oscillations. Each contribution has somewhat different dynamics, leading to complex frequency and mean fuel/air ratio dependencies. The present work extends these analyses to include stretch and non quasi-steady effects on the linear flame dynamics, as well as analysis of nonlinearities in flame response characteristics. Because the flame response is controlled by a superposition of multiple processes, each with a highly nonlinear dependence upon fuel/air ratio, the results are quite rich and the key nonlinearity mechanism varies with mean fuel/air ratio, frequency, and amplitude of excitation. In the quasi-steady framework, two key mechanisms leading to heat-release saturation have been identified. The first of these is the flame-kinematic mechanism, previously studied in the context of premixed flame response to flow oscillations and recently highlighted by Birbaud et al. (Combustion and Flame 154 (2008), 356–367). This mechanism arises due to fluctuations in flame position associated with the oscillations in flame speed. The second mechanism is due to the intrinsically nonlinear dependence of flame speed and mixture heat of reaction upon fuel/air ratio oscillations. This second mechanism is particularly dominant at perturbation amplitudes that cause the instantaneous stoichiometry to oscillate between lean and rich values, thereby causing non-monotonic variation of local flame speed and heat of reaction with equivalence ratio.     

Physical mechanism of ultrafast flame acceleration

We explain the physical mechanism of ultrafast flame acceleration in obstructed channels used in modern experiments on detonation triggering. It is demonstrated that delayed burning between the obstacles creates a powerful jetflow, driving the acceleration. This mechanism is much stronger than the classical Shelkin scenario of flame acceleration due to nonslip at the channel walls. The mechanism under study is independent of the Reynolds number, with turbulence playing only a supplementary role. The flame front accelerates exponentially; the analytical formula for the growth rate is obtained. The theory is validated by extensive direct numerical simulations and comparison to previous experiments.     

Violent folding of a flame front in a flame-acoustic resonance

The first direct numerical simulations of violent flame folding because of the flame-acoustic resonance are performed. Flame propagates in a tube from an open end to a closed one. Acoustic amplitude becomes extremely large when the acoustic mode between the flame and the closed tube end comes in resonance with intrinsic flame oscillations. The acoustic oscillations produce an effective acceleration field at the flame front leading to a strong Rayleigh-Taylor instability during every second half period of the oscillations. The Rayleigh-Taylor instability makes the flame front strongly corrugated with elongated jets of heavy fuel mixture penetrating the burnt gas and even with pockets of unburned matter separated from the flame front.     

Numerical study of strongly-nonlinear regimes of steady premixed flame propagation. The effect of thermal gas expansion and finite-front-thickness effects

Steady propagation of premixed flames in straight channels is studied numerically using the on-shell approach. A first numerical algorithm for solving the system of nonlinear integro-differential on-shell equations is presented. It is based on fixed-point iterations and uses simple (Picard) iterations or the Anderson acceleration method that facilitates separation of different solutions. Using these techniques, we scan the parameter space of the problem so as to study various effects governing formation of curved flames. These include the thermal gas expansion and the finite-front-thickness effects, namely flame stretch, curvature, and compression. In particular, flame compression is demonstrated to have a profound influence on the flame, strongly affecting the dependence of its propagation speed on the channel width b. Specifically, the solutions found exhibit a sharp increase of flame speed with channel width. Under weak flame compression, this increase commences at b/λ_c ≈ 2 ～ 3, where λ_c is the cutoff wavelength, but this ratio becomes significantly larger as the flame compression grows. The results obtained are also used to identify limitations of the analytical approach based on the weak-nonlinearity assumption, and to revise the role of noise in flame evolution.     

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

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

Effects of pressure rise rate on laminar flame speed under normal and engine-relevant conditions

Laminar flame speed (LFS) is one of the most important physicochemical properties of a combustible mixture. At normal and elevated temperatures and pressures, LFS can be measured using propagating spherical flames in a closed chamber. LFS is also used in certain turbulent premixed flame modelling for combustion in spark ignition engines. Inside the closed chamber or engine, transient pressure rise occurs during the premixed flame propagation. The effects of pressure rise rate (PRR) on LFS are examined numerically in this study. One-dimensional simulations are conducted for spherical flame propagation in a closed chamber. Detailed chemistry and transport are considered. Different values of PRR at the same temperature and pressure are achieved through changing the spherical chamber size. It is found that the effect of PRR on LFS is negligible under the normal and engine-relevant conditions considered in this study. This observation is then explained through the comparison between the unsteady and convection terms in the energy equation for a premixed flame.     

Exploring a critical diameter for thermo-acoustic instability of downward propagating flames in tubes

Thermo-acoustic oscillations are observed when a flame ignited at open end of a tube propagates towards the closed end due to interaction between unsteady heat release rate fluctuations from flame and acoustic fluctuations. In our past work, it was found that thermo-acoustic instability increases with decreasing diameter from 7.0 cm to 3.0 cm. A recent study in flame propagation in Hele–Shaw cells showed that thermo-acoustic instability is not observed for plate separation less than or equal to 0.4 cm. Thermoacoustic instabilities cannot be observed in very narrow tubes due to excessive damping from the wall. This opens up the possibility of a critical diameter where thermo-acoustic instability would be maximum. In this work we perform flame propagation experiments with diameter of combustion tube in the range 0.5 cm to 3 cm for a fixed length of 70.2 cm. It was found that thermo-acoustic parametric instability begins at lowest laminar burning velocity when the diameter is around 1.0 cm. This diameter is termed as critical diameter. Critical diameter is found to be independent of Lewis number of mixtures. Existence of a critical diameter is thus proved experimentally. Growth rates of primary instability increase with decreasing diameter and show a maximum around the critical diameter and decrease with further decrease in tube diameter. But, growth rates of secondary instability as well as maximum pressure fluctuation amplitude decreases continuously with decreasing diameter. Mechanisms responsible for these observations and existence of a critical diameter are clarified.     

Characterization of symmetric to non-symmetric flamefront transition in slender microchannels

The purpose of the present work is to analyze propagating two-dimensional flames confined in slender semi open channels, where the combustion process takes place towards the closed end. The study focuses on the calculation of the growth rate of the transition from symmetric to non-symmetric flames propagation identified by Jiménez et al. [1].The combustion cell is initially filled with a stoichiometric mixture of fuel and air at standard conditions. Ignition is induced close to the open end of the channel under planar and gaussian profiles in temperature and species mass fractions which activate a sustained combustion process. The gases inside the chamber, initially stagnant, are accelerated due to the heat generated in the chemical reactions, leading to the development of lateral boundary layers, so that the hot gases exit the channel following a Poiseuille velocity profile. This transverse flow velocity differences are accommodated by means of a symmetric tulip shape formed after a short initial transient.Acoustic waves generated in the ignition process, keep travelling along the channel, bouncing at the walls and interacting with the flame during all the combustion process. Additionally, the flame structure, curved by Darrieus-Landau instability, interacts with the pressure waves triggering small amplitude oscillations (primary oscillation mode), which under certain conditions can transition to higher amplitude oscillations (secondary mode).This transition is observed to be highly dependent not only on the cell geometry, but also in the initial conditions generated by the ignition procedure.The aim of this work is to improve the understanding of this process, complementing the work of Jiménez et al. [1], and to characterize the effect of the channel width in this transition.     

Numerical simulations of the flow field ahead of an accelerating flame in an obstructed channel

The development of the unburned gas flow field ahead of a flame front in an obstructed channel was investigated using large eddy simulation (LES). The standard Smagorinsky–Lilly and dynamic Smagorinsky–Lilly subgrid models were used in these simulations. The geometry is essentially two-dimensional. The fence-type obstacles were placed on the top and bottom surfaces of a square cross-section channel, equally spaced along the channel length at the channel height. The laminar rollup of a vortex downstream of each obstacle, transition to turbulence, and growth of a recirculation zone between consecutive obstacles were observed in the simulations. By restricting the simulations to the early stages of the flame acceleration and by varying the domain width and domain length, the three-dimensionality of the vortex rollup process was investigated. It was found that initially the rollup process was two-dimensional and unaffected by the domain length and width. As the recirculation zone grew to fill the streamwise gap between obstacles, the length and width of the computational domain started to affect the simulation results. Three-dimensional flow structures formed within the shear layer, which was generated near the obstacle tips, and the core flow was affected by large-scale turbulence. The simulation predictions were compared to experimental schlieren images of the convection of helium tracer. The development of recirculation zones resulted in the formation of contraction and expansion regions near the obstacles, which significantly affected the centerline gas velocity. Oscillations in the centerline unburned gas velocity were found to be the dominate cause for the experimentally observed early flame-tip velocity oscillations. At later simulation times, regular oscillations in the unburned streamwise gas velocity were not observed, which is contrary to the experimental evidence. This suggests that fluctuations in the burning rate might be the source of the late flame-tip velocity oscillations. The effect of the obstacle blockage ratio (BR) on the development of the unburned gas flow field was also investigated by varying the obstacle height. Simulation predictions show favorable agreement with the experimental results and indicate that turbulence production increases with increasing obstacle BR.     

Theoretical and experimental studies on mesoscale flame propagation and extinction

Mesoscale flame propagation and extinction of premixed flames in channels are investigated theoretically and experimentally. Emphasis is placed on the effect of wall heat loss and the wall–flame interaction via heat recirculation. At first, an analytical solution of flame speed in mesoscale channels is obtained. The results showed that channel width, flow velocity, and wall thermal properties have dramatic effects on the flame propagation and lead to multiple flame regimes and extinction limits. With the decrease in channel width, there exist two distinct flame regimes, a fast burning regime and a slow burning regime. The existence of the new flame regime and its extended flammability limit render the classical quenching diameter inapplicable. Furthermore, the results showed that at optimum conditions of flow velocity and wall thermal properties, mesoscale flames can propagate faster than the adiabatic flame. Second, numerical simulation with detailed chemistry demonstrated the existence of multiple flame regimes. The results also showed that there is a non-linear dependence of the flame speed on equivalence ratio. Moreover, it is shown that the Nusselt number has a significant impact on this non-linear dependence. Finally, the non-linear dependence of flame speed on equivalence ratio for both flame regimes is measured using a C₃H₈–air mixture. The results are in good agreement with the theory and numerical simulation.     

Effects of cryogenic temperature on premixed hydrogen/air flame propagation in a closed channel

As a carbon-free fuel, hydrogen has received significant attention recently since it can help enable low-carbon-economy. Hydrogen has very broad flammability range and very low minimum ignition energy, and thereby there are severe safety concerns for hydrogen transportation and utilization. Cryo-compressed hydrogen is popularly used in practice. Therefore, it is necessary to investigate the combustion properties of hydrogen at extremely low or cryogenic temperatures. This study aims to assess and interpret the effects of cryogenic temperature on premixed hydrogen/air flame propagation and acceleration in a thin closed channel. Different initial temperatures ranging from normal temperature (T₀ = 300 K) to cryogenic temperature (T₀ = 100 K) are considered. Both one- and two-dimensional hydrogen/air flames are investigated through transient simulations considering detailed chemistry and transport. It is found that when the initial temperature decreases from T₀ = 300 K to T₀ = 100 K, the expansion ratio and equilibrium pressure both increase substantially while the laminar flame speeds relative to unburned and burned gasses decrease moderately. The one-dimensional flame propagation is determined by laminar flame speed and thereby the combustion duration increases as the initial temperature decreases. However, the opposite trend is found to happen to two-dimensional flame propagation, which is mainly controlled by the flame surface area increase due to the no-slip side wall constraint and flame instability. Based on the change in flame surface area, three stages including the initial acceleration, steady burning and rapid acceleration are identified and investigated. It is demonstrated that the large expansion ratio and high pressure rise at cryogenic temperatures can significantly increase the flame surface area in early stage and promote both Darrieus-Landau instability (hydrodynamic instability) and Rayleigh-Taylor instability in later stage. These two instabilities can substantially increase the flame surface area and thereby accelerate flame propagation in hydrogen/air mixtures at cryogenic temperatures. The present study provides useful insights into the fundamental physics of hydrogen flames at extremely low temperatures, and is closely related to hydrogen safety.     

Edge flames stabilized in a non-premixed microcombustor

The dynamics of an edge flame confined in a non-premixed microcombustor model is studied numerically within the context of a diffusive-thermal model. Fuel and oxidizer, separated upstream by a thin plate, flow through a channel with a prescribed velocity. At the tip of the plate, the fuel and oxidizer mix and, when ignited, an edge flame is sustained at some distance from the plate. The objective in this work is to consider the effects of confinement, differential diffusion, and heat loss on the dynamics of an edge flame in a narrow channel. We consider a wide range of channel widths and allow for changing Lewis numbers, and both adiabatic conditions and heat losses along the channel walls. The results illustrate how the flame shape and standoff distance are affected by the channel width, by mixture composition through variations in Lewis numbers and by heat losses. Conditions for flame stabilization, flame oscillations and flame extinction or blowoff are predicted.     

Controlled detonation initiation in hypersonic flow

Experimental evidence of controlled detonation initiation and propagation in a hypersonic flow of premixed hydrogen-air is presented. This controlled detonation initiation is created in a hypersonic facility capable of producing a Mach 5 flow of hydrogen-air. Flow diagnostics such as high-speed schlieren and OH* chemiluminescence results show that a flame deflagration-to-detonation transition occurs as a combined result of turbulent flame acceleration and shock-focusing. The experimental results define three new distinct regimes in a Mach 5 premixed flow: deflagration-to-detonation transition (DDT), unsteady compressible turbulent flames, and shock-induced combustion. A two-dimensional implicit-LES (ILES) simulation, which solves the compressible, reactive Navier-Stokes equations on an adapting grid is conducted to provide additional insight into the local physical mechanism of detonation transition and propagation.     

Numerical simulations of flame propagation and DDT in obstructed channels filled with hydrogen–air mixture

We study flame acceleration and deflagration-to-detonation transition (DDT) in channels with obstacles using 2D and 3D reactive Navier–Stokes numerical simulations. The energy release rate for the stoichiometric H₂–air mixture is modeled by a one-step Arrhenius kinetics. Computations show that at initial stages, the flame and flow acceleration is caused by thermal expansion of hot combustion products. At later stages, shock–flame interactions, Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instabilities, and flame–vortex interactions in obstacle wakes become responsible for the increase of the flame surface area, the energy-release rate, and, eventually, the shock strength. Computations performed for different channel widths d with the distance between obstacles d and the constant blockage ratio 0.5 reproduce the main regimes observed in experiments: choking flames, quasi-detonations, and detonations. For quasi-detonations, both the initial DDT and succeeding detonation reignitions occur when the Mach stem, created by the reflection of the leading shock from the bottom wall, collides with an obstacle. As the size of the system increases, the time to DDT and the distance to DDT increase linearly with d². We also observe an intermediate regime of fast flame propagation in which local detonations periodically appear behind the leading shock, but do not reach it.     

Effects of fractal gating of potassium channels on neuronal behaviours

The classical model of voltage-gated ion channels assumes that according to a Markov process ion channels switch among a small number of states without memory, but a bunch of experimental papers show that some ion channels exhibit significant memory effects, and this memory effects can take the form of kinetic rate constant that is fractal. Obviously the gating character of ion channels will affect generation and propagation of action potentials, furthermore, affect generation, coding and propagation of neural information. However, there is little previous research on this series of interesting issues. This paper investigates effects of fractal gating of potassium channel subunits switching from closed state to open state on neuronal behaviours. The obtained results show that fractal gating of potassium channel subunits switching from closed state to open state has important effects on neuronal behaviours, increases excitability, rest potential and spiking frequency of the neuronal membrane, and decreases threshold voltage and threshold injected current of the neuronal membrane. So fractal gating of potassium channel subunits switching from closed state to open state can improve the sensitivity of the neuronal membrane, and enlarge the encoded strength of neural information.