Flame sheet dynamics of bluff-body stabilized flames during longitudinal acoustic forcing

Bluff-body stabilized flames are susceptible to combustion instabilities due to interactions between acoustics, vortical disturbances, and the flame. In order to elucidate these flow-flame interactions during an instability, an experimental and computational investigation of the flame-sheet dynamics of a harmonically excited flame was performed. It is shown that the flame dynamics are controlled by three key processes: excitation of shear layer instabilities by the axially oscillating flow, anchoring of the flame at the bluff body, and the kinematic response of the flame to this forcing. The near-field flame features are controlled by flame anchoring and the far-field by kinematic restoration. In the near-field, the flame response grows with downstream distance due to flame anchoring, which prevents significant flame movement near the attachment point. Theory predicts that this results in linear flame response characteristics as a function of perturbation amplitude, and a monotonic growth in magnitude of the flame-sheet fluctuations near the stabilization point, consistent with the experimental data. Farther downstream, the flame response reaches a maximum and then decays due to the dissipation of the vortical disturbances and action of flame propagation normal to itself, which acts to smooth out the wrinkles generated by the harmonic flow forcing. This behavior is strongly non-linear, resulting in significant variation in far-field flame-sheet response with perturbation amplitude.     

Flame balls dynamics in divergent channel

A three-dimensional reaction-diffusion model for lean low-Lewis-number premixed flames with radiative heat losses propagating in divergent channel is studied numerically. Effects of inlet gas velocity and heat-loss intensity on flame structure at low Lewis numbers are investigated. It is found that continuous flame front exists at small heat losses and the separate flame balls settled within restricted domain inside the divergent channel at large heat losses. It is shown that the time averaged flame balls coordinate may be considered as important characteristic analogous to coordinate of continuous flame stabilized in divergent channel.     

The role of radiative losses in self-extinguishing and self-wrinkling flames

We have developed a general theory of non-adiabatic premixed flames that is valid for flames of arbitrary shape that fully accounts for the hydrodynamic and diffusive-thermal processes, and incorporates the effects of volumetric heat losses. The model is used to describe aspects of experimentally observed phenomena of self-extinguishing (SEFs) and self-wrinkling flames (SWFs), in which radiative heat losses play an important role. SEFs are spherical flames that propagate considerable distances in sub-limit conditions before suddenly extinguishing. Our results capture many aspects of this phenomenon including an explicit determination of flame size and propagation speed at quenching. SWFs are hydrodynamically unstable flames in which cells spontaneously appear on the flame surface once the flame reaches a critical size. Our results yield expressions of the critical flame size at the onset of wrinkling and expected cell size beyond the stability threshold. The various possible burning regimes are mapped out in parameter space.     

The rate of expansion of spherical flames

In this paper we investigate the acceleration of the expansion of premixed spherical flames and evolution of the cellular patterns on their surfaces. An asymptotic model is used for the simulations and a spectral numerical algorithm is employed to study flames over large time intervals. Numerous numerical experiments indicate that for large enough time the acceleration of a two-dimensional expanding flame slows down but the expansion rate is still able to reach values significantly exceeding the burning rate of an exactly circular flame. The importance of the effect of forcing was also confirmed and the validity of simulations of sectors of circular flame fronts was studied in order to justify prospective use of the Fourier spectral model for three-dimensional spherical flames.     

AP/HTPB laminate propellant flame structure: Fuel-lean intrinsic instability

Two-dimensional laminate propellant flames of ammonium perchlorate (AP) and hydroxyl-terminated polybutadiene (HTPB) have been observed using infrared (IR) and ultraviolet (UV) emission and transmission imaging. Under fuel-lean conditions and at slightly elevated pressures (4 atm), intrinsic instability has been observed in the form of a leading-edge flame kernel whose location oscillates laterally about the central fuel binder layer. A mechanistic explanation for this behavior is described in terms of local gas-phase equivalence ratio, surface geometry, and gas–solid thermal coupling. The flame structure under these conditions is unique in having a leading-edge flame kernel that appears to be more spatially distinct from the trailing diffusion flame than under nonoscillatory conditions. Other results are reported, including gas-phase rotational and vibrational temperature estimates based on HCl emission imaging spectroscopy. These results add to a growing set of flame and burning surface observations being assembled for the purpose of comprehensive validation of multi-dimensional AP composite propellant computational combustion models.     

Diffusive-thermal instabilities of diffusion flames: onset of cells and oscillations

A comprehensive stability analysis of planar diffusion flames is presented within the context of a constant-density model. The analysis provides a complete characterization of the possible patterns that are likely to be observed as a result of differential and preferential diffusion when a planar flame becomes unstable. A whole range of physical parameters is considered, including the Lewis numbers associated with the fuel and the oxidizer, the initial mixture fraction, and the flow conditions. The two main forms of instability are cellular flames, obtained primarily in fuel-lean systems when the Lewis numbers are generally less than one, and planar pulsations, obtained in fuel-rich systems when the Lewis numbers are generally larger than one. The cellular instability is predominantly characterized by stationary cells of characteristic dimension comparable to the diffusion length, but smaller cells that scale on the reaction zone thickness are also possible near extinction conditions. The pulsating instability is characterized by planar oscillations normal to the flame sheet with a well-defined frequency comparable to the reciprocal of the diffusion time; high-frequency modes are also possible just prior to extinction. The analysis also alludes to other possible patterns, such as oscillating cellular structures, which result from competing modes of instability of comparable and/or disparate scales. The expected pattern depends of course on the underlying physical parameters. Consequently, stability boundaries have been identified for the onset of one or another form of the instability. The conditions for the onset of cellular and pulsating flames, as well as the predicted cell size and the frequency of oscillations, compare well with the experimental record.     

Influence of advective heat flux on extinction scalar dissipation rate and velocity in negative edge flames

A counterflow flame geometry, which has previously been experimentally shown to produce stable negative edge flames, was studied using numerical simulations. In this geometry, the flame edge is formed off the counterflow centreline owing to a local increase in scalar dissipation rate. Hot products from the stable nonpremixed flame on the centreline flow through the edge at velocities of ～ 1–5 m/s. The size of the counterflow burner and the gas flowrates are varied in the simulations to alter the flame strength and velocity at the flame edge. The advection of products through the edge is shown to extend the flame extinction to higher scalar dissipation rates than required for centreline extinction. For high velocities, the scalar dissipation rate required for flame extinction can be related to the centreline extinction value by considering only the effect of energy addition to the flame edge via advection. However, for lower edge flame velocities, the effects of increased thermal and species diffusion through the edge must also be included. Since the advection at the edge is a product of both the local velocity and temperature gradient, a single correlation between the scalar dissipation rate and the negative edge flame velocity does not exist.     

Influence of edge propagation on downward flame spread over three-dimensional PMMA samples

This systematic experimental study measures the effect of flame propagation along vertical edges on the overall downward spread of flames using Polymethyl Methacrylate (PMMA). Samples with a wide range of regular cross-sections – from triangular through octagonal – as well as irregular ones, have been used to test a large variation of internal angles. A MATLAB-based tool was used to calculate instantaneous spread rate for central and edge flames. The edge flame is shown to significantly enhance the spread rate, as much as five times, in respect to samples with no edges. This amplification is shown to depend primarily on the internal angle at the edge (the smaller the angle, the faster the flame) and fuel thickness, and not on other factors such as aspect ratio or cross-sectional area. Using a phenomenological argument, the edge propagation rate is correlated to the spread rate over an equivalent cylindrical fuel (the limiting shape with infinite edges) with an effective radius obtained from the geometry of the edges and the diffusion length scale of the solid phase. A formula for flame spread over cylindrical fuel from the literature is used to link the new results to existing models. Both thick and thin limits are shown to encompass the wide range of three-dimensional spread rate data within the effective radius (the independent variable), which can be determined from the known parameters. Based on these results, different types of cross-sectional areas can be sorted in the order of their inherent fire safety characteristics.     

Experimental research on combustion dynamics of medium-scale fire whirl

The medium-scale fire whirl was extensively investigated by experimental means, in order to establish correlations of the burning rate, flame height and flame temperature of fire whirl, and to clarify the difference between fire whirls and general pool fires. Experimental observations and data confirmed that a free burning fire whirl is a highly stable burning phenomenon with large quasi-steady periods. Burning rates of fire whirls depend on pool diameter similarly to those of general pool fires; however the transition turbulent burning occurs sooner as the pool diameter increases. The lip height seems to have little effect on the burning rate of fire whirls. The correlation was proposed to couple the height of fire whirl to the fire release rate and ambient circulation. It correlates the data from both this work and the literature. Radial temperature profiles in the continuous region of the fire whirl were confirmed to be hump-type, implying the existence of fuel-rich inner core. The pool diameter and heat release rate do not significantly affect the radial temperature profiles in non-dimensional radial coordinates. It was found that the fire plume of fire whirl involves three distinct zones just like that of pool fire, but with different normalized ranges. Fire whirls maintain a higher ratio of continuous flame height to the overall flame height, and also higher maximum centerline excess temperature in continuous flame region, as compared to general pool fires. It was further demonstrated that the fire whirl plume at its origin behaves like a turbulent jet with moderate swirling, and then tends to become buoyancy dominated downstream, with slight swirling. With an increase in dimensionless height adjusted by the plume origin, the plume centerline excess temperature decays rapidly and approaches the theoretical value of −5/3 for free buoyancy plume.     

Experimental measurements of two-dimensional planar propagating edge flames

The study of edge flames has received increased attention in recent years. This work reports the results of a recent study into two-dimensional, planar, propagating edge flames that are remote from solid surfaces (called here, “free-layer” flames, as opposed to layered flames along floors or ceilings). They represent an ideal case of a flame propagating down a flammable plume, or through a flammable layer in microgravity. The results were generated using a new apparatus in which a thin stream of gaseous fuel is injected into a low-speed laminar wind tunnel thereby forming a flammable layer along the centerline. An airfoil-shaped fuel dispenser downstream of the duct inlet issues ethane from a slot in the trailing edge. The air and ethane mix due to mass diffusion while flowing up towards the duct exit, forming a flammable layer with a steep lateral fuel concentration gradient and smaller axial fuel concentration gradient. We characterized the flow and fuel concentration fields in the duct using hot wire anemometer scans, flow visualization using smoke traces, and non-reacting, numerical modeling using COSMOSFloWorks. In the experiment, a hot wire near the exit ignites the ethane-air layer, with the flame propagating downwards towards the fuel source. Reported here are tests with the air inlet velocity of 25 cm/s and ethane flows of 967-1299 sccm, which gave conditions ranging from lean to rich along the centerline. In these conditions the flame spreads at a constant rate faster than the laminar burning rate for a premixed ethane-air mixture. The flame spread rate increases with increasing transverse fuel gradient (obtained by increasing the fuel flow rate), but appears to reach a maximum. The flow field shows little effect due to the flame approach near the igniter, but shows significant effect, including flow reversal, well ahead of the flame as it approaches the airfoil fuel source.     

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.     

Hydrodynamic and chemical scaling for blow-off dynamics of lean premixed flames stabilized on a meso-scale bluff-body

Direct numerical simulations were conducted to investigate the effect of two parameters, density ratio and laminar flame speed, on the conditions of the onset of local extinction and blow-off of lean premixed flames, stabilized on a meso-scale bluff-body in hydrogen-air and syngas–air mixtures. A total of six simulation cases were considered as isolated comparison of the two parametric effects of the fluid dynamic instability and flame time scale. For all cases under study, the general flame development towards the blow-off limit showed a sequence of five distinct modes, with possible cyclic patterns among the different modes for a range of velocity conditions. The onset of local extinction was observed during the asymmetric vortex shedding and vortex street mode. As the density ratio is decreased, the flow inunder reviewstability is promoted through the increased sinuous mode, and such behavior was properly scaled by the Strouhal number. Although the blow-off velocity is altered by the fluid dynamic effects, the condition for the onset of local extinction and blow-off was mainly dictated by the competition between flow residence time associated with the lateral flame motion and ignition delay of the local mixtures. Time scale analysis supported the validity of the findings across all the cases investigated.     

Cellular instabilities of expanding hydrogen/propane spherical flames at elevated pressures: theory and experiment

An experimental and theoretical investigation of the onset of cellular instabilities on spherically expanding flames in mixtures of hydrogen and propane in air at elevated pressures was conducted. Critical conditions for the onset of instability were measured and mapped out over a range of pressures and mixture compositions. An asymptotic theory of hydrodynamic and diffusional-thermal cell development on flames in mixtures comprised of two scarce fuels burning in air was also formulated. Predicted values of Peclet number, defined as the flame radius at the onset of instability normalized by the flame thickness, were shown to compare favorably with the experimentally measured values.     

Three-dimensional numerical simulations of cellular jet diffusion flames

Recent experimental investigations have demonstrated that the appearance of particular cellular states in circular non-premixed jet flames significantly depends on a number of parameters, including the initial mixture strength, reactant Lewis numbers, and proximity to the extinction limit (Damköhler number). For CO₂-diluted H₂/O₂ jet diffusion flames, these studies have shown that a variety of different cellular patterns or states can form. For given fuel and oxidizer compositions, several preferred states were found to co-exist, and the particular state realized was determined by the initial conditions. To elucidate the dynamics of cellular instabilities, circular non-premixed jet flames are modeled with a combination of three-dimensional numerical simulation and linear stability analysis (LSA). In both formulations, chemistry is described by a single-step, finite-rate reaction, and different reactant Lewis numbers and molecular weights are specified. The three-dimensional numerical simulations show that different cellular flames can be obtained close to extinction and that different states co-exist for the same parameter values. Similar to the experiments, the behavior of the cell structures is sensitive to (numerical) noise. During the transient blow-off process, the flame undergoes transitions to structures with different number of cells, while the flame edge close to the nozzle oscillates in the streamwise direction. For conditions similar to the experiments discussed, the LSA results reveal various cellular instabilities, typically with azimuthal wavenumber m = 1–6. Consistent with previous theoretical work, the propensity for the cellular instabilities is shown to increase with decreasing reactant Lewis number and Damköhler number.     

Flame stabilization with a tubular flame

A new technique to stabilize a flame in a high-velocity stream with use of a tubular flame has been proposed. To elucidate the validity of this technique, an experiment has been conducted by mounting a tubular flame burner on the nozzle. Flame stability limits and temperature distributions around the burner port have been determined, and experiments have been extended to the ducted combustion to measure pressure fluctuations and to analyze the burned gases. Results show that the tubular flame can successfully stabilize the main flame up to 130 m/s, which is the upper limit of the present supply facility. The main flame is well anchored at the exit of the nozzle, and the tubular flame efficiently supplies heat and radicals to the main flame. In the ducted combustion, the pressure fluctuations are reduced significantly. The exhaust gas analyses, however, indicate that an almost chemical equilibrium condition can be achieved at 50 m/s, but not at 90 and 130 m/s. Since the energy input relative to the main flame is just 6.1% at 130 m/s, the present tubular burner is not enough to burn all the unburned gas completely at high velocities, although the main flame can be anchored. The slit length and/or the slit width of the tubular flame burner should be larger to overcome this shortage. From the above results, it is concluded that the tubular flame has a potential for stabilizing a flame in a high-speed stream.     

Structure and dynamic response of radiative diffusion flames

Nitrogen-diluted hydrogen burning in air is modeled numerically using a constant density and one-step reaction model in a plane two-dimensional counterflow configuration. An optically thin assumption is used to investigate the effects of radiation on the dynamics, structure, and extinction of diffusion flames. While there exist dual steady-state extinction limits for the 1D radiative flame response, it is found that as the 1D radiative extinction point is approached the 1D low-stretch diffusion flame exhibits oscillatory response, even with sub-unity Lewis number fuel. These radiation-induced limit cycle oscillations are found to have increasing amplitude and decreasing frequency as the stretch rate is reduced. Flame oscillation eventually leads to permanent extinction at the stretch rate which is larger than the steady-state radiative extinction value. Along the 1D radiative response curve, the transition from 1D flame to 2D structure and the differences in the resulting 2D flame patterns are also examined using a variety of initial profiles, with special emphasis on the comparison of using the initial profiles with and without a flame edge. Similar to the previous studies on the high-stretch adiabatic edge flames using the same configuration, the high-stretch radiative flames are found to resist 1D blow-off quenching through various 2D structures, including propagating front and steady cellular flames for initial profiles with and without flame edges. For all initial profiles studied, the low-stretch radiative flames are also found to exhibit different 2D flame phenomena near the 1D radiative extinction limit, such as transient cellular structures, steady cellular structures, and pulsating ignition fronts. Although the results demonstrate the presence of low-stretch and high-stretch 2D bifurcation branches close to the corresponding 1D extinction limits irrespective of the initial profile used, particular 2D flame structures in certain stretch rate range are initial profile dependent. The existence of two-dimensional flame structures beyond the 1D steady-state radiative extinction limit suggests that the flammable range is expanded as compared to that predicted by the 1D model. Hence, multi-dimensional flame patterns need to be accounted for when determining the flammability limits for a given system.     

Combined experimental and computational study of laminar, axisymmetric hydrogen–air diffusion flames

We investigate the structure of two-dimensional, axisymmetric, laminar hydrogen–air flames in which a cylindrical fuel stream is surrounded by coflowing air, using laser-diagnostic and computational methods. Spontaneous Raman scattering and coherent anti-Stokes Raman scattering (CARS) are used to measure the distributions of major species and temperature. Computationally, we solve the governing conservation equations for mass, momentum, energy, and species, using detailed chemistry and transport. The fuel is diluted with nitrogen (1:1) to reduce heat transfer to the burner, to match the zero temperature gradient at the fuel exit. Three average fuel exit velocities are studied: 18, 27, and 50 cm/s. Comparisons of the measured and computed results are performed for radial profiles at a number of axial positions, and along the axial centerline. Peak major species mole fractions and temperatures are quantitatively predicted by the computations, and the axial species profiles are predicted to within the experimental uncertainty. In the radial profiles studied, base-case computations excluding thermal diffusion of light species were in excellent agreement with the measurements. While the addition of thermal diffusion led to some discrepancy with the measured results, the magnitude of the differences was no more than 25%. The computations predicted the axial centerline profiles from the burner exit to the maximum temperature well, though the experimental temperatures in the downstream mixing region decreased somewhat faster than the computed profiles. Radiative losses are seen to be negligible in these flames, and changes in transport properties and variations in initial flow velocities generally led to only modest changes in the axial profiles. The results also show that the detailed axial profiles of major species and temperature at different fuel jet velocities scale quantitatively with the jet velocity.