Flame structure and stabilization of lean-premixed sprays in a counterflow with low-volatility fuel

An experimental study was performed on the combustion of lean-premixed spays in a counterflow. n-Decane was used as a liquid fuel with low volatility. The flame structure and stabilization were discussed based on the flame-spread mechanism of a droplet array with a low-volatility fuel. The spray flame consisted of a blue region and a yellow luminous region. The flame spread among droplets and group-flame formation through the droplet interaction were observed on the premixed spray side, while envelope flames were also observed on the opposing airflow side. The blue-flame region consisted of premixed flames propagating in the mixture layer around each droplet, the envelope diffusion flames around each droplet, the lower parts of the group diffusion flame surrounding each droplet cluster, and the envelope flame around droplets passing through the group flame. The flame was stabilized within a specific range of the mean droplet diameter via a balance between the droplet velocity and the flame-spread rate of the premixed spray.     

Lift-off characteristics of triple flame with concentration gradient

The concentration gradient and uniform mean velocity of a triple flame in a mixing layer were studied using a multi-slot burner, which can stabilize the lift-off flame especially at a very small concentration gradient. Flame stabilization conditions were examined, and the lift-off heights of the triple flame were measured for methane and propane flames. A hot-wire anemometer was used to measure the velocity distributions. Mass spectroscopy (for methane) and Rayleigh scattering (for propane) were used to measure the concentration gradients. OH radical distribution was measured by laser-induced fluorescence (LIF), and in-stream velocity variation was measured with particle-image velocimetry (PIV). Maximum in-stream temperatures were measured using the coherent anti-Stokes Raman scattering (CARS) technique. Lift-off heights of triple flames have minimum values during the increase of concentration gradient, and the propagation velocity of triple flames reaches its maximum at a critical concentration gradient. This is caused by three factors: velocity distribution upstream, flammable limit of premixed gas, and reaction of diffusion flame. The critical concentration gradient, which maximizes the propagation velocity is suggested as a new criterion of transition from a premixed flame to a triple flame.     

Flame structure and flame spread rate over a solid fuel in partially premixed atmospheres

We have investigated the downward flame spread over a thin solid fuel. Hydrogen, methane, or propane, included in the gaseous product of pyrolysis reaction, is added in the ambient air. The fuel concentration is kept below the lean flammability limit to observe the partially premixing effect. Both experimental and numerical studies have been conducted. Results show that, in partially premixed atmospheres, both blue flame and luminous flame regions are enlarged, and the flame spread rate is increased. Based on the flame index, a so-called triple flame is observed. The heat release rate ahead of the original diffusion flame is increased by adding the fuel, and its profile is moved upstream. Here, we focus on the heat input by adding the fuel in the opposed air, which could be a direct factor to intensify the combustion reaction. The dependence of the flame spread rate on the heat input is almost the same for methane and propane/air mixtures, but larger effect is observed for hydrogen/air mixture. Since the deficient reactant in lean mixture is fuel, the larger effect of hydrogen could be explained based on the Lewis number consideration. That is, the combustion is surely intensified for all cases, but this effect is larger for lean hydrogen/air mixture (Le < 1), because more fuel diffuses toward the lean premixed flame ahead of the original diffusion flame. Resultantly, the pyrolysis reaction is promoted to support the higher flame spread rate.     

Flame stabilization in high pressure LOx/GH2 and GCH4 combustion

Planar laser induced fluorescence (PLIF) of OH is used to examine flame stabilization in high pressure cryogenic flames formed by injecting a central jet of low speed liquid oxygen surrounded by a high speed gaseous stream of hydrogen or methane. In the LO_x/GH₂ experiments injection conditions are transcritical as the chamber pressure is above critical but the temperature is below critical . In the LO_x/GCH₄ experiments the chamber pressure and LO_x injection temperature are below critical , . Hydrogen or methane are injected at room temperature LIF images delineate the flame edge in the injector nearfield. The two flames are stabilized in the vicinity of the liquid oxygen injector lip but the anchor point is found to lie closer to the lip in the LO_x/GH₂ case and its displacement from shot to shot is of a smaller amplitude than that corresponding to the LO_x/GCH₄ flame. Interpretation of these data is based on a previous analysis which indicates that stabilization is essentially controlled by a dimensionless group formed by comparing the lip thickness to the flame edge thickness Ψ = h_s/δ_f. It is found that Ψ slightly exceeds unity in the LO_x/GH₂ case essentially fulfilling the stability condition while Ψ < 1 in the LO_x/GCH₄ case. In this last situation the flame is thicker than the characteristic thickness h_s and it is therefore sensitive to the high speed methane stream. Anchoring is imperfect and the flame edge moves with the turbulent eddies shed from the lip. Global stabilization is achieved dynamically but the reactive layer is not well established and the large amplitude motion of the edge is a symptom of a possible lift-off. Theoretical estimates indicate that LO_x/GCH₄ flame stabilization requires a thicker lip size than the LO_x/GH₂ propellant couple.     

Flame dynamics

This lecture describes recent theoretical developments associated with the dynamics of flames, obtained primarily by exploiting the various temporal and length scales involved in the combustion process. In premixed flames the focus is on flame-flow interactions that occur during the nonlinear development of hydrodynamically unstable large-scale flames, or during the propagation of curved flames in two-dimensional channels. The second part of the paper deals with non-premixed and partially premixed flames, where the focus is on understanding the nature of diffusive-thermal instabilities including the effect of thermal expansion, and on stabilization mechanisms of edge flames, which possess characteristics of both premixed and diffusion flames. The results presented in this talk illustrate how simplified models, when analyzed to their extreme, yield predictions of qualitative nature with physical insight that have advanced our understanding of combustion. This insight can be used to guide the experimental efforts, explain observations and validate large-scale numerical simulations.     

A numerical study on the formation of diffusion flame islands in a turbulent hydrogen jet lifted flame

This paper presents a numerical study on the formation of diffusion flame islands in a hydrogen jet lifted flame. A real size hydrogen jet lifted flame is numerically simulated by the DNS approach over a period of about 0.5 ms. The diameter of hydrogen injector is 2 mm, and the injection velocity is 680 m/s. The lifted flame is composed of a stable leading edge flame, a vigorously turbulent inner rich premixed flame, and a number of outer diffusion flame islands. The relatively long-term observation makes it possible to understand in detail the time-dependent flame behavior in rather large time scales, which are as large as the time scale of the leading edge flame unsteadiness. From the observation, the following three findings are obtained concerning the formation of diffusion flame islands. (1) A thin oxygen diffusion layer is developed along the outer boundary of the lifted flame, where the diffusion flame islands burn in a rather flat shape. (2) When a diffusion flame island comes into contact with the fluctuating inner rich premixed flame, combustion is intensified due to an increase in the hydrogen supply by molecular diffusion. This process also works for the production of the diffusion flame islands in the oxygen diffusion layer. (3) When a large unburned gas volume penetrates into the leading edge flame, the structure of the leading edge flame changes. In this transformation process, a diffusion flame island comes near the leading edge flame. The local deficiency of oxygen plays an important role in this production process.     

Mechanistic Study of the Flame Retardancy of Epoxy Resin With a Novel Phosphorus and Silicon-Containing Flame Retardant

Abstract

An epoxy resin (EP) composite with a novel phosphorus and silicon-containing flame retardant (9,10-dihydro-9-oxa-10-phosphaphenanthrene-10-oxide-polyvinylsilicone-polyphenylaminosilicone, DOPO-V-PA) was prepared in this study. From the cone calorimeter measurements, it was confirmed that the EP/DOPO-V-PA (FREP) composite had relatively better flame retardancy than pure EP. The activation energy was calculated with the Kissinger and Ozawa–Flynn–Wall methods. The results showed that the activation energies of pure EP had slightly higher values than the FREP composite in the early and middle degradation stage (conversion ≤ 90% spaces), which indicated that the earlier degradation of the DOPO-V-PA at low temperature accelerated the degradation of the EP matrix. However, the activation energy values of the FREP were higher than those of pure EP in the final stage, which are attributed to the thermal stabilization effect of the DOPO-V-PA through the promoting of EP char formation.     

LES/CMC modelling of ignition and flame propagation in a non-premixed methane jet

The Large Eddy Simulation (LES) / Conditional Moment Closure (CMC) model with detailed chemistry is used for modelling spark ignition and flame propagation in a turbulent methane jet in ambient air. Two centerline and one off-axis ignition locations are simulated. We focus on predicting the flame kernel formation, flame edge propagation and stabilization. The current LES/CMC computations capture the three stages reasonably well compared to available experimental data. Regarding the formation of flame kernel, it is found that the convection dominates the propagation of its downstream edge. The simulated initial downstream and radial flame propagation compare well with OH-PLIF images from the experiment. Additionally, when the spark is deposited at off-centerline locations, the flame first propagates downstream and then back upstream from the other side of the stoichiometric iso-surface. At the leading edge location, the chemical source term is larger than others in magnitude, indicating its role in the flame propagation. The time evolution of flame edge position and the final lift-off height are compared with measurements and generally good agreement is observed. The conditional quantities at the stabilization point reflect a balance between chemistry and micro-mixing. This investigation, which focused on model validation for various stages of spark ignition of a turbulent lifted jet flame through comparison with measurements, demonstrates that turbulent edge flame propagation in non-premixed systems can be reasonably well captured by LES/CMC.     

Simulation of flame lift-off on a diesel jet using a generalized flame surface density modeling approach

A generalized flame surface density modelling approach is presented to simulate the transient ignition and flame stabilization of a diesel jet flame, for which experimental data are available. The approach consists of four submodels: a mixing model, a generalized flame surface density model, a generalized progress variable model, and a chemistry model. A database containing the laminar model reaction rates per unit generalized flame surface density is generated by solving the unsteady flamelet equations. The RANS-CFD code solves for the mean flame surface density and mean progress variable. The coupling of the models is done via the progress variable and the scalar dissipation rate. The proposed approach is found to be adapted to simulate such a lifted flame and yields good trend agreement for ignition delay and flame lift-off vs. liquid penetration. These first promising results are encouraging to further explore and to apply this method to a more industrial configuration such as a diesel engine.     

Dynamic behavior of a freely-propagating turbulent premixed flame under global stretch-rate oscillations

Dynamic features of a freely propagating turbulent premixed flame under global stretch rate oscillations were investigated by utilizing a jet-type low-swirl burner equipped with a high-speed valve on the swirl jet line. The bulk flow velocity, equivalence ratio and the nominal mean swirl number were 5 m/s, 0.80 and 1.23, respectively. Seven velocity forcing amplitudes, from 0.09 to 0.55, were examined with a single forcing frequency of 50 Hz. Three kinds of optical measurements, OH-PLIF, OH^* chemiluminescence and PIV, were conducted. All the data were measured or post-processed in a phase-locked manner to obtain phase-resolved information. The global transverse stretch rate showed in-phase oscillations centering around 60 (1/s). The oscillation amplitude of the stretch rate grew with the increment of the forcing amplitude. The turbulent flame structure in the core flow region varied largely in axial direction in response to the flowfield oscillations. The flame brush thickness and the flame surface area oscillated with a phase shift to the stretch rate oscillations. These two properties showed a maximum and minimum values in the increasing and decreasing stretch periods, respectively, for all the forcing amplitudes. Despite large variations in flame brush thickness at different phase angles, the normalized profiles collapse onto a consistent curve. This suggests that the self-similarity sustains in this dynamic flame. The global OH^* fluctuation response (i.e. response of global heat-release rate fluctuation) showed a linear dependency to the forcing velocity oscillation amplitudes. The flame surface area fluctuation response showed a linear tendency as well with a slope similar to that of the global OH^* fluctuation. This indicated that the flame surface area variations play a critical role in the global flame response.     

Premixed flame oscillations in obstructed channels with both ends open

An initially laminar premixed flame front accelerates extremely fast and may even trigger a detonation when propagating in a semi-open obstructed channel (one end of the channel is closed; the flame is ignited at the closed end and moves towards the open one). However, industrial and laboratory conduits oftentimes have both ends open, or vented, with a flame ignited at one of these ends. The latter constitutes the focus of the present work. Specifically, premixed flame propagation through a comb-shaped array of obstacles, in-built in a channel with both ends open, is studied by means of computational simulation of the reacting flow equations with fully-compressible hydrodynamics and an Arrhenius chemical kinetics. The parametric study includes various blockage ratios and spacing as well as the thermal expansion ratios, with oscillations of the burning rate observed in the majority of the cases, which conceptually differs from fast flame acceleration in semi-open channels. Such a difference is devoted to the fact that while the entire flame-generated jet-flow is pushed towards a single exit in a semi-open channel, in a channel with two ends open, this jet-flow is distributed between the upstream and downstream flows, thereby moderating flame propagation. The flame oscillations are nonlinear in all cases where they are observed. The oscillation period grows with the blockage ratio but decreases with the thermal expansion. The present results also support the recent experiments, modeling and theory of flames in obstructed channels with both ends open, which all yielded steady or quasi-steady flame propagation prior to the onset of flame acceleration. Indeed, the present oscillations can be treated as the fluctuations around a quasi-steady solution.     

Dynamic behavior of diffusion flame interacting with a large-scale vortex by laser imaging techniques

Instantaneous and simultaneous measurements of two-dimensional temperature and OH-LIF profiles by combining Rayleigh scattering with laser induced fluorescence (LIF) were demonstrated in a nitrogen-diluted hydrogen (H₂ 30% + N₂ 70%) laminar normal diffusion flame interacting with a large-scale vortex by oscillating central fuel flow or in an inverse diffusion flame by oscillating central airflow. The dynamic behavior of the diffusion flame extinction and reignition during the flame–vortex interaction processes was investigated. The results obtained are described as follows. (1) The width of the reaction zone decreases remarkably, and a decrease in flame temperature and OH-LIF is seen with increasing central airflow in an inverse diffusion flame. OH-LIF increases, and temperature does not change with increasing central fuel flow in a normal diffusion flame. The computations predict the experimental results well, and it is revealed that flame temperature characteristics result from the preferential diffusion of heat and species, which induces excess enthalpy or on enthalpy deficit, and an increase or decrease in H₂ mole fraction in the flame. (2) When a large velocity fluctuation is given to the central flow, the temperature and the OH-LIF at the reaction zone become thin at the convex and circumferential part of the vortex where a high temperature layer exists, and the temperature at the reaction zone is lowered in the inverse flame and the normal flame. (3) The width and temperature of the reaction zone interacting with the vortex recover quickly to that of the laminar steady flame after the vortex passing in the normal flame, but the recovery to that of the steady flame after the vortex passing is delayed in the inverse flame. (4) When a remarkably large velocity fluctuation is given to the central airflow in the inverse flame, thinning of temperature and reaction zone starts at the convex and circumferential part of the vortex, resulting in a and flame extinction completely occurs at the tail part of the vortex and makes the pair of edge flames. The outside edge flame reignites and connects with the upstream reaction zone. The inside edge flame finally extinguishes as the supply of fuel is interrupted by the outside edge flame.     

Compressible turbulent flame speeds of highly turbulent standing flames

In this paper we present the first measurement of turbulent burning velocities of a highly turbulent compressible standing flame induced by shock-driven turbulence in a Turbulent Shock Tube. High-speed schlieren, chemiluminescence, PIV, and dynamic pressure measurements are made to quantify flame–turbulence interaction for high levels of turbulence at elevated temperatures and pressure. Distributions of turbulent velocities, vorticity and turbulent strain are provided for regions ahead and behind the standing flame. The turbulent flame speed is directly measured for the high-Mach standing turbulent flame. From measurements of the flame turbulent speed and turbulent Mach number, transition into a non-linear compressibility regime at turbulent Mach numbers above 0.4 is confirmed, and a possible mechanism for flame generated turbulence and deflagration-to-detonation transition is established.     

Binary collision of a burning droplet and a non-burning droplet of xylene: Outcome regimes and flame dynamics

The binary collisions of a burning droplet and a non-burning droplet of xylene are experimentally investigated. The experimental parameters span an extensive range of Weber number and impact parameter, covering the collision outcome regimes of coalescence, reflexive separation, and stretching separation. A high-speed camera captures the temporal details of the collision process, involving flame spread, visible radiation, and flame distributions around droplets. For reflexive separation and stretching separation, the flame from the droplet spreads to the ligament, surrounding it during the interaction process, and then spreads around separated droplets and satellite droplets. Highly-interactive flames are formed in-between the droplets, with very sooty flames generated for most collisions. For the coalescence case, a swirling flame forms around the rotating coalesced droplet. For similar Weber numbers, visible flame radiation is compared for different collision regimes. The visible flame radiation changes more significantly for the reflexive and stretching separation cases than it does for the coalescence case. The change of the averaged visible flame radiation for reflexive separation and stretching separation is more than two times higher than that for coalescence. The map of three different collision regimes is plotted in the Weber number versus impact parameter domain and compared with available theoretical model predictions. Although the different outcomes of collision with the presence of flame can be well predicted by the model, using fluid properties determined by the averaged properties of the two droplets, the dynamics of the detailed processes involved in the collisions are very interesting and have strong implications on overall combustion behavior that go well beyond the mapped regimes.     

Modelling of lifted turbulent diffusion flames in a channel mixing layer by the flame hole dynamics

The partial quenching structure of turbulent diffusion flames in a turbulent mixing layer is investigated by the method of flame hole dynamics as an effort to develop a prediction model for the turbulent flame lift off. The essence of the flame hole dynamics is derivation of the random walk mapping, from the flame-edge theory, which governs expansion or contraction of the quenching holes initially created by the local quenching events. The numerical simulation for the flame hole dynamics is carried out in two stages. First, a direct numerical simulation is performed for a constant-density fuel–air channel mixing layer to obtain the background turbulent flow and mixing fields, from which a time series of two-dimensional scalar-dissipation-rate array is extracted. Subsequently, a Lagrangian simulation of the flame hole random walk mapping, projected to the scalar dissipation rate array, yields a temporally evolving turbulent extinction process and its statistics on partial quenching characteristics. In particular, the probability of encountering the reacting state, while conditioned with the instantaneous scalar dissipation rate, is examined to reveal that the conditional probability has a sharp transition across the crossover scalar dissipation rate, at which the flame edge changes its direction of propagation. This statistical characteristic implies that the flame edge propagation instead of the local quenching event is the main mechanism controlling the partial quenching events in turbulent flames. In addition, the conditional probability can be approximated by a heavyside function across the crossover scalar dissipation rate.     

Ignition of a sequential combustor: Evidence of flame propagation in the autoignitable mixture

Ignition of the second stage in a lab-scale sequential combustor is investigated experimentally. A fuel mixing section between jet-in-cross-flow injection and the second stage chamber allows the fuel and vitiated, hot cross-flow to partially mix upstream of the main heat release zone. The focus of the present work is on the transient ignition process leading to a stable flame in the second stage. High-speed OH-PLIF as well as OH chemiluminescence imaging is applied to obtain complementary planar and line-of-sight integrated information on the ignition. We find experimental evidence for the co-existence of two regimes dominating the chamber ignition, i.e. autoignition and flame propagation. As the mass flow of the dilution air injected downstream of the first stage is increased (i.e. mixing temperatures in the fuel mixing section are decreased), we transition from an autoignition to a flame propagation dominated regime. Hysteresis in the ignition behavior is observed indicating that the first stage in a sequential combustor may be operated at leaner conditions than required for ignition of the second stage. The time traces of integral heat release obtained simultaneously with a photomultiplier tube show distinct features depending on the dominating regime, which is important for high-pressure testing with limited optical access.     

Further studies of the reaction kernel structure and stabilization of jet diffusion flames

The structure of axisymmetric laminar jet diffusion flames of ethane, ethylene, acetylene, and propane in quasi-quiescent air has been studied numerically in normal earth gravity (1g) and zero gravity (0g). The time-dependent full Navier–Stokes equations with buoyancy were solved using an implicit, third-order accurate numerical scheme, including a C₃-chemistry model and an optically thin-media radiation model for heat losses. Observations of the flames were also made at the NASA Glenn 2.2-Second Drop Tower. For all cases of the fuels and gravity levels investigated, a peak reactivity spot, i.e., reaction kernel, was formed in the flame base, thereby holding a trailing diffusion flame. The location of the reaction kernel with respect to the burner rim depended inversely on the reaction-kernel reactivity or velocity. In the C₂ and C₃ hydrocarbon flames, the H₂–O₂ chain reactions were important at the reaction kernel, yet the CH₃ + O → CH₂O + H reaction, a dominant contributor to the heat-release rate in methane flames studied previously, did not outweigh other exothermic reactions. Instead of the C₁-route oxidation pathway in methane flames, the C₂ and C₃ hydrocarbon fuels dehydrogenated on the fuel side and acetylene was a major hydrocarbon fragment burning at the reaction kernel. The reaction-kernel correlations between the reactivity (the heat-release or oxygen-consumption rate) and the velocity, obtained previously for methane, were developed further for various fuels in more universal forms using variables related to local Damköhler numbers and Peclet numbers.     

A theoretical study of premixed turbulent flame development

Flame development in a statistically stationary and uniform, planar, one-dimensional turbulent flow is theoretically studied. A generalized balance equation for the mean combustion progress variable, which includes turbulent diffusion and pressure-driven transport terms, as well as the mean rate of product creation, is introduced and analyzed by invoking the sole assumption of a self-similar flame structure, well-supported by numerous experiments. The assumption offers the opportunity to simplify the problem by splitting the aforementioned partial differential equation into two ordinary differential equations, which separately model spatial variations of the progress variable and time variations of flame speed and thickness. The self-similar profile of the progress variable, obtained in numerous experiments, is theoretically predicted. Closures of the normalized pressure-driven transport term and mean rate of product creation are obtained. The closed balance equation shows that turbulent diffusion dominates during the initial stage of flame development, followed by the transition to counter-gradient transport in a sufficiently developed flame. A criterion of the transition is derived. The transition is promoted by the heat release and pressure-driven transport. Fully developed mean flame brush thickness and speed are shown to decrease when either density ratio or pressure-driven transport increases. Solutions for the development of the thickness are obtained. The development is accelerated by the pressure-driven transport and heat release.     

Impact of symmetry breaking on the Flame Transfer Function of a laminar premixed flame

This work presents a numerical study of the acoustic response of a laminar flame with tunable asymmetry. A V-shaped premixed flame is stabilised in the wake of a cylindrical flame holder that can be rotated. The configuration is symmetric when the flame holder is fixed but increasing its rotation rate breaks the symmetry of the flow. This configuration is submitted to acoustic forcing to measure the effect of rotation of the flame holder on the Flame Transfer Functions. It appears that the asymmetry of the two flame branches changes their respective time delays, resulting in interference in the global unsteady heat release rate fluctuations. Consequently, the Flame Transfer Function exhibits dips and bumps, which are studied via laminar Direct Numerical Simulation. Potential applications for the control of combustion instabilities are discussed.