A study of ignition and combustion of liquid hydrocarbon droplets in premixed fuel/air mixtures in a rapid compression machine

The combustion of two fuels with disparate reactivity such as natural gas and diesel in internal combustion engines has been demonstrated as a means to increase efficiency, reduce fuel costs and reduce pollutant formation in comparison to traditional diesel or spark-ignited engines. However, dual fuel engines are constrained by the onset of uncontrolled fast combustion (i.e., engine knock) as well as incomplete combustion, which can result in high unburned hydrocarbon emissions. To study the fundamental combustion processes of ignition and flame propagation in dual fuel engines, a new method has been developed to inject single isolated liquid hydrocarbon droplets into premixed methane/air mixtures at elevated temperatures and pressures. An opposed-piston rapid compression machine was used in combination with a newly developed piezoelectric droplet injection system that is capable of injecting single liquid hydrocarbon droplets along the stagnation plane of the combustion chamber. A high-speed Schlieren optical system was used for imaging the combustion process in the chamber. Experiments were conducted by injecting diesel droplet of various diameters (50 µm < d_o < 400 µm), into methane/air mixtures with varying equivalence ratios (0 < ϕ < 1.2) over a range of compressed temperatures (700 K < T_c < 940 K). Multiple autoignition modes was observed in the vicinity of the liquid droplets, which were followed by transition to propagating premixed flames. A computational model was developed with CONVERGE™, which uses a 141 species dual-fuel chemical kinetic mechanism for the gas phase along with a transient, analytical droplet evaporation model to define the boundary conditions at the droplet surface. The simulations capture each of the different ignition modes in the vicinity of the injected spherical diesel droplet, along with bifurcation of the ignition event into a propagating, premixed methane/air flame and a stationary diesel/air diffusion flame.     

Oscillatory cool flame combustion behavior of submillimeter sized n-alkane droplet under near limit conditions

This paper reports simulation results of oscillatory cool flame burning of an isolated, submillimeter sized n-heptane (n-C₇H₁₆) droplet in a selectively ozone (O₃) seeded nitrogen-oxygen (N₂-O₂) environments at atmospheric pressure. An evolutionary one-dimensional droplet combustion code encompassing relevant physics and detailed chemistry was employed to explore the roles of low-temperature chemistry, O₃ seeding, and dynamic flame structure on burning behaviors. For X_O2= 21% and a range of selective ozone seeding, near-quasi-steady cool flame burning is achieved directly (without requiring hot flame initiation and radiative extinction). Under low oxygen index conditions, but with significant O₃ seeding (X_O3 = 5%), a nearly quasi-steady cool flame is initially established that then transitions to a dynamically oscillating cool flame burning mode which continues until the droplet is completely consumed. It is found that the oscillation occurs as result of a initial depletion of fuel vapor-oxidizer layer evolving near the droplet surface and its dynamic re-establishment through liquid vaporization and vapor/oxidizer transport. A kinetic analysis indicates that the dynamic competition between the reaction classes- (a) degenerate chain branching and (b) chain termination/propagation - along with continuous fuel and oxygen leakage through the flame location contributes to an oscillatory burning phenomena of ever-increasing amplitude. Analysis based on single full-cycle of oscillatory burning shows that the reaction progression matrices (evolution of heat and species) for QOOH➔chain propagation/termination reactions (here, Q = C₇H₁₄-) directly scales with the gas phase temperature field. On the contrary, the QOOH➔degenerate branching reactions undergoes three distinct stages within the same oscillatory cycle. The coupled flame dynamics and kinetics suggest that in the oscillatory burning mode, kinetic processes dynamically cross through conditions characterizing the negative temperature coefficient (NTC) turnover temperature, separating low temperature and NTC kinetic regimes. In addition, a parametric study is conducted to determine the role of O₃ seeding level on the observed oscillation phenomena.     

Conditional statistics of inert droplet effects on turbulent combustion in reacting mixing layers

Direct numerical simulation (DNS) of turbulent reacting mixing layers laden with evaporating inert droplets is used to assess the droplet effects in the context of the conditional moment closure (CMC) for multiphase turbulent combustion. The temporally developing mixing layer has an initial Reynolds number of 1000 based on the vorticity thickness with more than 16 million Lagrangian droplets traced. An equivalent mixture fraction incorporating the inert vapour mass fractions clearly demonstrates the effects of vapour dilution on the mixture. Instantaneous fields and conditional statistics, such as the singly conditioned scalar dissipation rate, the gas temperature 〈 T _g|η 〉, conditional variance of the gas temperature 〈 T _g ^”2|η 〉 and conditional covariance between the fuel mass fraction and gas temperature 〈 Y _f ^” T _g ^”|η 〉 show considerable droplet effects. Comparison between the droplet-free and droplet-laden reacting mixing layer cases suggests significant extinction in the latter case. The resulting large conditional fluctuations around the conditional means contradict the basic assumption behind the first-order singly conditioned CMC. More sophisticated CMC approaches, such as doubly conditioned or second-order CMCs are, in principle, better able to incorporate the effects of evaporating droplets, but significant modelling challenges exist. The scalar dissipation rate doubly conditioned on the mixture fraction and a normalized gas temperature 〈 χ | η, ζ 〉 exemplifies the modelling complexity in the CMC of multiphase combustion.     

Three stage cool flame droplet burning behavior of n-alkane droplets at elevated pressure conditions: Hot,warm and cool flame

Transient, isolated n-alkane droplet combustion is simulated at elevated pressure for helium-diluent substituted-air mixtures. We report the presence of unique quasi-steady, three-stage burning behavior of large sphero-symmetric n-alkane droplets at these elevated pressures and helium substituted ambient fractions. Upon initiation of reaction, hot-flame diffusive burning of large droplets is initiated that radiatively extinguishes to establish cool flame burning conditions in nitrogen/oxygen “air” at atmospheric and elevated pressures. However, at elevated pressure and moderate helium substitution for nitrogen (X_He?>?20%), the initiated cool flame burning proceeds through two distinct, quasi-steady-state, cool flame burning conditions. The classical “Hot flame” (～1500?K) radiatively extinguishes into a “Warm flame” burning mode at a moderate maximum reaction zone temperature (～ 970?K), followed by a transition to a lower temperature (～765?K), quasi-steady “Cool flame” burning condition. The reaction zone (“flame”) temperatures are associated with distinctly different yields in intermediate reaction products within the reaction zones and surrounding near-field, and the flame-standoff ratios characterizing each burning mode progressively decrease. The presence of all three stages first appears with helium substitution near 20%, and the duration of each stage is observed to be strongly dependent on helium substitutions level between 20–60%. For helium substitution greater than 60%, the hot flame extinction is followed by only the lower temperature cool flame burning mode. In addition to the strong coupling between the diffusive loss of both energy and species and the slowly evolving degenerate branching in the low and negative temperature coefficient (NTC) kinetic regimes, the competition between the low-temperature chain branching and intermediate-temperature chain termination reactions control the “Warm” and “Cool” flame quasi-steady conditions and transitioning dynamics. Experiments onboard the International Space Station with n-dodecane droplets confirm the existence of these combustion characteristics and predictions agree favorably with these observations.     

Appearance of cool flame in flame spread over fuel droplets in microgravity

This research conducted microgravity experiments to investigate phenomena appearing around a droplet existing outside the flame-spread limit. n-Decane droplets are tethered at intersections of SiC fibers. The flame spreads to two- or three-interactive droplets to heat a droplet placed outside the flame-spread limit of the interactive droplets. The cool-flame appearance during the flame spread over droplets was detected using different methods. The droplet diameter was measured with a back illumination to evaluate the vaporization-rate constant and to judge whether the cool flame appears or not. The temperature around the droplet was measured by the thin-filament pyrometry using a near-infrared camera to detect the temperature rise due to cool-flame appearance. The infrared radiation distribution from the combustion products was measured using a mid-wave infrared camera to judge the cool-flame appearance. The results show that a cool flame appears around the droplet existing outside the hot-flame-spread limit and the vaporization completes with the cool flame if the heat input from the hot flame is sufficiently large. This type of flame spread is called hot-to-cool flame spread. The definition of flame spread should be extended considering the cool flame.     

Experimental and numerical investigation of ester droplet combustion: Application to butyl acetate

This paper presents an experimental and numerical study of the combustion of isolated n‑butyl acetate droplets in the standard atmosphere. Numerical simulations are reported using a model that incorporates unsteady gas and liquid transport, variable properties, and radiation. Three skeletal mechanisms of n‑butyl acetate, derived from a large detailed mechanism comprised of 819 species and 52,698 reactions, were used in the numerical simulations to evaluate the influence of the kinetic mechanism on burning. The reduced mechanisms comprised 212 species and 5413 reactions, 157 species and 3089 reactions, and 105 species and 1035 reactions. The numerical model did not include soot formation, though qualitatively mild sooting was noted only for droplets larger than 0.7 mm. The numerical predictions were in good agreement with experimental measurements of droplet and flame diameters. Flame extinction was numerically predicted which was attributed to a decrease of the characteristic diffusion time relative to the chemical time as droplet burned. Effects of initial droplet diameter on the evolution of maximum gas temperature (T_max) and peak mole fractions of CO₂ and CO are also examined numerically.     

The role of cool-flame chemistry in quasi-steady combustion and extinction of n-heptane droplets

Experiments on the combustion of large n-heptane droplets, performed by the National Aeronautics and Space Administration in the International Space Station, revealed a second stage of continued quasi-steady burning, supported by low-temperature chemistry, that follows radiative extinction of the first stage of burning, which is supported by normal hot-flame chemistry. The second stage of combustion experienced diffusive extinction, after which a large vapour cloud was observed to form around the droplet. In the present work, a 770-step reduced chemical-kinetic mechanism and a new 62-step skeletal chemical-kinetic mechanism, developed as an extension of an earlier 56-step mechanism, are employed to calculate the droplet burning rates, flame structures, and extinction diameters for this cool-flame regime. The calculations are performed for quasi-steady burning with the mixture fraction as the independent variable, which is then related to the physical variables of droplet combustion. The predictions with the new mechanism, which agree well with measured autoignition times, reveal that, in decreasing order of abundance, H₂O, CO, H₂O₂, CH₂O, and C₂H₄ are the principal reaction products during the low-temperature stage and that, during this stage, there is substantial leakage of n-heptane and O₂ through the flame, and very little production of CO₂ with no soot in the mechanism. The fuel leakage has been suggested to be the source of the observed vapour cloud that forms after flame extinction. While the new skeletal chemical-kinetic mechanism facilitates understanding of the chemical kinetics and predicts ignition times well, its predicted droplet diameters at extinction are appreciably larger than observed experimentally, but predictions with the 770-step reduced chemical-kinetic mechanism are in reasonably good agreement with experiment. The computations show how the key ketohydroperoxide compounds control the diffusion-flame structure and its extinction.     

Vaporization and combustion in three-dimensional droplet arrays

The quasi-steady vaporization and combustion of multiple-droplet arrays is studied numerically. Utilizing the Shvab–Zeldovich formulation, a transformation of the governing equations to a three-dimensional Laplace’s equation is performed, and the solution to Laplace’s equation is obtained numerically to find the effects of droplet interactions in symmetric, multiple-droplet arrays. Vaporization rates, flame surface shapes, and flame locations are found for different droplet array configurations and fuels. The number of droplets, the droplet arrangement within the arrays, and the droplet spacing within the arrays are varied to determine the effects of these parameters. Computations are performed for uniformly spaced three-dimensional arrays of up to 216 droplets, with center-to-center spacing ranging from 3 to 25 droplet radii. As a result of the droplet interactions, the number of droplets and relative droplet spacing significantly affect the vaporization rate of individual droplets within the array, and consequently the flame shape and location. For small droplet spacing, the individual droplet vaporization rate decreases below that obtained for an isolated droplet by several orders of magnitude. A similarity parameter which correlates vaporization rates with array size and spacing is identified. Individual droplet flames, internal group combustion, and external group combustion can be observed depending on the droplet geometry and boundary conditions.     

Measurements and modelling of oxy-fuel coal combustion

Oxy-fuel combustion is one of the most promising technologies to isolate efficiently and economically CO₂ emissions in coal combustion for the ready carbon sequestration. The high proportions of both H₂O and CO₂ in the furnace have complex impacts on flame characteristics (ignition, burnout, and heat transfer), pollutant emissions (NO_x, SO_x, and particulate matter), and operational concerns (ash deposition, fouling/slagging). In contrast to the existing literature, this review focuses on fundamental studies on both diagnostics and modelling aspects of bench- or lab-scale oxy-fuel combustion and, particularly, gives attention to the correlations among combustion characteristics, pollutant formation, and operational ash concerns. First, the influences of temperature and species concentrations (e.g., O₂, H₂O) on coal ignition, volatile combustion and char burning processes, for air- and oxy-firing, are comparatively evaluated and modelled, on the basis of data from optically-accessible set-ups including flat-flame burner, drop-tube furnace, and down-fired furnace. Then, the correlations of combustion-generated particulate/NO_x emissions with changes of combustion characteristics in both air and oxy-fuel firing modes are summarized. Additionally, ash deposition propensity, as well as its relation to the formation of fine particulates (i.e. PM_0.2, PM₁ and PM₁₀), for both modes are overviewed. Finally, future research topics are discussed. Fundamental oxy-fuel combustion research may provide an ideal alternative for validating CFD simulations toward industrial applications.     

Influence of nitrogen in aluminum droplet combustion

The influence of nitrogen on the aluminum droplet combustion under forced convection conditions has been studied. An aerodynamic levitation technique of millimetric size liquid droplets heated with a CO₂ laser has been adopted to characterize the combustion of aluminum droplets and, in particular, to observe the surface phenomena. The determination of the burning rate and of the droplet temperature in several atmospheres (H₂O/O₂, H₂O/Ar, H₂O/N₂, and air) has shown that they depend only on the nature and concentration of the oxidizers (O₂ and H₂O); a comparison of experiments in nitrogen and in argon containing mixtures demonstrated that N₂ did not influence the gas phase combustion. However, for nitrogen containing atmospheres we observed the formation of solid aluminum nitride (AlN) at the droplet surface after a latency time depending on the nitrogen pressure. AlN first interacts with the oxide cap producing an aluminum oxynitride, then completely covers the droplet, and finally prevents combustion. The existence of a latency time varying with the nitrogen pressure suggests that the AlN formation is controlled by heterogeneous kinetics. The phenomenon of oxide cap regression during combustion was also observed in all gases, and it is attributed to a chemical decomposition process of alumina by aluminum forming gaseous Al_xO_y species. Therefore, nitrogen effects are significant at the droplet surface rather than in the gas phase, and it is suggested that N₂ is probably one of the main species causing the manifestation of unsteady processes during aluminum droplet burning.     

Extinction characteristics of isolated n-alkane fuel droplets during low temperature cool flame burning in air

Large carbon number n-alkanes are a notable component in all real transportation fuels, and their chemical structure fosters substantial low temperature kinetic reactivity. Normal alkanes have been studied in various canonical configurations but rarely in systems with strong coupling between low temperature chemistry and transport for pure as well as for multi-component n-alkane mixtures. The Flame Extinguishment (FLEX) experiments onboard the International Space Station provided a unique platform for investigating low temperature multi-phase n-alkane and iso-alkane combustion. Among the many interesting phenomena experimentally observed, cool flame extinction can occur, accompanied by the concurrent formation of a surrounding cloud of condensed vapor. In this work we conduct numerical simulations of high and low temperature combustion of large, initially single-component n-heptane, n-decane and n-dodecane droplets. The role of initial droplet diameter, operating pressure, and n-alkyl carbon number on the extinction of hot and low temperature flames is investigated and compared against the available experimental data. While all three fuels exhibit similar hot flame behavior, cool flame activity increases with the carbon number, resulting in an increased cool flame temperature and decreased extinction diameter. Multi-cyclic “hot/cool flame transitions” are found in air as pressure is slightly increased above one atmosphere. The cyclic behaviors correspond to continuously varying hot and cool flame transitions across the high, low, and negative temperature coefficient (NTC) kinetic regimes. Further increase in pressure results in a second stage steady “Warm flame” transition. The extinction of hot and cool flame has a strong non-linear dependence on ambient pressure but as the hot flame extinction diameter increases with pressure the extinction diameter of the cool flame decreases. The computational results are compared with a recent asymptotic analysis of FLEX n-alkane cool flames.     

Droplet combustion in standing sound waves

Interaction between droplet combustion and acoustic oscillation is clarified. As the simplest model, an isolated fuel droplet is combusted in a standing sound wave. Apart from the conventional idea that oscillatory component of flow influences heat and mass transfer and promotes combustion, a new model that a secondary flow dominates combustion promotion is examined. The secondary flow, found by the authors in the previous work, is driven by acoustic radiation force due to Reynolds normal stress, and named as thermo-acoustic streaming. Since the force is described by the same equation as buoyancy, i.e., F = ΔρVg, the nature of the streaming is thought to be the same as natural convection. The flow patterns of the streaming are analyzed and its influence on burning rate of a droplet is predicted. Experimental investigation was mainly done with burning droplets located in the middle of node and anti-node of standing sound waves. This location realizes the strongest streaming. By varying sound pressure level, ambient pressure, and acoustic frequency, the strength of the streaming was controlled. Flame configuration including soot and burning rate were examined. Microgravity conditions were employed to clarify the influence of acoustic field through the streaming, since it is similar to and must be distinguished from natural convection. Experiments using microgravity conditions confirmed the new combustion promotion model and the way to quantify it. By introducing a new non-dimensional number Gr_a, that is the ratio of acoustic radiation force to viscosity, burning rate constants for various ambient and sound conditions are rearranged. As a result, it was found that the excess burning rate (k/k₀ − 1) is proportional to or , for weak sound and for strong sound, respectively.     

Low field63Cu NMR study of indirect nuclear spin-spin coupling in YBa2Cu3O6.9. A new approach to probe strong correlation effects in high-T c superconductors

The authors carried out a low field (6.3 kOe) NMR measurement of the Gaussian component of⁶³Cu nuclear spin-spin relaxation rate 1/T _2G at the planar Cu site in a high-T _c superconductor YBa₂Cu₃O_6.9 (T _c=92 K). They demonstrate that the results provide quantitative information concerning the static spin susceptibility χ′(q) at non-zero wave vector q. The detailed analysis of the data assuming a model susceptibility peaked at the corner of first Brillouin zone q=Q (Q=(π/a, θ/a),a: lattice constant) shows that χ′(q) satisfies a Curie-Weiss law around q=Q. Stoner enhancement over the calculated Lindhard function at q=Q is estimated to be of the order of ～10–20. They also demonstrate that combined analysis of 1/T ₁ and 1/T _2G allows one to judge whether anomalous shift of the low frequency spectral weight of χ″(Q,ω) to higher frequencies (i.e.pseudo gap) exists or not. Application of the method to YBa₂Cu₃O_6.9 revealed no appreciable pseudo gap in the material.     

Plasma supported combustion

Oxidation of molecular hydrogen and different hydrocarbons in stoichiometric mixtures with air and oxygen in the pulsed nanosecond discharges was studied at room temperature, and the detailed kinetics of the process has been numerically investigated. In the discharge afterglow, the reactions including electron-excited particles play a dominant role for the time up to 100 ns, ion–molecular reactions—for the time of microsecond range, and reactions including radicals mostly contribute for the time interval of several milliseconds. The principal role of processes with formation of excited components that support the development of the chain mechanism of oxidation has been shown. The spatial uniformity of the gas-mixture combustion initiated by a high-voltage nanosecond volume discharge is investigated at gas pressures of 0.3–2.4 atm and temperatures of 1000–2250 K. The self-ignition time and the time of discharge-induced ignition are determined. It is found that the discharge significantly (by 600 K) decreases the ignition temperature with very low energy in the discharge (10⁻² J/cm³). The influence of gas excitation by a pulsed nanosecond discharge with a high-voltage pulse amplitude up to 25 kV on the properties of a premixed propane–air flame has been investigated over a wide range of the equivalence ratios (0.4–5). It was experimentally found that the flame’s blow-off velocity increased more than twice at a discharge energy input less than 1% of the burner power. Efficient production of active radicals under the action of a barrier discharge has been observed. The increase in the flame’s propagation velocity is explained by the production of atomic oxygen in a discharge by the quenching of electronically excited molecular nitrogen N₂ and the dissociation of molecular oxygen on electron-impact. A numerical model has been developed, which describes the influence of pulsed electric discharges on the ignition, combustion, and flame propagation.     

Kinetics and extinction of non-premixed cool and warm flames of dimethyl ether at elevated pressure

The growing demand of clean and efficient propulsion and energy systems has sparked an interest in understanding low-temperature combustion at high pressure. Cool flame transition and extinction limits as well as oxygen concentration dependence at elevated pressures provide insights of the low-temperature and high-pressure fuel reactivity. A new experimental high-pressure counterflow burner platform was designed and developed to achieve the studies of high-pressure cool flames. Dimethyl ether (DME) was chosen to study its non-premixed cool flame in high-pressure counterflow burner at pressure up to 5 atm, perhaps for the first time. This paper investigates the effects of pressure on cool flame structure, extinction and transition limits, and oxygen concentration dependence as well as ozone assisted warm flames of DME in experiments and numerical simulations. The results show that the reignition transition from cool flame to hot flame occurs either with the decrease of the strain rate at a given fuel concentration and pressure or with the increase of fuel mole fraction or pressure at a given strain rate. Furthermore, it is shown that the higher pressure shifts the cool flame to higher strain rates and results in higher cool flame extinction strain rates. However, the existing kinetic model of DME fails in predicting the cool flame extinction limit at elevated pressures. Besides, the cool flame extinction limits are proportional to nth power of the oxygen concentration, [O₂]ⁿ, and the increase of pressure leads to stronger extinction limit dependence (larger n) on oxygen concentration. The present experiment and detailed kinetic analysis show clearly that increasing pressure promotes the low-temperature chemistry including the oxygen addition reactions. In addition, stable warm flame was first experimentally observed by using DME at elevated pressure with ozone sensitization.     

An analytical study of droplet combustion under microgravity: Quasi-steady transient approach

An analytical model based on an assumption of combined quasi-steady and transient behavior of the process is presented to exemplify the unsteady, sphero-symmetric single droplet combustion under microgravity. The model used in the present study includes an alternative approach of describing the droplet combustion as a process where the diffusion of fuel vapor residing inside the region between the droplet surface and the flame interface experiences quasi-steadiness while the diffusion of oxidizer inside the region between the flame interface and the ambient surrounding experiences unsteadiness. The modeling approach especially focuses on predicting; the variations of droplet and flame diameters with burning time, the effect of vaporization enthalpy on burning behavior, the average burning rates and the effect of change in ambient oxygen concentration on flame structure. The modeling results are compared with a wide range of experimental data available in the literature. It is shown that this simplified quasi-steady transient approach towards droplet combustion yields behavior similar to the classical droplet theory.     

On the influence of singlet oxygen molecules on characteristics of HCCI combustion: A numerical study

Mechanisms of homogeneous charge compression ignition (HCCI) combustion enhancement are investigated numerically when excited O₂(a ¹Δ_g) molecules are produced at different points in the compression stroke. The analysis is conducted with the use of an extended kinetic model involving the submechanism of nitric oxide formation in the presence of singlet oxygen O₂(a ¹Δ_g) or O₂(b ¹Σ_g ⁺) molecules in the methane-air mixture. It is demonstrated that the abundance of excited O₂(a ¹Δ_g) molecules in the mixture even in a small amounts intensifies the ignition and combustion and allows one to control the ignition event in the HCCI engine. Such a method of energy supply in the HCCI engine is much more effective in advancement of combustion timing than mere heating of the mixture, because it leads to acceleration of the chain-branching mechanism. The excitation of O₂ molecules to the a ¹Δ_g electronic state makes it possible to organise the successful combustion in the cylinder at diminished initial temperature of the mixture and increase the effective energy released during HCCI combustion. The advance in the value of this energy is much higher than the energy needed for the excitation of oxygen molecules. Moreover, in this case, the output concentration of NO and CO can be reduced significantly.     

Numerical study on flame spread of an n-decane droplet array in different temperature environment under microgravity

A series of numerical calculations of flame spread of an n-decane droplet array was conducted at different ambient temperatures (T_a = 300 and 573 K) for S/d₀ from 1.5 to 10, where S is the droplet interval and d₀ is the initial droplet diameter. The authors compared these numerical results with experimental results under similar conditions at different ambient temperatures for the first time in this study. Good qualitative agreement in flame spread behavior between numerical results and microgravity experiments is obtained. Flame spread mode changed with an increase in S/d₀. Also, appearance of the flame spread mode in a stepping-stone manner (Mode III in [Jpn. Soc. Mech. Eng. 68 (672) (2002) 2423]) in a normal temperature environment was verified by numerical calculations and microgravity experiments, although it was not predicted in the theoretical analysis. In addition, good qualitative agreement of flame spread rate V_f versus S/d₀ was obtained between numerical and experimental results, although numerical results were at least twice as large as experimental results. V_f had a maximum peak at a specific S/d₀ for a different ambient temperature. Employment of improved reaction model and consideration for thermal radiation heat transfer are expected to produce quantitatively better results. An increase in surface temperature of unburned droplets and the development of a flammable gas layer around the droplets were promoted in a high-temperature environment, due to an increase in heat transfer from ambient air to the droplet. As a result, V_f was increased by the higher ambient temperature, suggesting that ambient temperature plays a significant role both in the flame spread mode and the flame spread rate through promotion of a flammable gas layer around unburned droplets.     

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.