Transient soot dynamics in turbulent nonpremixed ethylene-air counterflow flames

The dynamics of soot formation in turbulent ethylene-air nonpremixed counterflow flames is studied using direct numerical simulation (DNS) with a semi-empirical soot model and the discrete ordinate method (DOM) as a radiation solver. Transient characteristics of soot behavior are studies by a model problem of flame interaction with turbulence inflow at various intensities. The interaction between soot and turbulence reveals that the soot volume fraction depends on the combined effects of the local conditions of flow, temperature, and fuel concentration, while the soot number density depends predominantly on the high temperature regions. Depending on the relative strength between mixing and reaction, the effects of turbulence on the soot formation lead to three distinct paths in deviating the data points away from the laminar flame conditions. It is found that turbulence has twofold effects of increasing the overall soot yield by generating additional flame volume and of reducing soot by dissipating soot pockets out of high-temperature regions. The relative importance between the two effects depends on the relative length scales of turbulence and flame, suggesting that a nonmonotonic response of soot yield to turbulence level may be expected in turbulent combustion.     

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.     

Synthesis of carbon nanotubes on Ni-alloy and Si-substrates using counterflow methane–air diffusion flames

We have conducted an experimental study to investigate the synthesis of multi-walled carbon nanotubes (CNTs) in counterflow methane–air diffusion flames, with emphasis on effects of catalyst, temperature, and the air-side strain rate of the flow on CNTs growth. The counterflow flame was formed by fuel (CH₄ or CH₄ + N₂) and air streams impinging on each other. Two types of substrates were used to deposit CNTs. Ni-alloy (60% Ni + 26% Cr + 14% Fe) wire substrates synthesized curved and entangled CNTs, which have both straight and bamboo-like structures; Si-substrates with porous anodic aluminum oxide (AAO) nanotemplates synthesized well-aligned, self-assembled CNTs. These CNTs grown inside nanopores had a uniform geometry with controllable length and diameter. The axial temperature profiles of the flow were measured by a 125 μm diameter Pt/10% Rh–Pt thermocouple with a 0.3 mm bead junction. It was found that temperature could affect not only the success of CNTs synthesis, but also the morphology of synthesized CNTs. It was also found, against previous general belief, that there was a common temperature region (1023–1073 K) in chemical vapor deposition (CVD) and counterflow diffusion flames where CNTs could be produced. CNTs synthesized in counterflow flames were significantly affected by air-side strain rate not through the residence time, but through carbon sources available for CNTs growth. Off-symmetric counterflow flames could synthesize high-quality CNTs because with this configuration carbon sources at the fuel side could easily diffuse across the stagnation surface to support CNTs growth. These results show the feasibility of using counterflow flames to synthesize CNTs for particular applications such as fabricating nanoscale electronic devices.     

Structure of partially premixed n-heptane–air counterflow flames

To avoid the complexities associated with the droplet/vapor transport and nonuniform evaporation processes, a fundamental investigation of liquid fuel combustion in idealized configurations is very useful. An experimental–computational investigation of prevaporized n-heptane nonpremixed and partially premixed flames established in a counterflow burner is described. There is a general agreement between various facets of our nonpremixed flame measurements and the literature data. The partially premixed flames are characterized by a double flame structure. This becomes more distinct as the strain rate decreases and partial premixing increases, which also increases the separation distance between the two reaction zones. The peak partially premixed flame temperature increases with increasing premixing of the fuel stream. The peak CO₂ and H₂O concentrations are relatively insensitive to partial premixing. The CO and H₂ peak concentrations on the premixed side increase as the fuel-side equivalence ratio decreases. These species are transported to the nonpremixed reaction zone where they oxidize. The C₂ species have peaks in the premixed reaction zone. The concentrations of olefins are ten times larger than those of the corresponding paraffins. The oxidizer is present in partially premixed flames throughout the combustion system and there are no regions characterized by simultaneous high temperature and high fuel concentration. As a result, pyrolysis reactions leading to soot formation are greatly diminished.     

Mechanism on oscillating lifted flames in nonpremixed laminar coflow jets

The oscillating lifted flame in a laminar nonpremixed nitrogen-diluted fuel jet is known to be a result of buoyancy, though the detailed physical mechanism of the initiation has not yet been properly addressed. We designed a systematic experiment to test the hypothesis that the oscillation is driven by competition between the positive buoyancy of flame and the negative buoyancy of a fuel stream heavier than the ambient air. The positive buoyancy was examined with various flame temperatures by changing fuel mole fraction, and the negative buoyancy was investigated with various fuel densities. The density of the coflow was also varied within a certain range by adding either helium or carbon dioxide to air, to study how it affected the positive and negative buoyancies at the same time. As a result, we found that the range of oscillation was well-correlated with the positive and the negative buoyancies; the former stabilized the oscillation while the latter triggered instability and became a source of the oscillation. Further measurements of the flow fields and OH radicals evidenced the important role of the negative buoyancy on the oscillation, detailing a periodic variation in the unburned flow velocity that affected the displacement of the flame.     

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.     

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.     

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.     

Characterization of the effects of pressure and hydrogen concentration on laminar burning velocities of methane–hydrogen–air mixtures

The aim of the present work was to characterize both the effects of pressure and of hydrogen addition on methane/air premixed laminar flames. The experimental setup consists of a spherical combustion chamber coupled to a classical shadowgraphy system. Flame pictures are recorded by a high speed camera. Global equivalence ratios were varied from 0.7 to 1.2 for the initial pressure range from 0.1 to 0.5 MPa. The mole fraction of hydrogen in the methane + hydrogen mixture was varied from 0 to 0.2. Experimental results were compared to calculations using a detailed chemical kinetic scheme (GRIMECH 3.0). First, the results for atmospheric laminar CH₄/air flames were compared to the literature. Very good agreements were obtained both for laminar burning velocities and for burned gas Markstein length. Then, increasing the hydrogen content in the mixture was found to be responsible for an increase in the laminar burning velocity and for a reduction of the flame dependence on stretch. Transport effects, through the reduction of the fuel Lewis number, play a role in reducing the sensitivity of the fundamental flame velocity to the stretch. Finally, when the pressure was increased, the laminar burning velocity decreased for all mixtures. The pressure domain was limited to 0.5 MPa due to the onset of instabilities at pressures above this value.     

Structure of propagating edge diffusion flames in hydrocarbon fuel jets

Direct numerical simulations with a C₃-chemistry model have been performed to investigate the transient behavior and internal structure of flames propagating in an axisymmetric fuel jet of methane, ethane, ethylene, acetylene, or propane in normal earth gravity (1g) and zero gravity (0g). The fuel issued from a 3-mm-i.d. tube into quasi-quiescent air for a fixed mixing time of 0.3 s before it was ignited along the centerline where the fuel–air mixture was at stoichiometry. The edge of the flame formed a vigorously burning peak reactivity spot, i.e., reaction kernel, and propagated through a flammable mixture layer, leaving behind a trailing diffusion flame. The reaction kernel broadened laterally across the flammable mixture layer and possessed characteristics of premixed flames in the direction of propagation and unique flame structure in the transverse direction. The reaction kernel grew wings on both fuel and air sides to form a triple-flame-like structure, particularly for ethylene and acetylene, whereas for alkanes, the fuel-rich wing tended to merge with the main diffusion flame zone, particularly methane. The topology of edge diffusion flames depend on the properties of fuels, particularly the rich flammability limit, and the mechanistic oxidation pathways. The transit velocity of edge diffusion flames, determined from a time series of calculated temperature field, equaled to the measured laminar flame speed of the stoichiometric fuel–air mixtures, available in the literature, independent of the gravity level.     

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.     

Effect of velocity gradient on propagation speed of tribrachial flames in laminar coflow jets

The effect of velocity gradient on the propagation speed of tribrachial flame edge has been investigated experimentally in laminar coflow jets for propane fuel. It was observed that the propagation speed of tribrachial flame showed appreciable deviations at various jet velocities in high mixture fraction gradient regime. From the similarity solutions, it was demonstrated that the velocity gradient varied significantly during the flame propagation. To examine the effect of velocity gradient, detail structures of tribrachial flames were investigated from OH LIF images and Abel transformed images of flame luminosity. It was revealed that the tribrachial point was located on the slanted surface of the premixed wing, and this slanted angle was correlated with the velocity gradient along the stoichiometric contour. The temperature field was visualized qualitatively by the Rayleigh scattering image. The propagation speed of tribrachial flame was corrected by considering the direction of flame propagation with the slanted angle and effective heat conduction to upstream. The corrected propagation speed of tribrachial flame was correlated well. Thus, the mixture fraction gradient together with the velocity gradient affected the propagation speed.     

Structure and propagation of two-dimensional,partially premixed,laminar flames in diesel engine conditions

We investigate the influence of inflow velocity (V_in) and scalar dissipation rate (χ) on the flame structure and stabilisation mechanism of steady, laminar partially premixed n-dodecane edge flames stabilised on a convective mixing layer. Numerical simulations were performed for three different χ profiles and several V_in (V_in = 0.2 to 2.5m/s). The ambient thermochemical conditions were the same as the Engine Combustion Network’s (ECN) Spray A flame, which in turn represents conditions in a typical heavy duty diesel engine. The results of a combustion mode analysis of the simulations indicate that the flame structure and stabilisation mechanism depend on V_in and χ. For low V_in the flame is attached. Increasing V_in causes the high-temperature chemistry (HTC) flame to lift-off, while the low-temperature chemistry (LTC) flame is still attached. A unique speed S_R associated with this transition is defined as the velocity at which the lifted height has the maximum sensitivity to changes in V_in. This transition velocity is negatively correlated with χ. Near $V_{in}=S_R$ a tetrabrachial flame structure is observed consisting of a triple flame, stabilised by flame propagation into the products of an upstream LTC branch. Further increasing the inlet velocity changes the flame structure to a pentabrachial one, where an additional HTC ignition branch is observed upstream of the triple flame and ignition begins to contribute to the flame stabilisation. At large V_in, the LTC is eventually lifted, and the speed at which this transition occurs is insensitive to χ. Further increasing V_in increases the contribution of ignition to flame stabilisation until the flame is completely ignition stabilised. Flow divergence caused by the LTC branch reduces the χ at the HTC branches making the HTC more resilient to χ. The results are discussed in the context of identification of possible stabilisation modes in turbulent flames.     

Experimental study of air dilution in oxy-liquid fuel flames

The influence of oxidizer dilution in oxy-liquid ethanol flames is experimentally investigated by using a coaxial air-assisted injector positioned in a vertical combustion chamber. This study accounts for the influence of a two-phase mode since two different injector geometries are used: for the first configuration, a vaporization mode is observed at nominal power in oxy conditions, while for the second one, a brush mode is observed. Dilution with air is applied by keeping oxidizer velocity constant. Flame structure is observed through CH emission: dilution leads to an increase in the flame diameter, and collective effects of two-phase combustion are encouraged. The effect of dilution on oxy flame stability is also studied: for a given oxygen mass fraction in the oxidizer, the oxidizer flow rate is increased until extinction occurs. Dilution leads to a less stable flame, which may be essentially explained by the decrease in laminar flame speed with dilution. For high oxidizer dilution levels, the change in flame structure might be another parameter to consider. Finally, species concentrations are measured using a standard gas sampling technique. NO and CO evolutions with dilution are different between both two-phase combustion regimes. An empirical approach based on thermal NO mechanism and CO oxidation reaction enables one to explain the evolutions for brush mode. For vaporization mode, the residence time in burned gases is also to be considered.     

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.     

Combustion dynamics of inverted conical flames

An inverted conical flame anchored on a central bluff-body in an unconfined burner configuration features a distinctive acoustic response. This configuration typifies more complex situations in which the thermo-acoustic instability is driven by the interaction of a flame with a convective vorticity mode. The axisymmetric geometry investigated in this article features a shear region between the reactive jet and the surrounding atmosphere. It exhibits self-sustained oscillations for certain operating conditions involving a powerful flame collapse phenomenon with sudden annihilation of flame surface area. This is caused by a strong interaction between the flame and vortices created in the outer jet shear layer, a process which determines the amplitude of heat release fluctuation and its time delay with respect to incident velocity perturbations. This process also generates an acoustic field that excites the burner and synchronizes the vortex shedding mechanism. The transfer functions between the velocity signal at the burner outlet and heat release are obtained experimentally for a set of flow velocities fluctuations levels. It is found that heat release fluctuations are a strong function of the incoming velocity perturbation amplitude and that the time delay between these two quantities is mainly determined by the convection of the large scale vortices formed in the jet shear layer. A model is formulated, which suitably describes the observed instabilities.     

A numerical study on extinction behaviour of laminar micro-diffusion flames

We conducted a numerical study on the fluid dynamic, thermal and chemical structures of laminar methane–air micro flames established under quiescent atmospheric conditions. The micro flame is defined as a flame on the order of one millimetre or less established at the exit of a vertically-aligned straight tube. The numerical model consists of convective–diffusive heat and mass transport with a one-step, irreversible, exothermic reaction with selected kinetics constants validated for near-extinction analyses. Calculations conducted under the burner rim temperature 300 K and the adiabatic burner wall showed that there is the minimum burner diameter for the micro flame to exist. The Damköhler number (the ratio of the diffusive transport time to the chemical time) was used to explain why a flame with a height of less than a few hundred microns is not able to exist under the adiabatic burner wall condition. We also conducted scaling analysis to explain the difference in extinction characteristics caused by different burner wall conditions. This study also discussed the difference in governing mechanisms between micro flames and microgravity flames, both of which exhibit similar spherical flame shape.     

Carbon nanotube synthesis on catalytic metal alloys in methane/air counterflow diffusion flames

Various morphologies of multi-walled carbon nanotubes (MWNTs) are grown catalytically on metal-alloy probes in counterflow diffusion flames using methane as fuel. Carbon nanotube (CNT) properties and morphologies are investigated as functions of local gas-phase temperatures, C-related species concentrations (e.g. C₂H₂, CO), sampling positions, C₂H₂ adding to the fuel, and metal-alloy compositions (i.e., Fe, Fe/Cr, Ni/Cu, Ni/Ti, Ni/Cr, Ni/Cr/Fe). MWNTs grow optimally in non-sooty regions of the flames. C₂H₂ addition is found to promote direct synthesis of vertically well-aligned MWNTs with uniform diameters from Ni/Cr/Fe and Ni/Ti alloys.