Particle formation in opposed-flow diffusion flames of ethylene: An experimental and numerical study

An experimental and numerical study on particles inception and growth is performed in opposed-flow diffusion flames of ethylene and air characterized by different sooting tendencies. Spectrally resolved UV-visible laser induced fluorescence, laser induced incandescence and laser light scattering measurements are used to characterize different classes of combustion-generated compounds. A detailed kinetic model accounting for both gas-phase and particle formation is used. Comparison between experimental results and numerical predictions gives a qualitative view of the mechanism of particle formation in opposed-flow flames.Particle inception is the result of both chemical growth and coagulation of aromatic compounds. In the region close to the flame front where the temperature is relatively high and radicals are abundant, the particle inception is due to a chemical growth mechanism by which aromatic molecules add aromatic radicals leading to the formation of biphenyl-like structures. The growth process continues as high-molecular mass aromatics are moved away from the flame zone towards the stagnation plane by the addition of acetylene and other aromatics forming particles of increasing sizes. Graphitization of these particles and thermal annealing lead to the formation of soot particles. At relatively lower temperatures, found across the stagnation plane, particles growth still occurs and it is mainly due to a process of physical coagulation of PAHs.The experimental and numerical results obtained in this work demonstrate and explain the sensibility of inception and growth of particles to radical concentration and temperature in opposed-flow flame configurations.     

Pressure and temperature dependence of soot in highly controlled counterflow ethylene diffusion flames

Soot volume fraction and dispersion index were measured by pyrometry in a series of highly controlled counterflow diffusion flames, with peak temperatures, T_max, spanning a few hundred degrees and pressure covering the 0.1–0.8 MPa range. An unprecedented level of control was implemented by selecting flames with a self-similar structure to ensure that the normalized temperature-time history experienced by the reactants was the same, regardless of pressure. The self-similarity was verified by suitably rescaling the transverse coordinate with respect to a characteristic diffusion length. At constant T_max, the soot volume fraction increases approximately by two orders of magnitude as the pressure is raised from 1 atm to 4 atm, and by one to two additional orders of magnitude with an additional doubling of the pressure to 8 atm. At constant pressure, the soot load spans two to three orders of magnitude and soot formation exhibits increased sensitivity to temperature as the pressure is raised. Soot inception occurs near the flame, with an increase in soot concentration that becomes steeper at higher T_max. The increase is accompanied by a decrease in the dispersion exponent that is suggestive of dehydrogenation and aging of the particles and is sharper at higher T_max. Soot experiences continuous growth in a monotonically decreasing temperature field until it is convected away radially at the stagnation plane, with essentially no opportunity for oxidation. Evidence of two distinct mechanisms for soot formation was found: the classic high temperature, high activation energy process affecting soot formed in the vicinity of the flame and followed by dehydrogenation; and a relatively low-temperature, zero activation energy process, associated with the increase in volume fraction at low-temperatures in proximity of the stagnation plane. The latter is tentatively attributed to dimerization of aromatics, as revealed by the concurrent increase in the dispersion index corresponding to an increase in the particle hydrogen content.     

The effects of dimethyl ether and ethanol on benzene and soot formation in ethylene nonpremixed flames

Soot volume fractions, C1-C12 hydrocarbon concentrations, and gas temperature were measured in ethylene/air nonpremixed flames with up to 10% dimethyl ether (CH₃OCH₃) or ethanol (CH₃CH₂OH) added to the fuel. The measurement techniques were laser-induced incandescence, photoionization mass spectroscopy, and thermocouples. Oxygenated hydrocarbons have been proposed as soot-reducing fuel additives, and nonpremixed flames are good laboratory-scale models of the fuel-rich reaction zones where soot forms in many full-scale combustion devices. However, addition of both dimethyl ether and ethanol increased the maximum soot volume fractions in the ethylene flames studied here, even though ethylene is a much sootier fuel than either oxygenate. Furthermore, dimethyl ether produced a larger increase in soot even though neat dimethyl ether flames produce less soot than neat ethanol flames. The detailed species measurements suggest that the oxygenates increase soot concentrations because they decompose to methyl radical, which promotes the formation of propargyl radical (C₃H₃) through C1 + C2 addition reactions and consequently the formation of benzene through propargyl self-reaction. Dimethyl ether has a stronger effect than ethanol because it decomposes more completely to methyl radical. Ethylene does not decompose to methyl, so its flames are particularly sensitive to this mechanism; the alkane-based fuels used in most practical fuels do decompose to methyl radical, so the mechanism will be much less important for practical devices.     

Experimental comparison of soot formation in turbulent flames of Diesel and surrogate Diesel fuels

Soot formation is compared in turbulent diffusion flames burning a commercial Diesel and two Diesel surrogates containing n-decane and α-methylnaphthalene. A burner equipped with a high-efficiency atomisation system has been specially designed and allows the stabilisation of liquid fuels flames with similar hydrodynamics conditions. The initial surrogate composition (70% n-decane, 30% α-methylnaphthalene) was previously used in the literature to simulate combustion in Diesel engines. In this work, a direct comparison of Diesel and surrogates soot tendencies is undertaken and relies on soot and fluorescent species mappings obtained respectively by Laser-Induced Incandescence (LII) at 1064 nm and Laser-Induced Fluorescence at 532 nm. LIF was assigned to soot precursors and mainly to high-number ring Polycyclic Aromatic Hydrocarbons (PAH). The initial surrogate was found to form 40% more soot than the tested Diesel. Consequently, a second surrogate containing a lower α-methylnaphthalene concentration (20%) has been formulated. That composition which presents a Threshold Soot Index (TSI) very close to Diesel one is also consistent with our Diesel composition that indicates a relatively low PAH content. The spatially resolved measurements of soot and fluorescent soot precursors are quite identical (in shape and intensity) in the Diesel and in the second surrogate flames. Furthermore the concordance of the LII temporal decays suggests that a similar growth of the primary soot particles has occurred for Diesel and surrogates. In addition, the comparison of the LII fluence curves indicates that physical/optical properties of soot contained in the different flames might be similar. The chemical composition present at the surface of soot particles collected in Diesel and surrogate flames has been obtained by laser-desorption ionisation time-of-flight mass spectrometry. An important difference is found between Diesel and surrogate samples indicating the influence of the fuel composition on soot content.     

Spatial scales of extinction and dissipation in the near field of non-premixed turbulent jet flames

Simultaneous high-resolution Rayleigh scattering imaging and planar laser-induced fluorescence (PLIF) of OH are combined to measure the dissipative scales associated with thermal mixing and the structure and scales of extinguished regions of the reaction zone. Measurements are performed throughout the near field (x/d = 5, 10, 15, 20) of two turbulent, non-premixed methane/hydrogen/nitrogen jet flames with Re = 15,200 and 22,800 (flames DLR-A and DLR-B of the TNF workshop). Locally extinguished regions are identified by discontinuities in the OH layers, and the extinction hole sizes are measured. For each flame, the probability density function of the hole sizes is very similar throughout the entire near field, with the most likely hole size being 1.9 mm in DLR-A and 1.1 mm in DLR-B. Extinction events are equally probable at all measurement locations in DLR-A. In the DLR-B flame, there is a progression from frequent extinction close to the nozzle to more continuous reaction zones further downstream. The approximate instantaneous location of the stoichiometric contour is determined using the OH-PLIF images, enabling statistical analysis of dissipative scales conditioned on rich and lean conditions. The widths of the thin, elongated structures that dominate the thermal dissipation field are measured. Statistics of this microscale are qualitatively similar in both flames, with the higher Reynolds number producing smaller scales throughout the flow field. For dissipation layers in rich regions, the layer widths increase significantly with increasing temperature, while on the lean side the layer widths decrease with increasing temperature.     

Numerical simulations of soot aggregation in premixed laminar flames

In this paper we make use of a detailed particle model and stochastic numerical methods to simulate the particle size distributions of soot particles formed in laminar premixed flames. The model is able to capture the evolution of mass and surface area along with the full structural detail of the particles. The model is validated against previous models for consistency and then used to simulate flames with bimodal and unimodal soot particle distributions. The change in morphology between the particles from these two types of flames provides further evidence of the interplay among nucleation, coagulation, and surface rates. The results confirm the previously proposed role of the strength of the particle nucleation source in defining the instant of transition from coalescent to fractal growth of soot particles.     

Modelling of turbulent scalar flux in turbulent premixed flames based on DNS databases

A transport equation for scalar flux in turbulent premixed flames was modelled on the basis of DNS databases. Fully developed turbulent premixed flames were obtained for three different density ratios of flames with a single-step irreversible reaction, while the turbulent intensity was comparable to the laminar burning velocity. These DNS databases showed that the countergradient diffusion was dominant in the flame region. Analyses of the Favre-averaged transport equation for turbulent scalar flux proved that the pressure related terms and the velocity–reaction rate correlation term played important roles on the countergradient diffusion, while the mean velocity gradient term, the mean progress variable gradient term and dissipation terms suppressed it. Based on these analyses, modelling of the combustion-related terms was discussed. The mean pressure gradient term and the fluctuating pressure term were modelled by scaling, and these models were in good agreement with DNS databases. The dissipation terms and the velocity–reaction rate correlation term were also modelled, and these models mimicked DNS well.     

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.     

A numerical study of auto-ignition in turbulent lifted flames issuing into a vitiated co-flow

This paper presents a numerical study of auto-ignition in simple jets of a hydrogen–nitrogen mixture issuing into a vitiated co-flowing stream. The stabilization region of these flames is complex and, depending on the flow conditions, may undergo a transition from auto-ignition to premixed flame propagation. The objective of this paper is to develop numerical indicators for identifying such behavior, first in well-known simple test cases and then in the lifted turbulent flames. The calculations employ a composition probability density function (PDF) approach coupled to the commercial CFD code, FLUENT. The in-situ-adaptive tabulation (ISAT) method is used to implement detailed chemical kinetics. A simple k–ε turbulence model is used for turbulence along with a low Reynolds number model close to the solid walls of the fuel pipe.

The first indicator is based on an analysis of the species transport with respect to the budget of convection, diffusion and chemical reaction terms. This is a powerful tool for investigating aspects of turbulent combustion that would otherwise be prohibitive or impossible to examine experimentally. Reaction balanced by convection with minimal axial diffusion is taken as an indicator of auto-ignition while a diffusive–reactive balance, preceded by a convective–diffusive balanced pre-heat zone, is representative of a premixed flame. The second indicator is the relative location of the onset of creation of certain radical species such as HO₂ ahead of the flame zone. The buildup of HO₂ prior to the creation of H, O and OH is taken as another indicator of autoignition.

The paper first confirms the relevance of these indicators with respect to two simple test cases representing clear auto-ignition and premixed flame propagation. Three turbulent lifted flames are then investigated and the presence of auto-ignition is identified. These numerical tools are essential in providing valuable insights into the stabilization behaviour of these flames, and the demarcation between processes of auto-ignition and premixed flame propagation.     

Effects of gas and soot radiation on soot formation in counterflow ethylene diffusion flames

Numerical study of soot formation in counterflow ethylene diffusion flames at atmospheric pressure was conducted using detailed chemistry and complex thermal and transport properties. Soot kinetics was modelled using a semi-empirical two-equation model. Radiation heat transfer was calculated using the discrete-ordinates method coupled with an accurate band model. The calculated soot volume fractions are in reasonably good agreement with the experimental results in the literature. The individual effects of gas and soot radiation on soot formation were also investigated.     

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.     

Axis switching in turbulent buoyant diffusion flames

Dynamics of buoyant diffusion flames from rectangular, square, and round fuel sources were investigated using direct numerical simulation (DNS). Fully three-dimensional simulations were performed employing high-order numerical methods and boundary conditions to solve governing equations for variable-density flow and finite-rate Arrhenius chemistry. Significant differences among the different cases were revealed in the vortex dynamics, entrainment rate, small-scale mixing, and consequently flame structures. Mixing and entrainment enhancement in non-circular flames in comparison with circular ones was explained using the Biot–Savart instability theory, which relates vortex dynamics to the local azimuthal curvature. An extension of the theory elucidated why rectangular flames entrain more efficiently and spread wider than square ones, although both configurations have corners. It also provided an explanation for the aspect ratio effects in the near field. In the far field, nonlinear effects were dominant and the general transport equations for vorticity were analyzed in detail. The corner effects and aspect ratio effects were shown to be augmented by the intricate interactions among vortex dynamics, combustion, and buoyancy through the various terms in the equations. The presence of corners in non-circular flames led to concentrated regions of fine-scale mixing and intense reactions centered around the corners. Moreover, the rectangular flames exhibited a different dynamic behavior from even the square one, by creating discrepancies in entrainment, mixing, and combustion between the minor and major axis directions. Increasing the aspect ratio exacerbated such directional discrepancies, and ultimately led to axis switching. It was the first time that axis switching was observed by DNS in a rectangular flame of aspect ratio 3, which raised further questions in combustion prediction and control. Finally, a unified explanation for corner and aspect ratio effects was given on the basis of the Biot–Savart instability theory and the vorticity transport equations.     

Experimental measurements of geometric properties of turbulent stratified flames

Combustion under stratified conditions is common in many systems. However, relatively little is known about the structure and dynamics of turbulent stratified flames. Two-dimensional imaging diagnostics are applied to premixed and stratified V-flames at a mean equivalence ratio of 0.77, and low turbulent intensity, within the corrugated flame range. The present results show that stratification affects the mean turbulent flame speed, structure and geometric properties. Stratification increases the flame surface density above the premixed flame levels in all cases, with a maximum reached at intermediate levels of stratification. The flame surface density (FSD) of stratified flames is higher than that of premixed flames at the same mean equivalence ratio. Under the present conditions, the FSD peaks at a stratification ratio around 3.0. The FSD curves for stratified flames are further skewed towards the product side. The distribution of flame curvature in stratified flames is broader and more symmetric relative to premixed flames, indicating an additional mechanism of curvature generation, which is not necessarily due to cusping. These experiments indicate that flame stratification affects the intrinsic behaviour of turbulent flames and suggest that models may need to be revised in the light of the current evidence.     

The stabilization mechanism and structure of turbulent hydrocarbon lifted flames

The structure and stabilization mechanism of turbulent lifted non-premixed hydrocarbon flames have been investigated using combined laser imaging techniques. The techniques include Rayleigh scattering, laser induced predissociation fluorescence of OH, LIF of PAH, LIF of CH₂O, and planar imaging velocimetry. The geometrical structure of multi-reaction zones and flow field at the stabilization region have been simultaneously measured in 16 hydrocarbon flames. The data reveal the existence of triple flame structure at the stabilization region of turbulent lifted flames. Increasing the jet velocity leads to an increase of the lift-off height and to a broadening of the lift-off region. Further analysis of the stabilization criterion at the lift-off height based on the premixed nature of triple-flame propagation and flow field data has been presented and discussed.     

Detailed measurements of local scalar front structures in stagnation-type turbulent premixed flames

Local scalar front structures of OH mole fraction, reaction progress variable, and its three-dimensional gradient have been measured in stagnation-type turbulent premixed flames. The reaction progress variable front is observed to change with increasing turbulence from parallel iso-scalar contours but reduced progress variable gradients, called the lamella-like front, to disrupted non-parallel iso-contours that deviate substantially from those of wrinkled laminar flamelets, called the non-flamelet front. This transition is attributed to the different scales of interaction between the flame internal structure and a spectrum of turbulence extending from the integral scale to the Kolmogorov scale. The lamella-like front pattern occurs when the length scales of interaction are smaller than the laminar flame thickness but the time scales are greater than the flame residence time. The non-flamelet front pattern occurs when the length scales of interaction are greater than the laminar flame thickness but the time scales are smaller than the flame residence time. This difference corresponds to the change of combustion regime from complex-strain flame front to turbulent flame front on a revised regime diagram. A correlation is also proposed for the turbulent flame brush thickness as a function of turbulent Reynolds number and heat release parameter. The heat release parameter is considered to arise from the non-passive effects of flame-surface wrinkling.