Soot formation and temperature field structure in co-flow laminar methane-air diffusion flames at pressures from 10 to 60 atm

The effects of pressure on soot formation and the structure of the temperature field were studied in co-flow methane-air laminar diffusion flames over a wide pressure range, from 10 to 60 atm in a high-pressure combustion chamber. The selected fuel mass flow rate provided diffusion flames in which the soot was completely oxidized within the visible flame envelope and the flame was stable at all pressures considered. The spatially resolved soot volume fraction and soot temperature were measured by spectral soot emission as a function of pressure. The visible (luminous) flame height remained almost unchanged from 10 to 100 atm. Peak soot concentrations showed a strong dependence on pressure at relatively lower pressures; but this dependence got weaker as the pressure is increased. The maximum conversion of the fuel’s carbon to soot, 12.6%, was observed at 60 atm at approximately the mid-height of the flame. Radial temperature gradients within the flame increased with pressure and decreased with flame height above the burner rim. Higher radial temperature gradients near the burner exit at higher pressures mean that the thermal diffusion from the hot regions of the flame towards the flame centerline is enhanced. This leads to higher fuel pyrolysis rates causing accelerated soot nucleation and growth as the pressure increases.     

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.     

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.     

The importance of thermal radiation transfer in laminar diffusion flames at normal and microgravity

The importance of radiation heat loss in laminar and turbulent diffusion flames at normal gravity has been relatively well recognized in recent years. There is currently lack of quantitative understanding on the importance of radiation heat loss in relatively small scale laminar diffusion flames at microgravity. The effects of radiation heat transfer and radiation absorption on the structure and soot formation characteristics of a coflow laminar ethylene/air diffusion flame at normal- and microgravity were numerically investigated. Numerical calculations were conducted using GRI-Mech 3.0 combustion chemistry without the NO_x mechanism and complex thermal and transport properties, an acetylene based soot formation model, and a statistical narrow-band correlated-k non-grey gas radiation model. Radiation heat transfer and radiation absorption in the microgravity flame were found to be much more important than their counterparts at normal gravity. It is important to calculate thermal radiation transfer accurately in diffusion flame modelling under microgravity conditions.     

TDLAS-based in situ measurement of absolute acetylene concentrations in laminar 2D diffusion flames

We report the first quantitative and calibration-free in situ C₂H₂ measurement in a flame environment using direct Tunable Diode Laser Absorption Spectroscopy(TDLAS). Utilizing a fiber-coupled Distributed Feedback diode laser near 1535 nm we measured spatially resolved, absolute C₂H₂ concentration profiles in a laminar non-premixed CH₄/air flame supported on a modified Wolfhard-Parker slot burner with N₂ purge slots to minimize end flames. We developed a wavelength tuning scheme combining laser temperature and current modulation to record with a single diode laser a mesh of 37 overlapping spectral windows and generate an ∼7 nm (30 cm⁻¹) wide, high-resolution absorption spectrum centered at 1538 nm. Experimental C₂H₂ spectra in a reference cell showed excellent agreement with simulations using HITRAN2004 data. The enhanced wavelength coverage was needed to establish correct C₂H₂ line identification and selection in the very congested high temperature flame spectra and led to the P17e line as the only candidate for in situ detection of C₂H₂ in the flame. We used highly efficient optical disturbance correction algorithms for treating transmission and background emission fluctuations in combination with a multiple Voigt line Levenberg-Marquardt fitting algorithm and Pt/Rh thermocouple measurements to achieve fractional optical resolutions of up to 4 × 10⁻⁵ OD (1σ) in the flame (T up to 2000 K). Temperature dependent C₂H₂ detection limits for the P17e line were 60 to 480 ppm. By translating the burner through the laser beam with a DC motor we obtained spatially resolved, absolute C₂H₂ concentration profiles along the flame sheet with 0.5 mm spatial resolution as measured with a knife edge technique. The TDLAS-based, transverse C₂H₂ concentration profiles without any scaling are in excellent agreement with published mass spectrometric C₂H₂ data for the same flame supported on a similar burner, thus validating our calibration-free TDLAS measurements.     

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.     

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.     

Constrained-temperature solutions of coflow laminar diffusion flames

We investigate the use of an experimental 2-D temperature profile to constrain detailed numerical solutions of a sooting coflow laminar diffusion flame. Experimentally, four optical diagnostic techniques are used to measure the two-dimensional temperature field in an ethylene-air coflow flame. This experimental temperature field is then used to impose the temperature in the solution process, thus obviating the need to solve the energy equation and, in particular, to incorporate costly models of radiative losses in the flame. Results are presented for a 40% ethylene-air flame on the Yale Coflow Burner. In the unconstrained solution of the complete set of governing equations, the location of maximum temperature is found along the flame wings, whereas the experimental temperature field has its maximum along the centerline. Similarly, the location of peak soot volume fraction migrates from along the flame wings in the unconstrained calculation, where soot surface growth processes dominate, to the centerline in the constrained case, where soot inception is the dominant condensed-phase formation mechanism. The distribution of soot in the constrained solution is much more consistent with experimental observations, and this fact illustrates how the validation of a soot sub-model may be complicated by the necessity of modeling distributed heat losses in the flame.     

Transport mechanisms controlling soot production inside a non-buoyant laminar diffusion flame

This study integrates new and existing numerical modeling and experimental observations to provide a consistent explanation to observations pertaining flame length and soot volume fractions for laminar diffusion flames. Integration has been attempted by means of scaling analysis. Emphasis has been given to boundary layer flames. For the experiments, ethylene is injected through a flat porous burner into an oxidizer flowing parallel to the burner surface. The oxidizer is a mixture of oxygen and nitrogen, flowing at various velocities. All experiments were conducted in microgravity to minimize the role of buoyancy in distorting the aerodynamics of the flames. A previous numerical study emphasizing fuel transport was extended to include the oxidizer flow. Fictitious tracer particles were used to establish the conditions in which fuel and oxidizer interact. This allowed establishing regions of soot formation and oxidation as well as relevant characteristic length and time scales. Adequate scaling parameters then allow to establish explanations that are consistent for different burner configurations as well as “open-tip” and “closed-tip” flames.     

Characteristics of microjet methane diffusion flames

Characteristics of microjet methane diffusion flames stabilized on top of the vertically oriented, stainless-steel tubes with an inner diameter ranging from 186 to 778 μ m are investigated experimentally, theoretically and numerically. Of particular interest are the flame shape, flame length and quenching limit, as they may be related to the minimum size and power of the devices in which such flames would be used for future micro-power generation. Experimental measurements of the flame shape, flame length and quenching velocity are compared with theoretical predictions as well as detailed numerical simulations. Comparisons of the theoretical predictions with measured results show that only Roper's model can satisfactorily predict the flame height and quenching velocity of microjet methane flames. Detailed numerical simulations, using skeletal chemical kinetic mechanism, of the flames stabilized at the tip of d = 186, 324 and 529 μ m tubes are performed to investigate the flame structures and the effects of burner materials on the standoff distance near extinction limit. The computed flame shape and flame length for the d = 186 μm flame are in excellent agreement with experimental results. Numerical predictions of the flame structures strongly suggest that the flame burns in a diffusion mode near the extinction limit. The calculated OH mass fraction isopleths indicate that different tube materials have a minor effect on the standoff distance, but influence the quenching gap between the flame and the tube.     

Head-on quenching of laminar premixed methane flames diluted with hot combustion products

Transient head-on quenching of laminar premixed methane flames diluted with hot combustion products is analyzed using full-chemistry 1D DNS. The impact of the dilution level, pressure and wall temperature on carbon monoxide (CO) emissions is investigated. Increasing dilution level and pressure reduce peak average near-wall CO concentrations, and reduce the near-wall CO reduction rate. However, the peak average near-wall CO and near-wall CO reduction rate increase with increasing wall temperature. Analysis of the species transport budget for CO near the wall before, during and after quenching indicates that there are conditions where diffusion is the dominant transport term. As a consequence, it may be possible to model the near-wall CO using only the integrated diffusion term within certain spatial regions. Dilution increases the size of these regions, whereas increasing pressure reduces this size.     

Use of laser-induced spark for studying ignition stability and unburned hydrogen escaping from laminar diluted hydrogen diffusion jet flames

Ignition and unburned hydrogen escaping from hydrogen jet diffusion flames diluted with nitrogen up to 70% were experimentally studied. The successful ignition locations were about 2/3 of the flame length above the jet exit for undiluted flames and moved much closer to the exit for diluted flames. For higher levels of dilution or higher flow rates, there existed a region within which a diluted hydrogen diffusion flame can be ignited and burns with a stable liftoff height. This is contrary to previous findings that pure and diluted hydrogen jet diffusion cannot achieve a stable lifted flame configuration. With liftoff, the flame is noisy and short with significant amount of unburned hydrogen escaping into the product gases. If ignition is initiated below this region, the flame propagates upstream quickly and attaches to the burner rim. Results from measurements of unburned hydrogen in the combustion products showed that the amount of unburned hydrogen increased as the nitrogen dilution level was increased. Thus, hydrogen diffusion flame diluted with nitrogen cannot burn completely.     

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.     

Calculating the soot particle size distribution function in turbulent diffusion flames using a sectional method

Soot formation in a turbulent jet diffusion flame is modeled using an unsteady flamelet approach in post-process. In the present work, we apply a detailed kinetic soot model with a sectional method, and study the evolution of the particle size distribution. Detailed information on the evolution of the soot particle size distribution function is acquired. It is found that the particle size distribution function is bimodal throughout the flame. The transition from the small to large particle size distributions is strongly influenced by surface growth and oxidation reactions. We find that large particles are most likely to be emitted from the flame.     

The effect of preferential diffusion on soot formation in a laminar ethylene/air diffusion flame

The influence of preferential diffusion on soot formation in a laminar ethylene/air diffusion flame was investigated by numerical simulation using three different transport property calculation methods. One simulation included preferential diffusion and the other two neglected preferential diffusion. The results show that the neglect of preferential diffusion or the use of unity Lewis number for all species results in a significant underprediction of soot volume fraction. The peak soot volume fraction is reduced from 8.0 to 2.0 ppm for the studied flame when preferential diffusion is neglected in the simulation. Detailed examination of numerical results reveals that the underprediction of soot volume fraction in the simulation neglecting preferential diffusion is due to the slower diffusion of some species from main reaction zone to PAH and soot formation layer. The slower diffusion of these species causes lower PAH formation rate and thus results in lower soot inception rate and smaller particle surface area. The smaller surface area further leads to smaller surface growth rate. In addition, the neglect of preferential diffusion also leads to higher OH concentration in the flame, which causes the higher specific soot oxidation rate. The lower inception rate, smaller surface growth rate and higher specific oxidation rate results in the lower soot volume fraction when preferential diffusion is neglected. The finding of the paper implies the importance of preferential diffusion for the modeling of not only laminar but maybe also some turbulent flames.     

Computational and experimental investigation of the interaction of soot and NO in coflow diffusion flames

A combined computational and experimental investigation that examines the relationship of soot formation and NO in coflow ethylene air diffusion flames is presented. While both NO and soot formation are often studied independently, there is a need to understand their coupled relationship as a function of system parameters such as fuel type, temperature and pressure. The temperature decrease due to radiative losses in systems in which significant soot is produced can affect flame length and other temperature-dependent processes such as the formation of NO. The results of a computational model that includes a sectional representation for soot formation with a radiation model are compared against laser-induced fluorescence measurements of NO. The sooting characteristics of these flames have been studied previously. Experimentally, a laser near 225.8 nm is used to excite the γ(0, 0) band in NO. Spectrally resolved fluorescence emission is imaged radially, for the (0, 0), (0, 1), (0, 2), (0, 3), and (0, 4) vibrational bands, at varying axial heights to create a two-dimensional image of NO fluorescence. A reverse quenching correction is applied to the computational results to determine an expected fluorescence signal for comparison with experimental results. Modeling results confirm that Fenimore NO is the dominant mechanism for NO production and suggest that for lightly sooting flames (peak soot volume fraction < 0.5 ppm), soot reduces only the Zeldovich NO formation (by a factor of two). For flames with increased soot levels (peak soot volume fraction ∼ 4 ppm), the model indicates not only that Zeldovich NO decreases by a factor of 2.5 through radiation loss, but that non-Zeldovich NO is reduced in the top center of the flame by about 30% through the oxidation of soot.     

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.